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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Orbital Evolution of Binary Black Holes in Active Galactic Nucleus Disks: A Disk Channel for Binary Black Hole Mergers?

We perform a series of high-resolution 2D hydrodynamical simulations of equal-mass binary black holes (BBHs) embedded in active galactic nucleus (AGN) accretion disks to study whether these binaries can be driven to merger by the surrounding gas. We find that the gravitational softening adopted for the BBH has a profound impact on this result. When the softening is less than 10% of the binary separation, we show that, in agreement with recent simulations of isolated equal-mass binaries, prograde BBHs expand in time rather than contract. Eventually, however, the binary separation becomes large enough that the tidal force of the central AGN disrupts them. Only when the softening is relatively large do we find that prograde BBHs harden. We determine through detailed analysis of the binary torque, that this dichotomy is due to a loss of spiral structure in the circum-single disks orbiting each black hole when the softening is a significant fraction of the binary separation. Properly resolving these spirals—both with high resolution and small softening—results in a significant source of binary angular momentum. Only for retrograde BBHs do we find consistent hardening, regardless of softening, as these BBHs lack the important spiral structure in their circum-single disks. This suggests that the gas-driven inspiral of retrograde binaries can produce a population of compact BBHs in the gravitational-wave-emitting regime in AGN disks, which may contribute a large fraction to the observed BBH mergers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Studies of localized mix and asymmetry caused by the stalk in inertial fusion experiments at OMEGA

In an inertially confined fusion experiment, a capsule of fusion fuel is compressed to high temperatures and densities using lasers, creating a hot plasma in which thermonuclear fusion reactions take place. To achieve efficient conversion of laser energy to fuel heating, the implosion must be as spherically symmetric as possible. One common source of asymmetry that degrades performance is the stalk and glue spot that hold the capsule in place before the experiment. The exact nature and magnitude of this perturbation had not been previously quantified. In this work, inertially confined fusion experiments were executed with an additional stalk on the capsule to measure and quantify the effect of the stalk and glue spot on the symmetry and performance. It is found that the stalk drives low-mode asymmetry by shadowing the laser drive on the part of the shell near the glue spot, enhancing the hot-spot velocity toward the stalk by 40–80 km/s and elongating the hot-spot along the stalk axis. It is also found that localized mix due to the additional stalk has a negligible effect on implosion performance. Two-dimensional radiation-hydrodynamic simulations using the xRAGE code reproduce the direction of the experimentally observed change to hot-spot velocity and implosion performance metrics, but not the observed elongation, suggesting that this feature results from an interaction between stalk perturbations and other preexisting asymmetry seeds, such as laser drive asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simultaneous inference of equation of state parameters and unknown data errors with uncertainty quantification via hierarchical Bayesian posterior maximization

Equations of state (EOSs) are a key component in running hydrodynamic simulations as they relate the thermodynamic states for the material. The Davis reactants EOS is commonly used for modeling high explosives (HEs), and the EOS model parameters are calibrated using material specific data. The calibrations are often performed with uncertainty quantification via Bayesian inference to account for uncertainty in the data and generate ensembles of likely parameters. However, there are relatively few HE data sets to use for calibration and many are historical and lack error information. In this work, we simultaneously calibrate the Davis reactants EOS model parameters and unknown data error terms for the high explosive PBX 9501. To quantify the uncertainty in the models and the data, we use a Bayesian framework for the calibration and compute the hierarchical Bayesian posterior distribution with both a posteriori maximization approach and Markov Chain Monte Carlo. In general, we find that, given our assumptions, the two approaches result in similar calibrated parameters, posterior covariance matrices, and insights about the parameters but that the posterior maximization requires far less computational resources.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

How the Galaxy–Halo Connection Depends on Large-scale Environment

We investigate the connection between galaxies, dark matter halos, and their large-scale environments at z = 0 with Illustris TNG300 hydrodynamic simulation data. We predict stellar masses from subhalo properties to test two types of machine learning (ML) models: explainable boosting machines (EBMs) with simple galaxy environment features and E(3)-invariant graph neural networks (GNNs). The best-performing EBM models leverage spherically averaged overdensity features on 3 Mpc scales. Interpretations via SHapley Additive exPlanations also suggest that in the context of the TNG300 galaxy–halo connection, simple spherical overdensity on ∼3 Mpc scales is more important than cosmic web distance features measured using the DisPerSE algorithm. Meanwhile, a GNN with connectivity defined by a fixed linking length, L, outperforms the EBM models by a significant margin. As we increase the linking length scale, GNNs learn important environmental contributions up to the largest scales we probe (L = 10 Mpc). We conclude that 3 Mpc distance scales are most critical for describing the TNG galaxy–halo connection using the spherical overdensity parameterization, but that information on larger scales, which is not captured by simple environmental parameters or cosmic web features, can further augment these models. Our study highlights the benefits of using interpretable ML algorithms to explain models of astrophysical phenomena, and the power of using GNNs to flexibly learn complex relationships directly from data while imposing constraints from physical symmetries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

