Panache volcanique

Un panache volcanique, panache éruptif, panache de cendres, nuage de cendres, colonne éruptive ou colonne volcanique est un nuage de gaz volcaniques chauds et de téphras (majoritairement des cendres), qui s'élève en altitude au-dessus d'un volcan au cours d'une éruption, notamment une éruption explosive.

Pour les articles homonymes, voir Panache.

Pour un article plus général, voir Panache (hydrodynamique).

Panache volcanique s'élevant au-dessus du Pinatubo aux Philippines au cours de son éruption plinienne de 1991.
Panache volcanique de petite taille s'élevant d'une bouche éruptive du Puʻu ʻŌʻō à Hawaii au cours d'une éruption hawaïenne en 2005.

Dans le cas des éruptions explosives, il arrive souvent que le panache volcanique, en montant et en se refroidissant, devienne plus dense que l'air ambiant et retombe sur les pentes du volcan (effondrement de panache) pour former des coulées pyroclastiques ou des déferlantes. Moins fréquemment, le panache reste moins dense que l'air ambiant et s'élève par convection pour former un pyrocumulonimbus.

Les panaches des éruptions les plus violentes, d'indice d'explosivité 4 à 8 (voire seulement 3 dans certains cas) atteignent la stratosphère et peuvent affecter le climat à moyen et long termes, jusqu'à provoquer un « hiver volcanique ». Les pyrocumulonimbus, troposphériques, perturbent le trafic aérien.

Caractéristiques

Même si tous les volcans émettent un panache volcanique lorsqu'ils entrent en éruption, les plus grands (plusieurs dizaines de kilomètres d'altitude) surviennent au cours d'éruptions pliniennes sur des volcans gris, les volcans rouges ne produisant parfois que des panaches volcaniques composés uniquement de gaz volcaniques et s'élevant peu en altitude (quelques dizaines de mètres parfois).

Modélisation analytique

Historique

Les travaux entrepris dans la description du comportement des panaches verticaux tels que les panaches volcaniques dans un milieu homogène ont été commencés par Rouse[1] qui a montré que le profil des propriétés d'écoulement est gaussien. Morton[2] a introduit l'hypothèse d'entraînement qui marquera une avancée majeure dans la modélisation théorique des panaches verticaux. Il a proposé une description simple de la nature du panache à travers un modèle théorique. Il a utilisé dans son approche l’hypothèse d’entraînement qui a été introduite précédemment par Taylor[3]. Cette hypothèse modélise le mécanisme d'aspiration du fluide ambiant dans le panache par une vitesse d’entrée perpendiculaire à l’axe du panache, et directement proportionnelle à la vitesse moyenne locale ascendante du panache.

Description du modèle d'entraînement de Morton
Modélisation numérique d'un panache vertical

Pour affiner les résultats sur le comportement du panache, Morton a, après analyse et développement, repris l’hypothèse de Taylor dans son étude. Il a supposé un profil uniforme (top-hat), i.e. une valeur constante à travers le panache et zéro ailleurs. Bien que les propriétés de l’écoulement telles que la vitesse et la densité suivent une distribution gaussienne le long de la zone de l'écoulement établi du panache, il est convenable de supposer une telle hypothèse pour simplifier l'écriture des équations. En fait, le passage du profil gaussien au profil top-hat peut être fait en effectuant une intégration le long de la section transversale du panache, i.e. obtenir une valeur moyenne sur la propriété voulue.

Description du modèle physique

Le modèle physique est représenté schématiquement sur la figure ci-contre. Il consiste en une éjection verticale d’un fluide léger turbulent dans un milieu plus dense homogène, depuis un orifice de rayon , à une vitesse initiale et une densité .

Il est considéré dans ce développement que le régime de l’écoulement est stationnaire et permanent et que le fluide est incompressible. La dispersion du panache dans un fluide homogène, est décrit mathématiquement au travers de diverses équations de conservation afin de prédire les trajectoires, les concentrations, les vitesses, les densités et le rayon en fonction de la distance verticale le long du panache. Il est à noter que cette étude concerne aussi bien le cas restrictif de l'approximation de Boussinesq que le cas général non-Boussinesq[incompréhensible].

Le nombre de Reynolds est suffisamment grand pour s'assurer que le régime d'écoulement est bien turbulent, ce qui nous permet d'affirmer que, ni le nombre de Reynolds, ni les propriétés moléculaires du fluide telles que la viscosité, n'entrent directement dans la détermination des propriétés du panache. Morton[4] précise la limite d'application de l'hypothèse d'entraînement : dans sa note, il indique que, pour les écoulements où l'équilibre n'est pas atteint, l'entraînement doit être lié non pas à la vitesse locale moyenne, mais plutôt au tenseur de Reynolds[incompréhensible].

Mise en équations

En considérant un élément de volume cylindrique, les équations de conservation de la masse, de la quantité de mouvement et de l'énergie sont écrites sous la forme suivante :

 ;
 ;
.

Selon l'hypothèse de Morton, la vitesse d'entraînement est proportionnelle à la vitesse locale du panache, ainsi :

 :

est la distance verticale ;
est la densité locale du panache ;
est la densité du milieu, considérée constante dans le cas homogène non stratifié ;
désigne la vitesse du panache, elle possède une seule composante suivant  ;
représente le rayon du panache qui délimite la frontière avec le milieu ;
désigne le coefficient d'entraînement ou appelé aussi le taux d'aspiration du fluide ambiant ;
est l'accélération de la pesanteur.

Michaux et al.[5] ont donné une formulation générale du problème, dans le cadre ou non de l'approximation de Boussinesq, ils ont introduit le rayon modifié , le déficit de densité adimensionnel et le coefficient d'entraînement corrigé définis respectivement par :

,

