# American Institute of Mathematical Sciences

2007, 7(3): 581-604. doi: 10.3934/dcdsb.2007.7.581

## Nonlinear three-dimensional simulation of solid tumor growth

 1 Department of Mathematics, University of California, Irvine, CA 92697, United States 2 School of Health Information Sciences, University of Texas Health Science Center, Houston, TX 77030, United States 3 Center for Mathematical and Computational Biology, Department of Mathematics, University of California, Irvine, CA 92697-3875 4 Department of Mathematics, Center for Mathematical and Computational Biology, University of California, Irvine, CA 92697, United States

Received  September 2006 Revised  January 2007 Published  February 2007

We present a new, adaptive boundary integral method to simulate solid tumor growth in 3-d. We use a reformulation of a classical model that accounts for cell-proliferation, apoptosis, cell-to-cell and cell-to-matrix adhesion. The 3-d method relies on accurate discretizations of singular surface integrals, a spatial rescaling and the use of an adaptive surface mesh. The discretized boundary integral equations are solved iteratively using GMRES and a discretized version of the Dirichlet-Neumann map, formulated in terms of a vector potential, is used to determine the normal velocity of the tumor surface. Explicit time stepping is used to update the tumor surface. We present simulations of the nonlinear evolution of growing tumors. At early times, good agreement is obtained between the results of a linear stability analysis and nonlinear simulations. At later times, linear theory is found to overpredict the growth of perturbations. Nonlinearity results in mode creation and interaction that leads to the formation of dimples and the tumor surface buckles inwards. The morphologic instability allows the tumor to increase its surface area, relative to its volume, thereby allowing the cells in the tumor bulk greater access to nutrient. This in turn allows the tumor to overcome the diffusional limitations on growth and to grow to larger sizes than would be possible if the tumor were spherical. Consequently, instability provides a means for avascular tumor invasion.
Citation: Xiangrong Li, Vittorio Cristini, Qing Nie, John S. Lowengrub. Nonlinear three-dimensional simulation of solid tumor growth. Discrete & Continuous Dynamical Systems - B, 2007, 7 (3) : 581-604. doi: 10.3934/dcdsb.2007.7.581
 [1] Yi Shi, Kai Bao, Xiao-Ping Wang. 3D adaptive finite element method for a phase field model for the moving contact line problems. Inverse Problems & Imaging, 2013, 7 (3) : 947-959. doi: 10.3934/ipi.2013.7.947 [2] Thomas Y. Hou, Pingwen Zhang. Convergence of a boundary integral method for 3-D water waves. Discrete & Continuous Dynamical Systems - B, 2002, 2 (1) : 1-34. doi: 10.3934/dcdsb.2002.2.1 [3] Thomas Y. Hou, Zuoqiang Shi. Dynamic growth estimates of maximum vorticity for 3D incompressible Euler equations and the SQG model. Discrete & Continuous Dynamical Systems - A, 2012, 32 (5) : 1449-1463. doi: 10.3934/dcds.2012.32.1449 [4] Michael V. Klibanov, Dinh-Liem Nguyen, Loc H. Nguyen, Hui Liu. A globally convergent numerical method for a 3D coefficient inverse problem with a single measurement of multi-frequency data. Inverse Problems & Imaging, 2018, 12 (2) : 493-523. doi: 10.3934/ipi.2018021 [5] Nikodem J. Poplawski, Abbas Shirinifard, Maciej Swat, James A. Glazier. Simulation of single-species bacterial-biofilm growth using the Glazier-Graner-Hogeweg model and the CompuCell3D modeling environment. Mathematical Biosciences & Engineering, 2008, 5 (2) : 355-388. doi: 10.3934/mbe.2008.5.355 [6] A. V. Fursikov. Stabilization for the 3D Navier-Stokes system by feedback boundary control. Discrete & Continuous Dynamical Systems - A, 2004, 10 (1&2) : 289-314. doi: 10.3934/dcds.2004.10.289 [7] Anna Kostianko, Sergey Zelik. Inertial manifolds for the 3D Cahn-Hilliard equations with periodic boundary conditions. Communications on Pure & Applied Analysis, 2015, 14 (5) : 2069-2094. doi: 10.3934/cpaa.2015.14.2069 [8] Mei Wang, Zilai Li, Zhenhua Guo. Global weak solution to 3D compressible flows with density-dependent viscosity and free boundary. Communications on Pure & Applied Analysis, 2017, 16 (1) : 1-24. doi: 10.3934/cpaa.2017001 [9] Fabrice Delbary, Kim Knudsen. Numerical nonlinear complex geometrical optics algorithm for the 3D Calderón problem. Inverse Problems & Imaging, 2014, 8 (4) : 991-1012. doi: 10.3934/ipi.2014.8.991 [10] Boyan Jonov, Thomas C. Sideris. Global and almost global existence of small solutions to a dissipative wave equation in 3D with nearly null nonlinear terms. Communications on Pure & Applied Analysis, 2015, 14 (4) : 1407-1442. doi: 10.3934/cpaa.2015.14.1407 [11] Yong Zhou. Remarks on regularities for the 3D MHD equations. Discrete & Continuous Dynamical Systems - A, 2005, 12 (5) : 881-886. doi: 10.3934/dcds.2005.12.881 [12] Hyeong-Ohk Bae, Bum Ja Jin. Estimates of the wake for the 3D Oseen equations. Discrete & Continuous Dynamical Systems - B, 2008, 10 (1) : 1-18. doi: 10.3934/dcdsb.2008.10.1 [13] Indranil SenGupta, Weisheng Jiang, Bo Sun, Maria Christina Mariani. Superradiance problem in a 3D annular domain. Conference Publications, 2011, 2011 (Special) : 1309-1318. doi: 10.3934/proc.2011.2011.1309 [14] Giovanny Guerrero, José Antonio Langa, Antonio Suárez. Biodiversity and vulnerability in a 3D mutualistic system. Discrete & Continuous Dynamical Systems - A, 2014, 34 (10) : 4107-4126. doi: 10.3934/dcds.2014.34.4107 [15] Niklas Hartung. Efficient resolution of metastatic tumor growth models by reformulation into integral equations. Discrete & Continuous Dynamical Systems - B, 2015, 20 (2) : 445-467. doi: 10.3934/dcdsb.2015.20.445 [16] 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 [17] Sadek Gala. A new regularity criterion for the 3D MHD equations in $R^3$. Communications on Pure & Applied Analysis, 2012, 11 (3) : 973-980. doi: 10.3934/cpaa.2012.11.973 [18] Eleonora Messina. Numerical simulation of a SIS epidemic model based on a nonlinear Volterra integral equation. Conference Publications, 2015, 2015 (special) : 826-834. doi: 10.3934/proc.2015.0826 [19] Luigi C. Berselli. An elementary approach to the 3D Navier-Stokes equations with Navier boundary conditions: Existence and uniqueness of various classes of solutions in the flat boundary case.. Discrete & Continuous Dynamical Systems - S, 2010, 3 (2) : 199-219. doi: 10.3934/dcdss.2010.3.199 [20] Jiahong Wu. Regularity results for weak solutions of the 3D MHD equations. Discrete & Continuous Dynamical Systems - A, 2004, 10 (1&2) : 543-556. doi: 10.3934/dcds.2004.10.543

2017 Impact Factor: 0.972