Box-Simulationen von rotierender Magnetokonvektion im flüssigen Erdkern
Numerische Simulationen der 3D-MHD Gleichungen sind mit Hilfe des Codes NIRVANA durchgeführt worden. Die Gleichungen für kompressible rotierende Magnetokonvektion wurden für erdähnliche Bedingungen numerisch in einer kartesischen Box gelöst.
Charakteristische Eigenschaften mittlerer Größen, wie der Turbulenz-Intensität oder der turbulente Wärmefluss, die durch die kombinierte Wirkung kleinskaliger Fluktuationen entstehen, wurden bestimmt. Die Korrelationslänge der Turbulenz hängt signifikant von der Stärke und der Orientierung des Magnetfeldes ab, und das anisotrope Verhalten der Turbulenz aufgrund von Coriolis- und Lorentzkraft ist für schnellere Rotation wesentlich stärker ausgeprägt. Die Ausbildung eines isotropen Verhaltens auf kleinen Skalen unter dem Einfluss von Rotation alleine wird bereits durch ein schwaches Magnetfeld verhindert. Dies resultiert in einer turbulenten Strömung, die durch die vertikale Komponente dominiert wird. In Gegenwart eines horizontalen Magnetfeldes nimmt der vertikale turbulente Wärmefluss leicht mit zunehmender Feldstärke zu, so dass die Kühlung eines rotierenden Systems verbessert wird. Der horizontale Wärmetransport ist stets westwärts und in Richtung der Pole orientiert. Letzteres kann unter Umständen die Quelle für eine großskalige meridionale Strömung darstellen, während erstes in globalen Simulationen mit nicht axialsymmetrischen Randbedingungen für den Wärmefluss von Bedeutung ist.
Die mittlere elektromotorische Kraft, die die Erzeugung von magnetischem Fluss
durch die Turbulenz beschreibt, wurde unmittelbar aus den Lösungen für Geschwindigkeit und Magnetfeld berechnet. Hieraus konnten die entsprechenden α-Koeffizienten hergeleitet werden. Aufgrund der sehr schwachen Dichtestratifizierung ändert der α-Effekt sein Vorzeichen nahezu exakt in der Mitte der Box. Der α-Effekt ist positiv in der oberen Hälfte und negativ in der unteren Hälfte einer auf der Nordhalbkugel rotierenden Box. Für ein starkes Magnetfeld ergibt sich zudem eine deutliche abwärts orientierte Advektion von magnetischem Fluss.
Ein Mean-Field Modell des Geodynamos wurde konstruiert, das auf dem α-Effekt basiert, wie er aus den Box-Simulationen berechnet wurde. Für eine äußerst beschränkte Klasse von radialen α-Profilen weist das lineare α^2-Modell
Oszillationen auf einer Zeitskala auf, die durch die turbulente Diffusionszeit
bestimmt wird. Die wesentlichen Eigenschaften der periodischen Lösungen werden präsentiert, und der Einfluss der Größe des inneren Kerns auf die Charakteristiken des kritischen Bereichs, innerhalb dessen oszillierende Lösungen auftreten, wurden untersucht. Reversals werden als eine halbe Oszillation interpretiert. Sie sind ein recht seltenes Ereignis, da sie lediglich dann stattfinden können, wenn das α-Profil ausreichend lange in dem periodische Lösungen erlaubenden Bereich liegt. Aufgrund starker Fluktuationen auf der konvektiven Zeitskala ist die Wahrscheinlichkeit eines solchen Reversals relativ klein. In einem einfachen nicht-linearen Mean-Field Modell mit realistischen Eingabeparametern, die auf den Box-Simulationen beruhen, konnte die Plausibilität des Reversal-Modells anhand von Langzeitsimulationen belegt werden. / Box-simulations of rotating magnetoconvection in the Earth's fluid core.
Simulations of the 3D MHD-equations have been performed using the code NIRVANA. The equations for compressible rotating magnetoconvection are numerically solved in a Cartesian box assuming conditions roughly suitable for the geodynamo.
The characteristics of averaged quantities like the turbulence intensity and the turbulent heat flux that are caused by the combined action of the small-scale fluctuations are computed. The correlation length of the turbulence significantly depends on the strength and orientation of the magnetic field and the anisotropic behavior of the turbulence intensity induced by Coriolis and Lorentz force is considerably more pronounced for faster rotation. The development of an isotropic behavior on the small scales -- as it is observed in pure rotating convection -- vanishes even for weak magnetic fields, which results in a turbulent flow that is dominated by the vertical component. In the presence of a horizontal magnetic field the vertical turbulent heat flux slightly increases with increasing field strength, so that cooling of the rotating system is facilitated. Horizontal transport of heat is always directed westwards and towards the poles. The latter might be a source of a large-scale meridional flow whereas the first would be important in global simulations in case of non-axisymmetric boundary conditions for the heat flux.
The mean electromotive force describing the generation of mean magnetic flux by turbulence in the rotating fluid is directly calculated from the simulations and the corresponding α-coefficients are derived. Due to the very weak density stratification, the α-effect changes its sign in the middle of the box. It is positive at the top and negative at the bottom of the convective instable layer. For strong magnetic fields we also find a clear downward advection of the mean magnetic field. Finally the quenching behavior of the α-effect in dependence of the imposed magnetic field strength is presented.
A geodynamo-model is constructed, which is based on an α-effect that has been computed from the box simulations. For a highly restricted class of radial α-profiles the linear α^2-model exhibits oscillating solutions on a timescale given by the turbulent diffusion time. The basic properties of the periodic solutions are presented and the influence of the inner core size on the characteristics of the critical range that allows for oscillating solutions is shown. Reversals are interpreted as half of such an oscillation. They are rather seldom events because they can only occur if the α-profile exists long enough within the small critical range that allows for periodic solutions. Due to strong fluctuations on the convective timescale the probability of such a reversal is very small. Finally, a simple non-linear mean-field model with reasonable input parameters based on the box-simulations demonstrates the plausibility of the presented theory with a long-time series of a (geo-)dynamo reversal sequence.
| Identifer | oai:union.ndltd.org:Potsdam/oai:kobv.de-opus-ubp:1360 |
| Date | January 2007 |
| Creators | Giesecke, André |
| Publisher | Universität Potsdam, Mathematisch-Naturwissenschaftliche Fakultät. Institut für Physik und Astronomie |
| Source Sets | Potsdam University |
| Language | German |
| Detected Language | English |
| Type | Text.Thesis.Doctoral |
| Format | application/pdf |
| Rights | http://creativecommons.org/licenses/by-nc-nd/2.0/de/ |
Page generated in 0.0022 seconds