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.

SourceResearch connectionFrom yearThrough year

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

At least 649 records · Page 36

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Geometric Scale-up Experiments on Fluidization of Geldart B Glass Beads

The objective of this work is to provide a valuable database from controlled experiments for validating computational models. Recently, coarse-grained techniques such as particle-in-cell (PIC) or coarse-grained discrete element modeling (DEM) have gained popularity due to their computational efficiency while modeling large-scale systems; however, the influence of model parameters and their sensitivities at different geometric scales and flow conditions remain to be analyzed. These datasets are critical for the multiphase flow research community to assess predictive capability of modeling techniques as well as elucidate the hydrodynamic behavior in these systems. This study performed fluidization experiments using three different test sections with internal diameters of 2.5, 4, and 6 in. The operating conditions, bed material, and range of flow velocities at the inlet were constant in all the units, which were not hydrodynamically scaled. Glass beads having a Sauter Mean Diameter of 332 μm were used. Superficial velocity was varied from 2.97 to 5.35 times the minimum fluidization velocity and the initial static bed height was 0.1524 m. The order in which the experiments were performed was randomized and replicates were included to provide uncertainty in measurements. Statistics of differential pressure and bed height from these tests were reported. Future plans include validating PIC methodology in the open-source software, MFiX (Multiphase Flow with Interphase Exchanges) using results from this study. This could further be extended to determine optimal model parameters using inverse techniques such as deterministic calibration or Bayesian inference.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2020

Geometric Scale-up Experiments on Fluidization of Geldart B Glass Beads

42 ENGINEERING↗

DOE OSTI2025

Implementation of tritium transport in a gas-liquid contactor CFD simulation of tritium extraction from lead-lithium in ANSYS fluent

Modifications to the Computational Fluid Dynamic (CFD) software ANSYS Fluent were done to quantify and characterize tritium transport in Gas-Liquid Contactors (GLCs). A double-slit, Ergun-like equation was employed for the porous media model, with Ergun coefficients validated with Sulzer’s Sulcol software. Tritium transport from PbLi within the GLC was verified against analytical models. The geometry of the CFD model was based on the MELODIE GLC experiment. The hydrodynamic CFD pressure drop results align well with SulCol estimations and fall between the predictions of the analytical Delft-Olujic and Billet & Schultes models. In terms of mass transfer efficiency, traditional mass transfer models showed a significant deviation from experimental results when using varying values of H solubility in PbLi. A saturation phenomenon occurred when utilizing high solubility values for hydrogen in PbLi. In conclusion, a modified film theory mass transfer coefficient, incorporating either the Delft-Olujic or Billet & Schultes wettability model, yielded CFD-predicted extraction efficiencies that closely matched experimental measurements.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

The effects of multispecies Hohlraum walls on stimulated Brillouin scattering, Hohlraum dynamics, and beam propagation

We report experiments and simulations have been conducted to investigate the efficacy of Ta 2 O 5 -lined Hohlraum walls at reducing stimulated Brillouin backscattering (SBS) as well as any subsequent effects on the Hohlraum dynamics and capsule implosions in indirect drive experiments at the National Ignition Facility. Using a 1.1 MJ 400 TW, 351 nm, shaped laser pulse, we measure a 5× reduction in SBS power in the peak of the pulse from the wall on the outer 50° cone beams. The SBS spectrum indicates a reduction in the high-Z spectral signature when using multispecies wall materials. Detailed hydrodynamic simulations were performed using different heat conduction models with flux limiters. Additional simulations were performed on the plasma maps using the 3D parallel paraxial code pF3D to compare backscatter powers between the pure Au and Ta 2 O 5 -lined Hohlraums. Further analysis, using hydrodynamically equivalent plasmas, shows that the SBS reduction is clearly a result of the added ion Landau damping caused by the oxygen ions and not from differences in plasma conditions. The experimental and simulation results also show an increase in the wall plasma expansion when using the Ta 2 O 5 liner leading to a 70% more oblate implosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Impact of toroidal rotation on the resistive ballooning modes in ASDEX Upgrade tokamak

