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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
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.↗
NASA NTRS1993

Initiating solar system formation through stellar shock waves

Isotopic anomalies in presolar grains and other meteoritical components require nucleosynthesis in stellar interiors, condensation into dust grains in stellar envelopes, transport of the grains through the interstellar medium by stellar outflows, and finally injection of the grains into the presolar nebula. The proximity of the presolar cloud to these energetic stellar events suggests that a shock wave from a stellar outflow might have initiated the collapse of an otherwise stable presolar cloud. We have begun to study the interactions of stellar shock waves with thermally supported, dense molecular cloud cores, using a three spatial dimension (3D) radiative hydrodynamics code. Supernova shock waves have been shown by others to destroy quiescent clouds, so we are trying to determine if the much smaller shock speeds found in, e.g., asymptotic giant branch (AGB) star winds, are strong enough to initiate collapse in an otherwise stable, rotating, solar-mass cloud core, without leading to destruction of the cloud.

Boss, A. P.↗
NASA NTRS1993

A mixed method Poisson solver for three-dimensional self-gravitating astrophysical fluid dynamical systems

A key ingredient in the simulation of self-gravitating astrophysical fluid dynamical systems is the gravitational potential and its gradient. This paper focuses on the development of a mixed method multigrid solver of the Poisson equation formulated so that both the potential and the Cartesian components of its gradient are self-consistently and accurately generated. The method achieves this goal by formulating the problem as a system of four equations for the gravitational potential and the three Cartesian components of the gradient and solves them using a distributed relaxation technique combined with conventional full multigrid V-cycles. The method is described, some tests are presented, and the accuracy of the method is assessed. We also describe how the method has been incorporated into our three-dimensional hydrodynamics code and give an example of an application to the collision of two stars. We end with some remarks about the future developments of the method and some of the applications in which it will be used in astrophysics.

Duncan, Comer↗
NASA NTRS1994

Axial focusing of impact energy in the Earth's interior: Proof-of-principle tests of a new hypothesis

A causal link between major impact events and global processes would probably require a significant change in the thermal state of the Earth's interior, presumably brought about by coupling of impact energy. One possible mechanism for such energy coupling from the surface to the deep interior would be through focusing due to axial symmetry. Antipodal focusing of surface and body waves from earthquakes is a well-known phenomenon which has previously been exploited by seismologists in studies of the Earth's deep interior. Antipodal focusing from impacts on the Moon, Mercury, and icy satellites has also been invoked by planetary scientists to explain unusual surface features opposite some of the large impact structures on these bodies. For example, 'disrupted' terrains have been observed antipodal to the Caloris impact basis on Mercury and Imbrium Basin on the Moon. Very recently there have been speculations that antipodal focusing of impact energy within the mantle may lead to flood basalt and hotspot activity, but there has not yet been an attempt at a rigorous model. A new hypothesis was proposed and preliminary proof-of-principle tests for the coupling of energy from major impacts to the mantle by axial focusing of seismic waves was performed. Because of the axial symmetry of the explosive source, the phases and amplitudes are dependent only on ray parameter (or takeoff angle) and are independent of azimuthal angle. For a symmetric and homogeneous Earth, all the seismic energy radiated by the impact at a given takeoff angle will be refocused (minus attenuation) on the axis of symmetry, regardless of the number of reflections and refractions it has experienced. Mantle material near the axis of symmetry will experience more strain cycles with much greater amplitude than elsewhere and will therefore experience more irreversible heating. The situation is very different than for a giant earthquake, which in addition to having less energy, has an asymmetric focal mechanism and a larger area. Two independent proof-of-principle approaches were used. The first makes use of seismic simulations, which are being performed with a realistic Earth model to determine the degree of focusing along the axis and to estimate the volume of material, if any, that experiences significant irreversible heating. The second involves two-dimensional hydrodynamic code simulations to determine the stress history, internal energy, and temperature rise as a function of radius along the axis.

Boslough, M. B.↗
NASA NTRS1995

Dimensional scaling for impact cratering and perforation

POD Associates have revisited the issue of generic scaling laws able to adequately predict (within better than 20 percent) cratering in semi-infinite targets and perforations through finite thickness targets. The approach used was to apply physical logic for hydrodynamics in a consistent manner able to account for chunky-body impacts such that the only variables needed are those directly related to known material properties for both the impactor and target. The analyses were compared and verified versus CTH hydrodynamic code calculations and existing data. Comparisons with previous scaling laws were also performed to identify which (if any) were good for generic purposes. This paper is a short synopsis of the full report available through the NASA Langley Research Center, LDEF Science Office.

Watts, Alan J.↗
NASA NTRS1994

X-ray clusters from a high-resolution hydrodynamic PPM simulation of the cold dark matter universe

A new three-dimensional hydrodynamic code based on the piecewise parabolic method (PPM) is utilized to compute the distribution of hot gas in the standard Cosmic Background Explorer (COBE)-normalized cold dark matter (CDM) universe. Utilizing periodic boundary conditions, 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, we find the X-ray-emitting clusters, compute the luminosity function at several wavelengths, the temperature distribution, and estimated sizes, as well as the evolution of these quantities with redshift. The results, which are compared with those obtained in the preceding paper (Kang et al. 1994a), may be used in conjuction with ROSAT and other observational data sets. Overall, the results of the two computations are qualitatively very similar with regard to the trends of cluster properties, i.e., how the number density, radius, and temeprature depend on luminosity and redshift. The total luminosity from clusters is approximately a factor of 2 higher using the PPM code (as compared to the 'total variation diminishing' (TVD) code used in the previous paper) with the number of bright clusters higher by a similar factor. The primary conclusions of the prior paper, with regard to the power spectrum of the primeval density perturbations, are strengthened: the standard CDM model, normalized to the COBE microwave detection, predicts too many bright X-ray emitting clusters, by a factor probably in excess of 5. The comparison between observations and theoretical predictions for the evolution of cluster properties, luminosity functions, and size and temperature distributions should provide an important discriminator among competing scenarios for the development of structure in the universe.

