Simulating complex tumor dynamics from avascular to vascular growth using a general level-set method

A comprehensive continuum model of solid tumor evolution and development is investigated in detail numerically, both under the assumption of spherical symmetry and for arbitrary two-dimensional growth. The level set approach is used to obtain solutions for a recently developed multi-cell transport model formulated as a moving boundary problem for the evolution of the tumor. The model represents both the avascular and the vascular phase of growth, and is able to simulate when the transition occurs; progressive formation of a necrotic core and a rim structure in the tumor during the avascular phase are also captured. In terms of transport processes, the interaction of the tumor with the surrounding tissue is realistically incorporated. The two-dimensional simulation results are presented for different initial configurations. The computational framework, based on a Cartesian mesh/narrow band level-set method, can be applied to similar models that require the solution of coupled advection-diffusion equations with a moving boundary inside a fixed domain. The solution algorithm is designed so that extension to three-dimensional simulations is straightforward.

[1]  R. D. Richtmyer,et al.  Difference methods for initial-value problems , 1959 .

[2]  G. Hedstrom,et al.  Numerical Solution of Partial Differential Equations , 1966 .

[3]  R. D. Richtmyer,et al.  Difference methods for initial-value problems , 1959 .

[4]  N. N. Yanenko,et al.  The Method of Fractional Steps , 1971 .

[5]  J. Folkman The vascularization of tumors. , 1976, Scientific American.

[6]  Burton I. Edelson,et al.  Global satellite communications , 1977 .

[7]  G. Smith,et al.  Numerical Solution of Partial Differential Equations: Finite Difference Methods , 1978 .

[8]  J. Sethian Curvature and the evolution of fronts , 1985 .

[9]  James A. Sethian,et al.  Numerical Methods for Propagating Fronts , 1987 .

[10]  P. Concus,et al.  Variational Methods for Free Surface Interfaces , 1987 .

[11]  Chia-Ven Pao,et al.  Nonlinear parabolic and elliptic equations , 1993 .

[12]  D. Chopp Computing Minimal Surfaces via Level Set Curvature Flow , 1993 .

[13]  R. LeVeque,et al.  A comparison of the extended finite element method with the immersed interface method for elliptic equations with discontinuous coefficients and singular sources , 2006 .

[14]  J. Sethian,et al.  A Fast Level Set Method for Propagating Interfaces , 1995 .

[15]  J. W. Thomas Numerical Partial Differential Equations: Finite Difference Methods , 1995 .

[16]  Wei Shyy,et al.  Computational Fluid Dynamics with Moving Boundaries , 1995 .

[17]  Baba C. Vemuri,et al.  Shape Modeling with Front Propagation: A Level Set Approach , 1995, IEEE Trans. Pattern Anal. Mach. Intell..

[18]  M. Chaplain,et al.  Modelling the role of cell-cell adhesion in the growth and development of carcinomas , 1996 .

[19]  J A Sethian,et al.  A fast marching level set method for monotonically advancing fronts. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[20]  Stanley Osher,et al.  A Hybrid Method for Moving Interface Problems with Application to the Hele-Shaw Flow , 1997 .

[21]  S. Osher,et al.  A Simple Level Set Method for Solving Stefan Problems , 1997, Journal of Computational Physics.

[22]  John A. Adam,et al.  General Aspects of Modeling Tumor Growth and Immune Response , 1997 .

[23]  James A. Sethian,et al.  The Fast Construction of Extension Velocities in Level Set Methods , 1999 .

[24]  J. A. Sethian,et al.  Fast Marching Methods , 1999, SIAM Rev..

[25]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods/ J. A. Sethian , 1999 .

[26]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods , 1999 .

[27]  L. Preziosi,et al.  Modelling and mathematical problems related to tumor evolution and its interaction with the immune system , 2000 .

[28]  L. Preziosi,et al.  ADVECTION-DIFFUSION MODELS FOR SOLID TUMOUR EVOLUTION IN VIVO AND RELATED FREE BOUNDARY PROBLEM , 2000 .

[29]  S Torquato,et al.  Simulated brain tumor growth dynamics using a three-dimensional cellular automaton. , 2000, Journal of theoretical biology.

[30]  Dantzig,et al.  Computation of dendritic microstructures using a level set method , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[31]  J. Sethian Evolution, implementation, and application of level set and fast marching methods for advancing fronts , 2001 .

[32]  D. Calhoun A Cartesian Grid Method for Solving the Two-Dimensional Streamfunction-Vorticity Equations in Irregular Regions , 2002 .

[33]  Luigi Preziosi,et al.  Cancer Modelling and Simulation , 2003 .

[34]  L. Preziosi,et al.  Modelling Solid Tumor Growth Using the Theory of Mixtures , 2001, Mathematical medicine and biology : a journal of the IMA.

[35]  Helen M. Byrne,et al.  A two-phase model of solid tumour growth , 2003, Appl. Math. Lett..

[36]  Marko Subasic,et al.  Level Set Methods and Fast Marching Methods , 2003 .

[37]  Ronald Fedkiw,et al.  Level set methods and dynamic implicit surfaces , 2002, Applied mathematical sciences.

[38]  V. Cristini,et al.  Nonlinear simulation of tumor growth , 2003, Journal of mathematical biology.

[39]  J. Sethian,et al.  FRONTS PROPAGATING WITH CURVATURE DEPENDENT SPEED: ALGORITHMS BASED ON HAMILTON-JACOB1 FORMULATIONS , 2003 .

[40]  Thomas S. Deisboeck,et al.  Simulating ‘structure–function’ patterns of malignant brain tumors , 2004 .

[41]  D L S McElwain,et al.  A history of the study of solid tumour growth: The contribution of mathematical modelling , 2004, Bulletin of mathematical biology.

[42]  Cosmina S. Hogeaa,et al.  Computational modeling of solid tumor evolution via a general Cartesian mesh / level set method , 2004 .

[43]  V. Cristini,et al.  Nonlinear simulation of tumor necrosis, neo-vascularization and tissue invasion via an adaptive finite-element/level-set method , 2005, Bulletin of mathematical biology.

[44]  J. Sethian,et al.  Implementation of the level set method for continuum mechanics based tumor growth models , 2005 .

[45]  Vasilios Alexiades,et al.  OVERCOMING THE STABILITY RESTRICTION OF EXPLICIT SCHEMES VIA SUPER-TIME-STEPPING , 2022 .