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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Mitigation of deceleration-phase Rayleigh–Taylor instability growth in inertial confinement fusion implosions

Rayleigh–Taylor growth during shell deceleration is one of the main limiting factors for target performance in inertial confinement fusion implosions. Here, using analytical scaling laws and hydrodynamic simulations, we show that such amplification can be mitigated by reducing the initial mass density in the central target region. The perturbation growth reduction is caused by a smaller hot-spot convergence ratio during deceleration, increased density scale length, and enhanced ablation stabilization. The required central density reduction can be achieved using the dynamic shell formation concept.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Convergence of small scale Ly α structure at high- z under different reionization scenarios

ABSTRACT The Ly α forest (LAF) at z > 5 probes the thermal and reionization history of the intergalactic medium (IGM) and the nature of dark matter, but its interpretation requires comparison to cosmological hydrodynamical simulations. At high-z, convergence of these simulations is more exacting since transmission is dominated by underdense voids that are challenging to resolve. With evidence mounting for a late end to reionization, small structures down to the sub-kpc level may survive to later times than conventionally thought due to the reduced time for pressure smoothing to impact the gas, further tightening simulation resolution requirements. We perform a suite of simulations using the Eulerian cosmological hydrodynamics code Nyx, spanning domain sizes of 1.25 − 10 h−1 Mpc and 5 − 80 h−1 kpc cells, and explore the interaction of these variables with the timing of reionization on the properties of the matter distribution and the simulated LAF at z = 5.5. In observable Ly α power, convergence within 10 per cent is achieved for k < 0.1 s km–1, but larger k shows deviation of up to 20 per cent. While a later reionization retains more small structure in the density field, because of the greater thermal broadening there is little difference in the convergence of LAF power between early (z = 9) and later (z = 6) reionizations. We conclude that at z ∼ 5.5, resolutions of 10 kpc are necessary for convergence of LAF power at k < 0.1 s km–1, while higher-k modes require higher resolution, and that the timing of reionization does not significantly impact convergence given realistic photoheating.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Lyman-α transmission properties of the intergalactic medium in the CoDaII simulation

ABSTRACT The decline in abundance of Lyman-α (Lyα) emitting galaxies at z ≳ 6 is a powerful and commonly used probe to constrain the progress of cosmic reionization. We use the CoDaII simulation, which is a radiation hydrodynamic simulation featuring a box of ∼94 comoving Mpc side length, to compute the Lyα transmission properties of the intergalactic medium (IGM) at z ∼ 5.8 to 7. Our results mainly confirm previous studies, i.e. we find a declining Lyα transmission with redshift and a large sightline-to-sightline variation. However, motivated by the recent discovery of blue Lyα peaks at high redshift, we also analyse the IGM transmission on the blue side, which shows a rapid decline at z ≳ 6 of the blue transmission. This low transmission can be attributed not only to the presence of neutral regions but also to the residual neutral hydrogen within ionized regions, for which a density even as low as $n_{\rm HI}\sim 10^{-9}\, \mathrm{cm}^{-3}$ (sometimes combined with kinematic effects) leads to a significantly reduced visibility. Still, we find that $\sim 1{{\ \rm per\ cent}}$ of sightlines towards M1600AB ∼ −21 galaxies at z ∼ 7 are transparent enough to allow a transmission of a blue Lyα peak. We discuss our results in the context of the interpretation of observations.

Gronke, Max (ORCID:000000032491060X)↗
DOE OSTI2022

SDSS-IV MaNGA: Cannibalism Caught in the Act—On the Frequency of Occurrence of Multiple Cores in Brightest Cluster Galaxies

Although it is generally accepted that massive galaxies form in a two-phased fashion, beginning with a rapid mass buildup through intense starburst activities followed by primarily dry mergers that mainly deposit stellar mass at outskirts, the late time stellar mass growth of brightest cluster galaxies (BCGs), the most massive galaxies in the universe, is still not well understood. Several independent measurements have indicated a slower mass growth rate than predictions from theoretical models. We attempt to resolve the discrepancy by measuring the frequency of BCGs with multiple cores, which serve as a proxy of the merger rates in the central region and facilitate a more direct comparison with theoretical predictions. Using 79 BCGs at z = 0.06–0.15 with integral field spectroscopic data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) project, we obtain a multiple-core fraction of 0.11 ± 0.04 at z ≈ 0.1 within an 18 kpc radius from the center, which is comparable to the value of 0.08 ± 0.04 derived from mock observations of 218 simulated BCGs from the cosmological hydrodynamical simulation IllustrisTNG. We find that most cores that appear close to the BCGs from imaging data turn out to be physically associated systems. Anchoring on the similarity in the multiple-core frequency between the MaNGA and IllustrisTNG, we discuss the mass growth rate of BCGs over the past 4.5 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The Odd Dark Matter Halos of Isolated Gas-rich Ultradiffuse Galaxies

