The influence of fluctuating conductivity on the coefficients known from the mean-field electrodynamics is considered. If the conductivity fluctuations are assumed as uncorrelated with the turbulent velocity field then only the effective magnetic diffusivity of the fluid is reduced and the decay time of a large-scale magnetic field is increased. If the fluctuations of conductivity and flow are correlated in a certain direction then an additional diamagnetic pumping effect results transporting magnetic field in opposite direction to the resistivity flux vector $\langle \eta'\vec{u}'\rangle$. Even for homogeneous turbulence fields in the presence of rotation an alpha effect appears. With the characteristic values of the outer core of the Earth or the solar convection zone, however, the dynamo number of the alpha effect never reaches supercritical values to operate as an $\alpha^2$-dynamo.

The possibility of explaining shear flow dynamos in terms of a magnetic shear-current effect (SC) is examined. A competing explanation is the incoherent $\alpha$-shear dynamo effect. Our primary diagnostics is the turbulent magnetic diffusivity tensor, in particular the off-diagonal diffusivity tensor component $\eta_{yx}$, a systematically negative sign of which would imply coherent dynamo action through SC. To be able to measure turbulent transport coefficients from systems with strong magnetic fluctuations, we present an extension of the test-field method that is capable of such measurements in the case where the pressure gradient term is dropped from the MHD equations, which is a nonlinear TFM (NLTFM). The hydrodynamic equation is related to Burger's equation and the resulting flows are referred to as magnetized burgulence. We use both kinetic and magnetic forcings, to mimic cases without and with simultaneous small-scale dynamo action. When the simplified MHD flow is forced kinetically, negative $\eta_{yx}$ values are obtained with exponential growth in both the radial and azimuthal magnetic field components. Using stochastic monochromatic magnetic forcing, the exponential growth is no longer seen, and NLTFM yields positive values of $\eta_{yx}$. In an attempt to recover the exponential growth seen in the kinetically forced case, we also employ an alternative forcing function with lowest $k$ vectors being removed - in this case the exponential growth is recovered, but the NLTFM results do not change dramatically. Our analysis of the dynamo numbers for the coherent SC and incoherent $\alpha$ and SC effects show that the incoherent effects are the main drivers of the dynamo instability in majority of cases. The coherent SC effect is mildly enhancing the dynamo action in the kinetically forced cases, while we find no evidence for magnetic SC in our simulations.

One interpretation of the activity and magnetism of late-type stars is that these both intensify with decreasing Rossby number up to a saturation level, suggesting that stellar dynamos depend on both rotation and convective turbulence. Some studies have claimed, however, that rotation alone suffices to parametrise this scaling adequately. Here, we tackle the question of the relevance of turbulence to stellar dynamos by including evolved, post main sequence stars in the analysis of the rotation-activity relation. These stars rotate very slowly compared with main sequence stars, but exhibit similar activity levels. We show that the two evolutionary stages fall together in the rotation-activity diagram and form a single sequence in the unsaturated regime in relation only to Rossby numbers derived from stellar models, confirming earlier preliminary results that relied on a more simplistic parametrisation of the convective turnover time. This mirrors recent results of fully convective M dwarfs, which likewise fall on the same rotation-activity sequence as partially convective solar-type stars. Our results demonstrate that turbulence plays a crucial role in driving stellar dynamos and suggest that there is a common turbulence-related dynamo mechanism explaining the magnetic activity of all late-type stars.

The $\Lambda$ effect describes a rotation-induced non-diffusive contribution to the Reynolds stress. It is commonly held responsible for maintaining the observed differential rotation of the Sun and other late-type stars. Here the sensitivity of the $\Lambda$ effect to small-scale magnetic fields and compressibility is studied by means of forced turbulence simulations either with anisotropic forcing in fully periodic cubes or in density-stratified domains with isotropic forcing. Effects of small-scale magnetic fields are studied in cases where the magnetic fields are self-consistently generated by a small-scale dynamo. The results show that small-scale magnetic fields lead to a quenching of the $\Lambda$ effect which is milder than in cases where also a large-scale field is present. The effect of compressibility on the $\Lambda$ effect is negligible in the range of Mach numbers from 0.015 to 0.8. Density stratification induces a marked anisotropy in the turbulence and a vertical $\Lambda$ effect if the forcing scale is roughly two times larger than the density scale height.

(abridged) Context: Turbulent diffusion of large-scale flows and magnetic fields play major roles in many astrophysical systems. Aims: Our goal is to compute turbulent viscosity and magnetic diffusivity, relevant for diffusing large-scale flows and magnetic fields, respectively, and their ratio, the turbulent magnetic Prandtl number, ${\rm Pm}_{\rm t}$, for isotropically forced homogeneous turbulence. Methods: We use simulations of forced turbulence in fully periodic cubes composed of isothermal gas with an imposed large-scale sinusoidal shear flow. Turbulent viscosity is computed either from the resulting Reynolds stress or from the decay rate of the large-scale flow. Turbulent magnetic diffusivity is computed using the test-field method. The scale dependence of the coefficients is studied by varying the wavenumber of the imposed sinusoidal shear and test fields. Results: We find that turbulent viscosity and magnetic diffusivity are in general of the same order of magnitude. Furthermore, the turbulent viscosity depends on the fluid Reynolds number (${\rm Re}$) and scale separation ratio of turbulence. The scale dependence of the turbulent viscosity is found to be well approximated by a Lorentzian. The results for the turbulent transport coefficients appear to converge at sufficiently high values of ${\rm Re}$ and the scale separation ratio. However, a weak decreasing trend is found even at the largest values of ${\rm Re}$. The turbulent magnetic Prandtl number converges to a value that is slightly below unity for large ${\rm Re}$ whereas for small ${\rm Re}$, we find values between 0.5 and 0.6. Conclusions: The turbulent magnetic diffusivity is in general consistently higher than the turbulent viscosity. The actual value of ${\rm Pm}_{\rm t}$ found from the simulations ($\approx0.9\ldots0.95$) at large ${\rm Re}$ and scale separation ratio is higher than any of the analytic predictions.

(abridged) Context: Convective motions overshooting to regions that are formally convectively stable cause extended mixing. Aims: To determine the scaling of overshooting depth ($d_{\rm os}$) at the base of the convection zone as a function of imposed energy flux ($\mathscr{F}_{\rm n}$) and to estimate the extent of overshooting at the base of the solar convection zone. Methods: Three-dimensional Cartesian simulations of compressible non-rotating convection with unstable and stable layers are used. The simulations use either a fixed heat conduction profile or a temperature and density dependent formulation based on Kramers opacity law. The simulations cover a range of almost four orders of magnitude in the imposed flux. Results: A smooth heat conduction profile (either fixed or through Kramers opacity law) leads to a relatively shallow power law with $d_{\rm os}\propto \mathscr{F}_{\rm n}^{0.08}$ for low $\mathscr{F}_{\rm n}$. A fixed step-profile of the heat conductivity at the bottom of the convection zone leads to a somewhat steeper dependency with $d_{\rm os}\propto \mathscr{F}_{\rm n}^{0.12}$. Experiments with and without subgrid-scale entropy diffusion revealed a strong dependence on the effective Prandtl number which is likely to explain the steep power laws as a function of $\mathscr{F}_{\rm n}$ reported in the literature. Furthermore, changing the heat conductivity artificially below the convection zone is shown to lead to substantial underestimation of overshooting depth. Conclusions: Extrapolating from the results obtained with smooth heat conductivity profiles suggest that the overshooting depth for the solar flux is of the order of $0.2$ pressure scale heights at the base of the convection zone which is two to four times higher than estimates from helioseismology.

Global and semi-global convective dynamo simulations of solar-like stars are known to show a transition from an anti-solar (fast poles, slow equator) to solar-like (fast equator, slow poles) differential rotation (DR) for increasing rotation rate. The dynamo solutions in the latter regime can exhibit regular cyclic modes, whereas in the former one, only stationary or temporally irregular solutions have been obtained so far. In this paper we present a semi-global dynamo simulation in the transition region, exhibiting two coexisting dynamo modes, a cyclic and a stationary one, both being dynamically significant. We seek to understand how such a dynamo is driven by analyzing the large-scale flow properties (DR and meridional circulation) together with the turbulent transport coefficients obtained with the test-field method. Neither an $\alpha\Omega$ dynamo wave nor an advection-dominated dynamo are able to explain the cycle period and the propagation direction of the mean magnetic field. Furthermore, we find that the $\alpha$ effect is comparable or even larger than the $\Omega$ effect in generating the toroidal magnetic field, and therefore, the dynamo seems to be $\alpha^2\Omega$ or $\alpha^2$ type. We further find that the effective large-scale flows are significantly altered by turbulent pumping.

Eclipsing time variations have been observed for a wide range of binary systems, including post-common-envelope binaries. A frequently proposed explanation, apart from the possibility of having a third body, is the effect of magnetic activity, which may alter the internal structure of the secondary star, particularly its quadrupole moment, and thereby cause quasi-periodic oscillations. Here we present two compressible non-ideal magneto-hydrodynamical (MHD) simulations of the magnetic dynamo in a solar mass star, one of them with three times the solar rotation rate ("slow rotator"), the other one with twenty times the solar rotation rate ("rapid rotator"), to account for the high rotational velocities in close binary systems. For the slow rotator, we find that both the magnetic field and the stellar quadrupole moment change in a quasi-periodic manner, leading to O-C (observed - corrected times of the eclipse) variations of ~0.025 s. For the rapid rotator, the behavior of the magnetic field as well as the quadrupole moment changes become considerably more complex, due to the less coherent dynamo solution. The resulting O-C variations are of the order 0.13 s. The observed system V471~Tau shows two modes of eclipsing time variations, with amplitudes of 151 s and 20 s, respectively. However, the current simulations may not capture all relevant effects due to the neglect of the centrifugal force and self-gravity. Considering the model limitations and that the rotation of V471 Tau is still a factor of 2.5 faster than our rapid rotator, it may be conceivable to reach the observed magnitudes.

