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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Focused wave interactions with floating structures: a blind comparative study

The paper presents results from the Collaborative Computational Project in Wave Structure Interaction (CCP-WSI) Blind Test Series 2. Without prior access to the physical data, participants, with numerical methods ranging from low-fidelity linear models to fully non-linear Navier-Stokes (NS) solvers, simulate the interaction between focused wave events and two separate, taut-moored, floating structures: a hemispherical-bottomed cylinder and a cylinder with a moonpool. The ‘blind’ numerical predictions for heave, surge, pitch and mooring load, are compared against physical measurements. Dynamic time warping is used to quantify the predictive capability of participating methods. In general, NS solvers and hybrid methods give more accurate predictions; however, heave amplitude is predicted reasonably well by all methods; and a WEC-Sim implementation, with CFD-informed viscous terms, demonstrates comparable predictive capability to even the stronger NS solvers. Large variations in the solutions are observed (even among similar methods), highlighting a need for standardisation in the numerical modelling of WSI problems.

42 ENGINEERING↗
DOE OSTI2025

3D Convective Urca Process in a Simmering White Dwarf

Abstract A proposed setting for thermonuclear (Type Ia) supernovae is a white dwarf that has gained mass from a companion to the point of carbon ignition in the core. In the early stages of carbon burning, called the simmering phase, energy released by the reactions in the core drive the formation and growth of a core convection zone. One aspect of this phase is the convective Urca process, a linking of weak nuclear reactions to convection, which may alter the composition and structure of the white dwarf. The convective Urca process is not well understood and requires 3D fluid simulations to properly model the turbulent convection, an inherently 3D process. Because the neutron excess of the fluid both sets and is set by the extent of the convection zone, the realistic steady state can only be determined in simulations with real 3D mixing processes. Additionally, the convection is relatively slow (Mach number less than 0.005) and thus a low Mach number method is needed to model the flow over many convective turnovers. Using the MAESTROeX low Mach number hydrodynamic software, we present the first full-star 3D simulations of the A = 23 convective Urca process, spanning hundreds of convective turnover times. Our findings on the extent of mixing across the Urca shell, the characteristic velocities of the flow, the energy-loss rates due to neutrino emission, and the structure of the convective boundary can be used to inform 1D stellar models that track the longer-timescale evolution.

Boyd, Brendan (ORCID:0000000254199751)↗
DOE OSTI2023

Fluctuating hydrodynamics and the Rayleigh–Plateau instability

The Rayleigh–Plateau instability occurs when surface tension makes a fluid column become unstable to small perturbations. At nanometer scales, thermal fluctuations are comparable to interfacial energy densities. Consequently, at these scales, thermal fluctuations play a significant role in the dynamics of the instability. These microscopic effects have previously been investigated numerically using particle-based simulations, such as molecular dynamics (MD), and stochastic partial differential equation–based hydrodynamic models, such as stochastic lubrication theory. In this paper, we present an incompressible fluctuating hydrodynamics model with a diffuse-interface formulation for binary fluid mixtures designed for the study of stochastic interfacial phenomena. An efficient numerical algorithm is outlined and validated in numerical simulations of stable equilibrium interfaces. We present results from simulations of the Rayleigh–Plateau instability for long cylinders pinching into droplets for Ohnesorge numbers of Oh = 0.5 and 5.0. Both stochastic and perturbed deterministic simulations are analyzed and ensemble results show significant differences in the temporal evolution of the minimum radius near pinching. Short cylinders, with lengths less than their circumference, were also investigated. As previously observed in MD simulations, we find that thermal fluctuations cause these to pinch in cases where a perturbed cylinder would be stable deterministically. Finally, we show that the fluctuating hydrodynamics model can be applied to study a broader range of surface tension–driven phenomena.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Validation of simulated wave energy converter responses to focused waves

