Drug release enhanced by temperature: An accurate discrete model for solutions in H3

In this paper we consider the coupling between two quasilinear diffusion equations: the diffusion coefficient of the first equation depends on its solution and the diffusion and convective coefficients of the second equation depend on the solution of the first one. This system can be used to describe the drug evolution in a target tissue when the drug transport is enhanced by heat. We study, from an analytical and a numerical viewpoints, the coupling of the heat equation with the drug diffusion equation. A fully discrete piecewise linear finite method is proposed to solve this system and its stability is studied. Assuming that the heat and the concentration are in H3 we prove that the method is second order convergent. Numerical experiments illustrating the theoretical results and the global qualitative behaviour of the coupling are also included.