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 973 records · Page 54

Violent mergers can explain the inflated state of some of the fastest stars in the Galaxy

A significant number of hypervelocity stars with velocities between 1500 − 2500 km s −1 have recently been observed. The only plausible explanation so far is that they were produced through thermonuclear supernovae in white dwarf binaries. Since these stars are thought to be surviving donors of Type Ia supernovae, a surprising finding was that these stars are inflated, with radii an order of magnitude higher than expected for Roche-lobe-filling donors. Recent attempts at explaining them have combined 3D hydrodynamical supernova explosion simulations with 1D stellar modelling to explain the impact of supernova shocks on runaway white dwarfs. However, only the hottest and most compact of those runaway stars can so far marginally be reproduced by detailed models of runaways from supernova explosions. In this and a companion paper, we introduce a new AREPO simulation of two massive CO white dwarfs that explode via a violent merger. During the merger, the primary white dwarf ignites when the secondary is on its last orbit and plunging towards the primary. In the corresponding aftermath, the core of the secondary white dwarf of 0.16 M ⊙ remains bound, moving at a velocity of ∼2800 km s −1 . We mapped this object into MESA and show that this runaway star can explain the observations of two hypervelocity stars that were dubbed D6-1 and D6-3 based on their original discovery motivated by the D6 scenario, though the violent merger scenario presented here is somewhat distinct from the D6 scenario.

Astronomy and AstroPhysics↗

Resonance production in Pb+Pb collisions at 5.02 TeV

The yields, mean transverse momenta, and flow of K* 0 , ρ 0 , Λ(1520) resonances provide an evidence of a late stage hadronic rescattering in ultrarelativistic central heavy ion collisions [1]. Using hydrodynamic + hadronic afterburner simulations of Pb+Pb collisions at 5.02 TeV we achieve a reasonable description of resonance yields and spectra as a function of collision centrality. We demonstrate that the measurements of Λ(1520)’s mean transverse momentum allow to constrain the unknown branching ratios of Σ* → Λ(1520)π decays. Hadronic dynamics leads to an enhanced Δ(1232) production in central collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The mass and galaxy distribution around SZ-selected clusters

ABSTRACT We present measurements of the radial profiles of the mass and galaxy number density around Sunyaev–Zel’dovich (SZ)-selected clusters using both weak lensing and galaxy counts. The clusters are selected from the Atacama Cosmology Telescope Data Release 5 and the galaxies from the Dark Energy Survey Year 3 data set. With signal-to-noise ratio of 62 (45) for galaxy (weak lensing) profiles over scales of about 0.2–20 h−1 Mpc, these are the highest precision measurements for SZ-selected clusters to date. Because SZ selection closely approximates mass selection, these measurements enable several tests of theoretical models of the mass and light distribution around clusters. Our main findings are: (1) The splashback feature is detected at a consistent location in both the mass and galaxy profiles and its location is consistent with predictions of cold dark matter N-body simulations. (2) The full mass profile is also consistent with the simulations. (3) The shapes of the galaxy and lensing profiles are remarkably similar for our sample over the entire range of scales, from well inside the cluster halo to the quasilinear regime. We measure the dependence of the profile shapes on the galaxy sample, redshift, and cluster mass. We extend the Diemer & Kravtsov model for the cluster profiles to the linear regime using perturbation theory and show that it provides a good match to the measured profiles. We also compare the measured profiles to predictions of the standard halo model and simulations that include hydrodynamics. Applications of these results to cluster mass estimation, cosmology, and astrophysics are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

Low-energy 23 Al β-delayed proton decay and 22 Na destruction in novae

The radionuclide 22 Na is a target of γ-ray astronomy searches, predicted to be produced during thermonuclear runaways driving classical novae. The 22 Na (p, γ) 23 Mg reaction is the main destruction channel of 22 Na during a nova, hence, its rate is needed to accurately predict the 22 Na yield. However, experimental determinations of the resonance strengths have led to inconsistent results. In this Rapid Communication, we report a measurement of the branching ratios of the 23 Al β-delayed protons as a probe of the key 204-keV (center-of-mass) 22 Na (p, γ) 23 Mg resonance strength. We report a factor of 5 lower branching ratio compared to the most recent literature value. The variation in 22 Na yield due to nuclear data inconsistencies was assessed using a series of hydrodynamic nova outburst simulations and has increased to a factor of 3.8, corresponding to a factor of ≈ 2 uncertainty in the maximum detectability distance. Finally, this is the first reported scientific measurement using the Gaseous Detector with Germanium Tagging system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

