Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Bayesian reconstruction of anisotropic flow fluctuations at fixed impact parameter

The cumulants of the distribution of anisotropic flow are measured accurately in Pb+Pb collisions at the LHC as a function of centrality classifiers (charged multiplicity and/or transverse energy). Using Bayesian inference, we reconstruct from these measurements the probability distribution of anisotropic flow in the ``theorists' frame'' where the impact parameter has a fixed magnitude and orientation, up to ∼70% centrality. The variation of flow fluctuations with impact parameter displays direct evidence of viscous damping, which is larger for higher Fourier harmonics, in line with expectations from hydrodynamics. We use intensive measures of non-Gaussian flow fluctuations, which have reduced dependence on centrality. Here, we infer from ATLAS data the magnitude of these intensive non-Gaussianities in each Fourier harmonic. They provide data-driven estimates of response coefficients to initial anisotropies, without resorting to any specific microscopic model of initial conditions. These estimates agree with viscous hydrodynamic calculations.

Bayesian methods↗
DOE OSTI2023

Pre-equilibrium evolution of conserved charges with initial conditions in the ICCING Monte Carlo event generator

Heavy-ion collisions can be well described through relativistic viscous hydrodynamics, but questions still remain when hydrodynamics is applicable because the initial state may begin very far from equilibrium. Thus, a pre-equilibrium evolution phase is used to bridge the gap between the initial state and hydrodynamics. KøMPøST is one such pre-equilibrium model that propagates the energy-momentum tensor by decomposing it into the background and fluctuations around that background, whose evolution is captured by Green's functions. We extend this formalism to include conserved charges and calculate the corresponding nonequilibrium Green's functions in the relaxation-time approximation. The ICCING algorithm initializes conserved charges in the initial state by sampling $g$ → $q$$\overline{q}$ splitting probabilities and is, thus, perfectly positioned to implement Green's functions for charge propagation. As a result, we show that this method alters the initial-state charge geometries and is applicable in central to mid-central collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Development, Verification, and Validation of an OpenFOAM-Based Solver for Modeling Inertial Fusion Energy Chambers

Our work seeks to introduce a computational tool tailored to the physics of inertial fusion energy chambers, in particular, those concepts based on thick liquid walls. In this approach, the structural materials are protected by several neutron mean-free-paths of renewable liquid and thus will be able to survive much longer than un-shielded walls, with virtually all structures lasting for the life of the plant and enabling the use of commercially available and qualified materials. The OpenFOAM-based solver named rhoCentralFoam has been used as a starting point. rhoCentralFoam belongs to the standard OpenFOAM solver toolset. It is a high-speed, explicit compressible flow solver with shock-capturing capability. While the main features have been retained, the solver had to be restructured to make use of tabular data for equations of states, a necessary addition to model the complex thermo-physical properties of ionized gasses. This entailed the need to change the independent state variables used by the solver, resulting in a new thermodynamic library and slightly different solution algorithm. Moreover, a radiation heat transfer model based on the P-1 approximation was added to the solver. The solver is verified against an analytical solution from the Sedov-Taylor-Neumann test problem to showcase the ability of the hydrodynamic solvers to handle strong shocks, whereas the P-1 model was verified using a simple one-dimensional problem with an analytical solution. Additionally, a validation case involving shock-wave propagation through jet array is presented, and the results are compared with experimental data from the open literature. Lastly, in order to showcase the utility of the solver for practical cases, we applied the refined solver to two representative scenarios: gas venting within the HYLIFE-II chamber and the compression of the gas following the partial ablation of the liquid wall.

Chamber dynamics↗
DOE OSTI2023

Exploration of cross-beam energy transfer mitigation constraints for designing an ignition-scale direct-drive inertial confinement fusion driver

The compression of direct-drive inertial confinement fusion (ICF) targets is strongly impacted by cross-beam energy transfer (CBET), a laser-plasma instability that limits ablation pressure by redirecting laser energy outward and that is projected to be mitigated by laser bandwidth. Here, we explore various CBET mitigation constraints to guide the design of future ICF facilities. First, we find that the flat, Gaussian, and Lorentzian spectral shapes have similar CBET mitigation properties, and a flat shape with nine spectral lines is a good surrogate for what can be obtained with other spectral shapes. Then, we conduct a comprehensive study across energy scales and ignition designs. 3D hydrodynamic simulations are used to derive an analytical model for the expected CBET mitigation as a function of laser and plasma parameters. From this model, we study the bandwidth requirements of conventional and shock ignition designs across four different energy scales and find that they require between 0.5 and 3±0.2% relative bandwidth. Further, best mitigation is achieved when the beam radius over critical radius R b /R c is kept low during the drive while the plasma temperature is kept high. In a steady state, we find that the bandwidth required to mitigate 85% of CBET scales as (R b /R c ) 2.15 L n –0.58 I 0.7 , where L n is the density scale length, and I the laser intensity. Finally, we find that the chamber beam port layout does not influence CBET mitigation. In the case of a driver using many monochromatic beamlets, we find that ~10 beamlets per port is required, with diminishing returns above ~20.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Code-to-code comparison and validation of the radiation-hydrodynamics capabilities of the FLASH code using a laboratory astrophysical jet

