Le modèle hydrodynamique MFIX

Les éruptions volcaniques sont des événements parfois dévastateurs dont il convient d’anticiper les phases paroxysmales afin d’en limiter les conséquences sur les populations ; la surveillance continue s’impose comme une composante majeure de gestion des risques. L’identification rapide des événements précurseurs et la prédiction à court terme d’une manière fiable constituent un ample domaine de recherche dans le but d’améliorer les systèmes d’alerte. De nombreux travaux ont démontré l’importance de l’écoute microsismique comme paramètre de caractérisation des mécanismes précurseurs aux éruptions, essentiellement explosives. L’interprétation des signatures des signaux acoustiques fournit également de précieuses indications sur le caractère effusif ou explosif d’une activité en augmentation. Cependant, la compréhension des complexes procédés de mécanique des fluides qui accompagnent les éruptions représentent un troisième axe de recherche que nous proposons d’explorer ici.

Les éléments volatiles saturant le magma sont plus ou moins solubles selon les conditions de pression et de température auxquelles ils sont soumis. Lorsque les conditions d’équilibre sont rompues, l’exsolution de ces volatiles sous forme de bulles résulte en la migration du magma vers la surface. Au cours de cette ascension, la croissance des bulles s’accélère et leur solubilité diminue, entrainant une diffusion du gaz limitée essentiellement par la viscosité du liquide silicaté. Cette vitesse d’ascension dépend donc de la capacité des bulles de gaz à croître au sein du magma. L’augmentation de la flottabilité résultant de la croissance de ces bulles fournit la force motrice de l’éruption ; une énergie considérable peut s’accumuler lorsque la surface refroidie se comporte comme une barrière et piège des gaz dont la pression peut se libérer violemment.

Les magmas montant lentement vers la surface dégazent d’une manière plutôt effusive alors que les magmas rapides entrent en éruption de manière plus explosive, souvent plus dangereuse pour les populations. Au cours de cette étude, nous avons utilisé la modélisation numérique de processus de mécanique des fluides dans le but d’améliorer la compréhension du régime effusif dont les dynamiques qui le gouverne sont encore insuffisamment connues. Le volcan Erebus, en Antarctique, présente un lac de cratère permanent dont l’activité est essentiellement effusive, accompagnée d’explosions stromboliennes de faible ampleur. Il représente donc une précieuse source de données que les scientifiques ont collectées depuis des décennies, véritable fenêtre ouverte sur le cœur d’un système magmatique.

L’idée centrale d’une approche multiphasique est basée sur la mécanique des milieux continus. Celle-ci considère qu’un système se compose d’un nombre suffisamment important de particules pour que les discontinuités puissent être aplanies et que ses différentes propriétés puissent être traitées comme étant continues . Ainsi, tout système peut être divisé en petites cellules appelées volumes finis , chacune d’elles étant orientée dans l’espace par rapport à chaque cellule voisine dans un système de coordonnées . Tout au long de cette étude, nous adoptons un signe positif de convention : les contraintes sont positives lorsque les forces extérieure et normale agissent dans la même direction par rapport au système de coordonnées x, y correspondant à E et vers le haut, respectivement. Chaque volume fini ou volume de référence  est composé de 1 à 2 phases   ; celle qui se présente sous forme de grains ou composés en petite quantité est appelée la phase dispersée ou granulaire (e.g. les cristaux) et la phase dans laquelle les grains sont dispersés est appelée la phase continue. Chaque phase granulaire est caractérisée par un diamètre et une densité uniques, et un coefficient d’élasticité.

Notre modélisation est caractéristique de modèles appelés continus dans lesquels les variables d’état d’une composante changent continuellement avec le temps (description Eulérienne) ; ceux-ci se différencient des modèles dit discrets dans lesquels les variables d’état d’une composante changent instantanément à des temps distincts (description Lagrangienne).

Les équations vectorielles du modèle hydrodynamique sont basées sur ces équations et sur plusieurs travaux développés depuis 1967. Pour résoudre les équations d’écoulement multiphasique, nous avons utilisé un code numérique basé sur MFIX Ver. 2.0, 2004 (Multiphase Flow with Interface eXchanges), un programme développé par le Laboratoire National Américain de Technologie Energétique du Département de l’Energie [Syamlal et al., 1993 ; Syamlal, 1994, 1998 ; Benyahia et al., 2008]. Ce programme permet de simuler la dynamique des fluides dont la complexité est telle qu’elle ne peut se résoudre à travers des méthodes numériques sans l’aide d’un ordinateur ; c’est-à-dire en utilisant la dynamique des fluides computationnelle (CFD). MFIX a la capacité de traiter jusqu’à 4 phases dispersées. Ces équations à dérivée partielle sont linéarisées et résolues à l’aide de techniques de ésolution linéaire   pour l’ensemble des cellules à un temps donné. Un compterendu plus détaillé des équations utilisées et des méthodes de discrétisation peut être trouvé dans la section du manuel de MFIX relative à la théorie et à la technique numériques [Syamlal, 1993; 1994].

Tenseur de contrainte de la phase granulaire

Nous nous sommes d’abord concentrés sur l’examen de la rhéologie globale de l’écoulement granulaire en utilisant les relations constitutives existantes dans MFIX, car elles constituent probablement l’approche la plus acceptable aujourd’hui pour modéliser les écoulements de particules [e.g., Valentine, 1994 ; Patiño, 1997 ; Dartevelle, 2004; Dartevelle et al., 2004; Dufek and Bergantz, 2005; Dartevelle and Valentine, 2007; Ruprecht et al., 2008; Dufek and Bachmann, 2010].

