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 55 records · Page 3

The extended atmosphere and circumstellar environment of the cool evolved star VX Sagittarii as seen by MATISSE

Context. VX Sgr is a cool, evolved, and luminous red star whose stellar parameters are difficult to determine, which affects its classification. Aims. We aim to spatially resolve the photospheric extent as well as the circumstellar environment. Methods. We used interferometric observations obtained with the MATISSE instrument in the L (3–4 μm), M (4.5–5 μm), and N (8–13 μm) bands. We reconstructed monochromatic images using the MIRA software. We used 3D radiation-hydrodynamics simulations carried out with CO5BOLD and a uniform disc model to estimate the apparent diameter and interpret the stellar surface structures. Moreover, we employed the radiative transfer codes OPTIM3D and RADMC3D to compute the spectral energy distribution for the L, M, and N bands, respectively. Results. MATISSE observations unveil, for the first time, the morphology of VX Sgr across the L, M, and N bands. The reconstructed images show a complex morphology with brighter areas whose characteristics depend on the wavelength probed. We measured the angular diameter as a function of the wavelength and showed that the photospheric extent in the L and M bands depends on the opacity through the atmosphere. In addition to this, we also concluded that the observed photospheric inhomogeneities can be interpreted as convection-related surface structures. The comparison in the N band yielded a qualitative agreement between the N-band spectrum and simple dust radiative transfer simulations. However, it is not possible to firmly conclude on the interpretation of the current data because of the difficulty in constraing the model parameters using the limited accuracy of our absolute flux calibration. Conclusions. MATISSE observations and the derived reconstructed images unveil the appearance of VX Sgr’s stellar surface and circumstellar environment across a very large spectral domain for the first time.

stars: atmospheres↗

Analysis of the early spectra and light curve of SN 1987A

Numerical modeling of supernova spectra, light curves, and hydrodynamics requires physical inputs, numerical techniques, approximations, and assumptions which must be thoroughly understood in order to study the details of supernova explosions. Here, we discuss some of these in the context of the early evolution of supernova 1987A. Gray radiation-hydrodynamics is used to calculate the bolometric light curve and the hydrodynamic evolution of the supernova. Synthetic spectra are then obtained for the resulting density and velocity structure. The spectrum calculations are performed using a special-relativistic treatment of the radiative transfer equation in the comoving frame, line blanketing by about 10(exp 5) spectral lines, and departures from local thermodynamic equilibrium (LTE) for H I, He I, Mg II, and Ca II. We find that we are able to simultaneously fit the early light curve and spectra reasonably well, using a progenitor model from Arnett (1991a), without fine-tuning the free parameters. Temperature structures and radiative equilibrium, non-LTE effects, homologous expansion, and mean opacities are discussed.

Hauschildt, Peter H.↗

Supercritical time-dependent accretion onto compact objects. I - Neutron stars

The time-dependent flow resulting from spherically symmetric, supercritical accretion onto a nonmagnetic neutron star is computed for accretion rates between 3 and 30. In addition, the fully time-dependent self-consistent set of coupled radiation-hydrodynamic equations governing such flow has been solved. Effects that have not previously been treated are taken into account, including separate energy equations for ions and electrons allowing for nonequilibrium, variable Eddington factors in the radiation transfer to account for anisotropy of the radiation field due to highly spherically extended envelopes, and convection of ions, electrons and photons in the accretion envelope. It is found that supercritical flow, for material falling from a radius of less than 10 to the 10th cm, is characterized by stable accretion with no evidence of stellar wind outflow at any phase of the evolution.

Klein, R. 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.↗

Shock-generated X-ray emission in radiatively driven winds - A model for Tau Scorpii

A one-dimensional radiation-hydrodynamics code is used to numerically investigate the structure and evolution of shocks in the winds of hot stars. Results are presented for the specific case of Tau Sco, a well-studied main-sequence B star for which there are X-ray data from the Einstien satellite's Solid State Spectrometer. A phenomenological radiative acceleration term and a mass-loss rate consistent with UV observations, are used to determine the time dependence of the temperatures within and X-ray emission from an isolated shock region. The driving acceleration leads to the formation of a two-component shock zone with 'forward' and 'reverse' shocks, each with their own characteristic temperature. A denser cold region forms between the two shocks, which could potentially account for the presence of narrow absorption features that are observed in the UV P Cygni profiles of many hot stars. The X-ray emission spectra from the shocks in the calculations are in good general agreement with two-temperature model fits to Einstein X-ray observations.

