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 343 records · Page 19

A hydrodynamic study of a slow nova outburst

The paper reports use of a Lagrangian implicit hydrodynamics computer code incorporating a full nuclear-reaction network to follow a thermonuclear runaway in the hydrogen-rich envelope of a 1.25 solar-mass white dwarf. In this evolutionary sequence the envelope was assumed to be of normal (solar) composition and the resulting outburst closely resembles that of the slow nova HR Del. In contrast, previous CNO-enhanced models resemble fast nova outbursts. The slow-nova model ejects material by radiation pressure when the high luminosity of the rekindled hydrogen shell source exceeds the local Eddington luminosity of the outer layers. This is in contrast to the fast nova outburst where ejection is caused by the decay of the beta(+)-unstable nuclei. Nevertheless, radiation pressure probably plays a major role in ejecting material from the fast nova remnants. Therefore, the sequence from slow to fast novae can be interpreted as a sequence of white dwarfs with increasing amounts of enhanced CNO nuclei in their hydrogen envelopes, although other parameters such as the white-dwarf mass and accretion rate probably contribute to the observed variation between novae.

Sparks, W. M.↗

The structure and appearance of winds from supercritical accretion disks. I - Numerical models

Equations for the structure and appearance of supercritical accretion disks and the radiation-driven winds which emanate from them are derived and solved by a steady-state hydrodynamic computer code with a relaxation technique used in stellar structure problems. The present model takes into account the mass of the accreting star, the total accretion rate, a generalization of the disk alpha parameter which accounts for heating by processes in addition to viscosity, and the ratio of the total luminosity to the Eddington luminosity. Solutions indicate that for accretion onto a hard-surfaced star, steady, optically thick winds result for even slightly supercritical accretion, and the object will appear as a supergiant star with a high mass loss rate and a nonblackbody spectrum. Winds from black hole accretion disks are expected to depend on the form of the accretion interior to the critical radius, possibly consisting of no ejection at all, a wind similar to that of a hard-surfaced star, or a column of material ejected from a hole in the accretion disk.

Meier, D. L.↗

Fragmentation of rotating protostellar clouds

With a three-dimensional hydrodynamic computer code, the behavior of rotating, isothermal gas clouds as they collapse from Jeans unstable configurations is examined in order to determine whether they are susceptible to fragmentation during the initial dynamic collapse phase of evolution. It is found that a gas cloud will not fragment unless (1) it begins collapsing from a radius much smaller than the Jeans radius (i.e., the cloud initially encloses many Jeans masses) and (2) irregularities in the cloud's initial structure (specifically, density inhomogeneities) enclose more than one Jeans mass of material. Instead of fragmenting, most of the models collapse to a ring configuration. The rings appear to be less susceptible to fragmentation from arbitrary perturbations in their structure than has previously been indicated in other work. Because the models, which include the effects of gas pressure, do not readily fragment during a phase of dynamic collapse, it is suggested that gas clouds in the galactic disk undergo fragmentation only during quasi-equilibrium phases of their evolution.

Tohline, J. E.↗

Numerical calculation of shock-induced initiation of detonations

Results of numerical calculations of the impact of steel cylinders and spheres on the plastic bonded high explosive PBX 9501 are described. The calculations were carried out by a reactive, multicomponent, two dimensional, Eulerian hydrodynamic computer code, 2DE. The 2DE computer code is a finite difference code that uses the donor acceptor cell method to compute mixed cell fluxes. The parameters in the Forest Fire burn model are developed from experiments where the induced shock approximates a plane wave and are applied, in this case, to a situation where the induced shock is a divergent wave with curvature that depends on the size and shape of the projectile. The calculated results are compared with results from experiments involving instrumented mock and live high explosives, with projectiles of varying size, shapes, and velocities.

Cort, G. E.↗

Criteria for fragmentation in a collapsing rotating cloud

