Idealized Global Models of Accretion Disks with Strong Toroidal Magnetic Fields

We present global magnetohydrodynamic (MHD) simulations of idealized accretion disks with a strong toroidal magnetic field using an equation of state that fixes the gas thermal scale height. The disk forms from the inflow of a rotating magnetized gas cloud with a toroidal magnetic field. We find that the system maintains a moderately strong mean azimuthal field in the midplane, with plasma- <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mi>β</mml:mi> <mml:mo>∼</mml:mo> <mml:mn>1</mml:mn> </mml:mrow> </mml:math> , trans-Alfv'enic fluctuations, and large accretion stresses <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mi>α</mml:mi> <mml:mo>∼</mml:mo> <mml:mn>0.1</mml:mn> </mml:mrow> </mml:math> . The azimuthal field in the disk is continuously escaping along the vertical direction but is also replenished via a local dynamo. The inflowing gas initially forms a strongly magnetized Keplerian disk with <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mi>β</mml:mi> <mml:mo>≪</mml:mo> <mml:mn>1</mml:mn> </mml:mrow> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mi>α</mml:mi> <mml:mo>≫</mml:mo> <mml:mn>1</mml:mn> </mml:mrow> </mml:math> . The disk gradually collapses from the inside out over <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mo>∼</mml:mo> <mml:mn>50</mml:mn> <mml:mo>−</mml:mo> <mml:mn>80</mml:mn> </mml:mrow> </mml:math> orbits to form a moderately magnetized disk with <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mi>β</mml:mi> <mml:mo>∼</mml:mo> <mml:mn>1</mml:mn> </mml:mrow> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mi>α</mml:mi> <mml:mo>∼</mml:mo> <mml:mn>0.1</mml:mn> </mml:mrow> </mml:math> . Radial advection of azimuthal magnetic field can maintain <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mi>β</mml:mi> <mml:mo>≲</mml:mo> <mml:mn>1</mml:mn> </mml:mrow> </mml:math> exterior to the circularization radius but not inside of it. Inclusion of a net initial vertical magnetic field can lead to an even more strongly magnetized disk midplane, consistent with previous work. When the gas thermal scale is not resolved ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mo>≲</mml:mo> <mml:mn>4</mml:mn> </mml:mrow> </mml:math> cells per thermal scale height), however, the disk remains highly magnetized with $ $. We discuss our results in the context of related shearing box simulations and other global disk simulations. The level of angular momentum transport found here is consistent with that inferred observationally in dwarf novae and X-ray transient outbursts, unlike simulations of weakly magnetized accretion disks.