Petri J. Käpylä

Cyclic Magnetic Activity due to Turbulent Convection in Spherical Wedge Geometry

We report on simulations of turbulent, rotating, stratified, magnetohydrodynamic convection in spherical wedge geometry. An initially small-scale, random, weak-amplitude magnetic field is amplified by several orders of magnitude in the course of the simulation to form oscillatory large-scale fields in the saturated state of the dynamo. The differential rotation is solar-like (fast equator), but neither coherent meridional poleward circulation nor near-surface shear layer develop in these runs. In addition to a poleward branch of magnetic activity beyond 50 degrees latitude, we find for the first time a pronounced equatorward branch at around 20 degrees latitude, reminiscent of the solar cycle.

MAGNETIC DIFFUSIVITY TENSOR AND DYNAMO EFFECTS IN ROTATING AND SHEARING TURBULENCE

The turbulent magnetic diffusivity tensor is determined in the presence of rotation or shear. The question is addressed whether dynamo action from the shear-current effect can explain large-scale magnetic field generation found in simulations with shear. For this purpose a set of evolution equations for the response to imposed test fields is solved with turbulent and mean motions calculated from the momentum and continuity equations. The corresponding results for the electromotive force are used to calculate turbulent transport coefficients. The diagonal components of the turbulent magnetic diffusivity tensor are found to be very close together, but their values increase slightly with increasing shear and decrease with increasing rotation rate. In the presence of shear, the sign of the two off-diagonal components of the turbulent magnetic diffusion tensor is the same and opposite to the sign of the shear. This implies that dynamo action from the shear-current effect is impossible, except perhaps for high magnetic Reynolds numbers. However, even though there is no alpha effect on the average, the components of the α tensor display Gaussian fluctuations around zero. These fluctuations are strong enough to drive an incoherent alpha-shear dynamo. The incoherent shear-current effect, on the other hand, is found to be subdominant.

Large-scale dynamos in turbulent convection with shear

Aims. To study the existence of large-scale convective dynamos under the influence of shear and rotation. Methods. Three-dimensional numerical simulations of penetrative compressible convection with uniform horizontal shear are used to study dynamo action and the generation of large-scale magnetic fields. We consider cases where the magnetic Reynolds number is either marginal or moderately supercritical with respect to small-scale dynamo action in the absence of shear and rotation. Our magnetic Reynolds number is based on the wavenumber of the depth of the convectively unstable layer. The effects of magnetic helicity fluxes are studied by comparing results for the magnetic field with open and closed boundaries. Results. Without shear no large-scale dynamos are found even if the ingredients necessary for the α-effect (rotation and stratification) are present in the system. When uniform horizontal shear is added, a large-scale magnetic field develops, provided the boundaries are open. In this case the mean magnetic field contains a significant fraction of the total field. For those runs where the magnetic Reynolds number is between 60 and 250, an additional small-scale dynamo is expected to be excited, but the field distribution is found to be similar to cases with smaller magnetic Reynolds number where the small-scale dynamo is not excited. In the case of closed (perfectly conducting) boundaries, magnetic helicity fluxes are suppressed and no large-scale fields are found. Similarly, poor large-scale field development is seen when vertical shear is used in combination with periodic boundary conditions in the horizontal directions. If, however, open (normal-field) boundary conditions are used in the x-direction, a large-scale field develops. These results support the notion that shear not only helps to generate the field, but it also plays a crucial role in driving magnetic helicity fluxes out of the system along the isocontours of shear, thereby allowing efficient dynamo action.

Reynolds stress and heat flux in spherical shell convection

Context. Turbulent fluxes of angular momentum and enthalpy or heat due to rotationally affected convection play a key role in determining differential rotation of stars. Their dependence on latitude and depth has been determined in the past from convection simulations in Cartesian or spherical simulations. Here we perform a systematic comparison between the two geometries as a function of the rotation rate. Aims. Here we want to extend the earlier studies by using spherical wedges to obtain turbulent angular momentum and heat transport as functions of the rotation rate from stratified convection. We compare results from spherical and Cartesian models in the same parameter regime in order to study whether restricted geometry introduces artefacts into the results. In particular, we want to clarify whether the sharp equatorial profile of the horizontal Reynolds stress found in earlier Cartesian models is also reproduced in spherical geometry. Methods. We employ direct numerical simulations of turbulent convection in spherical and Cartesian geometries. In order to alleviate the computational cost in the spherical runs, and to reach as high spatial resolution as possible, we model only parts of the latitude and longitude. The rotational influence, measured by the Coriolis number or inverse Rossby number, is varied from zero to roughly seven, which is the regime that is likely to be realised in the solar convection zone. Cartesian simulations are performed in overlapping parameter regimes. Results. For slow rotation we find that the radial and latitudinal turbulent angular momentum fluxes are directed inward and equatorward, respectively. In the rapid rotation regime the radial flux changes sign in accordance with earlier numerical results, but in contradiction with theory. The latitudinal flux remains mostly equatorward and develops a maximum close to the equator. In Cartesian simulations this peak can be explained by the strong banana cells. Their effect in the spherical case does not appear to be as large. The latitudinal heat flux is mostly equatorward for slow rotation but changes sign for rapid rotation. Longitudinal heat flux is always in the retrograde direction. The rotation profiles vary from anti-solar (slow equator) for slow and intermediate rotation to solar-like (fast equator) for rapid rotation. The solar-like profiles are dominated by the Taylor-Proudman balance.

Effects of enhanced stratification on equatorward dynamo wave propagation

