Source author record

James M Stone

James M Stone appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

3works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

3 published item(s)

preprint2011arXiv

Magnetothermal and magnetorotational instabilities in hot accretion flows

In a hot, dilute, magnetized accretion flow, the electron mean-free path can be much greater than the Larmor radius, thus thermal conduction is anisotropic and along magnetic field lines. In this case, if the temperature decreases outward, the flow may be subject to a buoyancy instability (the magnetothermal instability, or MTI). The MTI amplifies the magnetic field, and aligns field lines with the radial direction. If the accretion flow is differentially rotating, the magnetorotational instability (MRI) may also be present. Using two-dimensional, time-dependent magnetohydrodynamic simulations, we investigate the interaction between these two instabilities. We use global simulations that span over two orders of magnitude in radius, centered on the region around the Bondi radius where the infall time of gas is longer than the growth time of both the MTI and MRI. Significant amplification of the magnetic field is produced by both instabilities, although we find that the MTI primarily amplifies the radial component, and the MRI primarily the toroidal component, of the field, respectively. Most importantly, we find that if the MTI can amplify the magnetic energy by a factor $F_t$, and the MRI by a factor $F_r$, then when the MTI and MRI are both present, the magnetic energy can be amplified by a factor of $F_t \cdot F_r$. We therefore conclude that amplification of the magnetic energy by the MTI and MRI operates independently. We also find that the MTI contributes to the transport of angular momentum, because radial motions induced by the MTI increase the Maxwell (by amplifying the magnetic field) and Reynolds stresses. Finally, we find that thermal conduction decreases the slope of the radial temperature profile. The increased temperature near the Bondi radius decreases the mass accretion rate.

preprint2010arXiv

A Hybrid Godunov Method for Radiation Hydrodynamics

From a mathematical perspective, radiation hydrodynamics can be thought of as a system of hyperbolic balance laws with dual multiscale behavior (multiscale behavior associated with the hyperbolic wave speeds as well as multiscale behavior associated with source term relaxation). With this outlook in mind, this paper presents a hybrid Godunov method for one-dimensional radiation hydrodynamics that is uniformly well behaved from the photon free streaming (hyperbolic) limit through the weak equilibrium diffusion (parabolic) limit and to the strong equilibrium diffusion (hyperbolic) limit. Moreover, one finds that the technique preserves certain asymptotic limits. The method incorporates a backward Euler upwinding scheme for the radiation energy density and flux as well as a modified Godunov scheme for the material density, momentum density, and energy density. The backward Euler upwinding scheme is first-order accurate and uses an implicit HLLE flux function to temporally advance the radiation components according to the material flow scale. The modified Godunov scheme is second-order accurate and directly couples stiff source term effects to the hyperbolic structure of the system of balance laws. This Godunov technique is composed of a predictor step that is based on Duhamel's principle and a corrector step that is based on Picard iteration. The Godunov scheme is explicit on the material flow scale but is unsplit and fully couples matter and radiation without invoking a diffusion-type approximation for radiation hydrodynamics. This technique derives from earlier work by Miniati & Colella 2007. Numerical tests demonstrate that the method is stable, robust, and accurate across various parameter regimes.

preprint2010arXiv

The Effect of the Radial Pressure Gradient in Protoplanetary Disks on Planetesimal Formation

The streaming instability (SI) provides a promising mechanism for planetesimal formation because of its ability to concentrate solids into dense clumps. The degree of clumping strongly depends on the height-integrated solid to gas mass ratio Z in protoplanetary disks (PPDs). In this letter, we show that the magnitude of the radial pressure gradient (RPG) which drives the SI (characterized by $q\equivηv_K/c_s$, where $ηv_K$ is the reduction of Keplerian velocity due to the RPG and $c_s$ is the sound speed) also strongly affects clumping. We present local two-dimensional hybrid numerical simulations of aerodynamically coupled particles and gas in the midplane of PPDs. Magnetic fields and particle self-gravity are ignored. We explore three different RPG values appropriate for typical PPDs: $q=0.025, 0.05$ and 0.1. For each $q$ value, we consider four different particle size distributions ranging from sub millimeter to meter sizes and run simulations with solid abundance from Z=0.01 up to Z=0.07. We find that a small RPG strongly promotes particle clumping in that: 1) At fixed particle size distribution, the critical solid abundance $Z_{\rm crit}$ above which particle clumping occurs monotonically increases with $q$; 2) At fixed Z, strong clumping can occur for smaller particles when $q$ is smaller. Therefore, we expect planetesimals to form preferentially in regions of PPDs with a small RPG.