Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics simulations”

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 289 records · Page 16

Galaxy tracers in N-body simulations

Using the method of smoothed particle hydrodynamics, we have modeled the formation of a compact group of galaxies with sufficient resolution to trace galaxies. Radiative cooling allows the baryons to dissipate their thermal energy and collapse to overdensities characteristic of real galaxies. With their cross section greatly reduced, these galaxy tracers remain distinct during cluster formation while their dark matter halos merge. In addition, the number density, the mass distribution function, and even the morphology of these objects are similar to those of observed galaxies. A viable population of galaxy tracers can be unambiguously defined.

Summers, F. J.↗

Imposing a Lagrangian Particle Framework on an Eulerian Hydrodynamics Infrastructure in Flash

In many astrophysical simulations, both Eulerian and Lagrangian quantities are of interest. For example, in a galaxy cluster merger simulation, the intracluster gas can have Eulerian discretization, while dark matter can be modeled using particles. FLASH, a component-based scientific simulation code, superimposes a Lagrangian framework atop an adaptive mesh refinement Eulerian framework to enable such simulations. The discretization of the field variables is Eulerian, while the Lagrangian entities occur in many different forms including tracer particles, massive particles, charged particles in particle-in-cell mode, and Lagrangian markers to model fluid structure interactions. These widely varying roles for Lagrangian entities are possible because of the highly modular, flexible, and extensible architecture of the Lagrangian framework. In this paper, we describe the Lagrangian framework in FLASH in the context of two very different applications, Type Ia supernovae and galaxy cluster mergers, which use the Lagrangian entities in fundamentally different ways.

Dubey, A.↗

Wind flow over simulated structures for Martian and terrestrial atmospheres

A two-dimensional nonlinear hydrodynamic calculation is presented of a 50-m/sec wind flowing past a simulated Martian crater, and an attempt is made to reproduce a wind flow observed in a California valley. The computer program and boundary conditions used in the calculation are described, and streamlines of the flow over the simulated crater are plotted. The attempt to calculate the wind flow over Owens Valley, California, is made to test the validity of the results for the Martian wind. It is shown that the Martian calculation yields a quasi-steady flow similar to that produced in a wind tunnel simulation of a similar problem. The California wind could not be reproduced with a coarsely resolved calculation, leaving the validity of the Martian calculation and the wind tunnel simulation in doubt.

Jones, E. M.↗

The interplanetary electric field, cleft currents and plasma convection in the polar caps

The relationship between the pattern of plasma convection in the polar cleft and the dynamics of the interplanetary electric field (IEF) is examined theoretically. It is shown that owing to the geometrical properties of the magnetosphere, the East-West component of the IEF will drive field-aligned currents which connect to the ionosphere at points lying on either side of noon, while currents associated with the North-South component of the IEF will connect the two polar caps as sheet currents, also centered at 12 MLT. In order to describe the consequences of the Interplanetary Magnetic Field (IMF) effects upon high-latitude electric fields and convection patterns, a series of numerical simulations was carried out. The simulations were based on a solution to the steady-state equation of current continuity in a height-integrated ionospheric current. The simulations demonstrate that a simple hydrodynamical model can account for the narrow 'throats' of strong dayside antisunward convection observed during periods of southward interplanetary IMF drift, as well as the sunward convection observed during periods of strongly northward IMF drift.

Banks, P. 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.↗

Dynamics of binary-disk interaction. 1: Resonances and disk gap sizes

We investigate the gravitational interaction of a generally eccentric binary star system with circumbinary and circumstellar gaseous disks. The disks are assumed to be coplanar with the binary, geometrically thin, and primarily governed by gas pressure and (turbulent) viscosity but not self-gravity. Both ordinary and eccentric Lindblad resonances are primarily responsible for truncating the disks in binaries with arbitrary eccentricity and nonextreme mass ratio. Starting from a smooth disk configuration, after the gravitational field of the binary truncates the disk on the dynamical timescale, a quasi-equilibrium is achieved, in which the resonant and viscous torques balance each other and any changes in the structure of the disk (e.g., due to global viscous evolution) occur slowly, preserving the average size of the gap. We analytically compute the approximate sizes of disks (or disk gaps) as a function of binary mass ratio and eccentricity in this quasi-equilibrium. Comparing the gap sizes with results of direct simulations using the smoothed particle hydrodynamics (SPH), we obtain a good agreement. As a by-product of the computations, we verify that standard SPH codes can adequately represent the dynamics of disks with moderate viscosity, Reynolds number R approximately 10(exp 3). For typical viscous disk parameters, and with a denoting the binary semimajor axis, the inner edge location of a circumbinary disk varies from 1.8a to 2.6a with binary eccentricity increasing from 0 to 0.25. For eccentricities 0 less than e less than 0.75, the minimum separation between a component star and the circumbinary disk inner edge is greater than a. Our calculations are relevant, among others, to protobinary stars and the recently discovered T Tau pre-main-sequence binaries. We briefly examine the case of a pre-main-sequence spectroscopic binary GW Ori and conclude that circumbinary disk truncation to the size required by one proposed spectroscopic model cannot be due to Linblad resonances, even if the disk is nonviscous.

