Abstract:In this paper we propose and analyze a finite difference numerical scheme for the Poisson-Nernst-Planck equation (PNP) system. To understand the energy structure of the PNP model, we make use of the Energetic Variational Approach (EnVarA), so that the PNP system could be reformulated as a non-constant mobility H − 1 H^{-1} gradient flow, with singular logarithmic energy potentials involved. To ensure the unique solvability and energy stability, the mobility function is explicitly treated, while both the logarithmic and the electric potential diffusion terms are treated implicitly, due to the convex nature of these two energy functional parts. The positivity-preserving property for both concentrations, n n and p p , is established at a theoretical level. This is based on the subtle fact that the singular nature of the logarithmic term around the value of 0 0 prevents the numerical solution reaching the singular value, so that the numerical scheme is always well-defined. In addition, an optimal rate convergence analysis is provided in this work, in which many highly non-standard estimates have to be involved, due to the nonlinear parabolic coefficients. The higher order asymptotic expansion (up to third order temporal accuracy and fourth order spatial accuracy), the rough error estimate (to establish the ℓ ∞ \ell ^\infty bound for n n and p p ), and the refined error estimate have to be carried out to accomplish such a convergence result. In our knowledge, this work will be the first to combine the following three theoretical properties for a numerical scheme for the PNP system: (i) unique solvability and positivity, (ii) energy stability, and (iii) optimal rate convergence. A few numerical results are also presented in this article, which demonstrates the robustness of the proposed numerical scheme.

Well-Posedness For The Navier-Stokes-Nernst-Planck-Poisson System In Triebel-Lizorkin Space And Besov Space With Negative Indices

Well-posedness in Critical Spaces for the Density-dependent Incompressible Viscoelastic Fluid System

Global Well-Posedness for the Navier–Stokes–Nernst–Planck–Poisson System in Dimension Two

On the Well-Posedness and Decay Rates of Solutions to the Poisson–Nernst–Planck–Navier–Stokes System

Global Well-Posedness of A Dissipative System Arising in Electrohydrodynamics in Negative-Order Besov Spaces

Well-posedness of the viscous boussinesq system in besov spaces of negative regular index s = -1

OPTIMAL LARGE TIME BEHAVIOR OF THE COMPRESSIBLE BIPOLAR NAVIER-STOKES-POISSON SYSTEM

Well-posedness for the Euler-Nernst-Planck-Possion system in Besov spaces

On the Well-Posedness of 2-D Incompressible Navier-Stokes Equations with Variable Viscosity in Critical Spaces

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

Symmetry and Non-Existence of Positive Solutions for Pde System with Navier Boundary Conditions on A Half Space

Optimal Large Time Behavior of the Compressible Bipolar Navier–Stokes–Poisson System with Unequal Viscosities

Global Well-Posedness for Navier–Stokes Equations with Small Initial Value in B^0_n,∞(Ω)

Global Existence of Large Solutions for the 3D incompressible Navier--Stokes--Poisson--Nernst--Planck Equations

Global Well-Posedness for Navier–Stokes Equations with Small Initial Value in $${b^{0}_{n,\infty}(\omega)}$$ B N , ∞ 0 ( Ω )

Global well-posedness of the Navier–Stokes equations in homogeneous Besov spaces on the half-space

Global Well-Posedness and Long Time Decay of Fractional Navier-Stokes Equations in Fourier-Besov Spaces

Well-Posedness for the 3d Incompressible Nematic Liquid Crystal System in the Critical L-P Framework

ILL-POSEDNESS FOR THE COMPRESSIBLE NAVIER–STOKES EQUATIONS WITH THE VELOCITY IN FRAMEWORK

A Positivity-Preserving, Energy Stable and Convergent Numerical Scheme for the Poisson-Nernst-Planck System

Global Well-Posedness to the 3-D Incompressible Inhomogeneous Navier-Stokes Equations with a Class of Large Velocity