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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Critical lensing and kurtosis near a critical point in the QCD phase diagram in and out of equilibrium

In this work, we study the lensing effect of the QCD critical point on hydrodynamic trajectories and its consequences on the net-proton kurtosis $\kappa$ 4 . Including critical behavior by means of the BEST Collaboration equation of state (EoS), we first consider a scenario in equilibrium, then compare with hydrodynamic 0+1⁢D simulations with Bjorken expansion, including both shear and bulk viscous terms. Here, we find that, both in and out of equilibrium, the size and shape of the critical region directly affect if the signal will survive through the dynamical evolution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Future climate doubles the risk of hydraulic failure in a wet tropical forest

Summary Future climate presents conflicting implications for forest biomass. We evaluate how plant hydraulic traits, elevated CO 2 levels, warming, and changes in precipitation affect forest primary productivity, evapotranspiration, and the risk of hydraulic failure. We used a dynamic vegetation model with plant hydrodynamics (FATES‐HYDRO) to simulate the stand‐level responses to future climate changes in a wet tropical forest in Barro Colorado Island, Panama. We calibrated the model by selecting plant trait assemblages that performed well against observations. These assemblages were run with temperature and precipitation changes for two greenhouse gas emission scenarios (2086–2100: SSP2‐45, SSP5‐85) and two CO 2 levels (contemporary, anticipated). The risk of hydraulic failure is projected to increase from a contemporary rate of 5.7% to 10.1–11.3% under future climate scenarios, and, crucially, elevated CO 2 provided only slight amelioration. By contrast, elevated CO 2 mitigated GPP reductions. We attribute a greater variation in hydraulic failure risk to trait assemblages than to either CO 2 or climate. Our results project forests with both faster growth (through productivity increases) and higher mortality rates (through increasing rates of hydraulic failure) in the neo‐tropics accompanied by certain trait plant assemblages becoming nonviable.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Efficient Probabilistic Prediction and Uncertainty Quantification of Tropical Cyclone–Driven Storm Tides and Inundation

Abstract This study proposes and assesses a methodology to obtain high-quality probabilistic predictions and uncertainty information of near-landfall tropical cyclone–driven (TC-driven) storm tide and inundation with limited time and resources. Forecasts of TC track, intensity, and size are perturbed according to quasi-random Korobov sequences of historical forecast errors with assumed Gaussian and uniform statistical distributions. These perturbations are run in an ensemble of hydrodynamic storm tide model simulations. The resulting set of maximum water surface elevations are dimensionality reduced using Karhunen–Loève expansions and then used as a training set to develop a polynomial chaos (PC) surrogate model from which global sensitivities and probabilistic predictions can be extracted. The maximum water surface elevation is extrapolated over dry points incorporating energy head loss with distance to properly train the surrogate for predicting inundation. We find that the surrogate constructed with third-order PCs using elastic net penalized regression with leave-one-out cross validation provides the most robust fit across training and test sets. Probabilistic predictions of maximum water surface elevation and inundation area by the surrogate model at 48-h lead time for three past U.S. landfalling hurricanes (Irma in 2017, Florence in 2018, and Laura in 2020) are found to be reliable when compared to best track hindcast simulation results, even when trained with as few as 19 samples. The maximum water surface elevation is most sensitive to perpendicular track-offset errors for all three storms. Laura is also highly sensitive to storm size and has the least reliable prediction. Significance Statement The purpose of this study is to develop and evaluate a methodology that can be used to provide high-quality probabilistic predictions of hurricane-induced storm tide and inundation with limited time and resources. This is important for emergency management purposes during or after the landfall of hurricanes. Our results show that sampling forecast errors using quasi-random sequences combined with machine learning techniques that fit polynomial functions to the data are well suited to this task. The polynomial functions also have the benefit of producing exact sensitivity indices of storm tide and inundation to the forecasted hurricane properties such as path, intensity, and size, which can be used for uncertainty estimation. The code implementing the presented methodology is publicly available on GitHub.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Interplay between gain and loss in arrays of nonlinear plasmonic nanoparticles: toward parametric downconversion and amplification

With the help of a theoretical model and finite-difference time-domain (FDTD) simulations based on the hydrodynamic–Maxwell model, we examine the effect of difference-frequency generation (DFG) in an array of L-shaped metal nanoparticles (MNPs) characterized by intrinsic plasmonic nonlinearity. The outcomes of the calculations reveal the spectral interplay between gain and loss in the vicinity of the fundamental frequency of the localized surface plasmon resonances. Subsequently, we identify different array thicknesses and pumping regimes facilitating parametric amplification and spontaneous parametric downconversion. Our results suggest that the parametric amplification regime becomes feasible on a scale of hundreds of nanometers and spontaneous parametric downconversion on the scale of tens of nanometers, opening up new exciting opportunities for developing building blocks of photonic metasurfaces.