Abstract We analyze circular velocity profiles of seven ultradiffuse galaxies (UDGs) that are isolated and gas-rich. Assuming that the dark matter halos of these UDGs have a Navarro–Frenk–White (NFW) density profile or a Read density profile (which allows for constant-density cores), the inferred halo concentrations are systematically lower than the cosmological median, even as low as −0.6 dex (about 5 σ away) in some cases. Alternatively, similar fits can be obtained with a density profile that scales roughly as 1/ r 2 for radii larger than a few kiloparsecs. Both solutions require the radius where the halo circular velocity peaks ( R max ) to be much larger than the median expectation. Surprisingly, we find an overabundance of such large- R max halos in the IllustrisTNG dark-matter-only simulations compared to the Gaussian expectation. These halos form late and have higher spins compared to median halos of similar masses. The inner densities of the most extreme among these late-forming halos are higher than their NFW counterparts, leading to a ∼1/ r 2 density profile. However, the two well-resolved UDGs in our sample strongly prefer lower dark matter densities in the center than the simulated ones. Comparing to IllustrisTNG hydrodynamical simulations, we also find a tension in getting both low enough circular velocities and high enough halo mass to accommodate the measurements. Our results indicate that the gas-rich UDGs present a significant challenge for galaxy formation models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The development of thermodynamically consistent and physics-informed equation-of-state model through machine learning

Ab initio molecular dynamics (AIMD) simulations have become an important tool used in the construction of equations of state (EOS) tables for warm dense matter. Due to computational costs, only a limited number of system state conditions can be simulated, and the remaining EOS surface must be interpolated for use in radiation-hydrodynamic simulations of experiments. In this work, we develop a thermodynamically consistent EOS model that utilizes a physics-informed machine learning approach to implicitly learn the underlying Helmholtz free-energy from AIMD generated energies and pressures. The model, referred to as PIML-EOS, was trained and tested on warm dense polystyrene producing a fit within a 1% relative error for both energy and pressure and is shown to satisfy both the Maxwell and Gibbs–Duhem relations. In addition, we provide a path toward obtaining thermodynamic quantities, such as the total entropy and chemical potential (containing both ionic and electronic contributions), which are not available from current AIMD simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Beryllium–tungsten graded density inner shells in double shell capsules for improved hydrodynamic stability

The outer surface of the high-Z inner shell in the double shell configuration of inertial confinement fusion experiments experiences Rayleigh–Taylor instability growth during the implosion process due to inverted density and pressure gradients between a highly compressed foam interstitial layer and the accelerating dense inner shell. Graded density layers have long been known to reduce instability growth rates. In this study, we employ high-fidelity radiation hydrodynamic simulations to demonstrate this improved stability when grading beryllium into tungsten. We first characterize the response to L-band preheat of these layers using a newly calibrated radiation drive. While graded layer capsules suffer reduced performance (here, measured as DD neutron yield from a CD foam fuel) in 1D simulations due to reduced kinetic energy coupling and reduced fuel compression, they suffer less of a performance drop when 2D instabilities are accounted for. With the improved stability of graded layers, we explore the performance of capsules with larger fuel radii and thinner shells as a preliminary study to find new designs in which graded layers produce the highest yields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Constrained model calibration of grain structure dependent spall dynamics in shock-loaded tantalum

Here, we perform a gas gun experiment by shock loading tantalum samples of varying grain structures to assess the suitability of a numerical model for simulating spall behavior. The observed differences in spall strength, as well spallation and re-compression history, are not captured in uncalibrated hydrodynamic simulations. An optimization is performed on the Johnson spall model to determine the best parameters that fit the observed trends. Linear stability analysis is employed to motivate bounds on those parameters. Herein, optimized simulations agree well with the experimental results, reproducing pullback depth and recompression timescales across the different samples tested. Further, the observed pullback time of the single crystal sample was found to imply, via the stability analysis, a percolation threshold in good agreement with the theoretical value for a body centered cubic lattice. Therefore, the combined linear stability and percolation analysis shows promise and may be applied to other materials with diverse microstructures. Collectively, the findings demonstrate that the model is suitable for reproducing spall-induced free surface behavior across various microstructures, but also points to caution in using model coefficients for uncalibrated microstructures.

