Search NASA⌕ Search

SEARCH · Search NASA

Results for “Radiation Hydrodynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

Hydrodynamic effects in the atmosphere of variable stars

Numerical models of variable stars are established, using a nonlinear radiative transfer coupled hydrodynamics code. The variable Eddington method of radiative transfer is used. Comparisons are for models of W Virginis, beta Doradus, and eta Aquilae. From these models it appears that shocks are formed in the atmospheres of classical Cepheids as well as W Virginis stars. In classical Cepheids, with periods from 7 to 10 days, the bumps occurring in the light and velocity curves appear as the result of a compression wave that reflects from the star's center. At the head of the outward going compression wave, shocks form in the atmosphere. Comparisons between the hydrodynamic motions in W Virginis and classical Cepheids are made. The strong shocks in W Virginis do not penetrate into the interior as do the compression waves formed in classical Cepheids. The shocks formed in W Virginis stars cause emission lines, while in classical Cepheids the shocks are weaker.

Davis, C. G., Jr.↗

Theoretical mean colors for RR Lyrae variables

A hydrodynamically pulsating 0.6 solar mass model of a typical RR Lyrae variable was studied with a radiation transport-hydrodynamic computer program to predict theoretical T sub 3 and colors at many phases and to find the proper methods for getting mean colors and the consequent mean effective temperatures. The variable Eddington radiation approximation method was used with gray and with multifrequency absorption coefficients to represent the radiation flow in the outer optically thin layers. Comparison between observed and computed B-V colors indicate that these low Z population 2 models are reasonably accurate using King 1A composition opacities. The well known Oke, Giver, and Searle relation between B-V and T sub e reproduced. Mean colors were found by four different averaging methods. The method that gives a mean color and the mean T sub e closest to the nonpulsating model was the separate intensity means of B and V.

Davis, C. G.↗

Propagation of waves in a medium with high radiation pressure

The propagation and mutual transformation of acoustic and thermal waves are investigated in media with a high radiative pressure. The equations of hydrodynamics for matter and the radiative transfer equations in a moving medium in the Eddington approximation are used in the investigation. Model problems of waves in a homogeneous medium with an abrupt jump in opacity and in a medium of variable opacity are presented. The characteristic and the times of variability are discussed. Amplitude for the brightness fluctuations for very massive stars are discussed.

Bisnovatyy-Kogan, G. S.↗

Modeling Mg II During Solar Flares. I. Partial Frequency Redistribution, Opacity, and Coronal Irradiation

The Interface Region Imaging Spectrograph has routinely observed the flaring Mg II near-ultraviolet (NUV) spectrum, offering excellent diagnostic potential and a window into the location of energy deposition. A number of studies have forward-modeled both the general properties of these lines and specific flare observations. Generally these have forward-modeled radiation via post-processing of snapshots from hydrodynamic flare simulations through radiation transfer codes. There has, however, not been a study of how the physics included in these radiation transport codes affects the solution. A baseline setup for forward-modeling Mg II in flares is presented and contrasted with approaches that add or remove complexity. It is shown for Mg II that (1) partial frequency distribution (PRD) is still required during flare simulations despite the increased densities; (2) using full angle-dependent PRD affects the solution but takes significantly longer to process a snapshot; (3) including Mg I in non-LTE (NLTE) results in negligible differences to the Mg II lines but does affect the NUV quasi-continuum; (4) only hydrogen and Mg II need to be included in NLTE; (5) ideally the nonequilibrium hydrogen populations, with nonthermal collisional rates, should be used rather than the statistical equilibrium populations; (6) an atom consisting of only the ground state, h and k upper levels, and continuum level is insufficient to model the resonance lines; and (7) irradiation from a hot, dense flaring transition region can affect the formation of Mg II. We discuss modifications to the RH code allowing straightforward inclusion of the transition region and coronal irradiation in flares.

Kerr, Graham S.↗

Two-dimensional viscous accretion disk models. I - On meridional circulations in radiative regions

It has been suggested that thin Keplerian disks, like rotating stars, cannot attain hydrostatic and thermal equilibrium simultaneously if the angular velocity is a function of radius only. To address the problem of vertical structure and stability of thin accretion disks, we performed 2D axisymmetric calculations by solving the set of fully nonlinear hydrodynamic equations, including radiation transport. The results indicate that the variation of the rotational velocity with height is too small to sustain hydrostatic equilibrium. Instead, the variation of the radial (inflow) velocity with height suffices to establish a stationary state. In particular, we find that for low values of the viscosity parameter alpha, there is mass outflow in the central parts of the disk close to the equatorial plane, and inflow only near the surface. Only for modestly high values of alpha, the flow throughout the disk is directed inward. Thus, the flow within the disk is the result of viscous radial inflow plus a circulatory flow directed outward in the midplane and inward near the surface.

