Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamic simulation”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Hourglass control in staggered-grid hydrodynamics using virtual element stabilization techniques

Numerical simulations using the staggered-grid hydrodynamics (SGH) discretization suffer from hourglass instabilities. In this work, we develop a stabilization method to suppress the hourglass instabilities using techniques from the virtual element method (VEM). The stiffness matrix of the VEM consists of two terms: the consistency matrix which is rank deficient and the stability matrix. Here, we first show that in two dimensions and on general polygons, the stiffness matrix of the SGH is identical to the consistency matrix of the linear VEM for both the diffusion equation and the linear elasticity equation. These analyses explain the origin of the hourglass instabilities of the SGH discretization method, and establish a theoretical foundation for our proposed stabilization method by augmenting the stiffness matrix of the SGH discretization using the VEM stability matrix. Then, we present numerical examples using Lagrangian SGH simulations. The numerical experiments demonstrate that the proposed VEM stabilization method is effective at eliminating hourglass modes in the SGH discretization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A lumped particle direct simulation Monte-Carlo method combined with the collisional-radiative model for simulations of non-equilibrium laser-induced plasma plumes

Collisional plasma plumes induced by laser irradiation of material targets exhibit large variations in local density as well as ionization and excitation states, making purely hydrodynamic or kinetic simulations inaccurate or infeasible. To address this challenge and capture non-equilibrium effects in laser-induced plasma plumes at arbitrary degrees of ionization, we develop a hybrid computational approach that combines the kinetic direct simulation Monte Carlo (DSMC) method with a collisional-radiative model (CRM). This ℓDSMC-CRM approach utilizes a lumped particle method to represent minor fractions of excited ions in particle-based simulations and a special coarse-graining technique for atomic spectra and photoionization rates, ensuring numerical convergence at reduced computational cost. The hybrid approach is applied to simulate spatially homogeneous relaxation as well as one- and two-dimensional expansions of plasma plumes induced by irradiation of a copper target by a nanosecond laser pulse in a vacuum or background gas. The comparison with an equilibrium model, where local Saha-Boltzmann equilibrium is enforced, shows that the non-equilibrium effects play a dominant role. The equilibrium model can fail to predict the flow structure and strongly underestimate the degree of absorption of laser radiation by the plume. The ℓDSMC-CRM approach is validated against experimental data demonstrating reasonable agreement with the experimental electron density and temperature, while the equilibrium model is found to dramatically underestimate electron density and temperature. The flexibility of the ℓDSMC-CRM approach allows for its seamless integration into existing DSMC frameworks, making it a valuable tool for high-fidelity plasma modeling in laser-material interactions, laser-based manufacturing, and beyond.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Signatures of QCD conductivities in heavy-ion collisions

Dissipative processes are pivotal for understanding the hydrodynamic evolution of hot and dense quantum chromodynamics (QCD) matter created in relativistic nuclear collisions. The interplay of multiple conserved charges—net baryon, strangeness, and electric charge—is of particular interest. Here, we simulate the longitudinal hydrodynamic evolution with the three diffusion currents in a hydrodynamic model with a lattice-QCD-based equation of state, NEOS -4 D , and estimate rapidity distributions including diffusive corrections to the phase-space distribution in the presence of multiple charges, which ensure charge conservation at particlization. We determine the response of particle yields at midrapidity to changes in the diagonal and off-diagonal conductivities. Inversely, we find that most components of the conductivity matrix can be constrained experimentally using identified particle multiplicities at different collision energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Multiphase computational fluid dynamics modeling of reacting flows in absorption columns for carbon capture

First-principles derived computational fluid dynamics (CFD) simulations have been proposed as a fundamental tool for investigating solvent-based CO 2 absorption in packed columns due to their ability to accurately represent the underlying nonlinear, multiscale dynamics. Numerous studies have previously utilized such CFD simulations to investigate hydrodynamics of columns with structured and random packings by assessing the key hydrodynamic metrics such as the interfacial and wetted areas. While mapping such metrics for different conditions is essential to the optimization of absorption columns, it is not sufficient, as the CO 2 capture rate depends also on the coupled, nonlinear dynamics from the underlying chemical reaction kinetics, thermodynamics, and heat-transfer rates. In this work, we present detailed CFD simulation results augmented by incorporating the effects of interfacial physical mass transfer of CO 2 , heat release from chemical reaction kinetics, and thermophysical property variations from resulting temperature gradients. We demonstrate the applicability of the proposed approach in numerically assessing the performance of packed columns by evaluating key hydrodynamic quantities, CO 2 absorption rates, and temperature rise in a reference column with packings that are structurally similar to the Sulzer Mellapak™ 250.Y packing, for different solvent inflow velocities and temperatures. Predictions from simulation results are found to be consistent with the trends in experimental observations from the literature, suggesting that the predictive capabilities of the simulation framework can be leveraged to guide the future development of absorber-column designs and optimized process flowsheets.

Absorption columns↗
DOE OSTI2025

Exploring Flood Predictability in Taiwan through Coupled Atmospheric–Hydrological and High-Performance Hydrodynamic Models

Effective flood simulation capabilities can tremendously support early warning and disaster prevention. To examine the applicability of a fully physics-based and high-performance flood simulation and forecasting modeling framework for a flood-prone region in Taiwan, we conduct a numerical experiment that couples the Weather Research and Forecasting (WRF) Model, WRF-Hydrological modeling system (WRF-Hydro), and the Two-Dimensional Runoff Inundation Toolkit for Operational Needs (TRITON) to perform integrated rainfall, streamflow, and flood simulations. Furthermore, we first use the coupled WRF and WRF-Hydro (WWH) to predict rainfall and streamflow and then drive TRITON with the predicted streamflow hydrographs to simulate flood depth and inundation area. With the refined spatial resolution and parameterization, this framework can better predict rainfall with reasonable spatial patterns. Although WWH could overestimate the amount of rainfall in some areas, the uncertain rainfall–streamflow predictions produce reasonable flood maps able to pinpoint regions at risk of flooding. In terms of model efficiency, the graphics processing unit–based computation can yield a speed-up factor as high as ∼13 compared to the central processing unit–based computation, promoting the efficacy of the coupled modeling framework in practical real-time flood forecasting.

Coupled models↗
DOE OSTI2024

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Simulation of Binary-single Interactions in AGN Disk. I. Gas-enhanced Binary Orbital Hardening

Stellar-mass binary black hole (BBH) mergers within the accretion disks of active galactic nuclei may contribute to gravitational wave (GW) events detected by ground-based GW detectors. In particular, the interaction between a BBH and a single stellar-mass black hole (sBH), known as the binary-single interaction (BSI) process, can potentially lead to GW events with detectable nonzero eccentricity. Previous studies of the BSI process, which neglected the effects of gas, showed that BSIs contribute non-negligibly to GW events in a coplanar disk environment. In this work, we conduct a series of two-dimensional hydrodynamical and N-body simulations to explore the BSI in a gas environment by coupling REBOUND with Athena++. We perform 360 simulation runs, spanning parameters in disk surface density Σ 0 and impact parameter b. We find that the gas-induced energy dissipation within the three-body system becomes significant if the encounter velocity between the sBHs is sufficiently large (≫c s ). Our simulation results indicate that approximately half of the end states of the BSI are changed by gas. Furthermore, at higher gas density, the number of close encounters during the BSI process will increase, and the end-state BBHs tend to be more compact. Consequently, the presence of gas may shorten the GW merger timescale for end-state BBHs and increase the three-body merger rate.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The Effects of Material-Filled Voids on Detonation Wave Shape in Rubberized RDX Explosives

The sensitivity of explosives is affected by inhomogeneities within the material. This is evident in the increased shock sensitivity of explosives with slightly lower densities resulting from an increased number of hotspots. The influence of hotspots on explosive initiation has been well studied; however, few studies have been conducted on the effect of intermediate-sized voids (0.1–10 mm) on a propagating detonation wave. Cylindrical voids filled with air have been studied for diameters ranging from 0.3 mm to 0.8 mm for both 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) and 1,3,5-trinitro-1,3,5-triazinane (RDX)-based rubberized explosives. Continuing the investigation into single cylindrical voids, this study examined the effects of 0.5 mm diameter voids filled with different inert cylindrical metals on the detonation wave shape for an RDX-based rubberized explosive. The metals selected for experiments were 1066 aluminum, brass, copper, and tungsten. The propagation of the detonation wave was captured using a digital streak camera. Experimental results showed that the extent of detonation wave shaping was closely tied to the density differential between the bulk explosive and metal insert. Forty-four different filler materials, including non-metals, were simulated using a hydrodynamic code to further analyze material inclusion effects. The main factors hypothesized to be of interest were bulk sound speed, shock impedance, and filler material density. We found that the local detonation delay could be correlated fairly well to a ratio of bulk sound speed and density. Understanding the influence of material inclusions on detonation performance and wave shape allows for tailoring of detonations.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2026

Development of a coupled experimental–computational approach for engineering optimization of spout-fluidized bed particle coating systems

The design of spout-fluidized bed (SFB) coating systems for nuclear particle fuels typically relies on trial-and-error processes, comprising iterative and time-consuming coating deposition experiments and post-deposition characterization. At an engineering scale, this approach to guided SFB system design is inefficient, highlighting the need for streamlined experimental methodologies which can correlate fluidization conditions to downstream coating outcomes. In this study, we combine time-resolved particle image velocimetry (PIV) with CFD–DEM simulations to benchmark hydrodynamic behavior in a 3D spout-fluidized bed. By exploiting easily accessible optical measurements of particle motion at the bed wall and within the spouting region, we obtain quantitative velocity fields that can be directly compared with model predictions of the occluded bed region, without resorting to complex imaging and characterization techniques such as X-ray or magnetic resonance tomography. Experimental benchmarking reveals strong agreement between CFD–DEM and PIV in the spout and annulus regions, while discrepancies near the wall highlight areas for future model development. Here, the proposed integrated experimental–numerical framework will enable a direct connection between measured variables and numerically predicted fluidization performance of dense, surrogate nuclear particle fuel feedstock such that experimental SFB component design can be rapidly evaluated, informing design decisions for nozzle geometry and operating conditions. Future work will extend this framework by correlating quantified fluidization metrics across nozzle geometries and operating conditions with the resulting coating morphology, microstructure, and uniformity. Establishing these correlations will enable predictive links between hydrodynamic performance and coating quality, providing a rational, scalable basis for optimizing SFB design prior to coating deposition.

CFD/DEM↗
DOE OSTI2026

Violent mergers can explain the inflated state of some of the fastest stars in the Galaxy

A significant number of hypervelocity stars with velocities between 1500 − 2500 km s −1 have recently been observed. The only plausible explanation so far is that they were produced through thermonuclear supernovae in white dwarf binaries. Since these stars are thought to be surviving donors of Type Ia supernovae, a surprising finding was that these stars are inflated, with radii an order of magnitude higher than expected for Roche-lobe-filling donors. Recent attempts at explaining them have combined 3D hydrodynamical supernova explosion simulations with 1D stellar modelling to explain the impact of supernova shocks on runaway white dwarfs. However, only the hottest and most compact of those runaway stars can so far marginally be reproduced by detailed models of runaways from supernova explosions. In this and a companion paper, we introduce a new AREPO simulation of two massive CO white dwarfs that explode via a violent merger. During the merger, the primary white dwarf ignites when the secondary is on its last orbit and plunging towards the primary. In the corresponding aftermath, the core of the secondary white dwarf of 0.16 M ⊙ remains bound, moving at a velocity of ∼2800 km s −1 . We mapped this object into MESA and show that this runaway star can explain the observations of two hypervelocity stars that were dubbed D6-1 and D6-3 based on their original discovery motivated by the D6 scenario, though the violent merger scenario presented here is somewhat distinct from the D6 scenario.

