Unconditionally positivity preserving and energy dissipative schemes for Poisson--Nernst--Planck equations

Abstract

We develop a set of numerical schemes for the Poisson--Nernst--Planck equations. We prove that our schemes are mass conservative, uniquely solvable and keep positivity unconditionally. Furthermore, the first-order scheme is proven to be unconditionally energy dissipative. These properties hold for various spatial discretizations. Numerical results are presented to validate these properties. Moreover, numerical results indicate that the second-order scheme is also energy dissipative, and both the first- and second-order schemes preserve the maximum principle for cases where the equation satisfies the maximum principle.

24 pages, 10 figures

References in corpus (1)

• Ionic channels in biological membranes. Electrostatic analysis of a natural nanotube

Cited by in corpus (8)

• A structure-preserving, operator splitting scheme for reaction-diffusion equations with detailed balance
• High order spatial discretization for variational time implicit schemes: Wasserstein gradient flows and reaction-diffusion systems
• A dynamic mass transport method for Poisson-Nernst-Planck equations
• High-order space-time finite element methods for the Poisson-Nernst-Planck equations: Positivity and unconditional energy stability
• Structure-preserving numerical method for Maxwell-Ampère Nernst-Planck model
• Efficient numerical methods for the Navier-Stokes-Nernst-Planck-Poisson equations
• Local Structure-Preserving Relaxation Method for Charged Systems on Unstructured Meshes
• $Q$-tensor gradient flow with quasi-entropy and discretizations preserving physical constraints