Shah, Syed A.↗
DOE OSTI2021

Nyx: A Massively Parallel AMR Code for Computational Cosmology

Nyx is a highly parallel, adaptive mesh, finite-volume N-body compressible hydrodynamics solver for cosmological simulations. It has been used to simulate different cosmological scenarios with a recent focus on the intergalactic medium and Lyman alpha forest. Together, Nyx, the compressible astrophysical simulation code, Castro, and the low Mach number code MAESTROeX, make up the AMReX-Astrophysics Suite of open-source, adaptive mesh, performance-portable astrophysical simulation codes. Other examples of cosmological simulation research codes include Enzo, Enzo-P/Cello, RAMSES, ART, FLASH, Cholla, as well as Gadget, Gasoline, Arepo, Gizmo, and SWIFT.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

San Juan Islands Tidal Energy Characterization

The San Juan Islands are an archipelago with multiple tidal channels that produce strong tidal currents that could be harnessed for electricity generation. At present, there is limited electrical generation on the islands, with power primarily provided by a subsea cable connection with the mainland. As part of the C-MIST program, NOAA's CO-OPS program collected ADCP data at 46 stations in the San Juan Islands and adjacent waterways, from April to August of 2017. PNNL has recently developed a 3-D tidal hydrodynamic model for tidal energy resource characterization and assessment in the Salish Sea (Yang et al. 2021). Several tidal channels in the San Juan Islands were identified as top hotspots for potential tidal energy extraction. However, due to the complex geometry and inter-connected waterways, tidal currents around San Juan Islands also exhibit strong spatial and temporal asymmetry. Therefore, it is important to analyze the NOAA ADCP data to define metrics relevant to tidal energy, and further refine PNNL's tidal hydrodynamic model to accurately simulate currents in small tidal channels and account for sharp bathymetry gradients. The outcome of this technical assistance will allow developers to identify the promising opportunities for tidal energy in the San Juan Islands.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Physics-driven Explosions of Stripped High-mass Stars: Synthetic Light Curves and Spectra of Stripped-envelope Supernovae with Broad Light Curves

Stripped-envelope supernovae (SESNe) represent a significant fraction of core-collapse supernovae, arising from massive stars that have shed their hydrogen and, in some cases, helium envelopes. The origins and explosion mechanisms of SESNe remain a topic of active investigation. In this work, we employ radiative-transfer simulations to model the light curves and spectra of a set of explosions of single, solar-metallicity, massive Wolf–Rayet stars with ejecta masses ranging from 4 to 11 M ⊙ , which were computed from a turbulence-aided and neutrino-driven explosion mechanism. We analyze these synthetic observables to explore the impact of varying ejecta mass and helium content on observable features. We find that the light curve shape of these progenitors with high ejecta masses is consistent with observed SESNe with broad light curves but not the peak luminosities. The commonly used analytic formula based on rising bolometric light curves overestimates the ejecta mass of these high-initial-mass progenitor explosions by a factor of up to 2.6. In contrast, the calibrated method by Haynie et al., which relies on late-time decay tails, reduces uncertainties to an average of 20% within the calibrated ejecta mass range. Spectroscopically, the He I 1.083 μm line remains prominent even in models with as little as 0.02 M ⊙ of helium. However, the strength of the optical He I lines is not directly proportional to the helium mass but instead depends on a complex interplay of factors such as the 56 Ni distribution, composition, and radiation field. Thus, producing realistic helium features requires detailed radiative transfer simulations for each new hydrodynamic model.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1990

The modulation of a radar signal from the ocean surface due to slope and hydrodynamic effects

The modulation of the signal level of a tower-based radar scatterometer illuminating a small area of the ocean surface is predicted from direct measurements of the hydrodynamic conditions of the ocean surface and compared with X band radar measurements taken simultaneously. The radar backscatter is assumed to be due to Bragg resonance between the signal and the ripple waves. The effects of the slope of the ocean surface due to passing long waves are modeled by converting a measured wave height series into a slope time series. Laser-slope-gauge measurements of the spectrum of the ripple waves that ride upon the long ocean waves are used to predict the hydrodynamic effects. Comparisons with the measured radar signal show that the slope modulation alone gives a poor prediction of the modulation. Inclusion of the hydrodynamic effects in the simulation appears more promising, although it also is not entirely successful with this data set.

West, James C.↗
NASA NTRS1991

The origin of the moon and the single impact hypothesis. IV