X-ray-line coincidence photopumping in a potassium-chlorine mixed plasma

Exploiting the multiple long pulse capability and suite of x-ray diagnostics of the Orion laser, we have set out to explore line coincidence photopuming—the enhancement in population of an atomic level brought on by resonant absorption of x rays from a different emitting ion. Unlike previous work, the two ions are in the same plasma and so the experiment is an x-ray analog of the well-known Bowen resonance fluorescence mechanism that operates in astrophysical situations in the optical region. Herein, our measurements have shown enhanced fluorescence in a chlorine plasma, attributable to line coincident photopumping from co-mixed potassium ions. To detect this relatively low signal-to-noise phenomenon, the data from multiple shots are combined, and the statistical method of bootstrapping is used to assign a confidence value to the measured enhancement, resulting in an estimate of the enhancement of 39 ± 18 16 % compared to the null case, where no pumping occurs. The experimental results have been compared to coupled radiation-transport and radiation hydrodynamics simulations using the CRETIN code together with the NYM radiation hydrodynamics model and agreement has been found, with the simulations also predicting modest enhancement.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Betelgeuse as a Merger of a Massive Star with a Companion

We investigate the merger between a 16M ⊙ star, on its way to becoming a red supergiant (RSG), and a 4M ⊙ main-sequence companion. Our study employs three-dimensional hydrodynamic simulations using the state-of-the-art adaptive mesh refinement code Octo-Tiger. The initially corotating binary undergoes interaction and mass transfer, resulting in the accumulation of mass around the companion and its subsequent loss through the second Lagrangian point (L2). The companion eventually plunges into the envelope of the primary, leading to its spin-up and subsequent merger with the helium core. We examine the internal structural properties of the post-merger star, as well as the merger environment and the outflow driven by the merger. Our findings reveal the ejection of approximately ∼0.6 M ⊙ of material in an asymmetric and somewhat bipolar outflow. We import the post-merger stellar structure into the MESA stellar evolution code to model its long-term nuclear evolution. In certain cases, the post-merger star exhibits persistent rapid equatorial surface rotation as it evolves in the H–R diagram toward the observed location of Betelgeuse. These cases demonstrate surface rotation velocities of a similar magnitude to those observed in Betelgeuse, along with a chemical composition resembling that of Betelgeuse. In other cases, efficient rotationally induced mixing leads to slower surface rotation. This pioneering study aims to model stellar mergers across critical timescales, encompassing dynamical, thermal, and nuclear evolutionary stages.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

An experimentally informed design process for future inertial confinement fusion facilities

The achievement of ignition in the laboratory has renewed interest in defining the requirements for a future high-gain inertial confinement fusion (ICF) facility. Our best chance of predicting future ICF performance is with 3-D radiation hydrodynamic simulations that have been benchmarked against experimental data, but their high computational cost is prohibitive for use in practical design studies. We introduce a hierarchical approach where 3-D simulations are tuned to match experimental measurements and used to train 3-D degradation models in 1-D simulations allowing for accurate predictions over the entire OMEGA direct-drive database. A genetic algorithm was used in combination with the trained 1-D simulations to search for optimal direct-drive implosion designs at driver energies ranging from 20 kJ to 10 MJ. As the fidelity of 3-D codes improves, this approach will provide a viable experimentally informed tool for defining the next ICF facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Hydrodynamic expansion and near-infrared absorption of x-ray heated aluminum plasmas