Rotating stellar convection transports angular momentum towards the equator, generating the characteristic equatorial acceleration of the solar rotation while the radial flux of angular momentum is always inwards. New numerical box simulations for the meridional cross-correlation $\langle u_\theta u_\phi\rangle $, however, reveal the angular momentum transport towards the poles for slow rotation and towards the equator for fast rotation. The explanation is that for slow rotation a negative radial gradient of the angular velocity always appears, which in combination with a so-far neglected rotation-induced off-diagonal eddy viscosity term $\nu_\bot$ provides "antisolar rotation" laws with a decelerated equator. Similarly, the simulations provided positive values for the rotation-induced correlation $\langle u_r u_\theta\rangle $, which is relevant for the resulting latitudinal temperature profiles (cool or warm poles) for slow rotation and negative values for fast rotation. Observations of the differential rotation of slowly rotating stars will therefore lead to a better understanding of the actual stress-strain relation, the heat transport, and the underlying rotating convection.

Recent numerical work in helioseismology has shown that a periodically varying subsurface magnetic field leads to a fanning of the $f$-mode, which emerges from the density jump at the surface. In an attempt to model a more realistic situation, we now modulate this periodic variation with an envelope, giving thus more emphasis on localized bipolar magnetic structures in the middle of the domain. Some notable findings are: (i) compared to the purely hydrodynamic case, the strength of the $f$-mode is significantly larger at high horizontal wavenumbers $k$, but the fanning is weaker for the localized subsurface magnetic field concentrations investigated here than the periodic ones studied earlier; (ii) when the strength of the magnetic field is enhanced at a fixed depth below the surface, the fanning of the $f$-mode in the $k\omega$ diagram increases proportionally in such a way that the normalized $f$-mode strengths remain nearly the same in different such cases; (iii) the unstable Bloch modes reported previously in case of harmonically varying magnetic fields are now completely absent when more realistic localized magnetic field concentrations are imposed beneath the surface, thus suggesting that the Bloch modes are unlikely to be supported during most phases of the solar cycle; (iv) the $f$-mode strength appears to depend also on the depth of magnetic field concentrations such that it shows a relative decrement when the maximum of the magnetic field is moved to a deeper layer. We argue that detections of $f$-mode perturbations such as those being explored here could be effective tracers of solar magnetic fields below the photosphere before these are directly detectable as visible manifestations in terms of active regions or sunspots.

Context: Differential rotation in stars is driven by turbulent transport of angular momentum. Aims: To measure and parametrize the non-diffusive contribution to turbulent stress, known as $\Lambda$ effect, and its quenching as a function of rotation and magnetic field. Methods: Simulations of homogeneous, anisotropically forced turbulence in fully periodic cubes are used to extract the turbulent Reynolds and Maxwell stresses. Magnetic fields are introduced by imposing a uniform large-scale field on the system. Turbulent transport coefficients pertaining to the $\Lambda$ effect are obtained by fitting. Results: The numerical results agree qualitatively with analytic ones at slow rotation and low Reynolds numbers. This entails that vertical (horizontal) transport is downward (equatorward). The existence of a significant meridional $\Lambda$ effect is confirmed. Large-scale vorticity generation is found at rapid rotation when the Reynolds number exceeds a threshold value. The $\Lambda$ effect is severely quenched by large-scale magnetic fields due to the tendency of the Reynolds and Maxwell stresses to cancel each other. Rotational (magnetic) quenching of $\Lambda$ occurs at more rapid rotation (at lower field strength) in simulations than in theory. Conclusions: The current results largely confirm the earlier theoretical results but they also offer new insights: the non-negligible meridional $\Lambda$ effect possibly plays a role in the maintenance of meridional circulation in stars and the appearance of large-scale vortices raises the question of their effect on the angular momentum transport in stars. The results regarding magnetic quenching are consistent with the strong decrease of differential rotation in recent semi-global simulations and highlight the importance of including magnetic effects in differential rotation models. (abridged)

We consider the effect of a subadiabatic layer at the base of the convection zone on convection itself and the associated large-scale dynamos in spherical wedge geometry. We use a heat conduction prescription based on the Kramers opacity law which allows the depth of the convection zone to dynamically adapt to changes in the physical characteristics such as rotation rate and magnetic fields. We find that the convective heat transport is strongly concentrated toward the equatorial and polar regions in the cases without a substantial radiative layer below the convection zone. The presence of a stable layer below the convection zone significantly reduces the anisotropy of radial enthalpy transport. Furthermore, the dynamo solutions are sensitive to subtle changes in the convection zone structure. We find that the kinetic helicity changes sign in the deeper parts of the convection zone at high latitudes in all runs. This region expands progressively toward the equator in runs with a thicker stably stratified layer.

We test the sensitivity of hydrodynamic and magnetohydrodynamic turbulent convection simulations with respect to Mach number, thermal and magnetic boundary conditions, and the centrifugal force. We find that varying the luminosity, which also controls the Mach number, has only a minor effect on the large-scale dynamics. A similar conclusion can also be drawn from the comparison of two formulations of the lower magnetic boundary condition with either vanishing electric field or current density. The centrifugal force has an effect on the solutions, but only if its magnitude with respect to acceleration due to gravity is by two orders of magnitude greater than in the Sun. Finally, we find that the parameterisation of the photospheric physics, either by an explicit cooling term or enhanced radiative diffusion, is more important than the thermal boundary condition. In particular, runs with cooling tend to lead to more anisotropic convection and stronger deviations from the Taylor-Proudman state. In summary, the fully compressible approach taken here with the Pencil Code is found to be valid, while still allowing the disparate timescales to be taken into account.

We seek to understand the transition from nearly axisymmetric configurations at solar rotation rates to nonaxisymmetric configurations for rapid rotation using 3D numerical simulations of turbulent convection and considering rotation rates between 1 and 30 times the solar value. We find a transition from axi- to nonaxisymmetric solutions at around 1.8 times the solar rotation rate. This transition coincides with a change in the rotation profile from antisolar- to solar-like differential rotation with a faster equator and slow poles. In the solar-like rotation regime, the field configuration consists of an axisymmetric oscillatory field accompanied by an m=1 azimuthal mode (two active longitudes), which also shows temporal variability. At slow (rapid) rotation, the axisymmetric (nonaxisymmetric) mode dominates. The axisymmetric mode produces latitudinal dynamo waves with polarity reversals, while the nonaxisymmetric mode often exhibits a drift in the rotating reference frame and the strength of the active longitudes changes cyclically over time between the different hemispheres. Most of the obtained dynamo solutions exhibit cyclic variability either caused by latitudinal or azimuthal dynamo waves. In an activity-period diagram, the cycle lengths normalized by the rotation period form two different populations as a function of rotation rate or magnetic activity level. The slowly rotating axisymmetric population lies close to what is called the inactive branch in observations, while the rapidly rotating models are close to the superactive branch with a declining cycle to rotation frequency ratio with increasing rotation rate. We can successfully reproduce the transition from axi- to nonaxisymmetric dynamo solutions for high rotation rates, but high-resolution simulations are required to limit the effect of rotational quenching of convection at rotation rates above 20 times the solar value.

Using a recently developed two-scale formalism to determine the magnetic helicity spectrum (Brandenburg et al. 2017), we analyze synoptic vector magnetograms built with data from the Vector Spectromagnetograph (VSM) instrument on the \emph{Synoptic Optical Long-term Investigations of the Sun} (SOLIS) telescope during January 2010-July 2016. In contrast to an earlier study using only three Carrington rotations, our analysis includes 74 synoptic Carrington rotation maps. We recover here bihelical spectra at different phases of solar cycle~24, where the net magnetic helicity in the majority of the data is consistent with a large-scale dynamo with helical turbulence operating in the Sun. More than $20\%$ of the analyzed maps, however, show violations of the expected sign rule.

Small-scale dynamo action is often held responsible for the generation of quiet-Sun magnetic fields. We aim to determine the excitation conditions and saturation level of small-scale dynamos in non-rotating turbulent convection at low magnetic Prandtl numbers. We use high resolution direct numerical simulations of weakly stratified turbulent convection. We find that the critical magnetic Reynolds number for dynamo excitation increases as the magnetic Prandtl number is decreased, which might suggest that small-scale dynamo action is not automatically evident in bodies with small magnetic Prandtl numbers as the Sun. As a function of the magnetic Reynolds number (${\rm Rm}$), the growth rate of the dynamo is consistent with an ${\rm Rm}^{1/2}$ scaling. No evidence for a logarithmic increase of the growth rate with ${\rm Rm}$ is found.

We investigate dynamo action in global compressible solar-like convective dynamos in the framework of mean-field theory. We simulate a solar-type star in a wedge-shaped spherical shell, where the interplay between convection and rotation self-consistently drives a large-scale dynamo. To analyze the dynamo mechanism we apply the test-field method for azimuthally ($\phi$) averaged fields to determine the 27 turbulent transport coefficients of the electromotive force, of which six are related to the $\alpha$ tensor. This method has previously been used either in simulations in Cartesian coordinates or in the geodynamo context and is applied here for the first time to fully compressible simulations of solar-like dynamos. We find that the $\phi\phi$-component of the $\alpha$ tensor does not follow the profile expected from that of kinetic helicity. The turbulent pumping velocities significantly alter the effective mean flows acting on the magnetic field and therefore challenge the flux transport dynamo concept. All coefficients are significantly affected by dynamically important magnetic fields. Quenching as well as enhancement are being observed. This leads to a modulation of the coefficients with the activity cycle. The temporal variations are found to be comparable to the time-averaged values and seem to be responsible for a nonlinear feedback on the magnetic field generation. Furthermore, we quantify the validity of the Parker-Yoshimura rule for the equatorward propagation of the mean magnetic field in the present case.

