[论文解读] A positivity-preserving, energy stable and convergent numerical scheme for the Poisson-Nernst-Planck system

Chun Liu, Cheng Wang+3 more|arXiv (Cornell University)|Sep 17, 2020

Gas Dynamics and Kinetic Theory参考文献 35被引用 8

一句话总结

本文提出了一种泊松-能斯特-普朗克（PNP）系统的有限差分格式，该格式同时保证了唯一可解性、正性保持、能量稳定性和最优收敛速率。通过利用能量变分法（EnVarA），该格式将迁移率显式处理，对数/静电项则隐式处理，借助凸性及奇异对数势能，防止浓度爆破，实现了在非恒定迁移率下首次理论上的收敛性与稳定性。

ABSTRACT

In this paper we propose and analyze a finite difference numerical scheme for the Poisson-Nernst-Planck equation (PNP) system. To understand the energy structure of the PNP model, we make use of the Energetic Variational Approach (EnVarA), so that the PNP system could be reformulated as a non-constant mobility $H^{-1}$ gradient flow, with singular logarithmic energy potentials involved. To ensure the unique solvability and energy stability, the mobility function is explicitly treated, while both the logarithmic and the electric potential diffusion terms are treated implicitly, due to the convex nature of these two energy functional parts. The positivity-preserving property for both concentrations, $n$ and $p$, is established at a theoretical level. This is based on the subtle fact that the singular nature of the logarithmic term around the value of $0$ prevents the numerical solution reaching the singular value, so that the numerical scheme is always well-defined. In addition, an optimal rate convergence analysis is provided in this work, in which many highly non-standard estimates have to be involved, due to the nonlinear parabolic coefficients. The higher order asymptotic expansion (up to third order temporal accuracy and fourth order spatial accuracy), the rough error estimate (to establish the $\ell^\infty$ bound for $n$ and $p$), and the refined error estimate have to be carried out to accomplish such a convergence result. In our knowledge, this work will be the first to combine the following three theoretical properties for a numerical scheme for the PNP system: (i) unique solvability and positivity, (ii) energy stability, and (iii) optimal rate convergence. A few numerical results are also presented in this article, which demonstrates the robustness of the proposed numerical scheme.

研究动机与目标

• 为PNP系统开发一种数值格式，以保持正性与能量耗散等物理特性。
• 在非恒定迁移率和奇异对数势能条件下，建立有限差分格式的严格收敛性分析。
• 在一个统一的数值框架中整合三个关键理论性质——唯一可解性/正性、能量稳定性和最优收敛速率。
• 克服能量泛函中非线性抛物型系数与奇异对数项带来的挑战。

提出的方法

• 利用能量变分法（EnVarA）将PNP系统重表述为具有非恒定迁移率的H⁻¹梯度流形式。
• 显式处理迁移率函数，以确保离散系统的唯一可解性。
• 由于能量泛函具有凸性，对对数项和电势扩散项采用隐式时间格式。
• 通过高阶渐近展开（时间精度达三阶，空间精度达四阶）推导浓度的粗略ℓ∞范数界。
• 采用精细化误差估计，即使在非线性抛物型系数下，仍实现最优收敛速率。
• 利用对数势能在零附近奇异的特性，理论上保证数值解的正性。

实验结果

研究问题

• RQ1 能否为PNP系统构造一种有限差分格式，以同时保证唯一可解性与浓度的正性？
• RQ2 如何严格证明具有非恒定迁移率和奇异对数势能的格式的能量稳定性？
• RQ3 在如此复杂的非线性条件下，PNP系统数值格式的最优收敛速率是多少？
• RQ4 能否构建一个统一的理论框架，同时保证正性、能量稳定性和最优收敛速率？

主要发现

• 所提出的格式在离散层面实现了唯一可解性，并通过奇异对数势能防止数值爆破，从而保持了离子浓度n和p的正性。
• 通过能量泛函的凸性及对数项与静电项的隐式处理，证明了能量稳定性。
• 建立了最优收敛速率，其中高阶渐近展开达到三阶时间精度与四阶空间精度。
• 通过粗略误差估计推导出浓度的ℓ∞范数界，这对于处理非线性抛物型系数至关重要。
• 精细化误差估计证实了最优收敛速率，克服了非恒定迁移率与奇异势能带来的挑战。
• 数值结果验证了质量守恒、能量耗散与正性保持，与理论分析一致。