tels que dans le cas de l'approximation de Boussinesq et dans le cas non-Boussinesq.

désigne l'écart entre la densité du milieu et celle du panache.

Le système d'équations différentielles précédant devient ainsi :

Introduction de la fonction panache Γ

On introduit la fonction panache définie par :

Cette fonction n'est en fait qu'un nombre de Richardson normalisé, il traduit la compétition entre les effets de flottabilité et de quantité de mouvement.

Pour des faibles valeurs de , l'effet de la quantité de mouvement est prépondérant par rapport à celui de flottabilité, tandis que pour des fortes valeurs de , l'effet de flottabilité domine.

Forme d'un panache paresseux () avec présence d'un col (section minimale)
Forme d'un panache forcé ()
Forme d'un panache pur ()
Variations de la fonction panache en fonction de l'altitude pour un panache paresseux, forcé ou pur

Les équations différentielles qui régissent l'écoulement peuvent s'exprimer directement en fonction de  :

Ces dernières équations traduisent que le rayon modifié possède une valeur minimale (un col) et ceci correspond à , et que la vitesse atteint une valeur maximale à l'endroit où , autrement dit, il y a accélération de l'écoulement jusqu'à atteindre une valeur maximale, ceci a lieu juste après la section minimale.

L'équation définissant la fonction panache couplée avec ces dernières équations permet d'obtenir une équation différentielle en  :

L'avantage de l'introduction de la fonction panache est de pouvoir exprimer toutes les autres grandeurs uniquement en fonction de. Ainsi :

L'indice désigne la valeur de la variable au niveau de la source.

Il est ainsi possible de déterminer les valeurs de la vitesse maximale et de la section minimale :

Suivant la valeur de la fonction panache au niveau de la source, est la solution de :

avec :

Différentes configurations des panaches verticaux

L'introduction de la fonction panache permet de savoir la nature du panache suivant la valeur de .

  • Pour , la fonction est strictement décroissante, elle tend vers l'unité par valeurs supérieures lorsque tend vers l'infini. Ceci signifie que l'effet de flottabilité est dominant à proximité de la source. Le panache est dit paresseux (lazy).
  • Pour , la fonction est strictement croissante, elle tend de même vers l'unité par valeurs inférieures lorsque tend vers l'infini. Le comportement du jet domine au voisinage de la source et l'effet de flottabilité est négligeable devant celui de quantité de mouvement, mais plus loin de la source, les deux flux s'équilibrent, et la fonction panache tend vers l'unité. Le panache est dit dans ce cas forcé.
  • Pour , la fonction est constamment constante et valant l'unité, indépendamment de l'altitude , les deux flux sont toujours équilibrés et le panache est dit pur. C'est uniquement dans ce cas que des solutions analytiques exactes peuvent être retrouvées.

Solution semi-analytique

La solution du problème est donnée par :

avec :

tel que désigne la fonction bêta incomplète de paramètres et .

Il existe une autre formulation pour écrire ces deux équations précédentes d'une façon unique en remplaçant par son inverse dans le cas forcé, ceci va permettre d'effectuer une seule résolution à la fois :

 :

et pour le cas lazy et pour le cas forcé, vaut pour les deux cas .

Ainsi :

avec

Des solutions analytiques asymptotiques raccordées, déterminées par la méthode des perturbations ont été données par Candelier et al.[6].

Résolution dans le cas d'un panache pur

La fonction de panache vaut constamment l'unité, la résolution analytique donne directement les variations de , de et de en fonction de  :

Modélisation numérique

Les progrès de l'informatique et de la prévision météorologique permettent de maintenant mieux modéliser le déplacement de ces panaches (quand ils sont importants et tant qu'ils ne sont pas trop dilués)[7],[8], y compris pour une couche donnée d'altitude[9].

Notes et références

  1. (en)H. ROUSE, C. S. YIH and H. W. HUMPHREYS. Gravitational Convection from a Boundary Source. Tellus, 4: 201–210. doi: 10.1111/j.2153-3490.1952.tb01005.x (1952).
  2. (en)B. R. Morton, G. I. Taylor and J. S. Turner. Turbulent gravitational convection from maintained and instantaneous sources. Proc. R. Soc. Lond. A 234, 1-23 (1956).
  3. (en)G. I. Taylor. Dynamics of a mass of hot gas rising in air. US. Atomic Energy Commission MDDC 919. LADC 276 (1945).
  4. (en)B. R. Morton B. R. On similarity and the outer boundary conditions for flow in plumes. Tellus 20, 367-369 (1968).
  5. (en)G. Michaux, G. an O. Vauquelin. Solutions for turbulent buoyant plumes rising from circular sources. Phys. Fluids 20, 066601 (2008).
  6. (en)F. Candelier and O. Vauquelin. Mached asymptotic solutions for turbulent plumes. J. Fluid Mech. 699, 489-499 (2012).
  7. Exemple de modélisation prévisionnelle du déplacement du panache de l'Eyjafjöll : animation montrant le déplacement du panache pour la semaine du 16 au 21 avril 2010 (INERIS, consulté le 23 avril 2010).
  8. Modélisation du même panache de cendres - impact au sol (INERIS, consulté le 23 avril 2010).
  9. Exemple de modélisation du déplacement d'un nuage de cendres de l'Eyjafjöll à cinq kilomètres d'altitude du 19 au 21 avril 2010 (INERIS, consulté le 23 avril 2010).

Voir aussi

Articles connexes

  • Portail de la physique
  • Portail de l'analyse
  • Portail du volcanisme
  • Portail de l’environnement
Cet article est issu de Wikipedia. Le texte est sous licence Creative Commons - Attribution - Partage dans les Mêmes. Des conditions supplémentaires peuvent s'appliquer aux fichiers multimédias.