The potential for laser-produced plasmas to yield fundamental insights into high energy density physics (HEDP) and deliver other useful applications can sometimes be frustrated by uncertainties in modeling the properties and behavior of these plasmas using radiation-hydrodynamics codes. In an effort to overcome this and to corroborate the accuracy of the HEDP capabilities in the publicly available FLASH radiation-hydrodynamics code, we present detailed code-to-code comparisons between FLASH and the HYDRA code developed at Lawrence Livermore National Laboratory using previously published HYDRA simulations from Grava et al. [Phys. Rev. E 78, 016403 (2008)]. That study describes a laser experiment that produced a jet-like feature that the authors compare to astrophysical jets. Importantly, the Grava et al. [Phys. Rev. E 78, 016403 (2008)] experiment included detailed x-ray interferometric measurements of electron number densities and a time-integrated measurement of the soft x-ray spectrum. Despite markedly different methods for treating the computational mesh, and different equations of state and opacity models, the FLASH results resemble the results from HYDRA and, most importantly, the experimental measurements of electron density. Having validated the FLASH code in this way, we use the code to further investigate and understand the formation of the jet seen in the Grava et al. [Phys. Rev. E 78, 016403 (2008)] experiment and discuss its relation to the Wan et al. [Phys. Rev. E 55, 6293 (1997)] experiment at the NOVA laser.

Orban, Chris (ORCID:0000000312004538)↗
DOE OSTI2021

A statistical mechanics approach to macroscopic limits of car-following traffic dynamics

Herein we study the derivation of macroscopic traffic models from car-following vehicle dynamics by means of hydrodynamic limits of an Enskog-type kinetic description. We consider the superposition of Follow-the-Leader (FTL) interactions and relaxation towards a traffic-dependent Optimal Velocity (OV) and we show that the resulting macroscopic models depend on the relative frequency between these two microscopic processes. If FTL interactions dominate then one gets an inhomogeneous Aw–Rascle–Zhang model, whose (pseudo) pressure and stability of the uniform flow are precisely defined by some features of the microscopic FTL and OV dynamics. Conversely, if the rate of OV relaxation is comparable to that of FTL interactions then one gets a Lighthill–Whitham–Richards model ruled only by the OV function. We further confirm these findings by means of numerical simulations of the particle system and the macroscopic models. Unlike other formally analogous results, our approach builds the macroscopic models as physical limits of particle dynamics rather than assessing the convergence of microscopic to macroscopic solutions under suitable numerical discretisations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Measurement of beauty-strange meson production in Pb–Pb collisions at s NN = 5.02 TeV via non-prompt D s + mesons

The production yields of non-prompt ${D}^{+}_{s}$ mesons, namely ${D}^{+}_{s}$ mesons from beauty-hadron decays, were measured for the first time as a function of the transverse momentum ($p_T$) at midrapidity (|$y$| < 0.5) in central and semi-central Pb–Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s{NN}}$ = 5.02 TeV with the ALICE experiment at the LHC. The ${D}^{+}_{s}$ mesons and their charge conjugates were reconstructed from the hadronic decay channel ${D}^{+}_{s}$ → $φπ^+$, with $φ$ → K - K + , in the 4 < $p_T$ < 36 GeV/c and 2 < $p_T$ < 24 GeV/ $c$ intervals for the 0–10% and 30–50% centrality classes, respectively. The measured yields of non-prompt ${D}^{+}_{s}$ mesons are compared to those of prompt ${D}^{+}_{s}$ and non-prompt $D^0$ mesons by calculating the ratios of the production yields in Pb–Pb collisions and the nuclear modification factor $R_{AA}$. The ratio between the $R_{AA}$ of non-prompt ${D}^{+}_{s}$ prompt ${D}^{+}_{s}$ mesons, and that between the $R_{AA}$ of non-prompt ${D}^{+}_{s}$ and non-prompt $D^0$ mesons in central Pb–Pb collisions are found to be on average higher than unity in the 4 < $p_T$ < 12 GeV/ $c$ interval with a statistical significance of about 1.6σ and 1.7σ, respectively. The measured $R_{AA}$ ratios are compared with the predictions of theoretical models of heavy-quark transport in a hydrodynamically expanding QGP that incorporate hadronisation via quark recombination.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Multi-messenger dynamic imaging of laser-driven shocks in water using a plasma wakefield accelerator

