Scale-resolving simulations and data-driven modal analysis of turbulent transonic buffet cells on infinite swept wings

Dit onderzoek toont aan dat transsonische buffet op oneindig geschuinde vleugels voortkomt uit de superpositie van een quasi-tweedimensionale schokbeweging en door gemiddelde stroomscheiding gedreven driedimensionale instabiliteiten, waarbij de scheiding bij de schok een noodzakelijke voorwaarde is voor het ontstaan van dominante 3D-buffetdynamiek.

De kern

De Dans van de Schokgolf: Waarom Vliegtuigen Trillen op Transonische Snelheid

Stel je voor dat je vliegtuig vliegt op een snelheid die net onder het geluidsniveau ligt. Op dat moment gebeurt er iets vreemds op de vleugels: er ontstaat een onzichtbare muur van lucht, een schokgolf. Normaal gesproken is deze golf vrij stabiel, maar bij bepaalde omstandigheden begint hij te trillen, heen en weer te bewegen alsof hij een onrustige danser is. Dit fenomeen heet transonische buffet.

Het probleem? Deze trillingen kunnen het vliegtuig zo hard laten schudden dat het structureel beschadigd raakt. Vroeger dachten wetenschappers dat dit vooral een tweedimensionaal probleem was: de schokgolf bewoog alleen maar vooruit en achteruit, als een slinger. Maar nieuw onderzoek toont aan dat het veel complexer is.

Hier is wat David Lusher en Andrea Sansica van JAXA (het Japanse ruimtevaartagentschap) hebben ontdekt, vertaald naar begrijpelijke taal:

1. Het Probleem: Te Korte Vleugels in de Computer

Om dit te bestuderen, gebruiken wetenschappers supercomputers. Maar tot nu toe waren deze simulaties beperkt tot heel smalle vleugels (in de computerwereld).

• De Analogie: Stel je voor dat je probeert een danspartij te analyseren, maar je kijkt alleen door een smalle brievenbus. Je ziet misschien dat iemand met zijn voeten stampt (de schokgolf), maar je mist de hele dansvloer. Je ziet niet dat er ook mensen in cirkels om elkaar heen dansen (de 3D-buffetcellen).
• De Oplossing: Deze onderzoekers hebben een computermodel gebouwd met een vleugel die drie keer zo breed is als de lengte van het profiel. Ze hebben ook de vleugel "geschud" (een hoek gegeven, wat sweep of 'sweeping' wordt genoemd), zoals bij moderne passagiersvliegtuigen.

2. De Twee Soorten Dansers

Het onderzoek toont aan dat er twee verschillende soorten "dansers" (instabiliteiten) zijn die tegelijkertijd kunnen optreden:

• De Slinger (2D-modus): Dit is de klassieke schokgolf die vooruit en achteruit beweegt. Hij is als een oude, betrouwbare pendel. Deze beweegt altijd met dezelfde snelheid, ongeacht hoe de vleugel is geschud.
• De Dansende Cellen (3D-modus): Dit is het nieuwe, spannende deel. Wanneer de luchtstroming achter de schokgolf erg onrustig wordt (wat gebeurt als de vleugel steiler staat), ontstaan er bellen van losse lucht die over de vleugel heen rollen.
  • De Analogie: Denk aan een waslijn met natte lakens. Als de wind stil is, hangen ze recht. Maar als de wind opsteekt, beginnen de lakens te golven en te draaien. Deze "golven" (de buffetcellen) bewegen zijwaarts over de vleugel, van links naar rechts, met een eigen ritme dat sneller is dan de slinger.

3. De Magische Formule: De Hoek maakt het Verschil

De onderzoekers keken naar twee situaties:

1. Zachtjes vliegen (kleine hoek): De lucht stroomt vrij glad. De "slinger" (2D) is de baas. De "dansende cellen" (3D) zijn er nauwelijks, of ze zijn zo zwak dat je ze niet ziet. Het vliegtuig trilt, maar het is een voorspelbare trilling.
2. Steiler vliegen (grote hoek): De lucht stroomt chaotischer en gaat loslaten achter de schokgolf. Nu komen de dansende cellen echt tot leven.
   • Het effect van de hoek (Sweep): Als je de vleugel schuin zet (zoals bij een vliegtuig), verandert het gedrag van deze cellen drastisch. Ze beginnen niet meer stil te staan, maar rennen over de vleugel.
   • De snelheid: Hoe steiler de vleugel staat, hoe sneller deze cellen over de vleugel rennen en hoe hoger hun "toon" (frequentie) wordt. Het is alsof je een trommelstok harder slaat: het ritme gaat sneller.

