Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamic codes”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

X-ray clusters in a cold dark matter + lambda universe: A direct, large-scale, high-resolution, hydrodynamic simulation

A new, three-dimensional, shock-capturing, hydrodynamic code is utilized to determine the distribution of hot gas in a cold dark matter (CDM) + lambda model universe. Periodic boundary conditions are assumed: a box with size 85/h Mpc, having cell size 0.31/h Mpc, is followed in a simulation with 270(exp 3) = 10(exp 7.3) cells. We adopt omega = 0.45, lambda = 0.55, h identically equal to H/100 km/s/Mpc = 0.6, and then, from the cosmic background explorer (COBE) and light element nucleosynthesis, sigma(sub 8) = 0.77, omega(sub b) = 0.043. We identify the X-ray emitting clusters in the simulation box, compute the luminosity function at several wavelength bands, the temperature function and estimated sizes, as well as the evolution of these quantities with redshift. This open model succeeds in matching local observations of clusters in contrast to the standard omega = 1, CDM model, which fails. It predicts an order of magnitude decline in the number density of bright (h nu = 2-10 keV) clusters from z = 0 to z = 2 in contrast to a slight increase in the number density for standard omega = 1, CDM model. This COBE-normalized CDM + lambda model produces approximately the same number of X-ray clusters having L(sub x) greater than 10(exp 43) erg/s as observed. The background radiation field at 1 keV due to clusters is approximately the observed background which, after correction for numerical effects, again indicates that the model is consistent with observations.

Cen, Renyue↗
NASA NTRS1994

Synthesized spectra of turbulent clouds

A comparison has been made of the velocity structure of non-star-forming dense molecular clouds and the output of a 512(exp 3) hydrodynamic code for compressible gas displaying features of turbulence. A remarkable similarity exists between the interstellar line profiles and the velocity profiles found by integrating, with an optically thin assumption, through subsections of the simulation. This similarity is in terms of both the structure of individual spectra, showing skewness, multiple peaks, and wings, and also the variations of the spectra from point to point. It is possible to see the origin of the various velocity substructures and how they correlate with regions of high density and velocity divergence, or regions of high vorticity, by deprojecting individual spectra of the simulation. It is found that the interstellar medium spectra are best represented by the epoch of the simulation in which the energy spectrum is essentially Kolmogorov-like and initially developed shocks have been largely converted into vorticity structures. Further, it is found that the vorticity peaks ascribed to intermittency contribute to the non-Gaussian features and wings of the line profile, but at an intensity level that is weak compared to the contribution of the bulk of the flow.

Falgarone, E.↗
NASA NTRS1994

Convection, nucleosynthesis, and core collapse

We use a piecewise parabolic method hydrodynamics code (PROMETHEUS) to study convective burning in two dimensions in an oxygen shell prior to core collapse. Significant mixing beyond convective boundaries determined by mixing-length theory brings fuel (C-12) into the convective regon, causing hot spots of nuclear burning. Plumes dominate the velocity structure. Finite perturbations arise in a region in which O-16 will be explosively burned to Ni-56 when the star explodes; the resulting instabilities and mixing are likely to distribute Ni-56 throughout the supernova envelope. Inhomogeneities in Y(sub e) may be large enough to affect core collapse and will affect explosive nucleosynthesis. The nature of convective burning is dramatically different from that assumed in one-dimensional simulations; quantitative estimates of nucleosynthetic yields, core masses, and the approach to core collapse will be affected.

Bazan, Grant↗
NASA NTRS1994

Star-forming galaxy models: Blending star formation into TREESPH

