Viewpoint : graphene is thin, but not infinitely so

Atomically thin graphene is considered a prototypical 2D material, but high-pressure experiments now reveal the 3D nature of its mechanical properties

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

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

Mean flow instabilities of two-dimensional convection in strong magnetic fields

The interaction of magnetic fields with convection is of great importance in astrophysics. Two well-known aspects of the interaction are the tendency of convection cells to become narrow in the perpendicular direction when the imposed field is strong, and the occurrence of streaming instabilities involving horizontal shears. Previous studies have found that the latter instability mechanism operates only when the cells are narrow, and so we investigate the occurrence of the streaming instability for large imposed fields, when the cells are naturally narrow near onset. The basic cellular solution can be treated in the asymptotic limit as a nonlinear eigenvalue problem. In the limit of large imposed field, the instability occurs for asymptotically small Prandtl number. The determination of the stability boundary turns out to be surprisingly complicated. At leading order, the linear stability problem is the linearisation of the same nonlinear eigenvalue problem, and as a result, it is necessary to go to higher order to obtain a stability criterion. We establish that the flow can only be unstable to a horizontal mean flow if the Prandtl number is smaller than order , where B0 is the imposed magnetic field, and that the mean flow is concentrated in a horizontal jet of width in the middle of the layer. The result applies to stress-free or no-slip boundary conditions at the top and bottom of the layer

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.Comment: Submitted to ApJ

The Evolution of a Double Diffusive Magnetic Buoyancy Instability

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.Comment: 9 pages; 3 figures; accepted to appear in IAU symposium 27

Oscillations and secondary bifurcations in nonlinear magnetoconvection

Complicated bifurcation structures that appear in nonlinear systems governed by partial differential equations (PDEs) can be explained by studying appropriate low-order amplitude equations. We demonstrate the power of this approach by considering compressible magnetoconvection. Numerical experiments reveal a transition from a regime with a subcritical Hopf bifurcation from the static solution, to one where finite-amplitude oscillations persist although there is no Hopf bifurcation from the static solution. This transition is associated with a codimension-two bifurcation with a pair of zero eigenvalues. We show that the bifurcation pattern found for the PDEs is indeed predicted by the second-order normal form equation (with cubic nonlinearities) for a Takens-Bogdanov bifurcation with Z2 symmetry. We then extend this equation by adding quintic nonlinearities and analyse the resulting system. Its predictions provide a qualitatively accurate description of solutions of the full PDEs over a wider range of parameter values. Replacing the reflecting (Z2) lateral boundary conditions with periodic [O(2)] boundaries allows stable travelling wave and modulated wave solutions to appear; they could be described by a third-order system

Inverse cascade and symmetry breaking in rapidly-rotating Boussinesq convection

In this paper we present numerical simulations of rapidly-rotating Rayleigh-B\'enard 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 in the large-scale vortex motion