General Relativistic Magnetohydrodynamic Simulations of Black Hole Accretion Disks

Observations are providing increasingly detailed quantitative information about the accretion flows that power such high energy systems as X-ray binaries and active galactic nuclei. Analytic models of such systems must rely on assumptions such as regular flow geometry and a simple, parameterized stress. Global numerical simulations offer a way to investigate the basic physical dynamics of accretion flows without these assumptions. For black hole accretion studies one solves the equations of general relativistic magnetohydrodynamics. Magnetic fields are of fundamental importance to the structure and evolution of accretion disks because magnetic turbulence is the source of the anomalous stress that drives accretion. We have developed a three-dimensional general relativistic magnetohydrodynamic simulation code to evolve time-dependent accretion systems self-consistently. Recent global simulations of black hole accretion disks suggest that the generic structure of the accretion flow is usefully divided into five regimes: the main disk, the inner disk, the corona, the evacuated funnel, and the funnel wall jet. The properties of each of these regions are summarized.Comment: invited review at the conference "Stellar-mass, Intermediate-mass, and Supermassive Black Holes", held in Kyoto, Japan, Octorber 28-31, 2003, to be published in Progress of Theoretical Physics Supplemen

Numerical analysis of a downsized 2-stroke uniflow engine

In order to optimize the 2-stroke uniflow engine performance on vehicle applications, numerical analysis has been introduced, 3D CFD model has been built for the optimization of intake charge organization. The scavenging process was investigated and the intake port design details were improved. Then the output data from 3D CFD calculation were applied to a 1D engine model to process the analysis on engine performance. The boost system optimization of the engine has been carried out also. Furthermore, a vehicle model was also set up to investigate the engine in-vehicle performance

The Effect of Resistivity on the Nonlinear Stage of the Magnetorotational Instability in Accretion Disks

We present three-dimensional magnetohydrodynamic simulations of the nonlinear evolution of the magnetorotational instability (MRI) with a non-zero Ohmic resistivity. The properties of the saturated state depend on the initial magnetic field configuration. In simulations with an initial uniform vertical field, the MRI is able to support angular momentum transport even for large resistivities through the quasi-periodic generation of axisymmetric radial channel solutions rather than through the maintenance of anisotropic turbulence. Simulations with zero net flux show that the angular momentum transport and the amplitude of magnetic energy after saturation are significantly reduced by finite resistivity, even at levels where the linear modes are only slightly affected. This occurs at magnetic Reynolds numbers expected in low, cool states of dwarf novae, these results suggest that finite resistivity may account for the low and high angular momentum transport rates inferred for these systems.Comment: 8 figures, accepted for publication in Ap

BOSS Ultracool Dwarfs I: Colors and Magnetic Activity of M and L dwarfs

We present the colors and activity of ultracool (M7-L8) dwarfs from the Tenth Data Release of the Sloan Digital Sky Survey (SDSS). We combine previous samples of SDSS M and L dwarfs with new data obtained from the Baryon Oscillation Sky Survey (BOSS) to produce the BOSS Ultracool Dwarf (BUD) sample of 11820 M7-L8 dwarfs. By combining SDSS data with photometry from the Two Micron All Sky Survey and the Wide-Field Infrared Sky Explorer mission, we present ultracool dwarf colors from $i-z$ to $W2-W3$ as a function of spectral type, and extend the SDSS-2MASS-WISE color locus to include ultracool dwarfs. The $i-z$, $i-J$, and $z-J$ colors provide the best indication of spectral type for M7-L3 dwarfs. We also examine ultracool dwarf chromospheric activity through the presence and strength of H$\alpha$ emission. The fraction of active dwarfs rises through the M spectral sequence until it reaches $\sim$90% at spectral type L0. The fraction of active dwarfs then declines to 50% at spectral type L5; no H$\alpha$ emission is observed in the late-L dwarfs in the BUD sample. The fraction of active L0-L5 dwarfs is much higher than previously observed. The strength of activity declines with spectral type from M7 through L3, after which the data do not show a clear trend. Using one-dimensional chromosphere models, we explore the range of filling factors and chromospheric temperature structures that are consistent with H$\alpha$ observations of M0-L7 dwarfs. M dwarf chromospheres have a similar, smoothly varying range of temperature and surface coverage while L dwarf chromospheres are cooler and have smaller filling factors.Comment: 24 pages and 13 figures, submitted to AJ. A short video describing these results can be found at https://www.youtube.com/watch?v=wwX5WkuJCU

