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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

The formation phase of the solar nebula

Hydrodynamical calculations of the collapse of an axisymmetric, rotating protostellar cloud, with radiation transport and without magnetic fields, are presented. The collapse is assumed to start from a centrally condensed sphere of radius 5 x 10 to the 15th cm, a mean density of 4 x 10 to the -15th g/cu cm, a total mass of 1 solar mass, and a total angular momentum of 10 to the 53rd g sq cm per sec. The numerical grid is chosen to resolve the region of disk formation between 1 and 60 AU from the center. Frequency-dependent radiative transfer calculations show how the emergent spectrum of the structure depends upon viewing angle with respect to the rotation axis and how the observed isophotal contours should depend on wavelength and viewing angle. The central part of the protostar, interior to 1 AU, is not resolved numerically but is modeled approximately. At the end of the calculation, this region is found to have a mass of 0.6 solar mass and a ratio of rotational to gravitational energy of about 0.4, sufficiently large to be unstable to nonaxisymmetric perturbations. Although the disk is gravitationally stable according to the local Toomre criterion, the nonaxisymmetric structure in the center is likely to lead to angular momentum transport.

Bodenheimer, Peter↗
NASA NTRS2011

Constraining the Properties of the Eta Carinae System via 3-D SPH Models of Space-Based Observations: The Absolute Orientation of the Binary Orbit

The extremely massive (> 90 Stellar Mass) and luminous (= 5 x 10(exp 6) Stellar Luminosity) star Eta Carinae, with its spectacular bipolar "Homunculus" nebula, comprises one of the most remarkable and intensely observed stellar systems in the Galaxy. However, many of its underlying physical parameters remain unknown. Multiwavelength variations observed to occur every 5.54 years are interpreted as being due to the collision of a massive wind from the primary star with the fast, less dense wind of a hot companion star in a highly elliptical (e approx. 0.9) orbit. Using three-dimensional (3-D) Smoothed Particle Hydrodynamics (SPH) simulations of the binary wind-wind collision, together with radiative transfer codes, we compute synthetic spectral images of [Fe III] emission line structures and compare them to existing Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS) observations. We are thus able, for the first time, to tightly constrain the absolute orientation of the binary orbit on the sky. An orbit with an inclination of approx. 40deg, an argument of periapsis omega approx. 255deg, and a projected orbital axis with a position angle of approx. 312deg east of north provides the best fit to the observations, implying that the orbital axis is closely aligned in 3-D space with the Homunculus symmetry axis, and that the companion star orbits clockwise on the sky relative to the primary.

Madura, Thomas I.↗
NASA NTRS2022

Solar Flare Heating with Turbulent Suppression of Thermal Conduction

During solar flares plasma is typically heated to very high temperatures, and the resulting redistribution of energy via thermal conduction is a primary mechanism transporting energy throughout the flaring solar atmosphere. The thermal flux is usually modeled using Spitzer’s theory, which is based on local Coulomb collisions between the electrons carrying the thermal flux and those in the background. However, often during flares, temperature gradients become sufficiently steep that the collisional mean free path exceeds the temperature gradient scale size, so that thermal conduction becomes inherently non-local. Further, turbulent angular scattering, which is detectable in nonthermal widths of atomic emission lines, can also act to increase the collision frequency and so suppress the heat flux. Recent work by Emslie & Bian (2018) extended Spitzer’s theory of thermal conduction to account for both non-locality and turbulent suppression. We have implemented their theoretical expression for the heat flux (which is a convolution of the Spitzer flux with a kernel function) into the RADYN flare-modeling code and performed a parameter study to understand how the resulting changes in thermal conduction affect flare dynamics and hence the radiation produced. We find that models with reduced heat fluxes predict slower bulk flows, less intense line emission, and longer cooling times. By comparing features of atomic emission lines predicted by the models with Doppler velocities and nonthermal line widths deduced from a particular flare observation, we find that models with suppression factors between 0.3 to 0.5 relative to the Spitzer value best reproduce observed Doppler velocities across emission lines forming over a wide range of temperatures. Interestingly, the model that best matches observed nonthermal line widths has a kappa-type velocity distribution function

Solar flares↗
NASA NTRS1999

Reflight of the Solid Surface Combustion Experiment: Flame Radiation Near Extinction

