# American Institute of Mathematical Sciences

• Previous Article
The maximum surplus before ruin in a jump-diffusion insurance risk process with dependence
• DCDS-B Home
• This Issue
• Next Article
Comparison theorem and correlation for stochastic heat equations driven by Lévy space-time white noises
July  2019, 24(7): 3011-3035. doi: 10.3934/dcdsb.2018297

## Analysis and computation of some tumor growth models with nutrient: From cell density models to free boundary dynamics

 1 Department of Mathematics and Department of Physics, Duke University, Durham, NC 27708-0320, USA 2 Institute of Natural Sciences and Department of Mathematics, Shanghai Jiaotong University, Shanghai, 200240, China 3 School of Mathematics, University of Minnesota, Minneapolis, MN 55455, USA 4 Beijing International Center for Mathematical Research, Peking University, Beijing 100871, China

* Corresponding author: Zhennan Zhou

Received  February 2018 Revised  July 2018 Published  October 2018

Fund Project: J. Liu is partially supported by KI-Net NSF RNMS grant No.11-07444, NSF grant DMS-1812573 and NSF grant DMS-1514826. M. Tang is supported by Science Challenge Project No. TZZT2017-A3-HT003-F and NSFC 91330203. Z. Zhou is partially supported by RNMS11-07444 (KI-Net) and the start up grant from Peking University. L. Wang is partially supported by the start up grant from SUNY Buffalo and NSF grant DMS-1620135.

In this paper, we study a tumor growth equation along with various models for the nutrient component, including a in vitro model and a 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 $γ \to ∞$, 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.

Citation: Jian-Guo Liu, Min Tang, Li Wang, Zhennan Zhou. Analysis and computation of some tumor growth models with nutrient: From cell density models to free boundary dynamics. Discrete & Continuous Dynamical Systems - B, 2019, 24 (7) : 3011-3035. doi: 10.3934/dcdsb.2018297
##### References:

show all references

