Radiative transfers in earth atmospheric re-entry's hypersonic flows

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.Thermal protection design of spacecraft requires the knowledge of heat fluxes at the vehicle surface. For the kind of atmospheric entry considered in this work (earth reentry at 10-12 km/s or more), radiation significantly contributes to the surface heating. We present in this study, the development and the use of models and numerical tools to predict radiative transfers in the shock layers encountered in earth re-entry in nonequilibrium conditions and multi-dimensionnal geometry. First of all, we have formulated, in a line by line approach, the expressions of radiative properties of N2-O2 plasmas in nonequilibrium conditions. This formulation is suitable for a multi-temperature and/or a electronic state to state description of the thermal nonequilibrium. It has been used to simulate the test case FIRE II in order to determine radiative intensity at the stagnation point for four trajectory points. Results show a good agreement with flight data for different spectral measurement ranges form IR to UV, except for the trajectory point in strong nonequilibrium conditions. Calculations show moreover that VUV spectral range contributes significantly to the intensity at the wall and that taking into account chemical nonequilibrium is crucial. An approximate model of radiative properties has been developed on the basis of a statistical narrow band (SNB) model for optically thick (in our application) electronic molecular systems, of a box model for optically thin electronic molecular systems and continua, and of a line by line approach for atomic lines. The use of the hybrid model has required a suitable RTE formulation to take into account nonequilibrium and spectral correlations. Systematic validations of the SNB model for each contribution have been carried out. The hybrid model has then been validated on the test case FIRE II against line by line calculations. The hybrid model has been implemented in the ASTRE radiative solver, which is based on a Monte Carlo approach. The implementation has been validated in comparison with a ray tracing method on a tangent slab configuration. 3D radiative calculations have been then carried out on the test case FIRE II. Results show discrepancy of 10-15 % on the fluxes at the stagnation point with results obtained in the tangent slab approximation which is usually used in the litterature. The developed tools allow to obtain relatively well converged (5 %) results in tri-dimensional geometry (100000 cells) with a reasonable CPU time enabling to consider coupled calculations in futur works

Transferts radiatifs dans les écoulements hypersoniques de rentrée atmosphérique terrestre

Thermal protection design of spacecraft requires the knowledge of heat fluxes at the vehicle surface. For the kind of atmospheric entry considered in this work (earth reentry at 10-12 km/s or more), radiation significantly contributes to the surface heating. We present in this study, the development and the use of models and numerical tools to predict radiative transfers in the shock layers encountered in earth re-entry in nonequilibrium conditions and multi-dimensionnal geometry. First of all, we have formulated, in a line by line approach, the expressions of radiative properties of N2-O2 plasmas in nonequilibrium conditions. This formulation is suitable for a multi-temperature and/or a electronic state to state description of the thermal nonequilibrium. It has been used to simulate the test case FIRE II in order to determine radiative intensity at the stagnation point for four trajectory points. Results show a good agreement with flight data for different spectral measurement ranges form IR to UV, except for the trajectory point in strong nonequilibrium conditions. Calculations show moreover that VUV spectral range contributes significantly to the intensity at the wall and that taking into account chemical nonequilibrium is crucial. An approximate model of radiative properties has been developed on the basis of a statistical narrow band (SNB) model for optically thick (in our application) electronic molecular systems, of a box model for optically thin electronic molecular systems and continua, and of a line by line approach for atomic lines. The use of the hybrid model has required a suitable RTE formulation to take into account nonequilibrium and spectral correlations. Systematic validations of the SNB model for each contribution have been carried out. The hybrid model has then been validated on the test case FIRE II against line by line calculations. The hybrid model has been implemented in the ASTRE radiative solver, which is based on a Monte Carlo approach. The implementation has been validated in comparison with a ray tracing method on a tangent slab configuration. 3D radiative calculations have been then carried out on the test case FIRE II. Results show discrepancy of 10-15 % on the fluxes at the stagnation point with results obtained in the tangent slab approximation which is usually used in the litterature. The developed tools allow to obtain relatively well converged (5 %) results in tri-dimensional geometry (100000 cells) with a reasonable CPU time enabling to consider coupled calculations in futur works.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

