Engineering PapersSearch

Engineering topics

Klein, Richard I.

Publications and source records attributed to Klein, Richard I..

Accretion disk coronae in high-luminosity systems

We present the results of self-consistent models of Compton-heated accretion disk coronae. The models are calculated using a new method for computing monochromatic radiative transfer n two dimensions. The method splits the radiation into direct and scattered components. The direct radiation is computed by calculating the optical depth along rays, while transfer of the scattered radiation is approximated by flux-limited diffusion. The resulting code agrees with more accurate treatments to within 50%, and is highly efficient, making it practical for use in large hydrodynamic simulations. The coronal models are used to confirm the results of earlier work, and to extend it to higher luminosities. In contrast to earlier work, which found the outer disks to be shadowed by the inner corona at high luminosities, we find our results to form an almost continuous extension of the models at lower luminosities. This is due to the presence of multiply scattered radiation, which acts to partially offset the loss of direct radiation from the central source. Although the analytic methods derived at lower luminosities cannot be used to derive the coronal structure for L/L(sub Edd) approx. greater than 0.1, the results of the models are amenable to semiempirical fits. We also discuss possible observational consequences of the results for coronal veiling and line fluorescence from the disk.

Murray, Stephen D.

Shock interactions with magnetized interstellar clouds. 1: Steady shocks hitting nonradiative clouds

We study the interaction of a steady, planar shock with a nonradiative, spherical, interstellar cloud threaded by a uniform magnetic field. For strong shocks, the sonic Mach number scales out, so two parameters determine the evolution: the ratio of cloud to intercloud density, and the Alfven Mach number. We focus on the case with initial field parallel to the shock velocity, though we also present one model with field perpendicular to the velocity. Even with 100 zones per cloud radius, we find that the magnetic field structure converges only at early times. However, we can draw three conclusions from our work. First, our results suggest that the inclusion of a field in equipartition with the preshock medium can prevent the complete destruction of the cloud found in the field-free case recently considered by Klein, McKee, & Colella. Second, the interaction of the shock with the cloud can amplify the magnetic field in some regions up to equipartition with the post-shock thermal pressure. In the parallel-field case, the shock preferentially amplifies the parallel component of the field, creating a 'flux rope,' a linear structure of concentrated magnetic field. The flux rope dominates the volume of amplified field, so that laminar, rather than turbulent, amplification is dominant in this case. Third, the presence of the cloud enhances the production of X-ray and synchrotron emission. The X-ray emission peaks early, during the initial passage of the shock over the cloud, while the synchrotron emission peaks later, when the flow sweeps magnetic field onto the axis between the cloud and the main shock.

Low, Mordecai-Mark Mac

On the hydrodynamic interaction of shock waves with interstellar clouds. 1: Nonradiative shocks in small clouds

The interstellar medium (ISM) is inhomogeneous, with clouds of various temperatures and densities embedded in a tenuous intercloud medium. Shocks propagating through the ISM can ablate or destroy the clouds, at the same time significantly altering the properties of the intercloud medium. This paper presents a comprehensive numerical study of the simplest case of the interaction between a shock wave and a spherical cloud, in which the shock far from the cloud is steady and planar, and in which radiative losses, thermal conduction, magnetic fields, and gravitational forces are all neglected. As a result, the problem is completely specified by two numbers: the Mach number of the shock, M, and the ratio of the density of the cloud to that of the intercloud medium, Chi. For strong shocks we show that the dependence on M scales out, so the primary independent parameter is Chi. Variations from this simple case are also considered: the potential effect of radiative losses is assessed by calculations in which the ratio of specific heats in the cloud is 1.1 instead of 5/3; the effect of the initial shape of the cloud is studied by using a cylindrical cloud instead of a spherical one; and the role of the initial shock is determined by considering the case of a cloud embedded in a wind. Local adaptive mesh refinement techniques with a second-order, two-fluid, two-dimensional Godunov hydrodynamic scheme are used to address these problems, allowing heretofore unobtainable numerical resolution. Convergence studies to be described in a subsequent paper demonstrate that about 100 zones per cloud radius are needed for accurate results; previous calculations have generally used about a third of this number. The results of the calculations are analyzed in terms of global quantities which provide an overall description of te shocked cloud: the size and shape of the cloud, the mean density, the mean pressure, the mean velocity, the velocity dispersion, and the total circulation.

Klein, Richard I.

A new scheme for multidimensional line transfer. II - ETLA method in one dimension with application to iron K-alpha lines

The 1D version of the radiation hydrodynamics computer code ALTAIR is discribed. It is an implementation of the equivalent two-level atom method, and also incorporates a variety of other iterative techniques to achieve a fast, accurate solution of the problem of coupled atomic kinetics and radiative transfer. All the techniques are described in sufficient detail to enable their use by others. The application of the program is illustrated by a set of calculations of the iron X-ray emission spectrum from intercloud hot gas that may exist in the broad-line region of active galactic nuclei. These calculations demonstrate that such spectra may be a rich load to mine for information about the hot gas in AGNs, and also that the concept of 'the iron K-alpha line' inadequately depicts the spectrum.

Castor, John I.

Isothermal, Compton-heated coronae above accretion disks

The structure of Compton-heated coronae above accretion disks is studied here by using analytic and numerical approaches are used here to determine the direct and scattered radiation reaching the base of the corona for a range of central source luminosities. It is found that the outer region of the corona is unaffected by multiple scattering in the interior, provided that the luminosity of the central source is sufficient below the Eddington limit. How attenuation and scattering by the corona affects the strength of chromospheric emission lines is determined, as is the condition for which the irradiation due to the central source exceeds the locally generated flux from the disk. Finally, it is shown that the stability analysis for irradiated accretion disks of Tuchman et al. is not substantially altered by the corona.

Ostriker, Eve C.

Iron K-shell emission from NGC 1068

The X-ray iron line emission from NGC 1068 is modeled using the new multiline, multilevel, non-LTE radiative transport code Altair and a detailed atomic model for Ne-like through-stripped iron. The X-rays passing through the ionized gas induce iron K-alpha line emission. The atomic model was constructed to describe in detail the K-shell ionization and K-alpha line emission, as well as to calculate the ionization state properly. A greater equivalent width than previously predicted is found because the observed K-alpha line is produced not only by fluorescence but also by line scattering of the continuum into the line of sight. The K-alpha equivalent width and energy are functions not only of the ionization parameter, but also of the column depth and temperature. For a likely model of NGC 1068, it is found that the iron abundance is about twice solar, but that modifications of this model may permit a smaller abundance.

Band, David L.

Detailed non-LTE calculations of the iron emission from NGC 1068

The X-ray iron line emission from NGC 1068 observed by the Ginga satellite is modeled using the new multiline, multilevel, non-LTE radiative transport code ALTAIR and a detailed atomic model for Ne-like through stripped iron. The parameter space of the obscured type 1 Seyfert nucleus model for this object is studied. The equivalent width is greater than previously predicted. It is found that detailed radiative transfer can have a significant effect on the observed line flux both for the K alpha line and for the L-shell emission. The ionization of the iron increases with temperature. Therefore the K alpha equivalent width and energy is a function not only of the ionization parameter, but also of the column depth and temperature. For a likely model of NGC 1068 it is found that the iron abundance is about twice solar, but that modifications of this model may permit a smaller abundance.

Band, David L.