Astronomy and AstroPhysics↗
DOE OSTI2025

athenapk

The LANL fork of the AthenaPK project, which contains curvlinear coordinates The Code is a set of modifications to AthenaPK, an existing open source 3D hydrodynamics and magnetohydrodynamics simulation software for astrophysics, which enable simulation with uniform curvilinear geometries. Specifically it extends the existing Cartesian coordinates to uniform spherical coordinate and uniform cylindrical coordinate systems

Miller, Jonah [@LANL]↗
DOE OSTI2024

Benchmarking of hydrodynamic plasma waveguides for multi-GeV laser-driven electron acceleration

Hydrodynamic plasma waveguides initiated by optical field ionization have recently become a key component of multi-GeV laser wakefield accelerators. Here, we present the most complete and accurate experimental and simulation-based characterization to date, applicable to current multi-GeV experiments and future 100 GeV-scale laser plasma accelerators. Crucial to the simulations is the correct modeling of intense Bessel beam interaction with meter-scale gas targets, the results of which are used as initial conditions for hydrodynamic simulations. The simulations are in good agreement with our experiments measuring evolving plasma and neutral hydrogen density profiles using two-color short pulse interferometry, enabling realistic determination of the guided mode structure for application to laser-driven plasma accelerator design. Published by the American Physical Society 2024

Physics↗
DOE OSTI2026

Thermodynamically consistent incorporation of the Langmuir adsorption model into compressible fluctuating hydrodynamics

For a gas–solid interfacial system where chemical species undergo reversible adsorption, we develop a mesoscopic stochastic modeling method that simulates both gas-phase hydrodynamics and surface coverage dynamics by coupling the Langmuir adsorption model with compressible fluctuating hydrodynamics. To this end, we derive a thermodynamically consistent mass–energy update scheme that accounts for how the mass and energy variables in the gas and surface subsystems should be updated according to the changes in the number of molecules of each species in each subsystem due to adsorption and desorption events. By performing a stochastic analysis for the ideal Langmuir model and the full hydrodynamic system, we analytically confirm that our mass–energy update scheme captures thermodynamic equilibrium predicted by equilibrium statistical mechanics. We find that an internal energy correction term is needed, which is attributed to the difference in the mean kinetic energy of gas molecules colliding with the surface from that computed from the Maxwell–Boltzmann distribution. By performing an equilibrium simulation study for an ideal gas mixture of CO and Ar, with CO undergoing reversible adsorption, we validate our overall simulation method and implementation.

Adsorption↗
DOE OSTI2024

Strong Coupling of Hydrodynamics and Reactions in Nuclear Statistical Equilibrium for Modeling Convection in Massive Stars

We build on the simplified spectral deferred corrections (SDC) coupling of hydrodynamics and reactions to handle the case of nuclear statistical equilibrium (NSE) and electron/positron captures/decays in the cores of massive stars. Our approach blends a traditional reaction network on the grid with a tabulated NSE state from a very large, ${\mathcal O }(100)$ nuclei network. We demonstrate how to achieve second-order accuracy in the simplified-SDC framework when coupling NSE to hydrodynamics, with the ability to evolve the star on the hydrodynamics time step. We discuss the application of this method to convection in massive stars leading up to core collapse. We also show how to initialize the initial convective state from a 1D model in a self-consistent fashion. All of these developments are done in the publicly available Castro simulation code and the entire simulation methodology is fully GPU-accelerated.

Explosive nucleosynthesis↗
DOE OSTI2025

Experimental platforms for investigating feature-driven jets for HED mix model validation

