Abstract:Ion transport, often described by the Poisson--Nernst--Planck (PNP) equations, is ubiquitous in electrochemical devices and many biological processes of significance. In this work, we develop conservative, positivity-preserving, energy dissipating, and implicit finite difference schemes for solving the multi-dimensional PNP equations with multiple ionic species. A central-differencing discretization based on harmonic-mean approximations is employed for the Nernst--Planck (NP) equations. The backward Euler discretization in time is employed to derive a fully implicit nonlinear system, which is efficiently solved by a newly proposed Newton's method. The improved computational efficiency of the Newton's method originates from the usage of the electrostatic potential as the iteration variable, rather than the unknowns of the nonlinear system that involves both the potential and concentration of multiple ionic species. Numerical analysis proves that the numerical schemes respect three desired analytical properties (conservation, positivity preserving, and energy dissipation) fully discretely. Based on advantages brought by the harmonic-mean approximations, we are able to establish estimate on the upper bound of condition numbers of coefficient matrices in linear systems that are solved iteratively. The solvability and stability of the linearized problem in the Newton's method are rigorously established as well. Numerical tests are performed to confirm the anticipated numerical accuracy, computational efficiency, and structure-preserving properties of the developed schemes. Adaptive time stepping is implemented for further efficiency improvement. Finally, the proposed numerical approaches are applied to characterize ion transport subject to a sinusoidal applied potential.

A positivity-preserving, linear, energy stable and convergent numerical scheme for the Poisson–Nernst–Planck (PNP) system

Efficient, positive, and energy stable schemes for multi-D Poisson-Nernst-Planck systems

An unconditionally energy stable linear scheme for Poisson–Nernst–Planck equations

A Second Order Accurate, Positivity Preserving Numerical Method for the Poisson–Nernst–Planck System and Its Convergence Analysis

A positivity preserving scheme for Poisson-Nernst-Planck Navier-Stokes equations and its error analysis

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

A Linearized Local Conservative Mixed Finite Element Method for Poisson–Nernst–Planck Equations

Unconditional positivity-preserving and energy stable schemes for a reduced Poisson-Nernst-Planck system

Second-order, Positive, and Unconditional Energy Dissipative Scheme for Modified Poisson–Nernst–Planck Equations

Second-order, Positive, and Unconditional Energy Dissipative Scheme for Modified Poisson-Nernst-Planck Equations

A second order accurate numerical method for the Poisson-Nernst-Planck system in the energetic variational formulation

Structure-Preserving and Efficient Numerical Methods for Ion Transport

A linear, second-order accurate, positivity-preserving and unconditionally energy stable scheme for the Navier–Stokes–Poisson–Nernst–Planck system

A stabilized finite volume element method for solving Poisson–Nernst–Planck equations

A Positivity-Preserving and Free Energy Dissipative Hybrid Scheme for the Poisson-Nernst-Planck Equations on Polygonal and Polyhedral Meshes

Superconvergence analysis of finite element method for Poisson–Nernst–Planck equations

Nonconforming finite element method for coupled Poisson–Nernst–Planck equations

Unconditional superconvergence analysis of a structure-preserving finite element method for the Poisson-Nernst-Planck equations

Positivity Preserving and Mass Conservative Projection Method for the Poisson–Nernst–Planck Equation

Convergent discretizations for the Nernst–Planck–Poisson system

Convergence analysis of a positivity-preserving numerical scheme for the Cahn-Hilliard-Stokes system with Flory-Huggins energy potential