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 811 records · Page 45

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Wino and Real Minimal Dark Matter Excluded by Fermi Gamma-Ray Observations

We show that minimal, fermionic dark matter (DM) models in the $n$-dimensional representation of the weak force with zero hypercharge that make up 100% of the DM under the standard cosmological history are strongly excluded for $n < 9$. This includes the thermal wino, which we show is ruled out even allowing for DM core sizes up to $\sim$6.7 kpc with the preferred local DM density or $\sim$3.7 kpc in addition to the local DM density being half the preferred value, at less than $0.2$ GeV/cm$^3$. We reach these conclusions through dedicated searches with 14 years of Fermi gamma-ray data in the inner Galaxy between 30 GeV and 2 TeV for the continuum gamma-rays produced in the decays of unstable particles produced in DM annihilation and bound-state formation processes. We consider a variety of Milky Way DM profiles in our analyses, including those motivated by modern hydrodynamic cosmological simulations, and show that all the $n < 9$ minimal DM models are disfavored even under the most conservative assumptions for these density profiles. While wino, quintuplet ($n=5$), and $n = 7$ DM models are strongly disfavored by our analyses under the standard cosmology, we discuss how non-standard cosmological histories or DM sub-fractions could still allow for these particles to be realized in nature, with discovery opportunities at next-generation particle colliders and gamma-ray telescopes.

