An unconditionally energy stable linear scheme for Poisson–Nernst–Planck equations

Tian Qiao, Zhonghua Qiao, Shuyu Sun*, Shenggao Zhou

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

Abstract

This paper proposes a linear, unconditionally energy-stable scheme for the Poisson–Nernst–Planck (PNP) equations. Based on a gradient-flow formulation of the PNP equations, the energy factorization approach is applied to linearize the logarithm function at the previous time step, resulting in a linear semi-implicit scheme. Numerical analysis is conducted to illustrate that the proposed fully discrete scheme has desired properties at a discrete level, such as unconditional unique solvability, mass conservation, and energy dissipation. Numerical simulations verify that the proposed scheme, as expected, is first-order accurate in time and second-order accurate in space. Further numerical tests confirm that the proposed scheme can indeed preserve the desired properties. Applications of our numerical scheme to the simulations of electrolyte solutions demonstrate that, as a linear energy stable scheme of efficiency, it will be promising in simulating complicated transport phenomena of charged systems.

Original languageEnglish (US)
Article number115759
JournalJournal of Computational and Applied Mathematics
Volume443
DOIs
StatePublished - Jun 2024

Keywords

  • Electric double layer
  • Energy stability
  • Mass conservation
  • Poisson–Nernst–Planck equations

ASJC Scopus subject areas

  • Computational Mathematics
  • Applied Mathematics

Fingerprint

Dive into the research topics of 'An unconditionally energy stable linear scheme for Poisson–Nernst–Planck equations'. Together they form a unique fingerprint.

Cite this