Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics simulations”

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

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

At least 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Pre-equilibrium photons from the early stages of heavy-ion collisions

Abstract We use QCD kinetic theory to compute photon production in the chemically equilibrating Quark-Gluon Plasma created in the early stages of high-energy heavy-ion collisions. We do a detailed comparison of pre-equilibrium photon rates to the thermal photon production. We show that the photon spectrum radiated from a hydrodynamic attractor evolution satisfies a simple scaling form in terms of the specific shear viscosityη/sand entropy densitydS/dζ∼ (Tτ 1/3 ) 3/2 ∞ . We confirm the analytical predictions with numerical kinetic theory simulations. We use the extracted scaling function to compute the pre-equilibrium photon contribution in$$ \sqrt{s_{NN}} $$ s NN = 2.76 TeV 0–20% PbPb collisions. We demonstrate that our matching procedure allows for a smooth switching from pre-equilibrium kinetic to thermal hydrodynamic photon production. Finally, our publicly available implementation can be straightforwardly added to existing heavy ion models.

Physics↗
DOE OSTI2026

Plasma hydrodynamics from mean force kinetic theory

Mean force kinetic theory is used to evaluate the electrical conductivity, thermal conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity of a two-component plasma. The results are compared with molecular dynamics simulations. These simulations are made possible by assuming a repulsive Coulomb force for all interactions. Good agreement is found for all coefficients up to a Coulomb coupling strength of Γ ≈ 20. This is over 100-times larger than the coupling strength at which traditional theories break down. It is concluded that mean force kinetic theory is a promising framework for extending hydrodynamics to dense plasmas.

Electrical conductivity↗
DOE OSTI2021

Electroconvection near an ion-selective surface with Butler–Volmer kinetics

We study the effects of interfacial kinetics on the electro-hydrodynamics of ion transport near an ion-selective surface using a combination of linear stability analysis and numerical simulation. The finite kinetics of the electrolyte–electrode interface affects the ion transfer and electroconvection in many ways. On a surface of fixed topography, such as a metal surface of slow and stable ion deposition or covered by a polymer membrane, the finite kinetics reduces the current in one-dimensional ion diffusion/migration, increases the critical voltage for the onset of the electroconvective instability, changes the dynamics of the electroconvection and the overlimiting current, and enhances the lateral ion diffusion within the interfacial layer. The first three effects are indirectly caused by the reaction kinetics and can be characterized by an effective voltage difference across the liquid electrolyte. In comparison, the last effect is controlled by a direct interplay between kinetics and nonlinear electroconvection. Scaling laws for ion transport and features of electroconvection are proposed. We also analyse the linear stability of a surface which evolves under ion deposition and find that the finite kinetics decreases the growth rate of both electroconvective and morphological instabilities and therefore modifies the wavenumber of the most unstable mode.

Mechanics↗
DOE OSTI2025

Interpreting dynamic-compression experiments to uncover the time dependence of freezing: Application to gallium

Using pulsed-power magnetic field sources to compress gallium to gigapascal pressures on nanosecond timescales, we report here experiments on shockless dynamic compression of a liquid metal. Time-resolved velocimetry data reveal signatures of rapid freezing from a metastable liquid state, and we demonstrate that the kinetics of this nonequilibrium solidification can be accurately simulated with a computational modeling framework we have developed in previous studies, where classical nucleation theory is coupled with hydrodynamics. Notably, velocity traces in some of our experiments show evidence of a phase transition, while others do not, even though other types of evidence suggest that solidification may be occurring in all of them. We explain how predictions made by our models regarding the presence or absence of these phase-transition signatures motivated additional experiments that later confirmed the theoretical predictions. Our analysis shows that due to the rapid, quasi-isentropic nature of the loading path, our experiments were able to compress liquid gallium to metastable states that are undercooled below the equilibrium melt temperature by more than 300 K and exhibit pressures that approach five times the equilibrium melt pressure. The understanding gained in this study should form the basis for future dynamic-compression experiments aimed at interrogating melt curves at high pressures.

gallium↗
DOE OSTI2020

Collinear dual-pulse laser optical breakdown and energy deposition