Abstract In this work, we investigate the behavior of instabilities appearing between type-I edge localized modes (ELMs), with increasing neutral beam injection (NBI) power concomitant increase in toroidal rotation, and compare it to the modeling result of the linear magneto-hydrodynamic (MHD) code CASTOR3D. An injection of one NBI beam, increasing toroidal rotation, results in the mode slowing down from 12 kHz to 7 kHz, and its associated radial displacement decreases from 5 mm to 3.5 mm. In addition, modes shift radially outwards towards higher q , decreasing their poloidal mode numbers. The mode velocity is measured to be close to the E × B velocity with significant uncertainties. Through a set of CASTOR3D simulations with varying profiles, resistivity has been identified as the primary contributor to the growth rates. Only a small stabilizing effect due to toroidal rotation has been observed. While experimental results show a decrease of mode frequency with rotation, the opposite trend is observed in modeling. Reasons for discrepancies between modeling and experiment are discussed. Nevertheless, a main contributor to the mode frequency has been identified to be rotation velocity. CASTOR3D classifies modes as resistive ballooning modes as they do not appear unstable in ideal MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

X-ray observation of individual Ti-6Al-4V spherical powder particle impact in an in-situ operando laser directed energy deposition system

Understanding the particle capture mechanism in the laser directed energy deposition additive manufacturing process provides a foundation for improving productivity and reducing material loss. Capturing this phenomenon is challenging due to its highly transient time interval. To analyze the impact behavior of spherical Ti-6Al-4V powder, a custom laser directed energy deposition system was designed to control the deposition of individual particles. Synchrotron X-ray imaging at 24 kHz monitored the impact dynamics of individual powder particles. A multi-physics model based on the Smoothed-Particle Hydrodynamics scheme supported the analysis by providing the temperature and velocity fields of the molten pool. Results revealed that in laser directed energy deposition, surface tension forces dominated the powder–molten pool interaction. Hydrophobic Ti-6Al-4V powder particles had an equilibrium contact angle of 121 degrees with their molten state. Furthermore, Ti-6Al-4V powder particles required a high impact velocity (6 to 12 m/s) to transition from oscillation to submergence behavior. Lastly, oscillating powder particles that impacted near the location of the laser beam exhibited faster melting, highlighting the contribution of the laser beam to the melting mechanism. This work lays the foundation for investigations of powder particle impact in laser directed deposition and supports the validation of numerical models of powder–molten pool interactions.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Dark matter cores and cusps in spiral galaxies and their explanations

We compare proposed solutions to the core vs cusp issue of spiral galaxies, which has also been framed as a diversity problem, and demonstrate that the cuspiness of dark matter halos is correlated with the stellar surface brightness. We compare the rotation curve fits to the SPARC sample from a self-interacting dark matter (SIDM) model, which self-consistently includes the impact of baryons on the halo profile, and hydrodynamical N-body simulations with cold dark matter (CDM). The SIDM model predicts a strong correlation between the core size and the stellar surface density, and it provides the best global fit to the data. The CDM simulations without strong baryonic feedback effects fail to explain the large dark matter cores seen in low surface brightness galaxies. On the other hand, with strong feedback, CDM simulations do not produce galaxy analogs with high stellar and dark matter densities, and therefore they have trouble in explaining the rotation curves of high surface brightness galaxies. This implies that current feedback implementations need to be modified. We also explicitly show how the concentration-mass and stellar-to-halo mass relations together lead to a radial acceleration relation (RAR) in an averaged sense, and reiterate the point that the RAR does not capture the diversity of galaxy rotation curves in the inner regions. These results make a strong case for SIDM as the explanation for the cores and cusps of field galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Introducing the DREAMS Project: DaRk mattEr and Astrophysics with Machine Learning and Simulations

We introduce the DaRk mattEr and Astrophysics with Machine learning and Simulations (DREAMS) project, an innovative approach to understanding the astrophysical implications of alternative dark matter (DM) models and their effects on galaxy formation and evolution. The DREAMS project will ultimately comprise thousands of cosmological hydrodynamic simulations that simultaneously vary over DM physics, astrophysics, and cosmology in modeling a range of systems—from galaxy clusters to ultra-faint satellites. Such extensive simulation suites can provide adequate training sets for machine-learning-based analyses. This paper introduces two new cosmological hydrodynamical suites of warm dark matter (WDM), each comprising 1024 simulations generated using the arepo code. One suite consists of uniform-box simulations covering a ${(25\,{h}^{-1}\,{\rm{Mpc}})}^{3}$ volume, while the other consists of Milky Way zoom-ins with sufficient resolution to capture the properties of classical satellites. For each simulation, the WDM particle mass is varied along with the initial density field and several parameters controlling the strength of baryonic feedback within the IllustrisTNG model. We provide two examples, separately utilizing emulators and convolutional neural networks, to demonstrate how such simulation suites can be used to disentangle the effects of DM and baryonic physics on galactic properties. The DREAMS project can be extended further to include different DM models, galaxy formation physics, and astrophysical targets. In this way, it will provide an unparalleled opportunity to characterize uncertainties on predictions for small-scale observables, leading to robust predictions for testing the particle physics nature of DM on these scales.

