A <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" altimg="si31.svg"> <mml:mi>p</mml:mi> </mml:math> -adaptive polytopal discontinuous Galerkin method for high-order approximation of brain electrophysiology

• A novel p -adaptive algorithm for polytopal discontinuous Galerkin discretization of brain electrophysiology models. • A novel clustering algorithm to automatize the selection of elements for p -adaptive refinement/coarsening. • Extensive numerical assessment of the p -adaptive algorithm driven by different a-posteriori error estimators. • Computational modeling of brain electrophysiology during epileptic seizures within anatomically complex and heterogeneous sections of cortical geometries. 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.