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

Abstract

A gradient discretisation method (GDM) is an abstract setting that designs the unified convergence analysis of several numerical methods for partial differential equations and their corresponding models. In this paper, we study the GDM for anisotropic reaction diffusion problems, based on a general reaction term, with Neumann and Dirichlet boundary conditions. With natural regularity assumptions on the exact solution, the framework enables us to provide proof of the existence of weak solutions for the problem, and to obtain a uniform-in-time convergence for the discrete solution and a strong convergence for its discrete gradient. It also allows us to apply non conforming numerical schemes to the model on a generic grid; (the Crouzeix-Raviart scheme and the hybrid mixed mimetic (HMM) methods). Numerical experiments using the HMM method are performed to study the growth of glioma tumours in heterogeneous brain environment. The dynamics of their highly diffusive nature is also measured using the fraction anisotropic measure. The validity of the HMM is examined further using four different mesh types. The results indicate that the dynamics of the brain tumour is still captured by the HMM scheme, even in the event of a highly heterogeneous anisotropic case performed on the mesh with extreme distortions.