SEARCH · Search NASA

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.

SourceResearch connectionFrom yearThrough year

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

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Gas flows in elliptical galaxies

In preparation for the next generation of x ray telescopes, researchers have begun a program investigating the evolving x ray properties of elliptical galaxies. Their galaxy models consist of a modified King profile for the luminous portion of the galaxy and can include an isothermal dark halo comprising 90 percent of the total mass. The stellar population is assumed to form at a rate which decreases exponentially on a dynamical time scale with a Salpeter initial mass function. Stellar mass loss occurs instantaneously as stars evolve off the main sequence. All stars more massive than 8 solar mass produce type II supernovae, while less massive stars loss mass through a planetary nebulae. The evolving rate of type I supernovae is normalized to a fraction, gamma sub sn I, of Tammann's (1974) value. All of this information is then incorporated into a one-dimensional hydrodynamics code to determine the evolving dynamical state of the interstellar medium.

David, L. P.↗

NASA NTRS1990

Thermal instabilities in protogalactic clouds

The means by which a protogalaxy can fragment to form the first generation of stars and globular clusters remains an important problem in astrophysics. Gravitational instabilities grow on timescales too long to drive fragmentation before the background density grows by many orders of magnitude (see Murray and Lin 1989a, and references therein). Thermal instability provides a much more likely mechanism. After its initial collapse, a protogalactic cloud is expected to be shock heated to its virial temperature approx. 10(exp 6) K. Cooling by H and He+ below 10(exp 6) K has a negative slope, so that the cloud is subject to strong thermal instabilities. Density enhancements may then grow rapidly, fragmenting the protogalaxy as it cools to lower temperatures. The role of dynamical effects upon the growth of perturbations is considered here. The method used is similar to that used in Murray and Lin (1989a; see also the Erratum to appear September 15), which examined the growth of thermal instabilities with a one-dimensional Lagrangian hydrodynamics code, written for spherical symmetry. Perturbed regions therefore take the form of shells. The dynamical variables are integrated explicitly, while the temperature, ionization fraction, and molecular fraction are integrated implicitly, and account is taken for non-equilibrium values of these quantities.

Murray, S. D.↗

NASA NTRS1990

Dynamical experiments on models of colliding disk galaxies

Collisions between galaxies can induce large morphological changes in the participants and, in the case of colliding disk galaxies, bridges and tails are often formed. Observations of such systems indicate a wide variation in color (see Larson and Tinsley, 1978) and that some of the particpants are experiencing enhanced rates of star formation, especially in their central regions (Bushouse 1986, 1987; Kennicutt et al., 1987, Bushouse, Lamb, and Werner, 1988). Here the authors describe progress made in understanding some of the dynamics of interacting galaxies using N-body stellar dynamical computer experiments, with the goal of extending these models to include a hydrodynamical treatment of the gas so that a better understanding of globally enhanced star formation will eventually be forthcoming. It was concluded that close interactions between galaxies can produce large perturbations in both density and velocity fields. The authors measured, via computational experiments that represent a galaxy's stars, average radial velocity flows as large as 100 km/sec and 400 percent density increases. These can occur in rings that move outwards through the disk of a galaxy, in roughly homologous inflows toward the nucleus, and in off center, non-axisymmetric regions. Here the authors illustrate where the gas is likely to flow during the early stages of interaction and in future work they plan to investigate the fate of the gas more realistically by using an N-body/Smoothed Particle Hydrodynamics code to model both the stellar and gaseous components of a disk galaxy during a collision. Specifically, they will determine the locations of enhanced gas density and the strength and location of shock fronts that form during the interaction.

Gerber, Richard A.↗

NASA NTRS1991

Optimization techniques applied to passive measures for in-orbit spacecraft survivability

