Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics modeling”

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

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

At least 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Propagation of the 2014–2016 Northeast Pacific Marine Heatwave Through the Salish Sea

Effects and impacts of the Northeast Pacific marine heatwave of 2014–2016 on the inner coastal estuarine waters of the Salish Sea were examined using a combination of monitoring data and an established three-dimensional hydrodynamic and biogeochemical model of the region. The anomalous high temperatures reached the U.S. Pacific Northwest continental shelf toward the end of 2014 and primarily entered the Salish Sea waters through an existing strong estuarine exchange. Elevated temperatures up to + 2.3°C were observed at the monitoring stations throughout 2015 and 2016 relative to 2013 before dissipating in 2017. The hydrodynamic and biogeochemical responses to this circulating high-temperature event were examined using the Salish Sea Model over a 5-year window from 2013 to 2017. Responses of conventional water-quality indicator variables, such as temperature and salinity, nutrients and phytoplankton, zooplankton, dissolved oxygen, and pH, were evaluated relative to a baseline without the marine heatwave forcing. The simulation results relative to 2014 show an increase in biological activity (+14%, and 6% Δ phytoplankton biomass, respectively) during the peak heatwave year 2015 and 2016 propagating toward higher zooplankton biomass (+14%, +18% Δ mesozooplankton biomass). However, sensitivity tests show that this increase was a direct result of higher freshwater and associated nutrient loads accompanied by stronger estuarine exchange with the Pacific Ocean rather than warming due to the heatwave. Strong vertical circulation and mixing provided mitigation with only ≈+0.6°C domain-wide annual average temperature increase within Salish Sea, and served as a physical buffer to keep waters cooler relative to the continental shelf during the marine heatwave.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Constraining the P 30 ( p , γ ) S 31 Reaction Rate in ONe Novae via the Weak, Low-Energy, β -Delayed Proton Decay of Cl 31

The 30 P(p,γ) 31 S reaction plays an important role in understanding the nucleosynthesis of A ≥ 30 nuclides in oxygen-neon novae. The Gaseous Detector with Germanium Tagging was used to measure 31 Cl β-delayed proton decay through the key J π = 3/2 + , 260-keV resonance. The intensity $I$$^{260}_{βp}$ = $8.3$$^{+1.2}_{–0.9}$ × 10 –6 represents the weakest β-delayed, charged-particle emission ever measured below 400 keV, resulting in a proton branching ratio of Γ p /Γ = $2.5$$^{+0.4}_{–0.3}$ × 10 –4 . Here, by combining this measurement with shell-model calculations for Γγ and past work on other resonances, the total 30 P(p,γ) 31 S rate has been determined with reduced uncertainty. The new rate has been used in hydrodynamic simulations to model the composition of nova ejecta, leading to a concrete prediction of 30 Si: 28 Si excesses in presolar nova grains and the calibration of nuclear thermometers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

thornado-transport: Anderson- and GPU-accelerated nonlinear solvers for neutrino-matter coupling

Algorithms for neutrino-matter coupling in core-collapse supernovae (CCSNe) are investigated in the context of a spectral two-moment model, which is discretized in space with the discontinuous Galerkin method, integrated in time with implicit-explicit (IMEX) methods, and implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado). The model considers electron neutrinos and antineutrinos and tabulated opacities from Bruenn (1985), which includes neutrino-electron scattering and pair processes. The nonlinear system arising from implicit time discretization of the equations governing neutrino-matter coupling is iterated to convergence using Anderson-accelerated fixed-point methods, which avoid formation of Jacobians and inversion of dense linear systems. Numerical experiments show that, for a given tolerance, a nested iteration scheme which aims to reduce opacity evaluations can lower the computational cost. Our initial port to GPUs, using both OpenMP and OpenACC, shows an overall speedup of up to ~ 100× when compared to results using a single CPU core. These results indicate that the algorithms implemented in thornado are well-suited to GPU acceleration.

Laiu, Paul↗
DOE OSTI2022

Simulation of dislocation evolution in microparticle impacts over a wide range of impact velocities