Context: Convectively-driven flows play a crucial role in the dynamo processes that are responsible for producing magnetic activity in stars and planets. It is still not fully understood why many astrophysical magnetic fields have a significant large-scale component. Aims: Our aim is to investigate the dynamo properties of compressible convection in a rapidly rotating Cartesian domain, focusing upon a parameter regime in which the underlying hydrodynamic flow is known to be unstable to a large-scale vortex instability. Methods: The governing equations of three-dimensional nonlinear magnetohydrodynamics (MHD) are solved numerically. Different numerical schemes are compared and we propose a possible benchmark case for other similar codes. Results: In keeping with previous related studies, we find that convection in this parameter regime can drive a large-scale dynamo. The components of the mean horizontal magnetic field oscillate, leading to a continuous overall rotation of the mean field. Whilst the large-scale vortex instability dominates the early evolution of the system, it is suppressed by the magnetic field and makes a negligible contribution to the mean electromotive force that is responsible for driving the large-scale dynamo. The cycle period of the dynamo is comparable to the ohmic decay time, with longer cycles for dynamos in convective systems that are closer to onset. In these particular simulations, large-scale dynamo action is found only when vertical magnetic field boundary conditions are adopted at the upper and lower boundaries. Strongly modulated large-scale dynamos are found at higher Rayleigh numbers, with periods of reduced activity ("grand minima"-like events) occurring during transient phases in which the large-scale vortex temporarily re-establishes itself, before being suppressed again by the magnetic field.

We present numerical simulations of hydrodynamic overshooting convection in local Cartesian domains. We find that a substantial fraction of the lower part of the convection zone (CZ) is stably stratified according to the Schwarzschild criterion while the enthalpy flux is outward directed. This occurs when the heat conduction profile at the bottom of the CZ is smoothly varying, based either on a Kramers-like opacity prescription as a function of temperature and density or a static profile of a similar shape. We show that the subadiabatic layer arises due to nonlocal energy transport by buoyantly driven downflows in the upper parts of the CZ. Analysis of the force balance of the upflows and downflows confirms that convection is driven by cooling at the surface. We find that the commonly used prescription for the convective enthalpy flux being proportional to the negative entropy gradient does not hold in the stably stratified layers where the flux is positive. We demonstrate the existence of a non-gradient contribution to the enthalpy flux, which is estimated to be important throughout the convective layer. A quantitative analysis of downflows indicates a transition from a tree-like structure where smaller downdrafts merge into larger ones in the upper parts to a structure in the deeper parts where a height-independent number of strong downdrafts persist. This change of flow topology occurs when a substantial subadiabatic layer is present in the lower part of the CZ.

We focus on implementing and optimizing a sixth-order finite-difference solver for simulating compressible fluids on a GPU using third-order Runge-Kutta integration. Since graphics processing units perform well in data-parallel tasks, this makes them an attractive platform for fluid simulation. However, high-order stencil computation is memory-intensive with respect to both main memory and the caches of the GPU. We present two approaches for simulating compressible fluids using 55-point and 19-point stencils. We seek to reduce the requirements for memory bandwidth and cache size in our methods by using cache blocking and decomposing a latency-bound kernel into several bandwidth-bound kernels. Our fastest implementation is bandwidth-bound and integrates $343$ million grid points per second on a Tesla K40t GPU, achieving a $3.6 \times$ speedup over a comparable hydrodynamics solver benchmarked on two Intel Xeon E5-2690v3 processors. Our alternative GPU implementation is latency-bound and achieves the rate of $168$ million updates per second.

We explore the effects of an outer stably stratified coronal envelope on rotating turbulent convection, differential rotation, and large-scale dynamo action in spherical wedge models of the Sun. We solve the compressible magnetohydrodynamic equations in a two-layer model with unstable stratification below the surface, representing the convection zone, and a stably stratified coronal envelope above. The interface represents a free surface. We compare our model to models that have no coronal envelope. The presence of a coronal envelope is found to modify the Reynolds stress and the $\Lambda$ effect resulting in a weaker and non-cylindrical differential rotation. This is related to the reduced latitudinal temperature variations that are caused by and dependent on the angular velocity. Some simulations develop a near-surface shear layer that we can relate to a sign change in the meridional Reynolds stress term in the thermal wind balance equation. Furthermore, the presence of a free surface changes the magnetic field evolution since the toroidal field is concentrated closer to the surface. In all simulations, however, the migration direction of the mean magnetic field can be explained by the Parker--Yoshimura rule, which is consistent with earlier findings. A realistic treatment of the upper boundary in spherical dynamo simulations is crucial for the dynamics of the flow and magnetic field evolution.

(abidged) Context: Stellar convection zones are characterized by vigorous high-Reynolds number turbulence at low Prandtl numbers. Aims: We study the dynamo and differential rotation regimes at varying levels of viscous, thermal, and magnetic diffusion. Methods: We perform three-dimensional simulations of stratified fully compressible magnetohydrodynamic convection in rotating spherical wedges at various thermal and magnetic Prandtl numbers (from 0.25 to 2 and 5, respectively). Results: We find that the rotation profiles for high thermal diffusivity show a monotonically increasing angular velocity from the bottom of the convection zone to the top and from the poles toward the equator. For sufficiently rapid rotation, a region of negative radial shear develops at mid-latitudes as the thermal diffusivity is decreased. This coincides with a change in the dynamo mode from poleward propagating activity belts to equatorward propagating ones. Furthermore, the cyclic solutions disappear at the highest magnetic Reynolds numbers. The total magnetic energy increases with the magnetic Reynolds number in the range studied here ($5-151$), but the energies of the mean magnetic fields level off at high magnetic Reynolds numbers. The differential rotation is strongly affected by the magnetic fields and almost vanishes at the highest magnetic Reynolds numbers. In some of our most turbulent cases we find that two regimes are possible where either differential rotation is strong and mean magnetic fields relatively weak or vice versa. Conclusions: Our simulations indicate a strong non-linear feedback of magnetic fields on differential rotation, leading to qualitative changes in the behaviors of large-scale dynamos at high magnetic Reynolds numbers. Furthermore, we do not find indications of the simulations approaching an asymptotic regime where the results would be independent of diffusion coefficients.

Solar magnetic activity shows both smooth secular changes, such as the Grand Modern Maximum, and quite abrupt drops that are denoted as Grand Minima. Direct numerical simulations (DNS) of convection driven dynamos offer one way of examining the mechanisms behind these events. In this work, we analyze a solution of a solar-like DNS that has been evolved for roughly 80 magnetic cycles of 4.9 years, during which epochs of irregular behavior are detected. The emphasis of our analysis is to find physical causes for such behavior. The DNS employed is a semi-global (wedge) magnetoconvection model. For data analysis we use Ensemble Empirical Mode Decomposition (EEMD) and phase dispersion ($D^2$) methods. A special property of the DNS is the existence of multiple dynamo modes at different depths and latitudes. The dominant mode is solar-like. This mode is accompanied by a higher frequency mode near the surface and a low-frequency mode in the bottom of the convection zone. The overall behavior of the dynamo solution is very complex exhibiting variable cycle lengths, epochs of disturbed and even ceased surface activity, and strong short-term hemispherical asymmetries. Surprisingly, the most prominent suppressed surface activity epoch is actually a global magnetic energy maximum. We interpret the overall irregular behavior to be due to the interplay of the different dynamo modes showing different equatorial symmetries, especially the smoother part of the irregular variations being related to the variations of the mode strengths, evolving with different and variable cycle lengths. The abrupt low activity epoch in the dominant dynamo mode near the surface is related to a strong maximum of the bottom toroidal field strength, which causes abrupt disturbances especially in the differential rotation profile via the suppression of the Reynolds stresses.

We study the connection between spherical wedge and full spherical shell geometries using simple mean-field $\alpha^2$ dynamos. We solve the equations for a one-dimensional time-dependent mean-field dynamo to examine the effects of varying the polar angle $\theta_0$ between the latitudinal boundaries and the poles in spherical coordinates. We investigate the effects of turbulent magnetic diffusivity and $\alpha$ effect profiles as well as different latitudinal boundary conditions to isolate parameter regimes where oscillatory solutions are found. Finally, we add shear along with a damping term mimicking radial gradients to study the resulting dynamo regimes. We find that the commonly used perfect conductor boundary condition leads to oscillatory $\alpha^2$ dynamo solutions only if the wedge boundary is at least one degree away from the poles. Other boundary conditions always produce stationary solutions. By varying the profile of the turbulent magnetic diffusivity alone, oscillatory solutions are achieved with models extending to the poles, but the magnetic field is strongly concentrated near the poles and the oscillation period is very long. By introducing radial shear and a damping term mimicking radial gradients, we again see oscillatory dynamos, and the direction of drift follows the Parker--Yoshimura rule. Oscillatory solutions in the weak shear regime are found only in the wedge case with $\theta_0 = 1^\circ$ and perfect conductor boundaries. A reduced $\alpha$ effect near the poles with a turbulent diffusivity concentrated toward the equator yields oscillatory dynamos with equatorward migration and reproduces best the solutions in spherical wedges.

