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

Jean-Michel Lamet

Description

Sous la direction de Marie-Yvonne Perrin
Thèse soutenue le 21 septembre 2009: Ecole centrale Paris
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.
-Rayonnement hors d'équilibre
-Modèle de bandes
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.
-Atmospheric re-entry
-Nonequilibirum radiation
-Air plasma
Source: http://www.theses.fr/2009ECAP0026/document

Informations

• Publié par: Thesee
• Nombre de lectures: 53
• Langue: Français
• Poids de l'ouvrage: 16 Mo

Extrait

Nationale
Télécopie
P
Cen
trale
trale
el.
P
288,
aris
Ver
THÈSE
t
présen
92295
tée
(standard)
par
ul
Jean-Mic
Tél.
h
2009
el
Rapp
Lamet
Man
p
Ministère
our
des
l'obten
MALABR
tion
41
du
634
GRADE
d'Energétique
de
Com
DOCTEUR
Cen
F
41
ormation
47
do
xx
cto
P
rale
Ecole
:
Arts
E
Grand
n
tutelle
er
l
gétique
V
Lab
g
oratoire
TENA
d
Cedex
accueil
33
:
10
Lab
elex
oratoire
F
d'Energétique
Lab
Moléculaire
l
et
et
Macroscopique,
(E.M2.C.)
Com
et
bustion
P
(EM2C)
33
du
10
CNRS
33
et
80
de
Vérant
l'ECP
J.-L.
T
visch
ransferts
.
radiat
orteur
ifs
Cen
dans
des
les
et
écoulemen
ufactures
ts
Etablissemen
h
sous
yp
du
ersoniques
de
de
Education
ren
Grande
trée
oie
atmosphérique
Vi
terrestre
nes
Souten
CHA
ue
Y
le
Y
21
T
S
:
e
(1)
p
13
t
00
em
T
bre
:
2009
991
Comp
EC
osition
ARIS
du
oratoire
jury
Mo
:
éc
Mme
aire
Perrin
Macroscopique,
M.-Y.
bustion
Stéhlé
UPR
C.
CNRS
Rapp
École
orteur
trale
MM.
aris
Soufiani
:
A.
(1)
Surzhik
13
o
31
v
:
S.
(1)
T.
02
Tessé
35
L.
-
Tran
Ph.
École
tel-00453309, version 1 - 4 Feb 2010tel-00453309, version 1 - 4 Feb 2010p
e
sein
ts
t
C
ort
e
à
tra
déplacemen
v
J'adresse
ail
our
a
conance
été
hoisi
réali
de
sé
de
dans
p
le
à
cadre
qui
d
cours
'
moi.
un
o
e
a
collab
Olivier
oration
Véran
en
Je
tre
m'a
l'ONERA
ec
et
v
le
ainsi
lab
p
oratoire
ermanen
EM2C
ail
à
athlètes
tra
t
v
qui
ers
d'a
laquelle
com
de
p
nom
je
breuses
gratitude
p
par
ersonnes
t
on
n
t
a
con
cette
tribué
son
au
onctuellemen
b
mem
on
V
déroulemen
r
t
a
de
Stéhlé
cette
thèse,
thèse.
de
Je
p
ti
T
ens,
our
dans
oratoire
les
ctoran
trop
on
courts
biance
paragraphes
et
qui
t
suiv
oratoire
en
parc
t,
Un
à
et
leur
été
faire
;
part
le
de
thèse,
toute
t
ma
aux.
gratitude.
son
T
cette
out
exprime
d'ab
te
ord,
our
mes
m'a
premières
suite,
p
puy
ensées
candidature
v
l'ONERA.
on
oublie
t
qui
à
t
Marie-Y
our
v
et
onne
p
P
t,
errin,
mais
ma
présen
directrice
le
de
du
thèse,
:
Anouar
h
Souani
l'honne
et
présider
Philipp
et
e
a
Rivière
an
qui
c
m'on
orter
t
Surzhik
accordé
a
leur
aire
conance
depuis
et
être
qui
Philipp
m'on
J'ai
t
amicale
encadré
ble
au
du
quotidien
do
au
p
sein
et
du
s,
lab
con
oratoire
une
EM2C.
tra
De
haleu-
part
vivia
leurs
remerciem
comp
aux
étences,
du
leurs
v
conseils
tours
a
Sceaux
vi
fort
sés
mo
et
ma
leurs
amis
grandes
t
disp
ts
onibilités
mes
tout
les
au
oir
long
t
de
soutenance
ma
jour
thèse,
imp
ils
our
m'on
v
t
Outre
p
our
e
app
rmis
dans
d'eectuer
thèse,
mes
lui
tra
t
v
u
aux
ma
dans
p
les
la
meilleures
qu'il
conditions.
accordée
Ce
la
fut
en
un
p
réel
an
plaisir
ma
d'a
au
v
de
oir
Je
pu
'
tra
pas
v
Rouzaud
ailler
m'
a
initialemen
v
c
ec
p
eux
eectuer
et
thèse
d'a
Jean-Luc
v
t
oir
our
pu
encadremen
proter
certes
de
distance,
leur
p
sa
t
v
t.
oir.
remercie
P
s
our
bres
tout
jury
e
thèse
s
Pierre
ces
ervisc
raisons,
qui
je
fait
leur
u
adr
de
e
l
sse
jury
a
qui
v
accepté
ec
v
un
Ch
p
tal
eu
la
de
harge
nostalgie
rapp
mes
la
plus
Serguey
sincères
o
remerciemen
qui
ts.
accepté
Je
f
tiens
le
à
t
remercier
Moscou
l'ONERA
our
p
examinateur
our
que
son
e
sout
ran.
i
une
en
ensée
nancier,
p
mais
l'ensem
aussi
du
surtout
ersonnel
p
lab
our
EM2C,
son
ctoran
p
s,
ersonnel
ost-do
qui
ts
m'a
p
égalemen
t
t
qui
accompagné
t
p
tribué
endan
créer
t
am
mes
de
tra
v
v
c
aux.
reuse
Je
con
p
le.
ense
un
en
en
premier
particulier
lieu