Un milieu granulaire se comporte de façon très différente [Jaeger et al., 1996] en fonction du rapport ɛm/ɛs et de la manière dont la phase dispersée transmet le mouvement à l’intérieur d’elle-même  . Un ensemble de grains posés sur une table peut former un tas statique. Malgré des contraintes de cisaillement présentes dans le tas, le milieu reste immobile (stagnant) et se comporte donc comme un solide . Les théories décrivant ce régime sont gouvernées par la mécanique des solides, spécifiquement dans le champ de « l’élasticité ». Dans ce régime, le système est dominé par les interactions de contact entre les grains. Lorsque les grains sont en mouvement, comme dans un sablier par exemple, les particules interagissent à la fois par collision et contacts frictionnels de longue durée   ; lorsque les contacts entre les grains sont de type frictionnel ce régime est connu comme plastique. Porté à l’extrême, si l’on secoue énergiquement un tas de billes, le milieu devient très agité avec des particules bougeant dans tous les sens, se déplaçant par translation (cinétique) et interagissant par collision. Dans ce régime que l’on appelle collisionnel, le milieu ressemble à un gaz   ; ce régime est connu comme visqueux.

Table des matières

Introduction génèral
Chapitre I. Le modèle hydrodynamique MFIX
Résume
Les équations
Les équations de continuité ou de conservation de la masse
1.2. Les équations de bilan de la quantité de mouvement
1.2.1 Tenseur de déformation ou de contraintes
1.2.1.1. Tenseur de contrainte de la phase liquide
1.2.1.2. Tenseur de contrainte de la phase granulaire
1.2.1.2.1. Tenseur Plastique/Frictionelle
1.2.1.2.2. Tenseur Cinetique-Collisionelle
1.2.2. Interaction entre fluide et particules
1.3. Les équations de bilan d’energie
2. Conditions initiales et aux limites du système
3. Example d’application du modèle MFIX
3.1. Comment lancer un calcul et mettre en place l’expérience analogique MFIX ?
3.1.1. Le lancement d’une simulation vers le centre de calcul
3.1.2. Les paramètres numériques d’affinement
3.1.3. Les conditions aux limites
3.2. Comment traiter les données de sortie ? et quelques examples d’application
3.2.1. L’aire d’interface
3.2.2. Influence de la viscosité (µ) sur la croissance de l’aire d’interface
3.2.3. Influence du coefficient d’expansion thermique (α) sur la croissance de l’aire d’interface
3.2.4 Épaisseur theorique de la couche stagnante
4. Conclusions
Annexe 1
Annexe 2
Annexe 3
Annexe 4
Annexe 5
References
Chapitre II. Numerical simulations of convection crystal-bearing magmas: a case study of the magmatic system at Erebus, Antarctica
Abstract
1. Introduction
2. Physical model description
2.1. Governing equations
2.2. Thermal boundary conditions
2.3. Numerical considerations
3. Validation and verification of the physical model
3.1. Validation and verification of momentum and heat transfer for single-phase simulations
3.1.1. Comparison between numerical and experimental results
3.2. Verification of momentum coupling between melt and solid for bi-phase simulations
3.2.1. Comparison between steady-state, bi-phase flow and mixture theory
4. Erebus lava lake simulations
4.1. Physical parameters of the natural system
4.2. Idealized magmatic system
4.3. Overview of the simulations and numerical considerations
4.4. Results concerning how melt and crystals are treated in the simulations
4.4.1. Pure melt
4.4.2. Crystals as part of the melt
4.4.3. Crystals as a separate phase
4.4.4. Comparison of the three scenarios for representing phases
4.5. The influence of the conduit dimensions
4.6. Sensitivity to boundary conditions
4.6.1. Crystal settling in a closed system (no feeding)
4.6.2. Crystal settling in an open system (permanent feeding)
4.7. Analysis of salient model outputs: Behaviour of the lake surface
5. Discussion
6. Conclusions
Appendix A
Acknowledgements
References
Chapitre III. Numerical simulation of the effusive regime in bubbles-bearing magmas: a case study of the magmatic system at Erebus, Antarctica
Résume
Abstract
1. Introduction
2. Fluid dynamical model
2.1. Governing equations
2.1.1. Volatile exsolution
2.1.3. Sub-grid model for the expantion of the bubles
3. Erebus simulations
3.1. Physical parameters of the natural system
3.2. Idealized magmatic system
4. Results
4.1. Flow velocity
4.2. Reaction rates
4.3. Dissolved water
4.4. Gas volume fraction
4.5. Bubble diameter
4.6. Water release
4.7. Surface velocites at the top of the lake
4.8. Coupling of the mixture and gas velocities
5. Interpretations and implications for the effusive regime of Erebus
6. Main conclusions and future work
Appendix 1. Sub-grid model for gas expansion
Appendix 2. Diffusivity of dissolved water
Appendix 3. Viscosity model
Appendix 4. General notations and parameters of viscosity model
References
Conclusion

Cours gratuitTélécharger le document complet

Télécharger aussi :

Laisser un commentaire

Votre adresse e-mail ne sera pas publiée. Les champs obligatoires sont indiqués avec *