In this talk we propose and validate a multiscale model for a Poisson-Nernst-Planck (PNP) system, focusing on the correlated motion of positive and negative ions under the influence of a (potentially vibrating) trap. Specifically, we aim to model surface traps whose attraction range, of length delta, is much smaller than 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 meet the surface of the trap, they are reversibly adsorbed. As in our previous works, the effect of the attractive potential is replaced by a suitable boundary condition derived by mass conservation and asymptotic analysis. Since the potential range is delta, an asymptotic expansion is considered, and the boundary conditions are obtained by retaining the lowest order terms in the expansion. This model is particularly innovative as it addresses the interaction between two carriers-positive and negative ions-simultaneously, which is often neglected in traditional approaches that treat ion species independently. In the second part of the talk, we address the Coulomb interaction between carriers. For very short Debye lengths, of length epsilon, one can adopt the so called Quasi-Neutral limit which drastically simplifies the system, reducing it to a diffusion equation for a single carries with effective diffusion coefficient. This approach, while simplifying the mathematical model, can limit the scope of numerical simulations, as it may not capture the full range of behaviors near the Quasi-Neutral limit. The objective of this work is to develop an Asymptotic Preserving (AP) numerical scheme that remains stable as the Debye length tends to zero, with no stability restriction in the time step. Our AP scheme is designed to handle both regimes: the full system when the Debye length is small but non-negligible, and the Quasi-Neutral regime as the Debye length approaches zero. In our case, the numerical scheme we propose is also Asymptotic Accurate (AA). This means that, in addition to preserving stability, the scheme remains uniformly accurate in the Quasi-Neutral limit.