Le dimensionnement des protections thermiques des véhicules spatiaux nécessite la connaissance des flux de chaleur à la surface de l'engin. Pour le type de rentrée atmosphérique étudié dans ces travaux (rentrée terrestre à 10-12 km/s ou plus), le rayonnement contribue de manière significative à l'échauffement pariétal. Nous présentons dans cette étude, le développement et la mise en œuvre de modèles et d'outils numériques permettant de prédire les transferts radiatifs dans les couches de choc de rentrée terrestre hors d'équilibre en géométrie multi-dimensionnelle. Nous avons d'abord formulé dans une approche raie par raie l'expression des propriétés radiatives relatives à des plasmas de N2-O2 hors d'équilibre. Cette formulation est adaptée à une modélisation multi-température et/ou une modélisation par niveau électronique du déséquilibre thermique. Les expressions obtenues ont été mises en œuvre sur le cas d'épreuve FIRE II afin de déterminer les luminances au point d'arrêt pour quatre points de trajectoire. Les résultats montrent un bon accord global avec les données en vol dans les différents intervalles de mesure allant de l'IR à l'UV, hormis pour le point de rentrée le plus hors d'équilibre. Les calculs montrent également que la partie VUV du spectre contribue de manière très importante à la luminance à la paroi et que la prise en compte du déséquilibre chimique est primordiale. Un modèle approché de propriétés radiatives a ensuite été développé sur la base d'un modèle statistique à bandes étroites (MSBE) pour les systèmes électroniques moléculaires optiquement non minces (dans notre application), d'un modèle à coefficients d'émission et d'absorption moyennés par bande étroite pour les systèmes électroniques optiquement minces et les continua, et d'une approche raie par raie pour les transitions lié-lié atomiques. La mise en œuvre de cette modélisation hybride a nécessité une formulation adaptée de l'ETR pour tenir compte du déséquilibre et des corrélations spectrales. Des validations systématiques du MSBE, contribution par contribution, ont été réalisées. Le modèle hybride a ensuite été validé sur le cas d'épreuve FIRE II par rapport à des calculs raie par raie. Le modèle hybride a été implémenté dans le solveur radiatif ASTRE, code basé sur une approche de Monte Carlo. L'implémentation a été validée par rapport à une méthode de lancer de rayons sur une con_figuration plans tangents. Des calculs 3D de rayonnement ont ensuite été réalisés sur le cas d'épreuve FIRE II. Les résultats montrent des écarts de l'ordre de 10-15 % sur les flux au point d'arrêt par rapport aux résultats obtenus dans l'approximation des plans tangents habituellement utilisée dans la littérature. Les outils développés permettent d'obtenir des résultats relativement bien convergés (5 %) en géométrie tri-dimensionnelle (100000 cellules) en un temps CPU raisonnable permettant d'envisager des calculs couplés à terme.
Identifer | oai:union.ndltd.org:CCSD/oai:tel.archives-ouvertes.fr:tel-00453309 |
Date | 21 September 2009 |
Creators | Lamet, Jean-Michel |
Publisher | Ecole Centrale Paris |
Source Sets | CCSD theses-EN-ligne, France |
Language | French |
Detected Language | French |
Type | PhD thesis |
Page generated in 0.0047 seconds