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

Jie Ding,Shenggao Zhou

2023-09-07

Abstract:First-order energy dissipative schemes in time are available in literature for the Poisson-Nernst-Planck (PNP) equations, but second-order ones are still in lack. This work proposes novel second-order discretization in time and finite volume discretization in space for modified PNP equations that incorporate effects arising from ionic steric interactions and dielectric inhomogeneity. A multislope method on unstructured meshes is proposed to reconstruct positive, accurate approximations of mobilities on faces of control volumes. Numerical analysis proves that the proposed numerical schemes are able to unconditionally ensure the existence of positive numerical solutions, original energy dissipation, mass conservation, and preservation of steady states at discrete level. Extensive numerical simulations are conducted to demonstrate numerical accuracy and performance in preserving properties of physical significance. Applications in ion permeation through a 3D nanopore show that the modified PNP model, equipped with the proposed schemes, has promising applications in the investigation of ion selectivity and rectification. The proposed second-order discretization can be extended to design temporal second-order schemes with original energy dissipation for a type of gradient flow problems with entropy.

Numerical Analysis

What problem does this paper attempt to address?

The problem this paper attempts to address is the development of a new second-order time discretization method for solving the modified Poisson-Nernst-Planck (PNP) equations. Specifically, the method aims to:

1. **Achieve second-order accuracy in time**: Most existing energy-dissipative schemes can only achieve first-order time accuracy, while this study proposes a new second-order time discretization method.

2. **Ensure the positivity of the numerical solution**: Guarantee the positivity of the numerical solution without adding nonlinear artificial regularization terms, which is crucial for simulating physical quantities such as ion concentration.

3. **Maintain energy dissipation properties**: Ensure that the numerical scheme unconditionally maintains the original energy dissipation properties at the discrete level, which is an important characteristic of many physical systems.

4. **Preserve mass conservation and steady-state**: Ensure that the numerical scheme unconditionally maintains mass conservation and steady-state properties at the discrete level.

5. **Handle complex geometries**: By using a multi-slope method on unstructured grids, reconstruct the positivity and accurate mobility approximation on the control volume surfaces, thereby better handling complex geometries.

The achievement of these goals will contribute to more accurate simulations of ion transport processes, with significant application prospects in the study of ion selectivity and rectification phenomena in nanopores.