Distributed Roughness-Induced Transition on a Blunt Body at Mach 6: a Numerical Investigation

Deze studie presenteert de eerste directe numerieke simulatie van gedistribueerde ruwheid op een Mach 6 stompe cilinder, waarbij wordt onthuld dat de arrangement van de ruwheid het transitiemechanisme bepaalt door ofwel sinusoïdale streak-modi of 2D T-S golven te bevorderen, waarbij laatstgenoemde uniek wordt in stand gehouden door een interne akoestische feedbacklus die de noodzaak voor externe forcering elimineert.

De kern

Stel je een ruimtevaartuig voor dat met zes keer de snelheid van het geluid door de atmosfeer raast. Om de intense hitte te overleven, is het oppervlak bedekt met een speciaal materiaal dat langzaam wegbrandt (ablatie) om het schip te beschermen. Echter, terwijl dit materiaal wegbrandt, laat het geen perfect glad oppervlak achter, maar een bobbelige, ruwe textuur, een beetje zoals schuurpapier.

Dit artikel is een hoogwaardige computersimulatie die een eenvoudige maar cruciale vraag stelt: Hoe veranderen deze kleine bobbels op het oppervlak een gladde, geordende luchtstroom in een chaotische, turbulente luchtstroom?

Hier is de uitsplitsing van hun bevindingen met behulp van alledaagse analogieën:

1. De Opstelling: De "Schuurpapier" Cilinder

De onderzoekers bouwden een digitaal model van een stompe cilinder (zoals de neus van een raket) die met Mach 6 vliegt. In plaats van een glad oppervlak, bedekten ze het met kleine, kunstmatige "bobbels" (ruwheid) om de zandachtige textuur na te bootsen die door het brandende materiaal wordt achtergelaten.

Ze testten drie verschillende manieren om deze bobbels te rangschikken:

• Uitgelijnd: Zoals soldaten die in perfecte rijen en kolommen staan.
• Gestapeld (Staggered): Zoals een bakstenen muur, waarbij de bobbels in één rij verschoven zijn ten opzichte van de rij erachter.
• Willekeurig: Zoals kiezels die zonder patroon op een stoep zijn verspreid.

2. De Oude Theorie vs. De Nieuwe Ontdekking

Lange tijd dachten wetenschappers dat de overgang naar turbulentie werd veroorzaakt door een "langzame opbouw" van energie, vergelijkbaar met hoe een schommel steeds hoger wordt als je haar op het juiste moment een duwtje geeft. Dit wordt "transiënte groei" genoemd.

De Bevinding van het Papier:
De simulatie toonde aan dat deze "langzame opbouw"-theorie niet echt verklaart wat hier gebeurt. De bobbels op het oppervlak versterkten de energie niet alleen langzaam; ze fungeerden als destabilisatoren. Ze maakten de luchtstroom onmiddellijk instabiel en veranderden deze in een specifiek type golf dat zeer snel groeit.

Denk er zo over na: de oude theorie dacht dat de bobbels voorzichtig een dominosteen een zetje gaven om om te vallen. De nieuwe ontdekking laat zien dat de bobbels de domino eigenlijk een trap geven, waardoor hij onmiddellijk tegen de volgende klapt.

3. De Twee Soorten "Golven"

Afhankelijk van hoe de bobbels waren gerangschikt, reageerde de luchtstroom op twee verschillende manieren:

• De "Slang" (Sinuze Modus): Wanneer de bobbels uitgelijnd waren (perfecte rijen), begon de luchtstroom zijwaarts te wiebelen als een slang. Dit is een zeer specifieke, georganiseerde wiebel.
• De "Platte Golf" (Tollmien-Schlichting of T-S Golven): Wanneer de bobbels gestapeld of willekeurig waren, begon de luchtstroom te rimpelen in een plat, 2D-golfpatroon. Dit is een klassiek type golf dat meestal wordt gevonden in veel langzamere, laagfrequente luchtstromen, wat verrassend was om in deze hogesnelheidsomgeving te vinden.

Het Belangrijke Inzicht: De rangschikking van de bobbels bepaalde welke "dans" de lucht zou uitvoeren. De "slangendans" gebeurde bij uitgelijnde bobbels, terwijl de "platte golfdans" plaatsvond bij de andere configuraties.

4. Het "Haarspeld" Slotstuk

Zodra deze golven sterk genoeg waren, veroorzaakten ze de laatste fase van de crash. De constante "strepen" in de luchtstroom die door de bobbels werden gecreëerd (die lijken op lange, onzichtbare linten van langzame lucht) draaiden plotseling en knapten in haarspeldvortices.

Stel je een rubberen band voor die strak gespannen is. Plotseling draait deze en vormt een lus die op een haarspeld lijkt. Deze lussen zijn de geboorte van turbulentie. Zodra deze haarspelden ontstaan, breekt de gladde luchtstroom volledig af in chaos, en schiet de hitte op het oppervlak van het ruimtevaartuig dramatisch omhoog.

5. De "Echo-kamer" Verrassing

