Shearing Box Research Articles

ON VERTICALLY GLOBAL, HORIZONTALLY LOCAL MODELS FOR ASTROPHYSICAL DISKS

Disks with a barotropic equilibrium structure, for which the pressure is only a function of the density, rotate on cylinders in the presence of a gravitational potential, so that the angular frequency of such a disk is independent of height. Such disks with barotropic equilibria can be approximately modeled using the shearing box framework, representing a small disk volume with height-independent angular frequency. If the disk is in baroclinic equilibrium, the angular frequency does generally depend on height, and it is thus necessary to go beyond the standard shearing box approach. In this paper, we show that given a global disk model, it is possible to develop approximate models that are local in horizontal planes without an expansion in height with shearing-periodic boundary conditions. We refer to the resulting framework as the vertically global shearing box (VGSB). These models can be non-axisymmetric for globally barotropic equilibria but should be axisymmetric for globally baroclinic equilibria. We provide explicit equations for this VGSB which can be implemented in standard magnetohydrodynamic codes by generalizing the shearing-periodic boundary conditions to allow for a height-dependent angular frequency and shear rate. We also discuss the limitations that result from the radial approximations that are needed in order to impose height-dependent shearing periodic boundary conditions. We illustrate the potential of this framework by studying a vertical shear instability and examining the modes associated with the magnetorotational instability.

Local and global aspects of the linear MRI in accretion discs

We revisit the linear MRI in a cylindrical model of an accretion disk and uncover a number of attractive results overlooked in previous treatments. In particular, we elucidate the connection between local axisymmetric modes and global modes, and show that a local channel flow corresponds to the evanescent part of a global mode. In addition, we find that the global problem reproduces the local dispersion relation without approximation, a result that helps explain the success the local analysis enjoys in predicting global growth rates. MRI channel flows are nonlinear solutions to the governing equations in the local shearing box. However, only a small subset of MRI modes share the same property in global disk models, providing further evidence that the prominence of channels in local boxes is artificial. Finally, we verify our results via direct numerical simulations with the Godunov code RAMSES.

CHARACTERIZING THE MEAN-FIELD DYNAMO IN TURBULENT ACCRETION DISKS

The formation and evolution of a wide class of astrophysical objects is governed by turbulent, magnetized accretion disks. Understanding their secular dynamics is of primary importance. Apart from enabling mass accretion via the transport of angular momentum, the turbulence affects the long-term evolution of the embedded magnetic flux, which in turn regulates the efficiency of the transport. In this paper, we take a comprehensive next step towards an effective mean-field model for turbulent astrophysical disks by systematically studying the key properties of magnetorotational turbulence in vertically-stratified, isothermal shearing boxes. This allows us to infer emergent properties of the ensuing chaotic flow as a function of the shear parameter as well as the amount of net-vertical flux. Using the test-field method, we furthermore characterize the mean-field dynamo coefficients that describe the long-term evolution of large-scale fields. We simultaneously infer the vertical shape and the spectral scale dependence of these closure coefficients, with the latter describing non-local contributions to the turbulent electromotive force. Based on this, we infer a scale-separation ratio of about ten for the large-scale dynamo. We finally determine scaling properties of the mean-field dynamo coefficients. The relevant component of the dynamo {\alpha} effect is found to scale linearly with the shear rate, as is the corresponding turbulent diffusion, {\eta}. Together, these scalings allow us to predict, in a quantitative manner, the cycle period of the well-known butterfly diagram. This lends new support to the importance of the {\alpha}{\Omega} mechanism in determining the evolution of large-scale magnetic fields in turbulent accretion disks.

Motivation and challenge to capture both large-scale and local transport in next generation accretion theory

