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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Numerical Simulation of a Single Shock-Accelerated Particle

Particle drag models, which capture macro viscous and pressure effects, have been developed over the years for various flow regimes to enable cost effective simulations of particle-laden flows. The relatively recent derivation by Maxey and Riley has provided an exact equation of motion for spherical particles in a flow field based on the continuum assumption. Many models that have been simplified from these equations have provided reasonable approximations; however, the sensitivity of particle-laden flows to particle drag requires a very accurate model to simulate. To develop such a model, a 2D axisymmetric Navier-Stokes direct numerical simulation of a single particle in a transient, shock-driven flow field was conducted in the hydrocode FLAG. FLAG’s capability to run arbitrary Lagrangian-Eulerian hydrodynamics coupled with solid mechanic models makes it an ideal code to capture the physics of the flow field around and in the particle as it is shock-accelerated – a challenging regime to study. The goal of this work is twofold: to provide a validation for FLAG’s Navier-Stokes and heat diffusion solutions, and to provide a rationale for recent experimental particle drag measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Evidence for large baryonic feedback at low and intermediate redshifts from kinematic Sunyaev-Zel’dovich observations with ACT and DESI photometric galaxies

Recent advances in cosmological observations have provided an unprecedented opportunity to investigate the distribution of baryons relative to the underlying matter. In this work, we show that the gas is more extended than the dark matter, and the amount of baryonic feedback at $z \lesssim 1$ disfavors low-feedback models such as that of state-of-the-art hydrodynamical simulation IllustrisTNG compared with high-feedback models such as that of the original Illustris simulation. This has important implications for bridging the gap between theory and observations and understanding galaxy formation and evolution. Furthermore, a better grasp of the baryon-dark matter link is critical to future cosmological analyses, which are currently impeded by our limited knowledge of baryonic feedback. Here, we measure the kinematic Sunyaev-Zel'dovich (kSZ) effect from the Atacama Cosmology Telescope (ACT), stacked on the luminous red galaxy (LRG) sample of the Dark Energy Spectroscopic Instrument (DESI) imaging survey. This is the first analysis to use photometric redshifts for reconstructing galaxy velocities. Due to the large number of galaxies comprising the DESI imaging survey, this is the highest signal-to-noise stacked kSZ measurement to date: we detect the signal at 13$σ$, finding strong evidence that the gas is more spread out than the dark matter, as well as a preference for larger feedback compared to some commonly used state-of-the-art hydrodynamical simulations. Here, our work opens up the possibility of recalibrating large hydrodynamical simulations using the kSZ effect. In addition, our findings highlight the importance of properly accounting for baryonic feedback with future surveys such as LSST through direct probes such as the kSZ, and shed light on long-standing enigmas in astrophysics, such as the “missing baryon” problem.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Revisiting Bayesian constraints on the transport coefficients of QCD

Multistage models based on relativistic viscous hydrodynamics have proven successful in describing hadron measurements from relativistic nuclear collisions. These measurements are sensitive to the shear and the bulk viscosities of QCD and provide a unique opportunity to constrain these transport coefficients. Bayesian analyses can be used to obtain systematic constraints on the viscosities of QCD, through methodical model-to-data comparisons. In this manuscript, we discuss recent developments in Bayesian analyses of heavy ion collision data. Furthermore, we highlight the essential role of closure tests in validating a Bayesian analysis before comparison with measurements. We discuss the role of the emulator that is used as proxy for the multistage theoretical model. We use an ongoing Bayesian analysis of soft hadron measurements by the JETSCAPE Collaboration as context for the discussion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Design and evaluation of a dilute flow particle-to-air heat exchanger for energy storage applications

The use of inert and redox-active particles for high-temperature energy storage requires the development of components that can efficiently transfer energy to high-pressure working fluids like supercritical carbon dioxide (sCO 2 ). Dilute flow reactors can enable high working fluid outlet temperatures and minimal parasitic losses compared to moving packed bed and fluidized bed reactors. This research uses both computational and experimental methods to explore the design trade-offs and practical challenges of a novel component for transferring energy from dilute flows of hot, reduced metal oxide (MO x ) particles to sCO 2 in tubes. A discretized thermal resistance network model, which accounts for particle hydrodynamics, multi-mode heat transfer, and reaction equilibrium, guides the design of a prototype device. This device is experimentally tested with a surrogate heat transfer fluids and inert particle temperatures up to 400°C and a heat duty exceeding 1 kW. The data are used to validate the thermal hydraulic sub-models, allowing for the simulation of reacting particle scenarios. Under nominal design conditions, the flow rate of reactive particles is predicted to be 30% lower than that of inert particles for the same energy recovered, with over 70% of the stored particle energy transferred to the sCO 2 . Furthermore, these findings can inform the design of more efficient energy recovery reactors for particle-based systems and can be integrated into system-level concentrated solar power models with thermal storage to optimize operating conditions.