athenapk

The LANL fork of the AthenaPK project, which contains curvlinear coordinates The Code is a set of modifications to AthenaPK, an existing open source 3D hydrodynamics and magnetohydrodynamics simulation software for astrophysics, which enable simulation with uniform curvilinear geometries. Specifically it extends the existing Cartesian coordinates to uniform spherical coordinate and uniform cylindrical coordinate systems

Miller, Jonah [@LANL]↗

REYNO: A Reactive Hydrodynamics Modeling Suite

REYNO is a Python library that is primarily designed as a platform for the agile development and testing of novel reactive flow models. Equations of state and multi-step rate laws can be implemented using the provided class hierarchy with different thermodynamic closure conditions. The library includes a one-dimensional Lagrangian hydrodynamic solver to simulate multi-material reactive or inert problems in planar, cylindrical, or spherical coordinates. Custom time-dependent boundary conditions are also available for problems such as ramp compression. This document describes the different numerical methods and algorithms that are included in REYNO.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Catastrophic Cooling in Superwinds. III. Nonequilibrium Photoionization

Observations of some starburst-driven galactic superwinds suggest that strong radiative cooling could play a key role in the nature of feedback and the formation of stars and molecular gas in star-forming galaxies. These catastrophically cooling superwinds are not adequately described by adiabatic fluid models, but they can be reproduced by incorporating nonequilibrium radiative cooling functions into the fluid model. In this work, we have employed the atomic and cooling module maihem implemented in the framework of the flash hydrodynamics code to simulate the formation of radiatively cooling superwinds as well as their corresponding nonequilibrium ionization (NEI) states for various outflow parameters, gas metallicities, and ambient densities. We employ the photoionization program cloudy to predict radiation- and density-bounded photoionization for these radiatively cooling superwinds, and we predict UV and optical line emission. Our nonequilibrium photoionization models built with the NEI states demonstrate the enhancement of C iv, especially in metal-rich, catastrophically cooling outflows, and O vi in metal-poor ones.

Cooling flows↗

MeV–GeV Gamma-Ray Emission from SNR G327.1–1.1 Discovered by the Fermi-LAT

Abstract We report the discovery of MeV–GeV γ -ray emission by the Fermi-LAT positionally coincident with the TeV pulsar wind nebula (PWN) HESS J1554–550 within the host supernova remnant (SNR) G327.1–1.1. The γ -ray emission is point-like and faint but significant (>4 σ ) in the 300 MeV–2 TeV energy range. We report here the Fermi-LAT analysis of the observed γ -ray emission followed by a detailed multiwavelength investigation to understand the nature of the emission. The central pulsar powering the PWN within G327.1–1.1 has not been detected in any wave band; however, it is likely embedded within the X-ray nebula, which is displaced from the center of the radio nebula. The γ -ray emission is faint and therefore a pulsation search to determine if the pulsar may be contributing is not feasible. Prior detailed multiwavelength reports revealed an SNR system that is old, τ ∼ 18,000 yr, where the interaction of the reverse shock with the PWN is underway or has recently occurred. We find that the γ -ray emission agrees remarkably well with a detailed broadband model constructed in a prior report based on independent hydrodynamical and semianalytic simulations of an evolved PWN. We further investigate the physical implications of the model for the PWN evolutionary stage incorporating the new Fermi-LAT data and attempt to model the distinct particle components based on a spatial separation analysis of the displaced PWN counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗

Theory and modeling of atmospheric turbulence, part 2

Two dimensional geostrophic turbulence driven by a random force is investigated. Based on the Liouville equation, which simulates the primitive hydrodynamical equations, a group-kinetic theory of turbulence is developed and the kinetic equation of the scaled singlet distribution is derived. The kinetic equation is transformed into an equation of spectral balance in the equilibrium and non-equilibrium states. Comparison is made between the propagators and the Green's functions in the case of the non-asymptotic quasi-linear equation to prove the equivalence of both kinds of approximations used to describe perturbed trajectories of plasma turbulence. The microdynamical state of fluid turbulence is described by a hydrodynamical system and transformed into a master equation analogous to the Vlasov equation for plasma turbulence. The spectral balance for the velocity fluctuations of individual components shows that the scaled pressure strain correlation and the cascade transfer are two transport functions that play the most important roles.

Chen, C. M.↗

Liquid oxygen sloshing in Space Shuttle External Tank

