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

Observation of ionization trends in a laboratory photoionized plasma experiment at Z

We report experimental and modeling results for the charge state distribution of laboratory photoionized neon plasmas in the first systematic study over nearly an order of magnitude range of ionization parameter ξ ∝ F/N e . The range of ξ is achieved by flexibility in the experimental platform to adjust either the x-ray drive flux F at the sample or the electron number density N e or both. Experimental measurements of photoionized plasma conditions over such a range of parameters enable a stringent test of atomic kinetics models used within codes that are applied to photoionized plasmas in the laboratory and astrophysics. From experimental transmission data, ion areal densities are extracted by spectroscopic analysis that is independent of atomic kinetics modeling. The measurements reveal the net result of the competition between photon-driven ionization and electron-driven recombination atomic processes as a function of ξ as it affects the charge state distribution. Results from radiation-hydrodynamics modeling calculations with detailed inline atomic kinetics modeling are compared with the experimental results. There is good agreement in the mean charge and overall qualitative similarities in the trends observed with ξ but significant quantitative differences in the fractional populations of individual ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Model comparisons for two-temperature plasma equations of state

When a plasma is generated in the laboratory, energy is often deposited preferentially into either the electrons or the ions, giving rise to a quasiequilibrium state in which the two species, electrons and ions, are well described by two effective temperatures, T e and T i . Accurate hydrodynamic modeling of such a two-temperature plasma requires an equation of state that captures the relevant many-body physics without assuming a strict local thermodynamic equilibrium. Several models have been proposed within the literature, which extend conventional statistical approaches, each employing a different combination of assumptions for modifying the equilibrium equations. Here, in this work, we compare the predictions for several models, presenting derivations of the internal energy and pressure for each microscopic model within a unified framework so that the assumptions of each model may be more easily compared to one another. We find that for sufficiently weak coupling, all models agree with one another. However, as the coupling strength is increased, the disagreement between the models becomes more pronounced. Moreover, the relative sizes of the corrections predicted by each model depend on which species has the higher temperature, T e >T i vs T e < T i .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Embedding a critical point in a hadron to quark-gluon crossover equation of state

Lattice QCD simulations have shown unequivocally that the transition from hadrons to quarks and gluons is a crossover when the baryon chemical potential is zero or small. Many model calculations predict the existence of a critical point at a value of the chemical potential where current lattice simulations are unreliable. We show how to embed a critical point in a smooth background equation of state so as to yield the critical exponents and critical amplitude ratios expected of a transition in the same universality class as the liquid-gas phase transition and the three-dimensional Ising model. There are only two independent critical exponents; the relations α + 2β + γ = 2 and β(δ-1) = γ arise automatically, as does a relation between the two critical amplitudes. The resulting equation of state has parameters that may be inferred by hydrodynamic modeling of heavy-ion collisions in the Beam Energy Scan II at the BNL Relativistic Heavy Ion Collider or in experiments at other accelerators.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Slot-die-coating operability windows for polymer electrolyte membrane fuel cell cathode catalyst layers

Roll-to-roll (R2R) slot-die coating of polymer electrolyte membrane fuel cell (PEMFC) catalyst layers represents a scalable deposition method for producing 10–20 m 2 ·min -1 of catalyst-coated gas diffusion layers (GDLs). This high-throughput production technique will help lower the cost of PEMFC catalyst layers. The uniformity of the wet layer applied by slot die deposition is affected by process parameters such as substrate speed, vacuum pressure applied at the upstream meniscus, gap between the slot die lips and substrate, ink rheology, and other ink and substrate properties. The set of conditions for producing a defect-free coating with a dilute ink typically requires little to no upstream vacuum pressure, so suitable operating conditions can be found easily through trial and error and operator intuition. However, the higher viscosity of more concentrated inks dramatically shifts the range of settings that result in a homogeneous coating to higher vacuum levels, which are harder to find through hit or miss. A predictive model showing the range of operable conditions decreases material wastage inherent in experimentally searching for suitable parameters. In this study, the defect-free coating parameter window is explored experimentally and theoretically for two concentrations of PEFC cathode inks. Both a full capillary hydrodynamic model and a computationally cheaper viscocapillary model successfully predict the experimentally determined coating window within the experimental and model uncertainty limits for inks with 5.3 wt% and 12.0 wt% solids ink while maintaining the 0.1 mg Pt ·cm -2 Pt areal loading target. Here, this paper demonstrates a viable pathway for meeting the $30/kW net ultimate cost target of the United States Department of Energy (U.S. DOE) Hydrogen Fuel Cells Technologies Office (HFTO). The concentrated ink lowers the thermal energy and capital expenditure (CapEx) budget of the coating process by decreasing the amount of time, energy, and floorspace required for drying the coating.

25 ENERGY STORAGE↗
DOE OSTI2022

Embedding a Critical Point in a Hadron to Quark–Gluon Crossover Equation of State

It is shown how to embed a critical point in a smooth background equation of state so as to yield the critical exponents and critical amplitude ratios expected of a transition in the same universality class as the liquid–gas phase transition and the 3D Ising model. There are only two independent critical exponents; the relations α + 2β + $\gamma$ = 2 and β(δ - 1) = $\gamma$ arise automatically, as does a new relation between the two critical amplitudes. The resulting equation of state has parameters which may be inferred by hydrodynamic modeling of heavy-ion collisions in the Beam Energy Scan II at the Relativistic Heavy Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Explicit quantification of residence and flushing times in the Salish Sea using a sub-basin scale shoreline resolving model

The Salish Sea, located in the Pacific Northwest region of North America has complex currents and circulation features distributed over numerous interconnected deep basins with islands. Increased risk of exposure to oil spills and untreated wastewater from maritime emergencies and treatment plant failures have led to a need for quantifying residence and flushing characteristics at a sub-basin scale using a shoreline resolving hydrodynamic model. An unstructured grid model of the Salish Sea was developed using FVCOM with a ≈75-100m shoreline resolution. In addition to 7 tides and 23 salinity and temperature monitoring stations, an extensive currents data set from 135 stations collected over a span of three years was used for skill assessment and validation. Explicit forward computations were then conducted to define and quantify residence and flushing times in various sub-basins of interest using (a) Lagrangian particles and (b) Numerical/virtual dye experiments. The results in most basins show expected seasonal variability with longer flushing time associated with summer lower tides and lower freshwater inflows. However, contrary to expectation, flushing time is significantly longer in wintertime in large fjord-like basins such as Hood Canal (≈138 days), likely due to increased stratification and reduced mixing. The flushing time for the Puget Sound region of the Salish Sea is ≈115 days, while Georgia Basin is 240 days when analyzed as stand-alone basins with zero background concentrations. When examined as part of the flushing of the entire Puget Sound filled with virtual dye, the compounded flushing times for embedded sub-basins can be significantly longer in order of magnitudes and largely dictated by the flushing time of Puget Sound. The computed residence and flushing time scales tabulated over 36 sub-basins provide an improved understanding of water renewal in the system, informing pollution management actions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Tracking the environmental impacts of ecological engineering on coastal wetlands with numerical modeling and remote sensing

Coastal wetlands are the most valuable ecosystems on the earth but facing severe degradation and losses owing to climate change and anthropogenic activities. Many ecological engineering projects (EEP) have been conducted to mitigate the degradation of coastal wetlands. However, the geomorphological impacts of EEP on coastal wetlands have not been well documented. In this study, a method employed a process-based hydrodynamic model and remote sensing (RS) was developed to evaluate the impacts of EEP on the geomorphological change of a prototype Ramsar site. Results demonstrated that RS has great potential in improving the quality of bathymetry data for the numerical model with a decrease of RMSE from 0.52 m to 0.3 m. It also showed good capacity in trend detection of geomorphological change spatially. Results showed the Chongming Dongtan wetland experienced erosion with an annual rate of -0.035 m/yr from 2013 to 2016 after the implementation of EEP. The deposition rate changed significantly in the area within 200 m of the EEP. It is found that the EEP modified the composition of vegetation, sediment transportation, as well as substrate stability, affecting the geomorphological change of coastal wetlands. It is suggested that the EEP with moderate anthropogenic disturbance is a direct and effective way to recover the coastal habitats for waterbirds. However, the modification of the coastal wetland ecosystem by EEP will lead to the potential vulnerability to global climate change. Therefore, how to mitigate the advantages and disadvantages of the EEP is needed to be further studied to find a more sustainable way for coastal management.

geomorphological change, coastal wetlands, Ecologi↗
DOE OSTI2026

Assessing shellfish water exposure to fecal bacteria pollution in Salish Sea: three-dimensional modeling and implications for monitoring

Fecal bacteria (FB) contamination poses significant risks to shellfish safety and management in coastal and estuarine waters. Despite extensive pollution identification and correction efforts, FB contamination in shellfish-growing areas persists in the Salish Sea, highlighting the need to identify overlooked sources and better understand FB transport from riverine and shoreline inputs to shellfish beds. To address this, a high-resolution three-dimensional hydrodynamic model coupled with FB kinetics was developed and applied to a case study site in Salish Sea—Portage Bay—to simulate freshwater plume circulation, flushing dynamics, and bacterial transport. Daily FB loading from the major freshwater inflow—Nooksack River was generated by both linear interpolation and integrating a machine learning approach (XGBoost), trained on historical hydrological and meteorological data. The model successfully reproduced both the magnitude and seasonal variation of FB concentrations in Portage Bay for the year of 2021, demonstrating that simplified FB kinetics with first-order decay due to mortality was effective in this dynamic coastal environment with short flushing time. Model results identified the Nooksack River as the dominant far-field FB source, while scenario simulations showed that near-field coastal stormwater outfalls elevated local FB levels following rainfall, particularly under low-flow conditions. The XGBoost prediction provided comparable or superior accuracy to linear interpolation, particularly during periods of missing observational data, by capturing short-term variability and event-driven loading more effectively. Integrating data-driven riverine FB inputs with mechanistic coastal numerical modeling provides a robust framework for operational forecasting of shellfish bed exposure risk and supports adaptive monitoring and management of shellfish growing areas in the Salish Sea and similar coastal systems.

Salish Sea↗
DOE OSTI2023

Measurement of φ-meson production in Cu + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV and U + U collisions at $\sqrt{s_{NN}}$ = 193 GeV

The PHENIX experiment reports systematic measurements at the Relativistic Heavy Ion Collider of φ-meson production in asymmetric Cu + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV and in U + U collisions at $\sqrt{s_{NN}}$ = 193 GeV. Measurements were performed via the φ → K + K – decay channel at midrapidity |η| < 0.35. Features of φ- meson production measured in Cu + Cu, Cu + Au, Au + Au, and U + U collisions were found to not depend on the collision geometry, which was expected because the yields are averaged over the azimuthal angle and follow the expected scaling with nuclear-overlap size. The elliptic flow of the φ meson in Cu + Au, Au + Au, and U + U collisions scales with second-order-participant eccentricity and the length scale of the nuclear-overlap region (estimated with the number of participating nucleons). At moderate $p_T$ , φ-meson production measured in Cu + Au and U + U collisions is consistent with coalescence-model predictions, whereas at high $p_T$ the production is in agreement with expectations for in-medium energy loss of parent partons prior to their fragmentation. Finally, the elliptic flow for φ mesons measured in Cu + Au and U + U collisions is well described by a (2+1)-dimensional viscous-hydrodynamic model with specific-shear viscosity η/s = 1/4π.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Model Advancements to Enable Impact Analysis of Climate Change on Streamflow Temperature

With support from the Department of Energy’s Water Power Technologies Office, Pacific Northwest National Laboratory (PNNL) has developed new tools that incorporate cutting-edge climate and hydrological science capabilities to assess the potential long-term impacts of future climate conditions on unregulated streamflow and water temperature within watershed-river-reservoir systems. The objectives of this project were achieved by enhancing key hydrologic and hydrodynamic models and transferring them to a high-performance computing environment to provide a high-spatiotemporal resolution, multi-scale modeling framework. The new modeling framework has the potential to quantify risks related climate change impacts on runoff, unregulated streamflow, and water temperature. Initial development and demonstration of the modeling framework was conducted under historical and future climate conditions in the Columbia River Basin in the Pacific Northwest and the Connecticut River Basin in New England.

13 HYDRO ENERGY↗
DOE OSTI2021

Numerical Model of IProTech PIP WEC Device

