17.04.2023

Simulierte Schockwellen

Neuer Ansatz erleichtert Berechnung von stark turbulenten Strömungen.

Es muss nicht ein Hurrikan oder ein Tsunami sein – selbst ein dünner Wasserstrahl im Waschbecken löst eine physika­lische Schockwelle aus. Nun hat der Mathematiker Siddhartha Mishra von der Technischen Hochschule Zürich auf dem Supercomputer „Piz Daint“ am Schweizerischen Nationalen Hochleistungs­rechenzentrum CSCS einen Weg gefunden, die Schwierig­keiten beim Simulieren stark turbulenter Strömungen zu überwinden.

Abb.: Simulierte Schockwelle eines Tsunamis in Alaska: Das Bild zeigt ihre...
Abb.: Simulierte Schockwelle eines Tsunamis in Alaska: Das Bild zeigt ihre Ausbreitung 39 Sekunden nach dem Aufprall des Felsmaterials. (Bild: Sánchez-​Linares et al.)

Seien es Tornados, Lawinen, Überschwemmungen oder die Luftströme, die ein Flugzeug tragen: Starke Strömungen finden sich überall um uns herum. „Sie zu verstehen ist entscheidend, um technische Prozesse zu verbessern oder Natur­phänomene genauer vorhersagen“, sagt Siddhartha Mishra. Mishra interessiert sich vor allem für mächtige Phänomene, solche mit starken Turbulenzen oder Schockwellen. Zusammen mit seinem Team untersucht er solche Strömungs­dynamiken und entwickelt Algorithmen, um sie im Computer präziser und effizienter zu simulieren. Noch mehr jedoch versucht der Mathematiker, die grundlegenden Grenzen heutiger Strömungs­simulationen zu sprengen. 

Die eigentlichen mathe­matischen Gleichungen, um Flüssigkeiten zu beschreiben sind schon fast klassischer Herkunft: Der Schweizer Mathematiker Leonhard Euler hat sie vor fast 270 Jahren formuliert. Zwar berücksichtigte Euler damals noch nicht die Viskosität einer Strömung, dennoch lassen sich seine Differential­gleichungen heute für viele Flüsse anwenden, bei denen diese Eigenschaft vernachlässigbar ist. Auf dem Papier sehen die Gleichungen erstaunlich einfach aus – und sind doch unglaublich schwer zu lösen, einfach weil die Phänomene, die sie beschreiben, so komplex sind. „Wir haben die Gleichungen auch heute noch nicht komplett verstanden“, sagt Mishra.

Dieses Unverständnis lässt sich in der Praxis umgehen: Expertinnen und Experten für Fluiddynamik nutzen Algorithmen, um auf Super­computern wie dem „Piz Daint“ in Lugano Annäherungen an die Lösungen der Gleichungen zu berechnen. Solche numerischen Näherungs­methoden sind vielfach erfolgreich: Damit lassen sich etwa die tragenden Luftströmungen um Flugzeug­flügel oder Wasserläufe durch Turbinen simulieren und dadurch besser verstehen oder Wetterphänomene genauer vorhersagen. Mishras Team hat unter anderem verbesserte Algo­rithmen zur Simulation und Vorhersage von Tsunamis und Lawinen oder zur Modellierung von Supernovas entwickelt.

Doch streng genommen sind diese Methoden nicht mathematisch präzis. Für einige Phänomene sind sie an ihre Grenzen gestoßen. „Normaler­weise würde man erwarten, dass eine höhere Auflösung genauere Ergebnisse liefert“, sagt Mishra. Anders gesagt: Bezieht man in eine Berechnung mehr einzelne Punkte aus Zeit und Raum ein, sollte der Fehler eigentlich kleiner werden und die berechneten Näherungswerte sollten die Realität besser abbilden. Doch bei stark turbulenten Flüssig­keiten funktioniert das nicht, wie Mishra zeigte. Im Gegenteil: Bei hoher Auflösung erzielte Ergebnisse stimmen dann überhaupt nicht mit den Ergeb­nissen von schwächer aufgelösten Berechnungen überein – sie konvergieren nicht, sondern sehen komplett anders aus. „Das bedeutet auch, dass sich für solche Phänomene keine Vorhersagen berechnen lassen.“

