2025-11-24T04:19:17.762191

An Allen-Cahn tumor growth model with temperature

Gatti, Ipocoana, Miranville

In this paper, we propose a new non-isothermal Allen-Cahn (Ginzburg-Landau) model for tumor growth. After deriving it using a microforces approach, we study its well-posedness. In particular, we are able to prove the existence and uniqueness of a local and global-in-time solution to our PDE system.

academic

An Allen-Cahn tumor growth model with temperature

基本信息

论文ID: 2510.13283
标题: An Allen-Cahn tumor growth model with temperature
作者: Stefania Gatti, Erica Ipocoana, Alain Miranville
分类: math.AP (Analysis of PDEs)
发表时间: October 16, 2025
论文链接: https://arxiv.org/abs/2510.13283

摘要

本文提出了一个新的非等温Allen-Cahn（Ginzburg-Landau）肿瘤生长模型。通过使用微力方法推导该模型后，研究了其适定性问题。特别地，证明了PDE系统局部和全局时间解的存在性和唯一性。

研究背景与动机

问题背景

肿瘤生长建模的重要性：数学模型在癌症研究中发挥着越来越重要的作用，能够为医学界提供更好的理解和治疗癌症的支持，特别是在肿瘤生长机制、参数分析和长期行为研究方面。
相场模型的应用：相场（phase-field）模型在肿瘤生长研究中被广泛应用，通过序参量调节肿瘤演化，并由增殖、凋亡、趋化性等生物机制驱动。
温度效应的缺失：尽管已有大量肿瘤生长数学模型，但一个重要方面一直被忽视——温度对癌症的热效应。实验证实高温和低温都能导致肿瘤细胞的部分或完全破坏。

研究动机

填补研究空白：现有文献中缺乏考虑温度影响的Allen-Cahn肿瘤生长模型，仅有的非等温相场肿瘤模型是Cahn-Hilliard框架。
克服技术困难：在Cahn-Hilliard框架中，由于化学势项的存在，只能证明弱熵解的存在性，无法获得解的唯一性。Allen-Cahn系统能够克服类似困难。
热力学一致性：建立热力学一致的模型，能够包含温度对肿瘤生长的影响，同时实现不太弱解的全局存在性和唯一性。

核心贡献

提出新的非等温Allen-Cahn肿瘤生长模型：首次将温度效应纳入Allen-Cahn框架的肿瘤生长模型中。
严格的数学推导：使用Gurtin的微力方法从热力学第一、第二定律出发，严格推导出耦合的PDE系统。
完整的适定性理论：证明了系统局部和全局时间解的存在性和唯一性，克服了Cahn-Hilliard模型的技术困难。
解耦策略的成功应用：采用创新的解耦方法，结合Schauder不动点定理，成功处理了非线性耦合系统。
物理意义的保持：证明了解的非负性和有界性，符合生物学意义。

方法详解

任务定义

研究三维有界光滑区域Ω上的非等温Allen-Cahn肿瘤生长模型的适定性问题，系统包含三个演化方程：

$\begin{aligned} \phi_t - \Delta\phi + F'(\phi) - \theta &= (P\sigma - A)h(\phi) \quad \text{in } \Omega \times (0,T) \\ \theta_t - \text{div}[\kappa(\theta)\nabla\theta] - \phi_t^2 + \theta\phi_t &= 0 \quad \text{in } \Omega \times (0,T) \\ \sigma_t - \Delta\sigma &= -C\sigma h(\phi) + B(\sigma_B - \sigma) \quad \text{in } \Omega \times (0,T) \end{aligned}$

其中 $\phi$ 为序参量（肿瘤浓度）， $\theta$ 为绝对温度， $\sigma$ 为营养物浓度。

模型推导

1. 微力平衡方程

基于Gurtin的微力方法，建立内部微力平衡： $\text{div } \boldsymbol{\xi} + \pi + \gamma = 0$

其中 $\boldsymbol{\xi}$ 为微应力， $\pi$ 为内部微力， $\gamma$ 为外部微力。

2. 热力学关系

内能密度： $e = \psi + \theta s$ （Gibbs关系）
熵： $s = -\frac{\partial \psi}{\partial \theta}$
自由能： $\psi = \frac{\varepsilon}{2}|\nabla\phi|^2 + \frac{1}{\varepsilon}F(\phi) + f(\theta) - \theta\phi + N(\phi)$

3. 热力学定律应用

通过热力学第一定律和Clausius-Duhem不等式，推导出：

微应力关系： $\boldsymbol{\xi} = \varepsilon\nabla\phi$
本构关系： $\frac{1}{\theta}(\pi + \frac{\partial\psi}{\partial\phi}) = -\beta\frac{\partial\phi}{\partial t}$

技术创新点

1. 解耦策略

采用两步解耦方法：

第一步：固定温度 $\tilde{\theta}$ ，求解 $(\phi, \sigma)$ 系统
第二步：固定 $(\phi, \sigma)$ ，求解温度方程

2. 函数空间设计

温度空间： $\Theta_R(T) = \{\tilde{\theta} \in L^4(0,T;L^6(\Omega)) : \tilde{\theta} \geq 0, \|\tilde{\theta}\|_{L^4(0,T;L^6(\Omega))} \leq R\}$
相场-营养空间：适当的Sobolev空间组合

3. 关键估计技术

Moser迭代：获得序参量的 $L^\infty$ 估计
非线性Poincaré不等式：处理退化项
Young不等式和Hölder不等式：控制非线性项

实验设置

数学框架

区域： $\Omega \subset \mathbb{R}^3$ 为光滑有界区域
边界条件：齐次Neumann边界条件
初始条件：
- $\theta_0 \in L^{3q}(\Omega) \cap V$ ， $K(\theta_0) \in V$ ， $\theta_0 \geq 0$
- $\phi_0 \in W^{5/3,6}(\Omega)$ ， $\phi_0 \geq 0$
- $\sigma_0 \in L^\infty(\Omega)$ ， $0 \leq \sigma_0 \leq 1$

模型参数

物理参数： $P, A, C, B > 0$ （增殖率、凋亡率、消耗率、转移率）
热传导系数： $\kappa(\theta) = 1 + \theta^q$ ， $q \geq 2$
双阱势： $F(\phi) = \phi^2(1-\phi)^2$
调节函数： $h \in C^1(\mathbb{R})$ ，单调递增， $h(0) = 0$

$\theta \in C([0,T];H) \cap L^\infty(0,T;L^{3q}(\Omega) \cap V) \cap L^{4q}(0,T;L^{12q}(\Omega)) \cap H^1(0,T;H)$
$\phi \in C([0,T];V) \cap W^{1,4}(0,T;L^6(\Omega))$
$\sigma \in L^\infty(0,T;L^\infty(\Omega)) \cap L^2(0,T;V) \cap H^1(0,T;V')$