Spacecraft designers have always been concerned about the effects of meteoroid impacts on mission safety. The engineering solution to this problem has generally been to erect a bumper or shield placed outboard from the spacecraft wall to disrupt/deflect the incoming projectiles. Spacecraft designers have a number of tools at their disposal to aid in the design process. These include hypervelocity impact testing, analytic impact predictors, and hydrodynamic codes. Analytic impact predictors generally provide the best quick-look estimate of design tradeoffs. The most complete way to determine the characteristics of an analytic impact predictor is through optimization of the protective structures design problem formulated with the predictor of interest. Space Station Freedom protective structures design insight is provided through the coupling of design/material requirements, hypervelocity impact phenomenology, meteoroid and space debris environment sensitivities, optimization techniques and operations research strategies, and mission scenarios. Major results are presented.

Mog, Robert A.↗

NASA NTRS1990

Wave pressure in stellar atmospheres due to shock wave trains

Analytic expressions for the wave pressure of propagating shock wave trains in stellar atmospheres or winds are derived. Applications to weak shocks and stronger shocks with sawtooth profiles are discussed in detail. The shocks are treated as discontinuities. The results provide insight in the momentum balance of time-dependent stellar wind flows. The analytic expressions can be used as an independent test of hydrodynamic codes.

Gail, H.-P.↗

NASA NTRS1992

Colliding winds from early-type stars in binary systems

The dynamics of the wind and shock structure formed by the wind collision in early-type binary systems is examined by means of a 2D hydrodynamics code, which self-consistently accounts for radiative cooling, and represents a significant improvement over previous attempts to model these systems. The X-ray luminosity and spectra of the shock-heated region, accounting for wind attenuation and the influence of different abundances on the resultant level and spectra of X-ray emission are calculated. A variety of dynamical instabilities that are found to dominate the intershock region is examined. These instabilities are found to be particularly important when postshock material is able to cool. These instabilities disrupt the postshock flow and add a time variability of order 10 percent to the X-ray luminosity. The X-ray spectrum of these systems is found to vary with the nuclear abundances of winds. These theoretical models are used to study several massive binary systems, in particular V444 Cyg and HD 193793.

Stevens, Ian R.↗

NASA NTRS1992

Numerical models for the collapse and fragmentation of centrally condensed molecular cloud cores

The gravitational collapse and fragmentation of centrally condensed molecular cloud cores are investigated using a new hydrodynamical code. The numerical scheme is second-order accurate and uses explicit finite difference methods to advance the fluid variables on a 3D Cartesian grid. Two initial power-law density profiles, rho varies as r exp -1 and rho varies as r exp -2, are considered, as well as two initial density perturbations in the azimuthal coordinate theta, rho-i - rho(1 + a cos 2theta) where a = 0.1 and 0.5. Fragmentation is found to be possible in these centrally condensed cores if the initial conditions also include differential rotation. Models which collapse with initial uniform rotation do not produce fragments. If molecular cloud cores are indeed centrally condensed, as suggested by observations of star-forming regions and by studies of ambipolar diffusion, then differential rotation may be a mechanism for producing binary protostars during gravitational collapse.

Myhill, Elizabeth A.↗

NASA NTRS1992

Three-dimensional hydrodynamical simulations of colliding stars . III - Collisions and tidal captures of unequal-mass main-sequence stars

A 3D smooth-particle hydrodynamic code with 7000 particles is used to simulate collisions between lower main sequence (MS) stars whose masses differ by a factor of five. Collisions between point-mass intruders and MS stars are simulated to study how the physics of the encounters depends on the finite radius of the smaller MS stars. It is found that the two MS stars become gravitationally bound in a physical collision if their relative velocity at infinity is less than a critical velocity Vd. Vd decreases from 1000 km/s in a head-on collision to 150 km/s in a grazing one. If the less massive star is replaced by a point mass, Vd remains the same for grazing collision and tidal encounters but drops to about 600 km/s in head-on collisions. The collision and coalescence of two MS stars resets the nuclear clock of the coalesced star, so it contracts to the MS. Tidal capture of binary stars in globular clusters can only occur in encounters in which the closest approach of the two stars to their center of mass is less than 2.0 times the sum of their radii.

Benz, W.↗

NASA NTRS1992

Shock breakout in SN 1987A