We use x-ray pulses from dense argon plasmas at the Z Machine (Sandia National Laboratories) to generate hypersonic aluminum plasmas akin to material ejecta during proposed planetary defense missions, fusion reactor wall excursions, and other high-energy density processes. Near-infrared absorption is used to diagnose the controlled expansion of the plasmas through cylindrical cavities following their generation from x-ray heating of solid aluminum 7075 alloy. The data are compared to multidimensional radiation hydrodynamics simulations utilizing the ALEGRA multiphysics code, accounting for the dynamics of radiation scattering, material phase change, plasma expansion, thermal re-irradiation, and interactions with the cavity and with the infrared beams. To allow for accurate simulation, density functional theory is used to apply the Hagen–Rubens relation for the far-infrared and is adjoined with a detailed configuration accounting model using the Propaceos code, producing opacities spanning 10 −1 –10 4 eV photon energy for aluminum 7075 alloy, and in comparison with pure aluminum. The model is found to agree with experimental data in the higher-fluence regime when the Hagen–Rubens relation is applied. The ejected material, which is observed to travel up to 55 km/s, is comprised of a strongly ionized, non-LTE plasma front at ∼10 eV temperature followed by a weakly ionized LTE gas at higher density. The present findings lend some confidence to the broad-range equation of state and infrared opacity models for weakly ionized aluminum plasmas while demonstrating an approach to their future refinement, with potential application to astrophysical plasmas and other extreme processes.

Adiabatic process↗
DOE OSTI2025

Evaluating nonlocal heat transport in directly driven chromium spheres using x-ray spectroscopy

We report on experiments investigating heat transport in laser-generated plasmas using directly driven chromium spheres. The spheres are fielded at the OMEGA laser facility and are driven with laser intensities of 5×10 14 Wcm −2 . Plasma conditions in the corona and scattered light are measured experimentally and compared against predictions from two-dimensional (2D) radiation-hydrodynamic simulations using different heat transport models. Spectroscopic analysis of x-ray self-emission is used as an additional diagnostic. X-ray emission is integrated over a large region of the plasma, probing regions that are not observed by localized optical Thomson scattering. In particular, x-ray emission peaks near the plasma critical density, so emission from optically thin lines provides information on plasma conditions where nonlocal transport is most likely to be significant. Three common heat transport models are considered: local transport with flux limiters f = 0.15 and f = 0.03, and the nonlocal Schurtz–Nicolai–Busquet (SNB) model. Consistent with previous work, both the high-flux (f = 0.15) and SNB models show good agreement with experimentally measured plasma conditions in the corona despite overpredicting laser absorption, whereas the low-flux (f = 0.03) model fails to match any experimental data. Conditions inferred from x-ray self-emission line ratios support this conclusion during the period of laser peak power, although synthetic spectra for all models fail to match the experiment during the transient portions of the pulse. For these reasons, the low-flux model is again rejected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2008

Modeling Flows Around Merging Black Hole Binaries

Coalescing massive black hole binaries are produced by the merger of galaxies. The final stages of the black hole coalescence produce strong gravitational radiation that can be detected by the space-borne LISA. In cases in which the black hole merger takes place in the presence of gas and magnetic fields, various types of electromagnetic signals may also be produced. Modeling such electromagnetic counterparts of the final merger requires evolving the behavior of both gas and fields in the strong-field regions around the black holes. We have taken a first step towards this problem by mapping the flow of pressureless matter in the dynamic, 3-D general relativistic spacetime around the merging black holes. We report on the results of these initial simulations and discuss their likely importance for future hydrodynamical simulations.

Centrella, Joan↗
NASA NTRS2009

Massive Black Hole Mergers: Can we see what LISA will hear?

Coalescing massive black hole binaries are formed when galaxies merge. The final stages of this coalescence produce strong gravitational wave signals that can be detected by the space-borne LISA. When the black holes merge in the presence of gas and magnetic fields, various types of electromagnetic signals may also be produced. Modeling such electromagnetic counterparts requires evolving the behavior of both gas and fields in the strong-field regions around the black holes. We have taken a first step towards this problem by mapping the flow of pressureless matter in the dynamic, 3-D general relativistic spacetime around the merging black holes. We report on the results of these initial simulations and discuss their likely importance for future hydrodynamical simulations.

Centrella, Joan↗
DOE OSTI2025

Investigation of an octupole breathing mode of 208 Pb as a resolution to the elliptical-to-triangular azimuthal anisotropy puzzle in ultracentral relativistic heavy ion collisions

