Abstract:A second order accurate (in time) numerical scheme is proposed and analyzed for the Poisson–Nernst–Planck equation (PNP) system, reformulated as a non-constant mobility gradient flow in the Energetic Variational Approach (EnVarA). The centered finite difference is taken as the spatial discretization. Meanwhile, the highly nonlinear and singular nature of the logarithmic energy potentials has always been the essential difficulty to design a second order accurate scheme in time, while preserving the variational energetic structures. The mobility function is updated with a second order accurate extrapolation formula, for the sake of unique solvability. A modified Crank–Nicolson scheme is used to approximate the logarithmic term, so that its inner product with the discrete temporal derivative exactly gives the corresponding nonlinear energy difference; henceforth the energy stability is ensured for the logarithmic part. In addition, nonlinear artificial regularization terms are added in the numerical scheme, so that the positivity-preserving property could be theoretically proved, with the help of the singularity associated with the logarithmic function. Furthermore, an optimal rate convergence analysis is provided in this paper, in which the higher order asymptotic expansion for the numerical solution, the rough error estimate and refined error estimate techniques have to be included to accomplish such an analysis. This work combines the following theoretical properties for a second order accurate numerical scheme for the PNP system: (i) second order accuracy in both time and space, (ii) unique solvability and positivity, (iii) energy stability, and (iv) optimal rate convergence. A few numerical results are also presented.

Convergent discretizations for the Nernst–Planck–Poisson system

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

Adaptive finite element approximation for steady-state Poisson-Nernst-Planck equations

Discretizing the Fokker-Planck equation with second-order accuracy: a dissipation driven approach

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

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

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

A finite volume scheme for boundary-driven convection–diffusion equations with relative entropy structure

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

Convergence of an implicit Euler Galerkin scheme for Poisson-Maxwell-Stefan systems

Physics-based stabilized finite element approximations of the Poisson--Nernst--Planck equations

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

The Entropy Satisfying Discontinuous Galerkin Method for Fokker–Planck Equations

Consistency and convergence for a family of finite volume discretizations of the Fokker--Planck operator

A finite difference scheme for two-dimensional singularly perturbed convection-diffusion problem with discontinuous source term

Long-time behavior for discretization schemes of Fokker-Planck equations via couplings

Convergence of the numerical approximations and well-posedness: Nonlocal conservation laws with rough flux

Convergence of a finite volume scheme and dissipative measure-valued-strong stability for a hyperbolic-parabolic cross-diffusion system

Convergence of a continuous Galerkin method for hyperbolic-parabolic systems

A convergent finite-volume scheme for nonlocal cross-diffusion systems for multi-species populations

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