Accretion disc theory is less developed than stellar evolution theory although a similarly mature phenomenological picture is ultimately desired. While the interplay of theory and numerical simulations has amplified community awareness of the role of magnetic fields in angular momentum transport, there remains a long term challenge to incorporate the insights gained from simulations into improving practical models for comparison with observations. What has been learned from simulations that can lead to improvements beyond SS73 in practical models? Here, we emphasize the need to incorporate the role of non-local transport more precisely. To show where large-scale transport would fit into the theoretical framework and how it is currently missing, we review why the wonderfully practical approach of Shakura & Sunyaev (Astron. Astrophys., vol. 24, 1973, pp. 337–355, SS73) is necessarily a mean field theory, and one which does not include large-scale transport. Observations of coronae and jets, combined with the interpretation of results from shearing box simulations, of the magnetorotational instability (MRI) suggest that a significant fraction of disc transport is indeed non-local. We show that the Maxwell stresses in saturation are dominated by large-scale contributions and that the physics of MRI transport is not fully captured by a viscosity. We also clarify the standard physical interpretation of the MRI as it applies to shearing boxes. Computational limitations have so far focused most attention toward local simulations, but the next generation of global simulations should help to inform improved mean field theories. Mean field accretion theory and mean field dynamo theory should in fact be unified into a single theory that predicts the time evolution of spectra and luminosity from separate disc, corona and outflow contributions. Finally, we note that any mean field theory, including that of SS73, has a finite predictive precision that needs to be quantified when comparing the predictions to observations.

Magneto-rotational instability with the effect of cosmic rays

AbstractCosmic ray (CR) is an important component of the interstellar medium. It interacts with plasma via embedded magnetic irregularities. We study the magneto-rotational instability (MRI) in the presence of CRs. We analyse the effect of CRs including their diffusion using both linear stability analysis and MHD simulations. Two models are studied. (1) In the shearing box model, uniform initial state is adopted. Linear analysis shows that the growth rate of MRI is not sensitive to the value of CR diffusion coefficient. (2) In the differentially rotating cylinder model, the initial state is taken as a constant angular momentum polytropic disk threaded by weak uniform vertical magnetic field. Linear analysis shows that the growth rate of MRI becomes larger if the CR diffusion coefficient is larger. Both linear results are confirmed by MHD simulations.

GLOBAL PROPERTIES OF FULLY CONVECTIVE ACCRETION DISKS FROM LOCAL SIMULATIONS

We present an approach to deriving global properties of accretion disks from the knowledge of local solutions derived from numerical simulations based on the shearing box approximation. The approach consists of a two-step procedure. First a local solution valid for all values of the disk height is constructed by piecing together an interior solution obtained numerically with an analytical exterior radiative solution. The matching is obtained by assuming hydrostatic balance and radiative equilibrium. Although in principle the procedure can be carried out in general, it simplifies considerably when the interior solution is fully convective. In these cases, the construction is analogous to the derivation of the Hayashi tracks for protostars. The second step consists of piecing together the local solutions at different radii to obtain a global solution. Here we use the symmetry of the solutions with respect to the defining dimensionless numbers--in a way similar to the use of homology relations in stellar structure theory--to obtain the scaling properties of the various disk quantities with radius.

Angular momentum transport and large eddy simulations in magnetorotational turbulence: the small Pm limit

Angular momentum transport in accretion discs is often believed to be due to magnetohydrodynamic turbulence mediated by the magnetorotational instability. Despite an abundant literature on the MRI, the parameters governing the saturation amplitude of the turbulence are poorly understood and the existence of an asymptotic behavior in the Ohmic diffusion regime is not clearly established. We investigate the properties of the turbulent state in the small magnetic Prandtl number limit. Since this is extremely computationally expensive, we also study the relevance and range of applicability of the most common subgrid scale models for this problem. Unstratified shearing boxes simulations are performed both in the compressible and incompressible limits, with a resolution up to 800 cells per disc scale height. The latter constitutes the largest resolution ever attained for a simulation of MRI turbulence. In the presence of a mean magnetic field threading the domain, angular momentum transport converges to a finite value in the small Pm limit. When the mean vertical field amplitude is such that {\beta}, the ratio between the thermal and magnetic pressure, equals 1000, we find {\alpha}~0.032 when Pm approaches zero. In the case of a mean toroidal field for which {\beta}=100, we find {\alpha}~0.018 in the same limit. Both implicit LES and Chollet-Lesieur closure model reproduces these results for the {\alpha} parameter and the power spectra. A reduction in computational cost of a factor at least 16 (and up to 256) is achieved when using such methods. MRI turbulence operates efficiently in the small Pm limit provided there is a mean magnetic field. Implicit LES offers a practical and efficient mean of investigation of this regime but should be used with care, particularly in the case of a vertical field. Chollet-Lesieur closure model is perfectly suited for simulations done with a spectral code.