In flame spreading in quiescent and low-velocity opposing flows, effects of surface reradiation and flame radiation are important in establishing the spread rate and whether the flame, once ignited, survives to steady spread or extinguishes after a time long compared to the ignition event. A reflight of the Solid Surface Combustion Experiment (SSCE), supported by modelling, demonstrates that for thick, flat fuels, the ultimate fate of the flame is extinction rather than steady spread. A mismatch between the thermal scale in the gas, driven by radiation, and the species diffusion scale, driven by mass diffusion, develops such that the high temperature regions of the flame are ultimately located in a region to which oxygen cannot be supplied at a sufficient rate to sustain reaction, and extinction occurs. Results of the experiment conducted on Space Shuttle mission STS 85 on 9 August 1997 are reviewed. For the flat surface geometry, while the hydrodynamic phenomena associated with opposed-flow flame spread may be treated two dimensionally, the radiative effects are three dimensional, and so modelling the radiative processes, with the mismatch in dimensionality, is difficult. The cylindrical geometry at least one long compared to the radius, provides a configuration in which the radiative processes for spread in the axial direction are two dimensional, thus simplifying the modelling. The cylindrical geometry allows for the development of more sophisticated radiative models without the complication of dimensionality concerns, e.g., discrete transfer, which is discussed in detail by Bundy (1998). Additionally, the cylindrical geometry results in a "focussing" of the heat transfer to the surface and may allow for steady spread for radii that for thick fuels of the same half-thickness there is no steady spread.

Altenkirch, R. A.↗
NASA NTRS1991

Three-dimensional compressible hydrodynamic convection in the sun and stars

A pseudospectral code is used to obtain numerical solutions of 3D compressible hydrodynamic convection in a stratified medium where the boundaries are open and where radiative losses are strongest at the top. Results are presented for a series of three simulations in progressively deeper boxes. The code does not make the anelastic approximation, but follows sound waves explicitly. Spatial organization is found in the form of vertical 'stacking' of smaller structures on top of larger ones, with each structure being roughly Hp to a vertical extent. Compressibility effects are apparent in the density: the snapshots indicate that, on a horizontal plane, the density at certain points may be of an order of 10 percent larger or smaller than elsewhere on that plane.

Hossain, Murshed↗
NASA NTRS1994

Two-fluid simulations of galaxy formation

We investigate the formation of galaxies and larger structure with a simulation modeling two gravitationally coupled fluids representing dark matter and baryons. The baryon gas dynamics are calculated with a smoothed particle hydrodynamics (SPH) method, and the physics modeled includes thermal pressure, shock heating, and radiative cooling. We simulate a 16 Mpc periodic cube with 64(exp 3) particles in each fluid and 10% baryon mass fraction. We confirm, for the first time experimentally, disk formation as a natural consequence of hierarchical clustering in a large-scale cosmological environment. The majority of isolated galaxies exhibit centrifugally supported disks. A power-law relation between cold baryonic mass and maximum rotation velocity is found, M varies as nu(sub rot)(exp alpha) with alpha = 2.5 after correcting for differential numerical resolution. Both the spatial and velocity distributions of the simulated galaxies are biased with respect to the dark matter. A counts-in-cells analysis indicates that an unphysical degree of merging in the central cluster is likely responsible for the antibias signal in the correlation function. A robust, scale-dependent velocity bias is measured. The ratio of galaxy to dark matter pairwise velocity dispersions on a scale of 1 Mpc is 0.7. The amplitude is only mildly dependent on redshift or mass cutoff and scales with separation as r(exp 0.2). The degree to which these results depend on numerical parameters is discussed. Mass resolution plays a key role in controlling the resulting fraction of cold, dense baryons. The mass fraction associated with galaxies decreases by a factor of approximately greater than 3 when the mass per particle is increased by a factor 8. Photoionization and energy input from supernova will have to be included to determine more carefully the fraction of highly dissipated material and the characteristics of the stellar component of galaxies.

Evrard, August E.↗
NASA NTRS2001

The Nature of the Vela Supernova Remnant as Revealed by O VI and C IV Absorption

