[论文解读] A hybrid variational principle for the Keller-Segel system in $\mathbb R^2$

Adrien Blanchet, José A. Carrillo+4 more|arXiv (Cornell University)|Jul 21, 2014

Mathematical Biology Tumor Growth参考文献 22被引用 4

一句话总结

本文提出一种混合变分原理，用于构造二维Keller-Segel系统在总质量低于临界阈值 $8\pi$ 时的全局弱解。通过在Wasserstein乘积空间中利用梯度流结构，作者将Jordan-Kinderlehrer-Otto最小化方案推广至处理耦合的抛物-抛物系统，确保质量守恒且长期存在而不发生爆破。

ABSTRACT

We construct weak global in time solutions to the classical Keller-Segel system cell movement by chemotaxis in two dimensions when the total mass is below the well-known critical value. Our construction takes advantage of the fact that the Keller-Segel system can be realized as a gradient flow in a suitable functional product space. This allows us to employ a hybrid variational principle which is a generalisation of the minimising implicit scheme for Wasserstein distances introduced by Jordan, Kinderlehrer and Otto (1998).

研究动机与目标

• 在 $\mathbb{R}^2$ 中，当总质量低于临界值 $8\pi$ 时，建立抛物-抛物Keller-Segel系统全局弱解的存在性。
• 发展一种变分框架，以捕捉该系统在概率测度与 $H^1$ 函数乘积空间中的梯度流结构。
• 将Jordan、Kinderlehrer和Otto的最小化隐式格式推广至适用于具有耦合动力学的趋化模型的混合设置。
• 确保所构造的解保持质量守恒，并避免因非线性漂移与扩散作用导致的有限时间爆破。
• 提供一种严格的变分构造，所得自由能解具有改进的正则性与紧致性性质。

提出的方法

• 将Keller-Segel系统表述为在概率测度与 $H^1$ 函数乘积空间中的梯度流，使用自由能泛函 $\mathcal{E}$。
• 应用结合时间离散化与Wasserstein距离最小化的混合变分方案，以在每一步时间定义增量极小化子。
• 利用Matthes-McCann-Savare的流交换技术，推导离散方案中极小化子的正则性估计。
• 通过改进的Biler-Hebisch-Nadzieja不等式与Carleman型估计，建立离散能量与耗散的统一有界性。
• 通过命题C.1中向量场紧致性的应用，利用窄收敛与 $L^2$ 中弱收敛，证明插值解的紧致性。
• 在时间离散格式中取极限，得到满足连续性方程与输运方程分布意义下的全局弱解。

实验结果

研究问题

• RQ1 当总质量低于临界阈值 $8\pi$ 时，能否为二维抛物-抛物Keller-Segel系统构造全局弱解？
• RQ2 如何将Jordan-Kinderlehrer-Otto最小化方案推广至处理具有耦合抛物方程与非局部相互作用的系统？
• RQ3 Keller-Segel系统背后支撑的变分结构是什么，使其能够实现稳定的时间离散化与能量耗散控制？
• RQ4 从变分方案中涌现出哪些紧致性与正则性性质，以确保收敛至弱解？
• RQ5 混合变分原理如何在临界质量区域中保持质量守恒并防止有限时间爆破？

主要发现

• 当趋化敏感度 $\chi$ 低于 $8\pi$（即临界质量阈值）时，二维抛物-抛物Keller-Segel系统存在全局弱解。
• 所构造的解为自由能解，保持总质量守恒，且在时间上保持有界，避免有限时间爆破。
• 混合变分原理确保了离散能量耗散，并对密度及其对数导数的 $L^2$-范数提供了统一有界性。
• 通过改进的Biler-Hebisch-Nadzieja不等式与Carleman型估计，控制熵与对数矩，从而实现统一估计。
• 通过测度的窄收敛与向量场的弱收敛，建立了插值解的紧致性，从而实现向连续极限的过渡。
• 该方法构造的解在所有时刻均为概率测度，满足 $\rho(t,\cdot) \in \mathcal{P}(\mathbb{R}^2)$，确保物理一致性。