##### References:
Example 1: expanding disk with constant nutrient and initial data (5.55). Left: plot of solution at time $t = 0.5$ with different $\gamma = 20,\ 40, \ 80$. Here $\Delta r = 0.05$, and $\Delta t = 5\times 10^{-5}$ for $\gamma = 20, \ 40$ and $\Delta t = 2.5\times 10^{-5}$ for $\gamma = 80$. Right: comparison of the numerical solution with $\gamma = 80$ with the analytical solution (5.56) at different times $t = 0.0975$, $t = 0.2975$, $t = 0.4975$. Here the black solid curve is the numerical solution and the red dashed curve is the analytical solution
Example 2: a single annulus with constant nutrient and initial data (5.57). Left: plot of solution at time $t = 0.6$ with different $\gamma = 20,\ 40, \ 80$. Here $\Delta r = 0.05$, and $\Delta t = 2.5\times 10^{-5}$. Right: comparison of the numerical solution with $\gamma = 80$ with the analytical solution (5.58) at different times $t = 0.2494$, $t = 0.4994$, $t = 0.8$. Here we use $\Delta r = 0.025$ and $\Delta t = 6.25 \times 10^{-6}$. The black solid curve is the numerical solution and the red dashed curve is the analytical solution
Example 3: a double annulus with constant nutrient and initial data (5.59). Here we compare the numerical solution (black solid curve) and analytical solution (red dashed curve) at time $t = 0.2495$ (left) and $t = 0.6$ (right). Here we use $\Delta r = 0.025$ and $\Delta t = 5\times 10^{-6}$
Example 4: a 1D in vitro model with linear growing function. Left: plots of $n$ at time $t = 0.5$ with various $\gamma = 20, \ 40, \ 80$. The red curve is the analytical solution (5.61). Here $\Delta x = 0.05$ and $\Delta t = 2.5e-5$
Example 5: a 1D in vivo model with linear growing function. Left: plots of $n$ at time $t = 0.5$ with various $\gamma = 20, \ 40, \ 80$. The red curve is the analytical solution (5.61). Here $\Delta x = 0.05$ and $\Delta t = 2.5\times 10^{-5}$
A comparison of the front propagation speed for 1D in vitro model and in vivo model. The dots represent the position of the right boundary at each time, and the curves are computed via (3.37) and (3.39). Here $\Delta x = 0.05$, $\Delta t = 2.5 \times 10^{-5}$, $\gamma = 80$
A comparison of the front propagation speed in the 2D radial symmetric in vitro model and in vivo model. The dots represent the the position of the right boundary at each time, and the curve are computed via (3.43) and (3.44). Here $\Delta x = 0.05$, $\Delta t = 2.5\times 10^{-5}$, $\gamma = 80$
Plot of $n$ at four different times with initial data (5.63). From left to right, up to down, $t = 0$, $t = 0.0177$, $t = 0.0311$, $t = 0.05$
Plot of $n$ at four different times with initial data (5.64). From left to right, up to down, $t = 0$, $t = 0.0177$, $t = 0.0311$, $t = 0.05$
 [1] Hyung Ju Hwang, Youngmin Oh, Marco Antonio Fontelos. The vanishing surface tension limit for the Hele-Shaw problem. Discrete & Continuous Dynamical Systems - B, 2016, 21 (10) : 3479-3514. doi: 10.3934/dcdsb.2016108 [2] A. Bernardini, J. Bragard, H. Mancini. Synchronization between two Hele-Shaw Cells. Mathematical Biosciences & Engineering, 2004, 1 (2) : 339-346. doi: 10.3934/mbe.2004.1.339 [3] Nataliya Vasylyeva, Vitalii Overko. The Hele-Shaw problem with surface tension in the case of subdiffusion. Communications on Pure & Applied Analysis, 2016, 15 (5) : 1941-1974. doi: 10.3934/cpaa.2016023 [4] Elena Izquierdo-Kulich, José Manuel Nieto-Villar. Mesoscopic model for tumor growth. Mathematical Biosciences & Engineering, 2007, 4 (4) : 687-698. doi: 10.3934/mbe.2007.4.687 [5] Shihe Xu. Analysis of a delayed free boundary problem for tumor growth. Discrete & Continuous Dynamical Systems - B, 2011, 15 (1) : 293-308. doi: 10.3934/dcdsb.2011.15.293 [6] Yaodan Huang, Zhengce Zhang, Bei Hu. Bifurcation from stability to instability for a free boundary tumor model with angiogenesis. Discrete & Continuous Dynamical Systems - A, 2019, 39 (5) : 2473-2510. doi: 10.3934/dcds.2019105 [7] Didier Bresch, Thierry Colin, Emmanuel Grenier, Benjamin Ribba, Olivier Saut. A viscoelastic model for avascular tumor growth. Conference Publications, 2009, 2009 (Special) : 101-108. doi: 10.3934/proc.2009.2009.101 [8] Shihe Xu, Yinhui Chen, Meng Bai. Analysis of a free boundary problem for avascular tumor growth with a periodic supply of nutrients. Discrete & Continuous Dynamical Systems - B, 2016, 21 (3) : 997-1008. doi: 10.3934/dcdsb.2016.21.997 [9] Zejia Wang, Suzhen Xu, Huijuan Song. Stationary solutions of a free boundary problem modeling growth of angiogenesis tumor with inhibitor. Discrete & Continuous Dynamical Systems - B, 2018, 23 (6) : 2593-2605. doi: 10.3934/dcdsb.2018129 [10] 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 [11] Junde Wu, Shangbin Cui. Asymptotic behavior of solutions for parabolic differential equations with invariance and applications to a free boundary problem modeling tumor growth. Discrete & Continuous Dynamical Systems - A, 2010, 26 (2) : 737-765. doi: 10.3934/dcds.2010.26.737 [12] Shihe Xu, Meng Bai, Fangwei Zhang. Analysis of a free boundary problem for tumor growth with Gibbs-Thomson relation and time delays. Discrete & Continuous Dynamical Systems - B, 2018, 23 (9) : 3535-3551. doi: 10.3934/dcdsb.2017213 [13] T.L. Jackson. A mathematical model of prostate tumor growth and androgen-independent relapse. Discrete & Continuous Dynamical Systems - B, 2004, 4 (1) : 187-201. doi: 10.3934/dcdsb.2004.4.187 [14] J. Ignacio Tello. On a mathematical model of tumor growth based on cancer stem cells. Mathematical Biosciences & Engineering, 2013, 10 (1) : 263-278. doi: 10.3934/mbe.2013.10.263 [15] Hyun Geun Lee, Yangjin Kim, Junseok Kim. Mathematical model and its fast numerical method for the tumor growth. Mathematical Biosciences & Engineering, 2015, 12 (6) : 1173-1187. doi: 10.3934/mbe.2015.12.1173 [16] Mohammad A. Tabatabai, Wayne M. Eby, Karan P. Singh, Sejong Bae. T model of growth and its application in systems of tumor-immune dynamics. Mathematical Biosciences & Engineering, 2013, 10 (3) : 925-938. doi: 10.3934/mbe.2013.10.925 [17] Marcelo E. de Oliveira, Luiz M. G. Neto. Directional entropy based model for diffusivity-driven tumor growth. Mathematical Biosciences & Engineering, 2016, 13 (2) : 333-341. doi: 10.3934/mbe.2015005 [18] Kentarou Fujie. Global asymptotic stability in a chemotaxis-growth model for tumor invasion. Discrete & Continuous Dynamical Systems - S, 2020, 13 (2) : 203-209. doi: 10.3934/dcdss.2020011 [19] Jia-Feng Cao, Wan-Tong Li, Fei-Ying Yang. Dynamics of a nonlocal SIS epidemic model with free boundary. Discrete & Continuous Dynamical Systems - B, 2017, 22 (2) : 247-266. doi: 10.3934/dcdsb.2017013 [20] Wenzhen Gan, Peng Zhou. A revisit to the diffusive logistic model with free boundary condition. Discrete & Continuous Dynamical Systems - B, 2016, 21 (3) : 837-847. doi: 10.3934/dcdsb.2016.21.837

2018 Impact Factor: 1.008