Instabilités thermiques dans un fluide à seuil : de l'échelle du laboratoire à celle de la planète

Thermal instabilities in a yield-stress fluid : from the laboratory to the planetary scale

• Auteur / Autrice : Anna Massmeyer
• Direction : Anne Davaille
• Type : Thèse de doctorat
• Discipline(s) : Physique
• Date : Soutenance le 27/05/2013
• Etablissement(s) : Paris 11
• Ecole(s) doctorale(s) : Ecole doctorale Modélisation et Instrumentation en Physique, Energie, Géosciences et Environnement (Orsay, Essonne ; 2010-2015)
• Partenaire(s) de recherche : Laboratoire : Fluides, automatique et systèmes thermiques (Orsay, Essonne ; 1997-....)
• Jury : Président / Présidente : Andréa Tommasi
• Jury : Examinateurs / Examinatrices : Anne Davaille, Andréa Tommasi, Patrick Cordier, Chérif Nouar, Harold Auradou, Thibaut Divoux, Paul Tackley, Erika Di Giuseppe
• Jury : Rapporteurs / Rapporteuses : Patrick Cordier, Chérif Nouar

Mots clés

Mots clés contrôlés

Panaches mantelliques

Fluides complexes

Chaleur -- Convection

Mots clés libres

Convection

Carbopol

Diapirisme

Panache

Simulations numériques

Expériences de laboratoire

Modèle de Herschel-Bulkley

Contrainte seuil

Convection

Carbopol

Diapirism

Plume

Numerical simulations

Laboratory experiments

Herschel-Bulkley model

Yield-stress

Résumé

Des panaches sont connus pour migrer à travers le manteau ductile et quasi-Newtonien ; alors que les dikes se fracturent et se propagent dans la lithosphère solide. Cependant, la lithosphère est en fait visco-élastique. Afin de déterminer ce qui se passe dans ce cas complexe, nous avons réalisé une étude expérimentale et numérique sur le développement de panaches thermiques dans des solutions aqueuses de Carbopol, un gel de polymères formant un réseau continu d'éponges microscopiques. Ce fluide est rhéofluidifiant et présente un seuil de contrainte σ₀, de sorte que l'écoulement ne se produit que si la contrainte locale dépasse cette valeur critique σ₀. En dessous de cette valeur, le fluide agit comme un solide élastique. Les propriétés rhéologiques des solutions peuvent être systématiquement ajustées en variant la concentration de Carbopol. Le dispositif consiste en une source locale de chaleur de puissance constante placée au centre d'une cuve cubique. Selon la valeur du rapport entre la contrainte d'origine thermique et la contrainte seuil, Y₀, on peut observer trois régimes différents. A faible Y₀<Yc₁, aucun mouvement n'est détecté ; tandis que pour Yc₁<Y₀<Yc₂ une cellule se développe, puis évolue vers un panache pour Y₀>Yc₂. Nous montrons que les paramètres critiques (Yc₁,Yc₂) dépendent fortement de la géométrie du chauffage. Des mesures simultanées de température et de champs de vitesse montrent que la morphologie du panache ressemble à un doigt, contrairement à la forme de champignon rencontrée dans les fluides newtoniens. Utilisant des simulations numériques avec une description purement visqueuse, où la rhéologie du fluide est décrite par un modèle de Herschel-Bulkley régularisé, sont suffisantes pour rendre compte de la dynamique du panache. Une étude détaillée des paramètres indiquent que la dynamique du panache est gouvernée par la compétition entre la contrainte seuil, la contrainte induite par la flottabilité et les contraintes visqueuses. Nous avons identifié deux paramètres adimensionnés : le paramètre seuil Ψ comparant la contrainte induite par la flottabilité et la contrainte seuil, et le nombre de Bingham Bi comparant la contrainte seuil et les contraintes visqueuses. Un panache ne peut s'élever que si les deux paramètres sont supercritiques, i.e. la contrainte induite par la flottabilité et les contraintes visqueuses sont plus importantes que la contrainte seuil. Par conséquent, le panache peut s'arrêter avant d'atteindre la surface. Des lois d'échelles dans le conduit du panache ont été déterminées pour la vitesse, la température et la taille de la région cisaillée en régime permanent. Elles décrivent raisonnablement le comportement du conduit bien que seul l'effet rhéofluidifiant soit pris en compte. L'application de ces paramètres adimensionnés à la Terre contraignent significativement la limite de plasticité du manteau et de la lithosphère. La contrainte seuil maximale qui permet à une instabilité thermique de pénétrer dans la lithosphère ou le manteau supérieur est entre 100 kPa et 100 MPa, et elle dépend fortement de la taille et de l'anomalie de densité de l'intrusion.

Plumes are known to migrate through the ductile quasi-Newtonian mantle, while dikes fracture and propagate through the solid lithosphere. However, depending on the timescale, the lithosphere presents solid as well as viscous properties. To determine what happens in the complex case, where instabilities propagate through a visco-elastic matrix, we performed a combined study of laboratory experiments and numerical simulations. Here we investigate the development of thermal plumes in aqueous solutions of Carbopol, a polymer gel, forming a continuous network of micrometric sponges. This fluid is shear thinning and has a yield-stress σ₀, whereby flow occurs only if the local stress exceeds this critical value σ₀. Below this value, the fluid acts as an elastic solid. The rheological properties of the solutions can be systematically varied by varying the Carbopol concentration. The setup consists of a localized heat-source operated at constant power, placed at the centre of a square tank. Depending on the ratio of the thermally induced stresses and the yield stress, Y₀, three different regimes may be obtained. For low Y₀<Yc₁ no motion occurs, whereas for Yc₁<Y₀<Yc₂ a cell develops, that evolves into a plume for Y₀>Yc₂. We show that the critical parameters (Yc₁,Yc₂) strongly depend on the geometry of the heating. Combined temperature and velocity field measurements show that the morphology of the plume resembles a finger, contrary to the mushroom-like shape encountered in Newtonian fluids. Numerical simulations using a purely viscous description, where the rheology of the fluid is described by a regularized Herschel-Bulkley model, are sufficient to capture the plume dynamics. A detailed parametric study shows that the plume dynamics are governed by the interplay between yield stress, buoyancy induced stress and viscous stresses. We identify two non-dimensional parameters: the yield parameter Ψ comparing the buoyancy induced stress to the yield stress, and the Bingham number Bi comparing the yield stress to the viscous stresses. We show that a plume can rise only if both parameters are supercritical, i.e. if buoyancy induced stress and viscous stresses each overcome the yield stress. Therefore the plume may come to a halt before it reaches the surface. We propose scaling laws for the plume stem velocity, temperature and the size of the shear zone in the steady state. We show that the scaling laws describe the behaviour in the plume stem reasonably well, if the yield stress is neglected and only the shear thinning behaviour is taken into account. Applying the non-dimensional parameters to Earth places severe constraints on the strength of mantle and lithosphere. The maximum strength that allows for thermal instabilities to penetrate the lithosphere or upper mantle is in between 100 kPa and 100 MPa, and strongly depends on the size and buoyancy of the anomaly.