A Shakhov-based Bhatnagar-Gross-Krook model for polyatomic molecules and for atomic as well as polyatomic mixtures

Diese Arbeit erweitert das im Open-Source-Code PICLas implementierte Shakhov-Bhatnagar-Gross-Krook-Modell (SBGK) auf polyatomare Moleküle und deren Gemische unter Berücksichtigung von Nichtgleichgewichtszuständen der inneren Freiheitsgrade, wobei die Validierung durch Vergleich mit DSMC- und ESBGK-Simulationen eine hohe Genauigkeit, insbesondere bei der Erfassung von Stoßwellen, bestätigt.

Das Wesentliche

Stell dir vor, du versuchst, das Wetter vorherzusagen, aber nicht nur für einen ganzen Kontinent, sondern für winzige, unsichtbare Wolken aus Gas, die sich in der extremen Kälte des Weltraums oder in winzigen Mikrochips bewegen. Das ist die Aufgabe, die sich die Forscher von der Universität Stuttgart gestellt haben.

Hier ist eine einfache Erklärung ihrer Arbeit, ohne komplizierte Formeln:

1. Das Problem: Der "Verkehrsstau" der Computer

In der Welt der Gasströmungen gibt es zwei extreme Welten:

• Die dicke Welt (Kontinuum): Wie Luft, die wir atmen. Hier fliegen die Teilchen so dicht, dass sie sich wie eine Flüssigkeit verhalten. Dafür gibt es super schnelle Computerprogramme (CFD).
• Die leere Welt (Freie Moleküle): Wie im tiefen Weltraum. Hier fliegen die Teilchen so weit voneinander entfernt, dass sie sich wie einzelne Billardkugeln verhalten, die sich selten treffen. Dafür gibt es eine sehr genaue, aber extrem langsame Methode namens DSMC (Direct Simulation Monte Carlo).

Das Problem: Wenn du etwas simulieren willst, das von der dichten zur leeren Welt übergeht (z. B. ein Raumschiff, das in die Atmosphäre eintritt), musst du beide Methoden kombinieren. Aber das ist wie der Versuch, einen Elefanten und eine Ameise Hand in Hand tanzen zu lassen – die Methoden passen nicht zusammen und sind extrem rechenintensiv.

2. Die Lösung: Ein neuer "Tanzlehrer" (Der SBGK-Modell)

Die Forscher haben eine neue Methode entwickelt, die wie ein effizienter Tanzlehrer für Gaspartikel funktioniert. Sie nennen es Shakhov-BGK.

Stell dir vor, du hast eine riesige Menge an Partikeln (Atome und Moleküle), die durcheinander tanzen.

• Die alten Methoden (wie DSMC): Jeder Tänzer muss mit jedem anderen Tänzer kollidieren, um zu lernen, wie er sich bewegen soll. Das ist wie ein riesiges, chaotisches Fest, bei dem jeder mit jedem sprechen muss. Das dauert ewig.
• Die neue Methode (SBGK): Statt auf jede einzelne Kollision zu warten, sagt der "Tanzlehrer" (das Modell) einfach: "Hey, ihr alle, bewegt euch jetzt so, als hättet ihr gerade mit allen anderen gesprochen." Es ist eine Art Gedächtnis-Trick. Das Modell berechnet nicht jeden einzelnen Stoß, sondern schätzt den Durchschnitt sehr clever. Das ist viel schneller, fast so schnell wie die schnellen CFD-Methoden, aber fast so genau wie die langsame DSMC-Methode.

3. Die Besonderheit: Moleküle sind wie "Rucksäcke"

Das Besondere an dieser neuen Version ist, dass sie nicht nur einfache Atome (wie Helium) versteht, sondern auch komplexe Moleküle (wie Stickstoff oder Kohlendioxid).

• Atome sind wie einfache Kugeln. Sie haben nur eine Art von Energie: Sie fliegen einfach herum (Translation).
• Moleküle sind wie Kugeln mit einem Rucksack. Dieser Rucksack kann sich drehen (Rotation) und vibrieren (Vibration), wenn sie warm werden.

Frühere Modelle haben diesen Rucksack oft ignoriert oder falsch behandelt. Das neue Modell von Pfeiffer und Tuttas achtet genau darauf, wie dieser Rucksack Energie aufnimmt und abgibt. Es stellt sicher, dass die "Temperatur" im Rucksack (die innere Energie) und die "Temperatur" der Bewegung (die Flugtemperatur) korrekt berechnet werden, auch wenn das Gas nicht im Gleichgewicht ist (z. B. in einer Schockwelle).

4. Der Test: Der "Blitzschnelle Kegel"

Um zu beweisen, dass ihr neuer Tanzlehrer wirklich gut ist, haben die Forscher zwei Tests gemacht:

1. Der Supersonic-Couette-Flow: Stell dir zwei riesige Platten vor, die sich sehr schnell aneinander vorbeibewegen, mit Gas dazwischen. Sie haben gemischte Gase getestet (z. B. Helium und Argon). Das Ergebnis: Der neue Tanzlehrer hat die Temperaturverteilung fast perfekt vorhergesagt, genau wie die langsame, aber genaue DSMC-Methode.
2. Der 70°-Kegel (Hyperschall): Das ist wie ein Raumschiff, das mit extrem hoher Geschwindigkeit in die Atmosphäre fliegt. Vor dem Schiff entsteht eine Schockwelle – eine unsichtbare Mauer aus extrem heißer, komprimierter Luft.
   • Hier zeigte sich der große Vorteil: Das alte Modell (ESBGK) hat die Schockwelle etwas "verschwommen" dargestellt, als würde es den Moment des Aufpralls etwas zu früh oder zu spät erkennen.
   • Das neue SBGK-Modell hat die Schockwelle scharf und präzise eingefangen. Es hat genau gesehen, wann die Luft anfängt, sich zu erhitzen.

5. Warum ist das wichtig?

Stell dir vor, du willst ein neues Raumschiff für den Mars-Flug bauen. Der Mars hat eine Atmosphäre aus Kohlendioxid (ein komplexes Molekül mit Rucksack). Wenn du die Hitze beim Eintritt nicht genau berechnest, schmilzt das Schiff.

Mit diesem neuen Modell können Ingenieure:

• Schneller rechnen: Statt Tage oder Wochen dauert die Simulation nur Stunden.
• Genauer sein: Sie sehen die Schockwellen und Hitzezonen schärfer als mit den bisherigen schnellen Methoden.
• Alles simulieren: Sie können Mischungen aus verschiedenen Gasen (Atome und Moleküle) in einem einzigen Programm behandeln.

Zusammenfassend:
Die Forscher haben einen neuen, super-schnellen Algorithmus entwickelt, der wie ein genialer Dirigent wirkt. Er sorgt dafür, dass komplexe Gasgemische im Weltraum und in der Atmosphäre so realistisch simuliert werden können, als würde man jeden einzelnen Teilchen beobachten – aber ohne die immense Rechenzeit. Das ist ein großer Schritt für die Raumfahrt und die Entwicklung neuer Technologien.

Technische Zusammenfassung

Titel: Shakhov-BGK für polyatomare Gemische

Autoren: M. Pfeiffer und F. Tuttas (Institut für Raumfahrtsysteme, Universität Stuttgart)

---

1. Problemstellung

Numerische Simulationen von Strömungen im Weltraum, im Mikrobereich sowie in Vakuumtechnologien stellen erhebliche Herausforderungen dar, da sie große Dichtegradienten über verschiedene Strömungsregime hinweg (von der Kontinuumsströmung bis zum freien Molekülfluss) abdecken.

• DSMC-Limitierung: Die etablierte Direct Simulation Monte Carlo (DSMC)-Methode liefert zwar hochgenaue Lösungen für nicht-gleichgewichtige Strömungen, ist jedoch in den Übergangs- und Kontinuumsbereichen rechenintensiv.
• CFD-DSMC-Kopplung: Die Kopplung von DSMC mit klassischen CFD-Methoden ist aufgrund unterschiedlicher mathematischer Ansätze und des statistischen Rauschens der DSMC-Methode an den Schnittstellen schwierig.
• BGK-Alternativen: Partikelbasierte Kontinuumsmethoden wie das Bhatnagar-Gross-Krook (BGK)-Modell bieten eine effiziente Alternative, da sie den Boltzmann-Kollisionsintegral durch einen Relaxationsprozess approximieren. Dies erlaubt größere Zeitschritte und weniger Partikel als DSMC.
• Spezifisches Defizit: Bestehende BGK-Modelle (wie das Standard-BGK oder ESBGK) haben oft Schwierigkeiten, den korrekten Prandtl-Zahl-Wert zu reproduzieren oder die Stoßwellenstruktur (insbesondere bei hypersonischen Strömungen) präzise abzubilden. Zudem fehlen oft robuste Erweiterungen für komplexe Gemische aus Atomen und polyatomaren Molekülen unter Berücksichtigung nicht-gleichgewichtiger innerer Freiheitsgrade (Rotation, Vibration).

2. Methodik

Die Autoren erweitern die Implementierung des Shakhov-BGK (SBGK)-Modells im Open-Source-Partikelcode PICLas.

• Theoretisches Fundament:

  • Das BGK-Operator approximiert die Boltzmann-Gleichung durch eine Relaxation der Verteilungsfunktion $f_s$ gegen eine Zielverteilung $f_s^t$.
  • Polyatomare Moleküle: Die Zielverteilung wird faktorisiert in translatorische, rotatorische und vibronische Anteile. Die inneren Freiheitsgrade (DOF) werden quantisiert behandelt (harmonische Oszillatoren für Vibration).
  • Shakhov-Korrektur: Im Gegensatz zum Standard-BGK wird die Zielverteilung so modifiziert, dass der Wärmestromvektor korrekt abgebildet wird. Dies ermöglicht die Einstellung der korrekten Prandtl-Zahl ($Pr$) des Gemisches.
  • Einzelner Relaxationsterm: Das Modell verwendet für jede Spezies nur einen Relaxationsterm. Dies vereinfacht die Berechnung im Vergleich zu Multi-Relaxations-Modellen erheblich, da nur die Momente des Gesamtgemisches (nicht jeder einzelnen Spezies) für den Relaxationsterm benötigt werden.

