Complex mixing and magnetic field generation occurs within stellar interiors particularly where there is a strong shear flow. To obtain a comprehensive understanding of these processes, it is necessary to study the complex dynamics of shear regions. Due to current observational limitations, it is necessary to investigate the inevitable small-scale dynamics via numerical calculations. Here, we examine direct numerical calculations of a local model of unstable shear flows in a compressible polytropic fluid primarily in a two-dimensional domain, where we focus on determining how key parameters affect the global properties and characteristics of the resulting saturated turbulent phase. We consider the effect of varying both the viscosity and the thermal diffusivity on the non-linear evolution. Moreover, our main focus is to understand the global properties of the saturated phase, in particular estimating for the first time the spread of the shear region from an initially hyperbolic tangent velocity profile. We find that the vertical extent of the mixing region in the saturated regime is generally determined by the initial Richardson number of the system. Further, the characteristic quantities of the turbulence, i.e. typical length-scale and the root-mean-square velocity are found to depend on both the Richardson number, and the thermal diffusivity. Finally, we present our findings of our investigation into saturated flows of a `secular' shear instability in the low P\'eclet number regime with large Richardson numbers.

In stellar interiors shear flows play an important role in many physical processes. So far helioseismology provides only large-scale measurements, and so the small-scale dynamics remains insufficiently understood. To draw a connection between observations and three-dimensional DNS of shear driven turbulence, we investigate horizontally averaged profiles of the numerically obtained mean state. We focus here on just one of the possible methods that can maintain a shear flow, namely the average relaxation method. We show that although some systems saturate by restoring linear marginal stability this is not a general trend. Finally, we discuss the reason that the results are more complex than expected.

Shear flows are ubiquitous in astrophysical objects including planetary and stellar interiors, where their dynamics can have significant impact on thermo-chemical processes. Investigating the complex dynamics of shear flows requires numerical calculations that provide a long time evolution of the system. To achieve a sufficiently long lifetime in a local numerical model the system has to be forced externally. However, at present, there exist several different forcing methods to sustain large-scale shear flows in local models. In this paper we examine and compare various methods used in the literature in order to resolve their respective applicability and limitations. These techniques are compared during the exponential growth phase of a shear flow instability, such as the Kelvin-Helmholtz (KH) instability, and some are examined during the subsequent non-linear evolution. A linear stability analysis provides reference for the growth rate of the most unstable modes in the system and a detailed analysis of the energetics provides a comprehensive understanding of the energy exchange during the system's evolution. Finally, we discuss the pros and cons of each forcing method and their relation with natural mechanisms generating shear flows.

Shear flows have an important impact on the dynamics in an assortment of different astrophysical objects including accreditation discs and stellar interiors. Investigating shear flow instabilities in a polytropic atmosphere provides a fundamental understanding of the motion in stellar interiors where turbulent motions, mixing processes, as well as magnetic field generation takes place. Here, a linear stability analysis for a fully compressible fluid in a two-dimensional Cartesian geometry is carried out. Our study focuses on determining the critical Richardson number for different Mach numbers and the destabilising effects of high thermal diffusion. We find that there is a deviation of the predicted stability threshold for moderate Mach number flows along with a significant effect on the growth rate of the linear instability for small P\'eclet numbers. We show that in addition to a Kelvin-Helmholtz instability a Holmboe instability can appear and we discuss the implication of this in stellar interiors.

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.

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.

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.

Magnetic fields are known to reside in many astrophysical objects and are now believed to be crucially important for the creation of phenomena on a wide variety of scales. However, the role of the magnetic field in the bodies that we observe has not always been clear. In certain situations, the importance of a magnetic field has been over looked on the grounds that the large-scale magnetic field was believed to be too weak to play and important role in the dynamics. In this article I discuss some of the recent developments concerning magnetic fields in stars, planets and accretion discs. I choose to emphasise some of the situations where it has been suggested that weak magnetic fields may play a more significant role than previously thought. At the end of the article I list some of the questions to be answered in the future.

For more than a decade, the so-called shearing box model has been used to study the fundamental local dynamics of accretion discs. This approach has proved to be very useful because it allows high resolution and long term studies to be carried out, studies that would not be possible for a global disc. Localised disc studies have largely focused on examining the rate of enhanced transport of angular momentum, essentially a sum of the Reynolds and Maxwell stresses. The dominant radial-azimuthal component of this stress tensor is, in the classic Shakura-Sunayaev model, expressed as a constant alpha times the pressure. Previous studies have estimated alpha based on a modest number of orbital times. Here we use much longer baselines, and perform a cumulative average for alpha. Great care must be exercised when trying to extract numerical alpha values from simulations: dissipation scales, computational box aspect ratio, and even numerical algorithms all affect the result. This study suggests that estimating alpha becomes more, not less, difficult as computational power increases.

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.

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.