Fast Prediction of Aero-Thermal Performance During Hypersonic Vehicle Design

Technische Samenvatting: Snelle voorspelling van aero-thermische prestaties tijdens het ontwerp van hypersonische voertuigen

Probleemstelling
Het ontwerp van hypersonische voertuigen vereist een snelle schatting van aerodynamische krachten en oppervlaktewarmtefluxen tijdens de vroege ontwikkelingsfasen. Hoewel hoogwaardige Computational Fluid Dynamics (CFD)-simulaties gedetailleerde nauwkeurigheid bieden, zijn ze computationeel te kostbaar voor de uitgebreide ontwerpitereaties en trajectoptimalisaties die nodig zijn in de voorbereidende fasen. Traditionele engineeringmethoden, zoals technieken op basis van lokale oppervlaktehelling (tangent-wedge, gemodificeerde Newtoniaanse methode, Prandtl-Meyer), bieden snelheid maar missen vaak de getrouwheid bij het modelleren van viskeuze effecten. Specifiek gaan standaardbenaderingen er meestal van uit dat de grenslaag ofwel volledig laminair ofwel volledig turbulent is. De aanname van een volledig turbulente stroming, hoewel conservatief, leidt vaak tot een overschatting van de wrijvingsweerstand en de warmteflux aan het oppervlak. Dit resulteert in overgedimensioneerde thermische beschermingssystemen (TPS), een verhoogde voertuigmassa en een hoger brandstofverbruik. Bovendien bevatten weinig snelle engineeringtools momenteel effecten van grenslaagtransitie binnen lokale hellingskaders, ondanks de significante impact van transitie-stromingen op aerothermische belastingen.

Methodologie
Het artikel presenteert HipeX, een lichtgewicht aerothermodynamische analyse-tool die is ontworpen om de kloof te overbruggen tussen klassieke analytische methoden en hoogwaardige CFD. HipeX maakt gebruik van een hybride engineering-niveau benadering gebaseerd op lokale oppervlaktehellingsmethoden om inviscidische stromingseigenschappen te schatten en hanteert een geblend grenslaagmodel om viskeuze effecten en transitie aan te pakken.

• Inviscidische Stromingsmodellering: De solver verwerkt ongestructureerde oppervlaktemeshes (STL-formaat) en selecteert het passende stromingsmodel op basis van de lokale oppervlaktehelling ten opzichte van de invallende snelheidsvector:

  • Gemodificeerde Newtoniaanse Theorie: Toegepast voor losgelaten schokken (stompere lichamen) en grote deflectiehoeken.
  • Oblique Schokrelaties: Gebruikt voor aangelegde schokken en compressieregio's.
  • Prandtl-Meyer Expansie: Toegepast op stromingsexpansieregio's, met een veiligheidsmaatregel die de toepassing beperkt tot 50% van de theoretische maximale deflectiehoek om numerieke stabiliteit te garanderen.
  • Gaubeaud-formule: Gebruikt als empirische fallback wanneer Prandtl-Meyer oplossingen falen.
  • Atmosferische eigenschappen worden bepaald via standaardmodellen op basis van de vlieghoogte.

• Viskeuze en Thermische Modellering:

  • Compressibiliteit: Gecorrigeerd met de Eckert referentietemperatuurmethode, toepasbaar op zowel laminaire als turbulente regimes.
  • Transitiemodellering: Een kritisch Reynoldsgetal ($Re_{cr}$) wordt empirisch berekend als een functie van het getal van Mach om het begin van de transitie te bepalen. Het model gebruikt een intermittentiefactor ($\gamma_{it}$) variërend van 0 (laminair) tot 1 (turbulent) om de huidwrijvingscoëfficiënten van laminaire en turbulente stromingen over een gedefinieerd transitiegebied te mengen.
  • Berekening Warmteflux: De oppervlaktewarmteflux wordt afgeleid van de huidwrijvingscoëfficiënt met behulp van de Reynolds-analogie en het Stantongetal. Het model houdt rekening met zowel convectieve warmteoverdracht als radiatieve warmteoverdracht (uitgaande van een emissiviteit van 0,4), waarbij de netto warmteflux en de totale thermische last worden berekend door integratie over het bevochte oppervlak.

• Workflow: De tool accepteert 3D-geometrie en vluchtcondities (Mach, hoogte, invalshoek, zijwaartse slip). Het berekent de inviscidische eigenschappen, past viskeuze correcties toe, integreert krachten om aerodynamische coëfficiënten te bepalen en exporteert resultaten in VTK-formaat. De code ondersteunt parallelle uitvoering voor trajectanalyse.

