Abstract:Ion transport, often described by the Poisson--Nernst--Planck (PNP) equations, is ubiquitous in electrochemical devices and many biological processes of significance. In this work, we develop conservative, positivity-preserving, energy dissipating, and implicit finite difference schemes for solving the multi-dimensional PNP equations with multiple ionic species. A central-differencing discretization based on harmonic-mean approximations is employed for the Nernst--Planck (NP) equations. The backward Euler discretization in time is employed to derive a fully implicit nonlinear system, which is efficiently solved by a newly proposed Newton's method. The improved computational efficiency of the Newton's method originates from the usage of the electrostatic potential as the iteration variable, rather than the unknowns of the nonlinear system that involves both the potential and concentration of multiple ionic species. Numerical analysis proves that the numerical schemes respect three desired analytical properties (conservation, positivity preserving, and energy dissipation) fully discretely. Based on advantages brought by the harmonic-mean approximations, we are able to establish estimate on the upper bound of condition numbers of coefficient matrices in linear systems that are solved iteratively. The solvability and stability of the linearized problem in the Newton's method are rigorously established as well. Numerical tests are performed to confirm the anticipated numerical accuracy, computational efficiency, and structure-preserving properties of the developed schemes. Adaptive time stepping is implemented for further efficiency improvement. Finally, the proposed numerical approaches are applied to characterize ion transport subject to a sinusoidal applied potential.

NUMERICAL ANALYSIS OF NON-FICKIAN DIFFUSION PROBLEMS IN A POTENTIAL FIELD

Study of Non‐fickian Diffusion Problems with the Potential Field in the Cylindrical Co‐ordinate System

Effect of the Potential Field on Non-Fickian Diffusion Problems in a Sphere

Analysis of Non-Fickian Diffusion Problems in A Composite Medium

Study of Non-Linear Mesoscopic Diffusion with Potential Field Effect

Numerical Analysis for the Hyperbolic Heat Conduction Problem under a Pulsed Surface Disturbance

Quasi One-dimensional Non-passing Self-diffusion

A new approach in an analytical method for diffusion dynamics for the presence of delocalized sink in a potential well: Application to different potential curves

A provably stable numerical method for the anisotropic diffusion equation in confined magnetic fields

Numerical Solution of the Nonlocal Diffusion Equation on the Real Line.

Numerical analysis of a class of penalty discontinuous Galerkin methods for nonlocal diffusion problems

Structure-Preserving and Efficient Numerical Methods for Ion Transport

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

Numerical analysis of singularly perturbed parabolic reaction diffusion differential difference equations

An Asymptotic Preserving method for strongly anisotropic diffusion equations based on field line integration

Numerical Methods for Linear Diffusion Equations in the Presence of an Interface

Convergence of numerical schemes for the solution of partial integro-differential equations used in heat transfer

Intelligence Algorithm Research Of Convection Diffusion Equation And Numerical Simulation

Analyses of the contour integral method for time fractional subdiffusion-normal transport equation

Stability analysis and numerical results for some schemes discretising 2D nonconstant coefficient advection–diffusion equations