[論文レビュー] 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

ひとこと要約

本稿では、ポisson-ナーンスト・プランク（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系の数値スキームを構築し、正値性やエネルギー散逸といった物理的性質を保存すること。
• 非定数移動度および特異的対数ポテンシャル下での有限差分スキームに対して、厳密な収束解析を確立すること。
• 一意的可解性／正値性、エネルギー安定性、最適収束速度という3つの重要な理論的性質を、一つの数値フレームワークで統合すること。
• エネルギー関数形に特異的対数項を含む非線形放物型係数に起因する課題を克服すること。

提案手法

• エネルギー的変分的アプローチ（EnVarA）を用いて、PNP系を非定数移動度を持つH⁻¹勾配流として再定式化する。
• 離散系の一意的可解性を保証するために、移動度関数を明示的に取り扱う。
• 対数項および電気的ポテンシャル拡散項に凸性エネルギー構造があるため、陰的時間離散化を適用する。
• 濃度の粗いℓ∞バインディングを導出するために、高次漸近展開（時間方向3次、空間方向4次精度まで）を用いる。
• 非線形放物型係数に対処するために、精密な誤差推定を用いて最適収束率を達成する。
• ゼロ付近での対数ポテンシャルの特異性を活用して、数値解の正値性を理論的に保証する。

実験結果

リサーチクエスチョン

• RQ1 PNP系に対して、一意的可解性と濃度の正値性を両立させる有限差分スキームを構築可能か？
• RQ2 非定数移動度および特異的対数ポテンシャルを有するスキームに対して、エネルギー安定性を厳密に証明できるか？
• RQ3 このような複雑な非線形性下で、PNP系の数値スキームの最適収束速度は何か？
• RQ4 正値性、エネルギー安定性、最適収束速度を同時に保証する統一的理論枠組みを構築可能か？

主な発見

• 特異的対数ポテンシャルのおかげで、数値的爆発を防げるため、提案スキームは離散的に両方のイオン濃度nとpの正値性を保証する。
• エネルギー関数形の凸性と対数項・電気的ポテンシャル項の陰的取り扱いにより、エネルギー安定性が証明される。
• 時間方向3次、空間方向4次精度までの高次漸近展開を用いて、最適収束速度が確立される。
• 非線形放物型係数に対処するため、濃度のℓ∞バインディングは粗い誤差推定を通じて導出される。
• 精密な誤差推定により、非定数移動度および特異的ポテンシャルに起因する課題を克服し、最適収束速度が確認される。
• 数値結果により、質量保存、エネルギー散逸、正値性保存が離散的に確認され、理論的解析の妥当性が検証される。