Recent numerical simulations of dynamo action resulting from rotating convection have revealed some serious problems in applying the standard picture of mean field electrodynamics at high values of the magnetic Reynolds number, and have thereby underlined the difficulties in large-scale magnetic field generation in this regime. Here we consider kinematic dynamo processes in a rotating convective layer of Boussinesq fluid with the additional influence of a large-scale horizontal velocity shear. Incorporating the shear flow enhances the dynamo growth rate and also leads to the generation of significant magnetic fields on large scales. By the technique of spectral filtering, we analyse the modes in the velocity that are principally responsible for dynamo action, and show that the magnetic field resulting from the full flow relies crucially on a range of scales in the velocity field. Filtering the flow to provide a true separation of scales between the shear and the convective flow also leads to dynamo action; however, the magnetic field in this case has a very different structure from that generated by the full velocity field. We also show that the nature of the dynamo action is broadly similar irrespective of whether the flow in the absence of shear can support dynamo action.

We here report calculations of magnetic buoyancy instabilities of a sheared magnetic layer where two separate zones are unstable. The idea is to study the possible generation of large-scale helical structures which could then rise through a stellar convection zone and emerge at the surface to create active regions. The calculations shown here are a follow-up of the work of Favier et al. (2012) where the instability developed in a weakly magnetized atmosphere, consisting of a uniform field oriented in a different direction from the unstable layer below. Here, the top layer representing the atmosphere is itself unstable to buoyancy instabilities and thus quickly creates a more complex magnetic configuration with which the layer below will interact. We also find in this case that the accumulation of magnetic tension between the two unstable layers favors the creation of large-scale helical structures.

We consider the three-dimensional instability of a layer of horizontal magnetic field in a polytropic atmosphere where, contrary to previous studies, the field lines in the initial state are not unidirectional. We show that if the twist is initially concentrated inside the unstable layer, the modifications of the instability reported by several authors (see e.g. Cattaneo et al. (1990)) are only observed when the calculation is restricted to two dimensions. In three dimensions, the usual interchange instability occurs, in the direction fixed by the field lines at the interface between the layer and the field-free region. We therefore introduce a new configuration: the instability now develops in a weakly magnetised atmosphere where the direction of the field can vary with respect to the direction of the strong unstable field below, the twist being now concentrated at the upper interface. Both linear stability analysis and non-linear direct numerical simulations are used to study this configuration. We show that from the small-scale interchange instability, large-scale twisted coherent magnetic structures are spontaneously formed, with possible implications to the formation of active regions from a deep-seated solar magnetic field.

We perform idealised numerical simulations of magnetic buoyancy instabilities in a model of the solar tachocline. We introduce a simplified model of magnetic flux pumping in an upper layer (the convection zone), and study the effects of its inclusion on the evolution of buoyancy instabilities in a lower layer (the radiative interior). We study its effects on the instability of both a preconceived magnetic slab and of a shear-generated magnetic layer. In the former, we find that in the regime in which the downward pumping velocity is comparable with the Alfven speed of the magnetic layer, flux pumping is able to hold back the bulk of the magnetic field, with only small pockets of strong field able to rise into the upper layer. In simulations in which the magnetic layer is generated by shear, we find that the shear velocity is not necessarily required to exceed that of the pumping (therefore the kinetic energy of the shear is not required to exceed that of the overlying convection), for strong localised pockets of magnetic field to be produced which can rise into the upper layer. This is because magnetic flux pumping acts to store the field below the interface, allowing it to be amplified both by the shear, and by vortical fluid motions, until pockets of field can achieve sufficient strength to rise into the upper layer. In addition, we find that the interface between the two layers is a natural location for the production of strong vertical gradients in the magnetic field. If these gradients are sufficiently strong to allow the development of magnetic buoyancy instabilities, strong shear is not necessarily required to drive them (c.f. previous work by Vasil & Brummell). We find that the addition of magnetic flux pumping appears to be able to assist shear-driven magnetic buoyancy in producing strong flux concentrations that can rise into the convection zone from the radiative interior.