58 GEOSCIENCES↗
DOE OSTI2025

Resolving discrepancies in bang-time predictions for indirect-drive ICF experiments on the NIF: Insights from the Build-A-Hohlraum campaign

This study investigated discrepancies between measured and simulated x-ray drive in Indirect-Drive Inertial Confinement Fusion (ID-ICF) hohlraums at the National Ignition Facility. Despite advances in radiation-hydrodynamic simulations, a consistent “drive deficit” remains. Experimentally measured ID-ICF capsule bang-times are systematically 400–700 ps later than simulations predict. The Build-A-Hohlraum (BAH) campaign explored potential causes for this discrepancy by systematically varying hohlraum features, including laser entrance hole (LEH) windows, capsules, and gas fills. Overall, the agreement between simulated and experimental x-ray drive was found to be largely unaffected by these changes. The data allow us to exclude some hypotheses put forward to potentially explain the discrepancy. Errors in the local thermodynamic equilibrium (LTE) atomic modeling, errors in the modeling of LEH closure, and errors due to a lack of plasma species mix physics in simulations are shown to be inconsistent with our measurements. Instead, the data support the hypothesis that errors in NLTE emission modeling are a significant contributor to the discrepancy. X-ray emission in the 2–4 keV range is found to be approximately 30% lower than in simulations. This is accompanied by higher than predicted electron temperatures in the gold bubble region, pointing to errors in non-LTE modeling. Introducing an opacity multiplier of 0.87 on energy groups above 1.8 keV improves agreement with experimental data, reducing the bang-time discrepancy from 300 to 100 ps. These results underscore the need for refined NLTE opacity models to enhance the predictive power of hohlraum simulations.

Band emission↗
DOE OSTI2024

Nonplanar effects in simulations of laser-driven ejecta microjet experiments

Recent experiments of laser-driven ejecta microjets performed at OMEGA 60 reveal tortuous jets whereby the jets appear to deviate from their initial trajectory as they travel in vacuum. To understand these data, we perform two dimensional numerical simulations, considering different target thicknesses, pressures, and models of the drive conditions. In particular, modeling the finite laser spot size appears essential in reproducing qualitatively the non-planar shock observed in the experiment. Simulations capture jet deflection by accounting for a slight misalignment of the laser pointing with respect to the groove axis along with spatial variation of the laser pulse intensity. The principal physical mechanism appears to be that lateral momentum is imparted by release waves arising from the non-planar drive. The induced off-axis velocity is small in comparison to the jet axial velocity but integrates into a pronounced deflection over the course of the experiment. The analysis of jet axial and lateral mass distributions is found to be reproduced reasonably by the simulations. Simulated radiographs are in qualitative agreement with the experiments, though their differences point to potential shortcomings in modeling strictly three-dimensional experiments using two-dimensional hydrodynamic simulations. The simple analysis is able to explain part of the observed discrepancy in simulated vs experimental jet masses.

42 ENGINEERING↗
DOE OSTI2022

Role of self-generated magnetic fields in the inertial fusion ignition threshold

Magnetic fields spontaneously grow at unstable interfaces around hot-spot asymmetries during inertial confinement fusion implosions. Although difficult to measure, theoretical considerations and numerical simulations predict field strengths exceeding 5 kT in current National Ignition Facility experiments. Magnetic confinement of electrons then reduces the rate of hot-spot heat loss by >5%. We demonstrate this via magnetic post-processing of two-dimensional xRAGE hydrodynamic simulation data at bang time. We then derive a model for the self-magnetization, finding that it varies with the square of the hot-spot temperature and inversely with the areal density. The self-magnetized Lawson analysis then gives a slightly reduced ignition threshold. Time-dependent hot-spot energy balance models corroborate this finding, with the magnetic field quadrupling the fusion yield for near-threshold parameters. The inclusion of magnetized multi-dimensional fluid instabilities could further alter the ignition threshold and will be the subject of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Simulating Atomic Dark Matter in Milky Way Analogs