This paper describes a numerical simulation of the hydrodynamics within the liquid oxygen tank of the Space Shuttle External Tank during liftoff. Before liftoff, the tank is filled with liquid oxygen (LOX) to approximately 97 percent with the other 3 percent containing gaseous oxygen (GOX) and helium. During liftoff, LOX is drained from the bottom of the tank, and GOX is pumped into the tank's ullage volume. There is a delay of several seconds before the GOX reaches the tank which causes the ullage pressure to decrease for several seconds after liftoff; this pressure 'slump' is a common phenomenon in rocket propulsion. When four slosh baffles were removed from the tank, the ullage gas pressure dropped more rapidly than in all previous flights. The purpose of this analysis was to determine whether the removal of the baffles could have caused the increased pressure 'slump' by changing the LOX surface dynamics. The results show that the LOX surface undergoes very high vertical accelerations (up to 5 g) and, therefore, splashing almost certainly occurs. The number of baffles does not affect the surface if the structural motion is assumed; but, the number of baffles may affect the structural motion of the tank.

Kannapel, M. D.↗

A prediction of the gamma-ray flux from Nova Herculis 1991

UV observations made with the IUE satellite of Nova Her 1991 are presented and discussed. A distance of 3.4 kpc and reddening of 0.6 E(B-V) to Nova Her 1991 are determined by comparing its UV and optical outburst behavior with those of the extragalactic ONeMg nova LMC 1990 No. 1. Her 1991 was also detected by Rosat five days into its outburst, and it is argued that a small fraction of the X-rays may be the product of Comptonized 1.2 and 0.5 MeV emission from the decay of Na-22. If the mass of ejected shell is indeed as high as was estimated by Woodward et al. (1992), namely 0.0001 solar mass, then the present results predicted that Na-22 gamma-ray emission from the ejecta of Nova Her 1991 would still be detectable by GRO in the spring of 1992. The detection level was expected to be about 5 sigma in a two-week observation. Such a detection would strongly guide and constrain the hydrodynamic and nucleosynthesis simulations of nova outbursts on massive white dwarfs.

Starrfield, Sumner↗

Effects of convection electric field on upwelling and escape of ionospheric O(+)

A Monte Carlo code is used to explore the full effects of the convection electric field on distributions of upflowing O(+) ions from the cusp/cleft ionosphere. Trajectories of individual ions/neutrals are computed as they undergo multiple charge-exchange collisions. In the ion state, the trajectories are computed in realistic models of the magnetic field and the convection, corotation, and ambipolar electric fields. The effects of ion-ion collisions are included, and the trajectories are computed with and without simultaneous stochastic heating perpendicular to the magnetic field by a realistic model of broadband, low frequency waves. In the neutral state, ballistic trajectories in the gravitational field are computed. The initial conditions of the ions, in addition to ambipolar electric field and the number densities and temperatures of O(+), H(+), and electrons as a function of height in the cusp/cleft region were obtained from the results of Gombosi and Killeen (1987), who used a hydrodynamic code to simulate the time-dependent frictional-heating effects in a magnetic tube during its motion though the convection throat. The distribution of the ion fluxes as a function of height are constructed from the case histories.

Cladis, J. B.↗

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↗

LDEF data: Comparisons with existing models

The relationship between the observed cratering impact damage on the Long Duration Exposure Facility (LDEF) versus the existing models for both the natural environment of micrometeoroids and the man-made debris was investigated. Experimental data was provided by several LDEF Principal Investigators, Meteoroid and Debris Special Investigation Group (M&D SIG) members, and by the Kennedy Space Center Analysis Team (KSC A-Team) members. These data were collected from various aluminum materials around the LDEF satellite. A PC (personal computer) computer program, SPENV, was written which incorporates the existing models of the Low Earth Orbit (LEO) environment. This program calculates the expected number of impacts per unit area as functions of altitude, orbital inclination, time in orbit, and direction of the spacecraft surface relative to the velocity vector, for both micrometeoroids and man-made debris. Since both particle models are couched in terms of impact fluxes versus impactor particle size, and much of the LDEF data is in the form of crater production rates, scaling laws have been used to relate the two. Also many hydrodynamic impact computer simulations were conducted, using CTH, of various impact events, that identified certain modes of response, including simple metallic target cratering, perforations and delamination effects of coatings.

Coombs, Cassandra R.↗

Marshall Space Flight Center's Impact Testing Facility Capabilities