(abridged) Context: The mechanisms that cause the formation of sunspots are still unclear. Aims: We study the self-organisation of initially uniform sub-equipartition magnetic fields by highly stratified turbulent convection. Methods: We perform simulations of magnetoconvection in Cartesian domains that are $8.5$-$24$ Mm deep and $34$-$96$ Mm wide. We impose either a vertical or a horizontal uniform magnetic field in a convection-driven turbulent flow. Results: We find that super-equipartition magnetic flux concentrations are formed near the surface with domain depths of $12.5$ and $24$ Mm. The size of the concentrations increases as the box size increases and the largest structures ($20$ Mm horizontally) are obtained in the 24 Mm deep models. The field strength in the concentrations is in the range of $3$-$5$ kG. The concentrations grow approximately linearly in time. The effective magnetic pressure measured in the simulations is positive near the surface and negative in the bulk of the convection zone. Its derivative with respect to the mean magnetic field, however, is positive in the majority of the domain, which is unfavourable for the negative effective magnetic pressure instability (NEMPI). Furthermore, we find that magnetic flux is concentrated in regions of converging flow corresponding to large-scale supergranulation convection pattern. Conclusions: The linear growth of large-scale flux concentrations implies that their dominant formation process is tangling of the large-scale field rather than an instability. One plausible mechanism explaining both the linear growth and the concentrate on of the flux in the regions of converging flow pattern is flux expulsion. Possible reasons for the absence of NEMPI are that the derivative of the effective magnetic pressure with respect to the mean magnetic field has an unfavourable sign and that there may not be sufficient scale separation.

In light of new results, the one-dimensional mean-field dynamo model of Brandenburg & Kapyla (2007) with dynamical quenching and a nonlocal Babcock-Leighton alpha effect is re-examined for the solar dynamo. We extend the one-dimensional model to include the effects of turbulent downward pumping (Kitchatinov & Olemskoy 2011), and to combine dynamical quenching with shear. We use both the conventional dynamical quenching model of Kleeorin & Ruzmaikin (1982) and the alternate one of Hubbard & Brandenburg (2011), and confirm that with varying levels of non-locality in the alpha effect, and possibly shear as well, the saturation field strength can be independent of the magnetic Reynolds number.

We present results from four convectively-driven stellar dynamo simulations in spherical wedge geometry. All of these simulations produce cyclic and migrating mean magnetic fields. Through detailed comparisons we show that the migration direction can be explained by an $\alpha\Omega$ dynamo wave following the Parker--Yoshimura rule. We conclude that the equatorward migration in this and previous work is due to a positive (negative) $\alpha$ effect in the northern (southern) hemisphere and a negative radial gradient of $\Omega$ outside the inner tangent cylinder of these models. This idea is supported by a strong correlation between negative radial shear and toroidal field strength in the region of equatorward propagation.

(abridged) Context: Solar-like differential rotation is characterized by a rapidly rotating equator and slower poles. However, theoretical models and numerical simulations can result in a slower equator and faster poles when the rotation is slow. Aims: We study the critical rotational influence under which differential rotation flips from solar-like (fast equator, slow poles) to an anti-solar one (slow equator, fast poles). We estimate the non-diffusive ($\Lambda$ effect) and diffusive (turbulent viscosity) contributions to the Reynolds stress. Methods: We present the results of three-dimensional numerical simulations of mildly turbulent convection in spherical wedge geometry. Here we apply a fully compressible setup which would suffer from a prohibitive time step constraint if the real solar luminosity was used. We regulate the convective velocities by varying the amount of heat transported by thermal conduction, turbulent diffusion, and resolved convection. Results: Increasing the efficiency of resolved convection leads to a reduction of the rotational influence on the flow and a sharp transition from solar-like to anti-solar differential rotation for Coriolis numbers around 1.3. We confirm the recent finding of a large-scale flow bistability: contrasted with running the models from an initial condition with unprescribed differential rotation, the initialization of the model with certain kind of rotation profile sustains the solution over a wider parameter range. Conclusions: Our results may have implications for real stars that start their lives as rapid rotators implying solar-like rotation in the early main-sequence evolution. As they slow down, they might be able to retain solar-like rotation for lower Coriolis numbers before switching to anti-solar rotation. This could partially explain the puzzling findings of anti-solar rotation profiles for models in the solar parameter regime.

We compare simple analytical closure models of homogeneous turbulent Boussinesq convection for stellar applications with three-dimensional simulations. We use simple analytical closure models to compute the fluxes of angular momentum and heat as a function of rotation rate measured by the Taylor number. We also investigate cases with varying angles between the angular velocity and gravity vectors, corresponding to locating the computational domain at different latitudes ranging from the pole to the equator of the star. We perform three-dimensional numerical simulations in the same parameter regimes for comparison. The free parameters appearing in the closure models are calibrated by two fitting methods using simulation data. Unique determination of the closure parameters is possible only in the non-rotating case or when the system is placed at the pole. In the other cases the fit procedures yield somewhat differing results. The quality of the closure is tested by substituting the resulting coefficients back into the closure model and comparing with the simulation results. To eliminate the possibilities that the results obtained depend on the aspect ratio of the simulation domain or suffer from too small Rayleigh numbers we performed runs varying these parameters. The simulation data for the Reynolds stress and heat fluxes broadly agree with previous compressible simulations. The closure works fairly well with slow and fast rotation but its quality degrades for intermediate rotation rates. We find that the closure parameters depend not only on rotation rate but also on latitude. The weak dependence on Rayleigh number and the aspect ratio of the domain indicates that our results are generally valid

We study the generation of large-scale vortices in rotating turbulent convection by means of Cartesian direct numerical simulations. We find that for sufficiently rapid rotation, cyclonic structures on a scale large in comparison to that of the convective eddies, emerge, provided that the fluid Reynolds number exceeds a critical value. For slower rotation, cool cyclonic vortices are preferred, whereas for rapid rotation, warm anti-cyclonic vortices are favored. In some runs in the intermediate regime both types of cyclones co-exist for thousands of convective turnover times. The temperature contrast between the vortices and the surrounding atmosphere is of the order of five per cent. We relate the simulation results to observations of rapidly rotating late-type stars that are known to exhibit large high-latitude spots from Doppler imaging. In many cases, cool spots are accompanied with spotted regions with temperatures higher than the average. In this paper, we investigate a scenario according to which the spots observed in the temperature maps could have a non-magnetic origin due to large-scale vortices in the convection zones of the stars.

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.

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 3-7 times the solar value. The phase relation of the magnetic field is such that the toroidal field lags behind the radial field by about $\pi/2$, which can be explained by an oscillatory $\alpha^2$ dynamo caused by the sign change of the $\alpha$-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\pi$ azimuthal extent, reminiscent of the field configurations deduced from observations of rapidly rotating late-type stars.

We present results from compressible Cartesian convection simulations with and without imposed shear. In the former case the dynamo is expected to be of $\alpha^2\varOmega$ type which is generally expected to be relevant for the Sun, whereas the latter case refers to $\alpha^2$ dynamos which are more likely to occur in more rapidly rotating stars whose differential rotation is small. We perform a parameter study where the shear flow and the rotational influence are varied to probe the relative importance of both types of dynamos. Oscillatory solutions are preferred both in the kinematic and saturated regimes when the negative ratio of shear to rotation rates, $q\equiv -S/\varOmega$, is between 1.5 and 2, i.e., when shear and rotation are of comparable strengths. Other regions of oscillatory solutions are found with small values of $q$, i.e., when shear is weak in comparison to rotation, and in the regime of large negative $q$s, when shear is very strong in comparison to rotation. However, exceptions to these rules also appear so that for a given ratio of shear to rotation, solutions are non-oscillatory for small and large shear, but oscillatory in the intermediate range. Changing the boundary conditions from vertical field to perfect conductor ones changes the dynamo mode from oscillatory to quasi-steady. Furthermore, in many cases an oscillatory solution exists only in the kinematic regime whereas in the nonlinear stage the mean fields are stationary. However, the cases with rotation and no shear are always oscillatory in the parameter range studied here and the dynamo mode does not depend on the magnetic boundary conditions. The strengths of total and large-scale components of the magnetic field in the saturated state, however, are sensitive to the chosen boundary conditions.

We report on the results of four convective dynamo simulations with an outer coronal layer. The magnetic field is self-consistently generated by the convective motions beneath the surface. Above the convection zone, we include a polytropic layer that extends to 1.6 solar radii. The temperature increases in this region to $\approx8$ times the value at the surface, corresponding to $\approx1.2$ times the value at the bottom of the spherical shell. We associate this region with the solar corona. We find solar-like differential rotation with radial contours of constant rotation rate, together with a near-surface shear layer. This non-cylindrical rotation profile is caused by a non-zero latitudinal entropy gradient that offsets the Taylor--Proudman balance through the baroclinic term. The meridional circulation is multi-cellular with a solar-like poleward flow near the surface at low latitudes. In most of the cases, the mean magnetic field is oscillatory with equatorward migration in two cases. In other cases, the equatorward migration is overlaid by stationary or even poleward migrating mean fields.

In astrophysics, turbulent diffusion is often used in place of microphysical diffusion to avoid resolving the small scales. However, we expect this approach to break down when time and length scales of the turbulence become comparable with other relevant time and length scales in the system. Turbulent diffusion has previously been applied to the magneto-rotational instability (MRI), but no quantitative comparison of growth rates at different turbulent intensities has been performed. We investigate to what extent turbulent diffusion can be used to model the effects of small-scale turbulence on the kinematic growth rates of the MRI, and how this depends on angular velocity and magnetic field strength. We use direct numerical simulations in three-dimensional shearing boxes with periodic boundary conditions in the spanwise direction and additional random plane-wave volume forcing to drive a turbulent flow at a given length scale. We estimate the turbulent diffusivity using a mixing length formula and compare with results obtained with the test-field method. It turns out that the concept of turbulent diffusion is remarkably accurate in describing the effect of turbulence on the growth rate of the MRI. No noticeable breakdown of turbulent diffusion has been found, even when time and length scales of the turbulence become comparable with those imposed by the MRI itself. On the other hand, quenching of turbulent magnetic diffusivity by the magnetic field is found to be absent. Turbulence reduces the growth rate of the MRI in a way that is the same as microphysical magnetic diffusion.

