Positivity-preserving third order DG schemes for Poisson--Nernst--Planck equations
arXiv:2102.00101 · doi:10.1016/j.jcp.2021.110777
Abstract
In this paper, we design and analyze third order positivity-preserving discontinuous Galerkin (DG) schemes for solving the time-dependent system of Poisson--Nernst--Planck (PNP) equations, which has found much use in diverse applications. Our DG method with Euler forward time discretization is shown to preserve the positivity of cell averages at all time steps. The positivity of numerical solutions is then restored by a scaling limiter in reference to positive weighted cell averages. The method is also shown to preserve steady states. Numerical examples are presented to demonstrate the third order accuracy and illustrate the positivity-preserving property in both one and two dimensions.
7 figures, 16 tables
References in corpus (10)
- Ionic channels in biological membranes. Electrostatic analysis of a natural nanotube
- A discontinuous Galerkin method for nonlinear parabolic equations and gradient flow problems with interaction potentials
- 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
- A free energy satisfying discontinuous Galerkin method for one-dimensional Poisson--Nernst--Planck systems
- An energy preserving finite difference scheme for the Poisson-Nernst-Planck system
- Invariant-region-preserving DG methods for multi-dimensional hyperbolic conservation law systems, with an application to compressible Euler equations
- Structure-Preserving and Efficient Numerical Methods for Ion Transport
- An invariant-region-preserving limiter for DG schemes to isentropic Euler equations
- THIRD ORDER MAXIMUM-PRINCIPLE-SATISFYING DG SCHEMES Third Order Maximum-Principle-Satisfying DG schemes for Convection-Diffusion problems with Anisotropic Diffusivity DIFFUSIVITY