Kris Beckwith

Global Simulations of Accretion Disks. I. Convergence and Comparisons with Local Models

Grid-based magnetohydrodynamic (MHD) simulations have proven invaluable for the study of astrophysical accretion disks. However, the fact that angular momentum transport in disks is mediated by MHD turbulence (with structure down to very small scales) raises the concern that the properties of the modeled accretion disks are affected by the finite numerical resolution of the simulation. By implementing an orbital advection algorithm into the Athena code in cylindrical geometry, we have performed a set of global (but unstratified) Newtonian disk simulations extending up to resolutions previously unattained. We study the convergence of these models as a function of spatial resolution and initial magnetic field geometry. The usual viscosity parameter (?) or the ratio of thermal-to-magnetic pressure (?) is found to be a poor diagnostic of convergence, whereas the average tilt angle of the magnetic field in the (r, )-plane is a very good diagnostic of convergence. We suggest that this is related to the saturation of the MHD turbulence via parasitic modes of the magnetorotational instability. Even in the case of zero-net magnetic flux, we conclude that our highest resolution simulations (with 32 zones and 64 zones per vertical scale height) have achieved convergence. Our global simulations reach resolutions comparable to those used in local, shearing-box models of MHD disk turbulence. We find that the saturation predictors derived from local simulations correspond well to the instantaneous correlations between local flux and stress found in our global simulations. However, the conservation of magnetic flux implicit in local models is not realized in our global disks. Thus, the magnetic connectivity of an accretion disk represents physics that is truly global and cannot be captured in any ab initio local model.

Turbulence in the Outer Regions of Protoplanetary Disks. I. Weak Accretion with no Vertical Magnetic Flux

We use local numerical simulations to investigate the strength and nature of magnetohydrodynamic (MHD) turbulence in the outer regions of protoplanetary disks, where ambipolar diffusion is the dominant non-ideal MHD effect. The simulations include vertical stratification and assume zero net vertical magnetic flux. We employ a super time-stepping technique to ameliorate the Courant restriction on the diffusive time step. We find that in idealized stratified simulations, with a spatially constant ambipolar Elsasser number Am, turbulence driven by the magnetorotational instability (MRI) behaves in a similar manner as in prior unstratified calculations. Turbulence dies away for Am ? 1, and becomes progressively more vigorous as ambipolar diffusion is decreased. Near-ideal MHD behavior is recovered for Am ? 103. In the intermediate regime (10 ? Am ? 103) ambipolar diffusion leads to substantial increases in both the period of the MRI dynamo cycle and the characteristic scales of magnetic field structures. To quantify the impact of ambipolar physics on disk accretion, we run simulations at 30 AU and 100 AU that include a vertical Am profile based upon far-ultraviolet (FUV) ionized disk models. These models develop a vertically layered structure analogous to the Ohmic dead zone that is present at smaller radii. We find that, although the levels of surface turbulence can be strong (and consistent with constraints on turbulent line widths at these radii), the inferred accretion rates are at least an order of magnitude smaller than those observed in T Tauri stars. This discrepancy is very likely due to the assumption of zero vertical magnetic field in our simulations and suggests that vertical magnetic fields are essential for MRI-driven accretion in the outer regions of protoplanetary disks.

Turbulence in global simulations of magnetized thin accretion discs

We use a global magnetohydrodynamic simulation of a geometrically thin accretion disc to investigate the locality and detailed structure of turbulence driven by the magnetorotational instability (MRI). The model disc has an aspect ratio H/R≃ 0.07, and is computed using a higher order Godunov magnetohydrodynamics (MHD) scheme with accurate fluxes. We focus the analysis on late times after the system has lost direct memory of its initial magnetic flux state. The disc enters a saturated turbulent state in which the fastest growing modes of the MRI are well resolved, with a relatively high efficiency of angular momentum transport 〈〈α〉〉≈ 2.5 × 10−2. The accretion stress peaks at the disc mid-plane, above and below which exists a moderately magnetized corona with patches of superthermal field. By analysing the spatial and temporal correlations of the turbulent fields, we find that the spatial structure of the magnetic and kinetic energy is moderately well localized (with correlation lengths along the major axis of 2.5H and 1.5H, respectively), and generally consistent with that expected of homogeneous incompressible turbulence. The density field, conversely, exhibits both a longer correlation length and a long correlation time, results that we ascribe to the importance of spiral density waves within the flow. Consistent with prior results, we show that the mean local stress displays a well-defined correlation with the local vertical flux, and that this relation is apparently causal (in the sense of the flux stimulating the stress) during portions of a global dynamo cycle. We argue that the observed flux–stress relation supports dynamo models in which the structure of coronal magnetic fields plays a central role in determining the dynamics of thin-disc accretion.

A Second-order Godunov Method for Multi-dimensional Relativistic Magnetohydrodynamics

We describe a new Godunov algorithm for relativistic magnetohydrodynamics (RMHD) that combines a simple, unsplit second-order accurate integrator with the constrained transport (CT) method for enforcing the solenoidal constraint on the magnetic field. A variety of approximate Riemann solvers are implemented to compute the fluxes of the conserved variables. The methods are tested with a comprehensive suite of multi-dimensional problems. These tests have helped us develop a hierarchy of correction steps that are applied when the integration algorithm predicts unphysical states due to errors in the fluxes, or errors in the inversion between conserved and primitive variables. Although used exceedingly rarely, these corrections dramatically improve the stability of the algorithm. We present preliminary results from the application of these algorithms to two problems in RMHD: the propagation of supersonic magnetized jets and the amplification of magnetic field by turbulence driven by the relativistic Kelvin-Helmholtz instability (KHI). Both of these applications reveal important differences between the results computed with Riemann solvers that adopt different approximations for the fluxes. For example, we show that the use of Riemann solvers that include both contact and rotational discontinuities can increase the strength of the magnetic field within the cocoon by a factor of 10 in simulations of RMHD jets and can increase the spectral resolution of three-dimensional RMHD turbulence driven by the KHI by a factor of two. This increase in accuracy far outweighs the associated increase in computational cost. Our RMHD scheme is publicly available as part of the Athena code.

