We present a multi-layer mathematical model to describe the transdermal drug release from an iontophoretic system. The Nernst-Planck equation describes the basic convection-diffusion process, with the electric potential obtained by solving the Laplace's equation. These equations are complemented with suitable interface and boundary conditions in a multi-domain. The stability of the mathematical problem is discussed in different scenarios and a finite-difference method is used to solve the coupled system. Numerical experiments are included to illustrate the drug dynamics under different conditions.
Keywords: Nernst–Plank equations; drug release; finite-difference methods; iontophoresis.
© The authors 2016. Published by Oxford University Press on behalf of the Institute of Mathematics and its Applications. All rights reserved.