Discontinuous Galerkin Methods for Fisher-Kolmogorov Equation with Application to $\alpha$-Synuclein Spreading in Parkinson's Disease

Abstract

This spreading of prion proteins is at the basis of brain neurodegeneration. This paper deals with the numerical modelling of the misfolding process of $\alpha$-synuclein in Parkinson's disease. We introduce and analyze a discontinuous Galerkin method for the semi-discrete approximation of the Fisher-Kolmogorov (FK) equation that can be employed to model the process. We employ a discontinuous Galerkin method on polygonal and polyhedral grids (PolyDG) for space discretization, to accurately simulate the wavefronts typically observed in the prionic spreading and we prove stability and a priori error estimates. Next, we use a Crank-Nicolson scheme to advance in time. For the numerical verification of our numerical model, we first consider a manufactured solution, and then we consider a case with wavefront propagation in two-dimensional polygonal grids. Next, we carry out a simulation of $\alpha$-synuclein spreading in a two-dimensional brain slice in the sagittal plane with a polygonal agglomerated grid that takes full advantage of the flexibility of PolyDG approximation. Finally, we present a simulation in a three-dimensional geometry reconstructed from magnetic resonance images of a patient's brain.

arXiv admin note: text overlap with arXiv:2210.02272

References in corpus (1)

• Numerical Modelling of the Brain Poromechanics by High-Order Discontinuous Galerkin Methods

Cited by in corpus (9)

• Structure Preserving Polytopal Discontinuous Galerkin Methods for the Numerical Modeling of Neurodegenerative Diseases
• Exploring tau protein and amyloid-beta propagation: a sensitivity analysis of mathematical models based on biological data
• Discontinuous Galerkin approximations of the heterodimer model for protein-protein interaction
• Uncertainty Quantification for Fisher-Kolmogorov Equation on Graphs with Application to Patient-Specific Alzheimer Disease
• Numerical modelling of protein misfolding in neurodegenerative diseases: a computational study
• A coupled mathematical and numerical model for protein spreading and tissue atrophy, applied to Alzheimer's disease
• A polytopal discontinuous Galerkin method for the pseudo-stress formulation of the unsteady Stokes problem
• A structure-preserving LDG discretization of the Fisher-Kolmogorov equation for modeling neurodegenerative diseases
• A whole-brain model of amyloid beta accumulation and cerebral hypoperfusion in Alzheimer's disease