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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Stellar migration and chemical enrichment in the milky way disc: a hybrid model

ABSTRACT We develop a hybrid model of galactic chemical evolution that combines a multiring computation of chemical enrichment with a prescription for stellar migration and the vertical distribution of stellar populations informed by a cosmological hydrodynamic disc galaxy simulation. Our fiducial model adopts empirically motivated forms of the star formation law and star formation history, with a gradient in outflow mass loading tuned to reproduce the observed metallicity gradient. With this approach, the model reproduces many of the striking qualitative features of the Milky Way disc’s abundance structure: (i) the dependence of the [O/Fe]–[Fe/H] distribution on radius Rgal and mid-plane distance |z|; (ii) the changing shapes of the [O/H] and [Fe/H] distributions with Rgal and |z|; (iii) a broad distribution of [O/Fe] at sub-solar metallicity and changes in the [O/Fe] distribution with Rgal, |z|, and [Fe/H]; (iv) a tight correlation between [O/Fe] and stellar age for [O/Fe] > 0.1; (v) a population of young and intermediate-age α-enhanced stars caused by migration-induced variability in the Type Ia supernova rate; (vi) non-monotonic age–[O/H] and age–[Fe/H] relations, with large scatter and a median age of ∼4 Gyr near solar metallicity. Observationally motivated models with an enhanced star formation rate ∼2 Gyr ago improve agreement with the observed age–[Fe/H] and age–[O/H] relations, but worsen agreement with the observed age–[O/Fe] relation. None of our models predict an [O/Fe] distribution with the distinct bimodality seen in the observations, suggesting that more dramatic evolutionary pathways are required. All code and tables used for our models are publicly available through the Versatile Integrator for Chemical Evolution (VICE; https://pypi.org/project/vice).

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

The distribution of [α/Fe] in the Milky Way disc

ABSTRACT Using a sample of red giant stars from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) Data Release 16, we infer the conditional distribution $p([\alpha /{\rm Fe}]\, |\, [{\rm Fe}/{\rm H}])$ in the Milky Way disk for the α-elements Mg, O, Si, S, and Ca. In each bin of [Fe/H] and Galactocentric radius R, we model p([α/Fe]) as a sum of two Gaussians, representing ‘low-α’ and ‘high-α’ populations with scale heights $z_1=0.45\, {\rm kpc}$ and $z_2=0.95\, {\rm kpc}$, respectively. By accounting for age-dependent and z-dependent selection effects in APOGEE, we infer the [α/Fe] distributions that would be found for a fair sample of long-lived stars covering all z. Near the Solar circle, this distribution is bimodal at sub-solar [Fe/H], with the low-α and high-α peaks clearly separated by a minimum at intermediate [α/Fe]. In agreement with previous results, we find that the high-α population is more prominent at smaller R, lower [Fe/H], and larger |z|, and that the sequence separation is smaller for Si and Ca than for Mg, O, and S. We find significant intrinsic scatter in [α/Fe] at fixed [Fe/H] for both the low-α and high-α populations, typically ∼0.04-dex. The means, dispersions, and relative amplitudes of this two-Gaussian description, and the dependence of these parameters on R, [Fe/H], and α-element, provide a quantitative target for chemical evolution models and a test for hydrodynamic simulations of disk galaxy formation. We argue that explaining the observed bimodality will probably require one or more sharp transitions in the disk’s gas accretion, star formation, or outflow history in addition to radial mixing of stellar populations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

The hTERT core promoter forms three parallel G-quadruplexes

The structure of the 68 nt sequence with G-quadruplex forming potential within the hTERT promoter is disputed. One model features a structure with three stacked parallel G-quadruplex units, while another features an unusual duplex hairpin structure adjoined to two stacked parallel and antiparallel quadruplexes. We report here the results of an integrated structural biology study designed to distinguish between these possibilities. As part of our study, we designed a sequence with an optimized hairpin structure and show that its biophysical and biochemical properties are inconsistent with the structure formed by the hTERT wild-type sequence. By using circular dichroism, thermal denaturation, nuclear magnetic resonance spectroscopy, analytical ultracentrifugation, small-angle X-ray scattering, molecular dynamics simulations and a DNase I cleavage assay we found that the wild type hTERT core promoter folds into a stacked, three-parallel G-quadruplex structure. The hairpin structure is inconsistent with all of our experimental data obtained with the wild-type sequence. All-atom models for both structures were constructed using molecular dynamics simulations. These models accurately predicted the experimental hydrodynamic properties measured for each structure. We found with certainty that the wild-type hTERT promoter sequence does not form a hairpin structure in solution, but rather folds into a compact stacked three-G-quadruplex conformation.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Multiplicity dependence of $\pi$, K, and p production in pp collisions at $\sqrt{s}$ = 13 TeV