We have incorporated star-formation algorithms into a hybrid N-body/smoothed particle hydrodynamics code (TREESPH) in order to describe the star forming properties of disk galaxies over timescales of a few billion years. The models employ a Schmidt law of index n approximately 1.5 to calculate star-formation rates, and explicitly include the energy and metallicity feedback into the Interstellar Medium (ISM). Modeling the newly formed stellar population is achieved through the use of hybrid SPH/young star particles which gradually convert from gaseous to collisionless particles, avoiding the computational difficulties involved in creating new particles. The models are shown to reproduce well the star-forming properties of disk galaxies, such as the morphology, rate of star formation, and evolution of the global star-formation rate and disk gas content. As an example of the technique, we model an encounter between a disk galaxy and a small companion which gives rise to a ring galaxy reminiscent of the Cartwheel (AM 0035-35). The primary galaxy in this encounter experiences two phases of star forming activity: an initial period during the expansion of the ring, and a delayed phase as shocked material in the ring falls back into the central regions.

Mihos, J. Christopher↗
NASA NTRS1994

Hot gas in the cold dark matter scenario: X-ray clusters from a high-resolution numerical simulation

A new, three-dimensional, shock-capturing hydrodynamic code is utilized to determine the distribution of hot gas in a standard cold dark matter (CDM) model of the universe. Periodic boundary conditions are assumed: a box with size 85 h(exp -1) Mpc having cell size 0.31 h(exp -1) Mpc is followed in a simulation with 270(exp 3) = 10(exp 7.3) cells. Adopting standard parameters determined from COBE and light-element nucleosynthesis, sigma(sub 8) = 1.05, omega(sub b) = 0.06, and assuming h = 0.5, we find the X-ray-emitting clusters and compute the luminosity function at several wavelengths, the temperature distribution, and estimated sizes, as well as the evolution of these quantities with redshift. We find that most of the total X-ray emissivity in our box originates in a relatively small number of identifiable clusters which occupy approximately 10(exp -3) of the box volume. This standard CDM model, normalized to COBE, produces approximately 5 times too much emission from clusters having L(sub x) is greater than 10(exp 43) ergs/s, a not-unexpected result. If all other parameters were unchanged, we would expect adequate agreement for sigma(sub 8) = 0.6. This provides a new and independent argument for lower small-scale power than standard CDM at the 8 h(exp -1) Mpc scale. The background radiation field at 1 keV due to clusters in this model is approximately one-third of the observed background, which, after correction for numerical effects, again indicates approximately 5 times too much emission and the appropriateness of sigma(sub 8) = 0.6. If we have used the observed ratio of gas to total mass in clusters, rather than basing the mean density on light-element nucleosynthesis, then the computed luminosity of each cluster would have increased still further, by a factor of approximately 10. The number density of clusters increases to z approximately 1, but the luminosity per typical cluster decreases, with the result that evolution in the number density of bright clusters is moderate in this redshift range, showing a broad peak near z = 0.7, and then a rapid decline above redshift z = 3. Detailed computations of the luminosity functions in the range L(sub x) = 10(exp 40) - 10(exp 44) ergs/s in various energy bands are presented for both cluster central regions and total luminosities to be used in comparison with ROSAT and other observational data sets. The quantitative results found disagree significantly with those found by other investigators using semianalytic techniques. We find little dependence of core radius on cluster luminosity and a dependence of temperature on luminosity given by log kT(sub x) = A + B log L(sub x), which is slightly steeper (B = 0.38) than is indicated by observations. Computed temperatures are somewhat higher than observed, as expected, in that COBE-normalized CDM has too much power on the relevant scales. A modest average temperature gradient is found, with temperatures dropping to 90% of central values at 0.4 h(exp -1) Mpc and 70% of central values at 0.9 h(exp -1) Mpc. Examining the ratio of gas to total mass in the clusters normalized to Omega(sub B) h(exp 2) = 0.015, and comparing with observations, we conclude, in agreement with White (1991), that the cluster observations argue for an open universe.

Kang, Hyesung↗
NASA NTRS1995

Collapse and fragmentation of molecular cloud cores. 2: Collapse induced by stellar shock waves