Highly ionized gas, in particular C IV and O VI, is produced in the interstellar medium in regions with hot (T approx. 10(exp 6) K) X-ray emitting gas and at the boundaries where hot gas and cooler (T approx. 10(exp 4) K) gas interact. Supernova remnant shocks produce most of the hot gas in the ISM and, if they are in the correct range of speeds, should produce observable quantities of C IV and O VI absorption. In turn, the column densities of these ions are potentially powerful diagnostics of the shock speed and interstellar environment in which the SNR is evolving. With the advent of FUSE, the power of this diagnostic technique is now available. We have FUSE data toward 8 stars behind the Vela SNR, and have developed a data reduction and analysis method that produces reasonably reliable O VI column densities, in spite of the complexities of the FUSE spectra in this region. In order to gain insight into the observational results, the Vela SNR evolution was modelled using Piecewise Parabolic Method numerical hydrodynamics code. The code is 1-D and incorporates non-equilibrium ionization, radiative cooling, thermal conduction and magnetic pressure.

Lines, Nichols J.↗
NASA NTRS2011

X-Ray Spectra from MHD Simulations of Accreting Black Holes

We present new global calculations of X-ray spectra from fully relativistic magneto-hydrodynamic (MHO) simulations of black hole (BH) accretion disks. With a self consistent radiative transfer code including Compton scattering and returning radiation, we can reproduce the predominant spectral features seen in decades of X-ray observations of stellar-mass BHs: a broad thermal peak around 1 keV, power-law continuum up to >100 keV, and a relativistically broadened iron fluorescent line. By varying the mass accretion rate, different spectral states naturally emerge: thermal-dominant, steep power-law, and low/hard. In addition to the spectral features, we briefly discuss applications to X-ray timing and polarization.

Schnittman, Jeremy D.↗
NASA NTRS1994

Three-dimensional hydrodynamic Bondi-Hoyle accretion. 1: Code validation and stationary accretors

We investigate the hydrodynamics of three-dimensional classical Bondi-Hoyle accretion. Totally absorbing stationary spheres of varying sizes (from 10.0 down to 0.02 Bondi radii) accrete matter from a homogeneous and slightly perturbed medium, which is taken to be an ideal gas (gamma = 5/3 or 1.2). To accommodate the long-range gravitational forces, the extent of the computational volume is typically a factor of 100 larger than the radius of the accretor. We compare the numerical mass accretion rates with the theoretical predictions of Bondi, to assess the validity of the code. The hydrodynamics is modeled by the piecewise parabolic method. No energy sources (nuclear burning) or sinks (radiation, conduction) are included. The resolution in the vicinity of the accretor is increased by multiply nesting several (6-8) grids around the stationary sphere, each finer grid being a factor of 2 smaller spatially than the next coarser grid. This allows us to include a coarse model for the surface of the accretor (vacuum sphere) on the finest grid while at the same time evolving the gas on the coarser grids. The accretion rates derived numerically are in in very good agreement (to about 10% over several orders of magnitude) with the values given by Bondi for a stationary accretor within a hydrodynamic medium. However, the equations have to be changed in order to include the finite size of the accretor (in some cases very large compared to the sonic point or even to the Bondi radius).

Ruffert, Maximilian↗
NASA NTRS2010

Mapping and Modeling the Extended Winds of the Massive Interacting Binary, Eta Carinae

