An Energy-Stable S-SAV Finite Element Method for the Generalized Poisson-Nernst-Planck Equation

Authors:
Open Access
Publisher website
Google Scholar
Abstract
Designing structure-preserving numerical schemes for the generalized Poisson-Nernst-Planck (PNP) system is challenging due to its inherent strong nonlinearity and coupling. In this paper, we propose a class of efficient, unconditional energy-stable schemes based on the Stabilized Scalar Auxiliary Variable (S-SAV) framework combined with the finite element method. We construct both first-order (BE-S-SAV) and second-order (BDF2-S-SAV) fully discrete schemes. A distinguishing feature of our approach is the use of a linear decomposition strategy, which decouples the complex nonlinear system into a sequence of linear, constant-coefficient elliptic equations at each time step. This significantly reduces computational complexity by avoiding expensive nonlinear iterations. We provide rigorous theoretical proofs demonstrating that the proposed schemes are unconditionally energy stable and strictly preserve mass conservation. Numerical experiments satisfy the theoretical analysis, confirming optimal convergence rates and demonstrating robust preservation of mass conservation and modified energy stability in the tested regimes.
Funding Information
  • Natural Science Foundation of Xinjiang Uygur Autonomous Region (2025D01B197)
  • Research Foundation of China University of Petroleum–Beijing at Karamay (XQZX20230030)
  • Xinjiang University Doctoral Science and Technology Innovation Program (XJU2023BS026)
  • Tianyuan Fund for Mathematics of the National Natural Science Foundation of China (12326102)
  • 2025 National Undergraduate Innovation and Entrepreneurship Training Program of China (202519414009)

This publication has 40 references indexed in Scilit: