Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics modeling”

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 595 records · Page 33

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

Diagnostics of electron-heated solar flare models. III - Effects of tapered loop geometry and preheating

A series of hydrodynamic numerical simulations of nonthermal electron-heated solar flare atmospheres and their corresponding soft X-ray Ca XIX emission-line profiles, under the conditions of tapered flare loop geometry and/or a preheated atmosphere, is presented. The degree of tapering is parameterized by the magnetic mirror ratio, while the preheated atmosphere is parameterized by the initial upper chromospheric pressure. In a tapered flare loop, it is found that the upward motion of evaporated material is faster compared with the case where the flare loop is uniform. This is due to the diverging nozzle seen by the upflowing material. In the case where the flare atmosphere is preheated and the flare geometry is uniform, the response of the atmosphere to the electron collisional heating is slow. The upward velocity of the hydrodynamic gas is reduced due not only to the large coronal column depth, but also to the increased inertia of the overlying material. It is concluded that the only possible electron-heated scenario in which the predicted Ca XIX line profiles agree with the BCS observations is when the impulsive flare starts in a preheated dense corona.

Emslie, A. G.↗

Neutral atmospheres of comets - A distributed source model

It is argued that the typical nuclear region of a comet, under a variety of circumstances, consists of a central icy nucleus surrounded by an extended icy halo of grains forming a supplementary source for the observed atmospheric constituents. A complete hydrodynamic description of a neutral atmosphere corresponding to such a distributed source model when the predominant 'parent-molecule' is H2O is given, and numerical results corresponding to two different velocity distributions of the icy halo consistent with smaller and larger grains are computed. The applicability and relevance of these models are discussed.

Ip, W.-H.↗

A model for the recurrent flares in EXO 2030 + 375

It is shown that nonsteady hydrodynamical flows associated with mass and angular momentum capture by a neutron star during a mass ejection phase from a Be star can produce flares with remarkable resemblance to those observed during an outburst from the X-ray transient pulsar EXO 2030 + 375. To reproduce the recurrent time scale of the flares, the velocity of the outflowing matter is estimated to be about 550 km/s. Since the theoretical model requires that a transient disk circulating in one direction is followed by a transient disk circulating in the opposite direction, the time derivative of the pulse period is expected to change sign after each flare event.

Taam, Ronald E.↗

Constraining the Absolute Orientation of eta Carinae's Binary Orbit: A 3-D Dynamical Model for the Broad [Fe III] Emission

We present a three-dimensional (3-D) dynamical model for the broad [Fe III] emission observed in Eta Carinae using the Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS). This model is based on full 3-D Smoothed Particle Hydrodynamics (SPH) simulations of Eta Car's binary colliding winds. Radiative transfer codes are used to generate synthetic spectro-images of [Fe III] emission line structures at various observed orbital phases and STIS slit position angles (PAs). Through a parameter study that varies the orbital inclination i, the PA(theta) that the orbital plane projection of the line-of-sight makes with the apastron side of the semi-major axis, and the PA on the sky of the orbital axis, we are able, for the first time, to tightly constrain the absolute 3-D orientation of the binary orbit. To simultaneously reproduce the blue-shifted emission arcs observed at orbital phase 0.976, STIS slit PA = +38deg, and the temporal variations in emission seen at negative slit PAs, the binary needs to have an i approx. = 130deg to 145deg, Theta approx. = -15deg to +30deg, and an orbital axis projected on the sky at a P A approx. = 302deg to 327deg east of north. This represents a system with an orbital axis that is closely aligned with the inferred polar axis of the Homunculus nebula, in 3-D. The companion star, Eta(sub B), thus orbits clockwise on the sky and is on the observer's side of the system at apastron. This orientation has important implications for theories for the formation of the Homunculus and helps lay the groundwork for orbital modeling to determine the stellar masses.

Madura, T. I.↗

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↗

A model of Comet P/Giacobini-Zinner

A computer model of Comet P/Giacobini-Zinner is presented which contains photoprocesses, gas-phase chemical kinetics, energy balance, multifluid hydrodynamics with a transition to free molecular flow, and solar wind interaction. Recently the physics for electrons in the model has been improved by including electron impact ionization and dissociation and separately accounting for electron energetics. Electron heating and cooling mechanisms include photoprocesses, recombination processes, inelastic and elastic collisions with heavy molecules, and expansion cooling. The model incorporates an internally consistent interaction of the solar wind with the coma gas using the axisymmetric ideal fluid dynamic equations. The nuclear size and composition have been chosen to make the calculations relevant to the September 11, 1985 International Cometary Explorer (ICE) encounter with Comet P/Giacobini-Zinner. Model profiles of the temperature, velocity, and number density of the electrons are in good agreement with measurements along ICE's trajectory. The results indicate the probe passed through a region of the coma at the onset of the plasma tail.

Boice, D. C.↗

Far-infrared emission lines of CO and OH in the Orion-KL molecular shock