cosmological parameters↗

DOE OSTI2022

OC6 Phase Ia: CFD Simulations of the Free-Decay Motion of the DeepCwind Semisubmersible

Currently, the design of floating offshore wind systems is primarily based on mid-fidelity models with empirical drag forces. The tuning of the model coefficients requires data from either experiments or high-fidelity simulations. As part of the OC6 (Offshore Code Comparison Collaboration, Continued, with Correlation, and unCertainty (OC6) is a project under the International Energy Agency Wind Task 30 framework) project, the present investigation explores the latter option. A verification and validation study of computational fluid dynamics (CFD) models of the DeepCwind semisubmersible undergoing free-decay motion is performed. Several institutions provided CFD results for validation against the OC6 experimental campaign. The objective is to evaluate whether the CFD setups of the participants can provide valid estimates of the hydrodynamic damping coefficients needed by mid-fidelity models. The linear and quadratic damping coefficients and the equivalent damping ratio are chosen as metrics for validation. Large numerical uncertainties are estimated for the linear and quadratic damping coefficients; however, the equivalent damping ratios are more consistently predicted with lower uncertainty. Some difference is observed between the experimental and CFD surge-decay motion, which is caused by mechanical damping not considered in the simulations that likely originated from the mooring setup, including a Coulomb-friction-type force. Overall, the simulations and the experiment show reasonable agreement, thus demonstrating the feasibility of using CFD simulations to tune mid-fidelity models.

17 WIND ENERGY↗

DOE OSTI2026

Exploring composition mixing in kilonova ejecta with ray-by-ray simulations

Binary neutron star merger (BNSM) ejecta are considered a primary repository of r-process nucleosynthesis and a source of the observed heavy-element abundances. We implement composition mixing into ray-by-ray radiation-hydrodynamic simulations of BNSM ejecta, coupled with an online nuclear network (NN). We model mixing via a gradient-based mixing approximation that evolves simultaneously with the hydrodynamics. Here, we find that mixing occurs in regions where the electron fraction changes rapidly. While mixing smooths composition gradients in transition regions, it has a negligible impact on the heavy-element yields. This is because the primary r-process site (the equatorial ejecta) is initially homogeneous in free neutrons, leaving no strong gradients for mixing to act upon. In each angular ray, the abundances of the most produced elements are robust under mixing, while the less abundant ones are more affected. The total global abundances change only slightly from mixing, since each angular ray contributes its most abundant elements. Furthermore, the predicted kilonova light curves show only minor reddening, with differences below the detectability of state-of-the-art telescopes. In general, we do not observe significant effects from mixing in the time span of the r-process. Consequently, mixing only leads to minor variations in abundances and light curves in ray-by-ray simulations.

Explosive burning↗

DOE OSTI2020

Probing the Partonic Degrees of Freedom in High-Multiplicity p–Pb collisions at √ s NN = 5.02 TeV

We investigate the role of partonic degrees of freedom in high-multiplicity p–Pb collisions at √ s NN = 5.02 TeV carried out at the Large Hadron Collider (LHC) by studying the production and collective flow of identified hadrons at intermediate p T via the coalescence of soft partons from the viscous hydrodynamics (VISH 2 + 1) and hard partons from the energy loss model, linear Boltzmann transport (LBT). We find that combining these intermediate p T hadrons with the low p T hadrons from the hydrodynamically expanding fluid and high p T hadrons from the fragmentation of quenched jets, the resulting hydro-dynamics-coalescence-fragmentation model provides a nice description of measured p T spectra and differential elliptic flow v 2 (p T ) of pions, kaons, and protons over the p T range from 0 to 6 GeV. We further demonstrate the necessity of including the quark coalescence contribution to reproduce the experimentally observed approximate number of constituent quark scaling of hadron v 2 at intermediate p T . Furthermore, our results thus indicate the importance of partonic degrees of freedom and also hint at the possible formation of quark-gluon plasma in high-multiplicity p–Pb collisions at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