The isothermal collapse of a rotating protostar with pressure and self-gravity is calculated by a hydrodynamic computer code that treats three space dimensions. A number of initial conditions are tested to determine the conditions under which a cloud is unstable to fragmentation. It is shown that fragmentation does not proceed to a significant extent during the first free-fall time, but that in most cases fragmentation has occurred by 1.5 to 2 times the initial free-fall time. The extent to which the density in the fragments is enhanced over that in their surroundings on this time scale depends on the initial ratio of thermal to gravitational energy, but is only weakly dependent on the type and amplitude of the perturbation and on the initial rotational energy.

Bodenheimer, P.↗

Collapse and equilibrium of rotating, adiabatic clouds

A numerical hydrodynamics computer code analysis of the collapse and establishment of equilibrium of adiabatic gas clouds restricted to axial symmetry found that the clouds are originally uniform in density and rotation. The method can compare the dynamic collapse and approach to equilibrium with the data on incompressible uniformly rotating equilibrium clouds and on equilibrium structures of differentially rotating polytropes. It is concluded that the stellar formation theory indicates that the low alpha fragments produced at the termination of the dynamic isothermal collapse phase of interstellar clouds may undergo significant dynamic collapse in an adiabatic regime leading to transitory ring formation and additional fragmentation on a smaller scale.

Boss, A. P.↗

Numerical simulation of the magnetospheric gate model for X-ray bursters

A Lagrangian, fully implicit, one dimensional hydrodynamic computer code was used to investigate the evolution of a gas cloud impacting the surface of a 20 km, 1 Msub solar neutron star. This gas is initially at rest with respect to the surface of the neutron star, extends to 185 km above the surface, and is optically thick. The infall results in a burst which lasts about 0.1 seconds and reached a peak luminosity and effective temperature of 240,000 Lsub solar and 9 million; respectively. The burst was followed by a phase of oscillations with a period 0.2 seconds.

Starrfield, S. G.↗

Three-dimensional computer model for the atmospheric general circulation experiment

An efficient, flexible, three-dimensional, hydrodynamic, computer code has been developed for a spherical cap geometry. The code will be used to simulate NASA's Atmospheric General Circulation Experiment (AGCE). The AGCE is a spherical, baroclinic experiment which will model the large-scale dynamics of our atmosphere; it has been proposed to NASA for future Spacelab flights. In the AGCE a radial dielectric body force will simulate gravity, with hot fluid tending to move outwards. In order that this force be dominant, the AGCE must be operated in a low gravity environment such as Spacelab. The full potential of the AGCE will only be realized by working in conjunction with an accurate computer model. Proposed experimental parameter settings will be checked first using model runs. Then actual experimental results will be compared with the model predictions. This interaction between experiment and theory will be very valuable in determining the nature of the AGCE flows and hence their relationship to analytical theories and actual atmospheric dynamics.

Roberts, G. O.↗

The origin of the moon and the single impact hypothesis. III

Calculations of the single-impact hypothesis for the origin of the moon were performed using the smoothed particle hydrodynamics (SPH) code developed by Benz et al. (1986). Results are presented from calculations of a relatively low-level collision with an impactor mass in the range 6-8 x 10 to the 26th g. Several runs of the calculations are conducted for this mass range with variations in the SPH code, the equation of state, and the initial planetary models. The effects of these variations are compared. It is found that the orbiting mass is injected by gravitational torques.

Benz, W.↗

Tidal disruption of inviscid planetesimals

In view of previous efforts' demonstration that strongly dissipative planetesimals are immune to tidal disruption, an examination is presently conducted of the complementary case of inviscid planetesimals arising from collisions that are sufficiently energetic to entirely melt the resulting planetesimal and debris. The tidal disruption is numerically simulated by means of the smoothed particle hydrodynamics (SPH) code of Cameron and Benz (1991), concentrating on the tidal disruption of 0.01 earth-mass planetesimals passing by the earth with variations in the impact parameter at perigee and velocity at infinity. The SPH models show that tidal forces during a close encounter can efficiently convert orbital angular momentum into spin angular momentum, thereby initiating equatorial mass-shedding to inviscid planetesimals that have been spun up beyond the limit of rotational stability.