Kley, W.↗

Jet production in super-Eddington accretion disks

A two-dimensional, radiation-coupled, Newtonian hydrodynamic simulation is reported for a super-Eddington, mass accretion rate, M = 4 M(E) disk accretion flow onto a 3-solar mass pseudoblack hole. Near the disk midplane, convection cells effectively block the accretion flow, even though viscous heating maximizes there. Accretion predominantly occurs in a supersonic inflow which follows streamlines of approximately constant angular momentum. The optically thick inflow traps radiation so that 80 percent of the luminosity is absorbed by the black hole; the emergent power is sub-Eddington. An axial jet self consistently forms just outside a conical photosphere which bounds the accretion zone; radiation pressure accelerates the jet to about 10 to the 10th cm/s. The jet's mass efflux is only 0.4 percent of the total mass accretion rate.

Eggum, G. E.↗

Galactic absorption line coronae

Models are developed to determine if gaseous coronas around galaxies give rise to the absorption systems in quasar spectra. A two-dimensional hydrodynamic code with radiative cooling is used to study the dynamics and thermodynamics of a corona, and the position and rate at which clouds form. Coronal gas is found to rise from an extended gaseous disk in two models with base coronal temperatures of 10 to the 6th K and 0.50 x 10 to the 6th K. Column densities of 10 to the 19th sq cm out to a radius of 7 kcp are also achieved with a total gas mass of 2.0 x 10 to the 9th and 1.7 x 10 to the 9th solar masses, respectively. The velocity structure of the model absorption line systems may be reproduced, and clouds in pressure equilibrium with the hot coronas have densities consistent with basic observational and theoretical restrictions.

Bregman, J. N.↗

Evolution of the solar nebula. I - Nonaxisymmetric structure during nebula formation

Numerical solutions of the equations of hydrodynamics, gravitation, and radiative transfer in three spatial dimensions are used to model the formation and time evolution of the early solar nebula in order to learn whether or not gravitational torques between nonaxisymmetric structures in the solar nebula can transport angular momentum rapidly enough to produce nebula clearing on astronomically indicated (10 to the 5 to 10 to the 7 yr) time scales. The models involve solutions for the collapse of spherical clouds with assumed initial density and rotation profiles onto protosuns of variable mass. Most of the models assume uniform initial density and rotation, and have variations in the initial parameters of cloud mass, cloud rotation rate, and protosun mass which are chosen to simulate a range of possible phases of early solar nebula evolution. The models show little tendency for directly forming small numbers of giant gaseous protoplanets through gaseous gravitational instability.

Boss, Alan P.↗

Formation of the protosolar nebula

Theoretical models are discussed of the collapse of a dense molecular cloud core to form the protosolar nebula that produce the sun and the planet. The theoretical models use the equations of hydrodynamics, gravitation, and radiative transfer to follow the time evolution of a cloud collapsing under its own self-gravity. Both semi-analytical and fully numerical solutions (in two and three spatial dimensions) were calculated by several workers, One challenge is to find a set of initial conditions that will lead to the formation of a suitable protosolar nebula. Detailed results are shown for 2-D models, both with and without turbulent viscosity for redistributing angular momentum, and for 3-D models investigating the strength of gravitational torques associated with nonaxisymmetry produced during the collapse phase.

Tscharnuter, Werner M.↗

Two-dimensional viscous accretion disk models. II - On viscous overstability

Accretion disks have been shown to be susceptible to pulsational overstability when the viscosity coefficient increases upon compression. In this study, the stability of an axisymmetric accretion disk with vertical structure is considered. 2D axisymmetric calculations were performed by solving the set of fully nonlinear hydrodynamic equations, including radiation transport. The time evolution of a representative section (ring) of a thin radiative accretion disk is studied. The existence of the overstability in 2D is confirmed. It is found that the presence of periodic oscillations near the outer boundary of a real disk might account for some of the low-amplitude flickering observed in dwarf novae.

Kley, W.↗

Reionization in a cold dark matter universe: The feedback of galaxy formation on the intergalactic medium