Recently, Silvers, Vasil, Brummell, & Proctor (2009), using numerical simulations, confirmed the existence of a double diffusive magnetic buoyancy instability of a layer of horizontal magnetic field produced by the interaction of a shear velocity field with a weak vertical field. Here, we demonstrate the longer term nonlinear evolution of such an instability in the simulations. We find that a quasi two-dimensional interchange instability rides (or "surfs") on the growing shear-induced background downstream field gradients. The region of activity expands since three-dimensional perturbations remain unstable in the wake of this upward-moving activity front, and so the three-dimensional nature becomes more noticeable with time.

Numerical simulations are an important tool in furthering our understanding of turbulent dynamo action, a process that occurs in a vast range of astrophysical bodies. It is important in all computational work that comparisons are made between different codes and, if non-trivial differences arise, that these are explained. Kapyla et al (2010: MNRAS 402, 1458) describe an attempt to reproduce the results of Hughes & Proctor (2009: PRL 102, 044501) and, by employing a different methodology, they arrive at very different conclusions concerning the mean electromotive force and the generation of large-scale fields. Here we describe why the simulations of Kapyla et al (2010) are simply not suitable for a meaningful comparison, since they solve different equations, at different parameter values and with different boundary conditions. Furthermore we describe why the interpretation of Kapyla et al (2010) of the calculation of the alpha-effect is inappropriate and argue that the generation of large-scale magnetic fields by turbulent convection remains a problematic issue.

We investigate the interaction of a fluctuating alpha-effect with large-scale shear in a simple nonlinear 1-dimensional dynamo wave model. We firstly extend the calculations of Proctor (2007, MNRAS, 41, L39-L42) to include spatial variation of the fluctuations, and find that there can be a mechanism for magnetic field generation, even when the mean alpha is zero, provided the spatiotemporal spectrum of the fluctuations has an appropriate form. We investigate mean-field dynamo action when the new term arising from the fluctuations is non-zero, and present results concerning the stability and frequency of the solutions and parity selection in the nonlinear regime. The relation between the asymptotic theory and explicit simulation of a traditional mean-field model with a fluctuating function for the alpha-effect term is discussed.

A subcritical transition to turbulence in magnetized Keplerian shear flows is investigated using a statistical approach. Three-dimensional numerical simulations of the shearing box equations with zero net magnetic flux are employed to determine the transition from decaying to sustained turbulence as a function of the magnetic Reynolds number Rm. The results reveal no clear transition to sustained turbulence as the average lifetime of the transients grows as an exponential function of Rm, in accordance with a type-II supertransient law.

The Sun is a magnetic star whose cyclic activity is thought to be linked to internal dynamo mechanisms. A combination of numerical modelling with various levels of complexity is an efficient and accurate tool to investigate such intricate dynamical processes. We investigate the role of the magnetic buoyancy process in 2D Babcock-Leighton dynamo models, by modelling more accurately the surface source term for poloidal field. Methods. To do so, we reintroduce in mean-field models the results of full 3D MHD calculations of the non-linear evolution of a rising flux tube in a convective shell. More specifically, the Babcock-Leighton source term is modified to take into account the delay introduced by the rise time of the toroidal structures from the base of the convection zone to the solar surface. We find that the time delays introduced in the equations produce large temporal modulation of the cycle amplitude even when strong and thus rapidly rising flux tubes are considered. Aperiodic modulations of the solar cycle appear after a sequence of period doubling bifurcations typical of non-linear systems. The strong effects introduced even by small delays is found to be due to the dependence of the delays on the magnetic field strength at the base of the convection zone, the modulation being much less when time delays remain constant. We do not find any significant influence on the cycle period except when the delays are made artificially strong. A possible new origin of the solar cycle variability is here revealed. This modulated activity and the resulting butterfly diagram are then more compatible with observations than what the standard Babcock-Leighton model produces.

The coupled equations that describe the effect of large-scale magnetic and velocity fields on forced high-diffusivity magnetohydrodynamic flows are investigated through an extension of mean field electrodynamics. Our results generalise those of R\"adler & Brandenburg (2010), who consider a similar situation but assume that the effect of the Lorentz force on the momentum equation can be neglected. New mean coupling terms are shown to appear, which can lead to large-scale growth of magnetic and velocity fields even when the usual $\alpha$-effects are absent.