A non-equilibrium model for laser-generated plasmas is used to represent collinear dual nano-second-pulse interactions. In the specific case considered, the breakdown is initiated with an ultraviolet (UV) laser pulse at 266 nm, which pre-ionizes the gas, and it is followed by a near-infrared (NIR) pulse at 1064 nm, which deposits significant energy into the ionized mixture. The model is validated against corresponding experiments, and simulation results are interrogated to understand key features of the plasma-kernel dynamics and the post-discharge hydrodynamics. The hydrodynamics of the non-equilibrium plasma is governed by the two-temperature Navier–Stokes equations accounting for both multiphoton ionization and inverse bremsstrahlung. The interaction between the laser beam and the plasma is modeled based on the radiative transfer equation. The temporal and spatial offsets of the two pulses can generate ionization kernels with different topology and dynamics. It is shown that the UV pre-ionization pulse can tailor the plasma region, leading to a larger ionized volume than would occur for a single-pulse breakdown, and it increases the efficiency of the energy deposition for the following NIR discharge. Vorticity in the early post-discharge phase is generated via baroclinic torque from the misalignment of the radial gradient of density (sudden gas expansion) and the strong pressure gradient (initiated by energy deposition from the laser beam). It is predicted that the collinear dual pulse can be designed to specify the plasma kernel size, shape, and maximum temperature in the breakdown phase, and the initial post-breakdown vorticity and plasma-core decay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

An Assessment of the Influence of Uncertainty in Temporally Evolving Streamflow Forecasts on Riverine Inundation Modeling

Continental-scale river forecasting platforms forecast streamflow at reaches that can be used as boundary conditions to drive a local-scale flood inundation model. Uncertainty accumulated during this process stems not only from any part of the forecasting chain but can also be caused by the daily variations in weather forcing that keeps evolving as the event advances. This work aims to examine the influence of the evolving forecast streamflow on predicting the maximum inundation for extreme floods. A diagnostic case study was made on the basis of a hindcast of Hurricane Matthew striking the eastern U.S. in 2016. The U.S. National Water Model was one-way coupled to a hydrodynamic inundation model through a developed automated workflow. Although the river forcing has significantly mismatched hydrographs versus observations, the simulated peak water surface elevations and maximum extents were validated to be comparable with the observations, which indicates that the inundation model may not be sensitive to the inherited uncertainty from the weather forcing. Moreover, the uncertainty of the forecast streamflow time series caused only one order of magnitude fewer variations in inundation prediction; this dampening effect may become clearer for extreme events with large areas inundated. In addition, the forecast total volume of stream discharge appears to be an important metric for assessing the performance of river forcing for inundation mapping, as a linear correlation between the total volume and the accuracy of the predicted peak water surface elevation and maximum extent was found, with the coefficients of determination all above 0.8. Extra best-practice experience of running similar operational tasks demonstrated the tradeoff between the cost and accuracy gain.

flood, national water model, RIFT↗
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 OSTI2024

Performance characterization and modeling of an oscillating surge wave energy converter

Abstract Testing wave energy converters in the ocean could be expensive and complex, which necessitates the use of numerical modeling. However, accurately modeling the response of wave energy converters with high-fidelity simulations can be computationally intensive in the design stage where different configurations must be considered. Reduced-order models based on simplified equations of motion can be very useful in the design, optimization, or control of wave energy converters. Given the complex dynamics of wave energy converters, accurate representation, and evaluation of relative contributions by different forces are required. This effort is concerned with a performance characterization of the hydrodynamic response of an oscillating surge wave energy converter that is based on a reduced-order model. A state-space model is used to represent the radiation damping term. Morison’s representation of unsteady forces is used to account for the nonlinear damping. Wave tank tests are performed to validate simulations. A free response simulation is used to determine the coefficients of the state-space model. Torque-forced simulations are used to identify the coefficients of the nonlinear damping term for different amplitudes and wave frequencies. The impact of varying these coefficients on the response is investigated. An assessment of the capability of the model in predicting the hydrodynamic response under irregular forcing is performed. The results show that the maximum error is 3% when compared with high-fidelity simulations. It is determined that the nonlinear damping is proportional to the torque amplitude and its effects are more pronounced as the amplitude of the flap oscillations increases.

16 TIDAL AND WAVE POWER↗
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 OSTI2020

Simulation of coupled multiphase flow and geomechanics in porous media with embedded discrete fractures