Understanding dense matter hydrodynamics is critical for predicting plasma behavior in environments relevant to laser-driven inertial confinement fusion. Traditional diagnostic sources face limitations in brightness, spatiotemporal resolution, and in their ability to detect relevant electromagnetic fields. In this work, we present a dual-probe, multi-messenger laser wakefield accelerator platform combining ultrafast X-rays and relativistic electron beams at 1 Hz, to interrogate a free-flowing water target in vacuum, heated by an intense 200 ps laser pulse. This scheme enables high-repetition-rate tracking the evolution of the interaction using both particle types. Betatron X-rays reveal a cylindrically symmetric shock compression morphology assisted by low-density vapor, resembling foam-layer-assisted fusion targets. The synchronized electron beam detects time-evolving electromagnetic fields, uncovering charge separation and ion species differentiation during plasma expansion – phenomena not captured by photons or hydrodynamic simulations. We show that combining both probes provides complementary insights spanning kinetic to hydrodynamic regimes, highlighting the need for hybrid physics models to accurately predict fusion-relevant plasma behavior.

imaging techniques↗
DOE OSTI2021

A simple model of a strong shock driven by a spherical or cylindrical piston

A simple model of piston-driven spherical and cylindrical shocks is suggested. The model is based on a consistent use of two factors: (a) an almost uniform pressure across the shocked layer and, (b) continuous geometrical stretching of the surface elements of the expanding piston. It turns out that for a uniform pre-shock medium the gas between the piston and the shock behaves essentially as an incompressible fluid. An algebraic equation for the shock vs piston position is obtained. Detailed evaluation of the accuracy of the proposed solution shows that its accuracy is a few percent for the adiabatic index γ = 5/3 (as in ideal plasma). A closed-form solution describing enhancement of a weak ambient magnetic field by the shock is presented. The proposed model of piston-driven shocks goes beyond the classical self-similar solutions in that it: (1) naturally covers an early, non-asymptotic dynamics and its transition to asymptotic regime; (2) allows for smooth radial density variation of an ambient gas of the form of bumps, dimples or ramps between two constant values; and (3) allows for smooth temporal variation of piston velocity of the form of bumps, dimples, or ramps. Furthermore, this simple and versatile model provides some new insights into a classical hydrodynamical problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Two-dimensional parametric design of short pulse laser driven iron opacity experiments

High intensity short pulse lasers are commonly used to create high temperature and high density plasmas. These plasmas are used to study fundamental properties such as the equation of state and opacity. The targets often use small-diameter, thin buried layers of the materials of interest, surrounded by a tamper of low atomic number material such as plastic. Computational modeling is used to design and interpret experiments with short pulse lasers. Most of the modeling to date has assumed one-dimensional plane-parallel geometry. In this paper, the effects of radial gradients in the irradiation of thin planar targets are studied with one- and two-dimensional radiation/hydrodynamic simulations. It is found that the main effect of radial gradients is the averaging of plasma conditions and x-ray emission over the pattern of irradiation. Additionally, differences between one- and two-dimensional simulations arise because the plasma conditions and x-ray emission are, in general, nonlinear functions of the temperature of the plasma and thereby also nonlinear functions of the irradiation intensity. The differences increase with the ratio of the buried layer radius to the laser spot radius. The root mean square difference in the inferred iron L-shell opacity is less than 30% when the ratio is less than 1.0 and about 70% when the ratio is 1.5.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

200-fold increase in dynamic strength of platinum at 430 GPa without a phase transformation

Here, the high-pressure strength of solid platinum is studied at the National Ignition Facility via Rayleigh-Taylor ripple growth under ramp compression. Laser-driven experiments reached pressures up to 430 GPa, collecting velocimetry and radiography data. Despite no phase transformation, the observed ripple growth aligns with hydrodynamic simulations using an amplified Steinberg-Guinan model, resulting in a 200-fold increase in strength. Molecular dynamics simulations of platinum ramp compression yield similar von Mises stress at comparable high pressure and high strain rate conditions. These findings represent the highest-pressure strength measurements ever achieved for platinum and offer critical insight into its significant strengthening behavior under high pressure and strain rate conditions.