We consider the linear stability of two-dimensional nonlinear magnetohydrodynamic basic states to long-wavelength three-dimensional perturbations. Following Hughes & Proctor (2009a), the 2D basic states are obtained from a specific forcing function in the presence of an initially uniform mean field of strength $\mathcal{B}$. By extending to the nonlinear regime the kinematic analysis of Roberts (1970), we show that it is possible to predict the growth rate of these perturbations by applying mean field theory to \textit{both} the momentum and the induction equations. If $\mathcal{B}=0$, these equations decouple and large-scale magnetic and velocity perturbations may grow via the kinematic $\alpha$-effect and the AKA instability respectively. However, if $\mathcal{B} \neq 0$, the momentum and induction equations are coupled by the Lorentz force; in this case, we show that four transport tensors are now necessary to determine the growth rate of the perturbations. We illustrate these situations by numerical examples; in particular, we show that a mean field description of the nonlinear regime based solely on a quenched $\alpha$ coefficient is incorrect.

Previous theoretical work has speculated about the existence of double-diffusive magnetic buoyancy instabilities of a dynamically evolving horizontal magnetic layer generated by the interaction of forced vertically sheared velocity and a background vertical magnetic field. Here we confirm numerically that if the ratio of the magnetic to thermal diffusivities is sufficiently low then such instabilities can indeed exist, even for high Richardson number shear flows. Magnetic buoyancy may therefore occur via this mechanism for parameters that are likely to be relevant to the solar tachocline, where regular magnetic buoyancy instabilities are unlikely.

Motivated by the interface model for the solar dynamo, this paper explores the complex magnetohydrodynamic interactions between convective flows and shear-driven instabilities. Initially, we consider the dynamics of a forced shear flow across a convectively-stable polytropic layer, in the presence of a vertical magnetic field. When the imposed magnetic field is weak, the dynamics are dominated by a shear flow (Kelvin-Helmholtz type) instability. For stronger fields, a magnetic buoyancy instability is preferred. If this stably stratified shear layer lies below a convectively unstable region, these two regions can interact. Once again, when the imposed field is very weak, the dynamical effects of the magnetic field are negligible and the interactions between the shear layer and the convective layer are relatively minor. However, if the magnetic field is strong enough to favour magnetic buoyancy instabilities in the shear layer, extended magnetic flux concentrations form and rise into the convective layer. These magnetic structures have a highly disruptive effect upon the convective motions in the upper layer.

The transition to intermittent mean--field dynamos is studied using numerical simulations of isotropic magnetohydrodynamic turbulence driven by a helical flow. The low-Prandtl number regime is investigated by keeping the kinematic viscosity fixed while the magnetic diffusivity is varied. Just below the critical parameter value for the onset of dynamo action, a transient mean--field with low magnetic energy is observed. After the transition to a sustained dynamo, the system is shown to evolve through different types of intermittency until a large--scale coherent field with small--scale turbulent fluctuations is formed. Prior to this coherent field stage, a new type of intermittency is detected, where the magnetic field randomly alternates between phases of coherent and incoherent large--scale spatial structures. The relevance of these findings to the understanding of the physics of mean--field dynamo and the physical mechanisms behind intermittent behavior observed in stellar magnetic field variability are discussed.

We re-examine the nature of the turbulent magnetic diffusivity tensor of mean field electrodynamics and show that an inconsistency arises if it is calculated via consideration of time-independent magnetic fields. Specifically, the predicted growth rate of the mean magnetic field is, in general, incorrect. We describe how the traditional expansion procedure for the mean electromotive force should be extended, and illustrate the consistency of this approach by means of a perturbation analysis for a mean magnetic field varying on long spatial scales. Finally we examine the magnitude of this new contribution to the magnetic diffusion for a particular flow.

It is believed that some stars have two or more convection zones in close proximity near to the stellar photosphere. These zones are separated by convectively stable regions that are relatively narrow. Due to the close proximity of these regions it is important to construct mathematical models to understand the transport and mixing of passive and dynamic quantities. One key quantity of interest is a magnetic field, a dynamic vector quantity, that can drastically alter the convectively driven flows, and have an important role in coupling the different layers. In this paper we present the first investigation into the effect of an imposed magnetic field in such a geometry. We focus our attention on the effect of field strength and show that, while there are some similarities with results for magnetic field evolution in a single layer, new and interesting phenomena are also present in a three layer system.

