Design and analysis of a lightweight parallel adaptive scheme for the solution of the monodomain equation

Abstract : Numerical simulation of the nonlinear reaction-diffusion equations in computational electrocardiology requires locally high spatial resolution to capture the multiscale effects related to the electrical activation of the heart accurately, namely the strongly varying transmembrane potential. Here, we propose a novel lightweight adaptive algorithm which aims at ombining the plainness of structured meshes with the resolving capabilities of unstructered adaptive meshes. Our "patchwise adaptive" approach is based on locally structured mesh hierarchies which are glued along their interfaces by a nonconforming mortar element discretization. To further increase the overall efficiency, we keep the spatial meshes constant over suitable time windows in which error indicators are accumulated. This approach facilitates strongly varying mesh sizes in neighboring patches as well as in consecutive time steps. For the transfer of the dynamic variables between different spatial approximation spaces we compare the L2-projection and a local approximation. Finally, since an implicit-explicit time discretization is employed for stability reasons, we derive a spatial preconditioner which is tailored to the special structure of the patchwise adaptive meshes. We analyze the (parallel) performance and scalability of the resulting method by several examples from computational electrocardiology of different sizes. Additionally, we compare our method to a standard adaptive refinement strategy using unstructured meshes. As it turns out, our novel adaptive scheme provides a very good balance between reduction in degrees of freedom and overall (parallel) efficiency.
Type de document :
Article dans une revue
SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2014, 36 (2), pp.C163-C189. 〈10.1137/130912505〉
Liste complète des métadonnées

https://hal.inria.fr/hal-01024658
Contributeur : Mark Potse <>
Soumis le : mercredi 16 juillet 2014 - 14:55:08
Dernière modification le : jeudi 11 janvier 2018 - 06:23:41

Identifiants

Collections

Citation

Thomas Dickopf, Dorian Krause, Rolf Krause, Mark Potse. Design and analysis of a lightweight parallel adaptive scheme for the solution of the monodomain equation. SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2014, 36 (2), pp.C163-C189. 〈10.1137/130912505〉. 〈hal-01024658〉

Partager

Métriques

Consultations de la notice

250