A Three-Dimensional Tumor Growth Model and Its Boundary Instability

Abstract

In this paper, we investigate the instability of growing tumors by employing both analytical and numerical techniques to validate previous results and extend the analytical findings presented in a prior study by Feng et al. (Z Angew Math Phys 74:107, 2023). Building upon the insights derived from the analytical reconstruction of key results in the aforementioned work in one dimension and two dimensions, we extend our analysis to three dimensions. Specifically, we focus on the determination of boundary instability using perturbation and asymptotic analysis along with spherical harmonics. Additionally, we have validated our analytical results in a two-dimensional (2D) framework by implementing the Alternating Direction Implicit (ADI) method. Our primary focus has been on ensuring that the numerical simulation of the propagation speed aligns accurately with the analytical findings. Furthermore, we have matched the simulated boundary stability with the analytical predictions derived from the evolution function, which will be defined in subsequent sections of our paper. This alignment is essential for accurately determining the stability or instability of tumor boundaries.

Keywords

Tumor growth / Boundary instability / Asymptotic analysis / Spherical harmonics / Bessel functions

References

[1] AraujoRP, McElwainDS. A history of the study of solid tumour growth: the contribution of mathematical modelling. Bull. Math. Biol., 2004, 6651039-1091.

[2] BaskaranA, LowengrubJS, WangC, WiseSM. Convergence analysis of a second order convex splitting scheme for the modified phase field crystal equation. SIAM J. Numer. Anal., 2013, 5152851-2873.

[3] BidanCM, WangFM, DunlopJW. A three-dimensional model for tissue deposition on complex surfaces. Comput. Methods Biomech. Biomed. Eng., 2013, 16101056-1070.

[4] BrowningAP, MaclarenOJ, BuenzliPR, LanaroM, AllenbyMC, WoodruffMA, SimpsonMJ. Model-based data analysis of tissue growth in thin 3D printed scaffolds. J. Theor. Biol., 2021, 528. 110852

[5] ByrneH, AlarconT, OwenM, WebbS, MainiP. Modelling aspects of cancer dynamics: a review. Philos. Trans. R. Soc. A: Math. Phys. Eng. Sci., 2006, 36418431563-1578.

[6] ChenW, WangC, WangX, WiseSM. Positivity-preserving, energy stable numerical schemes for the Cahn-Hilliard equation with logarithmic potential. J. Comput. Phys. X, 2019, 3100031

[7] CristiniV, LowengrubJ, NieQ. Nonlinear simulation of tumor growth. J. Math. Biol., 2003, 46: 191-224.

[8] FalcóC, CohenDJ, CarrilloJA, BakerRE. Quantifying tissue growth, shape and collision via continuum models and Bayesian inference. J. R. Soc. Interface, 2023, 2020420230184.

[9] Feng, Y., He, Q., Liu, J.-G., Zhou, Z.: Rigorous derivation of a Hele-Shaw type model and its non-symmetric traveling wave solution. arXiv:2404.16353 (2024)

[10] FengY, TangM, XuX, ZhouZ. Tumor boundary instability induced by nutrient consumption and supply. Z. Angew. Math. Phys., 2023, 743107.

[11] FriedmanA, ReitichF. Symmetry-breaking bifurcation of analytic solutions to free boundary problems: an application to a model of tumor growth. Trans. Am. Math. Soc., 2001, 35341587-1634.

[12] GreenspanHP. Models for the growth of solid tumor as a problem by diffusion. Appl. Math. Comput., 1972, 30: 215-222

[13] GuY, ShenJ. Bound preserving and energy dissipative schemes for porous medium equation. J. Comput. Phys., 2020, 410. 109378

[14] HeinrichMA, AlertR, LaChanceJM, ZajdelTJ, KošmrljA, CohenDJ. Size-dependent patterns of cell proliferation and migration in freely-expanding epithelia. Elife, 2020, 9: 58945.

[15] HeinrichMA, AlertR, WolfAE, KošmrljA, CohenDJ. Self-assembly of tessellated tissue sheets by expansion and collision. Nat. Commun., 2022, 1314026.

[16] JacobsM, KimI, TongJ. Tumor growth with nutrients: regularity and stability. Commun. Am. Math. Soc., 2023, 304166-208.

[17] KimI, PožárN. Porous medium equation to Hele-Shaw flow with general initial density. Trans. Am. Math. Soc., 2018, 3702873-909.

[18] KimI, PožárN, WoodhouseB. Singular limit of the porous medium equation with a drift. Adv. Math., 2019, 349: 682-732.

[19] KimIC, PerthameB, SouganidisPE. Free boundary problems for tumor growth: a viscosity solutions approach. Nonlinear Anal., 2016, 138: 207-228.

[20] KimIC, TongJ. Interface dynamics in a two-phase tumor growth model. Interfaces Free Boundaries, 2021, 232191-304.

[21] LiuJ-G, TangM, WangL, ZhouZ. An accurate front capturing scheme for tumor growth models with a free boundary limit. J. Comput. Phys., 2018, 364: 73-94.

[22] LiuJ-G, WangL, ZhouZ. Positivity-preserving and asymptotic preserving method for 2D Keller-Segal equations. Math. Comput., 2018, 873111165-1189.

[23] LowengrubJS, FrieboesHB, JinF, ChuangY-L, LiX, MacklinP, WiseSM, CristiniV. Nonlinear modelling of cancer: bridging the gap between cells and tumours. Nonlinearity, 2009, 2311.

[24] QianY, WangC, ZhouS. Convergence analysis on a structure-preserving numerical scheme for the Poisson-Nernst-Planck-Cahn-Hilliard system. CSIAM Trans. Appl. Math., 2023, 4: 345-380.

[25] RooseT, ChapmanSJ, MainiPK. Mathematical models of avascular tumor growth. SIAM Rev., 2007, 492179-208.

[26] Tong, J., Zhang, Y.P.: Convergence of free boundaries in the incompressible limit of tumor growth models. arXiv:2403.05804 (2024)

[27] VázquezJLThe Porous Medium Equation: Mathematical Theory, 2006OxfordClarendon Press.

[28] WitelskiTP, BowenM. ADI schemes for higher-order nonlinear diffusion equations. Appl. Numer. Math., 2003, 452/3331-351.

Funding

National Key R&D Program of China(2021YFA1001200)

National Science Foundation of China Youth program(12101278)

Kunshan Shuangchuang Talent Program(kssc202102066)

RIGHTS & PERMISSIONS

Shanghai University