Reynolds stresses and mean fields generated by pure waves: applications to shear flows and convection in a rotating shell

A general reformulation of the Reynolds stresses created by two-dimensional waves breaking a translational or a rotational invariance is described. This reformulation emphasizes the importance of a geometrical factor: the slope of the separatrices of the wave flow. Its physical relevance is illustrated by two model systems: waves destabilizing open shear flows; and thermal Rossby waves in spherical shell convection with rotation. In the case of shear-flow waves, a new expression of the Reynolds–Orr amplification mechanism is obtained, and a good understanding of the form of the mean pressure and velocity fields created by weakly nonlinear waves is gained. In the case of thermal Rossby waves, results of a three-dimensional code using no-slip boundary conditions are presented in the nonlinear regime, and compared with those of a two-dimensional quasi-geostrophic model. A semi-quantitative agreement is obtained on the flow amplitudes, but discrepancies are observed concerning the nonlinear frequency shifts. With the quasi-geostrophic model we also revisit a geometrical formula proposed by Zhang to interpret the form of the zonal flow created by the waves, and explore the very low Ekman-number regime. A change in the nature of the wave bifurcation, from supercritical to subcritical, is found

Planetary Dynamos

The MAGSAT-program has added significantly to our knowledge of planetary magnetism. The accuracy of observations has been improved such that a reliable extrapolation of the magnetic field to the core surface is now much more feasible than it has been before, and the prospect of further MAGSAT missions raises the expectation that the time dependence of the geomagnetic field will be known with similar accuracy in the future. In the research support it has been attempted to develop dynamo theory with these applications in mind

Bistability and hysteresis of dipolar dynamos generated by turbulent convection in rotating spherical shells

Bistability and hysteresis of magnetohydrodynamic dipolar dynamos generated by turbulent convection in rotating spherical fluid shells is demonstrated. Hysteresis appears as a transition between two distinct regimes of dipolar dynamos with rather different properties including a pronounced difference in the amplitude of the axisymmetric poloidal field component and in the form of the differential rotation. The bistability occurs from the onset of dynamo action up to about 9 times the critical value of the Rayleigh number for onset of convection and over a wide range of values of the ordinary and the magnetic Prandtl numbers including the value unity

Quasi-geostrophic approximation of anelastic convection

The onset of convection in a rotating cylindrical annulus with parallel ends filled with a compressible fluid is studied in the anelastic approximation. Thermal Rossby waves propagating in the azimuthal direction are found as solutions. The analogy to the case of Boussinesq convection in the presence of conical end surfaces of the annular region is emphasised. As in the latter case, the results can be applied as an approximation for the description of the onset of anelastic convection in rotating spherical fluid shells. Reasonable agreement with three-dimensional numerical results published by Jones, Kuzanyan & Mitchell (J. Fluid Mech., vol. 634, 2009, pp. 291–319) for the latter problem is found. As in those results, the location of the onset of convection shifts outwards from the tangent cylinder with increasing number Nρof density scale heights until it reaches the equatorial boundary. A new result is that at a much higher number Nρ the onset location returns to the interior of the fluid shell

Magneto-inertial convection in rotating fluid spheres

The onset of convection in the form of magneto-inertial waves in a rotating fluid sphere permeated by a constant axial electric current is studied through a perturbation analysis. Explicit expressions for the dependence of the Rayleigh number on the azimuthal wavenumber are derived in the limit of high thermal diffusivity. Results for the cases of thermally infinitely conducting and of nearly thermally insulating boundaries are obtained.Comment: 10 pages, 5 figures, to be submitted for publicatio

Baroclinically-driven flows and dynamo action in rotating spherical fluid shells

The dynamics of stably stratified stellar radiative zones is of considerable interest due to the availability of increasingly detailed observations of Solar and stellar interiors. This article reports the first non-axisymmetric and time-dependent simulations of flows of anelastic fluids driven by baroclinic torques in stably stratified rotating spherical shells – a system serving as an elemental model of a stellar radiative zone. With increasing baroclinicity a sequence of bifurcations from simpler to more complex flows is found in which some of the available symmetries of the problem are broken subsequently. The poloidal component of the flow grows relative to the dominant toroidal component with increasing baroclinicity. The possibility of magnetic field generation thus arises and this paper proceeds to provide some indications for self-sustained dynamo action in baroclinically-driven flows. We speculate that magnetic fields in stably stratified stellar interiors are thus not necessarily of fossil origin as it is often assumed

Numerical simulations of rotating sunspots

A numerical model of idealized, axisymmetric, rotating sunspots is presented. The model contains a compressible plasma described by the nonlinear MHD equations, with density and temperature gradients simulating the upper layer of the sun’s convection zone. The solution forms a central flux tube in the cylindrical numerical domain, with convection cells pushing the magnetic field to the axis. When the numerical domain is rotated with a constant angular velocity, the umbra rotates as a rigid body while the surrounding convection cells show a swirling, vortical flow. As a result, the azimuthal velocity and magnetic field have their maximum values close to the flux tube, inside the innermost convection cell

Dynamo Effects Near The Transition from Solar to Anti-Solar Differential Rotation

Numerical MHD simulations play increasingly important role for understanding mechanisms of stellar magnetism. We present simulations of convection and dynamos in density-stratified rotating spherical fluid shells. We employ a new 3D simulation code for the solution of a physically consistent anelastic model of the process with a minimum number of parameters. The reported dynamo simulations extend into a "buoyancy-dominated" regime where the buoyancy forcing is dominant while the Coriolis force is no longer balanced by pressure gradients and strong anti-solar differential rotation develops as a result. We find that the self-generated magnetic fields, despite being relatively weak, are able to reverse the direction of differential rotation from anti-solar to solar-like. We also find that convection flows in this regime are significantly stronger in the polar regions than in the equatorial region, leading to non-oscillatory dipole-dominated dynamo solutions, and to concentration of magnetic field in the polar regions. We observe that convection has different morphology in the inner and at the outer part of the convection zone simultaneously such that organized geostrophic convection columns are hidden below a near-surface layer of well-mixed highly-chaotic convection. While we focus the attention on the buoyancy-dominated regime, we also demonstrate that conical differential rotation profiles and persistent regular dynamo oscillations can be obtained in the parameter space of the rotation-dominated regime even within this minimal model.Comment: Published in the Astrophysical Journa

Lagrangian Statistics of Navier-Stokes- and MHD-Turbulence

We report on a comparison of high-resolution numerical simulations of Lagrangian particles advected by incompressible turbulent hydro- and magnetohydrodynamic (MHD) flows. Numerical simulations were performed with up to $1024^3$ collocation points and 10 million particles in the Navier-Stokes case and $512^3$ collocation points and 1 million particles in the MHD case. In the hydrodynamics case our findings compare with recent experiments from Mordant et al. [1] and Xu et al. [2]. They differ from the simulations of Biferale et al. [3] due to differences of the ranges choosen for evaluating the structure functions. In Navier-Stokes turbulence intermittency is stronger than predicted by a multifractal approach of [3] whereas in MHD turbulence the predictions from the multifractal approach are more intermittent than observed in our simulations. In addition, our simulations reveal that Lagrangian Navier-Stokes turbulence is more intermittent than MHD turbulence, whereas the situation is reversed in the Eulerian case. Those findings can not consistently be described by the multifractal modeling. The crucial point is that the geometry of the dissipative structures have different implications for Lagrangian and Eulerian intermittency. Application of the multifractal approach for the modeling of the acceleration PDFs works well for the Navier-Stokes case but in the MHD case just the tails are well described.Comment: to appear in J. Plasma Phy