The combination HST/STIS high spatial and moderate spectral resolutions have revealed the massive interacting wind structure of Eta Carinae by forbidden lines of singly and doubly ionized elements. Throughout the 5.54-year period, lines of Fe++, Ne++, Ar++, S++ and N+ reveal the interacting wind structures, near critical electron densities of 10(exp 5) to 3 x 10(exp 7)cu cm, photoionized by the hot secondary, Eta Car B, Lines of Fe+ and Ni+ trace the denser (>10(exp 7)cu cm. less-ionized (< 8 eV) primary wind of Eta Car A as it wraps around the interacting binary stars. For 5 years of the 5.54 year period, the FUV radiation from Eta Car B escapes the orbital region, ionizing the boundaries of the expanding wind structures. But for three to six months, Eta Car B plunges into the primary wind approaching to within 1 to 2 AU, leading to cutoff of FUV and X-ray fluxes. The interacting wind structure, resolved out to 0.8", drops io ionization and then rebuilds as Eta Car B emerges from the primary wind envelope. Solid Particle Hydrodynamical(SPH) models have been developed extending out to 2000 AU and adapted to include FUV radiation effects of the winds. In turn, synthetic spectroimages of selected forbidden lines have been constructed and compared to the spectroimages recorded by the HST/STIS throughout 1998.0 to 2004.3, extending across the 1998 and 2003.5 minima. By this method, we show that the orbital axis of the binary system must bc within 15 degrees of the Homunculus axis of symmetry and that periastron occurs with Eta Car B passing on the far side of Eta Car B. This result ties the current binary orbit with the bipolar ejection with intervening skirt and leads to implications that the binary system influenced the mass ejection of the l840s and the lesser ejection of the 1890s.

Gull, Ted↗
NASA NTRS2009

Combined Modeling of Acceleration, Transport, and Hydrodynamic Response in Solar Flares: The Numerical Model - 1

Acceleration and transport of high-energy particles and fluid dynamics of atmospheric plasma are interrelated aspects of solar flares, but for convenience and simplicity they were artificially separated in the past. We present here self consistently combined Fokker-Planck modeling of particles and hydrodynamic simulation of flare plasma. Energetic electrons are modeled with the Stanford unified code of acceleration, transport, and radiation, while plasma is modeled with the Naval Research Laboratory flux tube code. We calculated the collisional heating rate directly from the particle transport code, which is more accurate than those in previous studies based on approximate analytical solutions. We repeated the simulation of Mariska et al. with an injection of power law, downward-beamed electrons using the new heating rate. For this case, a -10% difference was found from their old result. We also used a more realistic spectrum of injected electrons provided by the stochastic acceleration model, which has a smooth transition from a quasi-thermal background at low energies to a non thermal tail at high energies. The inclusion of low-energy electrons results in relatively more heating in the corona (versus chromosphere) and thus a larger downward heat conduction flux. The interplay of electron heating, conduction, and radiative loss leads to stronger chromospheric evaporation than obtained in previous studies, which had a deficit in low-energy electrons due to an arbitrarily assumed low-energy cutoff. The energy and spatial distributions of energetic electrons and bremsstrahlung photons bear signatures of the changing density distribution caused by chromospheric evaporation. In particular, the density jump at the evaporation front gives rise to enhanced emission, which, in principle, can be imaged by X-ray telescopes. This model can be applied to investigate a variety of high-energy processes in solar, space, and astrophysical plasmas.

STOCHASTIC ACCELERATION MODEL↗
NASA NTRS1984

Protostellar formation in rotating interstellar clouds. IV Nonisothermal collapse

Radiative transfer in the Eddington approximation is included in a multidimensional, self-gravitational, hydrodynamical computer code. Details of the numerical solution and thermodynamic relations are given. Comparison calculations with previous spherically symmetrical models of protostellar collapse are used to validate the basic approach and the artifices which allow the explicit hydrodynamics code to follow the accretion of gas onto a quasi-equilibrium core. A series of axisymmetric models is used to investigate the importance of rotation in collapsing clouds, as the initial amount of angular momentum is lowered, with an emphasis on the possible formation of rings. Rings readily form even in the nonisothermal regime except for very low initial angular momenta; even these clouds may experience ring formation prior to reaching stellar densities. The models imply that other effects (such as gravitational torques or turbulent viscosity) may be necesary to avoid binary formation and thus result in a presolar nebula consistent with the assumptions of either Safronov or Cameron.

Boss, A. P.↗
NASA NTRS1987

Advances in magnetospheric physics by numerical simulations - A critical analysis by an outsider

This is a selective, somewhat editorialized assessment, based on a literature survey, of computer simulation as it exists today in magnetospheric physics. Both large (MHD) and small scale dynamical simulations are described and considered from the perspective of what they are trying to do and with what success. Several specific problem areas where simulations are being carried out are called for commentary: global magnetospheric structure, magnetotail and magnetopause reconnection, Kelvin-Helmholtz instability, hydrodynamic expansion of ionospheric ions, electric double layers, heavy ion heating, and auroral kilometric radiation. It is concluded that simulation is a necessary tool for understanding magnetospheric physics and that significant progress has been made in simulation development. However, results should be evaluated knowing that many factors, some real and physical, others structural, may contribute to such output. A plea is made for greater intercalibration among different simulators working in parallel areas, so that facts can be distinguished from artifacts.

Birmingham, Thomas J.↗
NASA NTRS1991

Fragmentation in stellar collapse

A three dimensional Cartesian code was developed to model molecular cloud collapse and protostellar formation. The basic Eulerian scheme employs spatially second-order accurate, finite difference methods to advance the fluid variables in time. The hydrodynamic equations governing the collapse include the effects of self-gravity, rotation, and radiative transfer. Supplementary equations include Poisson's equation which relates the gravitational potential to the density distribution and a mean intensity equation derived with the Eddington approximation for radiative transfer. The gas pressures and internal energies are determined from equations of state derived for a molecular cloud composed of hydrogen, helium, and heavier elements.

Myhill, E. M.↗
NASA NTRS2021

Neutral–neutral synthesis of organic molecules in cometary comae

Remote and in situ observations of cometary gases have revealed the presence of a wealth of complex organic molecules, including carbon chains, alcohols, imines, and the amino acid glycine. Such chemical complexity in cometary material implies that impacts by comets could have supplied reagents for prebiotic chemistry to young planetary surfaces. However, the assumption that some of the molecules observed in cometary comae at millimetre wavelengths originate from ices stored inside the nucleus has not yet been proven. In fact, the comae of moderately-active comets reach sufficient densities within a few thousand kilometres of the nucleus for an active (solar radiation-driven) photochemistry to ensue. Here, we present results from our latest chemical-hydrodynamic models incorporating an updated reaction network, and show that the commonly-observed HC3N (cyanoacetylene) and NH2CHO (formamide) molecules can be efficiently produced in cometary comae as a result of two-body, neutral–neutral, gas-phase reactions involving well-known coma species. In the presence of a near-nucleus distributed source of CN (similar to that observed by the Rosetta spacecraft at comet 67P), we find that sufficient HC3N and NH2CHO can be synthesized to match the abundances of these molecules observed previously in Oort cloud comets. The precise coma origins of these (and other) complex organic molecules can be verified through radio interferometric mapping observations, for example using the Atacama Large Millimeter/submillimeter Array.

Comets↗
NASA NTRS1993

Hydrodynamics of Relativistic Fireballs

Many models of gamma-ray bursts involve a fireball, which is an optically thick concentration of radiation energy with a high ratio of energy density to rest mass. We examine analytically and numerically the evolution of a relativistic fireball. We show that, after an early rearrangement phase, most of the matter and energy in the fireball is concentrated within a narrow shell. The shell propagates at nearly the speed of light, with a frozen radial profile, and according to a simple set of scaling laws. The spectrum of the escaping radiation is harder at early times and softer later on. Depending on the initial energy-to-mass ratio, the final outcome of a fireball is either photons with roughly the initial temperature or ultrarelativistic baryons. In the latter case, the energy could be converted back to gamma-rays via interaction with surrounding material.

Piran, Tsvi↗
NASA NTRS1994

Three-dimensional hydrodynamic Bondi-Hoyle accretion. 2: Homogeneous medium at Mach 3 with gamma = 5/3

We investigate the hydrodynamics of three-dimensional classical Bondi-Hoyle accretion. Totally absorbing spheres of varying sizes (from 10 down to 0.01 accretion radii) move at Mach 3 relative to a homogeneous and slightly perturbed medium, which is taken to be an ideal gas (gamma = 5/3). To accommodate the long-range gravitational forces, the extent of the computational volume is 32(exp 3) accretion radii. We examine the influence of numerical procedure on physical behavior. The hydrodynamics is modeled by the 'piecewise parabolic method.' No energy sources (nuclear burning) or sinks (radiation, conduction) are included. The resolution in the vicinity of the accretor is increased by multiply nesting several (5-10) grids around the sphere, each finer grid being a factor of 2 smaller in zone dimension that the next coarser grid. The largest dynamic range (ratio of size of the largest grid to size of the finest zone) is 16,384. This allows us to include a coarse model for the surface of the accretor (vacuum sphere) on the finest grid, while at the same time evolving the gas on the coarser grids. Initially (at time t = 0-10), a shock front is set up, a Mach cone develops, and the accretion column is observable. Eventually the flow becomes unstable, destroying axisymmetry. This happens approximately when the mass accretion rate reaches the values (+/- 10%) predicted by the Bondi-Hoyle accretion formula (factor of 2 included). However, our three-dimensional models do not show the highly dynamic flip-flop flow so prominent in two-dimensional calculations performed by other authors. The flow, and thus the accretion rate of all quantities, shows quasi-periodic (P approximately equals 5) cycles between quiescent and active states. The interpolation formula proposed in an accompanying paper is found to follow the collected numerical data to within approximately 30%. The specific angular momentum accreted is of the same order of magnitude as the values previously found for two-dimensional flows.

Ruffert, Maximilian↗