This paper presents the measurements of π ± , K ± , p and $\bar{p}$ transverse momentum (p T ) spectra as a function of charged-particle multiplicity density in proton–proton (pp) collisions at $\sqrt{s}$ = 13TeV with the ALICE detector at the LHC. Such study allows us to isolate the center-of-mass energy dependence of light-flavour particle production. The measurements reported here cover a p T range from 0.1 to 20 GeV / c and are done in the rapidity interval | y| < 0.5. The p T -differential particle ratios exhibit an evolution with multiplicity, similar to that observed in pp collisions at $\sqrt{s}$ = 7 TeV, which is qualitatively described by some of the hydrodynamical and pQCD-inspired models discussed in this paper. Furthermore, the p T -integrated hadron-to-pion yield ratios measured in pp collisions at two different center-of-mass energies are consistent when compared at similar multiplicities. This also extends to strange and multi-strange hadrons, suggesting that, at LHC energies, particle hadrochemistry scales with particle multiplicity the same way under different collision energies and colliding systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Benchmark Solutions for Radiation Transport in Stochastic Media with Inhomogeneous Material Statistics

Accurately solving implicit Monte Carlo (IMC) thermal photon transport problems with mixed material cells is important in realistic applications. The production IMC package at LLNL treats mixed material cells arising from ALE remap and hydrodynamics using the same approximate model. The new Imp IMC thermal photon transport package currently under development has both a material interface reconstruction (MIR) algorithm and a Levermore-Pomraning (LP) stochastic medium algorithm for treating mixed material cells. Existing stochastic medium algorithms for treating mixed material cells in IMC lack a complete theoretical basis. The IMC LP algorithm implementation has been demonstrated to reproduce published deterministic LP solutions for the particular case of spatially homogeneous material statistics. Realistic simulations will include spatially inhomogeneous material statistics (material mean chord lengths). In a previous investigation, the LP-model for transport in binary stochastic media in rod geometry was generalized to accommodate spatially varying material chord lengths, i.e., the mixing statistics were allowed to be nonhomogeneous. Analytical solutions were obtained and used to produce a verifi cation suite for the Imp IMC Levermore-Pomraning implementation for different spatial variations of the chord lengths. However, the accuracy of the LP model when the mixing statistics are nonhomogeneous has not been assessed and leaves open the question of whether local accuracy is improved or further degraded when chord lengths are not uniform. This shortcoming is rectifi ed here by developing benchmark analytic solutions for transport in binary Markovian stochastic mixtures in rod geometry with nonhomogeneous mixing statistics, using spatially varying chord lengths considered in the previous investigation based on the LP model. Methods for sampling a nonhomogeneous Poisson process (NHPP) are first described and used to construct individual realizations of the binary mixtures in rod geometry. Analytic solutions are then obtained for the forward and backward directed fluxes on a given realization, now viewed as a deterministic medium with alternating layers of the two materials with known interface locations. Finally, material averaged scalar fluxes are obtained using these sampling schemes with spatially linear and quadratic chord lengths and used to assess the accuracy of the previously obtained LP-model results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1986

Hydrodynamics of internal solitons and a comparison of SIR-A and SIR-B data with ocean measurements

Large internal solitary waves have been observed by Shuttle SIR-A and SIR-B at locations in the Andaman Sea and the New York Bight. Satellite imagery and oceanographic measurements are used in conjunction with hydrodynamic interaction and electromagnetic scattering models to estimate the expected SAR image intensity modulations associated with the internal waves. There is reasonable agreement between the predicted and observed internal wave signatures.

Apel, J. R.↗

NASA NTRS1987

Electromagnetic jets from compact objects

The possibility that at least some astrophysical jets are initially electromagnetic in origin is examined. Subsequent pick-up of ionization would convert such electromagnetic jets into hydrodynamic jets. In such a model, relativistic outflow is formed into highly collimated beams simply through the interaction with the surrounding medium. Forming jets under such general circumstances is encouraging in view of the range of scales that appear to be involved. The overall properties of such jets are largely determined by a single dimensionless parameter: the characteristic electrostatic potential drop rewritten as a particle Lorentz factor. Consequently, the determination of any one observable, such as the total power output, also determines the particle energy scale, the electromagnetic field strengths, etc.

Michel, F. Curtis↗

NASA NTRS1987

Comparison of horizontal difference schemes for the shallow water equations on a sphere

The accuracy of horizontal difference schemes used in the hydrodynamics parts of General Circulation Models are compared by means of numerical experiments for the shallow water equations on a sphere. As expected, the phase lag of moving waves decreases as the order of accuracy of a scheme increases or as the grid resolution increases. Overall, Takano and Wurtele's partial fourth order energy and potential enstrophy conserving scheme on the C grid is most accurate. It is clearly superior to the other schemes for the Rossby-Haurwitz wave number 6 initial conditions for coarse grid resolution.