In fractured natural formations, the equations governing fluid flow and geomechanics are strongly coupled. Hydrodynamical properties depend on the mechanical configuration, and they are therefore difficult to accurately resolve using uncoupled methods. In recent years, significant research has focused on discretization strategies for these coupled systems, particularly in the presence of complicated fracture network geometries. In this work, we explore a finite-volume discretization for the multiphase flow equations coupled with a finite-element scheme for the mechanical equations. Fractures are treated as lower dimensional surfaces embedded in a background grid. Interactions are captured using the embedded discrete fracture model (EDFM) and the embedded finite element method (EFEM) for the flow and the mechanics, respectively. This nonconforming approach significantly alleviates meshing challenges. EDFM considers fractures as lower dimension finite volumes that exchange fluxes with the rock matrix cells. The EFEM method provides, instead, a local enrichment of the finite-element space inside each matrix cell cut by a fracture element. Both the use of piecewise constant and piecewise linear enrichments are investigated. They are also compared to an extended finite element approach. One key advantage of EFEM is the element-based nature of the enrichment, which reduces the geometric complexity of the implementation and leads to linear systems with advantageous properties. Synthetic numerical tests are presented to study the convergence and accuracy of the proposed method. It is also applied to a realistic scenario, involving a heterogeneous reservoir with a complex fracture distribution, to demonstrate its relevance for field applications.

58 GEOSCIENCES↗
DOE OSTI2022

Modeling pre-Exascale AMR Parallel I/O Workloads via Proxy Applications

The present work investigates the modeling of preexascale input/output (I/O) workloads of Adaptive Mesh Refinement (AMR) simulations through a simple proxy application. We collect data from the AMReX Castro framework running on the Summit supercomputer for a wide range of scales and mesh partitions for the hydrodynamic Sedov case as a baseline to provide sufficient coverage to the formulated proxy model. The non-linear analysis data production rates are quantified as a function of a set of input parameters such as output frequency, grid size, number of levels, and the Courant-Friedrichs-Lewy (CFL) condition number for each rank, mesh level and simulation time step. Linear regression is then applied to formulate a simple analytical model which allows to translate AMReX inputs into MACSio proxy I/O application parameters, resulting in a simple “kernel” approximation for data production at each time step. Results show that MACSio can simulate actual AMReX nonlinear “static” I/O workloads to a certain degree of confidence on the Summit supercomputer using the present methodology. The goal is to provide an initial level of understanding of AMR I/O workloads via lightweight proxy applications models to facilitate autotune data management strategies in anticipation of exascale systems.

Godoy, William↗
DOE OSTI2024

Numerical simulations of laser-driven experiments of ion acceleration in stochastic magnetic fields

We present numerical simulations used to interpret laser-driven plasma experiments at the GSI Helmholtz Centre for Heavy Ion Research. The mechanisms by which non-thermal particles are accelerated in astrophysical environments, e.g., the solar wind, supernova remnants, and gamma ray bursts, is a topic of intense study. When shocks are present, the primary acceleration mechanism is believed to be first-order Fermi, which accelerates particles as they cross a shock. Second-order Fermi acceleration can also contribute, utilizing magnetic mirrors for particle energization. Despite this mechanism being less efficient, the ubiquity of magnetized turbulence in the universe necessitates its consideration. Another acceleration mechanism is the lower-hybrid drift instability, arising from gradients of both density and magnetic field, which produce lower-hybrid waves with an electric field that energizes particles as they cross these waves. With the combination of high-powered laser systems and particle accelerators, it is possible to study the mechanisms behind cosmic-ray acceleration in the laboratory. In this work, we combine experimental results and high-fidelity three-dimensional simulations to estimate the efficiency of ion acceleration in a weakly magnetized interaction region. We validate the FLASH magneto-hydrodynamic code with experimental results and use OSIRIS particle-in-cell code to verify the initial formation of the interaction region, showing good agreement between codes and experimental results. We find that the plasma conditions in the experiment are conducive to the lower-hybrid drift instability, yielding an increase in energy ΔE of ~ 264 keV for 242 MeV calcium ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Fusion burn-propagation simulations using the collisional and radiative particle-in-cell code TRIFORCE

