A gradient Discretisation Method For Anisotropic Reaction Diffusion Models with applications to the dynamics of brain tumours

Yahya Alnashri,Hasan Alzubaidi

DOI: https://doi.org/10.48550/arXiv.2009.00139

2020-09-01

Abstract:A gradient discretisation method (GDM), Gradient schemes, Convergence analysis, Existence of weak solutions, Anisotropic reaction diffusion models, Dirichlet and Neumann boundary conditions, Non conforming finite element methods, Finite volume schemes, Hybrid mixed mimetic (HMM) method, Crouzeix--Raviart scheme, Brain tumour dynamics, Fractional anisotropy.

Numerical Analysis

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to develop a Gradient Discretisation Method (GDM) for anisotropic reaction - diffusion models and apply it to the study of brain tumor dynamics. Specifically, the author focuses on how to simulate and analyze the growth and diffusion processes of highly diffusive brain tumors (such as gliomas) in a heterogeneous brain environment by numerical methods.

### Main problem decomposition:

1. **Mathematical modelling and theoretical analysis**:
- Study anisotropic reaction - diffusion equations, whose general form is as follows:
\[
\frac{\partial \bar{c}(\mathbf{x}, t)}{\partial t}-\nabla\cdot(A(\mathbf{x})\nabla\bar{c}(\mathbf{x}, t)) = F(\bar{c}(\mathbf{x}, t)),\quad(\mathbf{x}, t)\in\Omega\times(0, T)
\]
where $\bar{c}(\mathbf{x}, t)$ represents the cancer cell density, $A(\mathbf{x})$ is the diffusion tensor, and $F(\bar{c}(\mathbf{x}, t))$ is the reaction term.
- Prove the existence of weak solutions and provide a unified convergence analysis framework applicable to multiple numerical methods (such as the finite - volume method, the mixed finite - element method, etc.).

2. **Application of numerical methods**:
- Apply non - conforming numerical methods (such as the Crouzeix - Raviart scheme and the Hybrid Mixed Mimetic (HMM) method) to solve the above - mentioned model and verify the effectiveness of these methods.
- Through numerical experiments, study the growth dynamics of brain tumors in a heterogeneous environment, especially use Fractional Anisotropy (FA) to measure their highly diffusive characteristics.

3. **Practical application background**:
- Glioma is a highly invasive brain tumor with a high recurrence rate, and it is difficult to be completely cured by existing treatment methods. Therefore, predicting the diffusion path of the tumor through mathematical modelling is helpful for guiding the formulation of treatment plans.
- Using Diffusion Tensor Imaging (DTI) data can better describe the influence of white matter fiber tracts on tumor diffusion, thereby more accurately simulating the non - uniform growth of the tumor.

### Summary:
This paper aims to provide a unified numerical analysis framework for anisotropic reaction - diffusion models by developing and applying the Gradient Discretisation Method, with particular attention to its application in brain tumor dynamics. This not only helps to understand the growth mechanism of tumors but also provides theoretical support for improving treatment strategies.