A p-adaptive polytopal discontinuous Galerkin method for high-order approximation of brain electrophysiology
Fuente:
arXiv
Salvato in:
| Autori principali: | , , |
|---|---|
| Natura: | Preprint |
| Pubblicazione: |
2025
|
| Soggetti: | |
| Accesso online: | |
| Tags: |
Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
|
| _version_ | 1866910923183095808 |
|---|---|
| author | Saglio, Caterina Beatrice Leimer Pagani, Stefano Antonietti, Paola F. |
| author_facet | Saglio, Caterina Beatrice Leimer Pagani, Stefano Antonietti, Paola F. |
| contents | 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. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2504_21657 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | A p-adaptive polytopal discontinuous Galerkin method for high-order approximation of brain electrophysiology Saglio, Caterina Beatrice Leimer Pagani, Stefano Antonietti, Paola F. Numerical Analysis 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. |
| title | A p-adaptive polytopal discontinuous Galerkin method for high-order approximation of brain electrophysiology |
| topic | Numerical Analysis |
| url | https://arxiv.org/abs/2504.21657 |