How to form planetesimals from mm-sized chondrules and chondrule aggregates

The size distribution of asteroids and Kuiper belt objects in the solar system is difficult to reconcile with a bottom-up formation scenario due to the observed scarcity of objects smaller than $\sim$100 km in size. Instead, planetesimals appear to form top-down, with large $100-1000$ km bodies forming from the rapid gravitational collapse of dense clumps of small solid particles. In this paper we investigate the conditions under which solid particles can form dense clumps in a protoplanetary disk. We use a hydrodynamic code to model the interaction between solid particles and the gas inside a shearing box inside the disk, considering particle sizes from sub-millimeter-sized chondrules to meter-sized rocks. We find that particles down to millimeter sizes can form dense particle clouds through the run-away convergence of radial drift known as the streaming instability. We make a map of the range of conditions (strength of turbulence, particle mass-loading, disk mass, and distance to the star) which are prone to producing dense particle clumps. Finally, we estimate the distribution of collision speeds between mm-sized particles. We calculate the rate of sticking collisions and obtain a robust upper limit on the particle growth timescale of $\sim$$10^5$ years. This means that mm-sized chondrule aggregates can grow on a timescale much smaller than the disk accretion timescale ($\sim$$10^6 - 10^7$ years). Our results suggest a pathway from the mm-sized grains found in primitive meteorites to fully formed asteroids. We speculate that asteroids may form from a positive feedback loop in which coagualation leads to particle clumping driven by the streaming instability. This clumping, in turn reduces collision speeds and enhances coagulation.} Future simulations should model coagulation and the streaming instability together to explore this feedback loop further.

Numerical simulations of the magnetorotational instability in protoneutron stars – I. Influence of buoyancy

The magneto-rotational instability (MRI) is considered to be a promising mechanism to amplify the magnetic field in fast rotating protoneutron stars. In contrast to accretion disks, radial buoyancy driven by entropy and lepton fraction gradients is expected to have a dynamical role as important as rotation and shear. We investigate the poorly known impact of buoyancy on the non-linear phase of the MRI, by means of three dimensional numerical simulations of a local model in the equatorial plane of a protoneutron star. The use of the Boussinesq approximation allows us to utilise a shearing box model with clean shearing periodic boundary conditions, while taking into account the buoyancy driven by radial entropy and composition gradients. We find significantly stronger turbulence and magnetic fields in buoyantly unstable flows. On the other hand, buoyancy has only a limited impact on the strength of turbulence and magnetic field amplification for buoyantly stable flows in the presence of a realistic thermal diffusion. The properties of the turbulence are, however, significantly affected in the latter case. In particular, the toroidal components of the magnetic field and of the velocity become even more dominant with respect to the poloidal ones. Furthermore, we observed in the regime of stable buoyancy the formation of long lived coherent structures such as channel flows and zonal flows. Overall, our results support the ability of the MRI to amplify the magnetic field significantly even in stably stratified regions of protoneutron stars.

HYSTERESIS BETWEEN DISTINCT MODES OF TURBULENT DYNAMOS

Nonlinear mean-field models of the solar dynamo show long-term variability, which may be relevant to different states of activity inferred from long-term radiocarbon data. This paper is aimed to probe the dynamo hysteresis predicted by the recent mean-field models of Kitchatinov \& Olemskoy (2010) with direct numerical simulations. We perform three-dimensional simulations of large-scale dynamos in a shearing box with helically forced turbulence. As initial condition, we either take a weak random magnetic field or we start from a snapshot of an earlier simulation. Two quasi-stable states are found to coexist in a certain range of parameters close to the onset of the large-scale dynamo. The simulations converge to one of these states depending on the initial conditions. When either the fractional helicity or the magnetic Prandtl number is increased between successive runs above the critical value for onset of the dynamo, the field strength jumps to a finite value. However, when the fractional helicity or the magnetic Prandtl number is then decreased again, the field strength stays at a similar value (strong field branch) even below the original onset. We also observe intermittent decaying phases away from the strong field branch close to the point where large-scale dynamo action is just possible. The dynamo hysteresis seen previously in mean-field models is thus reproduced by 3D simulations. Its possible relation to distinct modes of solar activity such as grand minima is discussed.