Relativistic heavy ion collisions provide a unique opportunity to probe the nuclear structure by taking an instantaneous snapshot of the colliding nuclei and converting it into momentum anisotropies of final emitted hadrons. A long-standing puzzle of too large a ratio of the elliptical-to-triangular (𝑣 2 -to-𝑣 3 ) anisotropies in ultracentral 208 Pb + 208 Pb collisions at the Large Hadron Collider (LHC) cannot be solved simply by hydrodynamic simulations with initial conditions containing the spherical or certain deformed shape of 208 Pb. Here, in this Letter, using the iEBE-VISHNU relativistic viscous hydrodynamic hybrid model simulations with the TR⁢ENTo initial condition, we show that a dynamic octupole deformation—a shape-breathing mode of 208 Pb —could potentially solve the 𝑣2-to-𝑣3 puzzle and simultaneously describe the 𝑣 3 ⁡{4} data measured in experiment. Our results highlight the unique capability of capturing transient collective properties of nuclei on yoctosecond (10 −24 s) timescales, unfeasible with low-energy nuclear reactions.

Xu, Hao-jie [Huzhou University, Zhejiang (China)] ↗
DOE OSTI2020

High mass and halo resolution from fast low resolution simulations

Generating mocks for future sky surveys requires large volumes and high resolutions, which is computationally expensive even for fast simulations. Here we try to develop numerical schemes to calibrate various halo and matter statistics in fast low resolution simulations compared to high resolution N-body and hydrodynamic simulations. For the halos, we improve the initial condition resolution and develop a halo finder "relaxed-FoF", where we allow different linking lengths for different halo mass and velocity dispersions. We show that our relaxed-FoF halo finder improves the common statistics, such as halo bias, halo mass function, halo auto power spectrum, cross correlation coefficient with the reference halo catalog, and halo-matter cross power spectrum. We also calibrate small-scale velocities of small halos to improve the power spectrum in redshift space. For the matter statistics, we incorporate the potential gradient descent (PGD) method into fast simulations to improve the matter distribution at nonlinear scales. By building a lightcone output, we show that the PGD method significantly improves the weak lensing convergence tomographic power spectrum. With these improvements FastPM is comparable to the high resolution full N-body simulation of the same mass resolution, with two orders of magnitude fewer time steps. These techniques can be used to improve the halo and matter statistics of FastPM simulations for mock catalogs of future surveys such as DESI and LSST.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The growth of intracluster light in XCS-HSC galaxy clusters from 0.1 < z < 0.5

ABSTRACT We estimate the intracluster light (ICL) component within a sample of 18 clusters detected in the XMM Cluster Survey (XCS) data using the deep (∼26.8 mag) Hyper Suprime-Cam Subaru Strategic Programme data release 1 i-band data. We apply a rest-frame μB = 25 mag arcsec−2 isophotal threshold to our clusters, below which we define light as the ICL within an aperture of RX,500 (X-ray estimate of R500) centred on the brightest cluster galaxy (BCG). After applying careful masking and corrections for flux losses from background subtraction, we recover ∼20 per cent of the ICL flux, approximately four times our estimate of the typical background at the same isophotal level (${\sim}5{{\ \rm per\ cent}}$). We find that the ICL makes up about ${\sim}24{{\ \rm per\ cent}}$ of the total cluster stellar mass on average (∼41 per cent including the flux contained in the BCG within 50 kpc); this value is well matched with other observational studies and semi-analytic/numerical simulations, but is significantly smaller than results from recent hydrodynamical simulations (even when measured in an observationally consistent way). We find no evidence for any links between the amount of ICL flux with cluster mass, but find a growth rate of 2–4 for the ICL between 0.1 < z < 0.5. We conclude that the ICL is the dominant evolutionary component of stellar mass in clusters from z ∼ 1. Our work highlights the need for a consistent approach when measuring ICL alongside the need for deeper imaging, in order to unambiguously measure the ICL across as broad a redshift range as possible (e.g. 10-yr stacked imaging from the Vera C. Rubin Observatory).

Furnell, Kate E.↗
DOE OSTI2021

Simulating highly eccentric common envelope jet supernova impostors

We conduct three-dimensional hydrodynamical simulations of eccentric common envelope jet supernova (CEJSN) impostors, i.e. a neutron star that crosses through the envelope of a red supergiant star on a highly eccentric orbit and launches jets as it accretes mass from the envelope. Because of numerical limitations, we apply a simple prescription where we inject the assumed jets’ power into two opposite conical regions inside the envelope. We find the outflow morphology to be very complicated, clumpy, and non-spherical, having a large-scale symmetry only about the equatorial plane. The outflow morphology can substantially differ between simulations that differ by their jets’ power. We estimate by simple means the light curve to be very bumpy, to have a rise time of one to a few months, and to slowly decay in about a year to several years. These eccentric CEJSN impostors will be classified as ‘gap’ objects, i.e. having a luminosity between those of classical novae and typical supernovae (termed also ILOTs for intermediate luminosity optical transients). We strengthen a previous conclusion that CEJSN impostors might account for some peculiar ILOTs, in particular those that might repeat over time-scales of months to years.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Evaluating quenching in cosmological simulations of galaxy formation with spectral covariance in the optical window

