[論文レビュー] A fully discrete variational scheme for solving nonlinear Fokker-Planck equations in higher space dimensions

Oliver Junge, Daniel Matthes+1 more|arXiv (Cornell University)|Sep 25, 2015

Geometric Analysis and Curvature Flows参考文献 11被引用数 11

ひとこと要約

本稿は、$L^2$-Wasserstein距離における勾配流構造と流れ写像による押し出し表現を活用することで、高次元空間における非線形Fokker-Planck方程式を完全離散的・変分的ラグランジュスキームで解く手法を提案する。このスキームはエントロピーを減少させ、質量を保存し、正則性仮定のもとで収束性を示す。2次元での数値的検証も行われている。

ABSTRACT

We introduce a novel spatio-temporal discretization for nonlinear Fokker-Planck equations on the multi-dimensional unit cube. This discretization is based on two structural properties of these equations: the first is the representation as a gradient flow of an entropy functional in the $L^2$-Wasserstein metric, the second is the Lagrangian nature, meaning that solutions can be written as the push forward transformation of the initial density under suitable flow maps. The resulting numerical scheme is entropy diminishing and mass conserving. Further, the scheme is weakly stable, which allows us to prove convergence under certain regularity assumptions. Finally, we present results from numerical experiments in space dimension $d=2$.

研究の動機と目的

• 高次元空間における非線形Fokker-Planck方程式の構造保存型数値スキームの開発。
• $L^2$-Wasserstein距離における勾配流構造と微分同相写像による流れを用いたラグランジュ的定式化の活用。
• 完全離散的枠組み内で、エントロピーの減少と質量保存といった重要な物理的性質の維持。
• 圧力関数、ポテンシャル、初期データの正則性仮定のもとで、スキームの弱安定性と収束性の確立。
• 微分同相写像の有限次元多様体上での最小化移動法を用いた2次元空間におけるスキームの数値的検証。

提案手法

• Fokker-Planck方程式を、$L^2$-Wasserstein距離におけるエントロピー汎関数の勾配流として定式化する。
• 解を時間発展する微分同相写像 $\mathfrak{T}_t$ による基準密度の押し出しとして表現する。
• 微分同相写像 $\mathfrak{T}_\Delta^n$ の時間発展を、段階的微分同相写像 $\mathfrak{t}_\Delta^m$ の合成として離散化する。
• 微分同相写像多様体上のエネルギー汎関数の勾配流に最小化移動法を適用し、時間離散化された作用汎関数を最小化する。
• 段階的微分同相写像に有限次元的アンザッツを用いることで、計算の実行可能性を確保するとともに、幾何的構造を保存する。
• 測度の変換とヤコビアン変換を介して、$L^2(\mathrm{d}x)$-勾配を微分同相写像多様体上の $L^2(\bar{u})$-勾配に変換する。

実験結果

リサーチクエスチョン

• RQ1 高次元空間における非線形Fokker-Planck方程式に対して、エントロピーと質量を保存する完全離散的変分スキームを構築可能か？
• RQ2 Fokker-Planck方程式のラグランジュ的構造を、空間的・時間的両方にわたって離散化しつつ、幾何的一致性を保てるか？
• RQ3 圧力関数および外部ポテンシャルの正則性仮定のもとで、このようなスキームの収束性はどのように振る舞うか？
• RQ4 微分同相写像多様体上での最小化移動法は、安定性と構造保存性をどのように保証するか？
• RQ5 2次元空間におけるスキームの数値的性能は、特にエントロピーの減少と質量保存の観点からどのように評価できるか？

主な発見

• 提案スキームはエントロピーを減少させ、質量を保存し、連続問題が持つ重要な物理的性質を継承する。
• スキームは弱安定性を示し、$\mathrm{P}$, $V$, および $u^0$ の正則性仮定のもとで収束性の証明が可能である。
• 離散的時間発展は、微分同相写像の有限次元多様体上での最小化移動法により構築され、計算の実行可能性が保証される。
• 微分同相写像多様体上のエネルギー汎関数の勾配は、押し出し測度とヤコビアン行列式を用いて明示的に特徴付けられる。
• 2次元空間における数値実験により、スキームが時間経過に伴うエントロピーの減少と質量保存を良好に維持していることが示された。
• 理論的解析により、$L^2(\mathrm{d}x)$-勾配が測度変換と押し出し構造を介して、微分同相写像多様体上の $L^2(\bar{u})$-勾配に対応することが確認された。