Asymptotic Preserving and Accurate scheme for Multiscale Poisson-Nernst-Planck (MPNP) system

Abstract: In this paper, we propose and validate a two-species Multiscale model for a Poisson-Nernst-Planck (PNP) system, focusing on the correlated motion of positive and negative ions under the influence of a trap. Specifically, we aim to model surface traps whose attraction range, of length delta, is much smaller then the scale of the problem. The physical setup we refer to is an anchored gas drop (bubble) surrounded by a diffusive flow of charged surfactants (ions). When the diffusing surfactants reach the surface of the trap, the anions are adsorbed. As in our previous works [11,6,9,4], the effect of the attractive potential is replaced by a suitable boundary condition derived by mass conservation and asymptotic analysis. The novelty of this work is the extension of the model proposed in [11], now incorporating the influence of both carriers - positive and negative ions - simultaneously, which is often neglected in traditional approaches that treat ion species independently. In the second part of the paper, we address the treatment of the Coulomb interaction between carriers. When the Debye length lambda_D (proportional to a small parameter epsilon) is very small, one can adopt the so-called Quasi-Neutral limit, which significantly simplifies the system, reducing it to a diffusion equation for a single carriers with effective diffusion coefficient [36,53]. This approach, while simplifying the mathematical model, does not capture the effects of non negligible values of epsilon. When the Debye length is small but not negligible, it may be very expensive to capture the small deviation from the Quasi-Neutral limit by standard methods in the literature. [...]