Righi, Gaia [Lawrence Livermore National Laborator↗
DOE OSTI2023

3D structure of anisotropic flow in small collision systems at energies available at the BNL Relativistic Heavy Ion Collider

Here, we present (3 + 1)-dimensional [(3 + 1)D] dynamical simulations of asymmetric nuclear collisions at the BNL Relativistic Heavy Ion Collider (RHIC). Employing a dynamical initial state model coupled to (3 + 1)D viscous relativistic hydrodynamics, we explore the rapidity dependence of anisotropic flow in the RHIC small system scan at 200 GeV center-of-mass energy. We calibrate parameters to describe central 3 He + Au collisions and make extrapolations to d + Au and p + Au collisions. Our calculations demonstrate that approximately 50% of the v 3 (p T ) difference between the measurements by the STAR and PHENIX Collaborations can be explained by the use of reference flow vectors from different rapidity regions. This emphasizes the importance of longitudinal flow decorrelation for anisotropic flow measurements in asymmetric nuclear collisions, and the need for (3 + 1)D simulations. We also present results for the beam energy dependence of particle spectra and anisotropic flow in d + Au collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

The Atacama Cosmology Telescope: Probing the baryon content of SDSS DR15 galaxies with the thermal and kinematic Sunyaev-Zel’dovich effects

We present high signal-to-noise measurements (up to 12$\sigma$) of the average thermal Sunyaev Zel'dovich (tSZ) effect from optically selected galaxy groups and clusters and estimate their baryon content within a 2.1$^\prime$ radius aperture. Sources from the Sloan Digital Sky Survey (SDSS) Baryon Oscillation Spectroscopic Survey (BOSS) DR15 catalog overlap with 3,700 sq. deg. of sky observed by the Atacama Cosmology Telescope (ACT) from 2008 to 2018 at 150 and 98 GHz (ACT DR5), and 2,089 sq. deg. of internal linear combination component-separated maps combining ACT and $\it{Planck}$ data (ACT DR4). The corresponding optical depths, $\bar{\tau}$, which depend on the baryon content of the halos, are estimated using results from cosmological hydrodynamic simulations assuming an AGN feedback radiative cooling model. We estimate the mean mass of the halos in multiple luminosity bins, and compare the tSZ-based $\bar{\tau}$ estimates to theoretical predictions of the baryon content for a Navarro-Frenk-White profile. We do the same for $\bar{\tau}$ estimates extracted from fits to pairwise baryon momentum measurements of the kinematic Sunyaev-Zel'dovich effect (kSZ) for the same data set obtained in a companion paper. We find that the $\bar{\tau}$ estimates from the tSZ measurements in this work and the kSZ measurements in the companion paper agree within $1\sigma$ for two out of the three disjoint luminosity bins studied, while they differ by 2-3$\sigma$ in the highest luminosity bin. The optical depth estimates account for one third to all of the theoretically predicted baryon content in the halos across luminosity bins. Potential systematic uncertainties are discussed. The tSZ and kSZ measurements provide a step towards empirical Compton-$\bar{y}$-$\bar{\tau}$ relationships to provide new tests of cluster formation and evolution models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Separated reactant mix width across diffusion-dominated and hydrodynamically dominated interface mix in inertial confinement fusion implosions

Diffusion-dominated mix in inertial confinement fusion (ICF) is characterized where the majority of the mix occurs in the immediate fuel-shell interface while hydrodynamic-dominated mix pulls shell material from farther away into the central fuel. A thin (150 nm) separated reactants ICF mix platform is highly sensitive to the amount of mix from the first micron of shell-fuel interface. This fine-spatial resolution platform has revealed that material mix in moderate convergence ( CR ∼ 12 ) ICF implosions is dominated by a diffusion mechanism. This technique has now been expanded across a set of OMEGA ICF implosions, observing an increase in mix width and mix amount for cooler, slower, and more compressive implosions. Hydrodynamic simulations require a buoyancy-drag mix model to capture the increasing mix width, suggesting a transition between these two mix mechanisms. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Supernova Hydrodynamics: A Lab-scale Study of the Blast-driven Instability Using High-speed Diagnostics

A novel experimental approach to study the blast-driven instability at a nondiffuse, gaseous interface with a density gradient is presented. Under Euler similarity, this approach enables study of dissipative-scale hydrodynamics relevant to many astrophysical and laboratory high energy density phenomena in a well-resolved manner. The instability is initiated by passing a Taylor–Sedov blast wave originating from a controlled detonation through a perturbed and stably stratified interface between two gases. The facility and driving blast wave are characterized to obtain repeatable conditions and capture large ensembles of time-resolved Mie scattering imaging that show consistent hydrodynamic development. We analyze the instability evolution between different gas pairs to demonstrate the wide range of development and turbulent behavior that may occur between different supernova layers. The mean evolution of the hydrodynamic instability is compared to a buoyancy–drag model that is frequently used to estimate perturbation growth in supernova mixing research. We propose a time delay to this model in order to reproduce the measured instability behavior and demonstrate model robustness in handling flows driven by a time-varying acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Pre-main-sequence Brackett Emitters in the APOGEE DR17 Catalog: Line Strengths and Physical Properties of Accretion Columns

Very young (t ≲ 10 Myr) stars possess strong magnetic fields that channel ionized gas from the interiors of their circumstellar disks to the surface of the star. Upon impacting the stellar surface, the shocked gas recombines and emits hydrogen spectral lines. To characterize the density and temperature of the gas within these accretion streams, we measure equivalent widths of Brackett (Br) 11–20 emission lines detected in 1101 APOGEE spectra of 326 likely pre-main-sequence accretors. For sources with multiple observations, we measure median epoch-to-epoch line strength variations of 10% in Br11 and 20% in Br20. We also fit the measured line ratios to predictions of radiative transfer models by Kwan & Fischer. We find characteristic best-fit electron densities of n e = 10 11 –10 12 cm -3 , and excitation temperatures that are inversely correlated with electron density (from T ~ 5000 K for n e ~ 10 12 cm -3 to T ~ 12,500 K at n e ~ 10 11 cm -3 ). These physical parameters are in good agreement with predictions from modeling of accretion streams that account for the hydrodynamics and radiative transfer within the accretion stream. We also present a supplementary catalog of line measurements from 9733 spectra of 4255 Brackett emission-line sources in the APOGEE Data Release 17 data set.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Optimizing Power Generation of a Bottomraised Oscillating Surge Wave Energy Converter Using a Theoretical Model

Preliminary sizing of an oscillating surge wave energy converter (OSWEC) is an iterative process that relies on knowledge of the relevant hydrodynamic coefficients for a given geometry. Often, the complex definition of the device geometry requires coefficients to be obtained through experiments or numerical boundary element solvers such as WAMIT. These techniques demand significant user and computational effort, therefore inhibiting the fine scale parametrization of object dimensions. In this study, a theoretical model, originally presented in Michele et al. (2016), is developed and demonstrated to efficiently optimize the power production for an OSWEC device (subjected to certain environmental conditions) with variations in device widths, heights, and distances from the seabed. Assuming negligible device thickness, the OSWEC motions are modelled as a bottom-raised two-dimensional flap in regular waves using potential flow theory formulated in elliptical coordinates. The solutions to this diffraction-radiation problem are obtained using Mathieu functions with appropriate boundary conditions. The resulting potentials are then used to derive frequency-dependent expressions for the added mass and radiation damping coefficients, along with wave excitation magnitude in surge, pitch, and coupled surge-pitch motions. Good agreement in hydrodynamic coefficient curves is shown between the theoretical model and the numerical results obtained from the boundary element-based program WAMIT. The theoretical model is then employed to maximize the time-averaged output power while maintaining or reducing the hinge reaction force, with variations in device dimensions, wave frequency and amplitude.

ACE↗
DOE OSTI2026

Modeling the Cosmological Lyman-𝛼 Forest at the Field Level

The distribution of absorption lines in the spectra of distant quasars, called the Lyman-𝛼 (Ly-𝛼) forest, is a unique probe of cosmology and the intergalactic medium at high redshifts and small scales. The statistical power of ongoing redshift surveys demands precise theoretical tools to model the Ly-𝛼 forest. We address this challenge by developing an analytic, perturbative forward model to predict the Ly-𝛼 forest at the field level for a given set of cosmological initial conditions. Our model shows a remarkable performance when compared with the Sherwood hydrodynamic simulations: it reproduces the Ly-𝛼 forest flux power spectrum, its cross-correlation with dark matter halos, and the one-point probability distribution function of both fields at the percent level down to scales of a few Mpc. Our work provides crucial tools that bridge analytic modeling on large scales with simulations on small scales, enabling field-level inference from Ly-𝛼 forest data and simulation-based priors for cosmological analyses. Furthermore, this is especially timely for realizing the full scientific potential of the Ly-𝛼 forest measurements by the dark energy spectroscopic instrument.

Cosmological parameters↗