Helical Magnetorotational Instability in Magnetized Taylor-Couette Flow

Hollerbach and Rudiger have reported a new type of magnetorotational instability (MRI) in magnetized Taylor-Couette flow in the presence of combined axial and azimuthal magnetic fields. The salient advantage of this "helical'' MRI (HMRI) is that marginal instability occurs at arbitrarily low magnetic Reynolds and Lundquist numbers, suggesting that HMRI might be easier to realize than standard MRI (axial field only). We confirm their results, calculate HMRI growth rates, and show that in the resistive limit, HMRI is a weakly destabilized inertial oscillation propagating in a unique direction along the axis. But we report other features of HMRI that make it less attractive for experiments and for resistive astrophysical disks. Growth rates are small and require large axial currents. More fundamentally, instability of highly resistive flow is peculiar to infinitely long or periodic cylinders: finite cylinders with insulating endcaps are shown to be stable in this limit. Also, keplerian rotation profiles are stable in the resistive limit regardless of axial boundary conditions. Nevertheless, the addition of toroidal field lowers thresholds for instability even in finite cylinders.Comment: 16 pages, 2 figures, 1 table, submitted to PR

Unstable magnetohydrodynamical continuous spectrum of accretion disks. A new route to magnetohydrodynamical turbulence in accretion disks

We present a detailed study of localised magnetohydrodynamical (MHD) instabilities occuring in two--dimensional magnetized accretion disks. We model axisymmetric MHD disk tori, and solve the equations governing a two--dimensional magnetized accretion disk equilibrium and linear wave modes about this equilibrium. We show the existence of novel MHD instabilities in these two--dimensional equilibria which do not occur in an accretion disk in the cylindrical limit. The disk equilibria are numerically computed by the FINESSE code. The stability of accretion disks is investigated analytically as well as numerically. We use the PHOENIX code to compute all the waves and instabilities accessible to the computed disk equilibrium. We concentrate on strongly magnetized disks and sub--Keplerian rotation in a large part of the disk. These disk equilibria show that the thermal pressure of the disk can only decrease outwards if there is a strong gravitational potential. Our theoretical stability analysis shows that convective continuum instabilities can only appear if the density contours coincide with the poloidal magnetic flux contours. Our numerical results confirm and complement this theoretical analysis. Furthermore, these results show that the influence of gravity can either be stabilizing or destabilizing on this new kind of MHD instability. In the likely case of a non--constant density, the height of the disk should exceed a threshold before this type of instability can play a role. This localised MHD instability provides an ideal, linear route to MHD turbulence in strongly magnetized accretion disk tori.Comment: 20 pages, 10 figures, accepted for publication in Astronomy & Astrophysic

Vortices in Thin, Compressible, Unmagnetized Disks

We consider the formation and evolution of vortices in a hydrodynamic shearing-sheet model. The evolution is done numerically using a version of the ZEUS code. Consistent with earlier results, an injected vorticity field evolves into a set of long-lived vortices, each of which has a radial extent comparable to the local scale height. But we also find that the resulting velocity field has a positive shear stress, . This effect appears only at high resolution. The transport, which decays with time as t^-1/2, arises primarily because the vortices drive compressive motions. This result suggests a possible mechanism for angular momentum transport in low-ionization disks, with two important caveats: a mechanism must be found to inject vorticity into the disk, and the vortices must not decay rapidly due to three-dimensional instabilities.Comment: 8 pages, 10 figures (high resolution figures available in ApJ electronic edition