We report the dynamic loading environment created by high velocity microparticle impacts is very difficult to model. The extremely high strain rates, large deformations, and high temperatures affect the flow stress of the material in ways that standard flow stress models cannot capture. Plastic deformation and strain hardening in metals is controlled by the motion of dislocations. Dislocations can be nucleated, stored as forest dislocations, or be annihilated as loading progresses. A comprehensive accounting of dislocation density change is needed to accurately describe rate and temperature dependent dislocation glide and evolution across the wide rage of loading conditions present in microparticle impact problems. Therefore, we implement and apply a newly proposed flow stress model (Hunter and Preston, 2015, 2022) that is then coupled with mobile and immobile dislocation density evolution equations. This model is implemented in Los Alamos National Laboratory’s hydrodynamics code, FLAG, to model copper-on-copper microparticle impacts. This new strength model allows for accurate simulation of particle rebound and flattening across a wide range of impact velocities.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Hydrodynamics in long-range interacting systems with center-of-mass conservation

In systems with a conserved density, the additional conservation of the center of mass (dipole moment) has been shown to slow down the associated hydrodynamics. At the same time, long-range interactions generally lead to faster transport and information propagation. Here, we explore the competition of these two effects and develop a hydrodynamic theory for long-range center-of-mass-conserving systems. Here we demonstrate that these systems can exhibit a rich dynamical phase diagram containing subdiffusive, diffusive, and superdiffusive behaviors, with continuously varying dynamical exponents. We corroborate our theory by studying quantum lattice models whose emergent hydrodynamics exhibit these phenomena.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Improved single-particle phase-space distributions for viscous fluid dynamic models of relativistic heavy ion collisions

Comparison of hydrodynamic calculations with experimental data inevitably requires a model for converting the fluid to particles. For shear viscous fluids, several different conversion models have been proposed, such as the quadratic Grad corrections, the Strickland-Romatschke (SR) ansatz, self-consistent shear corrections from linearized kinetic theory, or corrections from the relaxation time approach. We show, via comparison to nonlinear 2 → 2 kinetic theory evolution, that even for low specific shear viscosities 0.03 ≲ $\textit{η/s}$ ≲ 0.2 , these four models have varying accuracy and still significant errors in reproducing particle phase space densities from hydrodynamic fields. Moreover, we demonstrate that the overall reconstruction error of additive shear viscous $f = f_{eq} + δf$ models dramatically reduces by up to one order of magnitude, if one ensures through exponentiation that $\textit{f}$ is always positive. Further, we also illustrate how even more accurate viscous $\textit{δf}$ models can be constructed by incorporating the first time derivatives of hydrodynamic fields to account for the past evolution of the system. Such time derivatives are readily available in hydrodynamic simulations, though usually not included in the output. Though our comparisons are limited to a one-component massless system undergoing a (0 + 1)-dimensional (0 + 1D) longitudinal boost-invariant expansion, we expect that the improved models will be useful in 2 + 1D and 3 + 1D hydrodynamic studies as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Quantifying Nitrogen Bioextraction by Seaweed Farms – A Real-time Modeling-Monitoring Case Study in Hood Canal, WA

Anthropogenic nitrogen input is the leading cause of increasing marine eutrophication. Cultivating seaweeds at large scales may provide a viable solution for mitigating eutrophication especially when nitrogen extraction is achieved as an added environmental benefit of biofuel and food production. As the first step to promoting and implementing this solution, it is critical to quantify and assess the effectiveness of nitrogen removal by existing seaweed farms. Pacific Northwest National Laboratory and Hood Canal Mariculture Inc. conducted a case study by applying a coupled modeling-monitoring approach to quantify nitrogen bioextraction by a one-hectare commercial sugar kelp farm in Hood Canal, WA. Through this two-year effort, the project team conducted intensive field monitoring of water quality parameters at the farm site and subsequent laboratory analysis. Meanwhile, the team applied a coupled hydrodynamic-macroalgal growth model to simulate kelp growth and nitrogen bioextraction over the two kelp growing seasons. The model results showed an overall good comparison with field monitoring data and showed the effectiveness of nitrogen removal by kelp farms. The established modeling-monitoring approach can also be transferable for supporting broader research activities in macroalgal farming and eutrophication mitigation.

09 BIOMASS FUELS↗
DOE OSTI2022

An arbitrary Lagrange-Eulerian investigation of HRAM shallow jet pre-spurt formation and time sensitivities to impact plate dynamics

