No evidence of vorticity production from initially irrotational turbulent gravitational collapse

Technische Zusammenfassung: Kein Nachweis von Vortizitätsproduktion aus anfänglich irrotationaler turbulenter Gravitationskollaps

Problemstellung
Es wird weithin angenommen, dass Gravitationskollaps Turbulenz erzeugt, welche potenziell ein kleinskaliges Dynamo-Prinzip speisen könnte, um Magnetfelder unterhalb der Jeans-Skala zu verstärken. Kleinskalige Dynamos erfordern jedoch grundlegend vortikale (rotierende) Turbulenzen, während Gravitationskollaps primär irrotationale (potenzielle) Strömungen erzeugt. In Abwesenheit von Magnetfeldern oder externen Antriebskräften und unter Verwendung einer barotropen Zustandsgleichung (bei der Druck- und Dichtegradienten parallel verlaufen), kann der barokline Term keine Vortizität erzeugen. Folglich ist der einzige verfügbare Mechanismus zur Erzeugung von Vortizität aus einer anfänglich irrotationalen Strömung die Wechselwirkung mit der Viskosität. Frühere Studien haben nahegelegt, dass kollapsinduzierte Kompression einen Dynamo-Effekt antreiben könnte, es bleibt jedoch unklar, ob der Kollapsfluss selbst die notwendige Vortizität erzeugt oder ob beobachtete Dynamo-Effekte lediglich das Abklingen bereits existierender initialer Turbulenz sind.

Methodik
Die Autoren führen direkte numerische Simulationen (DNS) eines Gravitationskollapses unter Verwendung des Pencil Code durch, wobei sie sechswertige zentrierte Differenzen und ein Schema mit dritter Ordnung für die Zeitschritte anwenden. Die Studie nutzt ein periodisches kubisches Gebiet mit einer isentropen Zustandsgleichung ($\gamma = 5/3$) und einem modifizierten isothermen Bonnor-Ebert-Sphäroid als initiales Dichteprofil.

Wesentliche methodische Entscheidungen umfassen:

• Anfangsbedingungen: Das Geschwindigkeitsfeld wird im Fourier-Raum strikt irrotational konstruiert, mit zufälligen Phasen und einem spezifischen Spektralexponent. Es werden keine Magnetfelder einbezogen, um die Mechanismen der Vortizitätserzeugung zu isolenieren.
• Viskosität: Anstatt Schockviskosität (die die Ergebnisse verfälschen könnte) verwenden die Autoren eine konstante kinematische Viskosität ($\nu$).
• Auflösung und Parameter: Die Simulationen werden bei Auflösungen von $512^3$ bis $2048^3$ Gitterpunkten durchgeführt. Die Autoren variieren die initiale Geschwindigkeitsamplitude ($u_{ini}$) und die Viskosität ($\nu$), um numerisch zuverlässige Grenzwerte zu bestimmen.
• Analyse: Die Studie verfolgt die Entwicklung der Enstrophie ($\langle \omega^2 \rangle/2$) und zerlegt deren Zeitableitung in drei Terme: einen dynamo-ähnlichen Term ($T_{dyn}$), einen durch Viskosität und Dichtegradienten getriebenen Erzeugungs-Term ($T_{gen}$) und einen dissipativen Term ($T_{dis}$). Die Autoren analysieren zudem Geschwindigkeits-, Dichte- und Vortizitätsspektren und untersuchen die Skalierung der Vortizität mit der initialen Strömungsamplitude.

Hauptergebnisse

1. Viskositätsabhängigkeit: Die Produktion von Vortizität wird als ein Sekundäreffekt identifiziert, der durch die Verstärkung initialer Flussdivergenzen mittels Viskosität getrieben wird. Die Studie stellt fest, dass das Verhältnis von Vortizitätserzeugung zu Dissipation ($R_{gen}$) nicht proportional zur Viskosität ist; vielmehr ist die Vortizitätserzeugung in der frühen Phase vor dem Hauptkollaps am signifikantesten und nicht als Folge des Kollapsflusses selbst.
2. Skalierung mit den Anfangsbedingungen: Die quadratische Mittelwert-Vortizität ($\omega_{rms}$) skaliert mit dem Quadrat der initialen Geschwindigkeitsamplitude ($u_{ini}^2$), während die quadratische Mittelwert-Geschwindigkeitsdivergenz linear mit $u_{ini}$ skaliert. Diese Skalierung bleibt während der gesamten Simulation bestehen, was darauf hindeutet, dass der Gravitationskollaps die Vortizitätsproduktion nicht direkt antreibt. Stattdessen verstärkt der Kollaps die initialen Flussdivergenzen, welche die Viskosität dann in Vortizität umwandelt.
3. Numerische Zuverlässigkeit: Die Autoren stellen fest, dass die Ergebnisse unzuverlässig werden, wenn die Gitter-Reynolds-Zahl ($Re_{mesh}$) Werte zwischen 0,7 und 2,0 überschreitet, abhängig vom spezifischen Durchlauf. Im zuverlässigen Regime bleibt der dynamo-ähnliche Term ($T_{dyn}$) subdominant und trägt weniger als 20 % zur Enstrophieentwicklung bei.
4. Spektralanalyse: Das Vortizitätsspektrum zeigt einen $k^4$-Subinertialbereich, analog zum Verhalten magnetischer Felder in der MHD, doch die Autoren führen dies eher auf die Jeans-Instabilität zurück, die zufällige Geschwindigkeitsfluktuationen antreibt, als auf einen Dynamo-Prozess. Das Geschwindigkeitsspektrum zeigt einen steilen Inertialbereich ohne „Bottleneck“-Effekt.
5. Kritische Vortizität für den Dynamo: Bei der Überprüfung der kritischen magnetischen Vortizitäts-Reynolds-Zahl ($Pr_M Re_\omega$) für den Dynamo-Beginn bestätigen die Autoren eine Bifurkation um einen Wert von 300. Dies impliziert, dass rein irrotationale Strömungen, die ohne externe Antriebskräfte oder spezifische Anfangsbedingungen nicht genügend Vortizität aufrechterhalten können, nicht als Dynamos fungieren würden.

Bedeutung und Behauptungen
Die Arbeit kommt zu dem Schluss, dass Gravitationskollaps unter den Bedingungen einer barotropen Zustandsgleichung und anfänglich irrotationaler Turbulenz keine effiziente Vortizität produziert. Die beobachtete Vortizität ist ein Sekundäreffekt, der aus der viskosen Umwandlung initialer Flussdivergenzen resultiert, und nicht eine direkte Folge der Kollapsdynamik.

Die Autoren behaupten, dass ihre Ergebnisse bisherige Behauptungen über kollapsinduzierte kleinskalige Dynamo-Wirkung infrage stellen. Sie legen nahe, dass frühere Beobachtungen der Magnetfeldverstärkung in kollapsierenden Strömungen wahrscheinlich durch das Abklingen der initialen Turbulenz angetrieben wurden und nicht durch den Kollaps selbst. Die Studie betont die Notwendigkeit direkter numerischer Simulationen ohne Subgrid-Scale-Modellierung, um künstliche Vortizitätserzeugung zu vermeiden. Die Arbeit impliziert, dass das Erreichen eines selbsterhaltenden Dynamos in solchen Kollapsszenarien schwierig ist und dass der Kollapsfluss allein nicht ausreicht, um die notwendige vortikale Turbulenz für eine Dynamo-Wirkung zu erzeugen.