By incorporating a large-scale shear flow into turbulent rotating convection, we show that a sufficiently strong shear can promote dynamo action in flows that in the absence of shear do not act as dynamos. Our results are consistent with a dynamo driven by either the shear-current effect or by the interaction between a fluctuating $\alpha$-effect and the velocity shear; they are though inconsistent with either a classical $\alpha^2$ or $\alpha \omega$ mean field dynamo.

Identifying generic physical mechanisms responsible for the generation of magnetic fields and turbulence in differentially rotating flows is fundamental to understand the dynamics of astrophysical objects such as accretion disks and stars. In this paper, we discuss the concept of subcritical dynamo action and its hydrodynamic analogue exemplified by the process of nonlinear transition to turbulence in non-rotating wall-bounded shear flows. To illustrate this idea, we describe some recent results on nonlinear hydrodynamic transition to turbulence and nonlinear dynamo action in rotating shear flows pertaining to the problem of turbulent angular momentum transport in accretion disks. We argue that this concept is very generic and should be applicable to many astrophysical problems involving a shear flow and non-axisymmetric instabilities of shear-induced axisymmetric toroidal velocity or magnetic fields, such as Kelvin-Helmholtz, magnetorotational, Tayler or global magnetoshear instabilities. In the light of several recent numerical results, we finally suggest that, similarly to a standard linear instability, subcritical MHD dynamo processes in high-Reynolds number shear flows could act as a large-scale driving mechanism of turbulent flows that would in turn generate an independent small-scale dynamo.

Kilogauss-strength magnetic fields are often observed in intergranular lanes at the photosphere in the quiet Sun. Such fields are stronger than the equipartition field $B_e$, corresponding to a magnetic energy density that matches the kinetic energy density of photospheric convection, and comparable with the field $B_p$ that exerts a magnetic pressure equal to the ambient gas pressure. We present an idealised numerical model of three-dimensional compressible magnetoconvection at the photosphere, for a range of values of the magnetic Reynolds number. In the absence of a magnetic field, the convection is highly supercritical and is characterised by a pattern of vigorous, time-dependent, ``granular'' motions. When a weak magnetic field is imposed upon the convection, magnetic flux is swept into the convective downflows where it forms localised concentrations. Unless this process is significantly inhibited by magnetic diffusion, the resulting fields are often much greater than $B_e$, and the high magnetic pressure in these flux elements leads to their being partially evacuated. Some of these flux elements contain ultra-intense magnetic fields that are significantly greater than $B_p$. Such fields are contained by a combination of the thermal pressure of the gas and the dynamic pressure of the convective motion, and they are constantly evolving. These ultra-intense fields develop owing to nonlinear interactions between magnetic fields and convection; they cannot be explained in terms of ``convective collapse'' within a thin flux tube that remains in overall pressure equilibrium with its surroundings.

The mean electromotive force and alpha effect are computed for a forced turbulent flow using a simple nonlinear dynamical model. The results are used to check the applicability of two basic analytic ansatze of mean-field magnetohydrodynamics - the second order correlation approximation (SOCA) and the tau approximation. In the numerical simulations the effective Reynolds number Re is 2-20, while the magnetic Prandtl number varies from 0.1 to $10^{7}$. We present evidence that the $\tau$ approximation may be appropriate in dynamical regimes where there is a small-scale dynamo. Catastrophic quenching of the $\alpha$ effect is found for high $P_{m}$. Our results indicate that for high $P_{m}$ SOCA gives a very large value of the $\alpha$ coefficient compared with the ``exact'' solution. The discrepancy depends on the properties of the random force that drives the flow, with a larger difference occuring for $\delta$-correlated force compared with that for a steady random force.

We present results from simulations to examine how the separation between two convectively unstable layers affect their interaction. We show that two convectively unstable layers remain connected via the overshooting plumes even when they are separated by several pressure scale heights.