We report the finding of an azimuthal dynamo wave of a low-order (m=1) mode in direct numerical simulations (DNS) of turbulent convection in spherical shells. Such waves are predicted by mean field dynamo theory and have been obtained previously in mean-field models. Observational results both from photometry and Doppler imaging have revealed persistent drifts of spots for several rapidly rotating stars, but, although an azimuthal dynamo wave has been proposed as a possible mechanism responsible for this behavior, it has been judged as unlikely, as practical evidence for such waves from DNS has been lacking. The large-scale magnetic field in our DNS, which is due to self-consistent dynamo action, is dominated by a retrograde m=1 mode. Its pattern speed is nearly independent of latitude and does not reflect the speed of the differential rotation at any depth. The extrema of magnetic m=1 structures coincide reasonably with the maxima of m=2 structures of the temperature. These results provide direct support for the observed drifts being due to an azimuthal dynamo wave.

Using simulations of slowly rotating stratified turbulence, we show that the alpha effect responsible for the generation of astrophysical magnetic fields is proportional to the logarithmic gradient of kinetic energy density rather than that of momentum, as was previously thought. This result is in agreement with a new analytic theory developed in this paper for large Reynolds numbers. Thus, the contribution of density stratification is less important than that of turbulent velocity. The alpha effect and other turbulent transport coefficients are determined by means of the test-field method. In addition to forced turbulence, we also investigate supernova-driven turbulence and stellar convection. In some cases (intermediate rotation rate for forced turbulence, convection with intermediate temperature stratification, and supernova-driven turbulence) we find that the contribution of density stratification might be even less important than suggested by the analytic theory.

In this paper we first discuss observational evidence of longitudinal concentrations of magnetic activity in the Sun and rapidly rotating late-type stars with outer convective envelopes. Scenarios arising from the idea of rotationally influenced anisotropic convective turbulence being the key physical process generating these structures are then presented and discussed - such effects include the turbulent dynamo mechanism, negative effective magnetic pressure instability (NEMPI) and hydrodynamical vortex instability. Finally, we discuss non-axisymmetric stellar mean-field dynamo models, the results obtained with them, and compare those with the observational information gathered up so far. We also present results from a pure alpha-squared mean-field dynamo model, which show that time-dependent behavior of the dynamo solutions can occur both in the form of an azimuthal dynamo wave and/or oscillatory behavior related to the alternating energy levels of the active longitudes.

We present preliminary results from high resolution magneto-convection simulations where we find the formation of flux concentrations from an initially uniform magnetic field. We compute the effective magnetic pressure but find that the concentrations appear also in places where it is positive. The structures appear in roughly ten convective turnover times and live close to a turbulent diffusion time. The time scales are compatible with the negative effective magnetic pressure instability (NEMPI), although structure formation is not restricted to regions where the effective magnetic pressure is negative.

We present results of convective turbulent dynamo simulations including a coronal layer in a spherical wedge. We find an equatorward migration of the radial and azimuthal fields similar to the behavior of sunspots during the solar cycle. The migration of the field coexist with a spoke-like differential rotation and anti-solar (clockwise) meridional circulation. Even though the migration extends over the whole convection zone, the mechanism causing this is not yet fully understood.

We combine a convectively driven dynamo in a spherical shell with a nearly isothermal density-stratified cooling layer that mimics some aspects of a stellar corona to study the emergence and ejections of magnetic field structures. This approach is an extension of earlier models, where forced turbulence simulations were employed to generate magnetic fields. A spherical wedge is used which consists of a convection zone and an extended coronal region to $\approx1.5$ times the radius of the sphere. The wedge contains a quarter of the azimuthal extent of the sphere and $150\degr$ in latitude. The magnetic field is self-consistently generated by the turbulent motions due to convection beneath the surface. Magnetic fields are found to emerge at the surface and are ejected to the coronal part of the domain. These ejections occur at irregular intervals and are weaker than in earlier work. We tentatively associate these events with coronal mass ejections on the Sun, even though our model of the solar atmosphere is rather simplistic.

We investigate the effects of weakly and strongly stratified turbulent convection on the mean effective Lorentz force, and especially on the mean effective magnetic pressure. Earlier studies with isotropically forced non-stratified and stratified turbulence have shown that the contribution of the turbulence to the mean magnetic pressure is negative for mean horizontal magnetic fields that are smaller than the equipartition strength, so that the effective mean magnetic pressure that takes into account the turbulence effects, can be negative. Compared with earlier cases of forced turbulence with an isothermal equation of state, we find that the turbulence effect is similar to or even stronger in the present case of turbulent convection. This is argued to be due to the anisotropy of turbulence in the vertical direction. Another important difference compared with earlier studies is the presence of an evolution equation for the specific entropy. Mean-field modelling with entropy evolution indicates that the negative effective magnetic pressure can still lead to a large-scale instability which forms local flux concentrations, even though the specific entropy evolution tends to have a stabilizing effect when applied to a stably stratified (e.g., isothermal) layer. It is argued that this large-scale instability could be important for the formation of solar large-scale magnetic structures such as active regions.

We report on simulations of mildly turbulent convection in spherical wedge geometry with varying density stratification. We vary the density contrast within the convection zone by a factor of 20 and study the influence of rotation on the solutions. We demonstrate that the size of convective cells decreases and the anisotropy of turbulence increases as the stratification is increased. Differential rotation is found to change from anti-solar (slow equator) to solar-like (fast equator) at roughly the same Coriolis number for all stratifications. The largest stratification runs, however, are sensitive to changes of the Reynolds number. Evidence for a near-surface shear layer is found in runs with strong stratification and large Reynolds numbers.

Direct numerical simulations of isotropically forced homogeneous stationary turbulence with an imposed passive scalar concentration gradient are compared with an analytical closure model which provides evolution equations for the mean passive scalar flux and variance. Triple correlations of fluctuations appearing in these equations are described in terms of relaxation terms proportional to the quadratic correlations. Three methods are used to extract the relaxation timescales tau_i from direct numerical simulations. Firstly, we insert the closure ansatz into our equations, assume stationarity, and solve for tau_i. Secondly, we use only the closure ansatz itself and obtain tau_i from the ratio of quadratic and triple correlations. Thirdly we remove the imposed passive scalar gradient and fit an exponential decay law to the solution. We vary the Reynolds (Re) and P\'eclet (Pe) numbers while keeping their ratio at unity and the degree of scale separation and find for large Re fair correspondence between the different methods. The ratio of the turbulent relaxation time of passive scalar flux to the turnover time of turbulent eddies is of the order of three, which is in remarkable agreement with earlier work. Finally we make an effort to extract the relaxation timescales relevant for the viscous and diffusive effects. We find two regimes which are valid for small and large Re, respectively, but the dependence of the parameters on scale separation suggests that they are not universal.

We present a three-dimensional model of rotating convection combined with a simplified corona in spherical coordinates. The motions in the convection zone generate a large-scale magnetic field that gets sporadically ejected into the outer layers above. Our model corona is approximately isothermal, but it includes density stratification due to gravity.

We determine the timescales associated with turbulent diffusion and isotropization in closure models using anisotropically forced and freely decaying turbulence simulations and to study the applicability of these models. We compare the results from anisotropically forced three-dimensional numerical simulations with the predictions of the closure models and obtain the turbulent timescales mentioned above as functions of the Reynolds number. In a second set of simulations, turning the forcing off enables us to study the validity of the closures in freely decaying turbulence. Both types of experiments suggest that the timescale of turbulent diffusion converges to a constant value at higher Reynolds numbers. Furthermore, the relative importance of isotropization is found to be about 2.5 times larger at higher Reynolds numbers than in the more viscous regime.

In an earlier study, we reported on the excitation of large-scale vortices in Cartesian hydrodynamical convection models subject to rapid enough rotation. In that study, the conditions of the onset of the instability were investigated in terms of the Reynolds (Re) and Coriolis (Co) numbers in models located at the stellar North pole. In this study, we extend our investigation to varying domain sizes, increasing stratification and place the box at different latitudes. The effect of the increasing box size is to increase the sizes of the generated structures, so that the principal vortex always fills roughly half of the computational domain. The instability becomes stronger in the sense that the temperature anomaly and change in the radial velocity are observed to be enhanced. The model with the smallest box size is found to be stable against the instability, suggesting that a sufficient scale separation between the convective eddies and the scale of the domain is required for the instability to work. The instability can be seen upto the co-latitude of 30 degrees, above which value the flow becomes dominated by other types of mean flows. The instability can also be seen in a model with larger stratification. Unlike the weakly stratified cases, the temperature anomaly caused by the vortex structures is seen to depend on depth.

Using different analytical methods (the quasi-linear approach, the path-integral technique and tau-relaxation approximation) we develop a comprehensive mean-field theory for a pumping effect of the mean magnetic field in homogeneous non-rotating helical turbulence with imposed large-scale shear. The effective pumping velocity is proportional to the product of alpha effect and large-scale vorticity associated with the shear, and causes a separation of the toroidal and poloidal components of the mean magnetic field along the direction of the mean vorticity. We also perform direct numerical simulations of sheared turbulence in different ranges of hydrodynamic and magnetic Reynolds numbers and use a kinematic test-field method to determine the effective pumping velocity. The results of the numerical simulations are in agreement with the theoretical predictions.

Numerical simulations of the magnetorotational instability (MRI) with zero initial net flux in a non-stratified isothermal cubic domain are used to demonstrate the importance of magnetic boundary conditions.In fully periodic systems the level of turbulence generated by the MRI strongly decreases as the magnetic Prandtl number (Pm), which is the ratio of kinematic viscosity and magnetic diffusion, is decreased. No MRI or dynamo action below Pm=1 is found, agreeing with earlier investigations. Using vertical field conditions, which allow the generation of a net toroidal flux and magnetic helicity fluxes out of the system, the MRI is found to be excited in the range 0.1 < Pm < 10, and that the saturation level is independent of Pm. In the vertical field runs strong mean-field dynamo develops and helps to sustain the MRI.

