10.01.2015 Views

Simulation numérique directe de la turbulence en présence d ... - ISAE

Simulation numérique directe de la turbulence en présence d ... - ISAE

Simulation numérique directe de la turbulence en présence d ... - ISAE

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

<strong>Simu<strong>la</strong>tion</strong> <strong>numérique</strong> <strong>directe</strong> <strong>de</strong> <strong>la</strong> <strong>turbul<strong>en</strong>ce</strong> <strong>en</strong> prés<strong>en</strong>ce<br />

d’une paroi ab<strong>la</strong>table<br />

Résumé<br />

Lorsqu’une son<strong>de</strong> r<strong>en</strong>tre dans l’atmosphère à une vitesse hypersonique, le bouclier thermique<br />

reçoit plusieurs c<strong>en</strong>taines <strong>de</strong> MW.m −2 . Ce flux est absorbé par <strong>la</strong> paroi (composite C/C) grâce à<br />

<strong>de</strong>s réactions physico-chimiques (oxydation, sublimation, etc.). Les expéri<strong>en</strong>ces <strong>de</strong> jet <strong>de</strong> p<strong>la</strong>sma<br />

révèl<strong>en</strong>t que cette perte <strong>de</strong> matière s’accompagne <strong>de</strong> l’apparition <strong>de</strong> rugosités spécifiques lorsque<br />

le régime <strong>de</strong> l’écoulem<strong>en</strong>t est turbul<strong>en</strong>t. Dans ce contexte, cette thèse a pour objectif <strong>de</strong> caractériser,<br />

à <strong>la</strong> plus petite échelle <strong>de</strong> <strong>la</strong> <strong>turbul<strong>en</strong>ce</strong> (échelle <strong>de</strong> Kolmogorov), l’interaction <strong>en</strong>tre l’écoulem<strong>en</strong>t<br />

turbul<strong>en</strong>t et une paroi ab<strong>la</strong>table afin <strong>de</strong> mieux évaluer <strong>la</strong> formation <strong>de</strong> ces rugosités.<br />

La première phase <strong>de</strong> ce travail a ainsi consisté <strong>en</strong> <strong>la</strong> validation physique <strong>de</strong>s phénomènes turbul<strong>en</strong>ts<br />

simulés, d’abord <strong>en</strong> configuration périodique puis <strong>en</strong> prés<strong>en</strong>ce <strong>de</strong> surfaces <strong>de</strong> blocage<br />

pariétal. Pour ce<strong>la</strong>, nous avons développé et adapté à une configuration <strong>de</strong> prés<strong>en</strong>ce <strong>de</strong> paroi,<br />

une métho<strong>de</strong> <strong>de</strong> forçage spectral <strong>de</strong> <strong>la</strong> <strong>turbul<strong>en</strong>ce</strong> dont nous validons l’implém<strong>en</strong>tation. Pour<br />

ce<strong>la</strong>, un traitem<strong>en</strong>t statistique <strong>de</strong>s résultats a été é<strong>la</strong>boré pour pr<strong>en</strong>dre <strong>en</strong> compte le caractère<br />

aléatoire <strong>de</strong> <strong>la</strong> <strong>turbul<strong>en</strong>ce</strong>. L’étu<strong>de</strong> <strong>de</strong> l’évolution structurelle du matériau ab<strong>la</strong>té, couplé à un<br />

champ <strong>de</strong> vitesse turbul<strong>en</strong>t, démontre alors l’influ<strong>en</strong>ce du nombre <strong>de</strong> Reynolds <strong>de</strong> <strong>turbul<strong>en</strong>ce</strong> et<br />

<strong>de</strong> <strong>la</strong> taille <strong>de</strong>s structures porteuses d’énergie, sur <strong>la</strong> vitesse <strong>de</strong> récession <strong>de</strong> <strong>la</strong> paroi et les motifs<br />

rugueux observés.<br />

Mots clés : Ab<strong>la</strong>tion, DNS, Forçage, Turbul<strong>en</strong>ce<br />

Abstract<br />

Direct numerical simu<strong>la</strong>tion of <strong>turbul<strong>en</strong>ce</strong> near<br />

an ab<strong>la</strong>table wall<br />

During the atmospheric re-<strong>en</strong>try the heatshield of the probe suffers a significant overheating.<br />

The composite material un<strong>de</strong>rgoes an ab<strong>la</strong>tive process that consumes the heat flux by physicochemical<br />

reactions (sublimation, oxidation, etc.) and induces the wall disappearance. In the case<br />

of polycrystalline graphite, p<strong>la</strong>sma jet experim<strong>en</strong>ts reveal a structural surface roughness like a<br />

scalloped pattern wh<strong>en</strong> the flow is turbul<strong>en</strong>t. The main i<strong>de</strong>a of this study is to characterize the<br />

interaction betwe<strong>en</strong> a turbul<strong>en</strong>t flow and an ab<strong>la</strong>table material at the Kolmogorov scale (the<br />

smallest turbul<strong>en</strong>t scale) in or<strong>de</strong>r to analyze roughness appearance. First, we have validated the<br />

<strong>turbul<strong>en</strong>ce</strong> simu<strong>la</strong>tion for an infinite domain (periodic boundaries) and for the case of a solid<br />

wall. To this ext<strong>en</strong>t, we have <strong>de</strong>veloped a new way to force <strong>turbul<strong>en</strong>ce</strong> in the spectral space which<br />

is adapted to the solid wall configuration. To this ext<strong>en</strong>t, we have e<strong>la</strong>borated a statistic treatm<strong>en</strong>t<br />

of results as <strong>turbul<strong>en</strong>ce</strong> is an unpredictable ph<strong>en</strong>om<strong>en</strong>a. Th<strong>en</strong>, the evolution of the ab<strong>la</strong>table<br />

surface in interaction with a turbul<strong>en</strong>t flow proves that Reynolds number, wav<strong>en</strong>umber κ e and<br />

flow temperature have a strong influ<strong>en</strong>ce on wall recession velocity and roughness appearance.<br />

Keywords : Ab<strong>la</strong>tion, DNS, Forcing, Turbul<strong>en</strong>ce

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!