The ability to accurately model burn propagation in inertial confinement fusion plasmas is crucial for advancing fusion energy research. This work presents enhancements to the triforce hybrid fluid-kinetic multiphysics code, focusing on its kinetic half, which employs the particle-in-cell (PIC) method with Monte Carlo collisions (MCC). We use a moment-preserving collision model that mitigates numerical noise, particularly in spherical geometries where particle weights vary significantly. Additionally, we refine the treatment of inverse bremsstrahlung to account for electron–ion collision frequency reductions in degenerate plasmas and incorporate a blackbody radiation source to enable realistic photon injection. These improvements enable the simulation of 1-dimensional (1D) spherical fusion burn propagation in deuterium–tritium plasmas. Benchmark comparisons with the hydra radiation-hydrodynamics code confirm that triforce accurately captures the dynamics of hot-spot expansion and burn propagation, demonstrating sensitivity to ignition thresholds consistent with theoretical models. Findings show the ignition cliff to be less steep in our work compared to radiation-hydrodynamic modeling. These results highlight the role of kinetic effects in fusion ignition physics and underscore the necessity of hybrid fluid-kinetic models for advancing predictive capabilities in high-energy-density plasma systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Loads on a Jacket-Supported Wind Turbine During Hurricane Sandy Simulation

Abstract The dynamic response of a jacket-supported offshore wind turbine under coupled wind, wave, and current fields during Hurricane Sandy is the subject of this study. To illustrate the detailed procedure related to the response evaluation of a 5-MW offshore wind turbine with a jacket support structure, we consider a single site where the water depth is 50 m. Loads are computed using two simulation tools, fast and abaqus, with partial coupling. Aerodynamic loads on the turbine rotor are first evaluated using a wind turbine model in fast with a fixed base; then, these rotor aerodynamic loads are applied as point loads at the top of a model of a tower that is supported by a jacket structure in abaqus. On the basis of stochastic simulations, we discuss the applicability of the abaqus jacket modeling and describe the characteristics and comparisons of the aerodynamic and hydrodynamic effects on loads on the jacket members. Details related to the structural model and soil–pile interaction model employed in the analyses are also discussed.

Engineering↗
DOE OSTI2020

Computational Fluid Dynamics Analysis of a 12 MW Circulating Fluidized Bed Rise

A three-dimensional (3D) two fluid model (TFM) model of the 12 MW circulating fluidized bed (CFB) riser at Chalmers University was developed using the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite. The hydrodynamic behavior of the riser was validated by comparing axial pressure profiles from experimental measurements with the results from the simulations. Under some of the operating conditions where there was no external circulation of solids, the pressure profiles from the simulation matched well with the experimental measurements. For the cases where there was external circulation of solids, the pressure profiles with the monodisperse solid phase modeling underpredicted the solids concentration in the freeboard region. With the inclusion of polydispersity in the model, the predicted pressure profiles agreed well with the experimental measurements. A two-dimensional (2D) model of the riser was used to perform simulations of coal combustion. The predicted flue gas outlet compositions of oxygen and carbon dioxide compared reasonably to experiments. The experiments also reported some concentration of carbon monoxide which was overpredicted in the simulations. Additional experimental data such as axial oxygen concentration profile and pressure profiles were also compared with the experiment. It was observed that the oxygen concentration in the bottom bed from the simulations was higher than in the experiments, which suggested that most of the oxidation reactions occurred in the bottom bed. The concentration of volatiles and char in the simulations was lower in the bottom bed than in the experiments, which suggest that there is not sufficient penetration of the solid phase species into the bed.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2020

The Multiphysics on Advanced Platforms Project

