An energy-preserving discretization for the Poisson–Nernst–Planck equations
Journal of Computational Electronics
Poisson–Nernst–Planck equations, Energy conservation, Finite difference, Ion channels
The Poisson–Nernst–Planck (PNP) equations have recently been used to describe the dynamics of ion transport through biological ion channels besides being widely employed in semiconductor industry. This paper is about the design of a numerical scheme to solve the PNP equations that preserves exactly (up to roundoff error) a discretized form of the energy dynamics of the system. The proposed finite difference scheme is of second-order accurate in both space and time. Comparisons are made between this energy dynamics-preserving scheme and a standard finite difference scheme, showing a difference in satisfying the energy law. Numerical results are presented for validating the orders of convergence in both time and space of the new scheme for the PNP system.
Flavell, Allen; Kabre, Julienne; and Li, Xiaofan, "An energy-preserving discretization for the Poisson–Nernst–Planck equations" (2017). Mathematics Faculty Articles. 320.