Detailed models for the UV burst from SN 1987A are calculated using a two-temperature radiation hydrodynamics code. The sequence of events which occur during shock breakout is described, including the formation of a radiative precursor, the disappearance of the shock, and the growth of a thin, dense shell above the photosphere which gives rise to a new viscous, isothermal shock. Continuum spectra are calculated under the assumption that the emergent radiation is a dilute blackbody with a color temperature equal to the radiation temperature at the thermalization depth. This color temperature is 2-3 times higher than the effective temperature. Light-travel-time corrections have been included in the light curves and the spectra. Quantities most relevant to accurate modeling of the fluorescing ring around SN 1987A are highlighted.

Ensman, Lisa↗

NASA NTRS1992

Postcollapse hydrodynamics of SN 1987A - Two-dimensional simulations of the early evolution

The first few seconds of the explosion of SN 1987A are modeled here using a 2D cylindrical geometry smooth particle hydrodynamics code. The success of the explosion is determined to be sensitive to the duration of the infall, the timing of the bounce, and the subsequent neutrino heating. A semianalytical model for the initial structure of the collapsed object is used to present two simulations that differ by the mass that has been allowed to collapse into a neutron star prior to the bounce. In the case of a short initial infall, the explosion fails due to excessive cooling. For a longer initial infall, the cooling is less and the explosion is successful although relatively weak. It is shown that in this case, a successful explosion is brought about by the presence of an entropy gradient which, combined with the gravitational pull of the neutron star, leads to extremely strong instabilities. The critical importance of the global circulation for the success of the explosion is demonstrated.

Herant, Marc↗

NASA NTRS1992

Melting, vaporization, and energy partitioning for impacts on asteroidal and planetary objects

A three-dimensional smoothed particle hydrodynamics code was used to model normal and oblique impacts of silicate projectiles on asteroidal and planetary bodies. The energy of the system, initially in the kinetic energy of the impactor, is partitioned after impact into internal and kinetic energy of the impactor and the target body. These simulations show that, unlike the case of impacts onto a half-space, a significant amount of energy remains in the kinetic energy of the impacting body, as parts of it travel past the main planet and escape the system. This effect is greater for more oblique impacts, and for impacts onto the small planets. Melting and vaporization of both bodies were also examined. The amount of the target body melted was much greater in the case of smaller targets than for an impact of a similar scale on a larger body.

Smither, Catherine L.↗

NASA NTRS1993

Testing the pressure-confined Ly alpha cloud model

The Ly-alpha absorption line forest, seen in quasar spectra, is generally interpreted as being due to cosmologically distributed 'clouds' of primordial gas. Analyses of the observations reveal that the number distribution can be described by power laws: dN/dz is proportional to (1 + z)(sup gamma) and dN/dN(sub HI) is proportional to N(sub HI)(sup -beyda), where N(sub HI) is the HI column density. The typical values for power law indices range between 2 is approximately less than gamma is approximately less than 2.6 and 1.7 is approximately less than gamma is approximately less than 1.9. One model postulates that the Ly-alpha clouds are optically thin entities, photoionized by the background UV flux, J(sub nu) is proportional to (1 + z)(sup j), and confined by an adiabatically evolving intercloud medium (ICM): P(z) is proportional to (1 + z)(sup 5). Analytic studies of this model suggest that the ensuing Ly-alpha line statistics can account for the observations (in particular, the dN/dz and the dN/dN(sub HI) distributions) if the cloud mass spectrum is a power law dN/dN is proportional to M(sup -delta), delta is approximately 1.9, and j is approximately 4. One of the simplifying assumptions incorporated into these studies is the existence of a large mass range for the clouds at all epochs, the validity of which is questionable. The pressure-confined model is investigated using a 1-D spherically symmetric hydrodynamical code to simulate cloud evolution over the epoch 1.8 less than z less than 6. This enables us to relax many of the assumptions incorporated in the analytic studies.

Babul, Arif↗

NASA NTRS1993

Numerical modeling of the interstellar medium in galactic disks

