Leveraging unstructured grids for direct numerical simulations of wall turbulence

Dieser Beitrag stellt das {\eta}-Gitter vor, ein unstrukturiertes Gittergenerierungsframework für direkte numerische Simulationen von Wandturbulenz, das die Gittergrößen mit der lokalen Kolmogorov-Skala skaliert und dabei eine Genauigkeit erreicht, die mit konventionellen kartesischen Gittern vergleichbar ist, während es die Rechenkosten erheblich reduziert, insbesondere bei hohen Reynolds-Zahlen und über komplexen Rippelgeometrien.

Das Wesentliche

Stellen Sie sich vor, Sie versuchen zu simulieren, wie Luft über ein Auto strömt oder Wasser an einem Schiffsrumpf vorbeifließt. Um dies auf einem Computer genau zu berechnen, verwenden Wissenschaftler eine Technik namens Direkte Numerische Simulation (DNS). Denken Sie an DNS als die Erstellung eines riesigen, dreidimensionalen digitalen Mikroskops, das die Flüssigkeit (Luft oder Wasser) in Millionen winziger, unsichtbarer Würfel (ein Gitter) zerlegt. Der Computer berechnet dann, wie sich jeder einzelne Würfel bewegt und mit seinen Nachbarn interagiert.

Das Problem ist, dass die Strömung nahe einer Oberfläche (wie der Seite eines Schiffes) unglaublich chaotisch und detailliert ist. Um ein klares Bild zu erhalten, benötigen Sie eine enorme Anzahl dieser winzigen Würfel direkt neben der Oberfläche. Wenn Sie sich jedoch weiter von der Oberfläche entfernen, glättet sich das Chaos, und Sie benötigen keine so winzigen Würfel mehr.

Der alte Weg: Die „steife Mauer aus Ziegelsteinen"

Traditionell verwendeten Wissenschaftler ein kartesisches Gitter. Stellen Sie sich vor, Sie bauen eine Mauer aus identischen, steifen Ziegelsteinen.

• Das Problem: Um die winzigen Details nahe der Oberfläche zu erkennen, müssen Sie die Ziegelsteine am Boden sehr klein machen. Da diese Ziegelsteine jedoch steif sind und in einer geraden Linie verbunden sind, sind Sie gezwungen, diese gleichen winzigen Ziegelsteine bis zur Spitze Ihrer Mauer zu verwenden, auch dort, wo die Details nicht wichtig sind.
• Das Ergebnis: Sie landen bei einer Mauer aus Milliarden winziger Ziegelsteine, von denen die meisten unnötig sind. Dies macht die Computersimulation unglaublich langsam und teuer, wie der Versuch, jeden einzelnen Sandkorn an einem Strand zu zählen, nur um die Gezeiten zu messen.

Die neue Lösung: Das „intelligente, dehnbare Netz"

Dieser Artikel stellt eine neue Methode vor, die $\eta$-Gitter (Eta-Gitter) genannt wird. Anstelle von steifen Ziegelsteinen stellen Sie sich ein intelligentes, dehnbare Fischernetz vor.