4. Waarom is dit belangrijk?

Vroeger dachten ingenieurs dat ze alleen de "slinger" (de 2D-trilling) hoefden op te lossen. Maar dit onderzoek bewijst dat op echte vliegtuigen (die schuine vleugels hebben), de dansende cellen vaak de echte boosdoener zijn.

• Als de luchtstroming erg onrustig is, wordt deze 3D-dans zo sterk dat hij de oude 2D-slinger volledig overstemt.
• Dit verklaart waarom vliegtuigen in de praktijk soms heel anders reageren dan de simpele computermodellen voorspellen. De "slinger" verdwijnt bijna en de "dansende cellen" nemen het over.

Samenvatting in één zin

Dit onderzoek laat zien dat transonische buffet niet alleen een voor- en achterwaartse trilling is, maar een complex dansfeest waarbij, afhankelijk van de snelheid en de hoek van de vleugel, zijwaartse luchtstromen (de cellen) kunnen ontstaan die het vliegtuig in een ander, vaak gevaarlijker ritme laten trillen.

De les voor de toekomst: Als je een vliegtuig wilt ontwerpen dat veilig is, moet je niet alleen kijken naar de schokgolf zelf, maar ook naar hoe de luchtstroming eromheen "dansen" kan gaan, vooral als de vleugels schuin staan.

Technische samenvatting

Probleemstelling

Transsonische buffet is een vorm van zelfonderhoudende interactie tussen een schokgolf en een grenslaag (SBLI) die optreedt bij hoge subsonische Mach-getallen en matige tot hoge aanvalshoeken. Dit fenomeen veroorzaakt sterke, zelfopwekkende oscillaties in lift en weerstand, wat de structurele integriteit en het operationele vluchtenvelop van vliegtuigen kan bedreigen.

De literatuur onderscheidt twee instabiliteiten:

1. Tweedimensionale (2D) instabiliteit: Lage frequentie oscillaties van de hoofdschokgolf in de koordevlakrichting (Strouhal-getal $St \approx 0.05 - 0.1$).
2. Driedimensionale (3D) instabiliteit: Spanwijdte-gemoduleerde "buffet cells" (separatiebellen) die zich voortplanten langs de spanwijdte van de schokgolf bij hogere frequenties ($St \approx 0.2 - 0.4$).

Een groot probleem in het onderzoek is dat de meeste high-fidelity simulaties (zoals LES en DNS) beperkt blijven tot smalle vleugels met een lage aspectverhouding ($AR = L_z/c \approx 0.05 - 0.25$). Deze smalle domeinen zijn onvoldoende om de volledige 3D "buffet cell" instabiliteit te vangen, die een golflengte van ongeveer $1 - 1.5$ keer de koorde vereist. Daarnaast is de relatie tussen de 2D schokbeweging en de 3D separatie-instabiliteit, en de invloed van vleugelscheg (sweep), nog niet volledig gekwantificeerd in turbulente stromingen.

Methodologie

De auteurs hebben een reeks van tien zeer grote, high-fidelity Implicit Large Eddy Simulations (ILES) uitgevoerd op oneindige vleugels met een hoge aspectverhouding ($AR = 3$).

• Numerieke Oplossing: Gebruik werd gemaakt van de open-source solver OpenSBLI, een hoog-orde eind-differentie code die direct de ongestationaire compressibele Navier-Stokes vergelijkingen oplost zonder statistische gemiddelden of turbulentiemodellen (behalve de implicit filtering van de ILES methode).
• Rekenkracht: De simulaties omvatten ongeveer $8 \times 10^9$ mesh-punten per geval en draaiden op GPU- en CPU-clusters (JAXA en Fugaku), met een totale rekentijd van ongeveer 648.000 GPU-uren.
• Configuratie:
  • Profiel: NASA Common Research Model (CRM) vleugelprofiel.
  • Stromingscondities: $M_\infty = 0.72$, $Re_c = 0.5 \times 10^6$.
  • Variabelen: Twee aanvalshoeken ($\alpha = 5^\circ$ voor minimale gemiddelde separatie en $\alpha = 6^\circ$ voor grote gemiddelde separatie) en een reeks scheghoeken ($\lambda = 0^\circ, 5^\circ, 10^\circ, 15^\circ, 25^\circ, 35^\circ$).
• Data-analyse: De gegenereerde tijdsreeksen werden geanalyseerd met Spectrale Eigenwaarde Decompositie (SPOD). Deze data-gedreven methode werd gebruikt om dominante frequenties te identificeren en coherente stromingsstructuren te extraheren, met name om de 2D en 3D modi van elkaar te scheiden.