iProTech PIP wave energy converter (WEC) is a slack moored, single hull device with no moving parts in the water, joints or bearings. This submission includes data of the simulation, reports, and code for the iProTech PIP (WEC) project. The organization of the data included in the provided archive is detailed below and in the data description of the archive. The data teamer-iprotech-nrel folder includes and explains matlab and python code developed to hydrodynamically model the PIP WEC device in WEC-Sim. The subfolders cover the following steps: 1) report: explanatory information on device geometry 2) pip_mesher: python code to generate mesh panels from device profile data 3) wec-sim_models: matlab code to run WEC-Sim The data uploaded is a snapshot as of 11/02/2121 of code residing in a Github repository administered by David Ogden of NREL.

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↗
DOE OSTI2024

Flow resistance of phloem sieve plates revisited using an experimental model

Sieve plate resistance has been identified as the key to understanding the efficiency of phloem transport of carbohydrates in trees. These plates with small holes connecting sieve tubes are responsible for the largest resistance to flow in the phloem. Their structure determines how fast sugars can be transported through the phloem with certain pressure differences and what the limits for phloem transport in different plants are. Because experiments with the phloem are very challenging, our understanding of sieve plate resistance is mostly based on anatomical studies and hydrodynamic modeling of flow through sieve plates. These models calculate the resistance of the entire sieve tube–sieve plate system using the Hagen–Poiseuille flow resistance for the sieve cell lumen and a combination of the Hagen–Poiseuille resistance and Sampson flow resistance through the sieve pores. The resistance of the entire sieve plate is calculated by summing the Hagen–Poiseuille and Sampson flow resistances of each pore. To test the validity of this model formulation, an experimental model with aspect parameters similar to phloem sieve tubes was built using polyvinyl chloride (PVC) piping and plastic straws of different diameters and lengths. This system was used to measure flow rates and calculate flow resistance at Reynolds numbers 0.5–300. The results suggest that the current models may significantly overestimate the flow resistance caused by sieve plates and that the resistance might be better described by formulations used for perforated plates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Hydrological controls of a riparian wetland based on stable isotope data and model simulations

Isotopic evidence of groundwater and stream water is frequently used to investigate water exchanges with groundwater. Monthly sampling of rain, stream water, and groundwater was conducted at Tims Branch watershed in South Carolina for the oxygen and hydrogen stable isotope (δ 2 H and δ 18 O) measurement, as well as pH and oxidation–reduction potential (ORP). Together with a mass balance perspective, it was determined that it takes a few weeks to one month for groundwater in the hyporheic zone to fully exchange with stream water. From hydrodynamic modelling, we show that substantial (up to 70 %) groundwater exchange occurs at gaining and losing sites. Groundwater exfiltration, i.e. inflow into stream water, contributes up to 4 % to stream water, with the remainder from upstream exfiltration. A 2–4 % per day renewal rate of adjacent groundwater would indirectly indicate a groundwater residence time in the order of half a month to a full month (assuming either a well-mixed case or large dispersion rate in pulse flow case), in agreement with a greatly reduced variability of δ 2 H and δ 18 O of groundwater compared to stream water and rain. This reduced variability of stable isotope signal from groundwater confirms our hypothesis that riparian groundwater mixing at Tims Branch is more of a mixed type rather than a pulse flow type. As a result, a monthly time scale is sufficient for groundwater to become anoxic at exit points into stream water resulting in the episodic production of natural organic matter- and iron-rich flocs upon oxidation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Investigating laser ablated plume dynamics of carbon and aluminum targets

Recently acquired high-resolution images of nanosecond laser ablation plumes suggest a strong correlation between the internal plume structure and the type of material being ablated. However, the details of this relation are currently not well understood. In this work, we attempt to explore this correlation using a 2D radiation hydrodynamics model to study the dependence of internal plume structure formation on the ablation material. Here spatio-temporal emission maps and plume expansion velocities from experimental measurements are compared with the model predictions, including synthetic emission maps. The shape and expansion rate of an outer air plume region are found to be in good agreement for both carbon and aluminum, as are the inner material plume dynamics for carbon ablation. The largest disagreement is observed in the case of a polished aluminum target, where the chaotic inner plume features seen in the experimental images are not observed in the model. The possible physical mechanisms responsible for this discrepancy are discussed. This effort constitutes a continued development toward a predictive model of ablation plume dynamics and chemistry for various materials in extreme environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Survey of Multimaterial Treatments for Thermal Radiative Transfer