Bryan, Greg L.↗
NASA NTRS1994

A hydrodynamic approach to cosmology: The mixed dark matter cosmological scenario

We compute the evolution of spatially flat, mixed cold and hot dark matter models containing both baryonic matter and two kinds of dark matter. Hydrodynamics is treated with a highly developed Eulerian hydrodynamic code (see Cen 1992). A standard particle-mesh (PM) code is also used in parallel to calculate the motion of the dark matter components. We adopt the following parameters: h equivalent to (sub 0)/100 km/s Mpc(exp -1) = 0.5, OMEGA(sub C) = 0.3, and OMEGA(sub B) = 0.06, with amplitude of the perturbation spectrum fixed by the Cosmic Background Explorer Satellite (COBE) Dark Matter Radiation (DMR) measurements (Smoot et al. 1992) being sigma (sub 8) = 0.67. Four different boxes are simulated with box sizes of L = (64, 16, 4, 1) h(exp -1) Mpc, respectively, the two small boxes providing good resolution but little valid information due to absence of large-scale power. We use 128(exp 3) approximate 10(exp 6.3) baryonic cells, 128(exp .3) cold dark matter particles, and 2 x 128(exp 3) hot dark matter particles. In addition to the dark matter we follow separately six baryonic species (H, H(+), He, He(+), He(++), e(-)) with allowance for both (nonequilibrium) collisional and radiative ionization in every cell. The background radiation field is also followed in detail with allowance made for both continuum and line processes, to allow nonequilibrium heating and cooling processes to be followed in detail. The mean final Zeldovich-Sunyaev y parameter is estimated to be y Bar = (5.4 + or - 2.7) x 10(exp -7) below currently attainable observations, with a rms fluctuation of approximately delta bar y = (0.6 + or - 3.0) x 10(exp -7) on arcminute scales. The rate of galaxy formation peaks at an even later epoch (z approximate 0.3) than in the standard (OMEGA = 1, sigma sub 8 = 0.67) cold dark matter (CDM) model (z approximate 0.5) and, at a redshift of z = 4, is nearly a factor of 100 lower than for the CDM model with the same value of sigma sub 8. With regard to mass function, the smallest objects are stabilized against collapse by thermal energy: the mass-weighted mass spectrum has a broad peak in the vicinity of M(sub B) = 10(exp 9.5) solar mass with a reasonable fit to the Schechter luminosity function if the ratio of baryon mass to blue light is approximately 4. In addition, one very large PM simulation was made in a box with size (320 h(exp - 1) Mpc) containing 3 x 200(exp 3) = 10(exp 7.4) particles. Utilizing this simulation we find that the model yields a cluster mass function which is about a factor of 4 higher than observed, but a cluster-cluster correlation length marginally lower than observed, but that both are closer to observations than in the (COBE) normalized CDM model. The one-dimensional pairwise velocity dispersion is 605 + or - 8 km/s at 1/h separation, lower than that of the DCM model normalized to COBE, but still significant higher than observations (Davis & Peebles 1983). A plausible velocity bias b(sub v) = 0.8 + or - 0.1 on this scale will reduce but not remove the discrepancy. The velocity auto-correlat ion function has a coherence length of 40/h Mpc, which is somewhat lower than the observed counterpart. In all these respects the model would be improved by decreasing the cold fraction of the dark OMEGA(sub CDM)/ (OMEGA(sub CDM) + OMEGA(sub HDB). But formation of galaxies and clusters of galaxies is much later in this model than in COBE-normalized CDM, perhaps too late. To improve on these constraints a larger ratio of OMEGA(sub CDM)/ (OMEGA(sub CDM) + OMEGA(sub HDM)) is required than the value of 0.67 adopted here. It does not seem possible to find a value for this ratio which would satisfy all tests. Overall, the model is similar both on large and intermediate scales to the standard CDM model normalized to the same value of sigma(sub B), but the problem with regard to late formation of galaxies is more severe in this model than in that CDM model. Adding hot dark matter, significantly improves the ability of the COBE-normalized CDM scenario to fit existing observations, but the model is in fact not as good as the CDM model with the same sigma(sub 8) and is still probably unsatisfactory with regard to several critical tests.

Cen, Renyue↗
NASA NTRS1994

Acoustic waves in M dwarfs: Maintaining a corona

We use a time-dependent hydrodynamics code to follow the propagation of acoustic waves into the corona of an M dwarf star. An important qualitative difference between M dwarfs and stars such as the Sun is that the acoustic spectrum in M dwarfs is expected to peak at periods close to the acoustic cutoff P(sub A): this allows more effective penetration of waves into the corona. In our code, radiative losses in the photosphere, chromosphere, and corona are computed using Rosseland mean opacities, Mg II kappa and Ly alpha emission, and optically thin emissivities respectively. We find that acoustic heating can maintain a corona with a temperature of order 0.7-1 x 10(exp 6) K and a surface X-ray flux as large as 10(exp 5)ergs/sq cm/s. In a recent survey of X-rays from M dwarfs, some (20%-30%) of the stars lie at or below this limiting X-ray flux: we suggest that such stars may be candidates for acoustically maintained coronae.

Mullan, D. J.↗