Russell, Gary L.↗

NASA NTRS1989

Superbubble blowout dynamics

Multiple supernovae and stellar winds from OB associations carve large holes filled with hot gas in the galactic disk. These superbubbles sweep up H I into cold, thin, dense shells and eventually grow large enough to blow completely out of the galactic H I disk. When superbubbles blow out of the disk, they vent hot gas and supernova energy into the galactic corona. In this paper ZEUS, a two-dimensional hydrodynamics code, is used to model the blowout of a superbubble from exponential and Gaussian models for the vertical density stratification. The results are compared to those from the Kompaneets (thin-shell) approximation. It is found that this approximation works very well, and that most of the mass of the shell remains in the plane, with 5 percent of it accelerating upward. The venting of the hot gas and the stability of the shell depends strongly on the model of the density distribution. It is suggested that the low galactic halo actually consists of a froth of merged superbubbles.

Mac Low, Mordecai-Mark↗

NASA NTRS1990

X-rays from colliding stellar winds

A stellar wind from a massive OB or Wolf-Rayet star in a binary system will strike the surface or stellar wind of its companion, forming shocked gas that can radiate X-rays. The X-ray spectrum from the shocked winds will vary in a predictable way with orbital phase, owing to photoelectric absorption by the stellar winds. Detailed models are calculated for the hydrodynamics and X-ray emission from two such systems. In one of these systems (HD 165052), the winds are nearly identical in strength. In the other (V444 Cygni), the wind of the Wolf-Rayet star overwhelms and crushes that of its companion. The calculated X-ray luminosities agree fairly well with the observed values for HD 165052 and for V444 Cygni. These results can be scaled to other such systems.

Luo, Ding↗

NASA NTRS1991

Collisions of giant stars with compact objects - Hydrodynamical calculations

Smooth particle hydrodynamics are used here to model encounters between an 0.8 solar mass giant star and a 1.4 solar mass compact object on a parabolic orbit. It is found that all encounters with periastron distance r(p) less than about 2.5 stellar radii R(G) lead to significant disruption of the stellar envelope, including some mass loss. When r(p)/R(G) is less than about 0.5, the envelope is completely destroyed during the collision. About 40 percent of the gas escapes, while the rest forms a massive disk around the compact object, leaving behind the bare stellar core on a highly eccentric orbit. Subsequent close passages lead to continuous extraction of gas from the disk at a large rate, and the binary orbit may eventually stabilize while still retaining a large eccentricity. It is suggested that many globular cluster pulsars may have been spun up by accretion from a massive disk rather than from a binary companion, thereby providing a possible solution to the birthrate problem.

Rasio, Frederic A.↗

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

Models of the Cartwheel ring galaxy: Spokes and starbursts

Recent observations of this famous ring galaxy, including optical and near-infrared CCD surface photometry, and VLA radio continuum and 21 cm line mapping (Higdon 1992b, in prep.), have inspired a renewed modeling effort. Toomre's (1978, in The Large-scale Structure of the Universe, eds. Longair and Einasto) series of restricted three-body simulations demonstrated how the multiple rings could be produced in a nearly head-on galaxy collision. New models with a halo-dominated potential based on the 21 cm rotation curve are able to reproduce such details as the spacing between rings, ring widths, offset of the nucleus, and several kinematical features, thus providing strong support for the collisional theory. The new observations have shown there are little or no old stars in Cartwheel; it may consist almost entirely of gas and stars produced as a result of compression in the ring wave. To model this process Smooth Particle Hydrodynamics (SPH) simulations of the Cartwheel disk have been performed. Fixed gravitational potentials were used to represent the Cartwheel and a roughly 30 percent mass collision partner. The interaction dynamics was treated as in the usual restricted three-body approximation, and the effects of local self-gravity between disk particles were calculated. We are particularly interested in testing the theory that enhanced star formation in waves is the result of gravitational instability in the compressed region (see e.g. Kennicutt 1989, ApJ 344, 685). The gas surface density in a number of simulations was initialized to a value slightly below the threshold for local gravitational instability throughout most of the disk. The first ring wave produces relatively modest compressions (a factor of order a few), triggering instability in a narrow range of wavelengths. Self-gravity in the disk is calculated over a comparable range of scales. Simulations were run with isothermal, adiabatic, and adiabatic with radiative cooling characterized by a relatively short timescale. The isothermal approximation is good except in the vicinity of the strong second (inner) ring, and several snapshots from one case are shown in the figure below. Flocculent spiral segments are present before the collision, and these are compressed into dense knots in the ring wave. These knots are likely to be sites of vigorous star formation. In the strong rarefaction behind the outer ring most of the knots are radially stretched and sheared, giving rise to spoke-like features. A few dense knots are evidently very tightly bound, because they retain their coherence and are stretched relatively little through the rarefaction. This is in accord with evidence for continuing star formation in some spokes (Marcum, Appleton and Higdon 1992). The number and spacing of spokes is a direct function of the scale of the gravitational instability in the disk. Thus, the gravitational instability theory, together with the hypothesis that massive stars are only formed in dense knots of gas, can account for most of the distinct morphology of the Cartwheel.

