L. J. Silvers

Inverse cascade and symmetry breaking in rapidly-rotating Boussinesq convection

In this paper, we present numerical simulations of rapidly rotating Rayleigh-Benard convection in the Boussinesq approximation with stress-free boundary conditions. At moderately low Rossby number and large Rayleigh number, we show that a large-scale depth-invariant flow is formed, reminiscent of the condensate state observed in two-dimensional flows. We show that the large-scale circulation shares many similarities with the so-called vortex, or slow-mode, of forced rotating turbulence. Our investigations show that at a fixed rotation rate the large-scale vortex is only observed for a finite range of Rayleigh numbers, as the quasi-two-dimensional nature of the flow disappears at very high Rayleigh numbers. We observe slow vortex merging events and find a non-local inverse cascade of energy in addition to the regular direct cascade associated with fast small-scale turbulent motions. Finally, we show that cyclonic structures are dominant in the small-scale turbulent flow and this symmetry breaking persists ...

DOUBLE-DIFFUSIVE INSTABILITIES OF A SHEAR-GENERATED MAGNETIC LAYER

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.

Evolutionary graph theory revisited: when is an evolutionary process equivalent to the Moran process?

Evolution in finite populations is often modelled using the classical Moran process. Over the last 10 years, this methodology has been extended to structured populations using evolutionary graph theory. An important question in any such population is whether a rare mutant has a higher or lower chance of fixating (the fixation probability) than the Moran probability, i.e. that from the original Moran model, which represents an unstructured population. As evolutionary graph theory has developed, different ways of considering the interactions between individuals through a graph and an associated matrix of weights have been considered, as have a number of important dynamics. In this paper, we revisit the original paper on evolutionary graph theory in light of these extensions to consider these developments in an integrated way. In particular, we find general criteria for when an evolutionary graph with general weights satisfies the Moran probability for the set of six common evolutionary dynamics.

Long-term nonlinear behaviour of the magnetorotational instability in a localized model of an accretion disc

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.

Magnetic buoyancy instabilities in the presence of magnetic flux pumping at the base of the solar convection zone

We perform idealized numerical simulations of magnetic buoyancy instabilities in three dimensions, solving the equations of compressible magnetohydrodynamics in a model of the solar tachocline. In particular, we study the effects of including a highly simplified model of magnetic flux pumping in an upper layer (‘the convection zone’) on magnetic buoyancy instabilities in a lower layer (‘the upper parts of the radiative interior – including the tachocline’), to study these competing flux transport mechanisms at the base of the convection zone. The results of the inclusion of this effect in numerical simulations of the buoyancy instability of both a preconceived magnetic slab and a shear-generated magnetic layer are presented. In the former, we find that if we are in the regime that the downward pumping velocity is comparable with the Alfven speed of the magnetic layer, magnetic 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 localized 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 (cf. 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 up into the convection zone from the radiative interior.

Interactions between magnetohydrodynamic shear instabilities and convective flows in the solar interior

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.

Magnetic fields in astrophysical objects

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 overlooked on the grounds that the large-scale magnetic field was believed to be too weak to play an 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 emphasize 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.

The interaction of multiple convection zones in A-type stars

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 convection 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 affects the interaction. We go on to discuss the effect of varying the stiffness of the stable layer that lies between the unstable regions. We show that in A-type stars 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.

How can large-scale twisted magnetic structures naturally emerge from buoyancy instabilities?

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.

Three-Layer Magnetoconvection

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