Estimating the extent of glioblastoma invasion: Approximate stationalization of anisotropic advection–diffusion–reaction equations in the context of glioblastoma invasion

Glioblastoma Multiforme is a malignant brain tumor with poor prognosis. There have been numerous attempts to model the invasion of tumorous glioma cells via partial differential equations in the form of advection–diffusion–reaction equations. The patient-wise parametrization of these models, and their validation via experimental data has been found to be difficult, as time sequence measurements are mostly missing. Also the clinical interest lies in the actual (invisible) tumor extent for a particular MRI/DTI scan and not in a predictive estimate. Therefore we propose a stationalized approach to estimate the extent of glioblastoma (GBM) invasion at the time of a given MRI/DTI scan. The underlying dynamics can be derived from an instationary GBM model, falling into the wide class of advection-diffusion-reaction equations. The stationalization is introduced via an analytic solution of the Fisher-KPP equation, the simplest model in the considered model class. We investigate the applicability in 1D and 2D, in the presence of inhomogeneous diffusion coefficients and on a real 3D DTI-dataset.