14 SOLAR ENERGY↗
DOE OSTI2022

Dynamic responses of submerged components in advanced reactors: experimental and numerical studies

The seismic design of an advanced nuclear reactor must consider the interaction of vessel internal components with the surrounding coolant: fluid–structure interaction (FSI). Available analytical solutions for FSI of submerged components do not accommodate multiple-component, intense seismic inputs and complex reactor and internal geometries. Physical testing of reactor vessels and internals to inform seismic design is impractical and cost-prohibitive, leaving the use of verified and validated, robust numerical models as the only plausible option for analysis and design. Physical data that could be used for validating such numerical models for multi-component shaking are not available. This article describes an experimental program performed on a 6-degree-of-freedom earthquake simulator to generate data that could support validation of seismic FSI numerical models for submerged components in commercial finite element codes. A scaled model of a base-supported reactor vessel with simplified representations of submerged internals was tested to generate submerged-component response histories for a range of seismic inputs. The generated data were used to validate numerical models in the finite element code LS-DYNA. Numerical models were validated for calculation of hydrodynamic pressure responses on internals, in-water frequencies of internals, and acceleration responses of internals. The generated data and the analysis recommendations could aid engineering analysts designing submerged components and systems for seismic effects.

Engineering↗
DOE OSTI2020

Limitations to the ‘basic’ HOD model and beyond

ABSTRACT We use the IllustrisTNG cosmological, hydrodynamical simulations to test fundamental assumptions of the mass-based halo occupation distribution (HOD) approach to modelling the galaxy–halo connection. By comparing the clustering of galaxies measured in the 300 Mpc TNG box (TNG300) with that predicted by the standard (basic) HOD model, we find that, on average, the ‘basic’ HOD model underpredicts the real-space correlation function in the TNG300 box by ∼15 per cent on scales of $1 \,\,\lt\,\, r \,\,\lt\,\, 20 \ {\rm Mpc}\, h^{-1}$, which is well beyond the target precision demanded of next-generation galaxy redshift surveys. We perform several tests to establish the robustness of our findings to systematic effects, including the effect of finite box size and the choice of halo finder. In our exploration of ‘secondary’ parameters with which to augment the ‘basic’ HOD, we find that the local environment of the halo, the velocity dispersion anisotropy, β, and the product of the half-mass radius and the velocity dispersion, σ2Rhalfmass, are the three most effective measures of assembly bias that help reconcile the ‘basic’ HOD-predicted clustering with that in TNG300. In addition, we test other halo properties such as halo spin, formation epoch, and halo concentration. We also find that at fixed halo mass, galaxies in one type of environment cluster differently from galaxies in another. We demonstrate that a more complete model of the galaxy–halo connection can be constructed if we combine both mass and local environment information about the halo.

42 ENGINEERING↗
DOE OSTI2026

Cholla Galactic OutfLow Simulations (CGOLS)