ABSTRACT Cosmological hydrodynamical simulations provide valuable insights on galaxy evolution when coupled with observational data. Comparisons with real galaxies are typically performed via scaling relations of the observables. Here, we follow an alternative approach based on the spectral covariance in a model-independent way. We build upon previous work by Sharbaf et al. that studied the covariance of high-quality SDSS (Sloan Digital Sky Survey) continuum-subtracted spectra in a relatively narrow range of velocity dispersion ($\sigma \in [100,150]$ km s$^{-1}$). Here, the same analysis is applied to synthetic data from the eagle and IllustrisTNG100 simulations, to assess the ability of these runs to mimic real galaxies. The real and simulated spectra are consistent regarding spectral covariance, although with subtle differences that can inform the implementation of subgrid physics. Spectral fitting done a posteriori on stacks segregated with respect to latent space reveals that the first principal component (PC1) is predominantly influenced by the stellar age distribution, with an underlying age–metallicity degeneracy. Good agreement is found regarding star formation prescriptions but there is disagreement with active galactic nucleus (AGN) feedback, that also affects the subset of quiescent galaxies. We show a substantial difference in the implementation of the AGN subgrid prescriptions, regarding central black hole seeding, that could lead to the mismatch. Differences are manifest between these two simulations in the star formation histories stacked with respect to latent space. We emphasize that this methodology only relies on the spectral covariance to assess whether simulations provide a true representation of galaxy formation.

Sharbaf, Z. (ORCID:0009000450545946)↗
DOE OSTI2022

Improving Predictive Capability in REHEDS Simulations with Fast, Accurate, and Consistent Non-Equilibrium Material Properties

Predictive design of REHEDS experiments with radiation-hydrodynamic simulations requires knowledge of material properties (e.g. equations of state (EOS), transport coefficients, and radiation physics). Interpreting experimental results requires accurate models of diagnostic observables (e.g. detailed emission, absorption, and scattering spectra). In conditions of Local Thermodynamic Equilibrium (LTE), these material properties and observables can be pre-computed with relatively high accuracy and subsequently tabulated on simple temperature-density grids for fast look-up by simulations. When radiation and electron temperatures fall out of equilibrium, however, non-LTE effects can profoundly change material properties and diagnostic signatures. Accurately and efficiently incorporating these non-LTE effects has been a longstanding challenge for simulations. At present, most simulations include non-LTE effects by invoking highly simplified inline models. These inline non-LTE models are both much slower than table look-up and significantly less accurate than the detailed models used to populate LTE tables and diagnose experimental data through post-processing or inversion. Because inline non-LTE models are slow, designers avoid them whenever possible, which leads to known inaccuracies from using tabular LTE. Because inline models are simple, they are inconsistent with tabular data from detailed models, leading to ill-known inaccuracies, and they cannot generate detailed synthetic diagnostics suitable for direct comparisons with experimental data. This project addresses the challenge of generating and utilizing efficient, accurate, and consistent non-equilibrium material data along three complementary but relatively independent research lines. First, we have developed a relatively fast and accurate non-LTE average-atom model based on density functional theory (DFT) that provides a complete set of EOS, transport, and radiative data, and have rigorously tested it against more sophisticated first-principles multi-atom DFT models, including time-dependent DFT. Next, we have developed a tabular scheme and interpolation methods that compactly capture non-LTE effects for use in simulations and have implemented these tables in the GORGON magneto-hydrodynamic (MHD) code. Finally, we have developed post-processing tools that use detailed tabulated non-LTE data to directly predict experimental observables from simulation output.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Numerical simulations of bent, disrupted radio jets

We present preliminary results from three-dimensional hydrodynamical simulations designed to investigate the physics of jet bending and disruption. The specific scenario considered here involves a mildly supersonic jet crossing a contact discontinuity at the interface between the interstellar medium (ISM) and the intercluster medium (ICM) and then encountering a cross-wind in the ICM. The resultant morphologies show many of the features observed in radio sources including jet flaring, bending, and extended tails.

Loken, Chris↗