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.

SourceResearch connectionFrom yearThrough year

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

At least 595 records · Page 33

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Measurement and Modeling of Electron Temperature in Laboratory Photoionized Plasmas Relevant to Astrophysics

The electron temperature of photoionized plasmas characterizes the thermalization of photoelectrons, impacts the charge-state distribution, emissivity, and opacity through atomic recombination processes, and is needed to perform detailed comparisons with theory predictions. We discuss temperature measurements in laboratory photoionized plasmas and a comparison with model calculations done with several theory approximations and codes. These include a radiation-hydrodynamics simulation and two nonequilibrium heating and ionization models that tracked the evolution of the internal energy of the electrons. Furthermore, for the same physics model and X-ray flux time history, calculations were performed assuming steady-state or time-dependent conditions. The time history of steady-state results correlates with that of the X-ray flux, while that of the time-dependent cases does not, and it is qualitatively and quantitatively different from the steady-state case. Steady-state results significantly overestimated temperature measurements, while time-dependent results produced better approximations, which suggests the importance of transient effects in the experiment and also the need for time-resolved measurements.

Laboratory astrophysics↗

DOE OSTI2026

Dwarf Galaxy Number Counts within 25 Mpc: Predictions from Local Group Analogs in TNG50

The modern generation of wide-field galaxy surveys, such as LSST, Euclid, and Roman, will enable studies of dwarf galaxies (10 6 ≤ M * /M ⊙ ≤ 10 9 ) beyond the Local Group (LG) in unprecedented detail. Improved theoretical understanding of this population is necessary to guide these observations, since predictions in this regime are generally limited to specific environments like the LG. We present predictions for the population of dwarf galaxies from the TNG50 run of the IllustrisTNG suite of cosmological hydrodynamical simulations, focusing on the environments within 1 < D/Mpc < 25 of LG analogs at z = 0. In the simulated sample, there are ∼1000 and ∼12,000 dwarf galaxies within 10 and 25 Mpc, respectively. We compare our results with the 50 Mpc Galaxy Catalog and estimate that current observations are highly incomplete at low masses: for 10 6 ≤ M * /M ⊙ ≤ 10 7 (−13 ≲ M r ≲ −10), we find completeness fractions of ∼23% within 10 Mpc and ∼4% within 25 Mpc. The simulated galaxies below the completeness limits of the observations exist in a range of environments, with notable populations of field dwarfs at all distances and satellites around centrals with masses 10 8 ≲ M * /M ⊙ ≲ 10 11 within 10–25 Mpc. We find that there are ∼8 times more quiescent dwarf galaxies in the TNG50 sample than are currently cataloged. Our results suggest that upcoming observations should uncover a substantial population of dwarf galaxies, and that ≳15% of these will be red, currently quenched galaxies in the field.

