On Finite Volume Discretization of Infiltration Dynamics in Tumor Growth Models
This work addresses numerical accuracy issues for researchers simulating tumor growth with immune cell infiltration, but the contribution is incremental as it enhances existing finite volume methods.
The paper identifies that conventional finite volume methods fail to preserve constant solutions and incompressibility in tumor growth models with immune cell infiltration, and proposes a new framework with segregated flux computations that satisfies conservation laws, yielding more accurate results.
We address numerical challenges in solving hyperbolic free boundary problems described by spherically symmetric conservation laws that arise in the modeling of tumor growth due to immune cell infiltrations. In this work, we normalize the radial coordinate to transform the free boundary problem to a fixed boundary one, and utilize finite volume methods to discretize the resulting equations. We show that the conventional finite volume methods fail to preserve constant solutions and the incompressibility condition, and they typically lead to inaccurate, if not wrong, solutions even for very simple tests. These issues are addressed in a new finite volume framework with segregated flux computations that satisfy sufficient conditions for ensuring the so-called totality conservation law and the geometric conservation law. Classical first-order and second-order finite volume methods are enhanced in this framework. Their performance is assessed by various benchmark tests to show that the enhanced methods are able to preserve the incompressibility constraint and produce much more accurate results than the conventional ones.