skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: ON VERTICALLY GLOBAL, HORIZONTALLY LOCAL MODELS FOR ASTROPHYSICAL DISKS

Abstract

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 inmore » 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.« less

Authors:
;  [1]
  1. Niels Bohr International Academy, Niels Bohr Institute, Blegdamsvej 17, DK-2100, Copenhagen Ø (Denmark)
Publication Date:
OSTI Identifier:
22525343
Resource Type:
Journal Article
Resource Relation:
Journal Name: Astrophysical Journal; Journal Volume: 811; Journal Issue: 2; Other Information: Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
79 ASTROPHYSICS, COSMOLOGY AND ASTRONOMY; ACCRETION DISKS; ASTROPHYSICS; AXIAL SYMMETRY; BOUNDARY CONDITIONS; DENSITY; EQUILIBRIUM; INSTABILITY; MAGNETIC FIELDS; MAGNETOHYDRODYNAMICS; PERIODICITY; PLASMA; POTENTIALS; SHEAR

Citation Formats

McNally, Colin P., and Pessah, Martin E., E-mail: cmcnally@nbi.dk, E-mail: mpessah@nbi.dk. ON VERTICALLY GLOBAL, HORIZONTALLY LOCAL MODELS FOR ASTROPHYSICAL DISKS. United States: N. p., 2015. Web. doi:10.1088/0004-637X/811/2/121.
McNally, Colin P., & Pessah, Martin E., E-mail: cmcnally@nbi.dk, E-mail: mpessah@nbi.dk. ON VERTICALLY GLOBAL, HORIZONTALLY LOCAL MODELS FOR ASTROPHYSICAL DISKS. United States. doi:10.1088/0004-637X/811/2/121.
McNally, Colin P., and Pessah, Martin E., E-mail: cmcnally@nbi.dk, E-mail: mpessah@nbi.dk. 2015. "ON VERTICALLY GLOBAL, HORIZONTALLY LOCAL MODELS FOR ASTROPHYSICAL DISKS". United States. doi:10.1088/0004-637X/811/2/121.
@article{osti_22525343,
title = {ON VERTICALLY GLOBAL, HORIZONTALLY LOCAL MODELS FOR ASTROPHYSICAL DISKS},
author = {McNally, Colin P. and Pessah, Martin E., E-mail: cmcnally@nbi.dk, E-mail: mpessah@nbi.dk},
abstractNote = {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.},
doi = {10.1088/0004-637X/811/2/121},
journal = {Astrophysical Journal},
number = 2,
volume = 811,
place = {United States},
year = 2015,
month =
}
  • 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 resolutionmore » and initial magnetic field geometry. The usual viscosity parameter ({alpha}) or the ratio of thermal-to-magnetic pressure ({beta}) is found to be a poor diagnostic of convergence, whereas the average tilt angle of the magnetic field in the (r, {phi})-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.« less
  • We analyze a suite of global magnetohydrodynamic (MHD) accretion disk simulations in order to determine whether scaling laws for turbulence driven by the magnetorotational instability, discovered via local shearing-box studies, are globally robust. The simulations model geometrically thin disks with zero net magnetic flux and no explicit resistivity or viscosity. We show that the local Maxwell stress is correlated with the self-generated local vertical magnetic field in a manner that is similar to that found in local simulations. Moreover, local patches of vertical field are strong enough to stimulate and control the strength of angular momentum transport across much ofmore » the disk. We demonstrate the importance of magnetic linkages (through the low-density corona) between different regions of the disk in determining the local field, and suggest a new convergence requirement for global simulations-the vertical extent of the corona must be fully captured and resolved. Finally, we examine the temporal convergence of the average stress and show that an initial long-term secular drift in the local flux-stress relation dies away on a timescale that is consistent with turbulent mixing of the initial magnetic field.« less
  • 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 ismore » 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.« less
  • Oil/water displacement tests were conducted in preserved core material at reservoir pressure and at various temperatures ranging from room temperature to 522/sup 0/F (272/sup 0/C) to evaluate the effect of temperature on relative permeabilities. Both horizontally and vertically drilled plugs were tested to determine the influence of flow direction on relative permeabilities and on their temperature dependence. Results show that temperature influences both the endpoint saturations and the endpoint effective permeabilities. Irreducible water saturation increased with increasing temperature, whereas the dependence of residual oil saturation on temperature was less clear-cut. The residual oil saturation appeared to decrease up to anmore » optimum temperature, beyond which the trend reversed. The endpoint effective permeability to oil decreased with increasing temperature, while the endpoint effective permeability to water appeared to be independent of temperature. The endpoint effective permeability to oil was found to be from one to two orders of magnitude lower in vertically drilled plugs compared with horizontally drilled plugs. On the other hand, endpoint water permeability appeared to be similar in both directions.« less
  • An x-ray beam line BL-3A at the Photon Factory has been designed and constructed to have optical components of a collimating paraboloidal mirror, a sagittal-focusing double-crystal monochromator, and a focusing paraboloidal mirror. The evaluation at the (+,{minus},{minus}) setting of Si (111) showed that the vertical angular divergence was reduced from 13.4 to 9.5 (arcsec) by the collimating mirror. The vertical focusing with two mirrors gave the beam size of 0.2 mm (1/50). The horizontal beam size was reduced to 1/60 by sagittal focusing.