Focused waves are increasingly used as a design method to predict extreme loads for offshore structures, such as wave energy converters (WECs). However, the ability of computational methods to accurately simulate the prescribed focused waves and resulting wave structure interactions is not well validated. This study uses two computational fluid dynamic methods, WEC-Sim and STAR-CCM+, to evaluate the responses of two taut-moored WEC-like buoys to three different focused waves. WEC-Sim is a computationally efficient, mid-fidelity model, based on linearised potential flow theory. In contrast, STAR-CCM+ is a high-fidelity, three-dimensional, unsteady, Reynolds-averaged Navier–Stokes-based model. The two buoy geometries considered are a hemispherical-bottom cylinder and a truncated cylinder with a cylindrical moon pool. Three prescribed focused waves are generated in both codes, and the simulated displacements and mooring loads are obtained for each of the hydrodynamic bodies. The resulting STAR-CCM+ generated focused waves have approximately the same accuracy, in comparison with the analytic solution, as experimentally generated focused waves. The WEC-Sim simulated responses, when using viscous coefficients tuned to the measured displacements, are, on average, within 11.5% of experimentally measured values, whereas the STAR-CCM+ simulated responses are, on average, within 14.2% of the experimental values.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Cosmological Hydrodynamics at Exascale: A Trillion-Particle Leap in Capability

Resolving the most fundamental questions in cosmology requires simulations that match the scale, fidelity, and physical complexity demanded by next-generation sky surveys. To achieve the realism needed for this critical scientific partnership, detailed gas dynamics must be treated self-consistently with gravity for end-to-end modeling of structure formation. Exascale computing enables simulations that span survey-scale volumes while incorporating key astrophysical processes that shape complex cosmic structures. We present results from CRK-HACC, a cosmological hydrodynamics code built for extreme scalability. Using separation-of-scale techniques, GPU-resident tree solvers, in situ analysis pipelines, and multi-tiered I/O, CRK-HACCexecuted Frontier-E: a four trillion particle full-sky simulation, over an order of magnitude larger than previous efforts. The run achieved 513.1 PFLOPs peak performance, processing 46.6 billion particles per second and writing more than 100 PB of data in just over one week of runtime. Frontier-E marks a significant advance in predictive modeling for next-generation cosmological science.