GLOBAL SIMULATIONS OF PROTOPLANETARY DISKS WITH OHMIC RESISTIVITY AND AMBIPOLAR DIFFUSION

Protoplanetary disks are believed to accrete onto their central T Tauri star because of magnetic stresses. Recently published shearing box simulations indicate that Ohmic resistivity, ambipolar diffusion and the Hall effect all play important roles in disk evolution. In the presence of a vertical magnetic field, the disk remains laminar between 1-5au, and a magnetocentrifugal disk wind forms that provides an important mechanism for removing angular momentum. Questions remain, however, about the establishment of a true physical wind solution in the shearing box simulations because of the symmetries inherent in the local approximation. We present global MHD simulations of protoplanetary disks that include Ohmic resistivity and ambipolar diffusion, where the time-dependent gas-phase electron and ion fractions are computed under FUV and X-ray ionization with a simplified recombination chemistry. Our results show that the disk remains laminar, and that a physical wind solution arises naturally in global disk models. The wind is sufficiently efficient to explain the observed accretion rates. Furthermore, the ionization fraction at intermediate disk heights is large enough for magneto-rotational channel modes to grow and subsequently develop into belts of horizontal field. Depending on the ionization fraction, these can remain quasi-global, or break-up into discrete islands of coherent field polarity. The disk models we present here show a dramatic departure from our earlier models including Ohmic resistivity only. It will be important to examine how the Hall effect modifies the evolution, and to explore the influence this has on the observational appearance of such systems, and on planet formation and migration.

DUST TRANSPORT IN MRI TURBULENT DISKS: IDEAL AND NON-IDEAL MHD WITH AMBIPOLAR DIFFUSION

We study dust transport in turbulent protoplanetary disks using three-dimensional global unstratified magnetohydrodynamic (MHD) simulations including Lagrangian dust particles. The turbulence is driven by the magnetorotational instability (MRI) with either ideal or non-ideal MHD that includes ambipolar diffusion (AD). In ideal MHD simulations, the surface density evolution (except for dust that drifts fastest), turbulent diffusion, and vertical scale height of dust can all be reproduced by simple one-dimensoinal and/or analytical models. However, in AD dominated simulations which simulate protoplanetary disks beyond 10s of AU, the vertical scale height of dust is larger than previously predicted. To understand this anomaly in more detail, we carry out both unstratified and stratified local shearing box simulations with Lagrangian particles, and find that turbulence in AD dominated disks has very different properties (e.g., temporal autocorrelation functions and power spectra) than turbulence in ideal MHD disks, which leads to quite different particle diffusion efficiency. For example, MRI turbulence with AD has a longer correlation time for the vertical velocity, which causes significant vertical particle diffusion and large dust scale height. In ideal MHD the Schmidt numbers ($Sc$) for radial and vertical turbulent diffusion are $Sc_{r}\sim 1$ and $Sc_{z}\gtrsim 3$, but in the AD dominated regime both $Sc_{r}$ and $Sc_{z}$ are $\lesssim 1$. Particle concentration in pressure bumps induced by MRI turbulence has also been studied. Since non-ideal MHD effects dominate most regions in protoplanetary disks, our study suggests that modeling dust transport in turbulence driven by MRI with non-ideal MHD effects is important for understanding dust transport in realistic protoplanetary disks.

Dissipative effects on the sustainment of a magnetorotational dynamo in Keplerian shear flow