Cosmology and Nongalactic Astrophysics (astro-ph.C↗

DOE OSTI2022

A three-dimensional pore-scale model for redox flow battery electrode design analysis

A three-dimensional (3-D) pore-scale model has been developed to construct multiscale fibrous electrodes for redox flow batteries (RFB). New designs, such as biporous electrodes modify electrode structures by creating secondary pores on single carbon fiber to reduce internal battery resistance. Existing pore-scale models only resolve electrodes to the single carbon fiber scale and cannot incorporate recent multiscale electrode designs into numerical models. Our new model aims to bridge this gap and provide a tool to rapidly screen new electrode configurations. Two multiscale electrodes, laser-perforated and biporous electrodes, were investigated at varying operation conditions with the proposed 3-D pore-scale models. The laser-perforated electrode exhibits a reduced pressure drop, but follows the same permeability correlation compared to the corresponding pristine electrode. For the biporous electrode, the added specific surface area and faster reaction kinetics from the secondary pores are the most influential factors leading to improved battery efficiency. However, operating the biporous electrode in limiting current density conditions should be avoided due to the decreased mass transfer efficiency and a more significant voltage loss. Finally, we believe that our 3-D pore-scale model can accelerate the flow battery electrode design process and provide new insights into electrode geometry optimizations.

25 ENERGY STORAGE↗

DOE OSTI2024

A review on charged-particle transport modeling for laser direct-drive fusion

Inertial confinement fusion (ICF) with the laser-indirect-drive scheme has recently made a tremendous breakthrough recently after decades of intensive research effort. Taking this success to the next step, the ICF community is coming to a general consensus that laser direct-drive (LDD) fusion might be the viable way for enabling inertial fusion energy (IFE) and high-gain targets for other applications. Designing and understanding LDD fusion targets heavily rely on radiation-hydrodynamic code simulations, in which charged-particle transport plays an essential role in modeling laser-target energy coupling and bootstrap heating of fusion-produced α-particles. To better simulate charged-particle transport in LDD targets, over the past four decades the plasma physics community has advanced transport calculations from simple plasma physics models to sophisticated computations based on first-principles methods. In this review, we give an overview of the current status of charged-particle transport modeling for LDD fusion, including what challenges we still face and the possible paths moving forward to advance transport modeling for ICF simulations. We hope this review will provide a summary of exciting challenges to stimulate young minds to enter the field, facilitate further progress in understanding warm-dense matter physics, and ultimately bridge toward the success of reliable LDD fusion designs for IFE and other high-gain ICF applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Multibody for Everybody (M4E): A Symbolic Dynamics Modeling Tool with Applications in Simulation, Control, and Optimization

Developing the analytical model of a multibody system is often the initial step in control and optimization. The analytical model (equations of motion) describes a system’s time evolution under specified forcing conditions. Although developing these equations is easy for simple systems, this process becomes more complex for systems composed of multiple bodies. Deriving equations of motion for complex multibody systems requires specialized expertise in multibody dynamics, is time-consuming, and is susceptible to error. To address this issue, this paper presents an open-source, easy-to-use, systematic framework to derive symbolic equations of motion in both Python and MATLAB using the joint coordinate formulation. This formulation results in a set of ordinary differential equations that use the minimum set of coordinates needed to model a system. The symbolic representation provides better insight into the influence of design parameters on system performance, facilitates sensitivity analysis and parameter studies, and supports direct implementation of control and optimization routines. The tool enables numerical simulation for specified parameter sets, is modular for straightforward integration with other tools and libraries, and allows incorporation of hydrodynamics, mooring, and other external forces. The result is a reproducible, extensible pipeline for modeling, simulation, and design of complex multibody systems. The proposed tool is versatile and can be applied to domains such as robotics, control, and design. In addition, we integrated external libraries that provide capabilities for modeling offshore systems such as underwater robots and marine energy converters.

16 TIDAL AND WAVE POWER↗

DOE OSTI2022

Stream and ocean hydrodynamics mediate partial migration strategies in an amphidromous Hawaiian goby

Abstract Partial migration strategies, in which some individuals migrate but others do not, are widely observed in populations of migratory animals. Such patterns could arise via variation in migratory behaviors made by individual animals, via genetic variation in migratory predisposition, or simply by variation in migration opportunities mediated by environmental conditions. Here we use spatiotemporal variation in partial migration across populations of an amphidromous Hawaiian goby to test whether stream or ocean conditions favor completing its life cycle entirely within freshwater streams rather than undergoing an oceanic larval migration. Across 35 watersheds, microchemical analysis of otoliths revealed that most adult Awaous stamineus were freshwater residents (62% of n = 316 in 2009, 83% of n = 274 in 2011), but we found considerable variation among watersheds. We then tested the hypothesis that the prevalence of freshwater residency increases with the stability of stream flows and decreases with the availability of dispersal pathways arising from ocean hydrodynamics. We found that streams with low variation of daily discharge were home to a higher incidence of freshwater residents in each survey year. The magnitude of the shift in freshwater residency between survey years was positively associated with predicted interannual variability in the success of larval settlement in streams on each island based on passive drift in ocean currents. We built on these findings by developing a theoretical model of goby life history to further evaluate whether mediation of migration outcomes by stream and ocean hydrodynamics could be sufficient to explain the range of partial migration frequency observed across populations. The model illustrates that the proportion of larvae entering the ocean and differential survival of freshwater‐resident versus ocean‐going larvae are plausible mechanisms for range‐wide shifts in migration strategies. Thus, we propose that hydrologic variation in both ocean and stream environments contributes to spatiotemporal variation in the prevalence of migration phenotypes in A. stamineus . Our empirical and theoretical results suggest that the capacity for partial migration could enhance the persistence of metapopulations of diadromous fish when confronted with variable ocean and stream conditions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Toroidal modeling of runaway avalanche in DIII-D discharges

Abstract A toroidal modeling tool is developed to study the runaway electron (RE) avalanche production process in tokamak plasmas, by coupling the Rosenbluth–Putvinski avalanche model (Rosenbluth and Putvinski 1997 Nucl. Fusion 37 1355) with an n = 0 magneto-hydrodynamic (MHD) solver. Initial value numerical simulations are carried out for two DIII-D discharges with different plasma shapes (one near circular, and the other with high elongation). It is found that, assuming the same level of about 1% seed current level, the Rosenbluth–Putvinski model somewhat underestimates the RE plateau current for the circular-shaped plasma, as compared with that measured in DIII-D experiments. For an elongated, higher current plasma, simulations find strong runaway current avalanche production despite the lack of measured plateau RE current in experiments. A possible reason for this discrepancy is a lack of additional RE dissipation physics in the present two-dimensional model. Systematic scans of the plasma boundary shape, at fixed pre-disruption plasma current, find that the plasma elongation helps to reduce the RE avalanche production, confirming recent results obtained with an analytic model (Fülöp et al 2020 J. Plasma Phys. 86 474860101). The effect of the plasma triangularity (either positive or negative), on the other hand, has a minor effect. On the physics side, the avalanche process involves two competing mechanisms associated with the electric field. On the one hand, a stronger electric field produces a higher instantaneous avalanche growth rate. On the other hand, a fast growing RE current quickly reduces the fraction of the conduction current together with the electric field, and hence a faster dissipation of the poloidal flux. As a final result of these two dynamic processes, the runaway plateau current is not always the largest with the strongest initial electric field. These results lay the foundation for future self-consistent inclusion of the MHD dynamics and the RE amplification process.

Physics↗

DOE OSTI2020

The effect of jet-induced medium excitation on Λ/ k 0 s in jet in Pb-Pb collisions

In this work, we use CoLBT-hydro model for simultaneous event-by-event simulations of jet propagation and hydrodynamics evolution of the bulk medium including jet-induced medium excitation. The final reconstructed jet in heavy-ion collision include not only hadrons from medium response for the deposited energy and momentum of hard partons, but also hadrons from the fragmentation and recombination process of in-medium hard partons. We carry out the study with CoLBT-hydro of medium modification of the Λ/ k 0 s in jet in Pb-Pb collisions at √sNN=2.76TeV. We qualitatively describe the enhancement of the Λ/ k 0 s ratio in jet in the intermediate pT region due to the contribution from jet-induced medium flow.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Quantifying transport and electrocatalytic reaction processes in a gastight rotating cylinder electrode reactor via integration of Computational Fluid Dynamics modeling and experiments

Understanding the complexity of the multiple processes of mass, momentum, charge, and heat transport, and how these affect reaction kinetics at the electrode/electrolyte interface is one of the major challenges in the field of energy and catalysis. The rapid and rational scale-up of electrocatalytic systems to industrial scales require a detailed understanding of nonlinear transport-reaction processes, accessible only through the building of multi-physics models that capture with high fidelity the complexity of real-world devices. The gastight rotating cylinder electrode (RCE) reactor is a promising lab-scale tool that can decouple transport from intrinsic kinetics to generate data for first-principle models useful in the design of industrial, electrochemical reactors. Computational Fluid Dynamics (CFD) studies have previously been used to investigate the bulk flow in RCE reactors for simple corrosion and electroplating processes. However, the quantification of changes in local concentration within the viscous layer where catalysis takes place requires capturing the correct flow conditions inside the hydrodynamic boundary layer near the surface of the electrode. Further, this requires simulations with spatial resolution in the nm and μm scale and temporal resolutions between ms and s scales that are similar to the timescales for reactions on the electrode surface. In this study, experimental electrocatalysis is combined with CFD modeling to elucidate and parameterize the hydrodynamics in a gastight RCE reactor. CFD simulations of the electrochemical ferricyanide reduction reaction under mass transport limited conditions are used to evaluate the validity of the CFD model parameters by comparing calculated dimensionless mass transport descriptors to dimensionless correlations obtained experimentally. Justifications for assumptions and details of the simulation methods used in this study are presented to provide a detailed understanding of the effect that each model parameter has on the ability to accurately simulate electrocatalysis in RCE systems. The simulation methodology reported here is a first step towards the development of multi-scale models for the study of transport dependent electrocatalytic processes, such as the electrochemical transformation of CO 2 to fuels and chemicals.

42 ENGINEERING↗

DOE OSTI2025

High-redshift millennium and astrid galaxies in effective field theory at the field level

Effective field theory (EFT) modeling is expected to be a useful tool in the era of future higher-redshift galaxy surveys such as DESI-II and Spec-S5 due to its robust description of various large-scale structure tracers. However, large values of EFT bias parameters of higher-redshift galaxies could jeopardize the convergence of the perturbative expansion. Here, in this paper we measure the bias parameters and other EFT coefficients from samples of two types of star-forming galaxies in the state-of-the-art MilleniumTNG and astrid hydrodynamical simulations. Our measurements are based on the field-level EFT forward model that allows for precision EFT parameter measurements by virtue of cosmic variance cancellation. Specifically, we consider approximately representative samples of Lyman-break galaxies (LBGs) and Lyman-𝛼 emitters (LAEs) that are consistent with the observed (angular) clustering and number density of these galaxies at 𝑧 = 3. Reproducing the linear biases and number densities observed from existing LAE and LBG data, we find quadratic bias parameters that are roughly consistent with those predicted from the halo model coupled with a simple halo occupation distribution model. We also find nonperturbative velocity contributions (fingers of God) of a similar size for LBGs to the familiar case of luminous red galaxies. However, these contributions are quite small for LAEs despite their large satellite fraction values of up to ∼ 30%. Our results indicate that the effective momentum reach 𝑘 max at 𝑧 = 3 for LAEs (LBGs) will be in the range 0.3−0.6⁢ℎ Mpc −1 (0.2−0.8⁢ℎ Mpc −1 ), suggesting that EFT will perform well for high-redshift galaxy clustering. This work provides the first step toward obtaining realistic simulation-based priors on EFT parameters for LAEs and LBGs.

Sullivan, James M. [Massachusetts Inst. of Technol↗

DOE OSTI2023

Spark Channel Dynamics of Electrostatic Discharges (Rev. 1) [Slides]

To model spark discharges, this work implemented, verified, and validated: 1-D radial conservative Lagrangian hydrodynamic scheme (with heat conduction); Eddington/P1 radiative transfer approximation, time-implicit scheme; RLC circuit solver, time-implicit scheme (adaptable to other circuit types). Novel analytics EOS fit for air up to temperatures of ~150,000 K. Novel "seed" electron concept to initialize finite conductivity in spark without pre-heating the air in the channel, as typically used in spark modeling literature. Reasonable, if conservative (from an engineering viewpoint), agreement with experimental data for indirect (circuit) and direct (interferometry) measurements of spark. Many potential avenues of future work: Two-temperature dryodynamics (requires two-temperature EOS development); 2-D axisymmetric geometry for axial (z) variation, modeling delectric electrodes; Kinetic equations for time-dependent charged species populations.

42 ENGINEERING↗

DOE OSTI2022

A Hydro-particle-mesh Code for Efficient and Rapid Simulations of the Intracluster Medium

We introduce the cosmological HYPER code based on an innovative hydro-particle-mesh (HPM) algorithm for efficient and rapid simulations of gas and dark matter. For the HPM algorithm, we update the approach of Gnedin & Hui to expand the scope of its application from the lower-density intergalactic medium (IGM) to the higher-density intracluster medium (ICM). While the original algorithm tracks only one effective particle species, the updated version separately tracks the gas and dark matter particles, as they do not exactly trace each other on small scales. For the approximate hydrodynamics solver, the pressure term in the gas equations of motion is calculated using robust physical models. In particular, we use a dark matter halo model, ICM pressure profile, and IGM temperature–density relation, all of which can be systematically varied for parameter-space studies. We show that the HYPER simulation results are in good agreement with the halo model expectations for the density, temperature, and pressure radial profiles. Simulated galaxy cluster scaling relations for Sunyaev–Zel’dovich (SZ) and X-ray observables are also in good agreement with mean predictions, with scatter comparable to that found in hydrodynamic simulations. HYPER also produces lightcone catalogs of dark matter halos and full-sky tomographic maps of the lensing convergence, SZ effect, and X-ray emission. These simulation products are useful for testing data analysis pipelines, generating training data for machine learning, understanding selection and systematic effects, and for interpreting astrophysical and cosmological constraints.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Development of an x-ray radiography platform to study laser-direct-drive energy coupling at the National Ignition Facility

A platform has been developed to study laser-direct-drive energy coupling at the National Ignition Facility (NIF) using a plastic sphere target irradiated in a polar-direct-drive geometry to launch a spherically converging shock wave. To diagnose this system evolution, eight NIF laser beams are directed onto a curved Cu foil to generate Heα line emission at a photon energy of 8.4 keV. These x rays are collected by a 100-ps gated x-ray imager in the opposing port to produce temporally gated radiographs. The platform is capable of acquiring images during and after the laser drive launches the shock wave. A backlighter profile is fit to the radiographs, and the resulting transmission images are Abel inverted to infer radial density profiles of the shock front and to track its temporal evolution. Here, the measurements provide experimental shock trajectories and radial density profiles that are compared to 2D radiation-hydrodynamic simulations using cross-beam energy transfer and nonlocal heat-transport models.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Observational Signatures of Planets in Protoplanetary Disks: Temperature Structures in Spiral Arms

High-resolution imaging of protoplanetary disks has unveiled a rich diversity of spiral structure, some of which may arise from disk-planet interaction. Using 3D hydrodynamics with β cooling to a vertically stratified background, as well as radiative-transfer modeling, we investigate the temperature rise in planet-driven spirals. In rapidly cooling disks, the temperature rise is dominated by a contribution from stellar irradiation, 0.3%–3% inside the planet radius but always < 0.5% outside. When cooling time equals or exceeds dynamical time, however, this is overwhelmed by hydrodynamic PdV work, which introduces a ~10%–20% perturbation within a factor of ~2 from the planet’s orbital radius. We devise an empirical fit of the spiral amplitude Δ(T)=(M{sub p}/M{sub th}){sup c}(f{sub PdV}e{sup -t{sub a}{sub r}{sub m}/t{sub c}}+f{sub rad}) to take into account both effects. Where cooling is slow, we find also that temperature perturbations from buoyancy spirals—a strictly 3D, nonisothermal phenomenon—become nearly as strong as those from Lindblad spirals, which are amenable to 2D and isothermal studies. Our findings may help explain observed thermal features in disks like TW Hydrae and CQ Tauri, and underscore that 3D effects have a qualitatively important effect on disk structure.

42 ENGINEERING↗

DOE OSTI2021

Thermal effects mediating the flow induced by laser-induced optical breakdown

Here, a short focused laser pulse can generate a hot plasma, which expands then contracts and can eject a hot jet, the speed and direction of which is sensitive to the details of the plasma kernel. The coupling of thermal and chemical nonequilibrium plasma mechanisms with this hydrodynamic development is assessed with detailed simulations of a two-temperature, three-species plasma model. Time scales for electron recombination, thermal relaxation, and diffusion are compared to that of the plasma expansion to anticipate conditions in which these mechanisms might affect the vorticity generation that leads to the ultimate flow pattern. The effect of these mechanisms are analyzed through comparison with corresponding inert-gas and equilibrium models. Thermal-nonequilibrium effects are found to be weak due to rapid relaxation of the heavy-particle and electronic temperatures. In contrast, chemical equilibration occurs at a rate comparable to the expansion and thereby enhance both it and subsequent hydrodynamic mechanisms as the energy stored in ion formation during the preceding breakdown is released by electron recombination. Thermal conduction, enhanced by high-temperature free electrons, weakens the ejection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Computational Analysis of Different Sparging Systems and their Influence in the Fluid-Dynamic Behavior of Bubble Column Reactors

Bubble column bioreactors are being actively considered for gas fermentation applications, specifically for CO2 utilization, and sugars to fuels conversion. Their main advantages include good mass transfer without any moving parts and low-cost of operation and maintenance. However, the design and scale-up of such reactors is challenging specifically for carbon capture applications where a mixture of gases (e.g. CO2/CO/H2) with variable solubilities is used. The overall performance of scaled-up bioreactors (e.g., mass transfer rate) is largely affected by gas holdup, bubble size distribution (BSD), and multiphase hydrodynamics. We investigate the effect of gas sparger designs on the performance of these large-scale bioreactors using computational fluid dynamics simulations in this work, so as to improve CO2 conversion at scale. The gas distribution systems in bubble column reactors not only determines operational regime, but also affects the evolution of the BSD, which in turn influences interfacial mass transfer and ultimately the efficiency of the gas-liquid exchange process. In addition to the BSD, uniformity in gas sparging affects gas holdup and bubble residence time which constitute important metrics of performance in gas-liquid systems. In this work, we use computational models to simulate high fidelity representations of different sparger designs and their effect on the operation of a bubble column reactor. Four different types of spargers have been selected for the computational study (Fig. 1): ladder, multi-ring, single-ring and toroidal. Their effect on superficial velocity, gas holdup mixing efficiency, and BSD will be evaluated in this work. The model uses a multiphase Eulerian framework similar to [1] and include a composition of mixtures of H2/CO/CO2 gases, common in fermentation applications.

BIOMASS FUELS,MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Stochastic fluctuations in relativistic fluids: Causality, stability, and the information current

We develop a general formalism for introducing stochastic fluctuations around thermodynamic equilibrium which takes into account, for the first time, recent developments in the causality and stability properties of relativistic hydrodynamic theories. The method is valid for any covariantly stable theory of relativistic viscous fluid dynamics derived from a covariant maximum entropy principle. We illustrate the formalism with some applications, showing how it could be used to consistently introduce fluctuations in a model of relativistic heat diffusion and in conformally invariant Israel-Stewart theory in a general hydrodynamic frame. Furthermore, the latter example is used to study the hydrodynamic frame dependence of the symmetric two-point function of fluctuations of the energy-momentum tensor.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

A new effective theory for stochastic relativistic hydrodynamics

Thermal fluctuations are a fundamental feature of dissipative systems that are essential for understanding physics near the expected critical point of QCD and in small systems. When such fluctuations are modeled naively in relativistic systems, strange features can appear such as negative self-correlation functions. We construct an effective theory for nonlinear stochastic relativistic hydrodynamics that ensure a well-posed mathematical formulation. Using Crooks fluctuation theorem, we derive a symmetry of the effective action that incorporates fluctuations through a suitable free energy functional. For divergence type theories, the action can then be fully specified using a single vector generating current. The equations of motion obtained using this procedure are guaranteed to be flux conservative and symmetric hyperbolic when the dynamics is causal. This ensures that these equations are well-posed (for suitable initial data) and are in a form that can easily be simulated, including with Metropolis techniques.

Mullins, Nicki [University of Illinois at Urbana-C↗

DOE OSTI2025

Multiplicity dependence of (multi-)strange hadrons in oxygen-oxygen collisions at $\sqrt{s_{NN}}$ TeV using EPOS4 and AMPT

It is anticipated that the Large Hadron Collider (LHC) will collect data from oxygen-oxygen (O+O) collisions at a center-of-mass energy of $\sqrt{s_{NN}}$ = 7 TeV to explore the effects observed in high multiplicity proton-proton (p+p) and proton-lead (p+pb) collisions that closely related to lead-lead (Pb+Pb) collisions. Furthermore, these effects include azimuthal asymmetries in particle production, as well as variations in the abundances and momentum distributions across different hadron species, which are indicative of collective particle production mechanisms induced by the interactions in the presence of a QGP. The upcoming data on O+O collisions at the LHC are expected to constrain the model parameters and refine our understanding of theoretical models. In this work, the predicted transverse momentum (p T ) spectra, rapidity density distributions (dN/dy), particle yield ratios, and p T -differential ratios of (multi)strange hadrons produced in O+O collisions at $\sqrt{s_{NN}}$ = 7 TeV using AMPT and EPOS4 models are presented. AMPT focuses on preformed hadronic interactions, while EPOS4 incorporates a QGP phase. Stronger radial flow in EPOS4 as compared to AMPT is also observed. AMPT incorporates some flow effects, but the implementation of full hydrodynamic flow in EPOS4 appears to be significantly more effective in reproducing the existing experimental data. Both models predict the final state multiplicity overlap with p+p, p+pb, and Pb+Pb collisions.

Heavy-ion collisions↗