In 2015, the Lawrence Livermore National Laboratory started development of next-generation multiphysics simulation capabilities for the National Nuclear Security Administration under the Advanced Technologies Development and Mitigation (ATDM) element of the Advanced Simulation and Computing program in collaboration with the Exascale Computing Project (ECP). A key driver for this effort across the NNSA tri-lab was the emergence of advanced high performance computing (HPC) architectures based on heterogeneous compute capabilities, including GPU based systems, as part of the national drive toward exascale computing platforms at multiple Department of Energy (DOE) facilities. Developing a multiphysics code capable of meeting the various simulation needs of the NNSA as defined by the current generation of integrated codes (or ICs), initially developed as part of the Accelerated Strategic Computing Initiative (ASCI) program beginning in 1996, and able to scale to the current 100 petaflop class pre-exascale systems, as well the forthcoming exaflop class computers, is a daunting challenge. To accomplish this ambitious goal, LLNL has embraced two key themes: use of high-order numerical methods and a modular approach to code development. The LLNL next generation effort is organized under the Multi-Physics on Advanced Platforms Project (MAPP). A foundational component of MAPP is the Axom computer science (CS) toolkit which provides infrastructure for the development of modular, performance portable, multi-physics application codes. MARBL is a next-generation application code built on the Axom base to address the modeling needs of the high energy density physics (HEDP) community for simulating high-explosive, magnetic or laser driven experiments such as inertial confinement fusion (ICF), pulsed-power magneto-hydrodynamics (MHD), equation of state (EOS) and material strength studies as part of the NNSA’s stockpile stewardship program (SSP).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

ELG×LRG Distribution through Dark Matter Halo Dynamics

We investigate the clustering and halo occupation distribution (HOD) of DESI Y1 emission-line (ELGs) and luminous red (LRGs) galaxies at 0.8 < z < 1.1, including their cross-correlation (ELG×LRG), using the A BACUS S UMMIT suite and a new Halo Occupation Model (H OME ) for galaxy multitracers. This integrates intrahalo dynamics, halo exclusion, and quenching, bridging insights from hydrodynamical, HOD, abundance-matching, and semianalytic studies. Leveraging full phase-space information from the Uchuu N-body simulation, and sampling satellites from dark-matter particle positions via physically motivated prescriptions, Home reproduces the anisotropic clustering down to s = 200 h −1 kpc with unprecedented accuracy. Model parameters are inferred solely from two-point statistics using a two-level Bayesian framework, yielding high-fidelity ELG, LRG, and cross-reference catalogs. We find that satellite ELGs behave as incoherent flows within their parent halos, dominating the clustering below 4 h −1 Mpc. The HOD from the best-fit Home has the following properties: (i) 90.50% (85.91%) of ELGs (LRGs) are central galaxies without satellites, residing in halos of M vir ∼ 6.6 × 10 11 (1.2 × 10 13 ) h −1 M ⊙ ; (ii) the ELG×LRG cross-correlation is governed by central-central pairs and shaped by halo exclusion on 2–5 h −1 Mpc scales; (iii) 9.50% (14.09%) of ELGs (LRGs) are satellites, of which 1.09% (3.52%) inhabit halos with a central galaxy of the same species in a maximally conformal configuration, 7.02% (0.005%) orbit complementary hosts in a minimally conformal state, and 0.58% (10.57%) are orphans. The high sensitivity of Home precisely captures the dynamics of satellites in different host environments, opening a promising avenue for understanding systematics and the dynamical nature of dark matter, potentially distinguishing gravity models.

Favole, Ginevra [Universidad de La Laguna (Spain);↗
DOE OSTI2026

Second-order wave excitation forces in WEC-Sim/MOST: Implementation, experimental validation, and code-to-code comparison

Accurate prediction of second-order hydrodynamic loads is essential for floating bodies, including floating offshore wind turbines, wave energy converters, and hybrid wind–wave platforms. These nonlinear effects, arising from both sum- and difference-frequency forcing, are critical for capturing key response characteristics but remain challenging to model efficiently. In this work, we extend the open-source Wave Energy Converter Simulator / MATLAB for Offshore Simulation Tool by implementing second-order wave excitation forces, supporting both the full Quadratic Transfer Function formulation and the Newman approximation. The full Quadratic Transfer Function method is used for all code-to-code comparisons and experimental validation, while the Newman approximation is provided as a computationally lighter alternative. To benchmark the new capability, we perform a code-to-code comparison with OpenFAST and OrcaFlex. We then validate the enhanced model using wave-tank measurements of a 1:96 scale DeepCwind semi-submersible, showing that second-order effects are required to reproduce platform motions. The implementation employs a computationally efficient pre-computation strategy for second-order wave excitation forces, reducing simulation cost while maintaining engineering accuracy. Overall, this work advances the tool as an open-source and versatile tool for modelling floating offshore renewable-energy systems requiring second-order hydrodynamic fidelity.

17 WIND ENERGY↗