• Bestimmung der Transportgrößen:

  • Zur Berechnung von Viskosität und Wärmeleitfähigkeit für Gemische wird die erste Näherung basierend auf Kollisionsintegralen (Collision Integrals) verwendet, anstatt einfacher Mischungsregeln (wie Wilke's Regel), um höhere Genauigkeit zu gewährleisten.
  • Die VHS-Parameter (Variable Hard Sphere) für Kollisionen werden entweder direkt oder durch Mittelung für verschiedene Speziespaare bestimmt.

• Erhaltungseigenschaften:

  • Das Modell erfüllt die Erhaltungssätze für Masse, Impuls und Energie.
  • Die Relaxationsfrequenz und die effektiven Relaxationstemperaturen werden so gewählt, dass die Energieübertragung den Landau-Teller-Formen folgt.

3. Wichtige Beiträge

1. Erweiterung auf komplexe Gemische: Das SBGK-Modell wurde erstmals erfolgreich auf Mischungen aus Atomen und polyatomaren Molekülen (mit quantisierten Vibrationszuständen) angewendet.
2. Prandtl-Zahl-Modellierung: Herleitung des Modellparameters für die Prandtl-Zahl in Gemischen, der von den Massenanteilen und den inneren Freiheitsgraden abhängt.
3. Effizienzsteigerung: Durch den Einsatz eines einzelnen Relaxationsterms wird der Rechenaufwand im Vergleich zu Multi-Relaxations-Modellen signifikant reduziert, ohne die physikalische Genauigkeit für makroskopische Größen zu opfern.
4. Implementierung: Vollständige Integration in den Open-Source-Code PICLas mit Akzeptanz-Ablehnungs-Sampling (Acceptance-Rejection Sampling) für die neuen translatorischen Zustände.

4. Ergebnisse und Validierung

Das Modell wurde mit zwei Haupttestfällen validiert und mit DSMC sowie dem ESBGK (Ellipsoidal Statistical BGK)-Modell verglichen:

• Fall 1: Supersonische Couette-Strömung

  • Untersucht wurden reine $N_2$-Strömungen, ein Ar-He-Gemisch und ein $N_2$-$N$-Gemisch.
  • Ergebnis: Die SBGK-Ergebnisse zeigen eine sehr gute Übereinstimmung mit DSMC für Temperaturprofile.
  • Beobachtung: Bei großen Massenunterschieden (Ar-He) traten leichte Abweichungen auf, die auf die Approximation der effektiven Prandtl-Zahl zurückgeführt werden. Bei ähnlicheren Massenverhältnissen ($N_2$-$N$) war die Übereinstimmung nahezu perfekt.

• Fall 2: Hypersonische Strömung um einen 70°-abgerundeten Kegel

  • Getestet wurden: Reines $N_2$, ein $CO_2$-$N_2$-Gemisch (polyatomar-diatomar) und ein Dreikomponenten-Gemisch ($N_2$-$O_2$-$N$).
  • Schlüsselbefund (Stoßwelle): Das SBGK-Modell erfasst den Beginn der Stoßwelle und die Stoßstruktur (insbesondere den Anstieg der translatorischen Temperatur) präziser als das ESBGK-Modell. ESBGK neigt dazu, den Temperaturanstieg zu früh einzuleiten, was zu einer breiteren, weniger scharfen Stoßwelle führt.
  • Wärmeübergang: Sowohl SBGK als auch ESBGK zeigen hervorragende Übereinstimmung mit DSMC und experimentellen Daten für den Wärmestrom an der Kegelfläche im Nachlaufbereich.
  • Nicht-Gleichgewicht: Das Modell bildet die Nicht-Gleichgewichts-Zustände der Rotations- und Vibrationsenergien korrekt ab.

5. Bedeutung und Ausblick

• Genauigkeit vs. Effizienz: Das vorgestellte SBGK-Modell bietet einen optimalen Kompromiss zwischen der hohen Genauigkeit von DSMC (insbesondere in Stoßwellenregionen) und der Recheneffizienz von Kontinuumsmethoden.
• Anwendbarkeit: Es ist besonders geeignet für die Simulation von Wiedereintrittsströmungen (z. B. Mars-Eintritt mit $CO_2$-$N_2$-Gemischen) und anderen hochenergetischen Strömungen, wo innere Freiheitsgrade und Nicht-Gleichgewichtseffekte entscheidend sind.
• Zukunft: Als nächster Schritt ist die Implementierung chemischer Reaktionen in das Modell geplant, um reaktive Strömungen vollständig abbilden zu können.

Fazit: Die Arbeit demonstriert, dass das erweiterte Shakhov-BGK-Modell eine leistungsfähige, effiziente und physikalisch konsistente Methode zur Simulation von polyatomaren Gasgemischen darstellt, die insbesondere in der Erfassung von Stoßwellenstrukturen dem etablierten ESBGK-Modell überlegen ist.