Work

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

Year: 2025

Type: article

Abstract: Multiscale mathematical models have shown significant potential in computational brain electrophysiology. However, their practical implementation is still limited by the substantial computational costs associated with the brain’s rapid dynamics and complex geometries, which require exceedingly fine spatio-temporal resolution. In this paper, we propose a novel p -adaptive discontinuous Galerkin method on polytopal grids (PolyDG) coupled with Crank–Nicolson time stepping for the numerical discretization of a brain electrophysiology model consisting of the monodomain equation coupled with the Barreto–Cressman ionic model. The proposed p -adaptive strategy enhances local accuracy through dynamic, element-wise polynomial refinement and coarsening, guided by a-posteriori error estimators. To further enhance computational efficiency, we introduce a novel clustering algorithm that automatically and dynamically identifies the subset of mesh elements where p -adaptive updates are required. Comprehensive numerical experiments, including benchmark test cases and simulations of epileptic seizure activity in a sagittal section of the human brainstem, demonstrate the method’s ability to significantly reduce the global number of degrees of freedom while maintaining the accuracy necessary to resolve complex wavefront dynamics.

Source: Computer Methods in Applied Mechanics and Engineering

DOI: 10.1016/j.cma.2025.118249

Authors: Caterina B․ Leimer Saglio, Stefano Pagani, Paola F․ Antonietti

Institution: Politecnico di Milano

Language: English

Topic: Advanced Numerical Methods in Computational Mathematics

Subfield: Computational Mechanics

Field: Engineering

Domain: Physical Sciences

Open Access status: hybrid

Estimated APC paid: $4,830