Dark sector theories naturally lead to multicomponent scenarios for dark matter where a subcomponent can dissipate energy through self-interactions, allowing it to efficiently cool inside galaxies. We present the first cosmological hydrodynamical simulations of Milky Way analogs where the majority of dark matter is collisionless cold dark matter (CDM) but a subcomponent (6%) is strongly dissipative minimal atomic dark matter (ADM). The simulations, implemented in GIZMO and utilizing FIRE-2 galaxy formation physics to model the standard baryonic sector, demonstrate that the addition of even a small fraction of dissipative dark matter can significantly impact galactic evolution despite being consistent with current cosmological constraints. We show that ADM gas with roughly standard model–like masses and couplings can cool to form a rotating “dark disk” with angular momentum closely aligned with the visible stellar disk. The morphology of the disk depends sensitively on the parameters of the ADM model, which affect the cooling rates in the dark sector. The majority of the ADM gas gravitationally collapses into dark “clumps” (regions of black hole or mirror star formation), which form a prominent bulge and a rotating thick disk in the central galaxy. These clumps form early and quickly sink to the inner ∼kiloparsec of the galaxy, affecting the galaxy’s star formation history and present-day baryonic and CDM distributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Alignments between Galaxies and the Cosmic Web at z ~ 1–2 in the IllustrisTNG Simulations

Galaxy formation theories predict that galaxy shapes and angular momenta have nonrandom alignments with the cosmic web. This leads to so-called intrinsic alignment between pairs of galaxies, which is important to quantify as a nuisance parameter for weak lensing. We study galaxy–cosmic web alignment in the IllustrisTNG suite of hydrodynamical simulations at redshifts 1 and 2, finding that alignment trends are consistent with previous studies. However, we find that the magnitude of the spin alignment signal is ~2.4× weaker than seen in previous studies of the Horizon-AGN simulation, suggesting that this signal may have a significant dependence on subgrid physics. Based on IllustrisTNG, we then construct mock observational spectroscopic surveys that can probe shape–cosmic web alignment at z ~ 1–2, modeled on the low-z galaxy redshift and IGM tomography surveys on the upcoming Subaru Prime Focus Spectrograph Galaxy Evolution (PFS GE) survey. However, even over box sizes of L = 205 h -1 Mpc, we find that global anisotropies induce a sample variance in the 2D projected alignment signal that depend on the projected direction; this induces significant errors in the observed alignment. We predict a 5.3σ detection of IllustrisTNG's shape alignment signal at z ~ 1 from Subaru PFS GE, although a detection would be challenging at z ~ 2. However, a rough rescaling of the relative alignment signal strengths between the TNG and Horizon-AGN simulations suggests that PFS GE should be able to more easily constrain the latter's stronger signal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within ≲3σ of the mean of all simulations on a given grid (128 3 and 256 3 grid cells), where σ describes the statistical variation due to the flow’s time dependence. They also agree well with a 512 3 reference run. The 128 3 and 256 3 simulations agree within 9% and 4%, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Classical-quantum simulation of non-equilibrium Marshak waves

In the radiation hydrodynamic simulations used to design inertial confinement fusion (ICF) and pulsed power experiments, nonlinear radiation diffusion tends to dominate CPU time. This raises the interesting question of whether a quantum algorithm can be found for nonlinear radiation diffusion which provides a quantum speedup. Recently, such a quantum algorithm was introduced based on a quantum algorithm for solving systems of nonlinear partial differential equations (PDEs) which provides a quadratic quantum speedup. Here, we apply this quantum PDE (QPDE) algorithm to the problem of a non-equilibrium Marshak wave propagating through a cold, semi-infinite, optically thick target, where the radiation and matter fields are not assumed to be in local thermodynamic equilibrium. The dynamics is governed by a coupled pair of nonlinear PDEs which are solved using the QPDE algorithm, as well as two standard PDE solvers: (i) Python's py-pde solver; and (ii) the KULL ICF simulation code developed at Lawrence-Livermore National Laboratory. We compare the simulation results obtained using the QPDE algorithm and the standard PDE solvers and find excellent agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The impact of the $\text{WHIM}$ on the $\text{IGM}$ thermal state determined from the low- z Lyman $\alpha$ forest