Macfarlane, Joseph J.↗

Atomic Physics of Shocked Plasma in Winds of Massive Stars

High resolution diffraction grating spectra of X-ray emission from massive stars obtained with Chandra and XMM-Newton have revolutionized our understanding of their powerful, radiation-driven winds. Emission line shapes and line ratios provide diagnostics on a number of key wind parameters. Modeling of resolved emission line velocity profiles allows us to derive independent constraints on stellar mass-loss rates, leading to downward revisions of a factor of a few from previous measurements. Line ratios in He-like ions strongly constrain the spatial distribution of Xray emitting plasma, confirming the expectations of radiation hydrodynamic simulations that X-ray emission begins moderately close to the stellar surface and extends throughout the wind. Some outstanding questions remain, including the possibility of large optical depths in resonance lines, which is hinted at by differences in line shapes of resonance and intercombination lines from the same ion. Resonance scattering leads to nontrivial radiative transfer effects, and modeling it allows us to place constraints on shock size, density, and velocity structure

Physics↗

Near-critical spherical accretion by neutron stars - General relativistic treatment

Time-independent spherical accretion by a neutron star is studied using general relativistic radiation hydrodynamics. Numerical integrations of the flow equations are presented. These show that when the luminosity is sufficiently close to (but below) the Eddington limit, the flow velocity increases with decreasing radius far from the neutron star, reaches a maximum at an intermediate radius, and decreases at small radii. A large fraction of the binding energy of the flow is transferred to the radiation through scattering before the flow strikes the surface of the neutron star. Following Miller's treatment of accretion at luminosites near the Eddington limit (which neglected general relativistic effects), analytic approximations for the decelerating phase of the flow's velocity profile are derived. The dependence of the solutions on the variable Eddington factor prescription chosen to close the radiation moment equations is also examined.

Park, Myeong-Gu↗

The Atmospheric Response to High Nonthermal Electron Beam Fluxes in Solar Flares. I. Modeling the Brightest NUV Footpoints in the X1 Solar Flare of 2014 March 29

The 2014 March 29 X1 solar flare (SOL20140329T17:48) produced bright continuum emission in the far- and near-ultraviolet (NUV) and highly asymmetric chromospheric emission lines, providing long-sought constraints on the heating mechanisms of the lower atmosphere in solar flares. We analyze the continuum and emission line data from the Interface Region Imaging Spectrograph (IRIS) of the brightest flaring magnetic footpoints in this flare. We compare the NUV spectra of the brightest pixels to new radiative-hydrodynamic predictions calculated with the RADYN code using constraints on a nonthermal electron beam inferred from the collisional thick-target modeling of hard X-ray data from Reuven Ramaty High Energy Solar Spectroscopic Imager. We show that the atmospheric response to a high beam flux density satisfactorily achieves the observed continuum brightness in the NUV. The NUV continuum emission in this flare is consistent with hydrogen (Balmer) recombination radiation that originates from low optical depth in a dense chromospheric condensation and from the stationary beam-heated layers just below the condensation. A model producing two flaring regions (a condensation and stationary layers) in the lower atmosphere is also consistent with the asymmetric Fe II chromospheric emission line profiles observed in the impulsive phase.

Kowalski, Adam F.↗

Observations and Simulations of the Na I D-1 Line profiles in an M-Class Solar Flare

We study the temporal evolution of the Na I D1 line profiles in the M3.9 flare SOL2014-06-11T21:03UT, using observations at high spectral resolution obtained with the Interferometric Bidimensional Spectrometer instrumentation the Dunn Solar Telescope combined with radiative hydrodynamic simulations. Our results show a significant increase in the intensities of the line core and wings during the flare. The analysis of the line profiles from the flare ribbons reveals that the Na I D1 line has a central reversal with excess emission in the blue wing (blue asymmetry).We combine RADYN and RH simulations to synthesize Na I D1 line profiles of the flaring atmosphere and find good agreement with the observations. Heating with a beam of electrons modifies the radiation field in the flaring atmosphere and excites electrons from the ground state 3s 2S to the first excited state 3p 2P, which in turn modifies the relative population of the two states. The change in temperature and the population density of the energy states make the sodium line profile revert from absorption into emission. Furthermore, the rapid changes in temperature break the pressure balance between the different layers of the lower atmosphere, generating upflow/downflow patterns. Analysis of the simulated spectra reveals that the asymmetries of the Na I D1 flare profile are produced by the velocity gradients in the lower solar atmosphere.

