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.
Yuan et al. (Sat,) studied this question.