Shread, Evangela E. [California Institute of Techn↗

DOE OSTI2026

ODIN: Characterizing the Three-dimensional Structure of Two Protocluster Complexes at z = 3.1

We present a detailed study of the 3D morphology of two extended associations of multiple protoclusters at z = 3.1. These protocluster "complexes," designated COSMOS-z3.1-A and COSMOS-z3.1-C, are the most prominent overdensities of z = 3.1 Lyα emitters (LAEs) identified in the COSMOS field by the One-hundred-deg$^{2}$ DECam Imaging in Narrowbands survey. These protocluster complexes have been followed up with extensive spectroscopy from Keck, Gemini, and DESI. Using a probabilistic method that combines photometrically selected and spectroscopically confirmed LAEs, we reconstruct the 3D structure of these complexes on scales of ≈50 cMpc. We validate our reconstruction method using the IllustrisTNG300-1 cosmological hydrodynamical simulation and show that it consistently outperforms approaches relying solely on spectroscopic data. The resulting 3D maps reveal that both complexes are irregular and elongated along a single axis, emphasizing the impact of sightline on our perception of structure morphology. The complexes consist of multiple density peaks, 10 in COSMOS-z3.1-A and 4 in COSMOS-z3.1-C. The former is confirmed to be a proto-supercluster, similar to Hyperion at z = 2.4 but observed at an even earlier epoch. Multiple "tails" connected to the cores of the density peaks are seen, likely representing cosmic filaments feeding into these extremely overdense regions. The 3D reconstructions further provide strong evidence that Lyα blobs preferentially reside in the outskirts of the highest density regions. Descendant mass estimates of the density peaks suggest that COSMOS-z3.1-A and COSMOS-z3.1-C will evolve to become ultramassive structures by z = 0, with total masses log ( M / M ,⊙ ,) ≳ 15.3 , exceeding that of Coma.

Ramakrishnan, Vandana [Purdue U., West Lafayette] ↗

DOE OSTI2020

Ring Morphology with Dust Coagulation in Protoplanetary Disks

Tidal interactions between the embedded planets and their surrounding protoplanetary disks are often postulated to produce the observed complex dust substructures, including rings, gaps, and asymmetries. In this Letter, we explore the consequences of dust coagulation on the dust dynamics and ring morphology. Coagulation of dust grains leads to dust size growth that, under typical disk conditions, produces faster radial drifts, potentially threatening the dust ring formation. Utilizing 2D hydrodynamical simulations of protoplanetary disks that include a full treatment of dust coagulation, we find that if the planet does not open a gap quickly enough, the formation of an inner ring is impeded due to dust coagulation and subsequent radial drift. Furthermore, we find that a "buildup" of submillimeter-sized grains often appears in the dust emission at the outer edge of the dust disk.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Planet-induced Vortices with Dust Coagulation in Protoplanetary Disks

In this work, we study how the dust coagulation/fragmentation will influence the evolution and observational appearances of vortices induced by a massive planet embedded in a low-viscosity disk by performing global 2D high-resolution hydrodynamical simulations. Within the vortex, due to its higher gas surface density and steeper pressure gradients, dust coagulation, fragmentation, and drift (to the vortex center) are all quite efficient, producing dust particles ranging from 1 μm to ~1.0 cm, as well as an overall high dust-to-gas ratio (above unity). In addition, the dust size distribution is quite nonuniform inside the vortex, with the mass-weighted average dust size at the vortex center (~4.0 mm) being a factor of ~10 larger than other vortex regions. Both large (~millimeter) and small (tens of microns) particles contribute strongly to affect the gas motion within the vortex. As such, we find that the inclusion of dust coagulation has a significant impact on the vortex lifetime and the typical vortex lifetime is about 1000 orbits. Furthermore, after the initial gaseous vortex is destroyed, the dust spreads into a ring with a few remaining smaller gaseous vortices with a high dust concentration and a large maximum size (~millimeter). At late time, the synthetic dust continuum images for the coagulation case show as a ring inlaid with several hot spots at the 1.33 mm band, while only distinct hot spots remain at 7.0 mm.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Emergence of Microphysical Bulk Viscosity in Binary Neutron Star Postmerger Dynamics

Abstract In nuclear matter in isolated neutron stars, the flavor content (e.g., proton fraction) is subject to weak interactions, establishing flavor ( β -)equilibrium. However, there can be deviations from this equilibrium during the merger of two neutron stars. We study the resulting out-of-equilibrium dynamics during the collision by incorporating direct and modified Urca processes (in the neutrino-transparent regime) into general-relativistic hydrodynamics simulations with a simplified neutrino transport scheme. We demonstrate how weak-interaction-driven bulk viscosity in postmerger simulations can emerge and assess the bulk viscous dynamics of the resulting flow. We further place limits on the impact of the postmerger gravitational-wave strain. Our results show that weak-interaction-driven bulk viscosity can potentially lead to a phase shift of the postmerger gravitational-wave spectrum, although the effect is currently on the same level as the numerical errors of our simulation.

Astronomy & Astrophysics↗

DOE OSTI2025

Measuring the Mean Free Path of H i Ionizing Photons at 3.2 ≤ z ≤ 4.6 with DESI Y1 Quasars

The mean free path of ionizing photons in the intergalactic medium (IGM; λ $^{912}_{mpf}$ ) is a crucial quantity in modeling the ionization state of the IGM and the extragalactic ultraviolet background, and it is widely used in hydrodynamical simulations of galaxies and reionization. We construct the largest quasar spectrum data set to date—12,595 signal-to-noise ratio > 3 spectra—using the Y1 observations of the Dark Energy Spectroscopic Instrument, to make the most precise model-independent measurement of the mean free path at 3.2 ≤ z ≤ 4.6. By stacking the spectra in 17 redshift bins and modeling the Lyman-continuum profile, we get a redshift evolution λ $^{912}_{mpf}$ ∝ ( 1 + z ) − 4.27 at 2 ≤ z ≤ 5, which is much shallower than the previous estimate, λ $^{912}_{mpf}$ ∝ ( 1 + z ) − 5.4 . We then explore the sources of systematic bias, including the choice of the intrinsic quasar continuum, the consideration of Lyman-series opacity and Lyman-limit opacity evolution, and the definition of λ $^{912}_{mpf}$ . Combining our results with estimates of λ $^{912}_{mpf}$ at higher redshifts, we conclude at high confidence that the evolution in λ $^{912}_{mpf}$ steepens at z ≈ 5. We interpret this inflection as the transition from the end of H i reionization to a fully ionized plasma that characterizes the IGM of the past ∼10 billion yr.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Constraints on Nonthermal Pressure at Galaxy Cluster Outskirts from a Joint SPT and XMM-Newton Analysis

We present joint South Pole Telescope and XMM-Newton observations of eight massive galaxy clusters (0.8–2 × 10$^{15}$ M$_{⊙}$) spanning a redshift range of 0.16–0.35. Employing a novel Sunyaev–Zel’dovich + X-ray fitting technique, we effectively constrain the thermodynamic properties of these clusters out to the virial radius. The resulting best-fit electron density, deprojected temperature, and deprojected pressure profiles are in good agreement with previous observations of massive clusters. For the majority of the cluster sample (five out of eight clusters), the entropy profiles exhibit a self-similar behavior near the virial radius. We further derive hydrostatic mass, gas mass, and gas fraction profiles for all clusters up to the virial radius. Comparing the enclosed gas fraction profiles with the universal gas fraction profile, we obtain nonthermal pressure fraction profiles for our cluster sample at >0.5R$_{500}$, demonstrating a steeper increase between R$_{500}$ and R$_{200}$ that is consistent with the hydrodynamical simulations. Our analysis yields nonthermal pressure fraction ranges of 8%–28% (median: 15% ± 11%) at R$_{500}$ and 21%–35% (median: 27% ± 12%) at R$_{200}$. Notably, weak-lensing mass measurements are available for only four clusters in our sample, and our recovered total cluster masses, after accounting for nonthermal pressure, are consistent with these measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Dark matter decay signals in cosmic filaments

Cosmic filaments form the backbone of the cosmic web, yet their properties and evolution remain uncertain. Using the EAGLE hydrodynamical simulations, we investigate the dark matter density profiles in filaments and their implications for dark matter decay signals. We show that GeV-scale dark matter particles decaying into electron-positron pairs can produce detectable radio synchrotron emission. By leveraging stacked filament radio data, we place stringent constraints on the dark matter decay lifetime, improving existing limits by up to two orders of magnitude for strong filamentary magnetic fields.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Constraints on cosmology and baryonic feedback with joint analysis of Dark Energy Survey Year 3 lensing data and ACT DR6 thermal Sunyaev-Zel'dovich effect observations

We present a joint analysis of weak gravitational lensing (shear) data obtained from the first three years of observations by the Dark Energy Survey and thermal Sunyaev-Zel'dovich (tSZ) effect measurements from a combination of Atacama Cosmology Telescope (ACT) and Planck data. A combined analysis of shear (which traces the projected mass) with the tSZ effect (which traces the projected gas pressure) can jointly probe both the distribution of matter and the thermodynamic state of the gas, accounting for the correlated effects of baryonic feedback on both observables. We detect the shear$~\times~$tSZ cross-correlation at a 21$\sigma$ significance, the highest to date, after minimizing the bias from cosmic infrared background leakage in the tSZ map. By jointly modeling the small-scale shear auto-correlation and the shear$~\times~$tSZ cross-correlation, we obtain $S_8 = 0.811^{+0.015}_{-0.012}$ and $\Omega_{\rm m} = 0.263^{+0.023}_{-0.030}$, results consistent with primary CMB analyses from Planck and P-ACT. We find evidence for reduced thermal gas pressure in dark matter halos with masses $M < 10^{14} \, M_{\odot}/h$, supporting predictions of enhanced feedback from active galactic nuclei on gas thermodynamics. A comparison of the inferred matter power suppression reveals a $2-4\sigma$ tension with hydrodynamical simulations that implement mild baryonic feedback, as our constraints prefer a stronger suppression. Finally, we investigate biases from cosmic infrared background leakage in the tSZ-shear cross-correlation measurements, employing mitigation techniques to ensure a robust inference. Our code is publicly available on GitHub.

Pandey, S. [Johns Hopkins U.; Columbia U.] (ORCID:↗

DOE OSTI2023

Proton Number Cumulants and Correlation Functions from Hydrodynamics and the QCD Phase Diagram

We analyze the behavior of (net-)proton number cumulants in central collisions of heavy ions across a broad collision energy range by utilizing hydrodynamic simulations. The calculations incorporate essential non-critical contributions to proton fluctuations such as repulsive baryonic core and exact baryon-number conservation. The experimental data are consistent with non-critical physics at collision energies of √SNN ≥ 20 GeV. The data from the STAR and HADES collaborations at lower collision energies indicate an excess of (multi-)proton correlations over the non-critical reference. This observation is discussed in the context of different mechanisms, including the possibility of a critical point in the baryon-rich region of the QCD phase diagram.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Design of inertial fusion implosions reaching the burning plasma regime

One of the last remaining milestones in fusion research before reaching ignition is creating a burning plasma state, where alpha particles from deuterium-tritium (DT) fusion reactions redeposit their energy as the dominant source of heating in the plasma. The indirect-drive inertial confinement fusion approach at the National Ignition Facility (NIF) uses a laser-generated radiation cavity (hohlraum) to spherically implode DT fuel to high temperatures and densities in a central ”hot spot”. Here, we deliver more energy to the hot spot than ever before, while maintaining the extreme pressures required for inertial confinement, by increasing the size of the implosion compared to previous experiments. We develop more efficient hohlraums, to drive these larger implosions within NIF’s current laser energy and power capability and control symmetry by moving energy between laser beams and by changing the shape of the hohlraum. These designs resulted in record fusion powers of 1.5 petawatts, greater than the input power of the laser, and 170 kJ of fusion energy. Radiation hydrodynamics simulations show alpha particle heating as the dominant term in the hot spot energy balance, e.g. a burning plasma state. This work is expected to motivate future studies of burning plasmas and improve predictive capability by providing a benchmark for modeling used to understand the proximity to ignition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Experimental observations of detached bow shock formation in the interaction of a laser-produced plasma with a magnetized obstacle

The magnetic field produced by planets with active dynamos, like the Earth, can exert sufficient pressure to oppose supersonic stellar wind plasmas, leading to the formation of a standing bow shock upstream of the magnetopause, or pressure-balance surface. Scaled laboratory experiments studying the interaction of an inflowing solar wind analog with a strong, external magnetic field are a promising new way to study magnetospheric physics and to complement existing models, although reaching regimes favorable for magnetized shock formation is experimentally challenging. This paper presents experimental evidence of the formation of a magnetized bow shock in the interaction of a supersonic, super-Alfvenic plasma with a strongly magnetized obstacle at the OMEGA laser facility. The solar wind analog is generated by the collision and subsequent expansion of two counter- propagating, laser-driven plasma plumes. The magnetized obstacle is a thin wire, driven with strong electrical currents. Hydrodynamic simulations using the FLASH code predict that the colliding plasma source meets the criteria for bow shock formation. Spatially resolved, optical Thomson scat- tering measures the electron number density, and optical emission lines provide a measurement of the plasma temperature, from which we infer the presence of a fast magnetosonic shock far upstream of the obstacle. Proton images provide a measure of large-scale features in the magnetic field topology, and reconstructed path-integrated magnetic field maps from these images suggest the formation of a bow shock upstream of the wire and as a transient magnetopause. We compare features in the reconstructed fields to two-dimensional MHD simulations of the system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1975

Applications of Computer Graphics in Engineering

Various applications of interactive computer graphics to the following areas of science and engineering were described: design and analysis of structures, configuration geometry, animation, flutter analysis, design and manufacturing, aircraft design and integration, wind tunnel data analysis, architecture and construction, flight simulation, hydrodynamics, curve and surface fitting, gas turbine engine design, analysis, and manufacturing, packaging of printed circuit boards, spacecraft design.

Source record↗

NASA NTRS1984

Rapid expansion of polytropes

An extremely massive unevolved star may be idealized as an n = 3 polytropic sphere supported by radiation pressure. Such a polytrope is subject to explosive, supersonic expansion, whose later stages can be described by similarity solution of the dynamical equations. This picture is confirmed by numerical hydrodynamic simulation, which shows that an n = 3 polytrope, initially at rest near equilibrium, undergoes an explosive expansion which approaches the similarity solution at large times.

Barnes, A.↗

NASA NTRS1985

Jet production in super-Eddington accretion disks

A two-dimensional, radiation-coupled, Newtonian hydrodynamic simulation is reported for a super-Eddington, mass accretion rate, M = 4 M(E) disk accretion flow onto a 3-solar mass pseudoblack hole. Near the disk midplane, convection cells effectively block the accretion flow, even though viscous heating maximizes there. Accretion predominantly occurs in a supersonic inflow which follows streamlines of approximately constant angular momentum. The optically thick inflow traps radiation so that 80 percent of the luminosity is absorbed by the black hole; the emergent power is sub-Eddington. An axial jet self consistently forms just outside a conical photosphere which bounds the accretion zone; radiation pressure accelerates the jet to about 10 to the 10th cm/s. The jet's mass efflux is only 0.4 percent of the total mass accretion rate.

Eggum, G. E.↗

NASA NTRS1986

Jeans criterion in a turbulent medium

According to the classical Jeans analysis, all the molecular clouds of mass larger than a few 100 M(solar), size larger than about 1pc and kinetic temperature Tk less than 30K are gravitationally unstable. We have shown that in clouds supported by internal supersonic motions, local gravitational instabilities may appear within molecular clouds which are globally stable. The argument is threefold: (1) when the turbulent kinetic energy is included into the internal energy term, the virial equilibrium condition shows that molecular clouds such as those observed, which are gravitationally unstable according to the Jeans criterion, are indeed globally stable if supported by a turbulent velocity field of power spectrum steeper than 3; (2) 2D compressible hydrodynamical simulations show that a supersonic turbulent velocity field generates a turbulent pressure within clouds, the gradients of which stabilize the unstable scales (i.e., the largest scales and the cloud itself) against gravitational collapse; (3) an analysis similar to the Jeans approach but including the turbulent pressure gradient term, gives basically the same results as those given in (1). Clouds of mean density lower than a critical value are found to be stable even though more massive than their Jeans mass. In clouds of mean density larger than that critical value, the gravitational instability appears only over a range of scales smaller than the cloud size, the largest scales being stable. In practice, the observed mean densities are lower than this critical value: the observation of a small number of cores and stars of a few solar masses embedded in clouds of several hundred solar masses can only be understood in terms of small scale density fluctuations of large amplitude generated by the supersonic turbulence which would occasionally overtake the limit of gravitational stability.

Bonazzola, S.↗

NASA NTRS1986

Plasma turbulence and impulsive UV line emission in solar flares

Observations show that hard X-ray burst and UV lines rise and fall simultaneously on time scales of seconds. Hydrodynamic simulations of beam-heated atmospheres, based on collisional transport, however, produce only a gradual fall in UV emission, when the beam flux falls, due to the long time scale of conductive relaxation. It is suggested that this discrepancy might be explained by onset of plasma turbulence driven by the strong heat flux or by the beam return current going unstable. Such turbulence greatly reduces electrical and thermal conductivities. Fall in electrical conductivity reduces the hard X-ray flux by enhanced ohmic dissipation of the return current, while fall in thermal conductivity may cause the UV line to fall by reducing the transition region thickness.

Brown, John C.↗