[論文レビュー] Global solutions of a Keller--Segel system with saturated logarithmic sensitivity function

Qi Wang|arXiv (Cornell University)|Dec 1, 2013

Mathematical Biology Tumor Growth参考文献 14被引用数 3

ひとこと要約

本稿では、低濃度における特異性を回避するための正の定数 c > 0 を用いた飽和対数感受性関数 φ(v) = ln(v + c) を持つ完全な放物型 Keller–Segel きのこ移動系に対して、古典解の大域的存在を確立する。重み付きエネルギー推定と最大 Lp-正則性理論を用いて、拡散比 k = d2/d1 に依存する閾値内に化学走性係数 χ が存在する場合に大域的解が存在することを証明し、以前の研究で d1 = d2 を仮定していた条件を拡張する。

ABSTRACT

We study a Keller-Segel type chemotaxis model with a modified sensitivity function in a bounded domain $\Omega\subset \mathbb{R}^N$, $N\geq2$. The global existence of classical solutions to the fully parabolic system is established provided that the ratio of the chemotactic coefficient to the motility of cells is not too large.

研究の動機と目的

• v = 0 における ln v の特異性を回避するため、飽和対数感受性関数 φ(v) = ln(v + c) を有する放物型-放物型 Keller–Segel システムに対して、古典解の大域的存在を確立すること。
• 以前の研究で必要とされた等しい拡散係数（d1 = d2）の条件を一般化し、任意の正の拡散比 k = d2/d1 を許容すること。
• 化学走性係数 χ と拡散比 k が解の有限時刻爆発を防ぐ役割を解明すること。
• Winkler の手法を一般の正の係数 d1, d2, c1, c2 および非定数 k を持つ系に拡張すること。
• 先行研究における p = 1 に相当する境界状態の感受性関数 χ(v) ∝ 1/(v + c) に対する大域的解の存在という未解決問題に取り組むこと。

提案手法

• パラメータ χ = χ0/d1、k = d2/d1、α = c1/d1、β = c2/d1 を用いて、系を次元なし形式にスケーリングする。
• 最大 Lp-正則性推定と熱半群の滑らかさ性質を適用し、非線形化学走性項 χu∇v/(v + c) を制御する。
• 重み付き積分推定と改良された Hardy–Sobolev 不等式を用いて、u と ∇v の成長を制御する。
• χ ∈ (χ1, χ2) の条件下で、反復的 Lp-推定と Sobolev 嵌め込みを用いて u に対する事前 L∞ 界を確立し、k が大きい場合に有効である。
• µn+1 = (1/µn − (2/N − (χ − (k−1)/2)²/k)⁻¹)⁻¹ で定義される反復的スキーム µn を用いて、Lp-界を L∞ に伝搬する。
• 解の非負性と反応項 −c1v + c2u のおかげで v が一様有界であることに着目する。

実験結果

リサーチクエスチョン

• RQ1 化学走性係数 χ と拡散比 k = d2/d1 がどのような条件下で、飽和対数感受性関数を有する完全な放物型 Keller–Segel システムに対して大域的古典解をもたらすか？
• RQ2 技術的制限により以前は d1 = d2（つまり k = 1）を仮定していたが、この制限を超えて大域的解の結果を一般化できるか？
• RQ3 感受性関数 φ(v) = ln(v + c) の飽和パラメータ c が、低濃度 v における無限大の走性フラックスを引き起こす爆発を防ぐ役割を果たすか？
• RQ4 k が大きい場合に必要な下界 χ > χ1 は、この手法において技術的制限として不可避であるが、今後の解析で除去可能か？
• RQ5 大域的解が保証された場合でも、解は時間に一様有界であるか、それとも t に依存する有界性にとどまるか？

主な発見

• k ∈ (k1, k2) のとき χ ∈ (0, χ2)、k ∈ [k2, ∞) のとき χ ∈ (χ1, χ2) ならば、すべての t > 0 に対して大域的古典解が存在する。ここで k1, k2, χ1, χ2 は N と k の関数として定義される。
• 臨界閾値 χ2 = (k−1)/2 + √(2k/N) は k = 1 のとき √(2/N) に簡略化され、既存の文献と一致する。
• 解 (u, v) はすべての t > 0 に対して ∥u(·, t)∥L∞(Ω) ≤ C(t) および ∥v(·, t)∥L∞(Ω) ≤ C(t) を満たし、C(t) は t に依存するが、有限時刻で発散しない。
• 本手法により、以前の研究で技術的制限として要求されていた d1 = d2（つまり k = 1）の条件が明確に排除され、一般の正の拡散係数が許容される。
• k が大きい場合の下界 χ > χ1 は、このアプローチにおいて技術的に必要とされるが、将来的な技術で除去可能であると予想される。
• 解析により、飽和感受性関数 φ(v) = ln(v + c) が、χ が小さいとは限らない場合でも、v = 0 付近での走性フラックスを抑制することで爆発を効果的に防いでいることが確認された。