Boss, A. P.↗

Tidal disruption of inviscid protoplanets

Roche showed that equilibrium is impossible for a small fluid body synchronously orbiting a primary within a critical radius now termed the Roche limit. Tidal disruption of orbitally unbound bodies is a potentially important process for planetary formation through collisional accumulation, because the area of the Roche limit is considerably larger then the physical cross section of a protoplanet. Several previous studies were made of dynamical tidal disruption and different models of disruption were proposed. Because of the limitation of these analytical models, we have used a smoothed particle hydrodynamics (SPH) code to model the tidal disruption process. The code is basically the same as the one used to model giant impacts; we simply choose impact parameters large enough to avoid collisions. The primary and secondary both have iron cores and silicate mantles, and are initially isothermal at a molten temperature. The conclusions based on the analytical and numerical models are summarized.

Boss, Alan P.↗

A new scheme for multidimensional line transfer. II - ETLA method in one dimension with application to iron K-alpha lines

The 1D version of the radiation hydrodynamics computer code ALTAIR is discribed. It is an implementation of the equivalent two-level atom method, and also incorporates a variety of other iterative techniques to achieve a fast, accurate solution of the problem of coupled atomic kinetics and radiative transfer. All the techniques are described in sufficient detail to enable their use by others. The application of the program is illustrated by a set of calculations of the iron X-ray emission spectrum from intercloud hot gas that may exist in the broad-line region of active galactic nuclei. These calculations demonstrate that such spectra may be a rich load to mine for information about the hot gas in AGNs, and also that the concept of 'the iron K-alpha line' inadequately depicts the spectrum.

Castor, John I.↗

A comparative study between experimental results and numerical predictions of multi-wall structural response to hypervelocity impact

Over the last three decades, multiwall structures have been analyzed extensively, primarily through experiment, as a means of increasing the protection afforded to spacecraft structure. However, as structural configurations become more varied, the number of tests required to characterize their response increases dramatically. As an alternative, numerical modeling of high-speed impact phenomena is often being used to predict the response of a variety of structural systems under impact loading conditions. This paper presents the results of a preliminary numerical/experimental investigation of the hypervelocity impact response of multiwall structures. The results of experimental high-speed impact tests are compared against the predictions of the HULL hydrodynamic computer code. It is shown that the hypervelocity impact response characteristics of a specific system cannot be accurately predicted from a limited number of HULL code impact simulations. However, if a wide range of impact loadings conditions are considered, then the ballistic limit curve of the system based on the entire series of numerical simulations can be used as a relatively accurate indication of actual system response.

Schonberg, William P.↗

The LDCE Particle Impact Experiment as flown on STS-46

Many materials and techniques have been developed by the authors to sample the flux of particles in Low Earth Orbit (LEO). Though regular in-site sampling of the flux in LEO the materials and techniques have produced data which compliment the data now being amassed by the Long Duration Exposure Facility (LDEF) research activities. Orbital debris models have not been able to describe the flux of particles with d sub p less than or = 0.05 cm, because of the lack of data. Even though LDEF will provide a much needed baseline flux measurement, the continuous monitoring of micron and sub-micron size particles must be carried out. A flight experiment was conducted on the Space Shuttle as part of the LDCE payload to develop an understanding of the Spatial Density (concentration) as a function of size (mass) for particle sizes 1 x 10(exp 6) cm and larger. In addition to the enumeration of particle impacts, it is the intent of the experiment that hypervelocity particles be captured and returned intact. Measurements will be performed post flight to determine the flux density, diameters, and subsequent effects on various optical, thermal control and structural materials. In addition to these principal measurements, the Particle Impact Experiment (PIE) also provides a structure and sample holders for the exposure of passive material samples to the space environment, e.g., thermal cycling, and atomic oxygen, etc. The experiment will measure the optical property changes of mirrors and will provide the fluence of the ambient atomic oxygen environment to other payload experimenters. In order to augment the amount of material returned in a form which can be analyzed, the survivability of the experiment as well as the captured particles will be assessed. Using Sandia National Laboratory's hydrodynamic computer code CTH, hypervelocity impacts on the materials which comprise the experiments have been investigated and the progress of these studies are reported.