We have been developing detailed hydrodynamic models of the global interstellar medium in the hope of understanding the mass and volume occupied by various phases, as well as their structure and kinematics. In our model, the gas is modeled by one fluid while representative Pop 1 stars are modeled by a second fluid. The two fluids are coupled in that the gas forms into stars at a rate given by a Schmidt law while stellar mass loss returns matter into the gas phase (on a time scale of 100 Myr). Also, the stars heat the gas through stellar winds and the gas cools through optically thin radiation. The time behavior of these two fluids is studied in two spatial dimensions with the Eulerian finite difference numerical hydrodynamic code Zen. The two spatial dimensions are along the plane of a disk (x, total length of 2 kpc) and perpendicular to the disk (z, total height of +/- 15 kpc) and a galactic gravitational field in the z direction, typical of that at the solar circle, is imposed upon the simulation; self-gravity and rotation are absent. For the boundary conditions, outflow is permitted at the top and bottom of the grid (z = +/- 15 kpc) while periodic boundary conditions are imposed upon left and right sides of the grid. As initial conditions, we assumed a gaseous distribution like that seen for the H1 by earlier researchers, although the results are insensitive to the initial conditions. We have run simulations in which the heating due to stars, parameterized as a stellar wind velocity, a, is varied from low (a = 150 km/s), to intermediate (a = 300 km/s), to high (a = 600 km/s). Since the intermediate case is roughly equivalent to the Galactic energy injection rate from supernovae, this summary will concentrate on results from this simulation.

Rosen, A.↗

NASA NTRS1993

The stability of cosmic-ray-dominated shocks - A secondary instability

Previous studies in one spatial dimension showed that the precursors of cosmic-ray-dominated shocks are unstable against traveling acoustic disturbances. Here we report that a secondary, Rayleigh-Taylor type instability can exist in cosmic-ray-dominated media influenced by the acoustic instability. Using the local WKB approximation, the growth rate of the secondary instability is shown to be comparable to that of the one-dimensional acoustic instability itself in the cases we have considered. The nonlinear development has been followed numerically with a two-dimensional PPM hydrodynamics code that also incorporates the two-fluid cosmic-ray energy equation. We show that the secondary instability may cause the precursor and postshock flows to become highly turbulent.

Ryu, Dongsu↗

NASA NTRS1993

Numerical simulations of the convective flame in white dwarfs

A first step toward better understanding of the mechanism driving convective flames in exploding white dwarfs is presented. The propagation of the convective flame is examined using a two-dimensional implicit hydrodynamical code. The large scales of the instability are captured by the grid while the scales that are smaller than the grid resolution are approximated by a mixing-length approximation. It is found that largescale perturbations (of order of the pressure scale height) do grow significantly during the expansion, leading to a very nonspherical burning front. The combustion rate is strongly enhanced (compared to the unperturbed case) during the first second, but later the expansion of the star suppresses the flame speed, leading to only partial incineration of the nuclear fuel. Our results imply that large-scale perturbations by themselves are not enough to explain the mechanism by which convective flames are driven, and a study of the whole spectrum of relevant perturbations is needed. The implications of these preliminary results on future simulations, in the context of current models for Type Ia supernovae, are discussed.

Livne, Eli↗

NASA NTRS1993

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.↗

NASA NTRS1993

A hydrodynamic treatment of the tilted cold dark matter cosmological scenario

A standard hydrodynamic code coupled with a particle-mesh code is used to compute the evolution of a tilted cold dark matter (TCDM) model containing both baryonic matter and dark matter. Six baryonic species are followed, with allowance for both collisional and radiative ionization in every cell. The mean final Zel'dovich-Sunyaev y parameter is estimated to be (5.4 +/- 2.7) x 10 exp -7, below currently attainable observations, with an rms fluctuation of about (6.0 +/- 3.0) x 10 exp -7 on arcmin scales. The rate of galaxy formation peaks at a relatively late epoch (z is about 0.5). In the case of mass function, the smallest objects are stabilized against collapse by thermal energy: the mass-weighted mass spectrum peaks in the vicinity of 10 exp 9.1 solar masses, with a reasonable fit to the Schechter luminosity function if the baryon mass to blue light ratio is about 4. It is shown that a bias factor of 2 required for the model to be consistent with COBE DMR signals is probably a natural outcome in the present multiple component simulations.

Cen, Renyue↗

NASA NTRS1993

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.↗