2190 articles – 2574 Notices  [english version]
HAL : hal-00457731, version 4

Fiche détaillée  Récupérer au format
SIAM Journal on Scientific Computing 34, 1 (2012) 76-104
Versions disponibles :
New resolution strategy for multi-scale reaction waves using time operator splitting, space adaptive multiresolution and dedicated high order implicit/explicit time integrators
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
(31/01/2012)

In this paper, we tackle the numerical simulation of reaction-diffusion equations modeling multiscale reaction waves. This type of problems induces peculiar diffculties 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 which are spatially very localized. In a series of previous studies, the numerical analysis of operator splitting techniques has been conducted and such an approach has shown a great potential in the framework of reaction-diffusion and convection-diffusion-reaction systems. However, even if a firm theoretical background is available, an optimal strategy for high performance numerical simulation is still needed. In this paper, we introduce a new strategy for reaction-diffusion systems based on time operator splitting in the context of very localized and very stiff reaction fronts. It provides an optimal combination of adaptive spatial multiresolution, implicit resolution of reaction and explicit resolution of diffusion. The optimality is reached in terms of the choice of the operator splitting time step which, in the framework of self-similar reaction waves, allows a very good combination of the various dedicated solvers used in the proposed strategy. The computational effciency is then evaluated through the numerical simulation of configurations which were so far out of reach of standard methods in the field of nonlinear chemical dynamics for spiral waves and scroll waves as an illustration.
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/Analyse numérique

Mathématiques/Equations aux dérivées partielles

Sciences de l'ingénieur/Milieux fluides et réactifs
Reaction-diffusion – Adaptive multiresolution – Operator splitting – Multi-scale waves
Liste des fichiers attachés à ce document : 
PDF
Duarte_SISC081686.pdf(9.7 MB)