Arbitrary Lagrangian-Eulerian methods are a popular choice for hydrodynamic modeling in radiation (rad-hydro) simulations. Because these methods involve a relaxation step that moves the mesh relative to material boundaries, multimaterial spatial zones are generally present. Accurate treatments of these zones are needed to resolve various physical phenomena of interest for inertial confinement fusion applications. However, these codes are often paired with single-material, deterministic thermal radiative transfer (TRT) codes that are oblivious to the material compositions of each zone. These single-material TRT codes can only accept homogenized material properties (opacities, specific heats, etc.) from the hydrodynamic code and output homogenized solutions. After each TRT time step, the multimaterial hydrodynamic code must dehomogenize the quantities computed by the TRT package in order to update subzonal material temperatures. The process by which hydrodynamic codes perform this dehomogenization has not been well documented in previous literature, and the methods can vary significantly from code to code. The purpose of this paper is to document, study, and compare existing techniques used for rad-hydro simulations as well as present a new method with potentially promising results. We summarize several methods and give comparisons on infinite-medium problems as well a finite-medium problem for two of the methods.

42 ENGINEERING↗
DOE OSTI2024

Understanding spatial and temporal drivers of variation in tree hydraulic processes and their consequences for climate feedbacks (Final Technical Report)

This is the final technical report from the first phase of a project that changed institutions. The grant was titled “Understanding spatial and temporal drivers of variation in tree hydraulic processes and their consequences for climate feedbacks.” The overall objectives of this project were to (1) provide model‐compatible datasets of key plant hydraulic traits and status for model evaluation, parameterization and validation and (2) use these data to pinpoint ecosystem responses to a changing hydroclimate by addressing both long‐term climatic drying and episodic extreme droughts. We planned to address the objectives with three research activities to quantify plant responses to chronic water stress and episodic drought: (1) generate high frequency observations of soil and plant hydraulic data across different landscape positions at multiple sites, (2) quantify plant hydraulic trait plasticity in response to experimental soil moisture reduction in situ in two central hardwood forests, and (3) simulate the carbon consequences of incorporating plant hydrodynamics and plant acclimation to water stress in the DOE‐sponsored plant hydrodynamics model FATES‐HYDRO. As of the transfer of this project to another institution, we had made substantial progress on activities 1 and 2, and started activity 3.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Elliptic flow of charged particles at midrapidity relative to the spectator plane in Pb–Pb and Xe–Xe collisions

Measurements of the elliptic flow coefficient relative to the collision plane defined by the spectator neutrons v 2 {$Ψ_{SP}$} in collisions of Pb ions at center-of-mass energy per nucleon–nucleon pair $\sqrt{s_{NN}}=$ 2.76 TeV and Xe ions at $\sqrt{s_{NN}}=$ 5.44 TeV are reported. The results are presented for charged particles produced at midrapidity as a function of centrality and transverse momentum for the 5–70% and 0.2–6 GeV/c ranges, respectively. The ratio between v 2 {$Ψ_{SP}$} and the elliptic flow coefficient relative to the participant plane v 2 {4}, estimated using four-particle correlations, deviates by up to 20% from unity depending on centrality. This observation differs strongly from the magnitude of the corresponding eccentricity ratios predicted by the T R ENTo and the elliptic power models of initial state fluctuations that are tuned to describe the participant plane anisotropies. The differences can be interpreted as a decorrelation of the neutron spectator plane and the reaction plane because of fragmentation of the remnants from the colliding nuclei, which points to an incompleteness of current models describing the initial state fluctuations. A significant transverse momentum dependence of the ratio v 2 {$Ψ_{SP}$} /v 2 {4} is observed in all but the most central collisions, which may help to understand whether momentum anisotropies at low and intermediate transverse momentum have a common origin in initial state fluctuations. The ratios of v 2 {$Ψ_{SP}$} and v 2 {4} to the corresponding initial state eccentricities for Xe–Xe and Pb–Pb collisions at similar initial entropy density show a difference of (7.0 ± 0.9)% with an additional variation of +1.8% when including RHIC data in the T R ENTo parameter extraction. These observations provide new experimental constraints for viscous effects in the hydrodynamic modeling of the expanding quark–gluon plasma produced in heavy-ion collisions at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