The smoothed-particle hydrodynamics method developed to simulate the conditions under which a major planetary collision may have led to the formation of the moon is presently used in 41 simulations exploring the relevant parameter space. Two possible scenarios emerge: one in which the moon is formed directly in the collision or only after dissipation and evolution of the disk. Attention is given to the Giant Blowoff hypothesis, in which not only the primitive atmosphere but large amounts of mantle material flow over the surface of the orbiting disk and blow off in the forward direction. Large amounts of angular momentum can also be lost.

Cameron, A. G. W.↗
NASA NTRS1994

Merging neutron stars. 1. Initial results for coalescence of noncorotating systems

We present three-dimensional Newtonian simulations of the coalescence of two neutron stars, using a smoothed particle hydrodynamics (SPH) code. We begin the simulations with the two stars in a hard, circular binary, and have them spiral together as angular momentum is lost through gravitational radiation at the rate predicted by modeling the system as two point masses. We model the neutron stars as hard polytropes (gamma = 2.4) of equal mass, and investigate the effect of the initial spin of the two stars on the coalescence. The process of coalescence, from initial contact to the formation of an axially symmetric object, takes only a few orbital periods. Some of the material from the two neutron stars is shed, forming a thick disk around the central, coalesced object. The mass of this disk is dependent on the initial neutron star spins; higher spin rates result in greater mass loss and thus more massive disks. For spin rates that are most likely to be applicable to real systems, the central coalesced object has a mass of 2.4 solar mass, which is tantalizingly close to the maximum mass allowed by any neutron star equation of state for an object that is supported in part by rotation. Using a realistic nuclear equation of state, we estimate the temperature of the material after the coalescence. We find that the central object is at a temperature of approximately 10 MeV, while the disk is heated by shocks to a temperature of 2 to 4 MeV.

Davies, M. B.↗
NASA NTRS2002

Tethers as Debris: Hydrocode Simulation of Impacts of Tether Fragments on Planar Aerospace Materials

Tethers promise to find use in a variety of space applications. Despite being narrow objects, their great lengths result in them having large total areas, and so tethers are quite susceptible to being severed by orbital debris. Extensive work has been done designing tethers that resist severs by small debris objects, and hence have longer working lives. It is from this perspective that most recent work has considered the tether - debris question. The potential of intact tethers, or severed tether fragments, as debris to pose a significant collision risk to other spacecraft has been less well studied. Understanding the consequences of such encounters is important in assessing the risks to other spacecraft posed by tethers. In this paper I discuss the damage that two types of tethers may produce on planar aerospace materials, as revealed by hyper- velocity impact simulations using the SPHC hydrodynamic code. Tether types considered include a single nylon line and a complex design including metal wires. Target materials considered include the aluminum plates typically used in debris shielding, and solar panels.

Evans, Steven W.↗
NASA NTRS2003

Tethers as Debris: Hydrocode Simulation of Impacts of Polymer Tether Fragments on Aluminum Plates

Tethers promise to find use in a variety of space applications. Despite being narrow objects, their great lengths result in them having large total areas. Consequently, tethers are very susceptible to being severed by orbital debris. Extensive work has been done designing tethers that resist severs by small debris objects, in order to lengthen their working lives. It is from this perspective that most recent work has considered the tether - debris question. The potential of intact tethers, or severed tether fragments, as debris, to pose a significant collision risk to other spacecraft has been less well studied. Understanding the consequences of such collisions is important in assessing the risks tethers pose to other spacecraft. This paper discusses the damage that polymer tethers may produce on aluminum plates, as revealed by hypervelocity impact simulations using the SPHC hydrodynamic code.

Evans, Steven W.↗
NASA NTRSDate not supplied

Extending the Applicability of Thermal Protection System Ballistic Limit Equations Beyond the Testable Regime

All spacecraft use some sort of thermal insulation in their design, which is usually placed on the spacecraft’s exterior. As such, it is highly susceptible to high-speed impacts by meteoroids and orbital debris. These high-speed impacts can damage the insulation to a point where the protection it offers is below acceptable limits. It is, therefore, important to be able to characterize expected TPS damage levels stemming from such high-speed impacts. In a recent paper, we presented the results of a study that sought to characterize the high-speed impact damage that would be sustained by two thermal protection materials that have recently gained attention for possible use in interplanetary missions. As part of the analyses performed, ballistic limit equations for these materials were developed using recent high-speed impact test data, and then assessed by comparing their predictions against test data. In this paper, we extend the applicability of the ballistic limit equation for one of those thermal materials beyond the testable regime using the results of series of Smooth Particle Hydrodynamics Code (SPHC) impact simulations at velocities in excess of 10 km/s. When the predictions of the updated and expanded BLE were compared against the results of the SPHC runs, the predictions were seen to be consistent with the data: non-perforation data points were typically below the BLE curves, while perforation data points were typically above them.

Thermal Protection System↗
DOE OSTI2026

