Structure-preserving upwind DG scheme for a Cahn-Hilliard-Darcy model of tumor growth
Provides a provably stable and accurate numerical method for simulating tumor growth in fluid-saturated porous media, relevant for biomedical modeling.
Developed a structure-preserving numerical scheme for a Cahn-Hilliard-Darcy tumor growth model that guarantees mass conservation, pointwise bounds, and energy dissipation. Numerical experiments demonstrate robustness and show fluid effects on tumor evolution.
In this work, we develop a structure-preserving numerical scheme for a Cahn-Hilliard-Darcy model that describes tumor growth in a fluid-saturated porous medium. First, we derive a physically consistent model from the general framework proposed in [29] that guarantees mass conservation and pointwise bounds on the phase-field and nutrient variables, with a decreasing energy law. The resulting model couples the evolution of tumor cells via a Cahn-Hilliard equation with a diffusion equation for the nutrients thro chemotactic interactions and extends the model in [1] by introducing the effect of a surrounding fluid described by Darcy's law. Subsequently, we propose a fully discrete scheme that combines an upwind discontinuous Galerkin method in space and a convex splitting strategy in time, which inherits the fundamental properties of the continuous model: mass conservation, pointwise bounds and discrete energy law. Our theoretical analysis is accompanied by numerical experiments that demonstrate the robustness of the proposed scheme and show the influence of the surrounding fluid on the tumor evolution.