Frontiere, Nicholas [Argonne National Laboratory (↗
DOE OSTI2023

NED_Asteroid_Energy_Deposition

In the event of a potentially catastrophic asteroid impact, with sufficient warning time, deploying a nuclear device remains a powerful option for planetary defense if a kinetic impactor proves insufficient. Predicting the effectiveness of a potential nuclear deflection or disruption mission depends on accurate multiphysics simulations of the device's x-ray energy deposition into the asteroid and the resulting material ablation. These simulations span many orders of magnitude, require a variety of different complex physics packages, and are computationally expensive. Having an efficient and accurate way of modeling this system is necessary for exploring a mission's sensitivity to the asteroid's range of physical properties. To expedite future simulations, we present a completed x-ray energy deposition model developed using the radiation-hydrodynamics code Kull which can be used to initiate a nuclear mitigation mission hydrocode. The model spans a wide variety of possible mission initial conditions: four different asteroid-like materials (Silicon Dioxide, Forsterite, Iron, and Ice), two different source spectra (1 and 2 keV blackbodies), and then a broad range of radiation fluences (0.0001 to 1 kt per square meter), source durations (10 to 100 ns), and asteroid porosities (0 to 80 percent). Using blowoff momentum as the primary metric, the modelinitiated simulation results match the full radiation-hydrodynamics results to within 10 percent. Please reference the journal article: Burkey et al., X-Ray Energy Deposition Model for Simulating Asteroid Response to a Nuclear Planetary Defense Mitigation Mission, Planetary Science Journal, (2023) for more information.

Managan, RobertA↗
DOE OSTI2025

Staged Z-Pinch Modeling with Hydra and Chicago Codes (CRADA Final Report)

This was a collaborative effort between Lawrence Livermore National Security, LLC (“Contractor”), as manager and operator of Lawrence Livermore National Laboratory (“LLNL”) and Magneto Inertial Fusion Technologies, Inc. ("MIFTI" or "Participant"), to develop staged Z-pinch (SZP) modeling with HYDRA and CHICAGO codes. The contractor performed numerical simulations of the SZP to reproduce and validate Participant’s MACH2 simulations. Contractor used the state-of-the-art, radiation-magneto-hydrodynamic code HYDRA, as well as the hybrid particle-in-cell code CHICAGO. This work helped assess the physical and numerical uncertainties associated with the simulations and the overall fusion scheme. HYDRA is a well-established and validated code. This work has helped to independently assess the SZP concept to give confidence to the Participant and to the broader fusion community.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

GnarlyX: Eulerian multi-material hydrodynamics coupled to equation of state and hyperelastic, plastic constitutive models

GnarlyX is a new hydrocode for direct numerical simulations of the microstructural behavior of high explosives at the mesoscale. We summarize the computational framework for multidimensional, Eulerian multi-material hydrodynamics coupled to EOS and hyperelastic, plastic constitutive models. We present 1D verification tests of multi-material only and combined multimaterial and strength capabilities with comparisons to exact solutions of shock states resulting from an incident shock impacting the material interface of PBX9502 and aluminum. We show that GnarlyX performs well in capturing the resulting shock waves in examining numerical convergence with exact solutions. In later work, we will summarize the thermomechanics and multi-dimensional, parallel computing capabilities in GnarlyX with multi-dimensional verification tests.

42 ENGINEERING↗
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 OSTI2024

Comparative study of machine learning techniques for post-combustion carbon capture systems

Computational analysis of countercurrent flows in packed absorption columns, often used in solvent-based post-combustion carbon capture systems (CCSs), is challenging. Typically, computational fluid dynamics (CFD) approaches are used to simulate the interactions between a solvent, gas, and column's packing geometry while accounting for the thermodynamics, kinetics, heat, and mass transfer effects of the absorption process. These simulations can then be used explain a column's hydrodynamic characteristics and evaluate its CO 2 -capture efficiency. However, these approaches are computationally expensive, making it difficult to evaluate numerous designs and operating conditions to improve efficiency at industrial scales. In this work, we comprehensively explore the application of statistical ML methods, convolutional neural networks (CNNs), and graph neural networks (GNNs) to aid and accelerate the scale-up and design optimization of solvent-based post-combustion CCSs. We apply these methods to CFD datasets of countercurrent flows in absorption columns with structured packings characterized by several geometric parameters. We train models to use these parameters, inlet velocity conditions, and other model-specific representations of the column to estimate key determinants of CO 2 -capture efficiency without having to simulate additional CFD datasets. We also evaluate the impact of different input types on the accuracy and generalizability of each model. We discuss the strengths and limitations of each approach to further elucidate the role of CNNs, GNNs, and other machine learning approaches for CO 2 -capture property prediction and design optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

On the study of hydrodynamic instabilities in the presence of background magnetic fields in high-energy-density plasmas

Blast-wave-driven hydrodynamic instabilities are studied in the presence of a background B-field through experiments and simulations in the high-energy-density (HED) physics regime. In experiments conducted at the Laboratoire pour l’utilisation des lasers intenses (LULI), a laser-driven shock-tube platform was used to generate a hydrodynamically unstable interface with a prescribed sinusoidal surface perturbation, and short-pulse x-ray radiography was used to characterize the instability growth with and without a 10-T B-field. The LULI experiments were modeled in FLASH using resistive and ideal magnetohydrodynamics (MHD), and comparing the experiments and simulations suggests that the Spitzer model implemented in FLASH is necessary and sufficient for modeling these planar systems. These results suggest insufficient amplification of the seed B-field, due to resistive diffusion, to alter the hydrodynamic behavior. Although the ideal-MHD simulations did not represent the experiments accurately, they suggest that similar HED systems with dynamic plasma-β (=2μ0ρv2/B2) values of less than ∼100 can reduce the growth of blast-wave-driven Rayleigh–Taylor instabilities. These findings validate the resistive-MHD FLASH modeling that is being used to design future experiments for studying B-field effects in HED plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Impact of the numerical solution approach of a plant hydrodynamic model (v0.1) on vegetation dynamics

Abstract. Numerous plant hydrodynamic models have started to be implemented in vegetation dynamics models, reflecting the central role of plant hydraulic traits in driving water, energy, and carbon cycles, as well as plant adaptation to climate change. Different numerical approximations of the governing equations of the hydrodynamic models have been documented, but the numerical accuracy of these models and its subsequent effects on the simulated vegetation function and dynamics have rarely been evaluated. Using different numerical solution methods (including implicit and explicit approaches) and vertical discrete grid resolutions, we evaluated the numerical performance of a plant hydrodynamic module in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES-HYDRO version 0.1) based on single-point and global simulations. Our simulation results showed that when near-surface vertical grid spacing is coarsened (grid size >10 cm), the model significantly overestimates aboveground biomass (AGB) in most of the temperate forest locations and underestimates AGB in the boreal forest locations, as compared to a simulation with finer vertical grid spacing. Grid coarsening has a small effect on AGB in the tropical zones of Asia and South America. In particular, coarse surface grid resolution should not be used when there are large and prolonged water content differences among soil layers at depths due to long dry-season duration and/or well-drained soil or when soil evaporation is a dominant fraction of evapotranspiration. Similarly, coarse surface grid resolution should not be used when there is lithologic discontinuity along the soil depth. This information is useful for uncertainty quantification, sensitivity analysis, or the training of surrogate models to design the simulations when computational cost limits the use of ensemble simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Axisymmetric hydrodynamics in numerical relativity: treating coordinate singularity, artificial heating and modeling MHD instabilities

Two-dimensional axisymmetric simulations of binary neutron star (BNS) merger remnant are a cheap alternative to 3D simulations. To maintain realism for secular timescales, simulations must avoid accumulated errors from drifts in conserved quantities and artificial heating, and they must model turbulent transport in a way that remains plausible throughout the evolution. It is also crucial to avoid numerical artifacts due to the polar coordinate axis singularity. Methods that behave well near the axis often break flux-conservative form of the hydrodynamic equations, resulting in significant drifts in conserved quantities. We present a flux-conservative scheme that maintains smoothness near the axis without sacrificing conservative formulation of the equations or incurring drifts in conserved global quantities. We compare the numerical performance of different treatments of the hydrodynamic equations when evolving a hypermassive neutron star resembling the remnant of a BNS merger. These simulations demonstrate that the new scheme combines the axis smoothness of non-conservative methods with the mass and angular momentum conservation of other conservative methods on $\sim 10^2$ ms timescales of viscous and neutrino-driven evolution. Because fluid profiles remain smooth in the remnant interior, it is possible to remove artificial heating by evolving the entropy density. We show how physical heating and cooling terms can be easily calculated from source terms of the conservative evolution variables and demonstrate our implementation. Finally, we discuss and implement improvements to the effective viscosity scheme to better model the effect of magnetohydrodynamic instabilities as the remnant evolves.

axisymmetry↗
DOE OSTI2025

Demonstration of a code coupling framework for modeling beam-collimator impacts in the advanced photon source

The high-brightness beams being produced in current and future accelerators present new machine protection concerns with the potential for high-energy-density (HED) conditions ( >100 J/mm 3 ) in beam-intercepting components. Simulating HED conditions in accelerators requires utilizing a suite of physics codes for particle dynamics, particle-matter interactions, and hydrodynamics. This paper describes a method of coupling the codes elegant, fluka, and flash to simulate the effects of a rapid beam loss in the advanced photon source storage ring and the resulting interaction of the beam and collimators. This paper expands previous work [J. Dooling et al., Collimator irradiation studies at the advanced photon source, in Proceedings of the IBIC-2023 (2023), pp. 245–249] by introducing a definition of the evolving geometry of the collimator surface as well as providing methods for simulating the absorption of synchrotron radiation and tracking shower particles produced during beam strikes. We demonstrate this framework by simulating machine conditions of the APS ring before and after its recent upgrade. Simulation results are compared with observed damage to collimators and test samples taken from the APS ring.

Accelerator/storage ring control systems↗
DOE OSTI2025

60 years of science in ICF: from conception to scientific breakeven on the National Ignition Facility

The recent achievements of a burning plasma, fusion ignition, and scientific energy gain with deuterium-tritium (DT) fuel at Lawrence Livermore National Laboratory’s National Ignition Facility (NIF) represents a major milestone in the development of inertial confinement fusion (ICF) and all of fusion research. In these experiments, fuel pressures well in excess of hundreds of GBars were achieved in the compressed fuel, and robust alpha heating of the fuel, far in excess of the energy provided by the implosion, were demonstrated for the first time. These achievements occurred 60 years after the inception of ICF and the first laser demonstration, and were made possible by more than five decades of research at laser facilities around the world. Advances in laser technology both in wavelength and precision, motivated by improved understanding of laser-plasma interaction physics and the demands of targets; improvements in target fabrication inspired by the need to control and minimize hydrodynamic instabilities in the implosion; and multi-dimensional simulations and diagnostics have been critical to this achievement. This paper will summarize the scientific and technical advances, the surprises, and the challenges that had to be overcome to achieve these goals.

fusion↗
DOE OSTI2024

Pivotal Roles of Triple Screening-Topological, Electrostatic, and Hydrodynamic-On Dynamics in Semidilute Polyelectrolyte Solutions

For semidilute polyelectrolyte solutions, it is generally assumed that topological, electrostatic, and hydrodynamic interactions are screened (called triple screening). Despite a large body of research focused on polyelectrolyte solutions, the concept of triple screening has never been rigorously verified. In this work, we test the concept by probing concentration fluctuations in aqueous solutions containing a well-studied polyelectrolyte, sodium poly(styrenesulfonate) (NaPSS) with neutron scattering, theory, and molecular dynamics simulations. Neutron spin–echo (NSE) and small-angle neutron scattering (SANS) data from semidilute solutions of NaPSS are presented at different polymer and salt (NaCl) concentrations. A combined theory for structure (J. Chem. Phys. 105, 5183 (1996)) and dynamics (J. Chem. Phys. 107, 2619 (1997)), which captures effects of hydrodynamic, topological, and electrostatic screening, is used to interpret the experimental results. Here, the theory quantitatively predicts the decay rate obtained from the NSE measurements while capturing the shape and concentration dependencies of the polyelectrolyte peak observed in the SANS spectra. Detailed comparisons of the theory and the experiments reveal that the wavevector-dependent decay rate of concentration fluctuations in semidilute solutions of polyelectrolytes is dictated by the screening of hydrodynamic, topological, and electrostatic interactions. This conclusion is corroborated by coarse-grained molecular dynamics simulations, executed without any hydrodynamic interactions, which fail to capture the correct wavevector dependence of the decay rate. These results highlight that the theories based on the concept of triple screening provide a quantitative framework for predicting a relation between the structure and dynamics of polyelectrolyte solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

OpenFAST simulation of floating wind turbines with large heading change

Previous versions of OpenFAST, the physics-based wind turbine engineering and design tool developed by the National Renewable Energy Laboratory, were limited to small rotations of the floating platform. This prevented OpenFAST from being used to simulate important events, such as the loss of a mooring line, or specific floater concepts that might experience large platform yaw motion. To overcome this limitation, we modify the structural dynamics and hydrodynamics modules of OpenFAST to allow unrestricted platform yaw motion. We apply the improved version of OpenFAST to simulate the drifting of a floating wind turbine system after the loss of a mooring line. The results, including both global motion and the internal structural loads at selected locations, appear credible and consistent with expectations.

17 WIND ENERGY↗
DOE OSTI2021

Hydrodynamic description of D meson production in high-energy heavy-ion collisions *

The large values and constituent-quark-number scaling of the elliptic flow of low-p T D mesons imply that charm quarks, initially produced through hard processes, might be partially thermalized through strong interactions with quark-gluon plasma (QGP) in high-energy heavy-ion collisions. To quantify the degree of thermalization of low-p T charm quarks, we compare the D 0 meson spectra and elliptic flow from a hydrodynamic model to experimental data as well as transport model simulations. We use an effective charm chemical potential at the freeze-out temperature to account for the initial charm quark production from hard processes and assume that they are thermalized in the local comoving frame of the medium before freeze-out. D 0 mesons are sampled statistically from the freeze-out hyper-surface of the expanding QGP as described by the event-by-event (3+1)D viscous hydrodynamic model CLVisc. Both the hydrodynamic and transport models can describe the elliptic flow of D 0 mesons at P T < 3 GeV/c as measured in Au+Au collisions at $sqrts_{NN} = 200$ GeV. Though the experimental data on spectra are consistent with the hydrodynamic result at small p T ~ 1 GeV/c, they deviate from the hydrodynamic model at high transverse momentum, GeV/c. The diffusion and parton energy loss mechanisms in the transport model can describe the measured spectra reasonably well within the theoretical uncertainty. Finally, our comparative study indicates that charm quarks only approach local thermal equilibrium at small p T , even though they acquire sizable elliptic flow comparable to light-quark hadrons at both small and intermediate p T .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