Artymowicz, Pawel↗

The evolution of X-ray clusters in a cold plus hot dark matter universe

We present the first self-consistently computed results on the evolution of X-ray properties of galaxy clusters in a cold + hot dark matter (CHDM) model. We have performed a hydrodynamic plus N-body simulation for the COBE-compatible CHDM model with standard mass components: Omega(sub hot) = 0.3, Omega (sub cold) = 0.6 and Omega(sub baryon) = 0.1 (h = 0.5). In contrast with the CDM model, which fails to reproduce the observed temperature distribution function dN/dT (Bryan et al. 1994b), the CHDM model fits the observational dN/dT quite well. Our results on X-ray luminosity are less firm but even more intriguing. We find that the resulting X-ray luminosity functions at redshifts z = 0.0, 0.2, 0.4, 0.7 are well fit by observations, where they overlap. The fact that both temperatures and luminosities provide a reasonable fit to the available observational data indicates that, unless we are missing some essential physics, there is neither room nor need for a large fraction of gas in rich clusters: 10% (or less) in baryons is sufficient to explain their X-ray properties. We also see a tight correlation between X-ray luminosity and gas temperature.

Bryan, Greg L.↗

Comet Shoemaker-Levy 9: Impact on Jupiter and plume evolution

The impact of fragments of Comet Shoemaker-Levy 9 on Jupiter and the resulting vapor plume expansion are investigated by conducting three-dimensional numerical simulations using the smoothed particle hydrodynamics (SPH) method. An icy body, representing the cometary fragments, with a velocity of 60 km/sec and a diameter of 2 km can penetrate to 350 km below the 1-bar pressure level in the atmosphere. Most of the initial kinetic energy of the fragment is transferred to the atmosphere between 50 km and 300 km below the 1-bar pressure level. The shock-heated atmospheric gas in the wake is totally dissociated and partially ionized. Scaling our SPH results to other sizes indicates that fragments larger than approximately 100 m in diameter can penetrate to below the visible cloud decks. The energy deposited in the atmosphere is explosively released in the upward expansion of the resulting plume. The plume preferentially expands upward rather than horizontally due to the density gradient of the ambient atmosphere. It rises greater than or equal to 10(exp 2) km in approximately 10(exp 2) sec. Eventually the total atmospheric mass ejected to above 1 bar is greater than or equal to 40 times the initial mass of the impactor. The plume temperature at a radius approximately 10(exp 3) km is greater than 10(exp 3) K for 10(exp 3) sec for a 2-km fragment. We predict that impact-induced plumes will be observable with the remote sensing instruments of the Galileo spacecraft. As the impact site rotates into the view of Earth some 20 min after the impact, the plume expansion will be observable using the Hubble Space Telescope (HST) and from visible and infrared instruments on groundbased telescopes. The rising plume reaches approximately 3000 km altitude in approximately 10 min and will be visible from Earth.

Takata, Toshiko↗

Kinetic Modeling of Plasmaspheric Refilling

We have a new model of ion transport that we have applied to the problem of plasmaspheric flux tube refilling after a geomagnetic disturbance. Effects included are a time-varying ionospheric source, self-consistent Coulomb collisions, field-aligned electric field, and ion cyclotron wave heating. We see refilling rates similar to those of earlier observations and models, except when the electric field is included. In this case, the refilling rates are larger than previously predicted, yet still within the experimental limits. Wave interactions increase the trapped zone population and thus the density, especially during the few few hours of the refilling process. Results are compared with observations from DE/RIMS and Polar/TIDE, as well as hydrodynamic and kinetic particle model simulations.

Liemohn, M. W.↗

Modeling Electric Field Influences on Plasmaspheric Refilling

We have a new model of ion transport that we have applied to the problem of plasmaspheric flux tube refilling after a geomagnetic disturbance. This model solves the Fokker-Planck kinetic equation by applying discrete difference numerical schemes to the various operators. Features of the model include a time-varying ionospheric source, self-consistent Coulomb collisions, field-aligned electric field, hot plasma interactions, and ion cyclotron wave heating. We see refilling rates similar to those of earlier observations and models, except when the electric field is included. In this case, the refilling rates can be quite different that previously predicted. Depending on the populations included and the values of relevant parameters, trap zone densities can increase or decrease. In particular, the inclusion of hot populations near the equatorial region (specifically warm pancake distributions and ring current ions) can dramatically alter the refilling rate. Results are compared with observations as well as previous hydrodynamic and kinetic particle model simulations.

Liemohn, M. W.↗

Comparison of SPHC Hydrocode Results with Penetration Equations and Results of Other Codes

