2025/02/06 by Waixiang Cao, Cao, Waixiang, Yuzhe Qin +3 · 1 citation
Mathematics · Physics and Astronomy · #Advanced Mathematical Physics Problems #FOS: Mathematics #Gas Dynamics and Kinetic Theory #Numerical Analysis (math.NA) #Optical properties and cooling technologies in crystalline materials
paper · pdf · doi:10.48550/arxiv.2502.03892
openalex publication_date 2025/02/06 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
In this paper, we introduce and analyze a class of numerical schemes that demonstrate remarkable superiority in terms of efficiency, the preservation of positivity, energy stability, and high-order precision to solve the time-dependent Poisson-Nernst-Planck (PNP) system, which is as a highly versatile and sophisticated model and accommodates a plenitude of applications in the emulation of the translocation of charged particles across a multifarious expanse of physical and biological systems. The numerical schemes presented here are based on the energy variational formulation. It allows the PNP system to be reformulated as a non-constant mobility H-1 gradient flow, incorporating singular logarithmic energy potentials. To achieve a fully discrete numerical scheme, we employ a combination of first/second-order semi-implicit time discretization methods, coupled with either the k-th order direct discontinuous Galerkin (DDG) method or the finite element (FE) method for spatial discretization. The schemes are verified to possess positivity preservation and energy stability. Optimal error estimates and particular superconvergence results for the fully-discrete numerical solution are established. Numerical experiments are provided to showcase the accuracy, efficiency, and robustness of the proposed schemes.