[論文レビュー] Study of a chemo-repulsion model with quadratic production. Part I: Analysis of the continuous problem and time-discrete numerical schemes

Francisco Guillén‐González, María Ángeles Rodríguez‐Bellido+1 more|arXiv (Cornell University)|Mar 6, 2018

Mathematical Biology Tumor Growth参考文献 19被引用数 6

ひとこと要約

本稿は有界領域における2次生産を伴う化学的反発モデルを研究し、1次元および2次元で弱解の全域的存在および一意性を証明する。2つの無条件エネルギー安定な1次時間スキーム（後退オイラー法および線形化された結合スキーム）を導入し、可解性、収束性、誤差推定を示す。線形化スキームは正性保存を欠くが、収束性を示す。

ABSTRACT

We consider a chemo-repulsion model with quadratic production in a bounded domain. Firstly, we obtain global in time weak solutions, and give a regularity criterion (which is satisfied for $1D$ and $2D$ domains) to deduce uniqueness and global regularity. After, we study two cell-conservative and unconditionally energy-stable first-order time schemes: a (nonlinear and positive) Backward Euler scheme and a linearized coupled version, proving solvability, convergence towards weak solutions and error estimates. In particular, the linear scheme does not preserve positivity and the uniqueness of the nonlinear scheme is proved assuming small time step with respect to a strong norm of the discrete solution. This hypothesis is reduced to small time step in $nD$ domains ($n\le 2$) where global in time strong estimates are proved. Finally, we show the behavior of the schemes through some numerical simulations.

研究の動機と目的

• 有界領域における2次生産を伴う化学的反発モデルの全域的弱解の存在を確立すること。
• 1次元および2次元領域において、正則性基準の下で解の一意性および全域的正則性を証明すること。
• 非線形後退オイラー法および線形化された結合スキームの2つの無条件エネルギー安定な時間離散化スキームを設計・分析すること。
• 両スキームの可解性、弱解への収束性、および誤差推定を証明すること。
• 数値シミュレーションを通じてスキームの挙動を調査し、特にエネルギー安定性および収束性の特性に注目すること。

提案手法

• 2次生産 $ f(u) = u^2 $ を伴う化学的反発系を定式化し、系を分離するために補助変数 $ \boldsymbol{\sigma} = \nabla v $ を導入する。
• $ u $, $ \boldsymbol{\sigma} $, および $ v $ を含む再定式化された系 (3) を導出する。境界条件はノイマン型であり、初期データを含む。
• 時間離散化のための非線形後退オイラースキーム (104) を提案し、細胞保存則および無条件エネルギー安定性を保証する。
• 計算効率の高い代替手段として線形化された結合スキームを導入するが、正性保存は保証されない。
• Leray-Schauderの固定点定理を用いて、コンパクト性および連続性の議論により、時間離散スキームの解の存在を証明する。
• エネルギー推定および $ L^4 $-ノルムの有界性を用いて収束性および誤差推定を導出し、非線形スキームの一意性は強いノルムにおける小さな時間ステップ条件を要する。

実験結果

リサーチクエスチョン

• RQ1 2次生産を伴う化学的反発モデルは、有界領域で全域的弱解を有するか？
• RQ2 弱解の一意性および全域的正則性は、特に1次元および2次元でどのような条件下で成立するか？
• RQ3 本モデルに対して、無条件エネルギー安定な1次時間離散スキームを構築できるか？
• RQ4 提案されたスキーム、特に線形化されたバージョンの収束性および誤差推定は何か？
• RQ5 数値的に、エネルギー安定性、エネルギー減衰、解の分布の観点から、スキームはどのように振る舞うか？

主な発見

• 2次生産 $ f(u) = u^2 $ を伴う化学的反発モデルは、有界領域で全域的弱解を有する。
• 正則性基準の下で、1次元および2次元領域において解の一意性および全域的正則性が確立される。
• 非線形後退オイラースキームは、無条件エネルギー安定であり、細胞保存則を満たし、弱解へ収束し、誤差推定が得られる。
• 線形化された結合スキームは、無条件エネルギー安定であり、弱解へ収束するが、正性保存を欠く。
• 非線形スキームの一意性は、強いノルムにおける小さな時間ステップ条件のもとで証明され、$ n \leq 2 $ 次元領域では、全域的強い推定が得られるため、時間ステップの小さな条件に簡略化される。
• 数値的シミュレーションにより、両スキームのエネルギー安定性および収束性が確認され、理論的分析の妥当性が裏付けられる。