Belangrijkste Bijdragen

1. Geblend Transitiemodel: De primaire bijdrage is de integratie van een transitie-grenslaagmodel in een snelle, lokale hellingsmethode. Dit stelt de solver in staat om de verschuiving van laminaire naar turbulente stroming te schatten op basis van lokale Reynoldsgetallen en Machgetallen, in plaats van te vertrouwen op binaire laminaire/turbulente aannames.
2. Computationele Efficiëntie: HipeX is ontworpen voor snelle uitvoering en is in staat om de prestaties voor een volledig traject (ongeveer 300 punten) te berekenen op een standaard werkstation in ongeveer één minuut.
3. Validatie tegen Benchmarkgegevens: De tool is gevalideerd tegen CFD-data en experimentele referenties voor de HEXAFLY-INT configuratie (een blended lifting-body ontwerp voor Mach 7–8 cruise).

Resultaten
Validatie werd uitgevoerd met de HEXAFLY-INT configuratie over diverse Mach-getallen (M=2 tot M=7,5) en vluchttrajecten.

• Aerodynamische Coëfficiënten: Voor hoge Mach-getallen (M=6 en M=7) komen de voorspellingen van HipeX voor lift ($C_L$), weerstand ($C_D$) en drukcentrum ($X_{cp}$) nauw aan bij de CFD-data. De geblend transitieoplossing biedt over het algemeen nauwkeurigere resultaten dan puur laminaire of volledig turbulente aannames. Bij lagere Mach-getallen (M=4 en M=2) nemen de discrepanties toe, met name in de absolute $C_L$ en $C_D$ waarden, hoewel de lift-naar-weerstandverhouding ($L/D$) redelijk nauwkeurig blijft. De tool identificeert correct de grenzen van lokale hellingsmethoden wanneer de stroming afwijkt van sterk hypersonische condities.
• Thermische Lasten: In gevallen waar de stroming overwegend laminair is (hoge hoogte), levert de geblend oplossing resultaten op die vergelijkbaar zijn met het volledig laminaire geval. Echter, in transitie-regimes vangt het geblend model tussenliggende warmtefluxniveaus op die tussen de laminaire en turbulente extremen in liggen.
• TPS-Implicaties: De studie toont aan dat het aannemen van volledig turbulente stroming leidt tot een significante overschatting van de verhitting, vooral nabij de neus (leading edge). Door het laminaire of gedeeltelijk turbulente karakter van de grenslaag over grote delen van het traject te vangen, levert HipeX lagere, meer realistische schattingen van de warmteflux. Dit suggereert potentieel voor het verminderen van de TPS-massa zonder de veiligheidsmarges in gevaar te brengen.
• Transitie-onset: Vergelijkingen met de BLITZ-code laten vergelijkbare trends zien in de locatie van de transitie, hoewel HipeX de neiging heeft om een eerder begin van turbulentie te voorspellen. Het model heeft momenteel moeite om het voorkomen van relaminarisatie in bepaalde expansieregio's te voorkomen.

Betekenis en Claims
Het artikel positioneert HipeX als een robuuste basis voor vroege fase aero-thermische beoordeling en ontwerpoptimalisatie. De betekenis ligt in het bieden van een "gebalanceerde afweging" tussen ontwerpconservatisme en missie-efficiëntie. Door af te stappen van de overly conservatieve aanname van volledig turbulente stroming, stelt de tool in staat om nauwkeurigere dimensies voor thermische beschermingssystemen te bepalen, wat potentieel de voertuigmassa en het brandstofverbruik kan verminderen.

De auteurs erkennen dat HipeX de inherente beperkingen van lokale oppervlaktehellingsmethoden overneemt, specifificamente het onvermogen om schokinteracties, scheidingszones en sterke driedimensionale effecten te voorspellen. Daarnaast blijft de nauwkeurige bepaling van het begin van de transitie onder werkelijke vluchtcondities uitdagend vanwege de afhankelijkheden van oppervlakte-ruwheid en verstoringen in de vrije stroom. Desalniettende slaagt de tool erin de kloof te vullen tussen klassieke analytische methoden en hoogwaardige CFD, waardoor snelle evaluatie van volledige trajecten en grote verzamelingen geometrische configuraties mogelijk wordt, wat essentieel is voor het voorlopige ontwerp.