We study the coupled evolution of the intergalactic medium (IGM) and the emerging structure in the universe in the context of the cold dark matter (CDM) model, with a special focus on the consequences of imposing reionization and the Gunn-Peterson constraint as a boundary condition on the model. We have calculated the time-varying density of the IGM by coupling our detailed, numerical calculations of the thermal and ionization balance and radiative transfer in a uniform, spatially averaged IGM of H and He, including the mean opacity of an evolving distribution of gas clumps which correspond to quasar absorption line clouds, to the linearized equations for the growth of density fluctuations in both the gaseous and dark matter components in a CDM universe. We use the linear growth equations to identify the fraction of the gas which must have collapsed out at each epoch, an approach similar in spirit to the so-called Press-Schechter formalism. We identify the IGM density with the uncollapsed baryon fraction. The collapsed fraction is postulated to be a source of energy injection into the IGM, by radiation or bulk hydrodynamical heating (e.g., via shocks) or both, at a rate which is marginally enough to satisfy the Gunn-Peterson constraint at z less than 5. Our results include the following: (1) We find that the IGM in a CDM model must have contained a substantial fraction of the total baryon density of the universe both during and after its reionization epoch. (2) As a result, our previous conclusion that the observed Quasi-Stellar Objects (QSOs) at high redshift are not sufficient to ionize the IGM enough to satisfy the Gunn-Peterson constraint is confirmed. (3) We predict a detectable He II Gunn-Peterson effect at 304(1 + z) A in the spectra of quasars at a range of redshift z greater than or approx. 3, depending on the nature of the sources of IGM reionization. (4) We find, moreover, that a CDM model with high bias parameter b (i.e., b greater than or approx. 2) cannot account for the baryon content of the universe at z approximately 3 observed in quasar absorption line gas unless Omega (sub B) significantly exceeds the maximum value allowed by big bang nucleocynthesis. (5) For a CDM model with bias parameter within the allowed range of (lower) values, the lower limit to Omega(sub B) imposed by big bang nucleosynthesis (Omega(sub B) h(sup 2) greater than or equal to 0.01) combines with our results to yield the minimum IGM density for the CDM fodel. For CDM with b = 1 (Cosmic Background Explorer (COBE) normalization), we find Omega(sub IGM)(sup min) (z approximately 4) approx. equal 0.02-0.03, and Omega(sub IGM)(sup min)(z approximately 0) approx. equal 0.005-0.03, depending upon the nature of the sources of IGM reionization. (6) In general, we find that self-consistent reionization of the IGM by the collapsed baryon fraction has a strong effect on the rate of collapse. (7) As a further example, we show that the feedback effect on the IGM of energy release by the collapsed baryon fraction may explain the slow evolution of the observed comoving QSO number density between z = 5 and z = 2, followed by the sharp decline after z = 2.

Shapiro, Paul R.↗

Physics-Based Spectra of Accretion Disks around Black Holes

The purpose of this grant was to begin the process of deriving the light output of accretion disks around black holes directly from the actual processes that inject heat into the accreting matter, rather than from guessed dependences of heating rate on physical parameters. At JHU, the effort has focussed so far on models of accretion onto "intermediate mass black holes", a possible class of black holes, examples of which may have recently been discovered in nearby galaxies. There, Krolik and his student (Yawei Hui) have computed stellar atmospheres for uniformly-heated disks around this class of black holes. Their models serve two purposes: they are the very first serious attempts to compute the spectrum from accreting black holes in this mass range; and a library of such models can be used later in this program as contrasts for those computed on the basis of real disk dynamics. The output from these local disk calculations has also been successfully coupled to a program that applies the appropriate relativistic transformations and computes photon trajectories in order to predict the spectrum received by observers located at different polar angles. The principal new result of these calculations is the discovery of potentially observable ionization edges of H-like C and O at frequencies near the peak in flux from these objects. Most of the grant money at UCSB was spent on supporting graduate student Shane Davis. In addition. some money was spent on supporting two other students: Ari Socrates (now a Hubble Fellow at Princeton), and Laura Melling. Davis spent the year constructing stellar atmosphere models of accretion disks appropriate for the high/soft (thermal) state of black hole X-ray binaries. As with AGN models published previously by our collaboration with NASA support. our models include a complete general relativistic treatment of both the disk structure and the propagation of photons from the disk to a distant observer. They also include all important continuum opacity sources, including Compton scattering and bound-free opacity from abundant metal species. The principal new result is that bound-free opacity is very significant in altering the continuum spectral shape, resulting for example in quite different "color correction factors" compared to those predicted previously. In addition, the models predict a relationship between luminosity and inner disk temperature that is, for the first time, in accord with that observed. The primary purpose of the grant was to incorporate more realistic accretion disk physics, learned largely from simulations, into such spectral models. The Davis et al. paper includes consideration of a vertical dissipation profile computed from radiation magneto-hydrodynamic simulations of MRI turbulence by N. J. Turner (2004). So long as the disk is effectively thick, such dissipation profiles do not affect the predicted spectrum significantly. (More work needs to be done on these simulations, however.) A potentially more serious issue is that MRI turbulence produces substantial inhomogeneities, as do photon bubble instabilities. These inhomogeneities can affect the spectra by enhancing the effects of absorption opacity over scattering opacity. We have done some preliminary Monte Carlo calculations to explore these effects.