methods: numerical↗

Massive Stars and the Energy Balance of the Interstellar Medium: The Impact of an Isolated 60 M. Star - 1

We present results of numerical simulations carried out with a two-dimensional radiation hydrodynamics code in order to study the impact of massive stars on their surrounding interstellar medium. This first paper deals with the evolution of the circumstellar gas around an isolated 60 M. star. The interaction of the photo- ionized H II region with the stellar wind bubble forms a variety of interesting structures like shells, clouds, fingers, and spokes. These results demonstrate that complex structures found in H II regions are not necessarily relics from the time before the gas became ionized but may result from dynamical processes during the course of the H II region evolution. We have also analyzed the transfer and deposit of the stellar wind and radiation energy into the circumstellar medium until the star explodes as a supernova. Although the total mechanical wind energy supplied by the star is negligible compared to the accumulated energy of the Lyman continuum photons, the kinetic energy imparted to the circumstellar gas over the star s lifetime is 4 times higher than for a comparable windless simulation. Furthermore, the thermal energy of warm photoionized gas is lower by some 55%). Our results document the necessity to consider both ionizing radiation and stellar winds for an appropriate description of the interaction of OB stars with their circumstellar environment.

Freyer, Tim↗

Radiography in High Mass X-Ray Binaries: Micro-Structure of the Stellar Wind Through Variability of the Column Density

Context. In high mass X-ray binaries, an accreting compact object orbits a high mass star, which loses mass through a dense and in homogeneous wind. Aims. Using the compact object as an X-ray backlight, the time variability of the absorbing column density in the wind can be exploited in order to shed light on the micro-structure of the wind and obtain unbiased stellar mass-loss rates for high mass stars. Methods. We developed a simplified representation of the stellar wind where all the matter is gathered in spherical “clumps” that are radially advected away from the star. This model enables us to explore the connections between the stochastic properties of the wind and the variability of the column density for a comprehensive set of parameters related to the orbit and to the wind micro-structure, such as the size of the clumps and their individual mass. In particular, we focus on the evolution with the orbital phase of the standard deviation of the column density and of the characteristic duration of enhanced absorption episodes. Using the porosity length, we derive analytical predictions and compare them to the standard deviations and coherence time scales that were obtained. Results. We identified the favorable systems and orbital phases to determine the wind micro-structure. The coherence time scale of the column density is shown to be the self-crossing time of a single clump in front of the compact object. We thus provide a procedure to get accurate measurements of the size and of the mass of the clumps, purely based on the observable time variability of the column density. Conclusions. The coherence time scale grants direct access to the size of the clumps, while their mass can be deduced separately from the amplitude of the variability. We further show how monitoring the variability at superior conjunctions can probe the onset of the clump-forming region above the stellar photosphere. If the high column density variations in some high mass X-ray binaries are due to unaccreted clumps which are passing by the line-of-sight, this would require high mass clumps to reproduce the observed peak-to-peak amplitude and coherence time scales. These clump properties are marginally compatible with the ones derived from radiative-hydrodynamics simulations. Alternatively, the following components could contribute to the variability of the column density: larger orbital scale structures produced by a mechanism that has yet to be identified or a dense environment in the immediate vicinity of the accretor, such as an accretion disk, an outflow, or a spherical shell surrounding the magnetosphere of the accreting neutron star.

stars: mass loss↗

On the mass of the Cepheid X Cygni

Using a new non-Lagrangian hydrodynamic radiation transport code, models of the Cepheid X Cygni have been calculated and compared with new observations. Adopting an effective temperature of 5300 K, the luminosity was varied for stellar masses of 6, 8, and 9 solar masses in order to obtain a period of 16.4d. The results indicate that in order to find reasonable agreement between the light curves as calculated and observed, the mass used in the models must be nearly the evolutionary mass, not some 60% less as found in earlier work for Cepheids with periods from 7d to 10d. In addition, the new observations support the physical reality of the so-called 'artificial viscosity dip'.

Davis, C. G.↗

Grain processes in massive star formation

Observational evidence suggests that stars greater than 100 M(solar) exist in the Galaxy and Large Magellanic Cloud (LMC), however classical star formation theory predicts stellar mass limits of only approx. 60 M(solar). A protostellar accretion flow consists of inflowing gas and dust. Grains are destroyed as they are near the central protostar creating a dust shell or cocoon. Radiation pressure acting on the grain can halt the inflow of material thereby limiting the amount of mass accumulated by the protostar. We first consider rather general constraints on the initial grain to gas ratio and mass accretion rates that permit inflow. We further constrain these results by constructing a numerical model. Radiative deceleration of grains and grain destruction processes are explicitly accounted for in an iterative solution of the radiation-hydrodynamic equations. Findings seem to suggest that star formation by spherical accretion requires rather extreme preconditioning of the grain and gas environment.

