Approximation numérique et modélisation de l'ablation liquide

Numerical approximation and modelling of liquid ablation

• Auteur / Autrice : Simon Peluchon
• Direction : Luc Mieussens
• Type : Thèse de doctorat
• Discipline(s) : Mathématiques appliquées et calcul scientifique
• Date : Soutenance le 28/11/2017
• Etablissement(s) : Bordeaux
• Ecole(s) doctorale(s) : École doctorale Mathématiques et informatique (Talence, Gironde ; 1991-....)
• Partenaire(s) de recherche : Laboratoire : Centre d'études scientifiques et techniques d'Aquitaine (Le Barp, France)
• Jury : Président / Présidente : Mario Ricchiuto
• Jury : Examinateurs / Examinatrices : Héloïse Beaugendre, Gérard Gallice, Thierry Magin, Philippe Villedieu
• Jury : Rapporteurs / Rapporteuses : Rémi Abgrall, Christophe Chalons

Mots clés

Mots clés contrôlés

Mach, Nombre de

Écoulement diphasique

Analyse numérique

Mécanique des fluides

Mots clés libres

Écoulements diphasiques

Schémas de type Godunov

Méthode implicite

Écoulements bas Mach

Navier-Stokes compressible

Two-phase flow

Godunov-type schemes

Implicit method

Low Mach

Compressible Navier-Stokes

Résumé

Lors de sa rentrée dans l’atmosphère d’une planète, un engin spatial subit un échauffement important dû aux frottements des gaz atmosphériques sur la paroi. Cette élévation de température conduit à une dégradation physico-chimique du bouclier thermique de l’objet constitué de matériaux composites. Un composite est constitué de divers matériaux qui s’ablatent différemment. Dans cette thèse, nous nous intéressons essentiellement à la fusion d’un matériau durant sa phase de rentrée atmosphérique. Nous sommes donc en présence de trois phases : solide, liquide et gaz. Pour simuler ce phénomène, des méthodes numériques robustes ont été mises au point pour calculer l’écoulement diphasique compressible autour de l’objet. Le couplage entre le solide et l’écoulement fluide a aussi été étudié. Les méthodes numériques développées durant cette thèse sont basées sur une approche volumes finis. Une stratégie de décomposition d’opérateurs est utilisée pour résoudre le modèle diphasique à cinq équations avec les termes de dissipation modélisant l’écoulement fluide. L’idée principale de cette décomposition d’opérateurs est de séparer les phénomènes acoustiques et dissipatifs des phénomènes de transport. Un traitement implicite de l’étape acoustique est réalisé tandis que l’étape de transport est résolue explicitement. Le schéma semi-implicite global est alors très robuste, conservatif et préserve les discontinuités de contact. Les conditions d’interface entre les domaines fluide et solide sont déduites des bilans de masse et d’énergie à la paroi. Le front de fusion est suivi explicitement grâce à une formulation ALE des équations. La robustesse de l’approche et l’apport de la formulation semi-implicite sont finalement démontrés grâce à des expériences numériques mono et bidimensionnelles sur maillages curvilignes mobiles.

During atmospheric re-entry phase, a spacecraft undergoes a sudden increase of the temperature due to the friction of atmospheric gases. This rise drives to a physical-chemical degradation of the thermal protective system of the object made of composite material. A composite is made of several materials with ablates differently. In this thesis, we mainly focus on the melting of an object during its re-entry phase. Therefore there are three phases: solid, liquid and gas phases. In order to simulate this phenomenon, robust numerical methods have been developed to compute a compressible multiphase flow. The coupling strategy between the solid and the fluid have also been studied. Solvers developed in the present work are based on Finite Volume Method. A splitting strategy is used to compute compressible two-phase flows using the five-equation model with viscous and heat conduction effects. The main idea of the splitting is to separate the acoustic and dissipative phenomena from the transport one. An implicit treatment of the acoustic step is performed while the transport step is solved explicitly. The overall scheme resulting from this splitting operator strategy is very robust, conservative, and preserves contact discontinuities. The boundary interface condition between the solid and the multiphase flow is enforced by mass and energy balances at the wall. The melting front is tracked explicitly using an ALE formulation of the equations. The robustness of the approach and the interest of the semi-implicit formulation are demonstrated through numerical simulations in one and two dimensions on moving curvilinear grids.