High-energy-density (HED) systems, such as inertial confinement fusion (ICF), are susceptible to hydrodynamic instabilities that can significantly affect both experimental results and modeling predictions. Isolated features, such as fill tubes or divots in the capsule, can cause material to jet as a result of the compressive shock exciting the Richtmyer–Meshkov instability, and serve as one of the primary degradation mechanisms in ICF yield. Simulations of feature-driven jets and how they mix require extensive experimental validation, particularly for understanding to what degree the initial size and shape of a feature influence jet dynamics, and how much instability feeds through downstream layers. A better understanding of feature-driven jetting can improve our mix modeling capabilities and increase hydrodynamic simulation accuracy. This manuscript describes a series of experimental platforms fielded by Los Alamos National Laboratory as a part of the Mshock Omega 60 and ModCons Omega EP campaigns to explore feature-driven jetting. These platforms are designed to benchmark jet evolution and growth as a function of initial feature size and shape, investigate jet-layer interactions leading to instability feedthrough, and will be used to characterize jet-jet interactions resulting from clusters of features. In conclusion, preliminary results for both platforms are shown. The ModCons experiments are on-going, and a discussion of future work directions is included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Drivers and impacts of sediment deposition in Amazonian floodplains

Abstract The Amazon River carries enormous amounts of sediment from the Andes mountains, much of which is deposited in its floodplains. However, accurate quantification of the sediment sink at fine spatiotemporal scales is still challenging. Here, we present a high-resolution hydrodynamic-sediment model to simulate sediment deposition in a representative Amazon/Solimões floodplain. The process is found to be jointly driven by inundation, suspended sediment concentration in the Amazon River, and floodplain hydrodynamics and only weakly correlated with inundation level. By upscaling the sediment deposition rate (1.33 ± 0.24 kg m −2 yr −1 ), we estimate the trapping of 77.3 ± 13.9 Mt (or 6.1 ± 1%) of the Amazon River sediment by the Amazon/Solimões floodplains every year. Widespread deforestation would reduce the trapping efficiency of the floodplains over time, exacerbating downstream river aggradation. Additionally, we show that the deposition of sediment-associated organic carbon plays a minor role in fueling carbon dioxide and methane emissions in the Amazon.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Reaction-in-flight neutrons as a diagnostic for hydrodynamical mixing in double shell inertial confinement fusion capsules

We examine reaction-in-flight (RIF) neutrons as diagnostics for hydrodynamical mixing of high-Z shell material into the hotspot of double shell capsules that are designed for the National Ignition Facility. In particular, we consider the effects of different levels of mixing of tungsten shell material into the DT gas on RIF spectra. Using a set of 1D simulations from the radiation hydrodynamic code xRAGE to determine the temperature and density profiles of the mixed W-DT gas, we find that increasing the mass of mix systematically reduces the ratio of RIF neutrons to primary 14 MeV neutrons. The shape of the RIF spectrum also changes with mix, with the predicted spectrum softening in energy with increasing mix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Computational Analysis of the Effect of Structured Packing Design on Absorption Column Hydrodynamics for Post-Combustion Carbon Capture Applications

Solvent based post-combustion carbon capture technologies have a potential for reducing carbon emissions from fossil-fuel-fired power plants and industrial sources where CO2 emissions are inherently harder to mitigate, such as steel or cement industries. A prominent technology to achieve this is by retrofitting absorption columns to the existing infrastructure. While these systems have been among the less costly alternatives for carbon capture, they still impose a considerable energy penalty to the operation of power plants or industrial facilities. The optimization of CO2 capture rate in solvent-based absorption process is complex as it depends on several factors including CO2 solubility, solvent reaction kinetics and temperature effects on the solubility, reaction rates, surface tension, and thermophysical properties of the solvent and the flue gas. The overall heat and mass transfer also depends on the hydrodynamics, which in turn, is affected by the packing geometry. In the current work, we systematically quantify the effects the design of the structured packing has on the column hydrodynamics, by performing detailed CFD simulations for different geometrical configurations and operating conditions. We then obtain relationships between the key hydrodynamic metrics, such as liquid holdup, interfacial area, wetted area, and pressure drop to the parameters defining the packing geometries and identify new more effective packing designs for the given operating conditions.

Shah, Yash Girish↗