[論文レビュー] Convergence analysis of upwind type schemes for the aggregation equation with pointy potential

François Delarue, Frédéric Lagoutière+1 more|arXiv (Cornell University)|Sep 27, 2017

Mathematical Biology Tumor Growth参考文献 54被引用数 10

ひとこと要約

本稿では、鋭いポテンシャル（非滑らかで非局所的な相互作用力）を有する凝集方程式に対する上流有限体積スキームの収束解析を提示する。この力は不連続な速度場を引き起こす。Filippov流れ理論を用いて、測度値解の設定において Wasserstein 距離で 1/2 階の収束を証明し、滑らかでないポテンシャルを有する場合のこの種のスキームにおける最初の最適収束率を確立する。

ABSTRACT

A numerical analysis of upwind type schemes for the nonlinear nonlocal aggregation equation is provided. In this approach, the aggregation equation is interpreted as a conservative transport equation driven by a nonlocal nonlinear velocity field with low regularity. In particular, we allow the interacting potential to be pointy, in which case the velocity field may have discontinuities. Based on recent results of existence and uniqueness of a Filippov flow for this type of equations, we study an upwind finite volume numerical scheme and we prove that it is convergent at order 1/2 in Wasserstein distance. This result is illustrated by numerical simulations, which show the optimality of the order of convergence.

研究の動機と目的

• 非滑らかで鋭い相互作用ポテンシャルを有する凝集方程式に対する上流型有限体積スキームの収束を分析すること。
• 特異的で凸でないポテンシャルが引き起こす低正則性の速度場（不連続性を伴う）という課題に対処すること。
• 滑らかでない解が爆発するような状況を含め、測度値解フレームワークにおける数値スキームの厳密な収束速度を確立すること。
• 1次元および2次元における数値シミュレーションを通じて理論的収束次数の妥当性を検証すること。
• 従来の手法（例：Lax-Friedrichs）と比較して、本スキームのロバスト性と数値拡산の低減を示すこと。

提案手法

• 鋭いポテンシャルの勾配から導かれる非局所的かつ非滑らかな速度場によって駆動される保存則型輸送方程式として凝集方程式を定式化する。
• 不連続な速度場に対しても一意な流れを構築できるように、Filippovの微分包含理論を適用する。
• 低正則性仮定のもとで安定性と整合性を保証する上流有限体積スキームを構築し、輸送方程式を離散化する。
• 数値解と真の測度値解との間の収束を測る誤差指標として Wasserstein 距離を用いる。
• Filippov流れの安定性と離散的双対性の議論を活用し、Wasserstein 距離で 1/2 階の収束を確立する。
• 1次元および2次元での数値シミュレーションを実施し、コンパクトに台を持つ初期データを用い、径向的で鋭いポテンシャル（例：$W(x) = 1 - e^{-5|x|}$ および $W(x) = 5|x|$）を用いて理論的結果の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 相互作用ポテンシャルが鋭く、速度場が不連続である場合、凝集方程式に対する上流有限体積スキームの収束速度は何か？
• RQ2 低正則性仮定のもとで、測度値解の設定において Filippov 流れ理論を用いて解の存在と一意性を有効に正当化できるか？
• RQ3 提案された上流スキームは、非滑らかで非局所的な速度場の下でも安定性と収束性を維持できるか？
• RQ4 特に数値拡散の観点から、Lax-Friedrichs などの他の数値手法と比較して、本スキームの収束速度はどのように異なるか？
• RQ5 異なる鋭いポテンシャル（例：$1 - e^{-5|x|}$ と $5|x|$）は、凝集のダイナミクスおよび Dirac デルタの形成タイミングにどのように影響を与えるか？

主な発見

• 上流有限体積スキームは、Wasserstein 距離で真の測度値解へ 1/2 階の収束を示し、これが最適であることが証明された。
• 収束結果は、ポテンシャルが鋭い（すなわち、$ abla W$ が原点で不連続で、0 を除く点で Lipschitz 連続）という仮定のもとで成り立つ。この仮定は (A0)–(A3) を満たす。
• 1次元および2次元における数値シミュレーションにより、理論的収束次数が確認され、Lax-Friedrichs よりも少ない数値拡散を示した。
• 複数のバンプを有する初期データに対して、本スキームはクラスタの逐次的凝集が単一の Dirac デルタに合体する過程を正しく捉え、既知の爆発ダイナミクスと整合的である。
• ポテンシャルの選択は凝集ダイナミクスに顕著な影響を与える：$W(x) = 1 - e^{-5|x|}$ はより速く、コーナー付近に局在化した集中を引き起こすが、$W(x) = 5|x|$ はより均一に端部に集中し、Dirac デルタの形成が遅れる。
• 本スキームは、粒子法が実用的でない高次元でも、爆発後の挙動（Dirac マスの形成など）を効果的にシミュレートできた。