Een van de meest fascinerende ontdekkingen was hoe de turbulentie in eerste instantie ontstond voor de gestapelde en willekeurige gevallen.

Normaal gesproken denken wetenschappers dat je een externe "duw" nodig hebt (zoals een windvlaag of een trilling) om deze golven te starten. Maar de simulatie toonde iets zelfvoorzienends aan:

1. Turbulentie begint op één plek op de cilinder.
2. Omdat de lucht achter de schokgolf langzamer beweegt dan het geluid (subsonisch), reist het geluid van deze turbulentie terugwaarts stroomopwaarts als een echo.
3. Deze "echo" raakt het gladde deel van het oppervlak vóór de turbulentie en brengt de lucht daar in beweging, waardoor er nieuwe turbulentie ontstaat.
4. Dit creëert een feedbackloop: Turbulentie maakt geluid, het geluid reist terug, en het geluid creëert meer turbulentie.

Het is alsof een microfoon zijn eigen luidspreker opvangt en een gierend feedback-geluid produceert, maar in dit geval is de "gier" de lucht die turbulent wordt.

Samenvatting

Dit artikel gebruikte een supercomputer om de luchtstroom over een bobbelige, hogesnelheidscilinder te observeren. Ze ontdekten dat:

• Het patroon van de bobbels bepaalt precies hoe de lucht turbulent wordt.
• Het oude idee van "langzame energieopbouw" niet de hoofdschuldige is; in plaats daarvan destabiliseren de bobbels direct specifieke golven.
• Deze golven groeien totdat ze draaien in "haarspeldvormen", wat ervoor zorgt dat de lucht chaotisch wordt.
• In sommige gevallen creëert de turbulentie zijn eigen "echo" die terugreist om het proces opnieuw te starten, zonder dat daarvoor externe hulp nodig is.

Dit helpt ingenieurs begrijpen dat de kleine, willekeurige bobbels die achterblijven door het brandende hitteschild niet slechts kleine imperfecties zijn; ze zijn de primaire architecten van hoe en wanneer het oppervlak van een ruimtevaartuig gevaarlijk heet wordt.

Technische samenvatting

Technische Samenvatting: Gedistribueerde Ruwheid-geïnduceerde Transitie op een Stomp Lichaam bij Mach 6

Probleemstelling
De voorspelling van laminaire-turbulente transitie (LTT) over snelle, stompe geometrieën blijft een kritieke uitdaging voor het technisch ontwerp, met name voor voertuigen die de atmosfeer binnenkomen waarbij oppervlakte-ablatie zorgt voor gedistribueerde ruwheid. Hoewel er empirische correlaties bestaan (bijv. Hollis 2017, 2019, 2025), missen deze een duidelijke fysieke fundering. Eerdere theoretische inspanningen, zoals de Reshotko & Tumin (2004) correlatie, vertrouwen op mechanismen van transiënte groei maar kampen met beperkingen: ze gaan uit van een lineaire schaling van perturbaties met de ruwheidshoogte (geldig voor slechts kleine $k/\delta$), vereisen aanpassingsconstanten op basis van experimentele data, en hebben moeite om transitie te verklaren in gevallen waar de ruwheidshoogtes de dikte van de grenslaag overstijgen ($k/\delta > 1$). Bovendien simuleerden eerdere numerieke onderzoeken vaak ruwheidspatches in plaats van volledig gedistribueerde ruwheid, of slaagden ze er niet in om modale instabiliteitsversterking te observeren, waardoor de specifieke mechanismen die de transitie drijven in volledig gedistribueerde, hypersonische stromingen rond stompe lichamen onduidelijk bleven.

Methodologie
Deze studie presenteert de eerste Direct Numerical Simulation (DNS) van een stomme cilinder in Mach 6 cross-flow met ruwheidselementen die over het gehele oppervlak zijn gedistribueerd. De simulaties maken gebruik van de CHAMPS-solver om de samendrukbare Navier-Stokes vergelijkingen (CNSE) voor een ideaal gas te integreren. De geometrie is een cilinder met een diameter van 0,1524 m, gemodelleerd op basis van experimentele condities van Hollis (2017) (Run 38, Test 6975), met een vrije stroom Mach-getal van 5,98.

De ruwheid wordt gemodelleerd als een sinusvormige distributie met een vaste amplitude ($A \approx 35\%$ van de dikte van de grenslaag) en een golflengte die overeenkomt met nominale zandkorrelruwheid. Drie verschillende breedtevolgorde-faseringsconfiguraties werden onderzocht:

1. Gestapeld (Staggered): Rijen zijn $180^\circ$ uit fase.
2. Uitgelijnd (Aligned): Rijen zijn in fase ($0^\circ$ verschuiving).
3. Willekeurig (Random): Rijen hebben willekeurig verdeelde faseverschuivingen.

Om het transitieproces te analyseren, hanteerden de auteurs:

• DNS: Om de volledige niet-lineaire evolutie van de stroming te vangen.
• Lineaire Verstoringenvergelijkingen (LDE): Opgelost op de tijdgemiddelde gemiddelde stromingen uit de DNS om instabiliteitsmechanismen te isoleren middels puls-forcing.
• Lineaire Stabiliteitstheorie (LST): Toegepast op breedtegemiddelde profielen voor vergelijking.
• Spectrale Analyse: Fast Fourier Transforms (FFT) van snelheid- en drukvelden om dominante frequenties en golfgetallen te identificeren.

