A p-adaptive polytopal discontinuous Galerkin method for high-order approximation of brain electrophysiology

Keywords

Advanced Numerical Methods for Scientific Computing
Code:
22/2025
Title:
A p-adaptive polytopal discontinuous Galerkin method for high-order approximation of brain electrophysiology
Date:
Wednesday 30th April 2025
Author(s):
Leimer Saglio, C. B.; Pagani, S.; Antonietti P. F.
Download link:
Abstract:
Multiscale mathematical models have shown great promise in computational brain electrophysiology but are still hindered by high computational costs due to fast dynamics and complex brain geometries, requiring very fine spatio-temporal resolution. This paper introduces a novel p-adaptive discontinuous Galerkin method on polytopal grids (PolyDG) coupled with Crank–Nicolson time integration to approximate such models efficiently. The p-adaptive method enhances local accuracy via dynamic, element-wise polynomial refinement/de-refinement guided by a-posteriori error estimators. A novel clustering algorithm automatizes the selection of elements for adaptive updates, further improving efficiency. A wide set of numerical tests, including epileptic seizure simulations in a sagittal section of a human brain stem, demonstrate the method’s ability to reduce computational load while maintaining the accuracy of the numerical solution in capturing the dynamics of multiple wavefronts.