Observations of far infrared rotational emission lines which arise in the shocked gas associated with Orion-Kl are presented, including detections of the CO J = 34 yields 33, J = 31 yields 30, J = 26 yields 25, and OH sup 2 PI sub (3/2) J sup P = 7/2(-) yields 5/2(+) emission lines, as well as improved measurements of the CO J = 22 yields 21 and OH sup 2 PI sub (3/2) J = 5/2 yields 3/2 lines. These lines are observed to have velocity widths of Del V approx. 20 to 30 km/sec, somewhat less than either the 2 micro H sub 2 lines or the high velocity plateau component of the millimeter wave CO lines seen in this object. An H sub 2 column density of approx. 3 x 10 to the 21st power, a total mass of approx. 1 solar mass and characteristic temperature and density T approx. 750 K and approx. 2 x 10 to the 6th power per cu cm can be derived from the CO intensities. The density is too low by at least an order of magnitude for the observed infrared H sub 2 and far infrared CO emission to be accounted for by a purely hydrodynamic shock, and support is lent to hydromagnetic shock models. From the present measurements, the relative abundance of CO is estimated to be CO H sub 2 = 1.2 x .0001, corresponding to 20% of the cosmic abundance of C existing in the form of CO. The average relative abundance of OH in the shocked gas is O/H sub 2 or = 5 x 10 to the -7th power. An upper limit to the intensity of the HD J = 1 yields 0 line is used to derive an upper limit of tau or = 3 for the D/H relative abundance in the Orion cloud core.

Watson, D. M.↗

Far-infrared emission lines of CO and OH in the Orion-KL molecular shock

Observations of far infrared rotational emission lines which arise in the shocked gas associted with Orion-Kl are presented, including detections of the CO J = 34 yields 33, J = 31 yields 30, J = 26 yields 25, and OH sup 2 PI sub (3/2) J sup P = 7/2(-) yields 5/2(+) emission lines, as well as improved measurements of the CO J = 22 yields 21 and OH sup 2 PI sub (3/2) J = 5/2 yields 3/2 lines. These lines are observed to have velocity widths of Del V approx. 20 to 30 km/sec, somewhat less than either the 2 micro H sub 2 lines or the high velocity plateau component of the millimeter wave CO lines seen in this object. An H sub 2 column density of aprox. 3 x 10 to the 21st power, a total mass of approx. 1 solar mass and characteristic temperature and density T approx. 750 K and approx. 2 x 10 to the 6th power per cu cm can be derived from the CO intensities. The density is too low by at least an order of magnitude for the observed infrared H sub 2 and far infrared CO emission to be accounted for by a purely hydrodynamic shock, and support is lent to hydromagnetic shock models. From the present measurements, the relative abundance of CO is estimated to be CO H sub 2 = 1.2 x .0001, corresponding to 20 percent of the cosmic abundance of C existing in the form of CO. The average relative abundance of OH in the shocked gas is O/H sub 2 or = 5 x 10 to the -7th power. An upper limit to the intensity of the HD J - 1 yields 0 line is used to derive an upper limit of tau or = 3 for the D/H relative abundance in the Orion cloud core.

Watson, D. M.↗

MHD modelling of the earth's magnetosphere

Attention is given to the shock capturing technique, nonuniform grid system, and multiple time scale problem aspects of a global MHD model of the earth's magnetosphere that uses numerical methods. Because of the fast Alfven waves near the earth, the model is more difficult to solve than problems of hydrodynamic flow past bodies. The Rusanov (1962) scheme, which is a variant of the Lax (1954) method that maintains viscosity at a minimum value except where a large one is required, is used. The magnetosphere shape obtained by these means is noted.

Wu, C. C.↗

Remarks on the MHD problem of generic magnetospheres and magnetotails

An understanding of the plasma physics problem of MHD flow around a magnetized sphere or cylinder in different Reynolds number regimes is shown to be crucial to magnetospheric and magnetotail studies. A significant amount of knowledge can be acquired from experimental and computational hydrodynamic precedents. It is believed that highly symmetric laminar models of the magnetosphere and magnetotail are unlikely in view of the high Reynolds numbers involved.

Montgomery, D.↗

One-dimensional cloud fluid model for propagating star formation

The aim of this project was to study the propagation of star formation (SF) with a self-consistent deterministic model for the interstellar gas. The questions of under what conditions does star formation propagate in this model and what are the mechanisms of the propagation are explored. Here, researchers used the deterministic Oort-type cloud fluid model of Scalo and Struck-Marcell (1984, also see the review of Struck-Marcell, Scalo and Appleton 1987). This cloud fluid approach includes simple models for the effects of cloud collisional coalescence or disruption, collisional energy dissipation, and cloud disruption and acceleration as the result of young star winds, HII regions and supernovae. An extensive one-zone parameter study is presented in Struck-Marcell and Scalo (1987). To answer the questions above, researchers carried out one-dimensional calculations for an annulus within a galactic disk, like the so-called solar neighborhood of the galactic chemical evolution. In the calculations the left-hand boundary is set equal to the right hand boundary. The calculation is obviously idealized; however, it is computationally convenient to study the first order effects of propagating star formation. The annulus was treated as if it were at rest, i.e., in the local rotating frame. This assumption may remove some interesting effects of a supersonic gas flow, but was necessary to maintain a numerical stability in the annulus. The results on the one-dimensional propagation of SF in the Oort cloud fluid model follow: (1) SF is propagated by means of hydrodynamic waves, which can be generated by external forces or by the pressure generated by local bursts. SF is not effectively propagated via diffusion or variation in cloud interaction rates without corresponding density and velocity changes. (2) The propagation and long-range effects of SF depend on how close the gas density is to the critical threshold value, i.e., on the susceptibility of the medium.