• Wie es funktioniert: Die Autoren haben ein System entwickelt, bei dem die Größe der Maschen des Netzes automatisch basierend auf dem benötigten Detailgrad variiert.
  • Nahe der Oberfläche (Die „innere Schicht"): Das Netz hat sehr kleine, enge Maschen, um die winzigen, chaotischen Wirbel der Flüssigkeit einzufangen.
  • Weiter entfernt (Die „äußere Schicht"): Wenn die Flüssigkeit ruhiger wird, dehnt sich das Netz automatisch aus und macht die Maschen viel größer.
• Der geheime Bestandteil: Die Größe dieser Maschen basiert auf etwas, das als Kolmogorov-Skala (bezeichnet als $\eta$) bekannt ist. Denken Sie daran als die „kleinstmögliche Wirbel", die in der Flüssigkeit in einer bestimmten Höhe existieren kann. Das neue Gitter sagt einfach: „Machen Sie die Maschengröße gerade groß genug, um den kleinsten Wirbel in dieser spezifischen Höhe einzufangen, und nicht größer."

Warum dies eine große Sache ist

Die Autoren testeten dieses „intelligente Netz" auf zwei verschiedenen Arten von Computercodes (einer ähnlich der Spektralelementmethode, der andere der Finite-Volumen-Methode) und verglichen es mit der alten Methode der „steifen Ziegelsteine".

1. Genauigkeit: Die Ergebnisse waren nahezu identisch. Das „intelligente Netz" erfasste die Physik genauso gut wie die „steifen Ziegelsteine", mit weniger als 1 % Unterschied bei Schlüsselmessungen wie Reibung und Geschwindigkeit.
2. Massive Einsparungen: Hier passiert die Magie.
   • Bei glatten Oberflächen (wie einer flachen Wand) reduzierte das neue Gitter die Anzahl der erforderlichen „Ziegelsteine" (Gitterpunkte) bei hohen Geschwindigkeiten um etwa 90 %.
   • Bei rauen Oberflächen (wie einer Wand mit winzigen Rillen, sogenannten „Riblets", die entwickelt wurden, um den Widerstand zu verringern), waren die Einsparungen noch dramatischer – bis zu 97 % weniger Gitterpunkte.

Die Analogie der „gerillten Wand"

Um den Teil mit den Riblets zu verstehen, stellen Sie sich eine Wand vor, die mit winzigen, parallelen Rillen bedeckt ist (wie die Textur eines Golfballs oder der Haut eines Haies).

• Der alte Weg: Um dies zu simulieren, musste die Methode der steifen Ziegelsteine die Ziegelsteine überall winzig halten, da die Rillen das Gitter zwangen, bis weit nach oben fein zu sein. Es war wie der Versuch, jeden einzelnen Faden in einem Pullover zu zählen, auch die Teile weit entfernt vom Stoff.
• Der neue Weg: Das „intelligente Netz" weiß, dass die Strömung einige Zentimeter über diesen winzigen Rillen wieder glatt wird. Es dehnt die Maschen sofort über den Rillen aus und ignoriert die winzigen Details, die nicht mehr wichtig sind.

Das Fazit

Die Autoren haben einen Rahmen geschaffen, der wie eine intelligente Zoomlinse für Strömungssimulationen wirkt. Sie konzentriert ihre Rechenleistung genau dort, wo sie benötigt wird (nahe der Wand), und lockert sie dort, wo sie es nicht ist.

• Für glatte Wände: Sie skaliert den Aufwand viel langsamer, wenn die Simulation größer wird.
• Für raue Wände: Sie skaliert noch besser und macht Simulationen komplexer, widerstandsreduzierender Oberflächen auf Computern möglich, die dies zuvor nicht bewältigen konnten.

Kurz gesagt, sie haben einen Weg gefunden, die gleiche hochwertige Arbeit mit einem Bruchteil der Rechenleistung zu leisten und eine Aufgabe, die einen Supercomputer einen Monat kosten könnte, in eine zu verwandeln, die in wenigen Tagen erledigt werden kann.

Technische Zusammenfassung

Technische Zusammenfassung: Nutzung unstrukturierter Gitter für direkte numerische Simulationen von Wandturbulenz

Problemstellung
Die direkte numerische Simulation (DNS) von Wandturbulenz ist rechnerisch sehr aufwendig, hauptsächlich aufgrund der Notwendigkeit, die Kolmogorov-Längenskala ($\eta$) über das gesamte Rechengebiet aufzulösen. Herkömmliche Gitterlöser auf kartesischen Gittern, die weit verbreitet für glatte und nicht-glatte Oberflächen sind (z. B. mittels Immersed-Boundary-Methoden), erzwingen typischerweise eine feste oder langsam variierende Strömungsrichtung ($\Delta x^+$) und Querrichtung ($\Delta z^+$) des Gitters im gesamten Rechengebiet. Dieser Ansatz führt zu einer signifikanten Überauflösung in den äußeren Bereichen, wo die Kolmogorov-Skala wächst, und erzeugt so einen rechnerischen Engpass. Dies ist besonders ausgeprägt bei Simulationen mit komplexen Oberflächen wie Rillen (Riblets), wo das Gitter in der Nähe der Oberfläche extrem fein sein muss, um Mikro-Nuten aufzulösen, diese feine Auflösung jedoch in kartesischen Rahmenbedingungen unnötigerweise über die gesamte Dicke der Grenzschicht beibehalten wird.

Methodik
Die Autoren schlagen ein neues Gerüst zur Generierung unstrukturierter Gitter vor, das als $\eta$-Gitter bezeichnet wird. Das Kernprinzip besteht darin, die wandnormale ($\Delta y^+$) und die Querrichtung ($\Delta z^+$) der Gitterzellen proportional zur lokalen Kolmogorov-Skala ($\eta^+$) zu setzen, was eine Vergröberung des Gitters in den äußeren Strömungsbereichen ermöglicht, wo die Turbulenzskalen größer sind.

1. Gitterformulierung:

   • Innere Schicht: Eine dünne Schicht ($\sim 50$ viskose Einheiten) angrenzend an die Wand verwendet viskosskalierte Gittergrößen, die einer konventionellen DNS ähneln ($0.3 \lesssim \Delta y^+ \lesssim 4$, $\Delta z^+ \simeq 5$). Für nicht-glatte Oberflächen (z. B. Rillen) löst das Gitter die kleinste Oberflächenwellenlänge ($\ell^+$) mit $\ell^+/30 \lesssim \Delta y^+, \Delta z^+ \lesssim 4$ auf.
   • Äußerer Bereich: Oberhalb der inneren Schicht nehmen die Gittergrößen proportional zur lokalen Kolmogorov-Skala zu ($\Delta y^+ \simeq \Delta z^+ \simeq 2\eta^+$).
   • Skalierungsgesetze: Die Autoren verfeinern semi-empirische Anpassungen für $\eta^+$ sowohl im logarithmischen als auch im Nachlaufbereich für turbulente Kanalströmungen und turbulente Grenzschichten (TBL) ohne Druckgradienten (ZPG). Sie schlagen stückweise Potenzgesetz-Anpassungen für $\eta^+$ vor, die bei $y^+ \approx \delta^+/2$ übergehen.
   • Rillen-Anpassung: Für Rillen-Geometrien führt das Gerüst eine „Rillen-Unterschicht" ($\delta_r^+$) ein, die mit quadratischen Zellen gefüllt ist, um Sekundärströmungen und Kelvin-Helmholtz-Rollen aufzulösen, und die nahtlos in die $\eta$-basierte Skalierung darüber übergeht.

2. Implementierung:

   • Das Gerüst nutzt das Open-Source-Paket Gmsh zur Generierung von Mehrblock-unstrukturierten Gittern.
   • Die Gitter werden mit zwei unterschiedlichen Lösern getestet: SOD2D (Spektral-Elemente-Methode) und OpenFOAM (Finite-Volumen-Methode).
   • Die Simulationen umfassen turbulente Kanalströmungen und ZPG-TBLs über glatte Wände sowie verschiedene Rillen-Geometrien (dreieckige Nuten) bis zu einer Reibungs-Reynolds-Zahl $\delta_0^+ = 1000$.

Hauptbeiträge

• Gerüstentwicklung: Eine robuste Mehrblock-Strategie zur Generierung unstrukturierter Gitter, die den Gitterabstand dynamisch an die lokale Kolmogorov-Skala anpasst und sowohl für glatte als auch für komplexe raue Oberflächen anwendbar ist.
• Verfeinerte Skalierungsgesetze: Aktualisierte semi-empirische Anpassungen für die Kolmogorov-Skala im Nachlaufbereich von TBLs und Kanalströmungen, die die Genauigkeit der Gittergrößenbestimmung im äußeren Bereich verbessern.
• Löser-Unabhängigkeit: Der Nachweis, dass das $\eta$-Gitter-Gerüst konsistente Ergebnisse hoher Genauigkeit über verschiedene numerische Diskretisierungsschemata (SEM und FVM) hinweg liefert.
• Rillen-spezifische Formulierung: Eine maßgeschneiderte Gitterstrategie für Rillen, die die Physik der Unterschicht (Sekundärströmungen, KH-Rollen) auflöst und gleichzeitig eine schnelle Vergröberung im äußeren Strömungsbereich ermöglicht, wodurch die Einschränkungen kartesischer/IBM-Ansätze überwunden werden.

Ergebnisse

• Genauigkeit: Das $\eta$-Gitter liefert Ergebnisse mit weniger als 1 % Abweichung im Vergleich zu konventionellen kartesischen Gittern und Referenz-DNS-/Experimentaldaten. Zu den Metriken gehören Hautreibungskoeffizienten ($C_f$), mittlere Geschwindigkeitsprofile ($U^+$), turbulente Spannungen und Spektrogramme.
• Gittereffizienz:
  • Glatte Wände: Die Anzahl der Gitterpunkte ($N_\eta$) skaliert mit $\propto \delta_0^{+2.5}$, im Vergleich zu $\propto \delta_0^{+3.0}$ für kartesische Gitter mit hyperbolischer Tangenten-Streckung. Bei $\delta_0^+ = 6000$ benötigt das $\eta$-Gitter nur $\sim 10\%$ der Gitterpunkte eines Standard-kartesischen Gitters ($N_\eta/N_{Tanh} \simeq 0.1$).
  • Rillen: Der Effizienzgewinn ist noch ausgeprägter. Für drag-reduzierende dreieckige Rillen skaliert $N_\eta$ mit $\propto \delta_0^{+2.0}$, während kartesische Gitter mit $\propto \delta_0^{+3.0}$ skalieren. Bei $\delta_0^+ = 6000$ benötigt das $\eta$-Gitter nur $\sim 3\%$ der Gitterpunkte eines kartesischen Gitters ($N_\eta/N_{Tanh} \simeq 0.03$).
• Rillen-Physik: Das Gerüst erfasst erfolgreich komplexe Strömungsmechanismen über Rillen, einschließlich des virtuellen Ursprungs, Sekundärströmungen und Kelvin-Helmholtz-Rollen, wobei die Gitterkonvergenz mit moderaten Gitterparametern erreicht wird ($C_y=2.0, C_z=2.5, y_{in}^+=50$).
• TBL-Validierung: DNS von ZPG-TBLs über Rillen zeigte eine gute Übereinstimmung mit Experimentaldaten (Baron & Quadrio; Choi & Orchard) hinsichtlich des prozentualen Widerstandsabbaus und der Geschwindigkeitsprofile, trotz inhärenter Diskrepanzen zwischen Experiment und DNS in den Nachlaufbereichen.

Bedeutung
Die Arbeit behauptet, dass das $\eta$-Gitter-Gerüst eine „enorme Gittereinsparung" bietet und damit DNS von Wandturbulenz über komplexen Oberflächen bei hohen Reynolds-Zahlen rechnerisch machbar macht. Durch die Ausrichtung der Gitterauflösung auf die lokalen physikalischen Skalen (Kolmogorov-Skala) anstatt der Aufrechterhaltung fester Auflösungen reduziert die Methode die Rechenkosten (Gitterpunkte und Zeitschrittbegrenzungen) drastisch, ohne die Genauigkeit zu beeinträchtigen. Die Autoren betonen, dass dieser Ansatz die Durchführbarkeit von DNS von TBLs über komplexen Oberflächen bei Reynolds-Zahlen ermöglicht, die experimentellen Bedingungen entsprechen, eine Aufgabe, die zuvor durch die prohibitiven Kosten kartesischer Methoden behindert wurde. Die Arbeit validiert die Wirksamkeit unstrukturierter Gitter für hochpräzise Turbulenzsimulationen und geht damit über die Beschränkungen strukturierter kartesischer Gitter hinaus.