A phase-field model for viscoelastic compressible tumor growth

Growing tumors generate and respond to stress in their local microenvironment. Tissue re-arrangements can relax these mechanical stresses and make the tissue more fluid-like. Further, intricate coupling between mechanotransduction and biochemical signaling leads to complex patterns of growth. To predict the outcomes of these nonlinear interactions, we develop a phase-field model to simulate tumors growing into a surrounding medium taking into account their elastic and viscous properties as well as their compressibilities. We couple continuum modeling of the viscoelastic mechanics to the concentration of a diffusible growth-promoting nutrient in a mass conservative way. Additionally, the phase-field method is a stable and flexible way to describe the dynamics of arbitrarily shaped tumors. We demonstrate convergence of the phase-field model to a sharp interface model in radially symmetric geometries and can observe progression to stationary tumors. However, our results show that these stationary symmetric tumors are subject to symmetry-breaking instabilities in 2D and 3D driven by two primary mechanisms: (i) elastic buckling instabilities due to differential growth induced by the nutrient gradient and (ii) instabilities generated by apoptosis-related volumetric loss. Further, tissue fluidity and compressibility can lead to changes in tumor topologies. Our modeling framework provides a robust methodology for investigating how tissue mechanics and growth factor signaling influence the progression and invasive potential of solid tumors.