Global dynamo simulations solving the equations of magnetohydrodynamics (MHD) have been a tool of astrophysicists who try to understand the magnetism of the Sun for several decades now. During recent years many fundamental issues in dynamo theory have been studied in detail by means of local numerical simulations that simplify the problem and allow the study of physical effects in isolation. Global simulations, however, continue to suffer from the age-old problem of too low spatial resolution, leading to much lower Reynolds numbers and scale separation than in the Sun. Reproducing the internal rotation of the Sun, which plays a crucual role in the dynamo process, has also turned out to be a very difficult problem. In the present paper the current status of global dynamo simulations of the Sun is reviewed. Emphasis is put on efforts to understand how the large-scale magnetic fields, i.e. whose length scale is greater than the scale of turbulence, are generated in the Sun. Some lessons from mean-field theory and local simulations are reviewed and their possible implications to the global models are discussed. Possible remedies to some current issues of solar simulations are put forward.

Context. Earlier work has suggested that large-scale dynamos can reach and maintain equipartition field strengths on a dynamical time scale only if magnetic helicity of the fluctuating field can be shed from the domain through open boundaries. Aims. Our aim is to test this scenario in convection-driven dynamos by comparing results for open and closed boundary conditions. Methods. Three-dimensional numerical simulations of turbulent compressible convection with shear and rotation are used to study the effects of boundary conditions on the excitation and saturation of large-scale dynamos. Open (vertical-field) and closed (perfect-conductor) boundary conditions are used for the magnetic field. The shear flow is such that the contours of shear are vertical, crossing the outer surface, and are thus ideally suited for driving a shear-induced magnetic helicity flux. Results. We find that for given shear and rotation rate, the growth rate of the magnetic field is larger if open boundary conditions are used. The growth rate first increases for small magnetic Reynolds number, Rm, but then levels off at an approximately constant value for intermediate values of Rm. For large enough Rm, a small-scale dynamo is excited and the growth rate of the field in this regime increases as Rm^(1/2). Regarding the nonlinear regime, the saturation level of the energy of the total magnetic field is independent of Rm when open boundaries are used. In the case of perfect conductor boundaries, the saturation level first increases as a function of Rm, but then decreases proportional to Rm^(-1) for Rm > 30, indicative of catastrophic quenching. These results suggest that the shear-induced magnetic helicity flux is efficient in alleviating catastrophic quenching when open boundaries are used. The horizontally averaged mean field is still weakly decreasing as a function of Rm even for open boundaries.

(abridged) Context. Turbulent fluxes of angular momentum and heat due to rotationally affected convection play a key role in determining differential rotation of stars. Here we perform a systematic comparison between Cartesian and spherical geometries as a function of the rotation rate. Aims. We 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 particular, we want to clarify whether the sharp equatorial profile of the horizontal Reynolds stress found in earlier Cartesian models is reproduced in spherical models. Methods. We employ direct numerical simulations of turbulent convection. In order to reach as high spatial resolution as possible in the spherical runs, we model only parts of the latitude and longitude. The rotational influence, measured by the Coriolis 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 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'. The latitudinal heat flux is mostly equatorward for slow rotation but changes sign for rapid rotation. 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.

Angular momentum transport owing to hydrodynamic turbulent convection is studied using local three dimensional numerical simulations employing the shearing box approximation. We determine the turbulent viscosity from non-rotating runs over a range of values of the shear parameter and use a simple analytical model in order to extract the non-diffusive contribution (Lambda-effect) to the stress in runs where rotation is included. Our results suggest that the turbulent viscosity is of the order of the mixing length estimate and weakly affected by rotation. The Lambda-effect is non-zero and a factor of 2-4 smaller than the turbulent viscosity in the slow rotation regime. We demonstrate that for Keplerian shear, the angular momentum transport can change sign and be outward when the rotation period is greater than the turnover time, i.e. when the Coriolis number is below unity. This result seems to be relatively independent of the value of the Rayleigh number.

Convectively driven dynamos with rotation generating magnetic fields on scales large compared with the scale of the turbulent eddies are being reviewed. It is argued that such fields can be understood as the result of an alpha effect. Simulations in Cartesian domains show that such large-scale magnetic fields saturate on a time scale compatible with the resistive one, suggesting that the magnitude of the alpha effect is here still constrained by approximate magnetic helicity conservation. It is argued that, in the absence of shear and/or any other known large-scale dynamo effects, these simulations prove the existence of turbulent alpha^2-type dynamos. Finally, recent results are discussed in the context of solar and stellar dynamos.

A hypothesis for sunspot formation is the buoyant emergence of magnetic flux tubes created by the strong radial shear at the tachocline. In this scenario, the magnetic field has to exceed a threshold value before it becomes buoyant and emerges through the whole convection zone. We follow the evolution of a random seed magnetic field with the aim of study under what conditions it is possible to excite the dynamo instability and whether the dynamo generated magnetic field becomes buoyantly unstable and emerges to the surface as expected in the flux-tube context. We perform numerical simulations of compressible turbulent convection that include a vertical shear layer. Like the solar tachocline, the shear is located at the interface between convective and stable layers. We find that shear and convection are able to amplify the initial magnetic field and form large-scale elongated magnetic structures. The magnetic field strength depends on several parameters such as the shear amplitude, the thickness and location of the shear layer, and the magnetic Reynolds number ($\Rm$). Whenever the toroidal magnetic field reaches amplitudes greater a threshold value which is close to the equipartition value, it becomes buoyant and rises into the convection zone where it expands and forms mushroom shape structures. Some events of emergence, i.e. those with the largest amplitudes of the initial field, are able to reach the very uppermost layers of the domain. These episodes are able to modify the convective pattern forming either broader convection cells or convective eddies elongated in the direction of the field. However, in none of these events the field preserves its initial structure.

Hot luminous stars show a variety of phenomena in their photospheres and in their winds which still lack clear physical explanations at this time. Among these phenomena are non-thermal line broadening, line profile variability (LPVs), discrete absorption components (DACs), wind clumping and stochastically excited pulsations. Cantiello et al. (2009) argued that a convection zone close to the surface of hot, massive stars, could be responsible for some of these phenomena. This convective zone is caused by a peak in the opacity due to iron recombination and for this reason is referred as the "iron convection zone" (FeCZ). 3D MHD simulations are used to explore the possible effects of such subsurface convection on the surface properties of hot, massive stars. We argue that turbulence and localized magnetic spots at the surface are the likely consequence of subsurface convection in early type stars.

During their main sequence evolution, massive stars can develop convective regions very close to their surface. These regions are caused by an opacity peak associated with iron ionization. Cantiello et al. (2009) found a possible connection between the presence of sub-photospheric convective motions and small scale stochastic velocities in the photosphere of early-type stars. This supports a physical mechanism where microturbulence is caused by waves that are triggered by subsurface convection zones. They further suggest that clumping in the inner parts of the winds of OB stars could be related to subsurface convection, and that the convective layers may also be responsible for stochastic excitation of non-radial pulsations. Furthermore, magnetic fields produced in the iron convection zone could appear at the surface of such massive stars. Therefore subsurface convection could be responsible for the occurrence of observable phenomena such as line profile variability and discrete absorption components. These phenomena have been observed for decades, but still evade a clear theoretical explanation. Here we present preliminary results from 3D MHD simulations of such subsurface convection.

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 time scales 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.

Self-consistent convective dynamo simulations in wedge-shaped spherical shells are presented. Differential rotation is generated by the interaction of convection with rotation. Equatorward acceleration and dynamo action are obtained only for sufficiently rapid rotation. The angular velocity tends to be constant along cylinders. Oscillatory large-scale fields are found to migrate in the poleward direction. Comparison with earlier simulations in full spherical shells and Cartesian domains is made.

Using three-dimensional convection simulations it is shown that a sinusoidal variation of horizontal shear leads to a kinematic \alpha effect with a similar sinusoidal variation. The effect exists even for weak stratification and arises owing to the inhomogeneity of turbulence and the presence of impenetrable vertical boundaries. This system produces large-scale magnetic fields that also show a sinusoidal variation in the streamwise direction. It is argued that earlier investigations overlooked these phenomena partly because of the use of horizontal averaging and also because measurements of \alpha using an imposed field combined with long time averages give erroneous results. It is demonstrated that in such cases the actual horizontally averaged mean field becomes non-uniform. The turbulent magnetic diffusion term resulting from such non-uniform fields can then no longer be neglected and begins to balance the \alpha effect.

The predictive power of mean-field theory is emphasized by comparing theory with simulations under controlled conditions. The recently developed test-field method is used to extract turbulent transport coefficients both in kinematic as well as nonlinear and quasi-kinematic cases. A striking example of the quasi-kinematic method is provided by magnetic buoyancy-driven flows that produce an alpha effect and turbulent diffusion.

In this paper we make an effort to understand the interaction of turbulence generated by the magnetorotational instability (MRI) with turbulence from other sources, such as supernova explosions (SNe) in galactic disks. First we perform a linear stability analysis (LSA) of non-ideal MRI to derive the limiting value of Ohmic diffusion that is needed to inhibit the growth of the instability for different types of rotation laws. With the help of a simple analytical expression derived under first-order smoothing approximation (FOSA), an estimate of the limiting turbulence level and hence the turbulent diffusion needed to damp the MRI is derived. Secondly, we perform numerical simulations in local cubes of isothermal nonstratified gas with external forcing of varying strength to see whether the linear result holds for more complex systems. Purely hydrodynamic calculations with forcing, rotation and shear are made for reference purposes, and as expected, non-zero Reynolds stresses are found. In the magnetohydrodynamic calculations, therefore, the total stresses generated are a sum of the forcing and MRI contributions. To separate these contributions, we perform reference runs with MRI-stable shear profiles (angular velocity increasing outwards), which suggest that the MRI-generated stresses indeed become strongly suppressed as function of the forcing. The Maxwell to Reynolds stress ratio is observed to decrease by an order of magnitude as the turbulence level due to external forcing exceeds the predicted limiting value, which we interpret as a sign of MRI suppression. Finally, we apply these results to estimate the limiting radius inside of which the SN activity can suppress the MRI, arriving at a value of 14 kpc.