The magnetorotational (MRI) dynamo has long been considered one of the possible drivers of turbulent angular momentum transport in astrophysical accretion disks. However, various numerical results suggest that this dynamo may be difficult to excite in the astrophysically relevant regime of magnetic Prandtl number (Pm) significantly smaller than unity, for reasons currently not well understood. The aim of this article is to present the first results of an ongoing numerical investigation of the role of both linear and nonlinear dissipative effects in this problem. Combining a parametric exploration and an energy analysis of incompressible nonlinear MRI dynamo cycles representative of the transitional dynamics in large aspect ratio shearing boxes, we find that turbulent magnetic diffusion makes the excitation and sustainment of this dynamo at moderate magnetic Reynolds number (Rm) increasingly difficult for decreasing Pm. This results in an increase in the critical Rm of the dynamo for increasing kinematic Reynolds number (Re), in agreement with earlier numerical results. Given its very generic nature, we argue that turbulent magnetic diffusion could be an important determinant of MRI dynamo excitation in disks, and may also limit the efficiency of angular momentum transport by MRI turbulence in low Pm regimes.

FULLY CONVECTIVE MAGNETO-ROTATIONAL TURBULENCE IN LARGE ASPECT-RATIO SHEARING BOXES

We present a numerical study of turbulence and dynamo action in stratified shearing boxes with both finite and zero net magnetic flux. We assume that the fluid obeys the perfect gas law and has finite thermal diffusivity. The latter is chosen to be small enough so that vigorous convective states develop. The properties of these convective solutions are analyzed as the aspect ratio of the computational domain is varied and as the value of the mean field is increased. For the cases with zero net flux, we find that a well-defined converged state is obtained for large enough aspect ratios. In the converged state, the dynamo can be extremely efficient and can generate substantial toroidal flux. We identify solutions in which the toroidal field is mostly symmetric about the mid-plane and solutions in which it is mostly anti-symmetric. The symmetric solutions are found to be more efficient at transporting angular momentum and can give rise to a luminosity that is up to an order of magnitude larger than the corresponding value for the anti-symmetric states. In the cases with a finite net flux, the system appears to spend most of the time in the symmetric states.

ANALYSIS OF MAGNETOROTATIONAL INSTABILITY WITH THE EFFECT OF COSMIC-RAY DIFFUSION

We present the results obtained from linear stability analysis and 2.5-dimensional magnetohydrodynamic (MHD) simulations of the magnetorotational instability (MRI), including the effects of cosmic rays (CRs). We took into account of the CR diffusion along the magnetic field but neglect the cross-field-line diffusion. Two models are considered in this paper: shearing box model and differentially rotating cylinder model. We studied how MRI is affected by the initial CR pressure (i.e., energy) distribution. In the shearing box model, the initial state is uniform distribution. Linear analysis shows that the growth rate of MRI does not depend on the value of CR diffusion coefficient. In the differentially rotating cylinder model, the initial state is a constant angular momentum polytropic disk threaded by weak uniform vertical magnetic field. Linear analysis shows that the growth rate of MRI becomes larger if the CR diffusion coefficient is larger. Both results are confirmed by MHD simulations. The MHD simulation results show that the outward movement of matter by the growth of MRI is not impeded by the CR pressure gradient, and the centrifugal force which acts to the concentrated matter becomes larger. Consequently, the growth rate of MRI is increased. On the other hand, if the initial CR pressure is uniform, then the growth rate of the MRI barely depends on the value of the CR diffusion coefficient.

MAGNETOROTATIONAL INSTABILITY: NONMODAL GROWTH AND THE RELATIONSHIP OF GLOBAL MODES TO THE SHEARING BOX

We study the magnetorotational instability (MRI) using nonmodal stability techniques. Despite the spectral instability of many forms of the MRI, this proves to be a natural method of analysis that is well-suited to deal with the non-self-adjoint nature of the linear MRI equations. We find that the fastest growing linear MRI structures on both local and global domains can look very different to the eigenmodes, invariably resembling waves shearing with the background flow (shear waves). In addition, such structures can grow many times faster than the least stable eigenmode over long time periods, and be localized in a completely different region of space. These ideas lead -- for both axisymmetric and non-axisymmetric modes -- to a natural connection between the global MRI and the local shearing box approximation. By illustrating that the fastest growing global structure is well described by the ordinary differential equations (ODEs) governing a single shear wave, we find that the shearing box is a very sensible approximation for the linear MRI, contrary to many previous claims. Since the shear wave ODEs are most naturally understood using nonmodal analysis techniques, we conclude by analyzing local MRI growth over finite time-scales using these methods. The strong growth over a wide range of wave-numbers suggests that nonmodal linear physics could be of fundamental importance in MRI turbulence.

