[论文解读] 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

一句话总结

本文建立了具有饱和对数敏感函数 φ(v) = ln(v + c) 的完全抛物型 Keller–Segel 化学趋化系统经典解的全局存在性，其中 c > 0 可避免低化学浓度下奇点的出现。通过加权能量估计与最大 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) = ln(v + c) 的抛物-抛物 Keller–Segel 系统经典解的全局存在性，避免在 v = 0 处 ln v 的奇点问题。
• 通过允许任意正扩散比 k = d2/d1，推广此前要求扩散系数相等（d1 = d2）的结果。
• 分析趋化系数 χ 与扩散比 k 在防止解在有限时间内爆破中的作用。
• 将 Winkler 的技术扩展至具有通解正系数 d1, d2, c1, c2 及非恒定 k 的系统。
• 解决敏感函数 χ(v) ∝ 1/(v + c) 的边界情形下全局存在性的开放问题，该情形对应于先前研究中的 p = 1。

提出的方法

• 引入尺度变换，将系统化为无量纲形式，其中参数为 χ = χ0/d1，k = d2/d1，α = c1/d1，β = c2/d1。
• 应用最大 Lp-正则性估计与热半群的光滑性性质，以控制非线性趋化项 χu∇v/(v + c)。
• 利用加权积分估计与改进的 Hardy–Sobolev 不等式，控制 u 与 ∇v 的增长。
• 通过迭代 Lp-估计与 Sobolev 嵌入，建立 u 的先验 L∞ 有界性，利用当 k 较大时 χ ∈ (χ1, χ2) 的条件。
• 采用基于序列 µn 定义的迭代方案：µn+1 = (1/µn − (2/N − (χ − (k−1)/2)²/k)⁻¹)⁻¹，将 Lp 有界性传播至 L∞。
• 利用 v 的一致有界性，这是由于反应项 −c1v + c2u 及解的非负性所致。

实验结果

研究问题

• 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)，全局经典解对所有 t > 0 存在；当 k ∈ [k2, ∞) 时，若 χ ∈ (χ1, χ2)，全局经典解亦存在，其中 k1, k2, χ1, χ2 均以 N 与 k 表示。
• 临界阈值 χ2 = (k−1)/2 + √(2k/N) 在 k = 1 时简化为 √(2/N)，与文献中已知结果一致。
• 解 (u, v) 满足 ∥u(·, t)∥L∞(Ω) ≤ C(t) 与 ∥v(·, t)∥L∞(Ω) ≤ C(t)，对所有 t > 0 成立，其中 C(t) 依赖于 t 但不会在有限时间内爆破。
• 该方法成功消除了早期工作中存在的技术限制 d1 = d2（即 k = 1），从而允许一般正扩散系数。
• 对于大 k 值，χ > χ1 的下界在本方法中被视为技术上必要，但未来技术有望将其去除。
• 分析证实，饱和敏感函数 φ(v) = ln(v + c) 通过在 v = 0 附近抑制趋化通量，有效防止了爆破，即使 χ 并不很小亦然。