Reduced-Order CFD Modeling to Support Waste Loading Optimization in Hanford WTP Vitrification

The U.S. DOE Hanford Site stores over 56 million gallons of radioactive liquid tank waste that must be treated and immobilized for long-term disposal The Waste Treatment and Immobilization Plant (WTP) will vitrify this waste by feeding it into Joule-heated melters, where it is incorporated into a stable borosilicate glass Computational fluid dynamics (CFD) simulations of glass melters can provide insight into the maximum achievable waste loading under varying melter operating conditions Fully resolved VOF multiphase simulations were used as the reference model to capture bubble-driven convection in the melter, including bubble formation, rise behavior, and induced glass melt circulation Effective bubble column diameter and rise velocity were extracted from the resolved simulations, compared with empirical correlations, and refit across relevant viscosity and gas flow rate conditions Explicit gas–liquid interface tracking was replaced with a single-phase momentum source term model, enabling faster steady-state CFD simulations while preserving the dominant hydrodynamic effects of bubbling New empirical correlations were developed for effective bubble column diameter and bubble rise velocity by fitting resolved simulation data across expected melter viscosity and gas flow rate ranges, providing improved inputs for the momentum source term model compared with existing literature correlations The momentum source term model reduced fluid-domain mesh size by 89% and achieved an 8.4× computational speedup relative to resolved bubbling simulations The validated momentum source term approach enables prediction of process-relevant heat transfer behavior in the integrated melter model, including heat transfer from the molten glass to the cold cap, plenum, refractory walls, and surrounding structural regions under varying melter operating conditions

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2022

Evaluation of 3D Simulation Capabilities within xRAGE using a Suite of Test Problems

This presentation describes a 3D verification and validation suite for xRAGE. Three verification test problems for hydrodynamics are presented: the Kidder ball problem, the Verney shell problem, and a 5-material compression problem. There is excellent agreement between 2D and 3D XRAGE simulation results for these problems, and between the xRAGE results and the benchmark solutions. In addition, two 3D ICF test problems are presented, based on an OMEGA direct drive capsule experiment, and on a NIF indirect drive capsule experiment. These two ICF test problems are used to evaluate two hydrodynamic methods within xRAGE. In addition to hydrodynamics the simulations also include three temperature plasma physics, thermal conduction and radiation diffusion. When simulating either capsule in 3D the newer unsplit hydrodynamic method in xRAGE produces more vorticity relative to the older default method. For the indirect drive capsule the 3D simulations are in reasonable agreement with the experimental values of ion temperature and neutron production. Unstable 3D vortex rings are present near the DT ice/carbon interface which enhance the confinement of the DT ice.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

A two-dimensional relativistic (Gamma = 3.25) jet simulation

A numerical simulation of a fully relativistic jet in slab geometry is presented. The simulation uses a covariant divergence formulation of the equations of hydrodynamics. A Mach disk familiar from nonrelativistic simulations is observed. We also find relativistically narrow streaks in the pressure distribution near the aperture.

Van Putten, Maurice H. P. M.↗
NASA NTRS1991

Three-dimensional simulations of black hole tori

The global hydrodynamic nonaxisymmetric instabilities in thick constant angular momentum (l) accretion gas tori in orbit around a Schwarzschild black hole are investigated via 2D and 3D numerical simulations. A radially-wide torus is found to develop and m = 1 nonaxisymmetric density perturbation near the pressure maximum and a trailing spiral wave that extends through the outer part of the torus. This global wave transports angular momentum through the torus, causing the average angular momentum distribution to evolve slowly away from l = constant. The unstable mode drives an accretion flow from the torus into the black hole. The role of accretion is investigated by modeling a wide torus with an inner boundary at the cusp of the Schwarzschild effective potential. In such a torus, accretion is present throughout the evolution; only modest unstable mode growth is observed, and saturation occurs at low amplitude. Comparison simulations performed in the torus equatorial plane show qualitative similarities between the 2D and 3D system, although the 3D modes have a function dependence on height.

Hawley, John F.↗
DOE OSTI2021

Feasibility of the LAMMPS SPH Package for Simulating Bubble Motion in Vibrating Containers

The Smoothed Particle Hydrodynamics (SPH) package within LAMMPS is explored as a possible tool for simulating the motion of bubbles in a vibrating liquid-filled container. As an initial test case, the unphysical but computationally less intense situation of a two-dimensional single bubble rising in a quiescent liquid under the influence of gravity is considered herein. Although physically plausible behavior was obtained under certain conditions, this behavior depends strongly on the system parameters. Moreover, the large density ratio between the liquid and bubble requires extremely small timesteps, which make the simulations undesirably computationally expensive. Ultimately, it was determined that this method is not feasible for providing quantitatively accurate results for the desired application.

42 ENGINEERING↗