Efficient, positive, and energy stable schemes for multi-D Poisson-Nernst-Planck systems
arXiv:2001.08350
Abstract
In this paper, we design, analyze, and numerically validate positive and energy-dissipating schemes for solving the time-dependent multi-dimensional system of Poisson-Nernst-Planck (PNP) equations, which has found much use in the modeling of biological membrane channels and semiconductor devices. The semi-implicit time discretization based on a reformulation of the system gives a well-posed elliptic system, which is shown to preserve solution positivity for arbitrary time steps. The first order (in time) fully-discrete scheme is shown to preserve solution positivity and mass conservation unconditionally, and energy dissipation with only a mild time step restriction. The scheme is also shown to preserve the steady-state. For the fully second order (in both time and space) scheme with large time steps, solution positivity is restored by a local scaling limiter, which is shown to maintain the spatial accuracy. These schemes are easy to implement. Several three-dimensional numerical examples verify our theoretical findings and demonstrate the accuracy, efficiency, and robustness of the proposed schemes, as well as the fast approach to steady states.
arXiv comment32 pages, 3 tables, 4 figures
References in corpus (7)
- Diffuse-Charge Dynamics in Electrochemical Systems
- Cross-interaction drives stratification in drying film of binary colloidal mixtures
- A free energy satisfying finite difference method for Poisson--Nernst--Planck equations
- Numerical Solution of 3D Poisson-Nernst-Planck Equations Coupled with Classical Density Functional Theory for Modeling Ion and Electron Transport in a Confined Environment
- An energy preserving finite difference scheme for the Poisson-Nernst-Planck system
- Positive and free energy satisfying schemes for diffusion with interaction potentials
- Unconditional positivity-preserving and energy stable schemes for a reduced Poisson-Nernst-Planck system