# American Institute of Mathematical Sciences

## Cahn–Hilliard–Brinkman systems for tumour growth

 1 Fakultät für Mathematik, Universität Regensburg, 93040 Regensburg, Germany 2 Department of Mathematics, University of Trento, Trento, Italy

* Corresponding author: Harald Garcke

Received  March 2020 Revised  January 2021 Early access  March 2021

A phase field model for tumour growth is introduced that is based on a Brinkman law for convective velocity fields. The model couples a convective Cahn–Hilliard equation for the evolution of the tumour to a reaction-diffusion-advection equation for a nutrient and to a Brinkman–Stokes type law for the fluid velocity. The model is derived from basic thermodynamical principles, sharp interface limits are derived by matched asymptotics and an existence theory is presented for the case of a mobility which degenerates in one phase leading to a degenerate parabolic equation of fourth order. Finally numerical results describe qualitative features of the solutions and illustrate instabilities in certain situations.

Citation: Matthias Ebenbeck, Harald Garcke, Robert Nürnberg. Cahn–Hilliard–Brinkman systems for tumour growth. Discrete & Continuous Dynamical Systems - S, doi: 10.3934/dcdss.2021034
##### References:

show all references

##### References:
Typical situation for the formal asymptotic analysis
The tumour and healthy regions $\Omega_T$ and $\Omega_H$
Schematic sketch of the inner region close to $\Sigma(0)$
Initial tumour size for initial data $r$: A slightly perturbed sphere
Comparison of Cahn--Hilliard--Darcy and Cahn--Hilliard--Brinkman models: Tumour at time $t = 12$ for $\beta = 0.1$, left side for the CHD model, right side for the CHB model with $\eta = 10^{-5}$
Influence of different mobilities: Tumour at time $t = 9$ for $\eta = 10^{-5}$, $\beta = 0.1$ and $\alpha = \rho_S = 2$, but with different mobilities, left $m(\varphi) = \frac{1}{2}(1+\varphi)^2$, middle $m(\varphi) = \epsilon$, right $m(\varphi) = 10^{-3}\epsilon$
Influence of the adhesion parameter $\beta$: Evolution of the tumour with $m(\varphi) = \tfrac{1}{2}(1+\varphi)^2$ and $\eta = 0.1$, above for $\beta = 0.1$ at time $t = 1, 3, 6, 10$, below for $\beta = 0.01$ at time $t = 1, 1.5, 2, 2.5$
Influence of viscosiy I: Tumour and velocity for $\beta = 0.01$ at time $t = 2.5$, left for $\eta = 0.1$, right for $\eta = 100$, on top the tumour and below the velocity magnitude
Influence of viscosity II: Evolution of the tumour at time $t = 1, 3, 6, 10$ with $\beta = 0.1$, $\nu = 0$ and a no-slip boundary condition on the left boundary, on top for $\eta = 0.1$ and below for $\eta = 10$
Velocity profiles for different viscosities: The velocity magnitude at time $t = 10$ with $\beta = 0.1$, $\nu = 0$ and a no-slip boundary condition on the left boundary, left for $\eta = 0.1$, right for $\eta = 10$
Influence of viscosity contrast: Tumour at time $t = 10$ with $\beta = 0.1$, $\nu = 0$ and a no-slip b. c. on the left boundary, with $\eta_- = 0.01$, $\eta_+ = 1$; $\eta_- = 1$, $\eta_+ = 0.01$; $\eta_- = 0.01$, $\eta_+ = 10$; $\eta_- = 10$, $\eta_+ = 0.01$
Influence of initial profile I: Tumour at time $t = 0, 0.3, 1, 1.6$ with $\eta = 100$ and with the initial profile corresponding to $r_1$
Influence of initial profile II: Tumour at time $t = 0, 0.7, 1.2, 2.6$ with $\eta = 100$ and with the initial profile corresponding to $r_2$
Influence of initial profile III: Tumour at time $t = 0, 0.5, 1.2, 2.4$ with $\eta = 100$ and with the initial profile corresponding to $r_3$
Influence of initial profile IV: Tumour at time $t = 0, 0.3, 1.3, 2.3$ with $\eta = 0.01$ and with the initial profile corresponding to $r_4$
 [1] Andrea Signori. Penalisation of long treatment time and optimal control of a tumour growth model of Cahn–Hilliard type with singular potential. Discrete & Continuous Dynamical Systems, 2021, 41 (6) : 2519-2542. doi: 10.3934/dcds.2020373 [2] Pierluigi Colli, Gianni Gilardi, Danielle Hilhorst. On a Cahn-Hilliard type phase field system related to tumor growth. Discrete & Continuous Dynamical Systems, 2015, 35 (6) : 2423-2442. doi: 10.3934/dcds.2015.35.2423 [3] Andrea Signori. Optimal treatment for a phase field system of Cahn-Hilliard type modeling tumor growth by asymptotic scheme. Mathematical Control & Related Fields, 2020, 10 (2) : 305-331. doi: 10.3934/mcrf.2019040 [4] Pierluigi Colli, Gianni Gilardi, Elisabetta Rocca, Jürgen Sprekels. Asymptotic analyses and error estimates for a Cahn-Hilliard type phase field system modelling tumor growth. Discrete & Continuous Dynamical Systems - S, 2017, 10 (1) : 37-54. doi: 10.3934/dcdss.2017002 [5] Monica Conti, Stefania Gatti, Alain Miranville. A singular cahn-hilliard-oono phase-field system with hereditary memory. Discrete & Continuous Dynamical Systems, 2018, 38 (6) : 3033-3054. doi: 10.3934/dcds.2018132 [6] Cristina Pocci. On singular limit of a nonlinear $p$-order equation related to Cahn-Hilliard and Allen-Cahn evolutions. Evolution Equations & Control Theory, 2013, 2 (3) : 517-530. doi: 10.3934/eect.2013.2.517 [7] Elisabeth Logak, Chao Wang. The singular limit of a haptotaxis model with bistable growth. Communications on Pure & Applied Analysis, 2012, 11 (1) : 209-228. doi: 10.3934/cpaa.2012.11.209 [8] John R. King, Judith Pérez-Velázquez, H.M. Byrne. Singular travelling waves in a model for tumour encapsulation. Discrete & Continuous Dynamical Systems, 2009, 25 (1) : 195-230. doi: 10.3934/dcds.2009.25.195 [9] A. Jiménez-Casas. Invariant regions and global existence for a phase field model. Discrete & Continuous Dynamical Systems - S, 2008, 1 (2) : 273-281. doi: 10.3934/dcdss.2008.1.273 [10] Murat Uzunca, Ayşe Sarıaydın-Filibelioǧlu. Adaptive discontinuous galerkin finite elements for advective Allen-Cahn equation. Numerical Algebra, Control & Optimization, 2021, 11 (2) : 269-281. doi: 10.3934/naco.2020025 [11] Leonid Berlyand, Mykhailo Potomkin, Volodymyr Rybalko. Sharp interface limit in a phase field model of cell motility. Networks & Heterogeneous Media, 2017, 12 (4) : 551-590. doi: 10.3934/nhm.2017023 [12] Maciek D. Korzec, Hao Wu. Analysis and simulation for an isotropic phase-field model describing grain growth. Discrete & Continuous Dynamical Systems - B, 2014, 19 (7) : 2227-2246. doi: 10.3934/dcdsb.2014.19.2227 [13] Kelei Wang. The singular limit problem in a phase separation model with different diffusion rates $^*$. Discrete & Continuous Dynamical Systems, 2015, 35 (1) : 483-512. doi: 10.3934/dcds.2015.35.483 [14] Honghu Liu. Phase transitions of a phase field model. Discrete & Continuous Dynamical Systems - B, 2011, 16 (3) : 883-894. doi: 10.3934/dcdsb.2011.16.883 [15] Georgia Karali, Yuko Nagase. On the existence of solution for a Cahn-Hilliard/Allen-Cahn equation. Discrete & Continuous Dynamical Systems - S, 2014, 7 (1) : 127-137. doi: 10.3934/dcdss.2014.7.127 [16] T. Tachim Medjo. A Cahn-Hilliard-Navier-Stokes model with delays. Discrete & Continuous Dynamical Systems - B, 2016, 21 (8) : 2663-2685. doi: 10.3934/dcdsb.2016067 [17] Alain Miranville, Sergey Zelik. The Cahn-Hilliard equation with singular potentials and dynamic boundary conditions. Discrete & Continuous Dynamical Systems, 2010, 28 (1) : 275-310. doi: 10.3934/dcds.2010.28.275 [18] Xinfu Chen, G. Caginalp, Christof Eck. A rapidly converging phase field model. Discrete & Continuous Dynamical Systems, 2006, 15 (4) : 1017-1034. doi: 10.3934/dcds.2006.15.1017 [19] Dimitra Antonopoulou, Georgia Karali. Existence of solution for a generalized stochastic Cahn-Hilliard equation on convex domains. Discrete & Continuous Dynamical Systems - B, 2011, 16 (1) : 31-55. doi: 10.3934/dcdsb.2011.16.31 [20] Tania Biswas, Elisabetta Rocca. Long time dynamics of a phase-field model of prostate cancer growth with chemotherapy and antiangiogenic therapy effects. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021140

2020 Impact Factor: 2.425