To study the Reynolds stresses which describe turbulent momentum transport from turbulence affected by large-scale shear and rotation. Three-dimensional numerical simulations are used to study turbulent transport under the influences of large-scale shear and rotation in homogeneous, isotropically forced turbulence. We study three cases: one with only shear, and two others where in addition to shear, rotation is present. These cases differ by the angle (0 or 90\degr) the rotation vector makes with respect to the z-direction. Two subsets of runs are performed with both values of \theta where either rotation or shear is kept constant. When only shear is present, the off-diagonal stress can be described by turbulent viscosity whereas if the system also rotates, nondiffusive contributions (\Lambda-effect) to the stress can arise. Comparison of the direct simulations are made with analytical results from a simple closure model. We find that the turbulent viscosity is of the order of the first order smoothing result in the parameter regime studied and that for sufficiently large Reynolds numbers the Strouhal number, describing the ratio of correlation to turnover times, is roughly 1.5. This is consistent with the closure model based on the minimal tau-approximation which produces a reasonable fit to the simulation data for similar Strouhal numbers. In the cases where rotation is present, separating the diffusive and nondiffusive components of the stress turns out to be challenging but taking the results at face value, we can obtain nondiffusive contributions of the order of 0.1 times the turbulent viscosity. We also find that the simple closure model is able to reproduce most of the qualitative features of the numerical results provided that the Strouhal number is of the order of unity.

Estimates for the nonlinear alpha effect in helical turbulence with an applied magnetic field are presented using two different approaches: the imposed-field method where the electromotive force owing to the applied field is used, and the test-field method where separate evolution equations are solved for a set of different test fields. Both approaches agree for stronger fields, but there are apparent discrepancies for weaker fields that can be explained by the influence of dynamo-generated magnetic fields on the scale of the domain that are referred to as meso-scale magnetic fields. Examples are discussed where these meso-scale fields can lead to both drastically overestimated and underestimated values of alpha compared with the kinematic case. It is demonstrated that the kinematic value can be recovered by resetting the fluctuating magnetic field to zero in regular time intervals. It is concluded that this is the preferred technique both for the imposed-field and the test-field methods.

The existence of large-scale dynamos in rigidly rotating turbulent convection without shear is studied using three-dimensional numerical simulations of penetrative rotating compressible convection. We demonstrate that rotating convection in a Cartesian domain can drive a large-scale dynamo even in the absence of shear. The large-scale field contains a significant fraction of the total field in the saturated state. The simulation results are compared with one-dimensional mean-field dynamo models where turbulent transport coefficients, as determined using the test field method, are used. The reason for the absence of large-scale dynamo action in earlier studies is shown to be due to the rotation being too slow: whereas the alpha-effect can change sign, its magnitude stays approximately constant as a function of rotation, and the turbulent diffusivity decreases monotonically with increasing rotation. Only when rotation is rapid enough a large-scale dynamo can be excited. The one-dimensional mean-field model with dynamo coefficients from the test field results predicts reasonably well the dynamo excitation in the direct simulations. This result further validates the test field procedure and reinforces the interpretation that the observed dynamo is driven by a turbulent alpha-effect. This result demonstrates the existence of an alpha-effect and an alpha^2-dynamo with natural forcing.

(abridged) Aims: To 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 3D 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. 1D mean-field models are used with the derived turbulent transport coefficients to compare with direct simulations. Results: The alpha-effect increases monotonically as rotation increases. Turbulent diffusivity, eta_t, is proportional to the square of the turbulent vertical velocity. Whereas eta_t decreases approximately inversely proportional to the wavenumber of the field, the alpha-effect and turbulent pumping show a more complex behaviour. In the presence of shear and no rotation a small alpha-effect is induced which does not seem to show any consistent trend as a function of shear. If the shear is large enough, this small alpha 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 well below those in direct simulations suggesting that an incoherent alpha-shear dynamo may also act in them. If both rotation and shear are present, the alpha-effect is more pronounced. The combination of the shear-current and Omega x J-effects is also stronger than in the case of shear only, but subdominant to the alpha-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.

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 modelled in terms of an \alpha\Omega 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 \alpha-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 (Rm) 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 Rm.

Simulations of stochastically forced shear-flow turbulence in a shearing-periodic domain are used to study the spontaneous generation of large-scale flow patterns in the direction perpendicular to the plane of the shear. Based on an analysis of the resulting large-scale velocity correlations it is argued that the mechanism behind this phenomenon could be the mean-vorticity dynamo effect pioneered by Elperin, Kleeorin, and Rogachevskii in 2003 (Phys. Rev. E 68, 016311). This effect is based on the anisotropy of the eddy viscosity tensor. One of its components may be able to replenish cross-stream mean flows by acting upon the streamwise component of the mean flow. Shear, in turn, closes the loop by acting upon the cross-stream mean flow to produce stronger streamwise mean flows. The diagonal component of the eddy viscosity is found to be of the order of the rms turbulent velocity divided by the wavenumber of the energy-carrying eddies.

We study the dependence of turbulent transport coefficients, such as the components of the $\alpha$ tensor ($\alpha_{ij}$) and the turbulent magnetic diffusivity tensor ($\eta_{ij}$), on shear and magnetic Reynolds number in the presence of helical forcing. We use three-dimensional direct numerical simulations with periodic boundary conditions and measure the turbulent transport coefficients using the kinematic test field method. In all cases the magnetic Prandtl number is taken as unity. We find that with increasing shear the diagonal components of $\alpha_{ij}$ quench, whereas those of $\eta_{ij}$ increase. The antisymmetric parts of both tensors increase with increasing shear. We also propose a simple expression for the turbulent pumping velocity (or $\gamma$ effect). This pumping velocity is proportional to the kinetic helicity of the turbulence and the vorticity of the mean flow. For negative helicity, i.e. for a positive trace of $\alpha_{ij}$, it points in the direction of the mean vorticity, i.e. perpendicular to the plane of the shear flow. Our simulations support this expression for low shear and magnetic Reynolds number. The transport coefficients depend on the wavenumber of the mean flow in a Lorentzian fashion, just as for non-shearing turbulence.

We present local numerical models of accretion disk turbulence driven by the magnetorotational instability with varying shear rate. The resulting turbulent stresses are compared with predictions of a closure model in which triple correlations are modelled in terms of quadratic correlations. This local model uses five nondimensional parameters to describe the properties of the flow. We attempt to determine these closure parameters for our simulations and find that the model does produce qualitatively correct behaviour. In addition, we present results concerning the shear rate dependency of the magnetic to kinetic energy ratio. We find both the turbulent stress ratio and the total stress to be strongly dependent on the shear rate.

(abridged) Aims: 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. Methods: 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. 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 alpha-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, 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.

Aims: We determine the components of the $\Lambda$-effect tensor that quantifies the contributions to the turbulent momentum transport even for uniform rotation. Methods: Three-dimensional numerical simulations are used to study turbulent transport in triply periodic cubes under the influence of rotation and anisotropic forcing. Comparison is made with analytical results obtained via the so-called minimal tau-approximation. Results: In the case where the turbulence intensity in the vertical direction dominates, the vertical stress is always negative. This situation is expected to occur in stellar convection zones. The horizontal component of the stress is weaker and exhibits a maximum at latitude 30 degrees - regardless of how rapid the rotation is. The minimal tau-approximation captures many of the qualitative features of the numerical results, provided the relaxation time tau is close to the turnover time, i.e. the Strouhal number is of order unity.

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 alpha 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.

The motivation for considering distributed large scale dynamos in the solar context is reviewed in connection with the magnetic helicity constraint. Preliminary accounts of 3-dimensional direct numerical simulations (in spherical shell segments) and simulations of 2-dimensional mean field models (in spherical shells) are presented. Interesting similarities as well as some differences are noted.

Homogeneous anisotropic turbulence simulations are used to determine off-diagonal components of the Reynolds stress tensor and its parameterization in terms of turbulent viscosity and Lambda-effect. The turbulence is forced in an anisotropic fashion by enhancing the strength of the forcing in the vertical direction. The Coriolis force is included with a rotation axis inclined relative to the vertical direction. The system studied here is significantly simpler than that of turbulent stratified convection which has often been used to study Reynolds stresses. Certain puzzling features of the results for convection, such as sign changes or highly concentrated latitude distributions, are not present in the simpler system considered here.

Magnetic helicity effects are discussed in laboratory and astrophysical settings. First, dynamo action in Taylor-Green flows is discussed for different boundary conditions. However, because of the lack of scale separation with respect to the container, no large scale field is being produced and there is no resistively slow saturation phase as otherwise expected. Second, the build-up of a large scale field is demonstrated in a simulation where a localized magnetic eddy produces field on a larger scale if the eddy possesses a swirl. Such a set-up might be realizable experimentally through coils. Finally, new emerging issues regarding the connection between magnetic helicity and the solar dynamo are discussed. It is demonstrated that dynamos with a non-local (Babcock-Leighton type) alpha effect can also be catastrophically quenched, unless there are magnetic helicity fluxes.