Darum haben Mishra und sein Team nach einer Möglichkeit gesucht, diese Schwierig­keiten beim Simulieren stark turbulenter Ströme zu überwinden, und zwar mit statistischen Lösungen. Dazu haben die Forschenden das Problem randomisiert, also den Zufall berücksichtigt: Sie erzeugten in den untersuchten Strömungen viele winzige, zufällige Störungen – und analysierten dann das gemittelte Ergebnis. „Das ist die Grundlage für statis­tische Lösungen“, erklärt Mishra. „Wenn sich einzelne Messungen oder Experimente nicht konvergent verhalten, kann man sich stattdessen Durchschnitts­werte anschauen. Bei turbulenten Strömungen sind die Details das Problem. Deshalb sieht man mit Durchschnitts­werten mehr Struktur.“ Zusätzlich müsse man auch berück­sichtigen, dass die statistischen Eigen­schaften von Strömungen an verschiedenen Orten im Raum voneinander abhängig seien.

„Wir können solche statistischen Lösungen berechnen“, sagt Mishra. Tatsächlich konver­gierten die berechneten Mittelwerte nun mit höherer Auflösung, wie das Team erkannte. Das galt nicht nur für spezifische Größen, etwa die Dichte oder die Geschwin­digkeit einer Strömung, sondern auch für die statistischen Abweichungen dieser Größen und ihre Wahrscheinlichkeits­verteilung. „Bislang hat jede unserer Test­simulationen mit vereinfachten, zwei­dimensionalen Strömungen funktioniert.“ 

Jüngst haben die Forschenden ihren Code optimiert, sodass sie die Simulationen auf Piz Daint um mehr als das Zehnfache beschleunigen konnten. In den Simulationen wurden auf die Strömungen zunächst virtuell bestimmte Scherkräfte ausgeübt, um die Entstehung von Kelvin-Helmholz-Instabilitäten zu simulieren. Diese führen zu spezifischen Wirbeln, wie sie beispiels­weise im kräuselnden Rauch einer Kerze sichtbar sind, ebenso in bestimmten Wolkenformationen oder in Gasen in der Atmosphäre von Planeten. Auch hier zeigten die Ergebnisse: Betrachten die Forschenden einzelne Simulationen, ist keine Konvergenz zu erkennen. „Arbeiten wir aber mit den statistischen Eigen­schaften zahlreicher Simulationen, wie Mittelwerte, Abweichungen und Wahrscheinlichkeiten, konvergieren die Lösungen“, bestätigt Mishra.

Kürzlich führte sein Team ähnliche Simulationen eines anderen Phänomens durch, die Richtmeyer-Meshkov-Instabilität. Diese kennt man vor allem in der Astrophysik, sie entsteht, wenn eine Schockwelle auf die Grenzfläche zwischen zwei verschiedenen Strömungen trifft. „Eine solche Instabilität ist heftig, ähnlich wie eine Explosion“, erklärt Mishra. Wiederum zeigten die Simulationen und Analysen des Teams, dass die statis­tischen Lösungen vieler Simulationen mit zunehmender Auflösung konvergieren, wohingegen man mit einzelnen Simulationen ohne Konvergenz nicht weit kommt. „Durch die Computer­simulationen verstehen wir nun viel besser, was vor sich geht“, fasst Mishra zusammen. 

ETHZ / JOL

Weitere Infos

Sonderhefte

Physics' Best und Best of
Sonderausgaben

Physics' Best und Best of

Die Sonder­ausgaben präsentieren kompakt und übersichtlich neue Produkt­informationen und ihre Anwendungen und bieten für Nutzer wie Unternehmen ein zusätzliches Forum.

Virtuelle Jobbörse

Virtuelle Jobbörse
Eine Kooperation von Wiley-VCH und der DPG

Virtuelle Jobbörse

Innovative Unternehmen präsentieren hier Karriere- und Beschäftigungsmöglichkeiten in ihren Berufsfeldern.

Die Teilnahme ist kostenfrei – erforderlich ist lediglich eine kurze Vorab-Registrierung.

Meist gelesen

Themen