A Cahn‐Hilliard–type equation with application to tumor growth dynamics

We consider a Cahn-Hilliard–type equation with degenerate mobility and single-well potential of Lennard-Jones type. This equation models the evolution and growth of biological cells such as solid tumors. The degeneracy set of the mobility and the singularity set of the cellular potential do not coincide, and the absence of cells is an unstable equilibrium configuration of the potential. This feature introduces a nontrivial difference with respect to the Cahn-Hilliard equation analyzed in the literature. We give existence results for different classes of weak solutions. Moreover, we formulate a continuous finite element approximation of the problem, where the positivity of the solution is enforced through a discrete variational inequality. We prove the existence and uniqueness of the discrete solution for any spatial dimension together with the convergence to the weak solution for spatial dimension d=1. We present simulation results in 1 and 2 space dimensions. We also study the dynamics of the spinodal decomposition and the growth and scaling laws of phase ordering dynamics. In this case, we find similar results to the ones obtained in standard phase ordering dynamics and we highlight the fact that the asymptotic behavior of the solution is dominated by the mechanism of growth by bulk diffusion.

[1]  M. Gurtin,et al.  TWO-PHASE BINARY FLUIDS AND IMMISCIBLE FLUIDS DESCRIBED BY AN ORDER PARAMETER , 1995, patt-sol/9506001.

[2]  Charles M. Elliott,et al.  Numerical analysis of the Cahn-Hilliard equation with a logarithmic free energy , 1992 .

[3]  Jaemin Shin,et al.  Physical, mathematical, and numerical derivations of the Cahn–Hilliard equation , 2014 .

[4]  T. Hughes,et al.  Isogeometric analysis of the Cahn–Hilliard phase-field model , 2008 .

[5]  Ju Liu,et al.  Isogeometric analysis of the advective Cahn-Hilliard equation: Spinodal decomposition under shear flow , 2013, J. Comput. Phys..

[6]  F. Otto,et al.  Upper Bounds on Coarsening Rates , 2002 .

[7]  Krishna Garikipati,et al.  A discontinuous Galerkin method for the Cahn-Hilliard equation , 2006, J. Comput. Phys..

[8]  Johan van de Koppel,et al.  Phase separation explains a new class of self-organized spatial patterns in ecological systems , 2013, Proceedings of the National Academy of Sciences.

[9]  Charles M. Elliott,et al.  The Cahn–Hilliard equation with a concentration dependent mobility: motion by minus the Laplacian of the mean curvature , 1996, European Journal of Applied Mathematics.

[10]  N. Kenmochi,et al.  Subdifferential Operator Approach to the Cahn-Hilliard Equation with Constraint , 1995 .

[11]  Martin Rumpf,et al.  Nonnegativity preserving convergent schemes for the thin film equation , 2000, Numerische Mathematik.

[12]  G. Schimperna,et al.  Existence of solutions and separation from singularities for a class of fourth order degenerate parabolic equations , 2010, 1009.3148.

[13]  J. E. Hilliard,et al.  Free Energy of a Nonuniform System. I. Interfacial Free Energy , 1958 .

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

[15]  A. Karma,et al.  Evolution of nanoporosity in dealloying , 2001, Nature.

[16]  John W. Cahn,et al.  On Spinodal Decomposition , 1961 .

[17]  C. M. Elliott,et al.  On the Cahn-Hilliard equation with degenerate mobility , 1996 .

[18]  Sergey Zelik,et al.  The Cahn-Hilliard Equation with Logarithmic Potentials , 2011 .

[19]  Charles M. Elliott,et al.  On the Cahn-Hilliard equation , 1986 .

[20]  A. Bray Theory of phase-ordering kinetics , 1994, cond-mat/9501089.

[21]  J. Lowengrub,et al.  Quasi–incompressible Cahn–Hilliard fluids and topological transitions , 1998, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[22]  Harald Garcke,et al.  Finite Element Approximation of the Cahn-Hilliard Equation with Degenerate Mobility , 1999, SIAM J. Numer. Anal..

[23]  G. Grün Degenerate Parabolic Differential Equations of Fourth Order and a Plasticity Model with Non-Local Hardening , 1995 .

[24]  S. Tremaine,et al.  On the Origin of Irregular Structure in Saturn's Rings , 2002, astro-ph/0211149.

[25]  Simone Scacchi,et al.  A C1 Virtual Element Method for the Cahn-Hilliard Equation with Polygonal Meshes , 2015, SIAM J. Numer. Anal..

[26]  A. Friedman,et al.  Higher order nonlinear degenerate parabolic equations , 1990 .

[27]  Leonard M Sander,et al.  Generalized Cahn-Hilliard equation for biological applications. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[28]  Endre Süli,et al.  Discontinuous Galerkin Finite Element Approximation of the Cahn-Hilliard Equation with Convection , 2009, SIAM J. Numer. Anal..

[29]  L. Mcmaster Aspects of Liquid-Liquid Phase Transition Phenomena in Multicomponent Polymeric Systems , 1975 .

[30]  Andrea L. Bertozzi,et al.  Inpainting of Binary Images Using the Cahn–Hilliard Equation , 2007, IEEE Transactions on Image Processing.

[31]  Yan Xu,et al.  Local discontinuous Galerkin methods for the Cahn-Hilliard type equations , 2007, J. Comput. Phys..

[32]  Martin Z. Bazant,et al.  Phase Separation Dynamics in Isotropic Ion-Intercalation Particles , 2013, SIAM J. Appl. Math..

[33]  Clément Chatelain,et al.  Morphological changes in early melanoma development: influence of nutrients, growth inhibitors and cell-adhesion mechanisms. , 2011, Journal of theoretical biology.

[34]  Carlos J. García-Cervera,et al.  Ordering kinetics of a conserved binary mixture with a nematic liquid crystal component , 2014 .

[35]  M. Gurtin Generalized Ginzburg-Landau and Cahn-Hilliard equations based on a microforce balance , 1996 .

[36]  H. Frieboes,et al.  Three-dimensional multispecies nonlinear tumor growth--I Model and numerical method. , 2008, Journal of theoretical biology.

[37]  Amy Novick-Cohen,et al.  Upper bounds for coarsening for the degenerate Cahn-Hilliard equation , 2009 .

[38]  Franck Pigeonneau,et al.  A Hybrid High-Order Method for the Cahn-Hilliard problem in Mixed Form , 2015, SIAM J. Numer. Anal..