Numerical solution of a two dimensional tumour growth model with moving boundary

01/07/2020
by   Jerome Droniou, et al.
0

We consider a biphasic continuum model for avascular tumour growth in two spatial dimensions, in which a cell phase and a fluid phase follow conservation of mass and momentum. A limiting nutrient that follows a diffusion process controls the birth and death rate of the tumour cells. The cell volume fraction, cell velocity - fluid pressure system, and nutrient concentration are the model variables. A coupled system of a hyperbolic conservation law, a viscoelastic system, and a parabolic diffusion equation governs the dynamics of the model variables. The tumour boundary moves with the normal velocity of the outermost layer of cells, and this time-dependence is a challenge in designing and implementing a stable and fast numerical scheme. We recast the model into a form where the hyperbolic equation is defined on a fixed extended domain and retrieve the tumour boundary as the interface at which the cell volume fraction decreases below a threshold value. This procedure eliminates the need to track the tumour boundary explicitly and the computationally expensive re-meshing of the time-dependent domains. A numerical scheme based on finite volume methods for the hyperbolic conservation law, Lagrange P_2 - P_1 Taylor-Hood finite element method for the viscoelastic system, and mass-lumped finite element method for the parabolic equations is implemented in two spatial dimensions, and several cases are studied. We demonstrate the versatility of the numerical scheme in catering for irregular and asymmetric initial tumour geometries. When the nutrient diffusion equation is defined only in the tumour region, the model depicts growth in free suspension. On the contrary, when the nutrient diffusion equation is defined in a larger fixed domain, the model depicts tumour growth in a polymeric gel.

READ FULL TEXT

page 17

page 18

page 20

page 21

research
10/17/2019

Convergence analysis of a numerical scheme for a tumour growth model

We consider a one-spatial dimensional tumour growth model that consists ...
research
07/09/2023

Well posedness of fluid/solid mixture models for biofilm spread

Two phase solid-fluid mixture models are ubiquitous in biological applic...
research
10/05/2019

A Conservative Finite Element ALE Scheme for Mass-Conserving Reaction-Diffusion Equations on Evolving Two-Dimensional Domains

Mass-conservative reaction-diffusion systems have recently been proposed...
research
01/31/2023

Exploring numerical blow-up phenomena for the Keller-Segel-Navier-Stokes equations

The Keller-Segel-Navier-Stokes system governs chemotaxis in liquid envir...
research
09/18/2020

Point Forces in Elasticity Equation and Their Alternatives in Multi Dimensions

We consider several mathematical issues regarding models that simulate f...
research
07/10/2020

Conservation laws for free-boundary fluid layers

Time-dependent models of fluid motion in thin layers, subject to signed ...
research
01/27/2023

Laplacian regularized eikonal equation with Soner boundary condition on polyhedral meshes

In this paper, we propose a numerical algorithm based on a cell-centered...

Please sign up or login with your details

Forgot password? Click here to reset