Researcher profile

James Stone

James Stone contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 17 - UnverifiedVerification L1Unclaimed author
4works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

4 published item(s)

preprint2020arXiv

Transition region from turbulent to dead zone in protoplanetary disks: local shearing box simulations

The dynamical evolution of protoplanetary disks is of key interest for building a comprehensive theory of planet formation and to explain the observational properties of these objects. Using the magnetohydrodynamics code Athena++, with an isothermal shearing box setup, we study the boundary between the active and dead zone, where the accretion rate changes and mass can accumulate. We quantify how the turbulence level is affected by the presence of a non uniform ohmic resistivity in the radial-x direction that leads to a region of inhibited turbulence (or dead zone). Comparing the turbulent activity to that of ideal simulations, the turbulence inhibited area shows density fluctuations and magnetic activity at its boundaries, driven by energy injection from the active (ideal) zone boundaries. We find magnetic dissipation to be significantly stronger in the ideal regions, and the turbulence penetration through the boundary of the dead zone is determined by the value of the resistivity itself, through the ohmic dissipation process, though the thickness of the transition does not play a significant role in changing the dissipation. We investigate the 1D spectra along the shearing direction: magnetic spectra appear flat at large scales both in ideal as well as resistive simulations, though a Kolmogorov scaling over more than one decade persists in the dead zone, suggesting the turbulent cascade is determined by the hydrodynamics of the system: MRI dynamo action is inhibited where sufficiently high resistivity is present.

preprint2019arXiv

Global Radiation Magneto-hydrodynamic Simulations of Sub-Eddington Accretion Disks around Supermassive Black Holes

We use global three dimensional radiation magneto-hydrodynamic simulations to study the properties of inner regions of accretion disks around a 5\times 10^8 solar mass black hole with mass accretion rates reaching 7% and 20% of the Eddington value. This region of the disk is supported by magnetic pressure with surface density significantly smaller than the values predicted by the standard thin disk model but with a much larger disk scale height. The disks do not show any sign of thermal instability over many thermal time scales. More than half of the accretion is driven by radiation viscosity in the optically thin corona region for the lower accretion rate case, while accretion in the optically thick part of the disk is driven by the Maxwell and Reynolds stresses from MRI turbulence. Coronae with gas temperatures > 10^8 K are generated only in the inner \approx 10 gravitational radii in both simulations, being more compact in the higher accretion rate case. In contrast to the thin disk model, surface density increases with increasing mass accretion rate, which causes less dissipation in the optically thin region and a relatively weaker corona. The simulation results may explain the formation of X-ray coronae in Active Galactic Nuclei (AGNs), the compact size of such coronae, and the observed trend of optical to X-ray luminosity with Eddington ratio for many AGNs.

preprint2013arXiv

On the Thermal Stability of Radiation Dominated Accretion Disks

We study the long-term thermal stability of radiation dominated disks in which the vertical structure is determined self-consistently by the balance of heating due to dissipation of MHD turbulence driven by the magneto-rotational instability (MRI), and cooling due to radiation emitted at the photosphere. The calculations adopt the local shearing box approximation, and utilize the recently developed radiation transfer module in the Athena MHD code based on a variable Eddington tensor rather than an assumed local closure. After saturation of the MRI, in many cases the disk maintains a steady vertical structure for many thermal times. However, in every case in which the box size in the horizontal directions is at least one pressure scale height, fluctuations associated with MRI turbulence and dynamo action in the disk eventually trigger a thermal runaway which causes the disk to either expand or contract until the calculation must be terminated. During runaway, the dependence of the heating and cooling rates on total pressure satisfy the simplest criterion for classical thermal instability. We identify several physical reasons why the thermal runaway observed in our simulations differ from the standard alpha disk model, for example the advection of radiation contributes a non-negligible fraction to the vertical energy flux at the largest radiation pressure, most of the dissipation does not happen in the disk mid-plane, and the change of dissipation scale height with mid-plane pressure is slower than the change of density scale height. We discuss how and why our results differ from those published previously. Such thermal runaway behavior might have important implications for interpreting temporal variability in observed systems, but fully global simulations are required to study the saturated state before detailed predictions can be made.

preprint2012arXiv

Non-linear Evolution of Rayleigh-Taylor Instability in a Radiation Supported Atmosphere

The non-linear regime of Rayleigh-Taylor instability (RTI) in a radiation supported atmosphere, consisting of two uniform fluids with different densities, is studied numerically. We perform simulations using our recently developed numerical algorithm for multi-dimensional radiation hydrodynamics based on a variable Eddington tensor as implemented in Athena, focusing on the regime where scattering opacity greatly exceeds absorption opacity. We find that the radiation field can reduce the growth and mixing rate of RTI, but this reduction is only significant when radiation pressure significantly exceeds gas pressure. Small scale structures are also suppressed in this case. In the non-linear regime, dense fingers sink faster than rarefied bubbles can rise, leading to asymmetric structures about the interface. By comparing the calculations that use a variable Eddington tensor (VET) versus the Eddington approximation, we demonstrate that anisotropy in the radiation field can affect the non-linear development of RTI significantly. We also examine the disruption of a shell of cold gas being accelerated by strong radiation pressure, motivated by models of radiation driven outflows in ultraluminous infrared galaxies. We find that when the growth rate of RTI is smaller than acceleration time scale, the amount of gas that would be pushed away by the radiation field is reduced due to RTI.