2024/01/25 by Caterina Beatrice Leimer Saglio, Stefano Pagani, Saglio, Caterina Beatrice Leimer +5 · 2 citations
Chemistry · Computer Science · Physics and Astronomy · #Advanced Mathematical Modeling in Engineering #Advanced NMR Techniques and Applications #FOS: Mathematics #NMR spectroscopy and applications #Numerical Analysis (math.NA)
paper · pdf · doi:10.48550/arxiv.2401.14310
openalex publication_date 2024/01/25 · openalex created_date 2025/10/10 · openalex updated_date 2026/08/01
Epilepsy is a clinical neurological disorder characterized by recurrent and spontaneous seizures consisting of abnormal high-frequency electrical activity in the brain. In this condition, the transmembrane potential dynamics are characterized by rapid and sharp wavefronts traveling along the heterogeneous and anisotropic conduction pathways of the brain. This work employs the monodomain model, coupled with specific neuronal ionic models characterizing ion concentration dynamics, to mathematically describe brain tissue electrophysiology in grey and white matter at the organ scale. This multiscale model is discretized in space with the high-order discontinuous Galerkin method on polygonal and polyhedral grids (PolyDG) and advanced in time with a Crank-Nicolson scheme. This ensures, on the one hand, efficient and accurate simulations of the high-frequency electrical activity that is responsible for epileptic seizure and, on the other hand, keeps reasonably low the computational costs by a suitable combination of high-order approximations and agglomerated polytopal meshes. We numerically investigate synthetic test cases on a two-dimensional heterogeneous squared domain discretized with a polygonal grid, and on a two-dimensional brainstem in a sagittal plane with an agglomerated polygonal grid that takes full advantage of the flexibility of the PolyDG approximation of the semidiscrete formulation. Finally, we provide a theoretical analysis of stability and an a-priori convergence analysis for a simplified mathematical problem.