We study convective overshooting by means of local 3D convection calculations. Using a mixing length model of the solar convection zone (CZ) as a guide, we determine the Coriolis number (Co), which is the inverse of the Rossby number, to be of the order of ten or larger at the base of the solar CZ. Therefore we perform convection calculations in the range Co = 0...10 and interpret the value of Co realised in the calculation to represent a depth in the solar CZ. In order to study the dependence on rotation, we compute the mixing length parameters alpha_T and alpha_u relating the temperature and velocity fluctuations, respectively, to the mean thermal stratification. We find that the mixing length parameters for the rapid rotation case, corresponding to the base of the solar CZ, are 3-5 times smaller than in the nonrotating case. Introducing such depth-dependent alpha into a solar structure model employing a non-local mixing length formalism results in overshooting which is approximately proportional to alpha at the base of the CZ. Although overshooting is reduced due to the reduced alpha, a discrepancy with helioseismology remains due to the steep transition to the radiative temperature gradient. In comparison to the mixing length models the transition at the base of the CZ is much gentler in the 3D models. It was suggested recently (Rempel 2004) that this discrepancy is due to the significantly larger (up to seven orders of magnitude) input energy flux in the 3D models in comparison to the Sun and solar models, and that the 3D calculations should be able to approach the mixing length regime if the input energy flux is decreased by a moderate amount. We present results from local convection calculations which support this conjecture.

(abridged) Results from kinematic solar dynamo models employing alpha-effect and turbulent pumping from local convection calculations are presented. We estimate the magnitude of these effects to be around 2-3 m/s. The rotation profile of the Sun as obtained from helioseismology is applied. We obtain an estimate of the ratio of the two induction effects, C_alpha/C_Omega \approx 10^-3, which we keep fixed in all models. We also include a one-cell meridional circulation pattern having a magnitude of 10-20 m/s near the surface and 1-2 m/s at the bottom of the convection zone. The model essentially represents a distributed turbulent dynamo, as the alpha-effect is nonzero throughout the convection zone, although it concentrates near the bottom of the convection zone obtaining a maximum around 30 degrees of latitude. Turbulent pumping of the mean fields is predominantly down- and equatorward. We find that, when all these effects are included in the model, it is possible to correctly reproduce many features of the solar activity cycle, namely the correct equatorward migration at low latitudes and the polar branch at high latitudes, and the observed negative sign of B_r B_phi. Although the activity clearly shifts towards the equator in comparison to previous models due to the combined action of the alpha-effect peaking at midlatitudes, meridional circulation and latitudinal pumping, most of the activity still occurs at too high latitudes (between 5-60 degrees). Other problems include the relatively narrow parameter space within which the preferred solution is dipolar (A0), and the somewhat too short cycle lengths of the solar-type solutions. The role of the surface shear layer is found to be important only in the case where the alpha-effect has an appreciable magnitude near the surface.

Aims: The Strouhal number (St), which is a nondimensional measure of the correlation time, is determined from numerical models of convection. The Strouhal number arises in the mean-field theories of angular momentum transport and magnetic field generation, where its value determines the validity of certain widely used approximations, such as the first order smoothing (hereafter FOSA). More specifically, the relevant transport coefficients can be calculated by means of a cumulative series expansion if St < St_crit \approx 1. Methods: We define the Strouhal number as the ratio of the correlation and turnover times, which we determine separately, the former from the autocorrelation of velocity, and the latter by following test particles embedded in the flow. Results: We find that the Strouhal numbers are, generally, of the order of 0.1 to 0.4 which is close to the critical value above which deviations from FOSA become significant. Increasing the rotational influence tends to shorten both timescales in such a manner that St decreases. However, we do not find a clear trend as function of the Rayleigh number for the parameter range explored in the present study.

Aims. The alpha- and gamma-effects, which are responsible for the generation and turbulent pumping of large scale magnetic fields, respectively, due to passive advection by convection are determined in the rapid rotation regime corresponding to the deep layers of the solar convection zone. Methods. A 3D rectangular local model is used for solving the full set of MHD equations in order to compute the electromotive force (emf), E = <u x b>, generated by the interaction of imposed weak gradient-free magnetic fields and turbulent convection with varying rotational influence and latitude. By expanding the emf in terms of the mean magnetic field, E_i = a_ij <B_j>, all nine components of a_ij are computed. The diagonal elements of a_ij describe the alpha-effect, whereas the off-diagonals represent magnetic pumping. The latter is essentially the advection of magnetic fields by means other than the underlying large-scale velocity field. Comparisons are made to analytical expressions of the coefficients derived under the first-order smoothing approximation (FOSA). Results. In the rapid rotation regime the latitudinal dependence of the alpha-components responsible for the generation of the azimuthal and radial fields does not exhibit a peak at the poles, as is the case for slow rotation, but at a latitude of about 30 degrees. The magnetic pumping is predominantly radially down- and latitudinally equatorward as in earlier studies. The numerical results compare surprisingly well with analytical expressions derived under first-order smoothing, although the present calculations are expected to lie near the limits of the validity range of FOSA.

A two-dimensional mean field dynamo model is solved where magnetic helicity conservation is fully included. The model has a negative radial velocity gradient giving rise to equatorward migration of magnetic activity patterns. In addition the model develops longitudinal variability with activity patches travelling in longitude. These patches may be associated with active longitudes.

We study the mixing length concept in comparison to three-dimensional numerical calculations of convection with rotation. In a limited range, the velocity and temperature fluctuations are linearly proportional to the superadiabaticity, as predicted by the mixing length concept and in accordance with published results. The effects of rotation are investigated by varying the Coriolis number, Co = 2 Omega tau, from zero to roughly ten, and by calculating models at different latitudes. We find that \alpha decreases monotonically as a function of the Coriolis number. This can be explained by the decreased spatial scale of convection and the diminished efficiency of the convective energy transport, the latter of which leads to a large increase of the superadibaticity, \delta = \nabla - \nabla_ad as function of Co. Applying a decreased mixing length parameter in a solar model yields very small differences in comparison to the standard model within the convection zone. The main difference is the reduction of the overshooting depth, and thus the depth of the convection zone, when a non-local version of the mixing length concept is used. Reduction of \alpha by a factor of roughly 2.5 is sufficient to reconcile the difference between the model and helioseismic results. The numerical results indicate reduction of \alpha by this order of magnitude.

Three closely related stumbling blocks of solar mean field dynamo theory are discussed: how dominant are the small scale fields, how is the alpha effect quenched, and whether magnetic and current helicity fluxes alleviate the quenching? It is shown that even at the largest currently available resolution there is no clear evidence of power law scaling of the magnetic and kinetic energy spectra in turbulence. However, using subgrid scale modeling, some indications of asymptotic equipartition can be found. The frequently used first order smoothing approach to calculate the alpha effect and other transport coefficients is contrasted with the superior minimal tau approximation. The possibility of catastrophic alpha quenching is discussed as a result of magnetic helicity conservation. Magnetic and current helicity fluxes are shown to alleviate catastrophic quenching in the presence of shear. Evidence for strong large scale dynamo action, even in the absence of helicity in the forcing, is presented.

We determine the Strouhal number (hereafter St), which is essentially a nondimensional measure of the correlation time, from numerical calculations of convection. We use two independent methods to estimate St. Firstly, we apply the minimal tau-approximation (MTA) on the equation of the time derivative of the Reynolds stress. A relaxation time is obtained from which St can be estimated by normalising with a typical turnover time. Secondly, we calculate the correlation and turnover times separately, the former from the autocorrelation of velocity and the latter by following test particles embedded in the flow. We find that the Strouhal number is in general of the order of 0.1 to 1, i.e. rather large in comparison to the typical assumption in the mean-field theories that St << 1. However, there is a clear decreasing trend as function of the Rayleigh number and increasing rotation. Furthermore, for the present range of parameters the decrease of St does not show signs of saturation, indicating that in stellar convection zones, where the Rayleigh numbers are much larger, the Strouhal number may indeed be significantly smaller.

We study stellar convection using a local three-dimensional MHD model, with which we investigate the influence of rotation and large-scale magnetic fields on the turbulent momentum and heat transport. The former is studied by computing the Reynolds stresses, the latter by calculating the correlation of velocity and temperature fluctuations, both as functions of rotation and latitude. We find that the horisontal correlation, Q_(theta phi), capable of generating horisontal differential rotation, is mostly negative in the southern hemisphere for Coriolis numbers exceeding unity, corresponding to equatorward flux of angular momentum in accordance with solar observations. The radial component Q_(r phi) is negative for slow and intermediate rotation indicating inward transport of angular momentum, while for rapid rotation, the transport occurs outwards. Parametrisation in terms of the mean-field Lambda-effect shows qualitative agreement with the turbulence model of Kichatinov & R\"udiger (1993) for the horisontal part H \propto Q_(theta phi)/cos(theta), whereas for the vertical part, V \propto Q_(r phi)/sin(theta), agreement only for intermediate rotation exists. The Lambda-coefficients become suppressed in the limit of rapid rotation, this rotational quenching being stronger for the V component than for H. We find that the stresses are enhanced by the presence of the magnetic field for field strengths up to and above the equipartition value, without significant quenching. Concerning the turbulent heat transport, our calculations show that the transport in the radial direction is most efficient at the equatorial regions, obtains a minimum at midlatitudes, and shows a slight increase towards the poles. The latitudinal heat transport does not show a systematic trend as function of latitude or rotation.

The standard theory of the solar cycle in terms of an alpha-Omega dynamo hinges on a proper understanding of the nonlinear alpha effect. Boundary conditions play a surprisingly important role in determining the magnitude of alpha. For closed boundaries, the total magnetic helicity is conserved, and since the alpha effect produces magnetic helicity of one sign in the large scale field, it must simultaneously produce magnetic helicity of the opposite sign. It is this secondary magnetic helicity that suppresses the dynamo in a potentially catastrophic fashion. Open boundaries allow magnetic helicity to be lost. Simulations are presented that allow an estimate of alpha in the presence of open or closed boundaries, either with or without solar-like differential rotation. In all cases the sign of the magnetic helicity agrees with that observed at the solar surface (negative in the north, positive in the south), where significant amounts of magnetic helicity can be ejected via coronal mass ejections. It is shown that open boundaries tend to alleviate catastrophic alpha quenching. The importance of looking at current helicity instead of magnetic helicity is emphasized and the conceptual advantages are discussed.

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.