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 $H^{-1}$ 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.

Numerical methods for solution of the stochastic differential equations equivalent to the non-stationary Parker's transport equation

A stochastic method of solution of the Parker transport equation

Numerical Solution of Differential Equations by the Parker-Sochacki Method

Numerical solutions to time-fractional stochastic partial differential equations

Numerical approximation of singular-degenerate parabolic stochastic partial differential equations

Numerical analysis of the stochastic Stefan problem

Stochastic methods for solving high-dimensional partial differential equations

Efficient numerical methods for Anisotropic Diffusion of Galactic Cosmic Rays

Fractional Fokker-Planck Equation with Space and Time Dependent Drift and Diffusion

Numerical Solution of Stochastic Partial Differential Equations with Correlated Noise

A survey of numerical methods for stochastic differential equations

A stochastic solution with Gaussian stationary increments of the symmetric space-time fractional diffusion equation

Numerical Solution of Stochastic Differential Equations: Diffusion and Jump-Diffusion Processes

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

Numerical studies of a class of reaction–diffusion equations with Stefan conditions

A class of Galerkin schemes for time-dependent radiative transfer

Higher order time discretization method for a class of semilinear stochastic partial differential equations with multiplicative noise

Stabilized Time Series Expansions for High-Order Finite Element Solutions of Partial Differential Equations

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

Modeling Cosmic-Ray Transport: A CRPropa based stochastic differential equation solver

Numerical Methods of the Maxwell-Stefan Diffusion Equations and Applications in Plasma and Particle Transport