Belangrijkste Bijdragen

1. Schaal van de simulatie: Dit is een van de eerste studies die high-fidelity ILES/DNS toepast op oneindige vleugels met een aspectverhouding van $AR=3$, wat groot genoeg is om de volledige 3D buffet cell golflengte te accommoderen.
2. Koppeling van 2D en 3D modi: Het paper koppelt voor het eerst de intermitterende, quasi-stationaire 3D separatiemodus (geobserveerd op niet-geschegde vleugels) aan de reistende 3D buffet cell modus die wordt waargenomen op geschegde vleugels.
3. Rol van separatie: Het identificeert dat significante gemiddelde stromingsseparatie bij de schokgolf een noodzakelijke voorwaarde is voor het ontstaan van dominante 3D buffet dynamiek.
4. Invloed van scheg: Het kwantificeert systematisch hoe scheg de frequentie en aard van de 3D instabiliteit beïnvloedt, terwijl de 2D schokmodus grotendeels ongevoelig blijft voor scheg.

Resultaten

1. Geval met minimale separatie ($\alpha = 5^\circ$):

• De stroming gedraagt zich overwegend quasi-2D. De schokgolf oscilleert in de koordevlakrichting met een dominante frequentie van $St \approx 0.08$.
• Zelfs bij hoge scheghoeken ($\lambda = 25^\circ$) blijven de schokoscillaties spanwijdte-uniform.
• Er worden slechts zwakke, intermitterende separatiebelletjes waargenomen nabij de achterrand, maar deze interageren niet sterk met de hoofdschokgolf en veroorzaken geen macroscopische 3D buffet cells. De energie van eventuele 3D modi is hier minimaal (drie ordes van grootte lager dan de 2D modus).

2. Geval met grote separatie ($\alpha = 6^\circ$):

• Bij deze hoek is er een grote gemiddelde separatiebel bij de schokgolf. Hierdoor ontstaan duidelijke 3D buffet cells met een golflengte van $\lambda_z \approx 1 - 1.5c$.
• Invloed van scheg:
  • Bij lage scheg ($\lambda \le 10^\circ$) is de 3D modus nog zwak en verschijnt deze bij lage frequenties ($St \approx 0.02 - 0.1$).
  • Bij toenemende scheg ($\lambda = 15^\circ - 35^\circ$) verschuift de 3D modus monotoon naar hogere frequenties ($St = 0.16 - 0.35$).
  • De quasi-stationaire 3D modus van de niet-geschegde vleugel evolueert naar een spanwijdte-reistende modus op geschegde vleugels.
• Interactie: De 2D schokmodus ($St \approx 0.08$) blijft aanwezig en is ongevoelig voor scheg. Echter, bij hoge scheghoeken neemt de energie van de 3D modus toe tot een niveau dat vergelijkbaar is met of zelfs de 2D modus overtreft.
• Confectie: De convektiesnelheid van de buffet cells ($V_c$) volgt een schaalwet die lineair toeneemt met de tangens van de scheghoek ($V_c \propto \tan(\lambda)$), in overeenstemming met eerdere URANS-studies en theoretische modellen. De golflengte van de 3D modus blijft echter constant en wordt bepaald door de aspectverhouding van het domein, niet door de scheg.

Betekenis en Conclusie

De studie concludeert dat transsonische buffet op oneindige vleugels het resultaat is van een superpositie van twee co-existentie-instabiliteiten:

1. Een 2D schokbeweging (lage frequentie, ongevoelig voor scheg).
2. Een 3D separatie-gedreven instabiliteit (intermediaire frequentie, sterk afhankelijk van scheg).

De noodzakelijke voorwaarde voor het ontstaan van dominante 3D buffet dynamiek is een significante gemiddelde stromingsseparatie op de locatie van de schokgolf. Zonder deze separatie (zoals bij $\alpha = 5^\circ$) blijft de stroming quasi-2D, zelfs bij hoge aspectverhoudingen.

De bevindingen verklaren waarom de 2D schokmodus vaak moeilijk waar te nemen is op echte, eindige vleugels van vliegtuigen: de geometrische complexiteit en de inherente scheg van moderne vleugels bevorderen de 3D buffet cell instabiliteit, die bij hoge scheghoeken de energie van de 2D modus overneemt. Dit onderstreept het belang van het modelleren van 3D effecten en separatie voor het voorspellen en mitigeren van transsonische buffet in de luchtvaartindustrie.