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

一句话总结

本文研究了有界区域中具有二次生产项的化学排斥模型，证明了一维和二维情形下弱解的全局存在性与唯一性。提出了两种无条件能量稳定的首阶时间格式——向后欧拉格式与线性化耦合格式，证明了其可解性、收敛性及误差估计，其中线性格式虽不保持非负性，但仍表现出收敛性。

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.

研究动机与目标

• 建立具有二次生产项的化学排斥模型在有界区域中全局弱解的存在性。
• 在正则性准则下，证明一维与二维区域中解的唯一性与全局正则性。
• 设计并分析两种无条件能量稳定的首阶时间离散格式：非线性向后欧拉格式与线性化耦合格式。
• 证明两种格式的可解性、收敛至弱解的性质，并推导误差估计。
• 通过数值模拟研究格式的行为，特别关注能量稳定性与收敛性特征。

提出的方法

• 建立具有二次生产项 $ f(u) = u^2 $ 的化学排斥系统，引入辅助变量 $ \boldsymbol{\sigma} = \nabla v $ 以实现系统解耦。
• 推导包含 $ u $、$ \boldsymbol{\sigma} $ 与 $ v $ 的重构系统（3），并给出诺伊曼边界条件与初始数据。
• 提出非线性向后欧拉格式（104）进行时间离散化，确保细胞守恒与无条件能量稳定性。
• 引入线性化耦合格式作为计算高效的替代方案，尽管其不保持非负性。
• 应用 Leray-Schauder 不动点定理，通过紧致性与连续性论证，证明时间离散格式解的存在性。
• 利用能量估计与 $ L^4 $-范数有界性，推导收敛性与误差估计，其中非线性格式的唯一性需在强范数下满足小时间步长条件。

实验结果

研究问题

• RQ1 具有二次生产项的化学排斥模型在有界区域中是否存在全局弱解？
• RQ2 在何种条件下弱解具有唯一性与全局正则性，特别是在一维与二维情形下？
• RQ3 能否为该模型构造无条件能量稳定的首阶时间离散格式？
• RQ4 所提格式的收敛性与误差估计为何种形式，尤其是线性化版本？
• RQ5 从数值角度看，这些格式在稳定性、能量衰减与解分布特性方面表现如何？

主要发现

• 在有界区域中，具有二次生产项 $ f(u) = u^2 $ 的化学排斥模型存在全局弱解。
• 在正则性准则下，一维与二维区域中解的唯一性与全局正则性得以确立。
• 非线性向后欧拉格式具有无条件能量稳定性，保持细胞守恒，并以误差估计收敛至弱解。
• 线性化耦合格式具有无条件能量稳定性并收敛至弱解，尽管不保持非负性。
• 非线性格式的唯一性在强范数下需满足小时间步长条件，该条件在 $ nD $ 区域（$ n \leq 2 $）中可简化为小时间步长，因可导出全局强估计。
• 数值模拟验证了两种格式的能量稳定性与收敛行为，与理论分析一致。