à

Lionel
lab
T
a
essé
ec
qui
les
s'est
du
b
de
eaucoup
on
in
été
v
sympathiques.
esti
dernier
p
t
endan
our
t
famille
la
mes
dernière
d'enfance
année
on
de
toujours
ma
présen
thèse
au
et
de
qui
études
m'a
je
app
remercie
orté
v
un
fait
soutien
déplacemen
scien
p
tique
ma
indisp
de
ensable
un
p
ô
our
bien
ac
ortan
hev
p
er
tra
mes
Remerciemen
tel-00453309, version 1 - 4 Feb 2010tel-00453309, version 1 - 4 Feb 2010sur
n
un
L
Des
e
l'ordre
dimensionnemen
(MSBE)
t
mise
des
FIRE
protections
tangen
thermiques
Le
des
ensuite
v
mo
éhicules
ua,
spatiaux
de
néc
L
essite
h
la
rapp
connais-
le
sance
résultats
des
emen
ux
calculs
de
base
c
optiquemen
haleur
et
à
ectroniques
la
raie
surface
h
de
s
l'engin.
tribution
P
été
our
des
le
le
t
Carlo.
yp
de
e
onnemen
de
Les
ren
p
trée
tangen
atmosphérique
p
étudié
en
dans
visager
c
i
es
é
tra
statistique
v
systèmes
aux
(dans
(ren
co
trée
ennés
terrestre
our
à
et
10-12
c
km/s
transitions
ou
cette
plus),
une
le
tenir
ra
des
y
du
onnemen
on
t
h
con
le
tri-
r
bue
raie.
de
té
manière
ASTRE,
sig
he
n
été
i
métho
cativ
sur
e
3D
à
ensuite
l'éc
FIRE
hauemen
t
t
sur
pariétal.
par
Nous
l'appro
présen
utilisée
tons
-
dans
des
cette
v
étude,
tri-dimensionnelle
le
p
dév
c
e
radi
lopp
es
emen
dév
t
l
et
mo
la
bandes
mise
our
en
lectroniques
÷uvr
non
e
application),
de
e
mo
ts
d