We study the dynamics associated with the extension of turbulent convective motions from a convection zone (CZ) into a stable region (RZ) that lies below the latter. For that purpose, we have run a series of three-dimensional direct numerical simulations solving the Navier-Stokes equations under the Boussinesq approximation in a spherical shell geometry. We observe that the overshooting of the turbulent motions into the stably stratified region depends on three different parameters: the relative stability of the RZ, the transition width between the two, and the intensity of the turbulence. In the cases studied, these motions manage to partially alter the thermal stratification and induce thermal mixing, but not so efficiently as to extend the nominal CZ further down into the stable region. We find that the kinetic energy below the convection zone can be modeled by a half-Gaussian profile whose amplitude and width can be predicted a priori for all of our simulations. We examine different dynamical lengthscales related to the depth of the extension of the motions into the RZ, and we find that they all scale remarkably well with a lengthscale that stems from a simple energetic argument. We discuss the implications of our findings for 1D stellar evolution calculations.

The buoyant transport of magnetic fields from the solar interior towards the surface plays an important role in the emergence of active regions, the formation of sunspots and the overall solar dynamo. Observations suggest that toroidal flux concentrations often referred to as "flux tubes", rise from their region of initiation likely in the solar tachocline towards the solar surface due to magnetic buoyancy. Many studies have assumed the existence of such magnetic structures and studied the buoyant rise of an isolated flux tube in a quiescent, field-free environment. Here, motivated by their formation (Cline et al. 2003; Brummell et al. 2002), we relax the latter assumption and study the rise of a toroidal flux tube embedded in a large-scale poloidal background magnetic field. We find that the presence of the large-scale background field severely affects the dynamics of the rising tube. A relatively weak background field, as low as 6% of the tube strength, can destroy the rise of a tube that would otherwise rise in the absence of the background field. Surprisingly, the rise of tubes with one sign of the twist is suppressed by a significantly weaker background field than the other. This demonstrates a potential mechanism for the selection of the preferred helicity of rising and emerging tubes for the solar case that is commensurate with many features of the hemispherical rule.

We present a local but fully nonlinear model of the solar tachocline, using three-dimensional direct numerical simulations. The tachocline forms naturally as a statistically steady balance between Coriolis, pressure, buoyancy and Lorentz forces beneath a turbulent convection zone. Uniform rotation is maintained in the radiation zone by a primordial magnetic field, which is confined by meridional flows in the tachocline and convection zone. Such balanced dynamics has previously been found in idealised laminar models, but never in fully self-consistent numerical simulations.

We examine the dynamics associated with weakly compressible convection in a spherical shell by running 3D direct numerical simulations using the Boussinesq formalism [1]. Motivated by problems in astrophysics, we assume the existence of a finite adiabatic temperature gradient $\nabla T_{\rm{ad}}$ and use mixed boundary conditions for the temperature with fixed flux at the inner boundary and fixed temperature at the outer boundary. This setup is intrinsically more asymmetric than the more standard case of Rayleigh-B\'{e}nard convection in liquids between parallel plates with fixed temperature boundary conditions. Conditions where there is substantial asymmetry can cause a dramatic change in the nature of convection and we demonstrate that this is the case here. The flows can become pressure- rather than buoyancy- dominated leading to anomalous heat transport by upflows. Counter-intuitively, the background temperature gradient $\nabla\bar{T}$ can develop a subadiabatic layer (where $\boldsymbol{g}\cdot\nabla\bar{T}<\boldsymbol{g}\cdot\nabla T_{\rm{ad}}$, where $\boldsymbol{g}$ is gravity) although convection remains vigorous at every point across the shell. This indicates a high degree of non-locality.

Fingering convection (otherwise known as thermohaline convection) is an instability that occurs in stellar radiative interiors in the presence of unstable compositional gradients. Numerical simulations have been used in order to estimate the efficiency of mixing induced by this instability. However, fully three-dimensional (3D) computations in the parameter regime appropriate for stellar astrophysics (i.e. low Prandtl number) are prohibitively expensive. This raises the question of whether two-dimensional (2D) simulations could be used instead to achieve the same goals. In this work, we address this issue by comparing the outcome of 2D and 3D simulations of fingering convection at low Prandtl number. We find that 2D simulations are never appropriate. However, we also find that the required 3D computational domain does not have to be very wide: the third dimension need only contain a minimum of two wavelengths of the fastest-growing linearly unstable mode to capture the essentially 3D dynamics of small-scale fingering. Narrow domains, however, should still be used with caution since they could limit the subsequent development of any large-scale dynamics typically associated with fingering convection.

