[論文レビュー] Structure preserving schemes for nonlinear Fokker-Planck equations with anisotropic diffusion

Nadia Loy, Mattia Zanella|arXiv (Cornell University)|May 8, 2019

Advanced Thermodynamics and Statistical Mechanics被引用数 1

ひとこと要約

本稿では、非線形Fokker-Planck方程式に特徴的な異方性拡散を伴う高次構造保存型数値スキームを導入し、解の非負性、エントロピー散逸、および正確な定常状態近似を維持する。スキームは一時的状態では少なくとも2次精度を達成し、長時間においては高次精度を示し、集団行動や生命科学分野の異常拡散の安定なシミュレーションを可能にする。

ABSTRACT

In this work we propose novel numerical schemes for nonlinear Fokker-Planck-type equations with anisotropic diffusion matrix that preserve fundamental structural properties like non negativity of the solution, entropy dissipation and which guarantees an arbitrarily accurate approximation of the steady state of the problem. All the methods presented are at least second order accurate in the transient regimes and high order for large times. Applications of the schemes to models for collective phenomena and life sciences are considered, in these examples anomalous diffusion is often observed and must be taken into account in realistic models.

研究の動機と目的

• 非線形Fokker-Planck方程式に特徴的な異方性拡散を伴う数値スキームの開発を目的とし、解の非負性、エントロピー散逸、および正確な定常状態近似といった重要な構造的性質を保存すること。
• 集団現象や生命科学分野における異常拡散のモデル化という課題に取り組み、標準的なスキームが物理的整合性を保てない場合があることを補う。
• 一時的状態では少なくとも2次精度を達成し、長時間においては高次精度を達成することで、長時間にわたるシミュレーションの忠実性を向上させること。

提案手法

• 非負性とエントロピー構造の保存を目的として、半implicit時刻離散化と有限体積法による空間離散化を組み合わせた手法を採用する。
• 解変数に非線形変換を施し、非負性の保証と異方性拡散が存在する状況でのスキーム安定化を図る。
• フラックスの高次再構成技術を用いることで、時間方向に2次精度、空間方向に長時間において高次精度を達成する。
• 連続方程式の構造から導かれる離散的エントロピー不等式を通じて、エントロピー安定性を強制的に確保する。
• 構造的設計により、複雑な異方性拡散テンソルに対しても、任意の精度で正しい定常状態に収束するようにスキームを設計する。
• 集団ダイナミクスや生物学的系に現れる異常拡散を示すモデルへの応用を通じて、手法の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 非線形Fokker-Planck方程式に特徴的な異方性拡散を伴う数値スキームは、どのように解の非負性を保ちながら高次精度を達成できるか？
• RQ2 離散的構造保存技術として、連続問題と同一の熱力学的挙動を再現するエントロピー散逸を保証する手法は何か？
• RQ3 複雑な異方性拡散が存在する状況においても、どのようにして定常状態を任意の精度で近似できるか？
• RQ4 スキームは一時的状態および長時間状態の両方において、どの程度高次精度を維持できるか？
• RQ5 提案手法は、集団行動や生命科学分野のモデルにおける異常拡散を効果的にシミュレートできるか？

主な発見

• 提案されたスキームは、任意の異方性拡散行列に対しても解の非負性を保証し、物理的整合性を確保する。
• スキームは離散的エントロピー散逸法則を満たし、連続系の熱力学的挙動を模倣する。
• 一時的状態では2次精度が達成され、長時間においては高次精度が観察され、長期的シミュレーションの信頼性が向上する。
• 構造的設計のおかげで、時間刻み幅に依存せず、任意の精度で正しい定常状態に収束する。
• 集団運動や生命科学分野のモデルへの応用を通じて、異常拡散現象を安定的かつ高精度に捉えることの有効性が示された。