Engineering PapersSearch

Engineering topics

Vishniac, Ethan T.

Publications and source records attributed to Vishniac, Ethan T..

A Numerical Gamma-Ray Burst Simulation Using Three-Dimensional Relativistic Hydrodynamics: The Transition from Spherical to Jet-like Expansion

Utilizing 3D relativistic hydrodynamical calculations, we have examined the evolution of an expanding relativistic blob of gas intended to be representative of a jet associated with ejecta from an extremely energetic event such as a hypernova, that produces a gamma-ray burst (Aloy et al. 2000; Tan, Matzner, & McKee 2001; MacFadyen, Woosley, & Heger 2001, Zhang, Woosley, & Heger 2003, Zhang, Woosley, & MacFadyen 2003). Since these are the first such calculations applied to the blob during the time in which the afterglow radiation is produced, we have purposely kept them simple in an effort to concentrate on the most fundamental aspects of the physics. We restrict our attention to the transition from spherical to jetlike expansion that occurs during the time that the Lorentz factor becomes less than the reciprocal of the jet spreading angle. We have not yet attached specific numbers to our results. From the SRHD equations, one sees that the relevant quantities are the ratios of pressure to density, and of distance to time. If we specify either one of these two sets of numbers, the other one is also determined.

Cannizzo, John K.

The hard X-ray emission spectra from accretion columns in intermediate polars

We consider the hard (greater than 2 keV) X-ray emission from accretion columns in an intermediate polar system, GK Per, using a simple settling solution. The rate of photon emission per logarithmic energy interval can be fitted with a power law, E(exp -gamma), with gamma approximately 2.0, in agreement with observations. This index is only weakly dependent on the mass accretion rate, dot-M, for dot-M in the range of a few times 10(exp 16-18) g/s. The peak energy of the photon spectra (after photoelectric absorption) is expected to be E(sub p) approximately (5 keV) gamma(exp -1/3) (N(sub H)/10(exp 23)/sq cm)(exp 1/3) where N(sub H) is the hydrogen column density along the line of sight. The observed spectra of GK Per and possibly of V1223 Sgr suggest N(sub H) approximately 10(exp 23)/sq cm. This large N(sub H) may be due to partially ionized preshock column material. Alternatively, we also consider absorption by the cool outer parts of an accretion disk. In this case the photoelectric absorption depth in the disk is a sensitive function of inclination. For GK Per the required inclination is approximately 83 deg. For mass accretion rates larger than a critical rate of approximately 10(exp 18) g/s, X-ray emission from the column accretion is significantly affected by radiation drag. Although the mass accretion rate increases dramatically during outbursts, the observed hard (greater than 2 keV) X-ray luminosity will not rise proportionately. The slope and peak energy of the outburst spectra are only weakly affected. We conclude that the observed X-ray spectra can be explained by this simple analytic solution and that the production of hard X-rays from the accretion shock at the magnetic poles in the intermediate polars is in general agreement with the observations. However, since the X-ray emission and absorption depend on the mass accretion rate in a complicated manner, observed hard X-ray luminosities (greater than 2 keV) are not a good indicator of the mass accretion rate in the X-ray emitting regions.

Yi, Insu

Gravitational instabilities in a proto-planetary disk including the effects of magnetic fields

We investigate the gravitational instability of a thin, Keplerian protoplanetary disk including the effects of a largely azimuthal magnetic field. The model follows that of our previous work (Noh, Vishniac, & Cochran 1991) except for the inclusion of a magnetic field. The disk is assumed to consist of neutral and ionized gas and neutral dust which are coupled by gravity and friction. The growth rates and eigenfunctions are calculated numerically using nonaxisymmetric linear perturbation methods. The results show that the growth rate has a maximum at some intermediate azimuthal number m, but for each value of m it is reduced relative to the unmagnetized case. The effects of the magnetic field appear more strongly on small scales. As the strength of the equilibrium magnetic field increases the growth rates decrease, and the maximum instability occurs at a lower value on m due to the increasing magnetic pressure. The response of each component to the magnetic field is discussed using the behavior of the eigenfunctions in the radial direction. With the inclusion of the magnetic field, the effects of the ionization fraction and friction on the growth rates also appear to be important for high m modes. Increasing the ionization fraction or the friction suppresses instability, but only slightly changes the maximally unstable azimuthal scales. The enhanced growth rates due to a dust component for which thermal pressure is negligible are somewhat reduced by the inclusion of a magnetic field. The effects of different boundary conditions (reflecting and transmitting) on the growth rates are also shown.

Noh, Hyerim

An upper limit to the electrodynamic power output from a thin accretion disk around a black hole

We examine the electrodynamic power output from a thin accretion disk surrounding a Kerr black hole. We derive an upper limit for the total power. We find that this limit is only moderately sensitive to most of the relativistic correction factors. In fact, this limit is less sensitive to the change in the black hole spin than one would expect from the change in the inner radius of the accretion disk. The disk around a rapidly rotating black hole will produce more nonthermal radiation than a similar disk around a slowly rotating black hole of identical mass, but the difference is only a factor of approximately 1.7.

Park, Seok Jae

Local magnetohydrodynamic instabilities and the wave-driven dynamo in accretion disks

We consider the consequences of magnetic buoyancy and the magnetic shearing instability (MSI) on the strength and organization of the magnetic field in a thin accretion disk. We discuss a model in which the wave-driven dynamo growth rate is balanced by the dissipative effects of the MSI. As in earlier work, the net helicity is due to small advective motions driven by nonlinear interactions between internal waves. Assuming a simple model of the internal wave spectrum generated from the primary m = 1 internal waves, we find that the magnetic energy density saturates at about (H/r) exp 4/3 times the local pressure (where H is the disk thickness and r is its radius). On very small scales the shearing instability will produce an isotropic fluctuating field. For a stationary disk this is equivalent to a dimensionless 'viscosity' of about (H/r) exp 4/3. The vertical and radial diffusion coefficients will be comparable to each other. Magnetic buoyancy will be largely suppressed by the turbulence due to the MSI. We present a rough estimate of its effects and find that it removes magnetic flux from the disk at a rate comparable to that caused by turbulent diffusion.

