TY - JOUR
T1 - An iteration solver for the Poisson–Nernst–Planck system and its convergence analysis
AU - Liu, Chun
AU - Wang, Cheng
AU - Wise, Steven M.
AU - Yue, Xingye
AU - Zhou, Shenggao
N1 - Publisher Copyright:
© 2021 Elsevier B.V.
PY - 2022/5/1
Y1 - 2022/5/1
N2 - In this paper, we provide a theoretical analysis for an iteration solver to implement a finite difference numerical scheme for the Poisson-Nernst–Planck (PNP) system, based on the Energetic Variational Approach (EnVarA), in which a non-constant mobility H−1 gradient flow is formulated. In particular, the nonlinear and singular nature of the logarithmic energy potentials has always been the essential difficulty. In the numerical design, the mobility function is explicitly updated, for the sake of unique solvability analysis. The logarithmic and the electric potential diffusion terms, which come from the gradient of convex energy functional parts, are implicitly computed. The positivity-preserving property for all the concentrations, an unconditional energy stability, and the optimal rate error estimate have been established in a recent work. A modified Newton iteration for the nonlinear and logarithmic part, combined with a linear iteration for the electric potential part, is proposed to implement the given numerical scheme, in which a non-constant linear elliptic equation needs to be solved at each iteration stage. A theoretical analysis is presented in this article, and a linear convergence is proved for such an iteration, with an asymptotic error constant in the same order of the time step size. A numerical test is also presented in this article, which demonstrates the linear convergence rate of the proposed iteration solver.
AB - In this paper, we provide a theoretical analysis for an iteration solver to implement a finite difference numerical scheme for the Poisson-Nernst–Planck (PNP) system, based on the Energetic Variational Approach (EnVarA), in which a non-constant mobility H−1 gradient flow is formulated. In particular, the nonlinear and singular nature of the logarithmic energy potentials has always been the essential difficulty. In the numerical design, the mobility function is explicitly updated, for the sake of unique solvability analysis. The logarithmic and the electric potential diffusion terms, which come from the gradient of convex energy functional parts, are implicitly computed. The positivity-preserving property for all the concentrations, an unconditional energy stability, and the optimal rate error estimate have been established in a recent work. A modified Newton iteration for the nonlinear and logarithmic part, combined with a linear iteration for the electric potential part, is proposed to implement the given numerical scheme, in which a non-constant linear elliptic equation needs to be solved at each iteration stage. A theoretical analysis is presented in this article, and a linear convergence is proved for such an iteration, with an asymptotic error constant in the same order of the time step size. A numerical test is also presented in this article, which demonstrates the linear convergence rate of the proposed iteration solver.
UR - http://www.scopus.com/inward/record.url?scp=85121856690&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85121856690&partnerID=8YFLogxK
U2 - 10.1016/j.cam.2021.114017
DO - 10.1016/j.cam.2021.114017
M3 - Article
AN - SCOPUS:85121856690
SN - 0377-0427
VL - 406
JO - Journal of Computational and Applied Mathematics
JF - Journal of Computational and Applied Mathematics
M1 - 114017
ER -