[論文レビュー] On the hydrodynamical limit for a one dimensional kinetic model of cell aggregation by chemotaxis

François James, Nicolas Vauchelet|arXiv (Cornell University)|Oct 18, 2010

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

ひとこと要約

本稿は、細胞の走化性（チラクシス）がランダム運動を上回る1次元の運動論的モデルの流体力学的極限を確立し、保存則と楕円型方程式の結合系への収束を示している。巨視的速度が不連続で解が測度値をとるにもかかわらず、双対解を用いて収束を証明しており、数値シミュレーションにより、Dirac型構造への凝集と、キナーゲンチエントの勾配における特異性が確認されている。

ABSTRACT

The hydrodynamic limit of a one dimensional kinetic model describing chemotaxis is investigated. The limit system is a conservation law coupled to an elliptic problem for which the macroscopic velocity is possibly discontinuous. Therefore, we need to work with measure-valued densities. After recalling a blow-up result in finite time of regular solutions for the hydrodynamic model, we establish a convergence result of the solutions of the kinetic model towards solutions of a problem limit defined thanks to the flux. Numerical simulations illustrate this convergence result.

研究の動機と目的

• 走化性が無向運動を上回る状況下で、細胞凝集体の1次元運動論的モデルの流体力学的極限を導出すること。
• 密度が測度値をとり、キナーゲンチエント勾配が不連続である場合に、極限系における巨視的速度を定義する数学的課題に取り組むこと。
• 運動論的モデルの解が、キナーゲンチエントに対してスカラー保存則が結合された楕円型方程式を含む極限系に収束することを確立すること。
• キナーゲンチエント勾配における特異性の生成と、極限におけるDirac型細胞凝集体の出現を分析すること。

提案手法

• 小さなパラメータ ε = v₀/(φ₀x₀) を用いて運動論的モデルをスケーリングし、細胞の速度と回転率の比を表す。
• ε → 0 の極限において形式的に流体力学的極限を導出し、細胞密度の保存則とキナーゲンチエントの楕円型方程式の結合系が得られることを示す。
• 速度に依存する回転核から生じる不連続なフラックス関数 a(∂ₓS) を持つ保存則を扱うために、双対解を用いる。
• 双対解の存在を保証するため、フラックス関数 a(∂ₓS) に片側リプシッツ（OSL）条件を課す。
• 実際のパrameters（例：E. coli の速度 c = 20×10⁻⁶ m/s、初期状態はガウス分布）を用いた数値シミュレーションにより、収束性と凝集ダイナミクスを検証する。
• ε を変化させた際の細胞密度 ρ とキナーゲンチエント勾配 ∂ₓS の時間発展を分析し、ε → 0 に近づくにつれて凝集が速くなり、勾配に特異性が現れることを観察する。

実験結果

リサーチクエスチョン

• RQ1 走化性が無向運動を上回る1次元運動論的走化性モデルに対して、流体力学的極限を厳密に確立できるか？
• RQ2 解が測度値をとり、キナーゲンチエント勾配が不連続である場合に、極限系における巨視的速度をどのように定義できるか？
• RQ3 極限系は有限時間での爆発またはDirac型構造への凝集を示すか？ そして、その現象は数値的にどのように捉えられるか？
• RQ4 双対解は、不連続なフラックス関数 a(∂ₓS) を持つ極限保存則に対して適切な枠組みを提供できるか？
• RQ5 極限系におけるフラックスと、運動論的モデルから導出される巨視的速度との関係は何か？

主な発見

• 巨視的速度が不連続であっても、運動論的モデルの流体力学的極限は、キナーゲンチエントに対してスカラー保存則が結合された楕円型方程式の系に収束する。
• 極限系は測度値をとる密度と、キナーゲンチエント勾配におけるBV正則性を特徴とし、細胞凝集体付近で ∂ₓS に特異性が生じる。
• 数値シミュレーションにより、ε → 0 に近づくにつれて、細胞密度がドメイン中央にDirac型構造へ凝集することが確認され、ε が小さいほど凝集が速くなる。
• フラックス関数 a(∂ₓS) は片側リプシッツ（OSL）条件を満たしており、双対解の存在を保証するが、一意性は未解決のまま。
• キナーゲンチエント ∂ₓS はドメイン中央に急激なピークを示し、ε を小さくするにつれて特異性が形成されることを示している。
• ε = 10⁻⁴ の場合、明確な特異的勾配と急速な凝集が観察され、理論的収束が特異極限に一致することを支持する。