Quantifying the effect of turbulent magnetic diffusion on the growth rate of the magneto−rotational instability

In astrophysics, turbulent diffusion is often used in place of microphysical diffusion to avoid resolving the small scales. However, we expect this approach to break down when time and length scales of the turbulence become comparable with other relevant time and length scales in the system. Turbulent diffusion has previously been applied to the magneto-rotational instability (MRI), but no quantitative comparison of growth rates at different turbulent intensities has been performed. We investigate to what extent turbulent diffusion can be used to model the effects of small-scale turbulence on the kinematic growth rates of the MRI, and how this depends on angular velocity and magnetic field strength. We use direct numerical simulations in three-dimensional shearing boxes with periodic boundary conditions in the spanwise direction and additional random plane-wave volume forcing to drive a turbulent flow at a given length scale. We estimate the turbulent diffusivity using a mixing length formula and compare with results obtained with the test-field method. It turns out that the concept of turbulent diffusion is remarkably accurate in describing the effect of turbulence on the growth rate of the MRI. No noticeable breakdown of turbulent diffusion has been found, even when time and length scales of the turbulence become comparable with those imposed by the MRI itself. On the other hand, quenching of turbulent magnetic diffusivity by the magnetic field is found to be absent. Turbulence reduces the growth rate of the MRI in a way that is the same as microphysical magnetic diffusion.

CONVECTIVE OVERSTABILITY IN ACCRETION DISKS: THREE-DIMENSIONAL LINEAR ANALYSIS AND NONLINEAR SATURATION

Recently, Klahr & Hubbard (2014) claimed that a hydrodynamical linear overstability exists in protoplanetary disks, powered by buoyancy in the presence of thermal relaxation. We analyse this claim, confirming it through rigorous compressible linear analysis. We model the system numerically, reproducing the linear growth rate for all cases studied. We also study the saturated properties of the overstability in the shearing box, finding that the saturated state produces finite amplitude fluctuations strong enough to trigger the subcritical baroclinic instability. Saturation leads to a fast burst of enstrophy in the box, and a large-scale vortex develops in the course of the next $\approx$100 orbits. The amount of angular momentum transport achieved is of the order of $\alpha \approx 10^{-3}$, as in compressible SBI models. For the first time, a self-sustained 3D vortex is produced from linear amplitude perturbation of a quiescent base state.

On characterizing non-locality and anisotropy for the magnetorotational instability

The extent to which angular momentum transport in accretion discs is primarily local or non-local and what determines this is an important avenue of study for understanding accretion engines. Taking a step along this path, we analyse simulations of the magnetorotational instability (MRI) by calculating energy and stress power spectra in stratified isothermal shearing box simulations in several new ways. We divide our boxes in two regions, disc and corona where the disc is the MRI unstable region and corona is the magnetically dominated region. We calculate the fractional power in different quantities, including magnetic energy and Maxwell stresses and find that they are dominated by contributions from the lowest wave numbers. This is even more dramatic for the corona than the disc, suggesting that transport in the corona region is dominated by larger structures than the disc. By calculating averaged power spectra in one direction of k space at a time, we also show that the MRI turbulence is strongly anisotropic on large scales when analysed by this method, but isotropic on small scales. Although the shearing box itself is meant to represent a local section of an accretion disc, the fact that the stress and energy are dominated by the largest scales highlights that the locality is not captured within the box. This helps to quantify the intuitive importance of global simulations for addressing the question of locality of transport, for which similar analyses can be performed.

ON THE CONVERGENCE OF MAGNETOROTATIONAL TURBULENCE IN STRATIFIED ISOTHERMAL SHEARING BOXES

We consider the problem of convergence in stratified isothermal shearing boxes with zero net magnetic flux. We present results with the highest resolution to-date--up to 200 grid-point per pressure scale height--that show no clear evidence of convergence. Rather, the Maxwell stresses continue to decrease with increasing resolution. We propose some possible scenarios to explain the lack of convergence based on multi-layer dynamo systems.