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

Basic Information

  • Paper ID: 2510.13283
  • Title: An Allen-Cahn tumor growth model with temperature
  • Authors: Stefania Gatti, Erica Ipocoana, Alain Miranville
  • Classification: math.AP (Analysis of PDEs)
  • Publication Date: October 16, 2025
  • Paper Link: https://arxiv.org/abs/2510.13283

Abstract

This paper proposes a novel non-isothermal Allen-Cahn (Ginzburg-Landau) tumor growth model. After deriving the model using the microforce method, the well-posedness problem is investigated. In particular, the existence and uniqueness of local and global time solutions to the PDE system are established.

Research Background and Motivation

Problem Background

  1. Importance of Tumor Growth Modeling: Mathematical models play an increasingly important role in cancer research, providing the medical community with better understanding and support for cancer treatment, particularly in studying tumor growth mechanisms, parameter analysis, and long-term behavior.
  2. Application of Phase-Field Models: Phase-field models are widely applied in tumor growth research, regulating tumor evolution through order parameters driven by biological mechanisms such as proliferation, apoptosis, and chemotaxis.
  3. Neglected Temperature Effects: Despite extensive mathematical models for tumor growth, an important aspect has been overlooked—the thermal effects of temperature on cancer. Experiments confirm that both high and low temperatures can cause partial or complete destruction of tumor cells.

Research Motivation

  1. Filling Research Gaps: The existing literature lacks Allen-Cahn tumor growth models considering temperature effects; the only non-isothermal phase-field tumor models are within the Cahn-Hilliard framework.
  2. Overcoming Technical Difficulties: In the Cahn-Hilliard framework, due to the presence of chemical potential terms, only the existence of weak entropy solutions can be proven, preventing uniqueness of solutions. The Allen-Cahn system overcomes such difficulties.
  3. Thermodynamic Consistency: Establishing a thermodynamically consistent model that incorporates temperature effects on tumor growth while achieving global existence and uniqueness of non-weak solutions.

Core Contributions

  1. Proposing a Novel Non-Isothermal Allen-Cahn Tumor Growth Model: First incorporation of temperature effects into the Allen-Cahn framework for tumor growth modeling.
  2. Rigorous Mathematical Derivation: Using Gurtin's microforce method, the coupled PDE system is rigorously derived from the first and second laws of thermodynamics.
  3. Complete Well-Posedness Theory: Establishing existence and uniqueness of local and global time solutions, overcoming technical difficulties of Cahn-Hilliard models.
  4. Successful Application of Decoupling Strategy: Innovative decoupling methods combined with Schauder fixed-point theorem successfully handle the nonlinearly coupled system.
  5. Preservation of Physical Meaning: Proving non-negativity and boundedness of solutions, consistent with biological significance.

Detailed Methodology

Problem Formulation

Investigation of well-posedness for the non-isothermal Allen-Cahn tumor growth model on a three-dimensional bounded smooth domain Ω, with the system comprising three evolution equations:

\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}$$ where $\phi$ is the order parameter (tumor concentration), $\theta$ is absolute temperature, and $\sigma$ is nutrient concentration. ### Model Derivation #### 1. Microforce Balance Equation Based on Gurtin's microforce method, the internal microforce balance is established: $$\text{div } \boldsymbol{\xi} + \pi + \gamma = 0$$ where $\boldsymbol{\xi}$ is the microstress, $\pi$ is the internal microforce, and $\gamma$ is the external microforce. #### 2. Thermodynamic Relations - **Internal Energy Density**: $e = \psi + \theta s$ (Gibbs relation) - **Entropy**: $s = -\frac{\partial \psi}{\partial \theta}$ - **Free Energy**: $\psi = \frac{\varepsilon}{2}|\nabla\phi|^2 + \frac{1}{\varepsilon}F(\phi) + f(\theta) - \theta\phi + N(\phi)$ #### 3. Application of Thermodynamic Laws Through the first law of thermodynamics and the Clausius-Duhem inequality: - Microstress relation: $\boldsymbol{\xi} = \varepsilon\nabla\phi$ - Constitutive relation: $\frac{1}{\theta}(\pi + \frac{\partial\psi}{\partial\phi}) = -\beta\frac{\partial\phi}{\partial t}$ ### Technical Innovations #### 1. Decoupling Strategy A two-step decoupling approach is employed: - **Step One**: Fix temperature $\tilde{\theta}$, solve the $(\phi, \sigma)$ system - **Step Two**: Fix $(\phi, \sigma)$, solve the temperature equation #### 2. Function Space Design - Temperature space: $\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\}$ - Phase-field-nutrient space: Appropriate combinations of Sobolev spaces #### 3. Key Estimation Techniques - **Moser Iteration**: Obtaining $L^\infty$ estimates for the order parameter - **Nonlinear Poincaré Inequality**: Handling degenerate terms - **Young and Hölder Inequalities**: Controlling nonlinear terms ## Experimental Setup ### Mathematical Framework - **Domain**: $\Omega \subset \mathbb{R}^3$ is a smooth bounded domain - **Boundary Conditions**: Homogeneous Neumann boundary conditions - **Initial Conditions**: - $\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$ ### Model Parameters - **Physical Parameters**: $P, A, C, B > 0$ (proliferation, apoptosis, consumption, and transfer rates) - **Heat Conductivity Coefficient**: $\kappa(\theta) = 1 + \theta^q$, $q \geq 2$ - **Double-Well Potential**: $F(\phi) = \phi^2(1-\phi)^2$ - **Regulatory Function**: $h \in C^1(\mathbb{R})$, monotonically increasing, $h(0) = 0$ ## Experimental Results ### Main Theorems #### Theorem 3.1 (Local Existence and Uniqueness) For any $R > 0$, there exists $T = T(R) > 0$ such that when initial data satisfy the conditions and $\|\theta_0\|_{V \cap L^{3q}(\Omega)} + \|K(\theta_0)\|_V + \|\phi_0\|_{W^{5/3,6}(\Omega)} \leq R$, the system admits a unique local solution $(\theta, \phi, \sigma)$ satisfying: - $\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')$ #### Theorem 3.6 (Global Existence) The local solution can be extended to a global solution; that is, for any $T > 0$, the solution exists on the entire time interval $(0,T)$. ### Key Estimates #### 1. A Priori Estimates Through energy methods: $$\|\phi\|_{L^\infty(0,T;V) \cap L^2(0,T;H^2(\Omega))} + \|\sigma\|_{L^\infty(0,T;H) \cap L^2(0,T;V)} \leq m_0(R)m_1(T)$$ #### 2. Positivity Preservation - $\theta \geq 0$ a.e. in $Q_T$ - $\phi \geq 0$ a.e. in $Q_T$ - $0 \leq \sigma \leq 1$ a.e. in $Q_T$ #### 3. $L^\infty$ Estimates Proven via Moser iteration: $$\|\phi\|_{L^\infty(0,T;L^\infty(\Omega))} \leq m_0(R)m_1(T)$$ ## Related Work ### Allen-Cahn vs. Cahn-Hilliard Models - **Cahn-Hilliard Model**: Conserves mass but only yields weak solutions - **Allen-Cahn Model**: Does not conserve mass but yields strong solutions and uniqueness ### Non-Isothermal Phase-Field Models - **Fluid Mixtures**: Work by Eleuteri, Gatti, Schimperna, et al. - **Tumor Growth**: Ipocoana's Cahn-Hilliard model (weak solutions only) ### Allen-Cahn Models for Tumor Growth - Colli et al.: Prostate cancer model - Xu et al.: Vascularized tumor growth - This paper: First Allen-Cahn model incorporating temperature effects ## Conclusions and Discussion ### Main Conclusions 1. **Successfully established** the first non-isothermal Allen-Cahn tumor growth model 2. **Rigorously proved** existence and uniqueness of local and global solutions 3. **Preserved physical meaning**: Non-negativity and boundedness of solutions 4. **Overcame technical difficulties**: Advantages over Cahn-Hilliard models ### Limitations 1. **Simplified Assumptions**: Does not consider thermal control, necrosis, and other complex biological processes 2. **Parameter Restrictions**: Special form of heat conductivity coefficient $\kappa(\theta) = 1 + \theta^q$ 3. **Dimensional Restriction**: Only considers three-dimensional case 4. **Boundary Conditions**: Only handles homogeneous Neumann conditions ### Future Directions 1. **Biological Extensions**: Incorporating additional biological mechanisms (necrosis, immune response, etc.) 2. **Numerical Methods**: Developing efficient numerical algorithms 3. **Parameter Identification**: Determining model parameters from experimental data 4. **Optimal Control**: Studying temperature-controlled tumor treatment strategies ## In-Depth Evaluation ### Strengths 1. **Theoretical Rigor**: Complete mathematical derivation and rigorous well-posedness proofs 2. **Methodological Innovation**: Successful decoupling strategy and technical handling 3. **Physical Reasonableness**: Model with solid thermodynamic foundation 4. **Mathematical Sophistication**: Skillful application of Moser iteration, Schauder fixed-point theorem, etc. ### Weaknesses 1. **Missing Experimental Validation**: Lack of numerical experiments and biological verification 2. **Model Simplification**: Relatively simple biological assumptions 3. **Technical Limitations**: Restrictive requirements on certain parameters and conditions ### Impact 1. **Theoretical Contribution**: Establishing mathematical foundation for non-isothermal tumor growth models 2. **Methodological Value**: Decoupling strategy applicable to other coupled PDE systems 3. **Application Prospects**: Providing mathematical modeling support for temperature-based therapy ### Applicable Scenarios 1. **Theoretical Research**: PDE theory, phase-field model analysis 2. **Biomathematics**: Tumor growth modeling 3. **Medical Applications**: Mathematical modeling of thermotherapy and cryotherapy ## References The paper cites 35 important references, covering: - Phase-field theory foundations (Gurtin, Frémond, et al.) - Tumor growth models (Byrne, Wise, et al.) - PDE theory (Ladyzhenskaja, et al.) - Medical research on thermotherapy --- **Overall Assessment**: This is a high-quality mathematical analysis paper that achieves important theoretical breakthroughs, providing a rigorous mathematical foundation for non-isothermal tumor growth models. While lacking in experimental validation, its theoretical contributions and methodological innovations possess significant academic value.