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.

Dans un contexte industriel, l’utilisation de modèles diphasiques d’ordre réduit est nécessaire pour pouvoir effectuer des simulations numériques prédictives d’injection de combustible liquide dans les chambres de combustion automobiles et aéronautiques, afin de concevoir des équipements plus performants et moins polluants. Le processus d’atomisation du combustible, depuis sa sortie de l’injecteur sous un régime de phases séparées, jusqu’au brouillard de gouttelettes dispersées, est l’un des facteurs clés d’une combustion de bonne qualité. Aujourd’hui cependant, la prise en compte de toutes les échelles physiques impliquées dans ce processus nécessite une avancée majeure en termes de modélisation, de méthodes numériques et de calcul haute performance (HPC). Ces trois aspects sont abordés dans cette thèse. Premièrement, des modèles de mélange, dérivés par le principe variationnel de Hamilton et le second principe de la thermodynamique sont étudiés. Ils sont alors enrichis afin de pouvoir décrire des pulsations des interfaces au niveau de la sous-échelle. Des comparaisons avec des données expérimentales dans un contexte de milieux à bulles permettent de vérifier la cohérence physique des modèles et de valider la méthodologie. Deuxièmement, une stratégie de discrétisation est développée, basée sur une séparation d’opérateur, permettant la résolution indépendante de la partie convective des systèmes à l’aide de solveurs de Riemann approchés standards et les termes sources à l’aide d’intégrateurs d’équations différentielles ordinaires. Ces différentes méthodes répondent aux particularités des systèmes diphasiques compressibles, ainsi qu’au choix de l’utilisation de maillages adaptatifs (AMR). Pour ces derniers, une stratégie spécifique est développée : il s’agit du choix de critères de raffinement et de la projection de la solution d’une grille à une autre (plus fine ou plus grossière). Enfin, l’utilisation de l’AMR dans un cadre HPC est rendue possible grâce à la bibliothèque AMR p4est, laquelle a montré une excellente scalabilité jusqu’à plusieurs milliers de coeurs de calcul. Un code applicatif, CanoP, a été développé et permet de simuler des écoulements fluides avec des méthodes de volumes finis sur des maillages AMR. CanoP pourra être utilisé pour des futures simulations d’atomisation liquide. / In an industrial context, reduced-order two-phase models are used in predictive simulations of the liquid fuel injection in combustion chambers and help designing more efficient and less polluting devices. The combustion quality strongly depends on the atomization process, starting from the separated phase flow at the exit of the nozzle down to the cloud of fuel droplets characterized by a disperse-phase flow. Today, simulating all the physical scales involved in this process requires a major breakthrough in terms of modeling, numerical methods and high performance computing (HPC). These three aspects are addressed in this thesis. First, we are interested in mixture models, derived through Hamilton’s variational principle and the second principle of thermodynamics. We enrich these models, so that they can describe sub-scale pulsations mechanisms. Comparisons with experimental data in a context of bubbly flows enables to assess the models and the methodology. Based on a geometrical study of the interface evolution, new tracks are then proposed for further enriching the mixture models using the same methodology. Second, we propose a numerical strategy based on finite volume methods composed of an operator splitting strategy, approximate Riemann solvers for the resolution of the convective part and specific ODE solvers for the source terms. These methods have been adapted so as to handle several difficulties related to two-phase flows, like the large acoustic impedance ratio, the stiffness of the source terms and low-Mach issues. Moreover, a cell-based Adaptive Mesh Refinement (AMR) strategy is considered. This involves to develop refinement criteria, the setting of the solution values on the new grids and to adapt the standard methods for regular structured grids to non-conforming grids. Finally, the scalability of this AMR tool relies on the p4est AMR library, that shows excellent scalability on several thousands cores. A code named CanoP has been developed and enables to solve fluid dynamics equations on AMR grids. We show that CanoP can be used for future simulations of the liquid atomization.