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.

A Conservative, Free Energy Dissipating, and Positivity Preserving Finite Difference Scheme for Multi-Dimensional Nonlocal Fokker–Planck Equation

Positivity-preserving and energy-dissipating discontinuous Galerkin methods for nonlinear nonlocal Fokker-Planck equations

Nonstandard finite difference schemes for linear and non-linear Fokker–Planck equations

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

A Linear Finite Difference Scheme for the Two-Dimensional Nonlinear Schrödinger Equation with Fractional Laplacian

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

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

The Entropy Satisfying Discontinuous Galerkin Method for Fokker–Planck Equations

An optimization-based positivity-preserving limiter in semi-implicit discontinuous Galerkin schemes solving Fokker-Planck equations

A Linearly-Fitted Conservative (dissipative) Scheme for Efficiently Solving Conservative (dissipative) Nonlinear Wave PDEs

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

A Numerical Scheme for the Quantum Fokker-Planck-Landau Equation Efficient in the Fluid Regime

A Positivity-Preserving Finite Volume Scheme for Three-Temperature Radiation Diffusion Equations

On positivity preservation in nonlinear finite volume method for multi-term fractional subdiffusion equation on polygonal meshes

Numerical Analysis and Computation of the Finite Volume Element Method for the Nonlinear Coupled Time-Fractional Schrödinger Equations

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

Unconditionally Positivity and Boundedness Preserving Schemes for a FitzHugh–Nagumo Equation

A dual-pairing summation-by-parts finite difference framework for nonlinear conservation laws

A class of asymptotic-preserving schemes for the Fokker-Planck-Landau equation

Finite Element Scheme for Fokker-Planck Equation in Plasmas

A conservative, implicit solver for 0D-2V multi-species nonlinear Fokker-Planck collision equations