We analyse the role of a fluctuating alpha-effect in alpha-omega dynamo models, and show that there is a mechanism for magnetic field generation, valid at large scale separation, deriving from the interaction of mean shear and a fluctuating alpha-effect. It is shown that this effect can act as a dynamo even in the absence of a mean alpha-effect, and that the timescale for dynamo waves is strongly affected by the presence of fluctuations.

A-type stars have a complex internal structure with the possibility of multiple convection zones. If not sufficiently separated, such zones will interact through the convectively stable regions that lie between them. It is therefore of interest to ask whether the typical conditions that exist within such stars are such that these convections zones can ever be considered as disjoint. In this paper we present results from numerical simulations that help in understanding how increasing the distance between the convectively unstable regions are likely to interact through the stable region that separates them. This has profound implications for mixing and transport within these stars.

A three-dimensional nonlinear dynamo process is identified in rotating plane Couette flow in the Keplerian regime. It is analogous to the hydrodynamic self-sustaining process in non-rotating shear flows and relies on the magneto-rotational instability of a toroidal magnetic field. Steady nonlinear solutions are computed numerically for a wide range of magnetic Reynolds numbers but are restricted to low Reynolds numbers. This process may be important to explain the sustenance of coherent fields and turbulent motions in Keplerian accretion disks, where all its basic ingredients are present.

Direct numerical simulations of incompressible nonhelical randomly forced MHD turbulence are used to demonstrate for the first time that the fluctuation dynamo exists in the limit of large magnetic Reynolds number Rm>>1 and small magnetic Prandtl number Pm<<1. The dependence of the critical Rm_c for dynamo on the hydrodynamic Reynolds number Re is obtained for 1<Re<6700. In the limit Pm<<1, Rm_c is about three times larger than for the previously well established dynamo at large and moderate Prandtl numbers: Rm_c<200 for Re>6000 compared to Rm_c~60 for Pm>=1. Is is not as yet possible to determine numerically whether the growth rate of the magnetic energy is ~Rm^{1/2} in the limit Rm->infinity, as it should if the dynamo is driven by the inertial-range motions at the resistive scale.

This paper is a detailed report on a programme of simulations used to settle a long-standing issue in the dynamo theory and demonstrate that the fluctuation dynamo exists in the limit of large magnetic Reynolds number Rm>>1 and small magnetic Prandtl number Pm<<1. The dependence of the critical Rm_c vs. the hydrodynamic Reynolds number Re is obtained for 1<Re<6700. In the limit Pm<<1, Rm_c is ~3 times larger than for Pm>1. The stability curve Rm_c(Re) (and, it is argued, the nature of the dynamo) is substantially different from the case of the simulations and liquid-metal experiments with a mean flow. It is not as yet possible to determine numerically whether the growth rate is ~Rm^{1/2} in the limit Re>>Rm>>1, as should be the case if the dynamo is driven by the inertial-range motions. The magnetic-energy spectrum in the low-Pm regime is qualitatively different from the Pm>1 case and appears to develop a negative spectral slope, although current resolutions are insufficient to determine its asymptotic form. At 1<Rm<Rm_c, the magnetic fluctuations induced via the tangling by turbulence of a weak mean field are investigated and the possibility of a k^{-1} spectrum above the resistive scale is examined. At low Rm<1, the induced fluctuations are well described by the quasistatic approximation; the k^{-11/3} spectrum is confirmed for the first time in direct numerical simulations.

We demonstrate a close analogy between a viscoelastic medium and an electrically conducting fluid containing a magnetic field. Specifically, the dynamics of the Oldroyd-B fluid in the limit of large Deborah number corresponds to that of a magnetohydrodynamic (MHD) fluid in the limit of large magnetic Reynolds number. As a definite example of this analogy, we compare the stability properties of differentially rotating viscoelastic and MHD flows. We show that there is an instability of the Oldroyd-B fluid that is physically distinct from both the inertial and elastic instabilities described previously in the literature, but is directly equivalent to the magnetorotational instability in MHD. It occurs even when the specific angular momentum increases outwards, provided that the angular velocity decreases outwards; it derives from the kinetic energy of the shear flow and does not depend on the curvature of the streamlines. However, we argue that the elastic instability of viscoelastic Couette flow has no direct equivalent in MHD.