The standard scenario for low-mass star formation involves 'inside-out' collapse of a dense molecular cloud core following loss of magnetic field support through ambipolar diffusion. However, isotopic anomalies in presolar grains and meteoritical inclusions imply that the collapse of the presolar cloud may have been triggered by a stellar shock wave. This paper explores 'outside-in' collapse, that is, protostellar collapse initiated directly by the compression of quiescent dense cloud cores impacted by relatively slow stellar shock waves. A second-order accurate, gravitational hydrodynamics code has been used to study both the spherically symmetrical and three-dimensional evolution of initially centrally condensed, isothermal, self-gravitating, solar-mass cloud cores that are struck by stellar shock waves with velocities up to 25 km/s and postshock temperatures of 10 to 10,000 K. The models show that such mild shock waves do not completely shred and destroy the cloud, and that the dynamical ram pressure can compress the cloud to the verge of self-gravitational collapse. However, compression caused by a high postshock temperature is a considerably more effective means of inducing collapse. Shock-induced collapse produces high initial mass accretion rates (greater than 10(exp -4) solar mass/yr in a solar-mass cloud) that decline rapidly to much lower values, depending on the presence (approximately 10(exp -6) solar mass/yr) or absence (approximately 10(exp -8) to 10(exp -7) solar mass/yr) of an infinite reservoir of mass. Stellar mass accretion rates approximately 10(exp -7) solar mass/yr have been previously inferred from the luminosities of T Tauri stars; balanced mass accretion (stellar rate = envelope rate) at approximately 10(exp -7) solar mass/yr could then be possible if accretion occurs from a finite mass reservoir. Fluid tracers are used to determine what fraction of the stellar shock material is incorporated into the resulting protostellar object and disk; roughly half the impinging material is injected into the collapsing cloud core when there is a high postshock temperature. The models are consistent with a scenario where an AGB star wind triggered the collapse of the presolar cloud while injecting about 0.01 solar mass of matter derived from the AGB star envelope, as has been separately inferred on the basis of nucleosynthesis calculations.

Boss, Alan P.↗
NASA NTRS1996

Research activities in nuclear astrophysics and related areas

NASA/GRO grant NAG 5-2081, at the University of Chicago, has provided support for a broad program of theoretical research in nuclear astrophysics and related areas, with regard to gamma-ray and hard X-ray emission from classical nova explosions. This research emphasized the possible detection of 22Na gamma-ray line emission from nearby novae involving ONeMg white dwarfs, the detailed examination of 26Al production in novae, and the possible detection of the predicted early gamma ray emission from novae that arises from the decay of the short lived, positron emitting isotopes of CNO elements. Studies of nova related problems have consumed an increasing fraction of the Principal Investigator's research efforts over the past decade. Current research addresses problems associated with the standard model for the outbursts of the classical novae: the occurrence of thermonuclear runaways (TNR) in the accreted hydrogen rich envelopes on white dwarfs in close binary systems (see, e.g., the reviews by Truran 1982; and Shara 1989). Research in progress and planned for the next three years has three main objectives: (1) to gain an improved understanding of the early evolution of the light curves of, particularly, the fastest novae; (2) to gain an improved understanding of the relative importance of the various possible mechanisms of envelope hydrogen depletion (e.g. winds, common envelope driven mass loss, and nuclear burning) to the long term evolution of novae in outburst; and (3) to seek to provide a somewhat more definitive statement of the role of classical novae in nucleosynthesis. Our proposed 2-D studies of convection during the early phases of the TNR and our systematic attempt to incorporate an improved treatment of radiation hydrodynamics into the hydrodynamic code utilized in our calculations, are particularly relevant to the first of these objectives. Further 2-D studies of the effects of common envelope evolution are intended to provide more realistic constraints on the mass depletion mechanisms. Finally, detailed calculations of the thermonuclear history of the matter ejected in novae will be carried out for representative nova configurations involving both carbon-oxygen (CO) and oxygen-neon-magnesium (ONeMg) white dwarfs.

Source record↗
NASA NTRS1996

Interplanetary and Interstellar Dust Particles: Reentry Heating and Capture in Aerogel

Studies of interplanetary and interstellar dust can provide significant information on the evolution of the solar system or stars, respectively. However, for reliable analysis it is crucial to know how the particles have been modified during reentry (in the case of interplanetary dust particles, or [DP's) and impact into the capture medium. In the case of stratospheric capture, particles will be heated by atmospheric drag. Subsequent capture of the particles will result in heating, ablation, accretion of the capture medium and possible fragmentation. Modeling of these processes is a useful way of assessing their effects on the interpretation of the compositional data for these particles. Previous work on reentry heating has shown that heat diffusion alone cannot adequately account for temperature gradients observed in IDP's. In fact, for any reasonable thermal parameters, calculations show the particles to be nearly isothermal. Here we extend those studies to include phase transitions. These preliminary results are promising and show a significant temperature gradient for a 100 micron (diameter) particle. The actual capture of the particles in silica aerogel is being modeled using a comprehensive shock hydrodynamic code (called CTH). Various options of this code were investigated to attempt to make the most appropriate choice of methods of impact, equations of state, and processes of energy transfer from capture material to particle. The initial calculations with the code used only 'reasonable' estimates for the physical parameters of silica aerogel. Through the literature searches and personal contacts with the knowledgeable scientists, the best possible mechanical and thermal data have been made available for these computations.

Szydlik, Paul P.↗
NASA NTRS1997

SEDS Tether M/OD Damage Analyses

The Small Expendable Deployer System (SEDS) was designed to deploy an endmass at the end of a 20-km-long tether which acts as an upper stage rocket, and the threats from the meteoroid and orbital debris (M/OD) particle environments on SEDS components are important issues for the safety and success of any SEDS mission. However, the possibility of severing the tether due to M/OD particle impacts is an even more serious concern, since the SEDS tether has a relatively large exposed area to the M/OD environments although its diameter is quite small. The threats from the M/OD environments became a very important issue for the third SEDS mission, since the project office proposed using the shuttle orbiter as a launch platform instead of the second stage of a Delta II expendable rocket, which was used for the first two SEDS missions. A series of hyper-velocity impact tests were performed at the Johnson Space Center and Arnold Engineering Development Center to help determine the critical particle sizes required to sever the tether. The computer hydrodynamic code or hydrocode called CTH, developed by the Sandia National Laboratories, was also used to simulate the damage on the SEDS tether caused by both the orbital debris and test particle impacts. The CTH hydrocode simulation results provided the much needed information to help determine the critical particle sizes required to sever the tether. The M/OD particle sizes required to sever the tether were estimated to be less than 0.1 cm in diameter from these studies, and these size particles are more abundant in low-Earth orbit than larger size particles. Finally, the authors performed the M/OD damage analyses for the three SEDS missions; i.e., SEDS-1, -2, and -3 missions, by using the information obtained from the hypervelocity impact test and hydrocode simulations results.

Hayashida, K. B.↗
NASA NTRS2000

Cosmic Rays and Their Radiative Processes in Numerical Cosmology

A cosmological hydrodynamic code is described, which includes a routine to compute cosmic ray acceleration and transport in a simplified way. The routine was designed to follow explicitly diffusive, acceleration at shocks, and second-order Fermi acceleration and adiabatic loss in smooth flows. Synchrotron cooling of the electron population can also be followed. The updated code is intended to be used to study the properties of nonthermal synchrotron emission and inverse Compton scattering from electron cosmic rays in clusters of galaxies, in addition to the properties of thermal bremsstrahlung emission from hot gas. The results of a test simulation using a grid of 128 (exp 3) cells are presented, where cosmic rays and magnetic field have been treated passively and synchrotron cooling of cosmic ray electrons has not been included.

Ryu, Dongsu↗
NASA NTRS2001

Jet Stability and the Generation of Superluminal and Stationary Components

We present a numerical simulation of the response of an expanding relativistic jet to the ejection of a superluminal component. The simulation has been performed with a relativistic time-dependent hydrodynamical code from which simulated radio maps are computed by integrating the transfer equations for synchrotron radiation. The interaction of the superluminal component with the underlying jet results in the formation of multiple conical shocks behind the main perturbation. These trailing components can be easily distinguished because they appear to be released from the primary superluminal component instead of being ejected from the core. Their oblique nature should also result in distinct polarization properties. Those appearing closer to the core show small apparent motions and a very slow secular decrease in brightness and could be identified as stationary components. Those appearing farther downstream are weaker and can reach superluminal apparent motions. The existence of these trailing components indicates that not all observed components necessarily represent major perturbations at the jet inlet; rather, multiple emission components can be generated by a single disturbance in the jet. While the superluminal component associated with the primary perturbation exhibits a rather stable pattern speed, trailing components have velocities that increase with distance from the core but move at less than the jet speed. The trailing components exhibit motion and structure consistent with the triggering of pinch modes by the superluminal component. The increase in velocity of the trailing components is an indirect consequence of the acceleration of the expanding fluid, which is assumed to be relativistically hot; if observed, such accelerations would therefore favor an electron-positron (as opposed to proton rest mass) dominated jet.

Agudo, Ivan↗
NASA NTRS2002

Simulations of Moon-Forming Impacts

We have incorporated a newly upgraded version of the equation of state ANEOS that allows for the formation of molecular vapor into a smoothed-particle hydrodynamics code. Results of simulations of Moon-forming impacts with this new EOS are presented. Additional information is contained in the original extended abstract.

Canup, R. M.↗
NASA NTRS2002

The Disk Wind Model of the Broad Line Regions in Active Galactic Nuclei and Cataclysmic Variables

This is the final progress report for our Astrophysics Theory Program (NRA 97-OSS12) grant NAG5-7723. We have made considerable progress on incorporating photoionization calculations with a 2.5D hydrodynamical code to model disk winds in AGNs. Following up on our simultaneous broad band monitoring campaign of the type I Seyfert galaxy NGC 5548, we have investigated the constraints imposed on models of accretion in Seyfert galaxies by their optical, UV, and X-ray spectral energy distributions (SEDs). Using results from thermal Comptonization models that relate the physical properties of the hot inner accretion flow to the thermal reprocessing that occurs in the surrounding colder thin disk, we find that we can constrain the central black hole mass, accretion rate and size scale of the hot central flow. We have applied our model to observations of Seyfert galaxies NGC 3516, NGC 7469 and NGC 5548. Our mass and accretion rate estimates for these objects roughly agree with those found using other methods.

Begelman, Mitchell↗
NASA NTRS2002

Tethers as Debris: Hydrocode Simulation of Impacts of Tether Fragments on Planar Aerospace Materials

Tethers promise to find use in a variety of space applications. Despite being narrow objects, their great lengths result in them having large total areas, and so tethers are quite susceptible to being severed by orbital debris. Extensive work has been done designing tethers that resist severs by small debris objects, and hence have longer working lives. It is from this perspective that most recent work has considered the tether - debris question. The potential of intact tethers, or severed tether fragments, as debris to pose a significant collision risk to other spacecraft has been less well studied. Understanding the consequences of such encounters is important in assessing the risks to other spacecraft posed by tethers. In this paper I discuss the damage that two types of tethers may produce on planar aerospace materials, as revealed by hyper- velocity impact simulations using the SPHC hydrodynamic code. Tether types considered include a single nylon line and a complex design including metal wires. Target materials considered include the aluminum plates typically used in debris shielding, and solar panels.

Evans, Steven W.↗
NASA NTRS2003

Tethers as Debris: Hydrocode Simulation of Impacts of Polymer Tether Fragments on Aluminum Plates

Tethers promise to find use in a variety of space applications. Despite being narrow objects, their great lengths result in them having large total areas. Consequently, tethers are very susceptible to being severed by orbital debris. Extensive work has been done designing tethers that resist severs by small debris objects, in order to lengthen their working lives. It is from this perspective that most recent work has considered the tether - debris question. The potential of intact tethers, or severed tether fragments, as debris, to pose a significant collision risk to other spacecraft has been less well studied. Understanding the consequences of such collisions is important in assessing the risks tethers pose to other spacecraft. This paper discusses the damage that polymer tethers may produce on aluminum plates, as revealed by hypervelocity impact simulations using the SPHC hydrodynamic code.

Evans, Steven W.↗
NASA NTRS2004

The Interaction of Hot and Cold Gas in the Disk and Halo of Galaxies

Most of the thermal energy in the Galaxy and perhaps most of the baryons in the Universe are found in hot (log T approximately 5.5 - 7) gas. Hot gas is detected in the local interstellar medium, in supernova remnants (SNR), the Galactic halo, galaxy clusters and the intergalactic medium (IGM). In our own Galaxy, hot gas exists in large superbubbles up to several hundred pc in diameter that locally dominate the interstellar medium (ISM) and determine its thermal and dynamic evolution. While X-ray observations using ROSAT, Chandra and XMM have allowed us to make dramatic progress in mapping out the morphology of the hot gas and in understanding some of its spectral characteristics, there remain fundamental questions that are unanswered. Chief among these questions is the way that hot gas interacts with cooler phase gas and the effects these interactions have on hot gas energetics. The theoretical investigations we proposed in this grant aim to explore these interactions and to develop observational diagnostics that will allow us to gain much improved information on the evolution of hot gas in the disk and halo of galaxies. The first of the series of investigations that we proposed was a thorough exploration of turbulent mixing layers and cloud evaporation. We proposed to employ a multi-dimensional hydrodynamical code that includes non-equilibrium ionization (NEI), radiative cooling and thermal conduction. These models are to be applied to high velocity clouds in our galactic halo that are seen to have O VI by FUSE (Sembach et ai. 2000) and other clouds for which sufficient constraining observations exist.

Slavin, Jonathan↗
NASA NTRS2004

Foam on Tile Impact Modeling for the STS-107 Investigation

Following the breakup of the Space Shuttle Columbia during reentry a NASA/Contractor investigation team was formed to examine the probable damage inflicted on Orbiter Thermal Protection System elements by impact of External Tank insulating foam projectiles. The authors formed a working subgroup within the larger team to apply the Smooth Particle Hydrodynamics code SPHC to the damage estimation problem. Numerical models of the Orbiter's tiles and of the Tank's foam were constructed and used as inputs into the code. Material properties needed to properly model the tiles and foam were obtained from other working subgroups who performed tests on these items for this purpose. Two- and three-dimensional models of the tiles were constructed, including the glass outer layer, the main body of LI-900 insulation, the densified lower layer of LI-900, the Nomex felt mounting layer, and the Aluminum 2024 vehicle skin. A model for the BX-250 foam including porous compression, elastic rebound, and surface erosion was developed. Code results for the tile damage and foam behavior were extensively validated through comparison with Southwest Research Institute foam-on-tile impact experiments carried out in 1999. These tests involved small projectiles striking individual tiles and small tile arrays. Following code and model validation we simulated impacts of larger foam projectiles on the examples of tile systems used on the Orbiter. Results for impacts on the main landing gear door are presented in this paper, including effects of impacts at several angles, and of rapidly rotating projectiles. General results suggest that foam impacts on tiles at about 500 mph could cause appreciable damage if the impact angle is greater than about 20 degrees. Some variations of the foam properties, such as increased brittleness or increased density could increase damage in some cases. Rotation up to 17 rps failed to increase the damage for the two cases considered. This does not rule out other cases in which the rotational energy might lead to an increase in tile damage, but suggests that in most cases rotation will not be an important factor.

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