Marshall Space Flight Center s (MSFC) Impact Testing Facility (ITF) serves as an important installation for space and missile related materials science research. The ITF was established and began its research in spacecraft debris shielding in the early 1960s, then played a major role in the International Space Station debris shield development. As NASA became more interested in launch debris and in-flight impact concerns, the ITF grew to include research in a variety of impact genres. Collaborative partnerships with the DoD led to a wider range of impact capabilities being relocated to MSFC as a result of the closure of Particle Impact Facilities in Santa Barbara, California. The Particle Impact Facility had a 30 year history in providing evaluations of aerospace materials and components during flights through rain, ice, and solid particle environments at subsonic through hypersonic velocities. The facility s unique capabilities were deemed a "National Asset" by the DoD. The ITF now has capabilities including environmental, ballistic, and hypervelocity impact testing utilizing an array of air, powder, and two-stage light gas guns to accommodate a variety of projectile and target types and sizes. Numerous upgrades including new instrumentation, triggering circuitry, high speed photography, and optimized sabot designs have been implemented. Other recent research has included rain drop demise characterization tests to obtain data for inclusion in on-going model development. The current and proposed ITF capabilities range from rain to micrometeoroids allowing the widest test parameter range possible for materials investigations in support of space, atmospheric, and ground environments. These test capabilities including hydrometeor, single/multi-particle, ballistic gas guns, exploding wire gun, and light gas guns combined with Smooth Particle Hydrodynamics Code (SPHC) simulations represent the widest range of impact test capabilities in the country.

Finchum, Andy↗

Constraints on the Nature of CID-42: Recoil Kick or Supermassive Black Hole Pair?

The galaxy CXOC J100043.1+020637, also known as CID-42, is a highly unusual object. An apparent galaxy merger remnant, it displays signatures of both an inspiraling, kiloparsecscale active galactic nucleus (AGN) pair and of a recoiling AGN with a kick velocity approximately greater than 1300 km s(exp -1). Among recoiling AGN candidates, CID-42 alone has both spatial offsets (in optical and X-ray bands) and spectroscopic offsets. In order to constrain the relative likelihood of both scenarios, we develop models using hydrodynamic galaxy merger simulations coupled with radiative transfer calculations. Our gas-rich, major merger models are generally well matched to the galactic morphology and to the inferred stellar mass and star formation rate. We show that a recoiling supermassive black hole (SMBH) in CID-42 should be observable as an AGN at the time of observation. However, in order for the recoiling AGN to produce narrow-line emission, it must be observed shortly after the kick while it still inhabits a dense gaseous region, implying a large total kick velocity (vk approximately greater than 2000 km s(exp -1)). For the dual AGN scenario, an unusually large broad-line offset is required, and the best match to the observed morphology requires a galaxy that is less luminous than CID-42. Further, the lack of X-ray emission from one of the two optical nuclei is not easily attributed to an intrinsically quiescent SMBH or to a Compton-thick galactic environment. While the current data do not allow either the recoiling or the dual AGN scenario for CID-42 to be excluded, our models highlight the most relevant parameters for distinguishing these possibilities with future observations. In particular, high-quality, spatially-resolved spectra that can pinpoint the origin of the broad and narrow line features will be critical for determining the nature of this unique source.

Blecha, Laura↗

Delayed Development of Cool Plasmas in X-Ray Flares from the Young Sun-like Star κ1 Ceti

The Neutron star Interior Composition Explorer (NICER) X-ray observatory observed two powerful X-ray flares equivalent to superflares from the nearby young solar-like star κ 1 Ceti in 2019. NICER follows each flare from the onset through the early decay, collecting over 30 counts s −1 near the peak, enabling a detailed spectral variation study of the flare rise. The flare in September varies quickly in ∼800 s, while the flare in December has a few times longer timescale. In both flares, the hard-band (2–4 keV) light curves show typical stellar X-ray flare variations with a rapid rise and slow decay, while the soft X-ray light curves, especially of the September flare, have prolonged flat peaks. The time-resolved spectra require two temperature plasma components at kT ∼0.3–1 and ∼2–4 keV. Both components vary similarly, but the cool component lags by ∼200 s with a four to six times smaller emission measure (EM) compared to the hot component. A comparison with hydrodynamic flare loop simulations indicates that the cool component originates from X-ray plasma near the magnetic loop footpoints that mainly cools via thermal conduction. The time lag represents the travel time of the evaporated gas through the entire flare loop. The cool component has a several times smaller EM than its simulated counterpart, suggesting a suppression of conductive cooling, possibly by the expansion of the loop cross-sectional area or turbulent fluctuations. The cool component's time lag and EM ratio provide important constraints on the flare loop geometry.

main sequence stars↗