With dry bay fires persisting as a significant contributor to aircraft vulnerability despite chronicled developments in survivability technologies, an accurate fire prediction capability remains paramount for credible vulnerability assessments. Physics based modeling of the hydrodynamic ram (HRAM) fluid deposition process is a key component of such capability, wherein capturing the first instance of fluid spurt, referred to herein as shallow jet spurts, is a core focus. Such pre-spurts as they have been formerly identified have only been witnessed sporadically in HRAM spurt experiments. In this work, ALE3D, a first-principles multi-physics code was employed to model the shallow jet spurt phenomenon with spherical projectiles impacting water-filled tanks faced with aluminum panels, such that the underlying physics and sensitivities could be explored. Development and verification of the 2D axisymmetric model is described relative to trends observed in a prior experimental campaign. Results from the verified model suggest shallow jet spurts have at least a quadratic sensitivity to the fundamental vibrational mode of the impact plate across impact velocities 610–1829 m/s (2000–6000 ft/s). It is further explained how shallow jet spurts are arrested for impacts at the two extremes of plate rigidity. This research constitutes the first fluid-structure modeling of shallow jet spurts which future three-dimensional analyses will expound upon.

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

A Numerical Study of Hypoxia in Chesapeake Bay Using an Unstructured Grid Model: Validation and Sensitivity to Bathymetry Representation

A three-dimensional unstructured-grid hydrodynamic and water quality model (Semi-implicit Cross-scale Hydroscience Intergrated System Model-Integrated Compartment Model) is applied successfully for Chesapeake Bay. The model is validated with observations of salinity, chlorophyll-a, dissolved oxygen, nutrients, and phytoplankton productions from the year 1991 to 1995 for the mainstem and some major tributaries, based on multiple model skill scores. Model experiments are conducted to test the importance of having (1) an accurate representation of bathymetry to correctly predict hypoxia and other processes and (2) a high-resolution model grid for tributaries to correctly simulate water quality variables. Further, comparison with the model experiment results with bathymetry smoothing indicates that bathymetry smoothing, as commonly used for many systems, changes the stratification and lateral circulation pattern, resulting in more salt intrusion into shallow water regions, and an increase in the freshwater age. Consequently, a model with bathymetry smoothing can lead to an unrealistic prediction of the distribution of hypoxia and phytoplankton production. Local grid refinement shows significant improvement of model simulations on local stratification and water quality variables. Overall, the use of high-resolution unstructured grid model leads to a faithful representation of the complex geometry, and thus a seamless cross-scale capability for simulating water quality processes in the Bay including tributaries and tidal creeks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Disentangling the Effects of Vapor Pressure Deficit and Soil Water Availability on Canopy Conductance in a Seasonal Tropical Forest During the 2015 El Niño Drought

In land models used in Earth system models, vegetation water stress is heuristically related to soil water matric potential, and root water uptake in soil layers is simply proportional to root density to meet the transpiration demand. The efficiency of roots utilizing deeper soil water under water stress cannot be well represented, but it is of critical importance for transpiration regulation and carbon flux. To improve modeling of tropical forest response to drought, a plant hydrodynamic scheme (HYDRO), adapted from a trait-based model in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), has been implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) to represent the hydraulic dynamics and root water uptake strategy. The adapted model conforms to the big-leaf photosynthesis model in ELM. Vegetation water stress in the model is a function of prognostic leaf water potential and is used to regulate stomatal conductance. Point simulations were performed at two tropical forest sites in the Amazon and Panama. Model sensitivity to the numerical solver, plant hydraulic traits, soil properties, and atmospheric forcing has been evaluated. Our results showed that, even though water stress factors calculated at our study sites using HYDRO are comparable to that simulated by the heuristic water stress model, soil moisture prediction is improved with soil water redistributions through roots. Simulated results using HYDRO is not as sensitive to plant hydraulic traits measured for different tree species as to the soil properties and atmospheric forcing. Simplified numerical solver may cause inaccuracy as to where roots uptake water. This study suggests that plant hydrodynamic representation is critical for improved simulations of soil moisture under water stresses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Understanding the compound flood risk along the coast of the contiguous United States