Transport by meridional flows has significant consequences for stellar evolution, but is difficult to capture in global-scale numerical simulations because of the wide range of timescales involved. Stellar evolution models therefore usually adopt parameterizations for such transport based on idealized laminar or mean-field models. Unfortunately, recent attempts to model this transport in global simulations have produced results that are not consistent with any of these idealized models. In an effort to explain the discrepancies between global simulations and idealized models, we here use three-dimensional local Cartesian simulations of compressible convection to study the efficiency of transport by meridional flows below a convection zone in several parameter regimes of relevance to the Sun and solar-type stars. In these local simulations we are able to establish the correct ordering of dynamical timescales, although the separation of the timescales remains unrealistic. We find that, even though the generation of internal waves by convective overshoot produces a high degree of time dependence in the meridional flow field, the mean flow has the qualitative behavior predicted by laminar, "balanced" models. In particular, we observe a progressive deepening, or "burrowing", of the mean circulation if the local Eddington-Sweet timescale is shorter than the viscous diffusion timescale. Such burrowing is a robust prediction of laminar models in this parameter regime, but has never been observed in any previous numerical simulation. We argue that previous simulations therefore underestimate the transport by meridional flows.

We present a numerical study of dynamo action in a conducting fluid encased in a metallic spherical shell. Motions in the fluid are driven by differential rotation of the outer metallic shell, which we refer to as "the wall". The two hemispheres of the wall are held in counter-rotation, producing a steady, axisymmetric interior flow consisting of differential rotation and a two-cell meridional circulation with radial inflow in the equatorial plane. From previous studies, this type of flow is known to maintain a stationary equatorial dipole by dynamo action if the magnetic Reynolds number is larger than about 300 and if the outer boundary is electrically insulating. We vary independently the thickness, electrical conductivity, and magnetic permeability of the wall to determine their effect on the dynamo action. The main results are: (a) Increasing the conductivity of the wall hinders the dynamo by allowing eddy currents within the wall, which are induced by the relative motion of the equatorial dipole field and the wall. This processes can be viewed as a skin effect or, equivalently, as the tearing apart of the dipole by the differential rotation of the wall, to which the field lines are anchored by high conductivity. (b) Increasing the magnetic permeability of the wall favors dynamo action by constraining the magnetic field lines in the fluid to be normal to the wall, thereby decoupling the fluid from any induction in the wall. (c) Decreasing the wall thickness limits the amplitude of the eddy currents, and is therefore favorable for dynamo action, provided that the wall is thinner than the skin depth. We explicitly demonstrate these effects of the wall properties on the dynamo field by deriving an effective boundary condition in the limit of vanishing wall thickness.

We have performed numerical simulations of boundary-driven dynamos using a three-dimensional non-linear magnetohydrodynamical model in a spherical shell geometry. A conducting fluid of magnetic Prandtl number Pm=0.01 is driven into motion by the counter-rotation of the two hemispheric walls. The resulting flow is of von K\'arm\'an type, consisting of a layer of zonal velocity close to the outer wall and a secondary meridional circulation. Above a certain forcing threshold, the mean flow is unstable to non-axisymmetric motions within an equatorial belt. For fixed forcing above this threshold, we have studied the dynamo properties of this flow. The presence of a conducting outer wall is essential to the existence of a dynamo at these parameters. We have therefore studied the effect of changing the material parameters of the wall (magnetic permeability, electrical conductivity, and thickness) on the dynamo. In common with previous studies, we find that dynamos are obtained only when either the conductivity or the permeability is sufficiently large. However, we find that the effect of these two parameters on the dynamo process are different and can even compete to the detriment of the dynamo. Our self-consistent approach allow us to analyze in detail the dynamo feedback loop. The dynamos we obtain are typically dominated by an axisymmetric toroidal magnetic field and an axial dipole component. We show that the ability of the outer shear layer to produce a strong toroidal field depends critically on the presence of a conducting outer wall, which shields the fluid from the vacuum outside. The generation of the axisymmetric poloidal field, on the other hand, occurs in the equatorial belt and does not depend on the wall properties.

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.

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.

Meridional flows with velocities of a few meters per second are observed in the uppermost regions of the solar convection zone. The amplitude and pattern of the flows deeper in the solar interior, in particular near the top of the radiative region, are of crucial importance to a wide range of solar magnetohydrodynamical processes. In this paper, we provide a systematic study of the penetration of large-scale meridional flows from the convection zone into the radiative zone. In particular, we study the effects of the assumed boundary conditions applied at the convective-radiative interface on the deeper flows. Using simplified analytical models in conjunction with more complete numerical methods, we show that penetration of the convectively-driven meridional flows into the deeper interior is not necessarily limited to a shallow Ekman depth but can penetrate much deeper, depending on how the convective-radiative interface flows are modeled.