Transferts radiatifs dans les écoulements hypersoniques de rentrée atmosphérique terrestre

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.Thermal protection design of spacecraft requires the knowledge of heat fluxes at the vehicle surface. For the kind of atmospheric entry considered in this work (earth reentry at 10-12 km/s or more), radiation significantly contributes to the surface heating. We present in this study, the development and the use of models and numerical tools to predict radiative transfers in the shock layers encountered in earth re-entry in nonequilibrium conditions and multi-dimensionnal geometry. First of all, we have formulated, in a line by line approach, the expressions of radiative properties of N2-O2 plasmas in nonequilibrium conditions. This formulation is suitable for a multi-temperature and/or a electronic state to state description of the thermal nonequilibrium. It has been used to simulate the test case FIRE II in order to determine radiative intensity at the stagnation point for four trajectory points. Results show a good agreement with flight data for different spectral measurement ranges form IR to UV, except for the trajectory point in strong nonequilibrium conditions. Calculations show moreover that VUV spectral range contributes significantly to the intensity at the wall and that taking into account chemical nonequilibrium is crucial. An approximate model of radiative properties has been developed on the basis of a statistical narrow band (SNB) model for optically thick (in our application) electronic molecular systems, of a box model for optically thin electronic molecular systems and continua, and of a line by line approach for atomic lines. The use of the hybrid model has required a suitable RTE formulation to take into account nonequilibrium and spectral correlations. Systematic validations of the SNB model for each contribution have been carried out. The hybrid model has then been validated on the test case FIRE II against line by line calculations. The hybrid model has been implemented in the ASTRE radiative solver, which is based on a Monte Carlo approach. The implementation has been validated in comparison with a ray tracing method on a tangent slab configuration. 3D radiative calculations have been then carried out on the test case FIRE II. Results show discrepancy of 10-15 % on the fluxes at the stagnation point with results obtained in the tangent slab approximation which is usually used in the litterature. The developed tools allow to obtain relatively well converged (5 %) results in tri-dimensional geometry (100000 cells) with a reasonable CPU time enabling to consider coupled calculations in futur works.CHATENAY MALABRY-Ecole centrale (920192301) / SudocSudocFranceF

Narrow-band model for nonequilibrium air plasma radiation

International audienc

Numerical simulation of radiation in high altitude solid propellant rocket plumes

International audienceModels and computational schemes are presented to simulate the flow field and radiation in solid propellant rocket plumes at high altitude. Different models are developed to calculate plume radiation induced by both the gas and the dispersed phase of alumina particles. Simulations were made to reproduce a 110 km altitude rocket plume and comparisons with experimental data on the plume radiance are presented. Influence of different physical processes on the simulation of the plume flow field and radiation is discussed, such as alumina phase change, scattering and the coupling of radiation with the two-phase flow

CO2 IR radiation modelling with a multi-temperature approach in flows under vibrational nonequilibrium

International audienceThe aim of this paper is to develop a line by line model for CO2 vibrational nonequilibrium radiation and to investigate nonequilibrium effects in the case of some applications involving high temperature expanding flows. A vibrational specific collisional relaxation model is developed and is incorporated in a multi- temperature thermodynamic description of the gas mixture, in order to compute vibrational level populations along the flow. The HITEMP-2010 spectroscopic database is employed with a model for level energy splitting to provide line by line absorption and emission total and per-vibrational mode specific spectra. The specific spectra allow us to derive the radiative source terms to be used in the multi-temperature model if a coupled approach is required. The model is applied first to a simple conical nozzle flow and then to the plume of a high altitude solid propellant rocket engine. It is shown that, for the considered applications, the partial freeze of vibrational excitation in the expanding flow increases significantly the radiative intensity escaping from the mixture

Radiation modelling for air plasma in nonequilibrium conditions: application to earth atmospheric re−entry

International audienc

Modelling radiative properties of gas mixtures in nonequilibrium high-altitude rocket plumes

International audienceThis paper presents the development of a Statistical Narrow Band model (SNB) in a nonequilibrium vibrational state for HCl and CO molecules. The population densities of the energy levels are obtained by a multi-temperature approach to compute nonequilibrium Line By Line (LBL) spectra. The SNB parameters are obtained by fitting the curves of growth from the LBL approach by a least squares error minimization using a Newton method for pure Lorentz and Doppler broadening regimes. The model is tested in Voigt broadening regime using a mixing rule and agrees well with the LBL approach. Finally, spectral correlation problems between η_σ /κ_σ and κ_σ , where η_σ and κ_σ are the emission and absorption coefficients respectively, have been exhibited for CO_2