2185 articles – 2572 Notices  [english version]
HAL : hal-00590812, version 1

Fiche détaillée  Récupérer au format
Summer school on multiresolution and adaptive mesh refinement methods, E. Cancès, V. Louvet, M. Massot (Ed.) (2011) Vol. 34 277-290
New Resolution Strategy for Multi-scale Reaction Waves using Time Operator Splitting and Space Adaptive Multiresolution: Application to Human Ischemic Stroke
Max Duarte ( ) 1, Marc Massot 1, Stéphane Descombes 2, Christian Tenaud 3, Thierry Dumont 4, 5, Violaine Louvet 4, 5, Frédérique Laurent 1
PEPS Maths-ST2I project MIPAC Collaboration(s)
(22/12/2011)

We tackle the numerical simulation of reaction-diffusion equations modeling multi-scale reaction waves. This type of problems induces peculiar difficulties and potentially large stiffness which stem from the broad spectrum of temporal scales in the nonlinear chemical source term as well as from the presence of large spatial gradients in the reactive fronts, spatially very localized. In this paper, we introduce a new resolution strategy based on time operator splitting and space adaptive multiresolution in the context of very localized and stiff reaction fronts. Based on recent theoretical studies of numerical analysis, such a strategy leads to a splitting time step which is not restricted neither by the fastest scales in the source term nor by restrictive diffusive step stability limits, but only by the physics of the phenomenon. We aim thus at solving accurately complete models including all time and space scales of the phenomenon, considering large simulation domains with conventional computing resources. The efficiency of the method is evaluated through 2D and 3D numerical simulations of a human ischemic stroke model, conducted on a simplified brain geometry, for which a simple parallelization strategy for shared memory architectures was implemented, in order to reduce computing costs related to "detailed chemistry" features of the model.
1 :  Laboratoire d'Énergétique Moléculaire et Macroscopique, Combustion (EM2C)
CNRS : UPR288 – Ecole Centrale Paris
2 :  Laboratoire Jean Alexandre Dieudonné (JAD)
CNRS : UMR6621 – Université Nice Sophia Antipolis [UNS]
3 :  Laboratoire d'Informatique pour la Mécanique et les Sciences de l'Ingénieur [Orsay] (LIMSI)
CNRS : UPR3251 – Université Pierre et Marie Curie [UPMC] - Paris VI – Université Paris XI - Paris Sud
4 :  Institut Camille Jordan (ICJ)
CNRS : UMR5208 – Université Claude Bernard - Lyon I – Ecole Centrale de Lyon – Institut National des Sciences Appliquées (INSA) - Lyon
5 :  NUMED (ENS Lyon / UCB Lyon / Inria Grenoble Rhône-Alpes)
Université Claude Bernard - Lyon I – INRIA – École Normale Supérieure - Lyon – CNRS : UMR5669 – Unité de Mathématiques Pures et Appliquées
Mathématiques/Equations aux dérivées partielles

Informatique/Modélisation et simulation

Informatique/Calcul parallèle, distribué et partagé
Reaction-diffusion – Operator splitting – Adaptive multiresolution – Ischemic stroke – Parallel computing
Liste des fichiers attachés à ce document : 
PDF
Duarte_etal_ESAIM.pdf(381.3 KB)