The SPHC hydrodynamic code was used to simulate impacts of spherical aluminum projectiles on a single-wall aluminum plate and on a generic Whipple shield. Simulations were carried out in two and three dimensions. Projectile speeds ranged from 2 kilometers per second to 10 kilometers per second for the single-wall runs, and from 3 kilometers per second to 40 kilometers per second for the Whipple shield runs. Spallation limit results of the single-wall simulations are compared with predictions from five standard penetration equations, and are shown to fall comfortably within the envelope of these analytical relations. Ballistic limit results of the Whipple shield simulations are compared with results from the AUTODYN-2D and PAM-SHOCK-3D codes presented in a paper at the Hypervelocity Impact Symposium 2000 and the Christiansen formulation of 2003.

Evans, Steven W.↗

Meteoroid Risk Assessment of Lunar Habitat Concepts

Preliminary analyses have been performed of meteoroid impacts on lunar surface structures built with elements made from in situ materials. The SPHC hydrodynamic code was used to simulate impacts, and was validated against hypervelocity impact test results. Impact risk will affect necessary wall thicknesses, depth of burial, and placement of lunar base structures.

Evans, Steven W.↗

Low Mach number two-dimensional hydrodynamic turbulence - Energy budgets and density fluctuations in a polytropic fluid

Theory suggests that three distinct types of turbulence can occur in the low Mach number limit of polytropic flow: nearly incompressible flows dominated by vorticity, nearly pure acoustic turbulence dominated by compression, and flows characterized by near statistical equipartition of vorticity and compressions. Distinctions between these kinds of turbulence are investigated here by direct numerical simulation of two-dimensional compressible hydrodynamic turbulence. Dynamical scalings of density fluctuations, examination of the ratio of transverse to longitudinal velocity fluctuations, and spectral decomposition of the fluctuations are employed to distinguish the nature of these low Mach number solutions. A strong dependence on the initial data is observed, as well as a tendency for enhanced effects of compressibility at later times and at higher wave numbers, as suggested by theories of nearly incompressible flows.

Ghosh, S.↗

Experimental investigation of the hydrodynamic forces on the shroud of a centrifugal pump impeller

Fluid-induced forces acting on a rotating impeller are known to cause rotor-dynamic problems in turbomachines. The forces generated by leakage flow along the front shroud surface of a centrifugal turbomachine impeller play an important role among these fluid-induced forces. The present research was aimed to gain a better understanding of these shroud forces. An experimental apparatus was designed and constructed to simulate the impeller shroud leakage flow. Hydrodynamic forces and steady and unsteady pressure distributions on the rotating shroud were measured as functions of eccentricity, width of shroud clearance, face seal clearance and shaft rotating speed. The forces measured from the dynamometer and manometers agreed well. The hydrodynamic force matrices were found skew-symmetric and statically unstable. This is qualitatively similar to the result of previous hydrodynamic volute force measurements. Nondimensionalized normal and tangential forces decrease slightly as Reynolds number increases. As the width of the shroud clearance decreases and/or the eccentricity increases, the hydrodynamic forces increase nonlinearly. There was some evidence found that increased front seal clearance could reduce the radial shroud forces and the relative magnitude of the destabilizing tangential force. Subharmonic pressure fluctuations were also observed which may adversely affect the behavior of the rotor system.

Zhuang, Fei↗

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↗

Monte Carlo simulation of steady state shock structure including cosmic ray mediation and particle escape

Both hydrodynamic calculations (Drury and Volk, 1981, and Axford et al., 1982) and kinetic simulations imply the existence of thermal subshocks in high-Mach-number cosmic-ray-mediated shocks. The injection efficiency of particles from the thermal background into the diffusive shock-acceleration process is determined in part by the sharpness and compression ratio of these subshocks. Results are reported for a Monte Carlo simulation that includes both the back reaction of accelerated particles on the inflowing plasma, producing a smoothing of the shock transition, and the free escape of particles allowing arbitrarily large overall compression ratios in high-Mach-number steady-state shocks. Energy spectra and estimates of the proportion of thermal ions accelerated to high energy are obtained.

Ellison, D. C.↗

Thermocapillary Flows and Their Stability: Effects of Surface Layers and Contamination

The objective is to study quantitatively thermocapillary flows in two-dimensional and axisymmetric geometries in order to learn the impact of such flows on float-zone refining configurations. The work begins with the analytical study of steady flows in slots and zones to obtain predictions of surface curvature, flow and heat transport. It then examines the instability of such flows. Three new instability mechanisms have been identified and the results show that thermocapillary flows of high Prandtl number fluid become convectively unstable while those of low Prandtl number fluid become hydrodynamically unstable. The work includes numerical simulations of steady, high Marangoni number flows in unit aspect ratio boxes and asymptotic solutions in large boxes. Furthermore, the effect of surface contamination on the steady thermocapillary flows was examined. This fundamental study focused on fluid dynamical systems used to understand events in the flow of melts in float-zone geometries both on Earth and in microgravity environments.

Davis, S. H.↗

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