[論文レビュー] Unconditional positivity-preserving and energy stable schemes for a reduced Poisson-Nernst-Planck system

Hailiang Liu, Wumaier Maimaitiyiming|arXiv (Cornell University)|Sep 28, 2019

Quantum, superfluid, helium dynamics参考文献 30被引用数 6

ひとこと要約

本稿では、イオンチャネルをモデル化する低次元化されたポアソン-ナーンスト-プランク系に対する、2次精度で、正性を保証し、エネルギー安定性を満たす有限体積スキームを提示する。非対数形式のランダウ形式を用いた半implicit時間離散化により、いかなる時間刻みサイズでも、ポアソンソルバの選択に依存せず、イオン濃度の正性が保証される。また、ゼロフラックス境界条件下で質量保存則を満たし、自由エネルギーを散逸させる。

ABSTRACT

The Poisson-Nernst-Planck (PNP) system is a widely accepted model for simulation of ionic channels. In this paper, we design, analyze, and numerically validate a second order unconditional positivity-preserving scheme for solving a reduced PNP system, which can well approximate the three dimensional ion channel problem. Positivity of numerical solutions is proven to hold true independent of the size of time steps and the choice of the Poisson solver. The scheme is easy to implement without resorting to any iteration method. Several numerical examples further confirm the positivity-preserving property, and demonstrate the accuracy, efficiency, and robustness of the proposed scheme, as well as the fast approach to steady states.

研究の動機と目的

• 低次元ポアソン-ナーンスト-プランク系におけるイオン濃度の物理的正性を保証する数値スキームの開発を目的とする。
• 時間刻みサイズやポアソンソルバの選択にかかわらず、正性を非条件的に保証することを目的とする。
• ゼロフラックス境界条件下で質量保存則と自由エネルギー散逸を維持することを目的とする。
• イオン輸送シミュレーションに適した、単純で反復処理を必要とせず、容易に実装可能な数値手法を提供することを目的とする。
• 多様な初期条件および境界条件の下で、スキームの精度、頑健性、および定常状態への高速収束性を検証することを目的とする。

提案手法

• 正性を保証する離散化を可能にするために、低次元PNP系を非対数形式のランダウ形式に変換する。
• 2次精度を達成する半implicit有限体積法を時間方向に適用する。
• 非一様メッシュ上でのセル中心型有限体積法を空間方向に用いる。
• イオン濃度と電気的ポテンシャルの式を結合するために、非ステガラ式グリッド配置を採用する。
• 繰り返し補正を回避するため、フラックス再構築に基づく正性を保証する更新戦略を実装する。
• ディリクレ型およびロビン型境界条件（ゼロフラックスおよびポテンシャル反射境界条件を含む）とスキームを統合する。

実験結果

リサーチクエスチョン

• RQ1 低次元PNP系において、いかなる時間刻みサイズやポアソンソルバの選択に対しても、正性を保証する2次精度の数値スキームを構築可能か？
• RQ2 ゼロフラックス境界条件下で、質量保存則と自由エネルギー散逸が維持されるか？
• RQ3 初期データ、時間刻みサイズ、ポアソンソルバの選択に依存しないスキームの頑健性はいかがなものか？
• RQ4 異なるイオンチャネル構造における定常状態への収束挙動はどのように変化するか？
• RQ5 固定電荷密度および初期濃度プロファイルが、定常状態におけるイオン分布に与える影響は何か？

主な発見

• 提案スキームは、定理3.1で厳密に証明されたように、いかなる時間刻みサイズやポアソンソルバの選択に対しても、イオン濃度の正性を非条件的に保証する。
• 空間および時間方向で2次精度を達成し、数値的収束率は理論的期待値と一致する。
• 例4.4（初期データ式(4.8)）において、τ=10⁻³を用いて、2.7410時間単位（2741反復）で定常状態に到達し、残差誤差||ψⁿ−ψⁿ⁻¹||∞ ≈ 10⁻⁷となる。
• ランダム初期データを用いた場合、2.0440時間単位（2044反復）で定常状態に到達し、初期条件に依存しない頑健性が確認された。
• 例4.4の両方のデータセットにおいて同一のプロファイルが得られたことから、定常状態におけるイオン分布は初期データに依存しないことが確認された。
• 例4.5（ゼロフラックスおよびロビン境界条件）において、合計質量が保存され、時間経過に伴い自由エネルギーが単調に散逸することが確認され、定理3.2の妥当性が裏付けられた。