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.

Convergence of the Crank-Nicolson/newton Scheme for Nonlinear Parabolic Problem

Optimal Convergence of the Newton Iterative Crank–Nicolson Finite Element Method for the Nonlinear Schrödinger Equation

Fractional Crank-Nicolson Galerkin finite element analysis for coupled time-fractional nonlocal parabolic problem

Optimal $L^2$ error estimates of mass- and energy-conserved FE schemes for a nonlinear Schrödinger-type system

Conservative higher-order finite difference scheme for the coupled nonlinear Schrödinger equations

Second-order accurate numerical scheme with graded meshes for the nonlinear partial integrodifferential equation arising from viscoelasticity

Convergence analysis and error estimates for a second order accurate finite element method for the Cahn–Hilliard–Navier–Stokes system

A conservative relaxation Crank-Nicolson finite element method for the Schrödinger-Poisson equation

Unconditionally optimal H1-norm error estimates of a fast and linearized Galerkin method for nonlinear subdiffusion equations

Efficient second-order semi-implicit finite element method for fourth-order nonlinear diffusion equations

Convergence analysis of finite element approximations for a nonlinear second order hyperbolic optimal control problems

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

Superconvergence analysis of an anisotropic nonconforming FEM for parabolic equation

Improving Newton's method performance by parametrization: the case of Richards equation

Unconditional Superconvergence Analysis of a Crank–Nicolson Galerkin FEM for Nonlinear Schrödinger Equation

Nonconforming spectral element approximation for parabolic PDE with corner singularity

Error estimation for non-linear singularly perturbed reaction–diffusion parabolic problems via element-free Galerkin method

Two-grid reduced-dimension extrapolated method of Crank-Nicolson finite element solution coefficient vectors for nonlinear parabolic equation

Long-time Convergence of Numerical Approximations for Semilinear Parabolic Equations (Ⅰ)

High-order weighted compact nonlinear scheme for solving degenerate parabolic equations

A New Error Analysis of Crank–Nicolson Galerkin FEMs for a Generalized Nonlinear Schrödinger Equation