Long-term Nonlinear Behaviour of the Magnetorotational Instability in a Localised 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.Comment: 10 pages, 10 figures, 2 tables, accepted by MNRA

The Need for Greater Support in Academic Writing for PhD Students in Mathematics and Related Subjects

Within the United Kingdom (UK), the graduate student population in mathematics departments seeking to obtain the higher degree of Doctor of Philosophy (PhD) has become increasingly diverse as a result of a number of factors. This student body faces a variety of challenges that raise questions about what provisions universities should provide in order to give these students the greatest chance of a successful completion of their PhD programme. In this article I argue that universities should increase their provision for the development of the writing skills of PhD students and especially for those in mathematics and related disciplines. This is partially motivated by the new diversity of the graduate student body but also by the fact that undergraduate degrees in mathematics in the UK provide few opportunities for students to develop their writing skills. I argue for a centralized provision of support, either at the department or university level, to move the development of a student’s writing away from their PhD supervisor, which is not ideal

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

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 Alfvén 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

Mean flow evolution of saturated forced shear flows in polytropic atmospheres

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 threedimensional 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

Evolution and characteristics of forced shear flows in polytropic atmospheres: Large and small Péclet number regimes

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éclet number regime with large Richardson numbers

The effect of temperature-dependent viscosity and thermal conductivity on the onset of compressible convection

The linear equations of thermal convection in a compressible fluid with non-constant transport coefficients are derived. The criterion for the onset of convection is established, based on linear stability analysis, for a range of different temperature-dependent profiles of thermal conductivity and viscosity. Temperature-dependent transport coefficients are shown to lead to a more complex behaviour than their constant counterparts, and modifies the stability condition of the fluid. When the Rayleigh number is defined in terms of the mid-layer physical properties and the temperature gradient at the top is held constant, increasing the temperature-dependence of thermal conductivity is found to raise the critical Rayleigh number dramatically, as the convective disturbance is then concentrated mainly at the top of the layer. In contrast, for viscosity a more subtle effect on stability is identified

Shear instabilities in a fully compressible polytropic atmosphere

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 Peclet 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

Double-diffusive magnetic buoyancy instability in a quasi-two-dimensional Cartesian geometry

Magnetic buoyancy, believed to occur in the solar tachocline, is both an important part of large-scale solar dynamo models and the picture of how sunspots are formed. Given that in the tachocline region the ratio of magnetic diffusivity to thermal diffusivity is small it is important, for both the dynamo and sunspot formation pictures, to understand magnetic buoyancy in this regime. Furthermore, the tachocline is a region of strong shear and such investigations must involve structures that become buoyant in the double-diffusive regime which are generated entirely from a shear flow. In a previous study, we have illustrated that shear-generated doublediffusive magnetic buoyancy instability is possible in the tachocline. However, this study was severely limited due to the computational requirements of running three-dimensional magnetohydrodynamic simulations over diffusive time-scales. A more comprehensive investigation is required to fully understand the double-diffusive magnetic buoyancy instability and its dependency on a number of key parameters; such an investigation requires the consideration of a reduced model. Here we consider a quasi-two-dimensional model where all gradients in the x direction are set to zero. We show how the instability is sensitive to changes in the thermal diffusivity and also show how different initial configurations of the forced shear flow affect the behaviour of the instability. Finally, we conclude that if the tachocline is thinner than currently stated then the double-diffusive magnetic buoyancy instability can more easily occur