Maag, Carl R.↗

TSS tether cable meteoroid/orbital debris damage analysis

This report summarizes the damage analyses performed on the tether cable used for the tethered satellite system (TSS), for the damage that could be caused by meteoroid or orbital debris impacts. The TSS consists of a tethered satellite deployer and a tethered satellite. The analytical studies were performed at Marshall Space Flight Center (MSFC) with the results from the following tests: (1) hypervelocity impact tests to determine the 'critical' meteoroid particle diameter, i.e., the maximum size of a meteoroid particle which can impact the tether cable without causing 'failure'; (2) electrical resistance tests on the damaged and undamaged tether cable to determine if degradation of current flow occurred through the damaged tether cables; and (3) tensile load tests to verify the load carrying capability of the damaged tether cables. Finally, the HULL hydrodynamic computer code was used to simulate the hypervelocity impact of the tether cable by particles at velocities higher than can be tested, to determine the extent of the expected tether damage.

Hayashida, K. B.↗

The formation process of the He I lambda 10830 line in cool giant stars

The Final Report on the formation process of the He I lambda 10830 line in cool giant stars is presented. The research involves observing a sample of cool giant stars with ROSAT. These stars were selected from the list of bright stars which display He I lambda 10830 in absorption or emission and lie on the cool side of the coronal dividing line. With measured x ray fluxes or upper limits measured by the Position Sensitive Proportional Counter (PSPC), the role x rays play in the formation of this important line was investigated using the non-LTE radiative transfer code PANDORA. Hydrodynamic calculations were performed to investigate the contributions of acoustic wave heating in the formation of this line as well.

Luttermoser, Donald G.↗

A study of angular momentum loss in binaries using the free Lagrange method

The evolution of a binary star system depends greatly on the angular momentum losses in the system brought about by gravitational radiation and mass outflow (e.g., evaporating winds and magnetic braking) from the secondary component of the binary. Using a three-dimensional hydrodynamic fluid code based on the free Lagrange method, we study the loss of specific angular momentum from a binary system due to an evaporative wind from the companion of a millisecond pulsar. We consider binaries of different mass ratios and winds of different initial velocities and in particular attempt to model the system PSR 1957+20. We are in the process of incorporating the effect of the radiation force from the pulsar and the magnetic field of the companion on the mass outflow. The latter effect would also enable us to study magnetic braking in cataclysmic variables and low-mass X-ray binaries.

Rajasekhar, A. M.↗

Shoemaker-Levy 9 and the tidal disruption of comets

The break-up of Periodic Comet Shoemaker-Levy 9 into multiple pieces following its grazing encounter with Jupiter in July 1992 can be used to study tidally-induced fracture in comets. This spectacular event allows us not only to set limits on the size, strength and density of Shoemaker-Levy 9 itself, but provides invaluable guidance to numerical modeling of such encounters. In an extensive treatment of tidal breakup which assumed self-gravitating, homogeneous, perfectly elastic bodies, Dobrovolskis derived simple analytical expressions for the tidally-induced surface and central stresses. Both can be cast in such a way that Poisson's ratio is the only material dependent constant entering these expressions. Whether both surface and central failure must be initiated as a criterion for breakup, or either one of them is sufficient, remains a subject of disagreement. To resolve this debate, we model the details of cometary breakup using a three-dimensional Smooth Particle Hydrodynamics (SPH) code modified to simulate fracture in small solid objects. At lower stresses associated with brittle failure, we use a rate-dependent strength based on the nucleation of incipient flaws whose number density is given by a Weibull distribution.

Benz, W.↗