We present results from simulations of rotating magnetized turbulent convection in spherical wedge geometry representing parts of the latitudinal and longitudinal extents of a star. Here we consider a set of runs for which the density stratification is varied, keeping the Reynolds and Coriolis numbers at similar values. In the case of weak stratification, we find quasi-steady dynamo solutions for moderate rotation and oscillatory ones with poleward migration of activity belts for more rapid rotation. For stronger stratification, the growth rate tends to become smaller. Furthermore, a transition from quasi-steady to oscillatory dynamos is found as the Coriolis number is increased, but now there is an equatorward migrating branch near the equator. The breakpoint where this happens corresponds to a rotation rate that is about three to seven times the solar value. The phase relation of the magnetic field is such that the toroidal field lags behind the radial field by about ?/2, which can be explained by an oscillatory ?2 dynamo caused by the sign change of the ?-effect about the equator. We test the domain size dependence of our results for a rapidly rotating run with equatorward migration by varying the longitudinal extent of our wedge. The energy of the axisymmetric mean magnetic field decreases as the domain size increases and we find that an m = 1 mode is excited for a full 2? azimuthal extent, reminiscent of the field configurations deduced from observations of rapidly rotating late-type stars.

Non-Fickian diffusion and tau approximation from numerical turbulence

Evidence for non-Fickian diffusion of a passive scalar is presented using direct simulations of homogeneous isotropic turbulence. The results compare favorably with an explicitly time-dependent closure model based on the tau approximation. In the numerical experiments three different cases are considered: (i) zero mean concentration with finite initial concentration flux, (ii) an initial top hat profile for the concentration, and (iii) an imposed background concentration gradient. All cases agree in the resulting relaxation time in the tau approximation relating the triple correlation to the concentration flux. The first-order smoothing approximation is shown to be inapplicable.

Alpha effect and turbulent diffusion from convection

Aims. We study turbulent transport coefficients that describe the evolution of large-scale magnetic fields in turbulent convection. Methods. We use the test field method, together with three-dimensional numerical simulations of turbulent convection with shear and rotation, to compute turbulent transport coefficients describing the evolution of large-scale magnetic fields in mean-field theory in the kinematic regime. We employ one-dimensional mean-field models with the derived turbulent transport coefficients to examine whether they give results that are compatible with direct simulations. Results. The results for the α-effect as a function of rotation rate are consistent with earlier numerical studies, i.e. increasing magnitude as rotation increases and approximately cos θ latitude profile for moderate rotation. Turbulent diffusivity, ηt, is proportional to the square of the turbulent vertical velocity in all cases. Whereas ηt decreases approximately inversely proportional to the wavenumber of the field, the α-effect and turbulent pumping show a more complex behaviour with partial or full sign changes and the magnitude staying roughly constant. In the presence of shear and no rotation, a weak α-effect is induced which does not seem to show any consistent trend as a function of shear rate. Provided that the shear is large enough, this small α-effect is able to excite a dynamo in the mean-field model. The coefficient responsible for driving the shear-current effect shows several sign changes as a function of depth but is also able to contribute to dynamo action in the mean-field model. The growth rates in these cases are, however, well below those in direct simulations, suggesting that an incoherent α-shear dynamo may also act in the simulations. If both rotation and shear are present, the α-effect is more pronounced. At the same time, the combination of the shear-current and Ω × J-effects is also stronger than in the case of shear alone, but subdominant to the α-shear dynamo. The results of direct simulations are consistent with mean-field models where all of these effects are taken into account without the need to invoke incoherent effects.

Convective dynamos in spherical wedge geometry

Astronomy Unit, School of Mathematical Sciences, Queen Mary University of London, Mile End Road, London E14NS, United KingdomReceived 2009 Sep 9, accepted 2009 Nov 16Published online 2009 Dec 30Key words Sun: magnetic fields – magnetohydrodynamics (MHD)Self-consistent convective dynamo simulations in wedge-shaped spherical shells are presented. Differential rotation isgenerated by the interaction of convection with rotation. Equatorward acceleration and dynamo action are obtained onlyfor sufficiently rapid rotation. The angular velocity tends to be constant along cylinders. Oscillatory large-scale fields arefound to migrate in the poleward direction. Comparison with earlier simulations in full spherical shells and Cartesiandomains is made.

Turbulent Dynamos with Shear and Fractional Helicity

Dynamo action owing to helically forced turbulence and large-scale shear is studied using direct numerical simulations. The resulting magnetic field displays propagating wave-like behavior. This behavior can be modeled in terms of an αΩ dynamo. In most cases super-equipartition fields are generated. By varying the fraction of helicity of the turbulence the regeneration of poloidal fields via the helicity effect (corresponding to the α-effect) is regulated. The saturation level of the magnetic field in the numerical models is consistent with a linear dependence on the ratio of the fractional helicities of the small and large-scale fields, as predicted by a simple nonlinear mean-field model. As the magnetic Reynolds number (Re M ) based on the wavenumber of the energy-carrying eddies is increased from 1 to 180, the cycle frequency of the large-scale field is found to decrease by a factor of about 6 in cases where the turbulence is fully helical. This is interpreted in terms of the turbulent magnetic diffusivity, which is found to be only weakly dependent on the Re M .

OSCILLATORY MIGRATING MAGNETIC FIELDS IN HELICAL TURBULENCE IN SPHERICAL DOMAINS

We present direct numerical simulations of the equations of compressible magnetohydrodynamics in a wedge-shaped spherical shell, without shear, but with random helical forcing which has negative (positive) helicity in the northern (southern) hemisphere. We find a large-scale magnetic field that is nearly uniform in the azimuthal direction and approximately antisymmetric about the equator. Furthermore, the large-scale field in each hemisphere oscillates on nearly dynamical timescales with reversals of polarity and equatorward migration. Corresponding mean-field models also show similar migratory oscillations with a frequency that is nearly independent of the magnetic Reynolds number. This mechanism may be relevant for understanding equatorward migration seen in the solar dynamo.