Titus, Timothy N.↗

Impact cratering calculations

In the course of carrying out finite difference calculations, it was discovered that for large craters, a previously unrecognized type of crater (diameter) growth occurred which was called lip wave propagation. This type of growth is illustrated for an impact of a 1000 km (2a) silicate bolide at 12 km/sec (U) onto a silicate half-space at earth gravity (1 g). The von Misses crustal strength is 2.4 kbar. The motion at the crater lip associated with this wave type phenomena is up, outward, and then down, similar to the particle motion of a surface wave. It is shown that the crater diameter has grown d/a of approximately 25 to d/a of approximately 4 via lip propagation from Ut/a = 5.56 to 17.0 during the time when rebound occurs. A new code is being used to study partitioning of energy and momentum and cratering efficiency with self gravity for finite-sized objects rather than the previously discussed planetary half-space problems. These are important and fundamental subjects which can be addressed with smoothed particle hydrodynamic (SPH) codes. The SPH method was used to model various problems in astrophysics and planetary physics. The initial work demonstrates that the energy budget for normal and oblique impacts are distinctly different than earlier calculations for silicate projectile impact on a silicate half space. Motivated by the first striking radar images of Venus obtained by Magellan, the effect of the atmosphere on impact cratering was studied. In order the further quantify the processes of meteor break-up and trajectory scattering upon break-up, the reentry physics of meteors striking Venus' atmosphere versus that of the Earth were studied.

Ahrens, Thomas J.↗

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↗

CFD Study of the J2-X LOX Two-Bladed Inducer in the Presence of a Circumferential Groove

Results are presented of a computational fluid dynamics (CFD) study done in support of water flow experiments of the J-2X Oxidizer Turbopump (OTP) 2-bladed inducer with a circumferential groove that were conducted at Marshall Space Flight Center (MSFC). Sub-scale water flow testing results indicate that the circumferential groove greatly reduces synchronous cavitation and subsequent bearing loads at a minimal performance cost, but the energy reappears as high order cavitation (HOC) that spans a relatively large frequency range. Thus, HOC may have implications for the full-scale OTP inducer in terms of reduced structural margin at higher mode frequencies. Simulations using the LOCI-Stream CFD program were conducted in order to explore the root physical cause of the HOC. It was found that the axial recirculation pattern in the circumferential groove causes high-swirl fluid to interact with the nearly-axial incoming fluid just above the inducer blades. The high-shear interface between the fluids is Kelvin-Helmholtz unstable, resulting in trains of low pressure regions or 'pearls' forming near the upstream edge of the groove. When the pressure in these regions becomes low enough and they get cut by the blade leading edge, HOC is thought to occur. Although further work is required, the numerical models indicate that the root cause of HOC is hydrodynamic. That is, the pearls are always present, even when cavitation is not. Comparisons to ongoing water flow experiments will be discussed, as well as predictions for the full-scale OTP inducer.

Coker, Robert↗

Aggregates: The Fundamental Building Blocks of Planetesimals?

The initial accretion of primitive asteroids (meteorite parent bodies) from freely-floating nebula particles remains problematic. Traditional growth-by-sticking models in turbulent nebulae encounter a formidable "meter-size barrier" due to both drift and destruction, or even a mmtocmsize "bouncing" barrier. Even if growth by sticking could somehow breach these barriers (perhaps if the actual sticking or strength is larger than current estimates, which are based on pure ice or pure silicate), turbulent nebulae present further obstacles through the 1-10km size range. On the other hand, nonturbulent nebulae form large asteroids too quickly to explain long spreads in formation times, or the dearth of melted asteroids. Thus, the intensity of nebula turbulence (or "alpha ") is critical to the entire process. Theoretical understanding of nebula turbulence continues to evolve; while recent models of MRI (magneticallydriven) turbulence favor lowornoturbulence environments, purely hydrodynamic turbulence is making a comeback with three recently discovered mechanisms generating turbulence of moderate which do not rely on magnetic fields at all.

Cuzzi, J. N.↗

Intracranial Effects of Microgravity: A Prospective Longitudinal MRI Study

Purpose of Research: To investigate the intracranial effects of microgravity by measuring combined changes in intracranial volumetric parameters, pituitary morphologic structure, and aqueductal cerebrospinal fluid (CSF) hydrodynamics relative to spaceflight and to establish a comprehensive model of recovery after return to Earth.

Larry A Kramer↗