TY - JOUR
T1 - Numerical study of cancer cell invasion dynamics using adaptive mesh refinement: the urokinase model
AU - Kolbe, Niklas
AU - Katuchova, Jana
AU - Sfakianakis, Nikolaos
AU - Hellmann, Nadja
AU - Lukacova-Medvidova, Maria
PY - 2014/8/1
Y1 - 2014/8/1
N2 - In the present work we investigate the chemotactically and
proteolytically driven tissue invasion by cancer cells. The model
employed is a system of advection-reaction-diffusion equations that
features the role of the serine protease urokinase-type plasminogen
activator. The analytical and numerical study of this system constitutes
a challenge due to the merging, emerging, and travelling concentrations
that the solutions exhibit. Classical numerical methods applied to this
system necessitate very fine discretization grids to resolve these
dynamics in an accurate way. To reduce the computational cost without
sacrificing the accuracy of the solution, we apply adaptive mesh
refinement techniques, in particular h-refinement. Extended numerical
experiments exhibit that this approach provides with a higher order,
stable, and robust numerical method for this system. We elaborate on
several mesh refinement criteria and compare the results with the ones
in the literature. We prove, for a simpler version of this model,
$L^\infty$ bounds for the solutions, we study the stability of its
conditional steady states, and conclude that it can serve as a test case
for further development of mesh refinement techniques for cancer
invasion simulations.
AB - In the present work we investigate the chemotactically and
proteolytically driven tissue invasion by cancer cells. The model
employed is a system of advection-reaction-diffusion equations that
features the role of the serine protease urokinase-type plasminogen
activator. The analytical and numerical study of this system constitutes
a challenge due to the merging, emerging, and travelling concentrations
that the solutions exhibit. Classical numerical methods applied to this
system necessitate very fine discretization grids to resolve these
dynamics in an accurate way. To reduce the computational cost without
sacrificing the accuracy of the solution, we apply adaptive mesh
refinement techniques, in particular h-refinement. Extended numerical
experiments exhibit that this approach provides with a higher order,
stable, and robust numerical method for this system. We elaborate on
several mesh refinement criteria and compare the results with the ones
in the literature. We prove, for a simpler version of this model,
$L^\infty$ bounds for the solutions, we study the stability of its
conditional steady states, and conclude that it can serve as a test case
for further development of mesh refinement techniques for cancer
invasion simulations.
KW - Mathematics - Numerical Analysis
KW - Quantitative Biology - Cell Behavior
KW - 92B05
KW - 35Q92
KW - 65M08
KW - 65M50
M3 - Article
SN - 0096-3003
VL - 273
JO - Applied Mathematics and Computation
JF - Applied Mathematics and Computation
ER -