Abstract. Compound flooding is a type of flood event caused by multiple flood drivers. The associated risk has usually been assessed using statistics-based analyses or hydrodynamics-based numerical models. This study proposes a compound flood (CF) risk assessment (CFRA) framework for coastal regions in the contiguous United States (CONUS). In this framework, a large-scale river model is coupled with a global ocean reanalysis dataset to (a) evaluate the CF exposure related to the coastal backwater effects on river basins, and (b) generate spatially distributed data for analyzing the CF hazard using a bivariate statistical model of river discharge and storm surge. The two kinds of risk are also combined to achieve a holistic understanding of the continental-scale CF risk. The estimated CF risk shows remarkable inter- and intra-basin variabilities along the CONUS coast with more variabilities in the CF hazard over the US west and Gulf coastal basins. Different risk assessment methods present significantly different patterns in a few key regions such as the San Francisco Bay area, the lower Mississippi River, and Puget Sound. Our results highlight the need to weigh different CF risk measures and avoid using single statistics-based or hydrodynamics-based CFRAs. Uncertainty sources in these CFRAs include the use of gauge observations, which cannot account for the flow physics or resolve the spatial variability of risks, and underestimations of the flood extremes and the dependence of CF drivers in large-scale models, highlighting the importance of understanding the CF risks for developing a more robust CFRA.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Viscosity Measurements in Extreme Conditions

Radiation-hydrodynamics simulations are critical to the NNSA complex for understanding high energy density physics and for addressing problems in national security science. However, current rad-hydro codes typically do not account for viscosity, which generally leads to large uncertainties in problems involving mixing of materials under shock loading. This report details progress made towards experimental measurements of fluid viscosity at high pressures and temperatures achieved by shock wave compression techniques, as well as complementary material model development and hydrodynamic simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Characterization of thermal transport and evolution of Au plasma in ICF experiments by Thomson scattering

Here, this paper demonstrates the capability of optical Thomson scattering (OTS) to measure thermal transport, and provides support to radiation hydrodynamic and kinetic simulations of electron thermal transport and plasma evolution. OTS theory and plasma simulations are applied to the interpretation of experimental measurements of laser-produced plasma from spherical gold targets on the OMEGA laser facility. The dynamical form factor, S($\vec{k}$, ω), of electron density fluctuations that is used in the fitting of Thomson scattering spectra includes ion–ion collisions and effects of non-Maxwellian distribution functions. OTS measurements and their interpretation are consistent with the nonlocal transport model in radiation hydrodynamic simulations as well as with kinetic simulations in the second half of the probe pulse duration. In particular, the reversal of heat transport during cooling is observed to be consistent with simulations, while some discrepancies are noted during the initial heating of the Au targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Effects of alpha-ion stopping on ignition and ignition criteria in inertial confinement fusion experiments

With the advent of ignited plasmas at the National Ignition Facility (NIF), alpha physics has become a driving factor in theoretical understanding and experimental behavior. In this communication, we explore aspects of direct alpha-ion heating through comparison of the consequences from the one-fluid and two-fluid models in the hydrodynamic approach. We show that the case with all alpha energy deposited in electrons raises the ignition criteria by ∼4 keV or ∼0.2 g/cm2 in the hotspot relative to the case with all alpha energy deposited in ions. In the case of the recently ignited NIF implosion, 30% of the 3.5 MeV α energy is deposited into the DT fuel ions, for which there is negligible difference between the one-fluid and two-fluid ignition criteria. However, changes in the ion stopping fraction through profile effects and alternate stopping power models could lead to ignition curve shifts of ∼1 keV.

Physics↗
DOE OSTI2023

Pinwheel Experiment Confirmatory Shot 1 Pre-Shot Report

The purpose of this experimental series is to validate the reactive burn model in the PAGOSA hydrodynamic code. It will also help validate the EOS and Strength models for the aluminum and tantalum cylinders. C4 high explosive (HE) is used to drive an aluminum cylinder into an aluminum sandwiched PBX-9501. The follow-on experiment will drive the aluminum cylinder into a tantalum sandwiched PBX-9501. The experiment will provide time of arrival (TOA) data for the burn front in the HE through the metal cylinders using Fiber Bragg Gratings (FBG). Photon Doppler Velocimetry (PDV) will provide shock wave TOA and the velocity of the inward moving cylinder. The PAGOSA code will be used to predict TOA and cylinder wall velocity using PDV and tracers.

42 ENGINEERING↗