Wolfire, M. G.↗

A fast operator perturbation method for the solution of the special relativistic equation of radiative transfer in spherical symmetry

A fast method for the solution of the radiative transfer equation in rapidly moving spherical media, based on an approximate Lambda-operator iteration, is described. The method uses the short characteristic method and a tridiagonal approximate Lambda-operator to achieve fast convergence. The convergence properties and the CPU time requirements of the method are discussed for the test problem of a two-level atom with background continuum absorption and Thomson scattering. Details of the actual implementation for fast vector and parallel computers are given. The method is accurate and fast enough to be incorporated in radiation-hydrodynamic calculations.

Hauschildt, P. H.↗

Non-LTE radiating acoustic shocks and Ca II K2V bright points

We present, for the first time, a self-consistent solution of the time-dependent 1D equations of non-LTE radiation hydrodynamics in solar chromospheric conditions. The vertical propagation of sinusoidal acoustic waves with periods of 30, 180, and 300 s is calculated. We find that departures from LTE and ionization recombination determine the temperature profiles of the shocks that develop. In LTE almost all the thermal energy goes into ionization, so the temperature rise is very small. In non-LTE, the finite transition rates delay the ionization to behind the shock front. The compression thus goes into thermal energy at the shock front leading to a high temperature amplitude. Further behind the shock front, the delayed ionization removes energy from the thermal pool, which reduces the temperature, producing a temperature spike. The 180 s waves reproduce the observed temporal changes in the calcium K line profiles quite well. The observed wing brightening pattern, the violet/red peak asymmetry and the observed line center behavior are all well reproduced. The short-period waves and the 5 minute period waves fail especially in reproducing the observed behavior of the wings.

Carlsson, Mats↗

An interpretation of the multipeaked structure in X-ray bursts

Some X-ray bursts from sources such as 4U/MXB 1636 - 53, 4U 1608 - 52, and possibly MXB 1727 - 34, display a multipeaked structure in the bolometric flux that is difficult to reconcile with the standard thermonuclear flash model or with the effects due to photospheric expansion and contraction. Understanding this phenomenon is important primarily because it might be due to a variable Thomson optical depth along a line of sight threading the burst-induced accretion disk corona. We explore the conditions under which these burst profiles might be observed, in light of the somewhat surprising recent results from the radiative-hydrodynamic simulation of the coronal structure. We find that double- (and triple-) peaked bursts are produced at inclination angles in the range 40-65 deg. Interestingly, the orbital information for 4U/MXB 1636 - 53 and 4U 1608 - 52 suggests i is less than 70 deg, in these sources. The orbital inclination of MXB 1728 - 34 is not yet known.

Melia, Fulvio↗

Collapse and fragmentation of molecular cloud cores. I - Moderately centrally condensed cores

3D calculations of the collapse of moderately centrally condensed molecular cloud cores with varied thermal and rotational energies are presented. The calculations are carried out using a newly developed and tested second-order accurate radiative hydrodynamics code. Because of the use of a second-order accurate numerical scheme and initial clouds that resemble both observed prolate molecular cloud cores and magnetically supported clouds at the initiation of the dynamic collapse phase, the new models provide a superior estimate of the likelihood of fragmentation as a mechanism for binary star formation.

Boss, Alan P.↗

Midplane temperatures in the solar nebula

Cosmochemical analyses of meteorites imply that maximum temperatures in the inner solar nebula were on the order of 1300 K, yet standard viscous accretion disk models predict much lower midplane temperatures (approx. 300 K at 2 AU to 3 AU) in a minimum mass nebula. A second-order accurate radiative hydrodynamics code has been used to construct models of the late-phase solar nebula appropriate for low-mass star formation (M is approximately 10(exp -6) to 10(exp -5) solar-M yr(exp -1). For a minimum mass (0.02 solar-M) nebula and a solar-mass protostar, the new models show that compressional heating due to mass accretion onto the nebula and subsequent vertical contraction of the nebula are sufficient to lead to midplane temperatures T(sub m) greater than 1400 K at 1 AU and T(sub m) greater than 1000 K at 2.5 AU.

Boss, A. P.↗