Local simulations of instabilities in relativistic jets - I. Morphology and energetics of the current-driven instability

We present the results of a numerical investigation of current-driven instability in magnetized jets. Utilizing the well-tested, relativistic magnetohydrodynamic code Athena, we construct an ensemble of local, co-moving plasma columns in which initial radial force balance is achieved through various combinations of magnetic, pressure, and rotational forces. We then examine the resulting flow morphologies and energetics to determine the degree to which these systems become disrupted, the amount of kinetic energy amplification attained, and the non-linear saturation behaviors. Our most significant finding is that the details of initial force balance have a pronounced effect on the resulting flow morphology. Models in which the initial magnetic field is force-free deform, but do not become disrupted. Systems that achieve initial equilibrium by balancing pressure gradients and/or rotation against magnetic forces, however, tend to shred, mix, and develop turbulence. In all cases, the linear growth of current-driven instabilities is well-represented by analytic models. CDI-driven kinetic energy amplification is slower and saturates at a lower value in force-free models than in those that feature pressure gradients and/or rotation. In rotating columns, we find that magnetized regions undergoing rotational shear are driven toward equipartition between kinetic and magnetic energies. We show that these results are applicable for a large variety of physical parameters, but we caution that algorithmic decisions (such as choice of Riemann solver) can affect the evolution of these systems more than physically motivated parameters.

Recollimation boundary layers in relativistic jets

We study the collimation of relativistic hydrodynamic jets by the pressure of an ambient medium in the limit where the jet interior has lost causal contact with its surroundings. For a jet with an ultrarelativistic equation of state and external pressure that decreases as a power of spherical radius, p / r � , the jet interior will lose causal contact when η > 2. However, the outer layers of the jet gradually collimate toward the jet axis as long as η < 4, leading to the formation of a shocked boundary layer. Assuming that pressure-matching across the shock front determines the shape of the shock, we study the resulting structure of the jet in two ways: first by assuming that the pressure remains constant across the shocked boundary layer and looking for solutions to the shock jump equations, and then by constructing self-similar boundary-layer solutions that allow for a pressure gradient across the shocked layer. We demonstrate that the constant-pressure solutions can be characterized by four initial parameters that determine the jet shape and whether the shock closes to the axis. We show that self-similar solutions for the boundary layer can be constructed that exhibit a monotonic decrease in pressure across the boundary layer from the contact discontinuity to the shock front, and that the addition of this pressure gradient in our initial model generally causes the shock front to move outwards, creating a thinner boundary layer and decreasing the tendency of the shock to close. We discuss trends based on the value of the pressure power-law index η.

Quantifying energetics and dissipation in magnetohydrodynamic turbulence

We perform a suite of two- and three-dimensional magnetohydrodynamic (MHD) simulations with the Athena code of the non-driven Kelvin-Helmholtz instability in the subsonic, weak magnetic field limit. Focusing the analysis on the non-linear turbulent regime, we quantify energy transfer on a scale-by-scale basis and identify the physical mechanisms responsible for energy exchange by developing the diagnostic known as spectral energy transfer function analysis. At late times when the fluid is in a state of MHD turbulence, magnetic tension mediates the dominant mode of energy injection into the magnetic reservoir, whereby turbulent fluid motions twist and stretch the magnetic field lines. This generated magnetic energy turbulently cascades to smaller scales, while being exchanged backwards and forwards with the kinetic energy reservoir, until finally being dissipated. Incorporating explicit dissipation pushes the dissipation scale to larger scales than if the dissipation were entirely numerical. For scales larger than the dissipation scale, we show that the physics of energy transfer in decaying MHD turbulence is robust to numerical effects.

As a Matter of Force—Systematic Biases in Idealized Turbulence Simulations

Many astrophysical systems encompass very large dynamical ranges in space and time, which are not accessible by direct numerical simulations. Thus, idealized subvolumes are often used to study small-scale effects including the dynamics of turbulence. These turbulent boxes require an artificial driving in order to mimic energy injection from large-scale processes. In this Letter, we show and quantify how the autocorrelation time of the driving and its normalization systematically change properties of an isothermal compressible magnetohydrodynamic flow in the sub- and supersonic regime and affect astrophysical observations such as Faraday rotation. For example, we find that

Emergent mesoscale phenomena in magnetized accretion disc turbulence

\delta

RESISTIVITY-DRIVEN STATE CHANGES IN VERTICALLY STRATIFIED ACCRETION DISKS

-in-time forcing with a constant energy injection leads to a steeper slope in kinetic energy spectrum and less efficient small-scale dynamo action. In general, we show that shorter autocorrelation times require more power in the acceleration field, which results in more power in compressive modes that weaken the anticorrelation between density and magnetic field strength. Thus, derived observables, such as the line-of-sight magnetic field from rotation measures, are systematically biased by the driving mechanism. We argue that