Krolik, Julian H.↗

Acoustic radiation from plates excited by flow noise.

Theoretical formulation of noise radiation from thin, rectangular, flexible plates excited by turbulent wall pressure fluctuations. An extensive parametric study including surface roughness and both plate and flow parameters is presented. From the results, presented graphically, it is possible to predict the spectrum of the total radiated power provided that the power spectral density of the wall-pressure fluctuations, the flow velocity, the plate dimensions and the plate modal damping parameters are known. Along with the parametric study, an estimate of the dependence of the noise radiation on plate boundary conditions based on a consideration of the wavenumber-frequency selectivity of the plate is presented. The parametric study is based on the assumption that the wave number/frequency spectrum of the wall-pressure fluctuations can be described by a Fourier transform of a Corcos model of the cross spectral density. The results indicate that at frequencies above the hydrodynamic-coincidence frequency the radiated power is dependent on the low-wave-number spectrum of the excitation. The dependence on the high-wave-number spectrum, however, is limited to a narrow band of frequencies below the hydrodynamic-coincidence frequency.

Aupperle, F. A.↗

Calculations of the early evolution of Jupiter

The evolution of the protoplanet Jupiter is followed, using a hydrodynamic computer code with radiative energy transport. Jupiter is assumed to have formed as a subcondensation in the primitive solar nebula at a density just high enough for gravitational collapse to occur. The initial state has a density of 0.0015 nanograms per cu cm and a temperature of 43 K; the calculations are carried to an equilibrium state where the central density reaches 0.5 g per cu cm and the central temperature reaches 25,000 K. During the early part of the evolution the object contracts in quasi-hydrostatic equilibrium; later on hydrodynamic collapse occurs, induced by the dissociation of hydrogen molecules. After dissociation is complete, the planet regains hydrostatic equilibrium with a radius of a few times the present value. Further evolution beyond this point is not treated here; however the results are consistent with the existence of a high-luminosity phase shortly after the planet settles into its final quasi-static contraction.-

Bodenheimer, P.↗

Planetary geology, stellar evolution and galactic cosmology

Field studies of selected basalt flows in the Snake River Plain, Idaho, were made for comparative lunar and Mars geological investigations. Studies of basalt lava tubes were also initiated in Washington, Oregon, Hawaii, and northern California. The main effort in the stellar evolution research is toward the development of a computer code to calculate hydrodynamic flow coupled with radiative energy transport. Estimates of the rotation effects on a collapsing cloud indicate that the total angular momentum is the critical parameter. The study of Paschen and Balmer alpha lines of positronium atoms in the center of a galaxy is mentioned.

Source record↗

Evolution of rotating interstellar clouds. I - Numerical techniques

A method is described for the numerical calculation of the hydrodynamic evolution of a self-gravitating configuration in two space dimensions with assumed axial symmetry. The calculation is formulated in cylindrical coordinates with respect to a moving Eulerian grid and is solved using explicit hydrodynamics combined with implicit radiative transfer. The physics included is appropriate for calculation of the collapse of a rotating protostellar cloud. The gravitational field is obtained by means of an alternating-direction iterative technique. Numerical tests to demonstrate the correctness of the method are presented for special cases.

Black, D. C.↗

The Eddington limit and supercritical accretion. I - Time-independent calculations

Spherically symmetric, steady state accretion of an ionized hydrogen plasma onto a neutron star is considered for accretion rates which exceed a critical rate at which the Eddington limiting luminosity is produced. The coupled hydrodynamic and frequency integrated, radiative transfer equations are solved for accretion rates up to 10 times the nominal limit. Steady state solutions are presented that imply a multiplicity of different luminosity solutions for a single accretion rate in this 'supercritical' regime.

Burger, H. L.↗