Vergleichende Analyse kompressibler Löser-Schemata für die Aerodynamik unterexpandierter Düsen mit Schlieren-Validierung

Warum schnelle Gasstrahlen wichtig sind

Von Raketendüsen über Industrieabführungen bis hin zu Vulkanausbrüchen beeinflussen schnelle Gasstrahlen, wie Schub, Lärm und Wärme in die umgebende Luft übertragen werden. Ihr Verhalten genau vorherzusagen ist anspruchsvoll, weil sich Wellen aus komprimiertem und expandiertem Gas zu komplexen Mustern formen, die herkömmliche Rechenwerkzeuge schwer nachverfolgen können. Diese Studie stellt eine praktische Frage: Kann ein preiswerterer Typ von numerischem Löser, der üblicherweise für langsamere Strömungen eingesetzt wird, diese anspruchsvollen Strahlprobleme ohne Genauigkeitsverlust bewältigen?

Zwei verschiedene Wege, dieselbe Strömung zu lösen

Ingenieure nutzen häufig Computational Fluid Dynamics (CFD), um Entwürfe vor dem Bau zu prüfen. Für schnelle Gasströmungen galt lange ein „dichtebasierter“ Löser als Standard, weil er auf schnelle Druck- und Dichteänderungen in der Nähe von Schockwellen abgestimmt ist. Ein „druckbasierter“ Löser dagegen wird meist für langsame, nahezu inkompressible Strömungen verwendet und ist dafür wirtschaftlicher. In den letzten Jahren berichteten einige Forscher, dass druckbasierte Löser auch bei höheren Geschwindigkeiten gut funktionieren können; ihre Leistungsfähigkeit bei stark unterexpandierten Strahlen – wenn Gas die Düse mit deutlich höherem Druck als die Umgebung verlässt – war jedoch bislang unklar.

Eine Laborstrahl aufbauen und unsichtbare Wellen sichtbar machen

Um beide Ansätze fair zu prüfen, bauten die Autoren ein einfaches, aber sorgfältig kontrolliertes Experiment auf. Druckluft wurde aus einem zylindrischen Vorratstank durch eine kleine konvergente Düse in ruhende Raumluft freigesetzt. Durch Variation des Versorgungsdrucks stellten sie vier Druckverhältnisse zwischen Tank und Atmosphäre ein, beginnend knapp über dem „geblockten“ Zustand, bei dem das Gas die Düse mit Schallgeschwindigkeit verlässt, bis hin zu deutlich unterexpandierten Strahlen mit ausgeprägten Schockmustern. Um diese sonst unsichtbaren Strukturen zu sehen, nutzten sie eine Schlieren-Optik: einen hellen Punktleuchter, einen konkaven Spiegel, eine Rasierklinge und eine Digitalkamera, so angeordnet, dass kleine Dichteänderungen als helle und dunkle Bänder sichtbar werden. Diese Aufnahmen lieferten eine visuelle Karte der stromabwärts der Düse entstehenden Schockzellen.

Die Computermodelle auf die Probe stellen

Parallel dazu wurden Düse und Strömungsbedingungen im CFD-Code ANSYS Fluent rekonstruiert. Die Forschenden verwendeten für beide Löser identische Netze, Randbedingungen und Turbulenzmodelle und veränderten nur das zugrundeliegende Lösungsverfahren. Sie verfolgten Druck, Geschwindigkeit und Turbulenz entlang und um den Strahl und überprüften, wie schnell jeder Löser zu einer stabilen Lösung konvergierte. Beide Werkzeuge reproduzierten Schlüsselmerkmale: die Bildung eines potenziellen Kerns im geblockten Zustand, das Auftreten diamantförmiger Schockzellen, sobald der Strahl unterexpandierte, und die zunehmende Länge und Stärke dieser Zellen mit steigendem Druckverhältnis. Der erste starke Schock, als Mach-Scheibe bekannt, erschien an nahezu denselben Positionen in beiden Simulationen und in veröffentlichten experimentellen Daten, mit nur wenigen Prozent Unterschied in der Spitzengeschwindigkeit.

Wo die Löser übereinstimmen und wo sie sich unterscheiden

Entlang der Strahlmitte lieferten die beiden Löser nahezu identische Oszillationen in Druck und Mach-Zahl, die das wiederkehrende Muster von Kompression und Expansion widerspiegeln. Unterschiede zeigten sich weiter stromabwärts, wo der Strahl auf Unterschallgeschwindigkeiten abgesenkt war und Turbulenz dominierte. Dort tendierte der dichtebasierte Löser dazu, höhere turbulente kinetische Energie vorherzusagen als der druckbasierte Löser, was mit bekannten Einschränkungen dichtebasierter Verfahren bei niedrigen Geschwindigkeiten übereinstimmt. Trotz dessen lieferten beide Löser nahezu denselben Massenstrombeiwert (discharge coefficient), ein Maß dafür, wie effektiv die Düse den Massenstrom passieren lässt, mit Unterschieden von weniger als 0,2 Prozent. In Bezug auf den Rechenaufwand erreichte der druckbasierte Löser das gleiche strenge Konvergenzziel in etwa 22 Prozent kürzerer CPU-Zeit.

Was das für reale Konstruktionen bedeutet

Für Ingenieure, die Hochgeschwindigkeitsstrahlen aus relativ einfachen Düsen simulieren müssen, liefert diese Studie beruhigende Erkenntnisse. Ein druckbasierter Löser, obwohl ursprünglich für langsamere, inkompressible Strömungen entwickelt, kann die wesentlichen Schockmuster und Leistungskennwerte unterexpandierter Überschallstrahlen mit vergleichbarer Genauigkeit wie ein spezialisierterer dichtebasierter Löser erfassen und dabei weniger Rechenzeit benötigen. Die Autoren weisen darauf hin, dass dichtebasierte Verfahren in manchen Bereichen nach wie vor Vorteile haben, insbesondere dort, wo sehr starke Schocks und komplexe Chemie beteiligt sind. Für viele praktische Düsenprobleme kann jedoch der kostengünstigere druckbasierte Ansatz eine verlässliche Wahl sein und es den Konstrukteuren ermöglichen, mehr Varianten ohne prohibitive Simulationskosten zu erkunden.

Zitation: Alsaedi, S.S., Al-Sadawi, L.A., Al-Haddad, L.A. et al. Comparative analysis of compressible solver schemes for underexpanded jet aerodynamics with Schlieren validation. Sci Rep 16, 15724 (2026). https://doi.org/10.1038/s41598-026-44651-w

Schlüsselwörter: unterexpandierter Strahl, Überschalldüse, Schlieren-Aufnahme, CFD-Löser, Schockwellen