Vishniac, Ethan T.

Disk instability and the time-dependent X-ray emission from the intermediate polar GK Persei

The correlation between the disk instability model for the 1981-1989 optical outbursts of the intermediate polar GK Per and the accompanying X-ray emission is examined, and the self-consistency of the combined optical-X-ray model is investigated. Special attention is given to the nature of the transition in the X-ray emission due to the time-dependent accretion rates in the simple column accretion model. The large variation in the efficiency of hard X-ray production is explained.

Yi, Insu

The dynamic instability of adiabatic blast waves

Adiabatic blastwaves, which have a total energy injected from the center E varies as t(sup q) and propagate through a preshock medium with a density rho(sub E) varies as r(sup -omega) are described by a family of similarity solutions. Previous work has shown that adiabatic blastwaves with increasing or constant postshock entropy behind the shock front are susceptible to an oscillatory instability, caused by the difference between the nature of the forces on the two sides of the dense shell behind the shock front. This instability sets in if the dense postshock layer is sufficiently thin. The stability of adiabatic blastwaves with a decreasing postshock entropy is considered. Such blastwaves, if they are decelerating, always have a region behind the shock front which is subject to convection. Some accelerating blastwaves also have such region, depending on the values of q, omega, and gamma where gamma is the adiabatic index. However, since the shock interface stabilizes dynamically induced perturbations, blastwaves become convectively unstable only if the convective zone is localized around the origin or a contact discontinuity far from the shock front. On the other hand, the contact discontinuity of accelerating blastwaves is subject to a strong Rayleigh-Taylor instability. The frequency spectra of the nonradial, normal modes of adiabatic blastwaves have been calculated. The results have been applied to the shocks propagating through supernovae envelopes. It is shown that the metal/He and He/H interfaces are strongly unstable against the Rayleigh-Taylor instability. This instability will induce mixing in supernovae envelopes. In addition the implications of this work for the evolution of planetary nebulae is discussed.

Ryu, Dongsu

Expansion-cooled Lyman-alpha clouds

It is shown that recent observations by Pettini et al. (1990) which indicate that low-N H I Ly-alpha forest lines have small velocity widths and that the velocity widths are positively correlated with N (H I) can be understood as the result of adiabatic cooling of expanding clouds. It is argued that expansion cooling can efficiently lower temperatures and velocity widths of diffuse ionized clouds, and that this trend of diminishing temperature and velocity width in a wide range of plausible cloud models is consistent with double-quasar data. Expansion can provide a natural explanation for the steep z-evolution of the cloud numbers.

Duncan, Robert C.

A model for the distribution of the intergalactic medium

The evolution and distribution of the intergalactic medium (IGM) in a universe dominated by cold dark matter with Omega(0) = 1 and h = 0.5 are investigated. Galaxies form and eject energy into the IGM from z about 20 up to the present, and the distribution of the IGM is dominated by large connected structures. The power spectrum and two-point correlation function of the IGM show a suppressed growth due to the energy injected from galaxies and the mass subtraction to form galaxies. The high-temperature regions of the IGM correspond to the low-density regions and the low-temperature regions correspond to the high-density regions. The temperature of the IGM increases from z = 1 to z = 0, while the prsssure decreases. The present temperature distribution shows a peak at about 10 million K. The mass fraction of the IGM with temperature below 100,000 K is negligible, indicating almost all the hydrogen is ionized.

Ryu, Dongsu

A model for the distribution of dark matter, galaxies, and the intergalactic medium in a cold dark matter-dominated universe

The spatial distribution of the cold-dark-matter (CDM) and baryonic components of CDM-dominated cosmological models are characterized, summarizing the results of recent theoretical investigations. The evolution and distribution of matter in an Einstein-de Sitter universe on length scales small enough so that the Newtonian approximation is valid is followed chronologically, assuming (1) that the galaxies, CDM, and the intergalactic medium (IGM) are coupled by gravity, (2) that galaxies form by taking mass and momentum from the IGM, and (3) that the IGM responds to the energy input from the galaxies. The results of the numerical computations are presented in extensive graphs and discussed in detail.

Ryu, Dongsu

The distribution of dark matter, galaxies, and the intergalactic medium in a cold dark matter dominated universe

The evolution and distribution of galaxies and the intergalactic medium (IGM) have been studied, along with collisionless dark matter in a Universe dominated by cold dark matter. The Einstein-deSitter universe with omega sub 0 = 1 and h = 0.5 was considered (here h = H sub 0 bar 100/kms/Mpc and H sub 0 is the present value of the Hubble constant). It is assumed that initially dark matter composes 90 pct and baryonic matter composes 10 pct of total mass, and that the primordial baryonic matter is comprised of H and He, with the abundance of He equal to 10 pct of H by number. Galaxies are allowed to form out of the IGM, if the total density and baryonic density satisfy an overdensity criterion. Subsequently, the newly formed galaxies release 10 to the 60th ergs of energy into the IGM over a period of 10 to the 8th years. Calculations have been performed with 32 to the 3rd dark matter particles and 32 to the 3rd cells in a cube with comoving side length L = 9.6/h Mpc. Dark matter particles and galaxies have been followed with an N-body code, while the IGM has been followed with a fluid code.

Ryu, Dongsu