Belangrijkste Resultaten

• Gemiddelde Stroming en Streaks: Alle configuraties genereerden stromingsstreaks in de stroomrichting via een lift-up mechanisme, gedreven door tegenrotende vortexparen (CVP's). De uitgelijnde configuratie produceerde de sterkste stationaire streaks door een consistente fasering die vloeistof met een lage massaflux hoger in de grenslaag bracht. De gestapelde en willekeurige configuraties produceerden zwakkere, complexere streak-patronen. In tegenstelling tot de voorspellingen van de optimale transiënte groeitheorie, vertoonden de streak-amplitudes een lage relatieve groei (minder dan een factor 2) en bleven ze door het hele domein niet-lineair, wat het nut van lineaire transiënte groei als voorspellend instrument voor dit regime uitdaagt.

• Instabiliteitsmechanismen: Het transitieproces werd gedreven door convectieve instabiliteiten, maar het specifieke type modus hing af van de ruwheidsfasering:

  • Uitgelijnd geval: Gedomineerd door een fundamentele sinusoïdale streak-modus. De Tollmien-Schlichting (T-S) golven werden zwak gedestabiliseerd.
  • Gestapelde en Willekeurige gevallen: Gedomineerd door 2D T-S golven (varicose modus). Deze configuraties vertoonden een significante destabilisatie van T-S golven, wat leidde tot snelle exponentiële groei.
  • Breakdown: In alle gevallen, zodada de dominante instabiliteitsgolven een voldoende amplitude bereikten, triggerden zij de stationaire streaks om hairpin-vortices te vormen, wat leidde tot breakdown en turbulentie.

• Transitielocaties: De DNS voorspelde transitielocaties ($\Theta_{tr}$) aanzienlijk verder stroomafwaarts (ca. $11^\circ$–$17^\circ$) dan de meest geavanceerde empirische correlaties (die $6^\circ$–$13^\circ$ voorspelden). Het uitgelijnde geval vertoonde de vroegste transitie onder de ruwheidscasussen vanwege de benodigde forcing om LTT te initiëren, terwijl de gestapelde en willekeurige gevallen later transiteerden.

• Feedbackmechanisme: Een nieuwe bevinding in de gestapelde en willekeurige gevallen was de identificatie van een feedbackloop. Turbulentie gegenereerd in het subsonische deel van de cilinder verstrooide akoestische golven die stroomopwaarts voortplantten. Deze trage akoestische golven exciteerden T-S golven aan de tegenoverliggende zijde van de cilinder (over de stagnatie-lijn heen) bij hun neutrale punt. Dit mechanisme, dat geen externe forcing vereist, houdt het transitieproces in stand.

Betekenis en Claims
De auteurs beweren dat dit werk een nieuw fysiek begrip biedt van ruwheid-geïnduceerde transitie op hypersonische stompe lichamen, waarbij de huidige afhankelijkheid van transiënte groeitheorie wordt uitgedaagd. Belangrijkste bijdragen zijn:

1. Verwerping van Transiënte Groei als Primaire Drijver: De studie toont aan dat voor volledig gedistribueerde ruwheid met $k/\delta \approx 0,35$, de niet-lineaire ontwikkeling van streaks en de lage amplificatieverhoudingen die worden waargenomen, de optimale transiënte groeitheorie een onvoldoende voorspeller maken van het werkelijke transitieproces.
2. Destabilisatie van Modale Instabiliteiten: Het artikel stelt vast dat gedistribueerde ruwheid effectief modale instabiliteiten (specifiek T-S golven) kan destabiliseren, zelfs op stompe lichamen waar de gladde wand-stroming typisch eigenmode-stabiel is. De ruwheidsconfiguratie bepaalt of de transitie wordt gedreven door streak-modi of T-S golven.
3. Ontdekking van een Akoestische Feedbackloop: De identificatie van een zelfonderhoudend feedbackmechanisme waarbij door turbulentie gegenereerde akoestiek T-S golven aan de andere zijde van het lichaam exciteert, biedt een potentiële verklaring voor transitie in subsonische regio's van hypersonische stromingen, een fenomeen dat voorheen onverklaarbaar was door lineaire mechanismen.
4. Gevoeligheid voor Fasering: De resultaten benadrukken de extreme gevoeligheid van LTT voor de relatieve fasering van ruwheidselementen, wat suggereert dat eenvoudige topologieën het gedrag van complexere, willekeurige ruwheidspatronen kunnen begrenzen.

De auteurs concluderen dat hoewel het voorgestelde akoestische feedbackmechanisme verdere rigoureuze bewijsvoering vereist, de blootstelling van deze instabiliteitsmechanismen via hoogwaardige DNS wijst op een fundamentele verschuiving in het begrip van de verstoringen die verantwoordelijk zijn voor LTT op hypersonische stompe lichamen.