Struck-Marcell, Curtis↗

NASA NTRS1992

Origin of noble gases in the terrestrial planets

Current models of the origin of noble gases in the terrestrial planets are reviewed. Primary solar system volatile sources and processes are examined along with the current data base on noble gases and its applications to evolutionary processing. Models of atmospheric evolution by hydrodynamic escape are addressed.

Pepin, Robert O.↗

NASA NTRS1994

Static and dynamic force/moment measurements in the Eidetics water tunnel

Water tunnels have been utilized in one form or another to explore fluid mechanics and aerodynamics phenomena since the days of Leonardo da Vinci. Water tunnel testing is attractive because of the relatively low cost and quick turn-around time to perform flow visualization experiments and evaluate the results. The principal limitation of a water tunnel is that the low flow speed, which provides for detailed visualization, also results in very small hydrodynamic (aerodynamic) forces on the model, which, in the past, have proven to be difficult to measure accurately. However, the advent of semi-conductor strain gage technology and devices associated with data acquisition such as low-noise amplifiers, electronic filters, and digital recording have made accurate measurements of very low strain levels feasible. The principal objective of this research effort was to develop a multi-component strain gage balance to measure forces and moments on models tested in flow visualization water tunnels. A balance was designed that allows measuring normal and side forces, and pitching, yawing and rolling moments (no axial force). The balance mounts internally in the model and is used in a manner typical of wind tunnel balances. The key differences between a water tunnel balance and a wind tunnel balance are the requirement for very high sensitivity since the loads are very low (typical normal force is 0.2 lbs), the need for water proofing the gage elements, and the small size required to fit into typical water tunnel models.

Carlos J Suarez↗

NASA NTRS1996

Interplanetary and Interstellar Dust Particles: Reentry Heating and Capture in Aerogel

Studies of interplanetary and interstellar dust can provide significant information on the evolution of the solar system or stars, respectively. However, for reliable analysis it is crucial to know how the particles have been modified during reentry (in the case of interplanetary dust particles, or [DP's) and impact into the capture medium. In the case of stratospheric capture, particles will be heated by atmospheric drag. Subsequent capture of the particles will result in heating, ablation, accretion of the capture medium and possible fragmentation. Modeling of these processes is a useful way of assessing their effects on the interpretation of the compositional data for these particles. Previous work on reentry heating has shown that heat diffusion alone cannot adequately account for temperature gradients observed in IDP's. In fact, for any reasonable thermal parameters, calculations show the particles to be nearly isothermal. Here we extend those studies to include phase transitions. These preliminary results are promising and show a significant temperature gradient for a 100 micron (diameter) particle. The actual capture of the particles in silica aerogel is being modeled using a comprehensive shock hydrodynamic code (called CTH). Various options of this code were investigated to attempt to make the most appropriate choice of methods of impact, equations of state, and processes of energy transfer from capture material to particle. The initial calculations with the code used only 'reasonable' estimates for the physical parameters of silica aerogel. Through the literature searches and personal contacts with the knowledgeable scientists, the best possible mechanical and thermal data have been made available for these computations.

Szydlik, Paul P.↗

NASA NTRS1997

Quantitative UV spectroscopy of early O stars on the Magellanic Clouds: The determination of the stellar metallicities

Ultraviolet spectra of four O stars in the Magellanic Clouds obtained with the faint object spectrograph of the Hubble Space Telescope are analyzed with respect to their metallicity. The metal abundances are derived from the stellar parameters and the mass loss rate with a two step procedure: hydrodynamic radiation-driven wind models with metallicity as a free parameter are constructed to fit the observed wind momentum rate and thus yield a dynamical metallicity, and synthetic spectra are computed for different metal abundances and compared to the observed spectra in order to obtain a spectroscopic metallicity.

Haser, Stefan M.↗

NASA NTRS2003

Studies of Solar Wind Interaction and Ionospheric Processes at Venus and Mars

This is the final report summarizing the work done during the last three years under NASA Grant NAG5-8946. Our efforts centered on a systematic development of a new generation of three dimensional magneto-hydrodynamic (MHD) numerical code, which models the interaction processes of the solar wind or fast flowing magnetospheric plasma with 'non-magnetic' solar system bodies (e.g. Venus, Mars, Europa, Titan). We have also worked on a number of different, more specific and discrete studies, as various opportunities arose. In the next few pages we briefly summarize these efforts.

Bogan, Denis↗