At z ≲ 1, shock heating caused by large-scale velocity flows and possibly violent feedback from galaxy formation, converts a significant fraction of the cool gas (T ~ 10 4 K) in the intergalactic medium (IGM) into warm–hot phase (WHIM) with T > 10 5 K, resulting in a significant deviation from the previously tight power-law IGM temperature–density relationship, T = T 0 (P/$\bar{p}$) γ-1 ⁠. This study explores the impact of the WHIM on measurements of the low-z IGM thermal state, [T 0 , γ], based on the b–NH 1 distribution of the Ly α forest. Exploiting a machine learning-enabled simulation-based inference method trained on Nyx hydrodynamical simulations, we demonstrate that [T 0 , γ] can still be reliably measured from the b–NH 1 distribution at z = 0.1, notwithstanding the substantial WHIM in the IGM. To investigate the effects of different feedback, we apply this inference methodology to mock spectra derived from the IllustrisTNG and Illustris simulations at z = 0.1. The results suggest that the underlying [T 0 , γ] of both simulations can be recovered with biases as low as |Δlog(T 0 /K)| ≲ 0.05 dex, |Δγ| ≲ 0.1, smaller than the precision of a typical measurement. Given the large differences in the volume-weighted WHIM fractions between the three simulations (Illustris 38 percent, IllustrisTNG 10 percent, and Nyx 4 per cent), we conclude that the b–N H1 distribution is not sensitive to the WHIM under realistic conditions. Finally, we investigate the physical properties of the detectable Ly α absorbers, and discover that although their T and Δ distributions remain mostly unaffected by feedback, they are correlated with the photoionization rate used in the simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Agora : Multicomponent simulation for cross-survey science

ABSTRACT Next-generation cosmological surveys will observe large portions of the sky, with significant overlap between them. Multi-wavelength observations will enable us to analyse the same large-scale structure from different angles using a variety of tracers and astrophysical effects. The complex interplay between these observables calls for a model that can accurately and coherently describe their collective behaviour, posing a challenge that can only be met through the use of simulations. In this work, a suite of simulated extragalactic skies is presented, including maps and/or catalogues of cosmic microwave background (CMB) lensing, thermal and kinetic Sunyaev–Zel’dovich (tSZ/kSZ) effects, cosmic infrared background (CIB), radio sources, galaxy overdensity and galaxy weak lensing. Each of these probes is implemented in the lightcone using halo catalogues and/or particles from the Multidark-Planck2 ($\small {MDPL2}$) N-body simulation, and the modelling is calibrated using hydrodynamic simulations and publicly available data. The auto and cross-spectra of the individual probes, as well as the cross-spectra between the observables, are shown to be consistent with theoretical models and measurements from data. The simulation is shown to have a wide range of applications, including forecasting, pipeline testing, and evaluating astrophysical biases in cross-correlation measurements. It is further demonstrated that the simulation products produced in this work have sufficient accuracy to recover the input cosmology when subjected to a full cosmological analysis and are ready for application in real-world analyses for ongoing and future surveys. The simulation products presented in this work can be accessed at: https://yomori.github.io/agora/index.html.

Astronomy & Astrophysics↗
DOE OSTI2020

Detectability of collective neutrino oscillation signatures in the supernova explosion of a 8.8 M ⊙ star

In order to investigate the impact of collective neutrino oscillations (CNOs) on the neutrino signal from a nearby supernova, we perform three-flavor neutrino oscillation simulations employing the multiangle effect. The background hydrodynamic model is based on the neutrino hydrodynamic simulation of a 8.8 M ⊙ progenitor star. In this work, we find that CNO commences after some 100 ms post bounce. Before this, CNO is suppressed by matter-induced decoherence. In the inverted mass hierarchy, the spectrum of $\bar{ν}_e$ becomes softer after the onset of CNO. To evaluate the detectability of this modification, we define a hardness ratio between the number of high energy neutrino events and low energy neutrino events adopting a fixed critical energy. We show that Hyper-Kamiokande (HK) can distinguish the effect of CNO for supernova distances out to ~ 10 kpc. On the other hand, for the normal mass hierarchy, the spectrum of ν e becomes softer after the onset of CNO, and we show that DUNE can distinguish this feature for supernova distances out to ~ 10 kpc . More work is necessary to optimize the best value of critical energy for maximum sensitivity. We also show that if the spectrum of $\bar{ν}_e$ in HK becomes softer due to CNO, the spectrum of $ν_e$ in DUNE becomes harder, and vice versa. These synergistic observations in $\bar{ν}_e$ and $ν_e$, by HK and DUNE, respectively, will be an intriguing opportunity to test the occurrence of CNO.

79 ASTRONOMY AND ASTROPHYSICS↗