A Parametric Study of the SASI Comparing General Relativistic and Nonrelativistic Treatments

We present numerical results from a parameter study of the standing accretion shock instability (SASI), investigating the impact of general relativity (GR) on the dynamics. Using GR hydrodynamics with GR gravity, and nonrelativistic (NR) hydrodynamics with Newtonian gravity, in an idealized model setting, we vary the initial radius of the shock, and by varying its mass and radius in concert, the proto-neutron star compactness. We investigate four compactnesses expected in a post-bounce core-collapse supernova (CCSN). We find that GR leads to a longer SASI oscillation period, with ratios between the GR and NR cases as large as 1.29 for the highest-compactness suite. We also find that GR leads to a slower SASI growth rate, with ratios between the GR and NR cases as low as 0.47 for the highest-compactness suite. We discuss implications of our results for CCSN simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

First observation of increased DT yield over prediction due to addition of hydrogen

Here, in a number of reported instances, implosions utilizing fuel mixtures have resulted in anomalously low fusion yields below those predicted by radiation-hydrodynamics simulations. Inter-species ion diffusion has been suggested as a possible cause of the observed yield degradation in fuel mixture implosions. An experimental platform utilizing hydro-equivalent deuterium–tritium (DT), deuterium–tritium–hydrogen (DTH), and deuterium-tritium-helium3 (DT 3 He) capsule implosions was developed to determine whether the inter-species ion diffusion theory may describe the resulting fuel mixture implosion behavior. The implosion experiments were performed at the Omega laser facility. X-ray images and shell areal density diagnostics results show that the hydro-equivalent three capsules (DT, DTH, and DT 3 He) have similar compression behavior. However, nuclear yield deviation was observed from the scaling determined using a fusion yield formula. In the DT 3 He mixture, a reduced yield of a factor of 0.65 ± 0.13 was observed, which is similar to a yield reduction observed in D 3 He mixture by Rygg et al. (i.e., Rygg effect). In contrast, in the DTH mixture, a factor of 1.17 ±0.15 yield increase was observed, which we named the inverse Rygg effect. The yield increase observed in the DTH mixture is consistent with the inter-species ion diffusion theory where lighter H diffuses away from the core and concentrated DT in the core produces higher yield. An inter-species ion diffusion model, the Zimmerman–Paquette–Kagan–Zhdanov model, implemented in a Lagrangian radiation-hydrodynamics fluid code, was also used to analyze the present data, without the need to assume hydrodynamic equivalence of the capsules, but it does not completely explain the DTH or DT 3 He capsules although its effects are in the correct direction. Simulation-based Bayesian inference was used in the latter analysis to quantify the uncertainty in the numerical simulations. The simulation-based analysis resulted in an inferred Rygg-effect yield decrease factor of 0.91 ± 0.02 for the DT 3 He mixture, and an inferred inverse-Rygg yield increase factor of 1.21 ± 0.04 for the DTH mixture, based on simulations ignoring ion diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Coupled aero-hydro-geotech real-time hybrid simulation of offshore wind turbine monopile structures

Real-time hybrid simulation (RTHS) divides a structural system into an analytical and experimental substructure. The former is based on a well-established analytical model while the latter consists of a physical model in the laboratory, for which there is not a well-established analytical model. This paper extends real-time hybrid simulation to monopile-type Offshore Wind Turbines (OWTs) to enable the investigation of their behavior considering the response of pile foundations under operational and more severe conditions. The embedded foundation and surrounding soil of the OWT are modeled physically in a soil box in the laboratory while the remaining parts of the system and loading are modeled analytically. The program OpenFAST, developed by the National Renewable Energy Laboratory (NREL), is linked to the RTHS coordinator to determine the hydrodynamic and aerodynamic loads acting on the OWT, along with modeling the dynamics of the electric power generation equipment and associated controller for the OWT. The RTHS framework along with its initial implementation and validation are described in this paper. RTHSs of a 5 MW OWT subjected to operational and more severe conditions are performed to experimentally validate the framework. The framework offers a realistic approach to investigate the behavior of OWT structures supported on monopiles. Furthermore, this approach accounts for the coupled response of the OWT structure with its foundation, while experimentally capturing the nonlinearities of the soil-foundation interaction in real-time.

