[论文解读] Analysis and computation of some tumor growth models with nutrient: from cell density models to free boundary dynamics

Jian‐Guo Liu, Min Tang+2 more|arXiv (Cornell University)|Feb 2, 2018

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

一句话总结

本文分析了具有营养动力学的肿瘤生长模型，将受压力定律 $ p(n) = n^\gamma $ 控制的细胞密度模型与不可压缩极限下 $ \gamma \to \infty $ 时的自由边界 Hele-Shaw 模型联系起来。研究建立了分析估计，推导了肿瘤前缘传播的基准解，并通过一种保守且保持非负性的数值格式验证了从细胞密度模型到自由边界动力学的收敛性。

ABSTRACT

In this paper, we study the tumor growth equation along with various models for the nutrient component, including the \emph{in vitro} model and the \emph{in vivo} model. At the cell density level, the spatial availability of the tumor density $n$ is governed by the Darcy law via the pressure $p(n)=n^γ$. For finite $γ$, we prove some a priori estimates of the tumor growth model, such as boundedness of the nutrient density, and non-negativity and growth estimate of the tumor density. As $γ ightarrow \infty$, the cell density models formally converge to Hele-Shaw flow models, which determine the free boundary dynamics of the tumor tissue in the incompressible limit. We derive several analytical solutions to the Hele-Shaw flow models, which serve as benchmark solutions to the geometric motion of tumor front propagation. Finally, we apply a conservative and positivity preserving numerical scheme to the cell density models, with numerical results verifying the link between cell density models and the free boundary dynamical models.

研究动机与目标

• 建立肿瘤细胞密度模型在营养依赖性增殖下的严格分析估计——包括营养物质的有界性、非负性及增殖控制。
• 形式化证明当 $ \gamma \to \infty $ 时，具有 $ p(n) = n^\gamma $ 的细胞密度模型收敛于 Hele-Shaw 类型的自由边界问题，代表肿瘤组织的不可压缩极限。
• 推导极限 Hele-Shaw 模型的解析解，作为肿瘤前缘传播的基准解。
• 开发并验证一种保守且保持非负性的数值格式，准确捕捉从细胞密度模型到自由边界动力学的过渡。

提出的方法

• 通过由压力 $ p(n) = n^\gamma $ 驱动的退化扩散方程建模肿瘤细胞密度 $ n(x,t) $，其中营养消耗由 $ -\Delta c + \psi(n)c = 0 $ 在体外和体内条件下描述。
• 使用渐近分析推导 $ \gamma \to \infty $ 时的 Hele-Shaw 极限，得到一个自由边界问题，其中 $ p_\infty $ 满足 $ -\Delta p_\infty = G(c) $ 在 $ D_\infty(t) $ 内，且在 $ \partial D_\infty(t) $ 上 $ p_\infty = 0 $。
• 推导体外和体内模型中径向对称情况的解析解，为前缘传播速度提供精确基准。
• 实施一种保守且保持非负性的有限差分格式，数值求解有限 $ \gamma $ 下的细胞密度模型，确保稳定性和物理一致性。
• 通过比较不同 $ \gamma $ 下的前缘传播速度和空间分布，验证数值解向 Hele-Shaw 极限的收敛性。
• 对 $ n $ 和 $ c $ 使用 Neumann 和 Dirichlet 边界条件，通过网格细化确保前缘追踪和营养分布的准确性。

实验结果

研究问题

• RQ1 在营养影响的模型中，肿瘤细胞密度的有界性、非负性及增殖估计如何依赖于压力指数 $ \gamma $？
• RQ2 在 $ \gamma \to \infty $ 的极限下，具有 $ p(n) = n^\gamma $ 的细胞密度模型与 Hele-Shaw 自由边界模型之间是否存在严格的数学联系？
• RQ3 在径向对称的体外和体内设定下，Hele-Shaw 模型是否存在解析解？它们如何预测前缘传播速度？
• RQ4 一种保守且保持非负性的数值格式在多大程度上能准确再现 $ \gamma $ 增大时从细胞密度动力学到自由边界运动的过渡？

主要发现

• 对于有限 $ \gamma $，肿瘤细胞密度 $ n $ 保持非负，且以受控速率增长，营养浓度均匀有界。
• 当 $ \gamma \to \infty $ 时，细胞密度模型在形式上收敛于 Hele-Shaw 自由边界问题，其中肿瘤前缘以速度 $ v = -\nabla p_\infty $ 演化。
• 体外和体内模型中径向对称情况的解析解表明，在长时间极限下，体外模型中的肿瘤前缘传播速度是体内模型的两倍。
• 当 $ \gamma = 80 $ 时的数值模拟显示，细胞密度模型与极限 Hele-Shaw 解之间具有极佳的一致性，验证了前缘动力学的收敛性。
• 体外模型中的前缘传播速度由 $ R(t) \sim \sqrt{2t} $ 预测，而体内模型中为 $ R(t) \sim \sqrt{t} $，与数值结果一致。
• 该数值格式成功保持了非负性并守恒质量，即使在高 $ \gamma $ 条件下也能实现稳定且准确的前缘追踪，但当 $ \gamma $ 较小且密度超过 1 时精度会下降。