These datasets contain full hydro-field snapshots from the galactic outflow simulations in the CGOLS suite, models I-V. The datasets were generated using the Cholla hydrodynamics code (https://github.com/cholla-hydro/cholla); descriptions of the models are in the associated publications (Schneider & Robertson 2018, ApJ; Schneider et al. 2018, ApJ; Schneider et al. 2020, ApJ; and Schneider & Mao, 2024, ApJ). Each hdf5 dataset is numbered according to the simulation time of the snapshot, in Myr. Fields include density, x momentum, y momentum, z momentum, total energy, and thermal energy (for models I - III), as well as a passive scalar field (models IV and V). 2 dimensional density and temperature projections, as well as slices along each midplane are also included if they exist.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Future Response of Coastal Wetlands to Environmental Stresses: Algorithm Comparison of Numerical Models

Coastal wetlands are a critical component of the earth system that strongly influence the global water and biogeochemical cycles. They are also likely important sentinel of climate change. Because eco-geomorphological processes have long been recognized to be important for coastal wetland survival under accelerated sea-level rise (SLR), many eco-geomorphology models have been developed to assess the impact of climate change on coastal wetlands. Although these models differ substantially in complexity and numerical methods, few studies have investigated the algorithm-level uncertainties in these models. In this study, we developed a multiple-algorithm model framework of coastal wetlands that represents coastal hydrodynamics (such as water level, significant wave height and bottom shear stress) and four eco-geomorphological processes: mineral accretion, organic matter accretion, storm surge erosion and landward migration. We validated the model at three representative coastal wetland sites (Venice Lagoon, Plum Island Estuary and Hunter Estuary) for hydrodynamics, mineral accretion and organic matter accretion. Through model-data comparison, we showed that the model can well capture the dynamics of hydrodynamical and eco-geomorphological conditions in the study sites. Importantly, analysis of the multiple-algorithm simulations suggests that differences in the process representation of mineral and organic matter accretion may contribute to the recent contradicting predictions of coastal wetland evolution under accelerated SLR.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Viscous flow properties and hydrodynamic diameter of phenothiazine-based redox-active molecules in different supporting salt environments

We report viscous flow properties of a redox-active organic molecule, N-(2-(2-methoxyethoxy)ethyl)phenothiazine (MEEPT), a candidate for non-aqueous redox flow batteries, and two of its radical cation salts. A microfluidic viscometer enabled the use of small sample volumes in determining viscosity as a function of shear rate and concentration in the non-aqueous solvent, acetonitrile, both with and without supporting salts. All solutions tested show Newtonian behavior over shear rates of up to 30 000 s−1, which was rationalized by scaling arguments for the diffusion-based relaxation time of a single MEEPT molecule without aggregation. Neat MEEPT is flowable but with a large viscosity (412 mPa⋅s at room temperature), which is ∼1000 times larger than that of acetonitrile. MEEPT solutions in acetonitrile have low viscosities; at concentrations up to 0.5 M, the viscosity increases by less than a factor of two. From concentration-dependent viscosity measurements, molecular information was inferred from intrinsic viscosity (hydrodynamic diameter) and the Huggins coefficient (interactions). Model fit credibility was assessed using the Bayesian Information Criterion. It is found that the MEEPT and its charged cations are “flowable” and do not flocculate at concentrations up to 0.5 M. MEEPT has a hydrodynamic diameter of around 8.5 Å, which is almost insensitive to supporting salt and state of charge. This size is comparable to molecular dimensions of single molecules obtained from optimized structures using density functional theory calculations. The results suggest that MEEPT is a promising candidate for redox flow batteries in terms of its viscous flow properties.

Wang, Yilin (ORCID:0000000316102995)↗
DOE OSTI2022

Three-dimensional hot-spot x-ray emission tomography from cryogenic deuterium–tritium direct-drive implosions on OMEGA

A three-dimensional model of the hot-spot x-ray emission has been developed and applied to the study of low-mode drive asymmetries in direct-drive inertial confinement fusion implosions on OMEGA with cryogenic deuterium–tritium targets. Here, the steady-state model assumes an optically thin plasma and the data from four x-ray diagnostics along quasi-orthogonal lines of sight are used to obtain a tomographic reconstruction of the hot spot. A quantitative analysis of the hot-spot shape is achieved by projecting the x-ray emission into the diagnostic planes and comparing this projection to the measurements. The model was validated with radiation-hydrodynamic simulations assuming a mode-2 laser illumination perturbation resulting in an elliptically shaped hot spot, which was accurately reconstructed by the model using synthetic x-ray images. This technique was applied to experimental data from implosions in polar-direct-drive illumination geometry with a deliberate laser-drive asymmetry, and the hot-spot emission was reconstructed using spherical-harmonic modes of up to ℓ = 3. A 10% stronger drive on the equator relative to that on the poles resulted in a prolate-shaped hot spot at stagnation with a large negative A 2,0 coefficient of A 2,0 = –0.47 ± 0.03, directly connecting the modal contribution of the hot-spot shape with the modal contribution in laser-drive asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The benefits of diligence: how precise are predicted gravitational wave spectra in models with phase transitions?

Models of particle physics that feature phase transitions typically provide predictions for stochastic gravitational wave signals at future detectors and such predictions are used to delineate portions of the model parameter space that can be constrained. The question is: how precise are such predictions? Uncertainties enter in the calculation of the macroscopic thermal parameters and the dynamics of the phase transition itself. We calculate such uncertainties with increasing levels of sophistication in treating the phase transition dynamics. Currently, the highest level of diligence corresponds to careful treatments of the source lifetime; mean bubble separation; going beyond the bag model approximation in solving the hydrodynamics equations and explicitly calculating the fraction of energy in the fluid from these equations rather than using a fit; and including fits for the energy lost to vorticity modes and reheating effects. The lowest level of diligence incorporates none of these effects. We compute the percolation and nucleation temperatures, the mean bubble separation, the fluid velocity, and ultimately the gravitational wave spectrum corresponding to the level of highest diligence for three explicit examples: SMEFT, a dark sector Higgs model, and the real singlet-extended Standard Model (xSM). In each model, we contrast different levels of diligence in the calculation and find that the difference in the final predicted signal can be several orders of magnitude. Our results indicate that calculating the gravitational wave spectrum for particle physics models and deducing precise constraints on the parameter space of such models continues to remain very much a work in progress and warrants care.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Turbulence modeling to aid tidal energy resource characterization in the Western Passage, Maine, USA

Numerical models combined with field measurements are regularly used to characterize tidal energy resources at potential energetic sites. However, most existing works only focus on the tidal hydrodynamic characteristics, and turbulence parameters are often not reported because of the lack of high-quality turbulence measurements and the limitations of numerical models in resolving turbulent eddies. In this study, we used FVCOM - a hydrostatic primitive equation (HPE) model - to characterize the tidal energy resource in the Western Passage, Maine, USA, by taking care of the essential macro-scale turbulence properties. We observed an excellent model performance using the Mellor-Yamada Level 2.5 Turbulence Model; estimating the spatial and vertical distribution of the turbulent kinetic energy and intensity added a new perspective to the site ranking for tidal energy converter (TEC) deployments. In addition, we also examined the role of channel geometry and bathymetry, such as headlands and underwater sills, in enhancing turbulent eddies around potential TEC siting locations. Ultimately, the detailed analysis of the turbulent flow characteristics has changed the site-ranking results and demonstrated that the regional-scale HPE models could be used for the relative understanding of more or less turbulent sites for a refined resource assessment.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Modeling of a spatially resolved ion temperature diagnostic for inertial confinement fusion

Here, the performance of modern laser-driven inertial confinement fusion (ICF) experiments is degraded by contamination of the deuterium–tritium (DT) fuel with high-Z material during compression. Simulations suggest that this mix can be described by the ion temperature distribution of the implosion, given that such contaminants deviate in temperature from the surrounding DT plasma. However, existing neutron time-of-flight (nTOF) diagnostics only measure the spatially integrated ion temperature. This paper describes the techniques and forward modeling used to develop a novel diagnostic imaging system to measure the spatially resolved ion temperature of an ICF implosion for the first time. The technique combines methods in neutron imaging and nTOF diagnostics to measure the ion temperature along one spatial dimension at yields currently achievable on the OMEGA laser. A detailed forward model of the source and imaging system was developed to guide instrument design. The model leverages neutron imaging reconstruction algorithms, radiation hydrodynamics and Monte Carlo simulations, optical ray tracing, and more. The results of the forward model agree with the data collected on OMEGA using the completed diagnostic. The analysis of the experimental data is still ongoing and will be discussed in a separate publication.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Hydrological Controls on Riverbed Methane Emissions: A Numerical Investigation of Hydrodynamic and Ebullitive Mechanisms from Site to Basin Scales

Methane production and emission from riverbed sediments constitute a significant yet underexplored component of the river methane budget. Despite their importance, a mechanistic and quantitative understanding of these processes under dynamic flow conditions remains elusive, particularly at the basin scale. In this study, we investigate the multiscale (site-to-basin) mechanisms governing methane emissions from riverbeds, focusing on the interplay between biogeochemical processes and hydrological dynamics. We developed a reactive transport model that integrates methane production with hydrodynamic and ebullitive pathways and coupled it with a basin-scale groundwater flow model to simulate methane emissions across scales. By incorporating key factors such as vertical hydrological exchange flow (VHEF), particulate organic carbon availability, sediment properties, and temperature sensitivity, our model captures the spatiotemporal variability of methane ebullition at both the site and basin scales. Our results reveal that methane ebullition is strongly modulated by VHEF, with lower flow rates promoting methane accumulation and subsequent ebullition. Sensitivity analyses underscore nonlinear responses of methane fluxes to environmental drivers. Furthermore, these findings emphasize the critical roles of sediment composition, hydrological exchange, and environmental conditions in shaping methane dynamics. Overall, this work advances our understanding of riverine methane emissions and offers insights for guiding climate models and mitigation strategies.

Fluxes↗
DOE OSTI2024

A hybrid finite volume method and smoothed particle hydrodynamics approach for efficient and accurate blast simulations

Modeling strong shock waves in fluids remains a persistent challenge in computational physics. Essential to research efforts in industry and defense, numerous methods have been devised to improve the accuracy and efficiency of shock simulations. A novel, hybrid Finite Volume Method (FVM)-Smoothed Particle Hydrodynamics (SPH) approach is capable of further improving efficiency and retaining accuracy by exploiting the favorable characteristics of each respective method. This hybrid approach is presented for shock capturing in compressible fluids. The Python framework Pyro2 is employed to simulate a coarse FVM mesh, while the Python framework PySPH is utilized to model the fluid in regions with high gradients through SPH particles. The performance of the hybrid FVM-SPH scheme, compared to the individual FVM and SPH methods, is assessed in 1 kt and 10 kt blast simulations. Our results indicate that the hybrid approach offers higher computational efficiency than SPH while preserving its accuracy and characteristics. The hybrid approach had a relative speedup of 11.3x and 22.3x over the FVM and SPH approaches for the 1 kt simulation and a relative speedup of 14.7x and 20.9x over the FVM and SPH approaches for the 10 kt simulation. The hybrid SPH algorithm enables future compressible fluid simulations with more extensive capabilities than grid-based methods alone, presenting potential applications in modeling fluid-structure interactions and solid deformation and fracturing in blast simulations.

Myers, Conner↗
DOE OSTI2024

Development of spatially and temporally resolved electron density measurements for the assessment of radiation hydrodynamics simulations of laboratory X-ray photoionized plasmas

The photoionized plasma gas cell experiment is an established platform we use to make at-parameter (ξ > > 1 ergs cm s –1 ) measurements of plasma properties with application to high-energy astrophysical systems. We model the experiments with 1D radiation hydrodynamics simulations using the HELIOS-CR code to inform our understanding and assist in the interpretation of results. The simulations predict that the bulk of the plasma is in a quasi-uniform and hydrodynamically unperturbed state throughout the duration of the experiment. To evaluate this prediction, we introduced a photonic Doppler velocimetry (PDV) diagnostic to measure spatially and temporally resolved plasma electron density. The initial measurements were successful but had limitations that made model-data comparisons challenging. To address this, we re-designed the gas cell PDV diagnostic and doubled the number of measurement locations to sample across two thirds of the depth of the cell. We also present a comparison of the results from the upgraded PDV diagnostic to the HELIOS-CR simulations for the first time. As a result, the experimental data confirms the prediction of an unperturbed region in the bulk of the plasma but reveals discrepancies in the time evolution and spatial distribution of the simulated electron density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Number of constituent quark scaling of elliptic flows in high multiplicity p-Pb collisions at $\sqrt{s_{NN}}$=5.02TeV

In this work, we briefly summarize our recent study on the number of constituent quark (NCQ) scaling of hadron elliptic flows in high multiplicity p-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. With the inclusion of hadron production via the quark coalescence model at intermediate p T , the viscous hydrodynamics at low p T , and jet fragmentation at high p T , our Hydro – Coal – Frag model provides a nice description of the p T -spectra and differential elliptic flow v 2 (p T ) of pions, kaons and protons over the pT range from 0 to 6 GeV. Our results demonstrate that including the quark coalescence is essential for reproducing the observed approximate NCQ scaling of hadron v 2 at intermediate p T in experiments, indicating strongly the existence of partonic degrees of freedom and the formation of quark-gluon plasma in high multiplicity p-Pb collisions at the LHC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