17 WIND ENERGY↗

DOE OSTI2025

Collision-Energy Dependence in Heavy-Ion Collisions from Nonlinear QCD Evolution

We explore the effects of including the energy dependence determined from evolution equations within the color glass condensate framework on observables in ultrarelativistic heavy-ion collisions. This amounts to integrating the JIMWLK evolution equations into the impact parameter-dependent glasma model, which is then coupled to viscous relativistic hydrodynamics. This methodology allows for a systematic representation of nuclei at specific Bjorken-x values, which are probed at different center-of-mass energies of the collision and rapidities of final state particles. Comparing the methodology to the conventional impact parameter-dependent glasma model, we find significant effects on multiplicity distributions and particle spectra, especially in smaller collision systems at the highest center-of-mass energies. Our results highlight the importance of incorporating nonlinear QCD evolution in the description of heavy-ion collisions at varying center-of-mass energies, as the precise extraction of transport coefficients will be affected. This Letter establishes a robust framework for understanding the quark gluon plasma and nuclear structure at high energy, integrating small-x physics into the initial conditions of heavy-ion collisions.

Mäntysaari, Heikki↗

DOE OSTI2020

Constraining physical models at gigabar pressures

High-energy-density (HED) experiments in convergent geometry are able to test physical models at pressures beyond hundreds of millions of atmospheres. The measurements from these experiments are generally highly integrated and require unique analysis techniques to procure quantitative infor mation. This work describes a methodology to constrain the physics in convergent HED experiments by adapting the methods common to many other fields of physics. As an example, a mechanical model of an imploding shell is constrained by data from a thin-shelled direct-drive exploding-pusher experiment on the OMEGA Laser System using Bayesian inference, resulting in the reconstruction of the shell dynamics and energy transfer during the implosion. The model is tested by analyzing synthetic data from a 1-D hydrodynamics code and is sampled using a Markov chain Monte Carlo to generate the posterior distributions of the model parameters. As a result, the goal of this work is to demonstrate a general methodology that can be used to draw conclusions from a wide variety of HED experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Constraining Physical Models at Gigabar Pressures

High-energy-density (HED) experiments in convergent geometry are able to test physical models at pressures beyond hundreds of millions of atmospheres. The measurements from these experiments are generally highly integrated and require unique analysis techniques to procure quantitative information. This work describes a methodology to constrain the physics in convergent HED experiments by adapting the methods common to many other fields of physics. As an example, a mechanical model of an imploding shell is constrained by data from a thin-shelled direct-drive exploding-pusher experiment on the OMEGA Laser System using Bayesian inference, resulting in the reconstruction of the shell dynamics and energy transfer during the implosion. The model is tested by analyzing synthetic data from a 1-D hydrodynamics code and is sampled using a Markov chain Monte Carlo to generate the posterior distributions of the model parameters. The goal of this work is to demonstrate a general methodology that can be used to draw conclusions from a wide variety of HED experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Coupling solidification kinetics with phase-behavior computations in hydrodynamic simulations of high-pressure, dynamic-compression processes

In this study, we report a numerical scheme to integrate models for the kinetics of solidification processes together with phase-behavior computations in the context of continuum-scale hydrodynamic simulations. The objective of the phase-behavior computations is to determine the pressure and temperature, given the following three sets of inputs: (1) an appropriate equation of state to describe our system, (2) the phase fraction(s) produced by the kinetic models, (3) and the volume and internal energy obtained by solving the conservation equations that govern the hydrodynamic behavior. The kinetics are assumed to be governed by the Kolmogorov–Johnson–Mehl–Avrami equation, and the nucleation and growth rates that enter into that equation are functions of the pressure and temperature produced by the phase-behavior computations. Our formulation allows for the fluid and solid phases to be at different temperatures (thermal nonequilibrium) and pressures (arising from surface-tension-induced Laplace contributions). The formulation is presented in a fairly general setting that is independent of any particular material, although we demonstrate it in some examples related to high-energy-density science applications where materials are rapidly compressed to pressures exceeding several gigapascals in less than a microsecond. We conclude with a critical evaluation of our approach and provide suggestions for future work to improve the predictive capabilities and generality of the models.

42 ENGINEERING↗