Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamic simulation”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

A Recent Impact Origin of Saturn’s Rings and Mid-Sized Moons

We simulate the collision of precursor icy moons analogous to Dione and Rhea as a possible origin forSaturn’s remarkably young rings. Such an event could have been triggered a few hundred million years ago by resonant instabilities in a previous satellite system. Using high-resolution SPH simulations, we find that this kind of impact can produce a wide distribution of massive objects and scatter material throughout the system. This includes the direct placement of pure-ice ejecta onto orbits that enterSaturn’s Roche limit, which could form or rejuvenate rings. In addition, fragments and debris of rock and ice totalling more than the mass of Enceladus can be placed onto highly eccentric orbits that would intersect with any precursor moons orbiting in the vicinity of Mimas, Enceladus, or Tethys. This could prompt further disruption and facilitate a collisional cascade to distribute more debris for potential ring formation, the re-formation of the present-day moons, and to evolve into an eventual cratering population of planeto-centric impactors.

Saturnian satellites↗
NASA NTRS2021

Crucial Factors for Lyα Transmission in the Reionizing Intergalactic Medium: Infall Motion, H II Bubble Size, and Self-shielded Systems

Using the CoDa II simulation, we study the Lyα transmissivity of the intergalactic medium (IGM) during reionization. At z>6, a typical galaxy without an active galactic nucleus fails to form a proximity zone around itself due to the overdensity of the surrounding IGM. The gravitational infall motion in the IGM makes the resonance absorption extend to the red side of Lyα, suppressing the transmission up to roughly the circular velocity of the galaxy. In some sight lines, an optically thin blob generated by a supernova in a neighboring galaxy results in a peak feature, which can be mistaken for a blue peak. Redward of the resonance absorption, the damping-wing opacity correlates with the global IGM neutral fraction and the UV magnitude of the source galaxy. Brighter galaxies tend to suffer lower opacity because they tend to reside in larger H II regions, and the surrounding IGM transmits redder photons, which are less susceptible to attenuation, owing to stronger infall velocity. The H II regions are highly nonspherical, causing both sight-line-to-sight-line and galaxy-to-galaxy variation in opacity. Also, self-shielded systems within H II regions strongly attenuate the emission for certain sight lines. All these factors add to the transmissivity variation, requiring a large sample size to constrain the average transmission. The variation is largest for fainter galaxies at higher redshift. The 68% range of the transmissivity is similar to or greater than the median for galaxies with M(sub UV) ≽ -21 at z ≽ 7, implying that more than a hundred galaxies would be needed to measure the transmission to 10% accuracy.

Reionization↗
NASA NTRS1976

A three-dimensional, time-dependent model of Mobile Bay

A three-dimensional, time-variant mathematical model for momentum and mass transport in estuaries was developed and its solution implemented on a digital computer. The mathematical model is based on state and conservation equations applied to turbulent flow of a two-component, incompressible fluid having a free surface. Thus, bouyancy effects caused by density differences between the fresh and salt water, inertia from thare river and tidal currents, and differences in hydrostatic head are taken into account. The conservation equations, which are partial differential equations, are solved numerically by an explicit, one-step finite difference scheme and the solutions displayed numerically and graphically. To test the validity of the model, a specific estuary for which scaled model and experimental field data are available, Mobile Bay, was simulated. Comparisons of velocity, salinity and water level data show that the model is valid and a viable means of simulating the hydrodynamics and mass transport in non-idealized estuaries.

Pitts, F. H.↗
NASA NTRS1979

The interaction between giant gaseous protoplanets and the primitive solar nebula

The manner in which a giant gaseous protoplanet becomes embedded in the primitive solar nebula determines surface boundary conditions which must be used in studying the evolution of such objects. On the one hand, if the system resembles a contact binary system, then the envelope of the protoplanet should approach the entropy of the surrounding nebula. On the other hand angular momentum transfer by resonance and tidal effects between the nebula and the protoplanet may cause the nebula to exhibit a zone of avoidance near the protoplanet, thus inhibiting exchange of material. This problem has been studied with a computer program developed by D. N. C. Lin which simulates disk hydrodynamics by particle motions with dissipation. These studies suggest that for expected values of the protoplanet/protosun mass ratios, significant inhibition of mass exchange is likely, so that it is a reasonable next step to undertake protoplanet evolution studies with the imposition of minimum protoplanet surface temperatures.

Cameron, A. G. W.↗
NASA NTRS1980

Modeling of vortex generated sound in solid propellant rocket motors

There is considerable evidence based on both full scale firings and cold flow simulations that hydrodynamically unstable shear flows in solid propellant rocket motors can lead to acoustic pressure fluctuations of significant amplitude. Although a comprehensive theoretical understanding of this problem does not yet exist, procedures were explored for generating useful analytical models describing the vortex shedding phenomenon and the mechanisms of coupling to the acoustic field in a rocket combustion chamber. Since combustion stability prediction procedures cannot be successful without incorporation of all acoustic gains and losses, it is clear that a vortex driving model comparable in quality to the analytical models currently employed to represent linear combustion instability must be formulated.

Flandro, G. A.↗
NASA NTRS1989

Direct numerical simulations of turbulent convection with a variable gravity and Keplerian rotation

Thermal convection was proposed as a possible mechanism for generation and maintenance of turbulence in the inner accretion disk regime of the primordial solar nebula. It is of fundamental interest to design experiments with the basic physical features of the solar nebula conditions cannot be produced in the laboratory, numerical simulations of hydrodynamic flows, which have been very successful in describing aerodynamic flows, can be suitable modified to provide experimental data for solar nebula modelling. The goals are to modify an extant, proven hydrodynamics code with the most important features of the solar nebula and other thin accretion disks: bouyancy terms to generate convection, internal heating representing the release of gravitational potential energy, a variable gravity linearly proportional the the distance from the vertical midplane due to centrifugal balance, rapid rotation with axis aligned with gravity, and Keplerian rotational shear; to determine the effect that these features have on the turbulent convection by introducing them individually and to determine the cumulative nature of the turbulent convection for accretion disk conditions; and to model the convection and the turbulence. In this manner, prior solar nebula models can be tested and their deficiencies rectified.

Cabot, William H.↗
NASA NTRS1993

Theoretical studies of a molecular beam generator

The following is a proposed baseline model that is being developed for the simulation of hydrodynamic generator, which can be converted at a later date to magnetohydrodynamic MHD thruster by adding the necessary electric and magnetic fields. The following development will include the electric and magnetic terms, however, the initial computer program will not include these terms. The analysis that follows is a one species, single temperature model constructed over a domain that is defined by an enclosed region.

Heinbockel, John H.↗
NASA NTRS1993

The Kelvin-Helmholtz instability in photospheric flows - Effects of coronal heating and structure

A series of hydrodynamic numerical simulations has been used to investigate the nonlinear evolution of driven, subsonic velocity shears under a range of typical photospheric conditions. These calculations show that typical photospheric flows are susceptible to the Kelvin-Helmholtz instability (KHI), with rapid nonlinear growth times that are approximately half of a typical granule lifetime. The KHI produces vortical structures in intergranule lanes comparable to a typical fluxule radius; this is precisely the correct scale for maximum power transfer to the corona.

Karpen, Judith T.↗
NASA NTRS1993

An acoustic charge transport imager for high definition television applications

This report covers: (1) invention of a new, ultra-low noise, low operating voltage APD which is expected to offer far better performance than the existing volume doped APD device; (2) performance of a comprehensive series of experiments on the acoustic and piezoelectric properties of ZnO films sputtered on GaAs which can possibly lead to a decrease in the required rf drive power for ACT devices by 15dB; (3) development of an advanced, hydrodynamic, macroscopic simulator used for evaluating the performance of ACT and CTD devices and aiding in the development of the next generation of devices; (4) experimental development of CTD devices which utilize a p-doped top barrier demonstrating charge storage capacity and low leakage currents; (5) refinements in materials growth techniques and in situ controls to lower surface defect densities to record levels as well as increase material uniformity and quality.

Hunt, William D.↗
NASA NTRS1994

Tidal disruption of Periodic Comet Shoemaker-Levy 9 and a constraint on its mean density

The apparent tidal disruption of Periodic Comet Shoemaker-Levy 9 (1993e) during a close encounter within approximately 1.62 planetary radii of Jupiter can be used along with theoretical models of tidal disruption to place an upper bound on the density of the predisruption body. Depending on the theoretical model used, these upper bounds range from rho(sub c) less than 0.702 +/- 0.080 g/cu cm for a simple analytical model calibrated by numerical smoothed particle hydrodynamics (SPH) simulations to rho(sub c) less than 1.50 +/- 0.17 g/cu cm for a detailed semianalytical model. The quoted uncertainties stem from an assumed uncertainty in the perijove radius. However, the uncertainty introduced by the different theoretical models is the major source of error; this uncertainty could be eliminated by future SPH simulations specialized to cometary disruptions, including the effects of initially prolate, spinning comets. If the SPH-based upper bound turns out to be most appropriate, it would be consistent with the predisruption body being a comet with a relatively low density and porous structure, as has been asserted previously based on observations of cometary outgassing. Regardless of which upper bound is preferable, the models all agree that the predisruption body could not have been a relatively high-density body, such as an asteroid with rho approximately = 2 g/cu cm.

Boss, Alan P.↗
NASA NTRS1995

Cooling of solar flares plasmas. 1: Theoretical considerations

Theoretical models of the cooling of flare plasma are reexamined. By assuming that the cooling occurs in two separate phase where conduction and radiation, respectively, dominate, a simple analytic formula for the cooling time of a flare plasma is derived. Unlike earlier order-of-magnitude scalings, this result accounts for the effect of the evolution of the loop plasma parameters on the cooling time. When the conductive cooling leads to an 'evaporation' of chromospheric material, the cooling time scales L(exp 5/6)/p(exp 1/6), where the coronal phase (defined as the time maximum temperature). When the conductive cooling is static, the cooling time scales as L(exp 3/4)n(exp 1/4). In deriving these results, use was made of an important scaling law (T proportional to n(exp 2)) during the radiative cooling phase that was forst noted in one-dimensional hydrodynamic numerical simulations (Serio et al. 1991; Jakimiec et al. 1992). Our own simulations show that this result is restricted to approximately the radiative loss function of Rosner, Tucker, & Vaiana (1978). for different radiative loss functions, other scaling result, with T and n scaling almost linearly when the radiative loss falls off as T(exp -2). It is shown that these scaling laws are part of a class of analytic solutions developed by Antiocos (1980).

Cargill, Peter J.↗
NASA NTRS1996

Latitudinal Dependence of the Radial IMF Component: Coronal Imprint

Measurements by Ulysses have confirmed that there is no significant gradient with respect to heliomagnetic latitude in the radial component, B(sub r,) of the interplanetary magnetic field. In the corona, the plasma, beta is much less than 1, except directly above streamers, so longitudinal and latitudinal gradients in field strength will relax due to the transverse magnetic pressure gradient force as the solar wind carries magnetic flux away from the Sun. This happens quickly enough so that the field is essentially uniform by 5 - 10 solar radius, apparently remaining so as it is carried to beyond 1 AU. Here, we illustrate the coronal relaxation with a qualitative physical argument and by reference to a detailed Magneto HydroDynamics (MHD) simulation.

Suess, S. T.↗
NASA NTRS2010

3D Modeling of Forbidden Line Emission in the Binary Wind Interaction Region of Eta Carinae

We present recent work using three-dimensional (3D) Smoothed Particle Hydrodynamics (SPH) simulations to model the high ([Fe III], [Ar III], [Ne III] and [S III]) and low ([Fe II], [Ni II]) ionization forbidden emission lines observed in Eta Carinae using the HST/STIS. These structures are interpreted as the time-averaged, outer extensions of the primary wind and the wind-wind interaction region directly excited by the FUV of the hot companion star of this massive binary system. We discuss how analyzing the results of the 3D SPH simulations and synthetic slit spectra and comparing them to the spectra obtained with the HST/STIS helps us determine the absolute orientation of the binary orbit and helps remove the degeneracy inherent to models based solely on the observed RXTE X-ray light curve. A key point of this work is that spatially resolved observations like those with HST/STIS and comparison to 3D models are necessary to determine the alignment or misalignment of the orbital angular momentum axis with the Homunculus, or correspondingly, the alignment of the orbital plane with the Homunculus skirt.

Madura, Thomas↗
NASA NTRS2010

Constraining the Properties of the Eta Carinae System via 3-D SPH Models of Space-Based Observations: The Absolute Orientation of the Binary Orbit

The extremely massive (> 90 Solar Mass) and luminous (= 5 x 10(exp 6) Solar Luminosity) star Eta Carinae, with its spectacular bipolar "Homunculus" nebula, comprises one of the most remarkable and intensely observed stellar systems in the galaxy. However, many of its underlying physical parameters remain a mystery. Multiwavelength variations observed to occur every 5.54 years are interpreted as being due to the collision of a massive wind from the primary star with the fast, less dense wind of a hot companion star in a highly elliptical (e approx. 0.9) orbit. Using three-dimensional (3-D) Smoothed Particle Hydrodynamics (SPH) simulations of the binary wind-wind collision in Eta Car, together with radiative transfer codes, we compute synthetic spectral images of [Fe III] emission line structures and compare them to existing Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS) observations. We are thus able, for the first time, to constrain the absolute orientation of the binary orbit on the sky. An orbit with an inclination of i approx. 40deg, an argument of periapsis omega approx. 255deg, and a projected orbital axis with a position angle of approx. 312deg east of north provides the best fit to the observations, implying that the orbital axis is closely aligned in 3-1) space with the Homunculus symmetry axis, and that the companion star orbits clockwise on the sky relative to the primary.

Madura, Thomas I.↗
NASA NTRS2011

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

Constraining the Properties of the Eta Carinae System via 3-D SPH Models of Space-Based Observations: The Absolute Orientation of the Binary Orbit

The extremely massive (> 90 Stellar Mass) and luminous (= 5 x 10(exp 6) Stellar Luminosity) star Eta Carinae, with its spectacular bipolar "Homunculus" nebula, comprises one of the most remarkable and intensely observed stellar systems in the Galaxy. However, many of its underlying physical parameters remain unknown. Multiwavelength variations observed to occur every 5.54 years are interpreted as being due to the collision of a massive wind from the primary star with the fast, less dense wind of a hot companion star in a highly elliptical (e approx. 0.9) orbit. Using three-dimensional (3-D) Smoothed Particle Hydrodynamics (SPH) simulations of the binary wind-wind collision, together with radiative transfer codes, we compute synthetic spectral images of [Fe III] emission line structures and compare them to existing Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS) observations. We are thus able, for the first time, to tightly constrain the absolute orientation of the binary orbit on the sky. An orbit with an inclination of approx. 40deg, an argument of periapsis omega approx. 255deg, and a projected orbital axis with a position angle of approx. 312deg east of north provides the best fit to the observations, implying that the orbital axis is closely aligned in 3-D space with the Homunculus symmetry axis, and that the companion star orbits clockwise on the sky relative to the primary.

Madura, Thomas I.↗
NASA NTRS2013

Effects of the LBV Primary's Mass-loss Rate on the 3D Hydrodynamics of eta Carinae's Colliding Winds

At the heart of eta Carinae's spectacular "Homunculus" nebula lies an extremely luminous (L(sub Total) greater than approximately 5 × 10(exp 6) solar luminosity) colliding wind binary with a highly eccentric (e approximately 0.9), 5.54-year orbit (Figure 1). The primary of the system, a Luminous Blue Variable (LBV), is our closest (D approximately 2.3 kpc) and best example of a pre-hypernova or pre-gamma ray burst environment. The remarkably consistent and periodic RXTE X-ray light curve surprisingly showed a major change during the system's last periastron in 2009, with the X-ray minimum being approximately 50% shorter than the minima of the previous two cycles1. Between 1998 and 2011, the strengths of various broad stellar wind emission lines (e.g. Halpha, Fe II) in line-of-sight (l.o.s.) also decreased by factors of 1.5 - 3 relative to the continuum2. The current interpretation for these changes is that they are due to a gradual factor of 2 - 4 drop in the primary's mass-loss rate over the last approximately 15 years1, 2. However, while a secular change is seen for a direct view of the central source, little to no change is seen in profiles at high stellar latitudes or reflected off of the dense, circumbinary material known as the "Weigelt blobs"2, 3. Moreover, model spectra generated with CMFGEN predict that a factor of 2 - 4 drop in the primary's mass-loss rate should lead to huge changes in the observed spectrum, which thus far have not been seen. Here we present results from large- (plus or minus 1620 AU) and small- (plus or minus 162 AU) domain, full 3D smoothed particle hydrodynamics (SPH) simulations of eta Car's massive binary colliding winds for three different primary-star mass-loss rates (2.4, 4.8, and 8.5 × 10(exp -4) solar mass/yr). The goal is to investigate how the mass-loss rate affects the 3D geometry and dynamics of eta Car's optically-thick wind and spatially-extended wind-wind collision (WWC) regions, both of which are known sources of observed X-ray, optical, UV, and near-IR emission and absorption. We use two domain sizes in order to better understand how the primary's mass-loss rate influences the various observables that form at different length scales. The 3D simulations provide information important for helping constrain η Car's recent mass-loss history and future state.

lavestigate↗