Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow physics”

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,423 records · Page 79

Measurements of the Elliptic and Triangular Azimuthal Anisotropies in Central He 3 + Au , d + Au and p + Au Collisions at s N N = 200 GeV

The elliptic (v 2 ) and triangular (v 3 ) azimuthal anisotropy coefficients in central 3He + Au, d + Au, and p + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV are measured as a function of transverse momentum (P T ) at midrapidity (|n| < 0.9) via the azimuthal angular correlation between two particles both at |n| < 0.9 while the v 2 (P T ) values depend on the colliding systems, the v 3 (P T ) values are system independent within the uncertainties, suggesting an influence on eccentricity from subnucleonic fluctuations in these small-sized systems. Furthermore, these results also provide stringent constraints for the hydrodynamic modeling of these systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Viscosities of the Baryon-Rich Quark-Gluon Plasma from Beam Energy Scan Data

Here, this work presents the first Bayesian inference study of the (3 + 1)D dynamics of relativistic heavy-ion collisions and quark-gluon plasma viscosities using an event-by-event (3 + 1)D hydrodynamics + hadronic transport theoretical framework and data from the Relativistic Heavy Ion Collider Beam energy scan program. Robust constraints on initial state nuclear stopping and the baryon chemical potential-dependent shear viscosity of the produced quantum chromodynamic (QCD) matter are obtained. The specific bulk viscosity of the QCD matter is found to exhibit a preferred maximum around $\sqrt{s_{NN}}$ = 19.6 GeV . This result allows for the alternative interpretation of a reduction (and/or increase) of the speed of sound relative to that of the employed lattice-QCD based equation of state for net baryon chemical potential μ B ~0.2(0.4) GeV .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Disentangling Sources of Momentum Fluctuations in Xe+Xe and Pb+Pb Collisions with the ATLAS Detector

High-energy nuclear collisions create a quark-gluon plasma, whose initial condition and subsequent expansion vary from event to event, impacting the distribution of the eventwise average transverse momentum [𝑃⁡([𝑝 T ])]. Disentangling the contributions from fluctuations in the nuclear overlap size (geometrical component) and other sources at a fixed size (intrinsic component) remains a challenge. This problem is addressed by measuring the mean, variance, and skewness of 𝑃⁡([𝑝 T ]) in 208 Pb + 208 Pb and 129 Xe + 129 Xe collisions at $\sqrt{S{NN}}$ = 5.02 and 5.44 TeV, respectively, using the ATLAS detector at the LHC. All observables show distinct features in ultracentral collisions, which are explained by a suppression of the geometrical component as the overlap area reaches its maximum. These results demonstrate a new technique to separate geometrical and intrinsic fluctuations, providing constraints on initial conditions and properties of the quark-gluon plasma, such as the speed of sound.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Upper limit on the chiral magnetic effect in isobar collisions at the Relativistic Heavy-Ion Collider

The chiral magnetic effect (CME) is a phenomenon that arises from the QCD anomaly in the presence of an external magnetic field. The experimental search for its evidence has been one of the key goals of the physics program of the Relativistic Heavy-Ion Collider. The STAR Collaboration has previously presented the results of a blind analysis of isobar collisions ( Ru 44 96 + Ru 44 96 , Zr 40 96 + Zr 40 96 ) in the search for the CME. The isobar ratio ( Y ) of CME-sensitive observable, charge separation scaled by elliptic anisotropy, is close to but systematically larger than the inverse multiplicity ratio, the naive background baseline. This indicates the potential existence of a CME signal and the presence of remaining nonflow background due to two- and three-particle correlations, which are different between the isobars. In this postblind analysis, we estimate the contributions from those nonflow correlations as a background baseline to Y , utilizing the isobar data as well as Heavy Ion Jet Interaction Generator simulations. This baseline is found consistent with the isobar ratio measurement, and an upper limit of 10% at 95% confidence level is extracted for the CME fraction in the charge separation measurement in isobar collisions at s NN = 200 GeV. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Evidence for Similar Collectivity of High Transverse-Momentum Particles in 𝑝-Pb and Pb-Pb Collisions

Charged hadron elliptic anisotropies (𝑣 2 ) are presented over a wide transverse momentum (𝑝 T ) range for proton-lead (𝑝 ⁢Pb) and lead-lead (PbPb) collisions at nucleon-nucleon center-of-mass energies of 8.16 and 5.02 TeV, respectively. The data were recorded by the CMS experiment and correspond to integrated luminosities of 186 and 0.607 nb −1 for the 𝑝⁢Pb and PbPb systems, respectively. A four-particle cumulant analysis is performed using subevents separated in pseudorapidity to effectively suppress noncollective effects. At high 𝑝 T (𝑝 T > 8 GeV), significant positive 𝑣 2 values that are similar between 𝑝 ⁢Pb and PbPb collisions at comparable charged particle multiplicities are observed. This observation suggests a common origin for the multiparticle collectivity for high-𝑝 T particles in the two systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Experimental search for the chiral magnetic effect in relativistic heavy-ion collisions: A perspective

The chiral magnetic effect (CME) refers to generation of the electric current along a magnetic field in a chirally imbalanced system of quarks. The latter is predicted by quantum chromodynamics to arise from quark interaction with nontrivial topological fluctuations of the vacuum gluonic field. The CME has been actively searched for in relativistic heavy-ion collisions, where such gluonic field fluctuations and a strong magnetic field are believed to be present. The CME-sensitive observables are unfortunately subject to a possibly large non-CME background, and firm conclusions on a CME observation have not yet been reached. In this perspective, we review the experimental status and progress in the CME search, from the initial measurements more than a decade ago to the dedicated program of isobar collisions in 2018 and the release of the isobar blind analysis result in 2022 to intriguing hints of a possible CME signal in Au + Au collisions, and discuss future prospects of a potential CME discovery in the anticipated high-statistic Au + Au collision data at the Relativistic Heavy-Ion Collider by 2025. We hope such a perspective will help sharpening our focus on the fundamental physics of the CME and steer its experimental search.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Sample Complexity of Power System State Estimation Using Matrix Completion

In this paper, we propose an analytical framework to quantify the amount of data samples needed to obtain accurate state estimation in a power system - a problem known as sample complexity analysis in computer science. Motivated by the increasing adoption of distributed energy resources into the distribution-level grids, it becomes imperative to estimate the state of distribution grids in order to ensure stable operation. Traditional power system state estimation techniques mainly focus on the transmission network which involve solving an overdetermined system and eliminating bad data. However, distribution networks are typically underdetermined due to the large number of connection points and high cost of pervasive installation of measurement devices. In this paper, we consider the recently proposed state-estimation method for underdetermined systems that is based on matrix completion. In particular, a constrained matrix completion algorithm was proposed, wherein the standard matrix completion problem is augmented with additional equality constraints representing the physics (namely power-flow constraints). We analyze the sample complexity of this general method by proving an upper bound on the sample complexity that depends directly on the properties of these constraints that can lower number of needed samples as compared to the unconstrained problem. To demonstrate the improvement that the constraints add to state estimation, we test the method on a 141-bus distribution network case study and compare it to the traditional least squares minimization state estimation method.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A modern concept of Lagrangian hydrodynamics

Here, we offer a modern interpretation of Lagrangian hydrodynamics as employed in Lagrangian simulations of compressible fluid flow. Our main result is to show that artificial viscosity, traditionally viewed as a numerical artifice to control unphysical oscillations in flows with shocks, actually represents a physical process and is necessary to derive accurate simulations in any compressible flow. We begin by reviewing the origins of two numerical devices, artificial viscosity and finite-volume methods. We proceed to construct a mathematical (PDE) model that incorporates those numerics and in which a new length scale, the observer, arises representing the discretization. Associated with that length scale, there are new inviscid fluxes that are the artificial viscosity as first formulated by Richtmyer and an artificial heat flux postulated by Noh but typically not included in Lagrangian codes. We discuss the connection of our results to bivelocity hydrodynamics. We conclude with some speculation as to the direction of future developments in multidimensional Lagrangian codes as computers get faster and have larger memories.

97 MATHEMATICS AND COMPUTING↗

TEAL

TEAL is a financial performance calculator plugin for the RAVEN code, framework, resolving around the computation of Net Present Value and associated financial metrics. TEAL can make use of inflation rates, taxation, escalation factors, capital expenditure economy of scale scaling factors. The unique feature of TEAL is the capability to be linked with RAVEN external models and build corresponding cash flows using the variables computed by those external models. In addition to be able to use the capability to generate cash flows derived from complex physical models generated by RAVEN, another distinctive feature of TEAL is the capability to provide financial risk/probabilistic metrics that can empower RAVEN to perform optimization/analysis driven by financial risk augmentations. Optimization, robust optimization, parametric studies, large parallel simulations, sensitivity analysis, data mining, etc. are just some of the capabilities that can be leveraged.

Alfonsi, Andrea↗

Implicit SPH v. 1.0

ISPH code is developed to solve multi-physics meso-scale flow problems using implicit SPH method. In particular, the code can provides solutions for incompressible, multi phase flow and electro-kinetic flows.

Kim, Kyungjoo↗

GEFSv12 High- and Low-Skill Day-10 Tornado Forecasts

Abstract On average, modern numerical weather prediction forecasts for daily tornado frequency exhibit no skill beyond day 10. However, in this extended-range lead window, there are particular model cycles that have exceptionally high forecast skill for tornadoes because of their ability to correctly simulate the future synoptic pattern. Here, model initial conditions that produced a more skillful forecast for tornadoes over the United States were exploited while also highlighting potential causes for low-skill cycles within the Global Ensemble Forecasting System, version 12 (GEFSv12). There were 88 high-skill and 91 low-skill forecasts in which the verifying day-10 synoptic pattern for tornado conditions revealed a western U.S. thermal trough and an eastern U.S. thermal ridge, a favorable configuration for tornadic storm occurrence. Initial conditions for high skill forecasts tended to exhibit warmer sea surface temperatures throughout the tropical Pacific Ocean and Gulf of Mexico, an active Madden–Julian oscillation, and significant modulation of Earth-relative atmospheric angular momentum. Low-skill forecasts were often initialized during La Niña and negative Pacific decadal oscillation conditions. Significant atmospheric blocking over eastern Russia—in which the GEFSv12 overforecast the duration and characteristics of the downstream flow—was a common physical process associated with low-skill forecasts. This work helps to increase our understanding of the common causes of high- or low-skill extended-range tornado forecasts and could serve as a helpful tool for operational forecasters. Significance Statement This research provides a framework for the anticipation of a more (or less) skillful 10-day tornado forecast in an operational numerical weather prediction system. High-skill forecasts were associated with substantial tropical convection and warm sea surface temperature throughout the Pacific Ocean and Gulf of Mexico, whereas the underlying cause of low-skill forecasts were typically associated with a blocking anticyclone over eastern Russia. These findings are important because they permit increased or decreased confidence in a long-range forecast of tornado occurrence based on a dynamical prediction system.

Meteorology & Atmospheric Sciences↗

Relationship between Anisotropy in Soil Hydraulic Conductivity and Saturation

Anisotropy in unsaturated hydraulic conductivity is saturation-dependent. Accurate characterization of soil anisotropy is very important in simulating flow and contaminant (e.g., radioactive nuclides in Hanford) transport. A recently developed tensorial connectivity-tortuosity (TCT) concept describes the hydraulic conductivity tensor of the unsaturated anisotropic soils as the product of a scalar variable, the symmetric connectivity tortuosity tensor, and the hydraulic conductivity tensor at saturation. In this study, the TCT model is used to quantify soil anisotropy in unsaturated hydraulic conductivity. The TCT model can describe different types of soil anisotropy; e.g., the anisotropy coefficient, C, can be monotonically increase or decrease with saturation and can vary from greater than unity to less than unity and vice versa. Soil anisotropy is independent of soil water retention properties and can be characterized by the ratio of the saturated hydraulic conductivities and the difference of the tortuosity-connectivity coefficients in two directions. ln(C) is linearly proportional to ln(Se) with Se being the effective saturation. The log-linear relationship between C and Se allows the saturation-dependent anisotropy to be determined using linear regression with the measurements of the directional hydraulic conductivities at a minimum of two water content levels, of which one may be at full saturation. The model was tested using measurements of directional hydraulic conductivities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Real-time POD-CFD Wind-Load Calculator for PV Systems

The primary objective of this project is to create an accurate web-based real-time wind-load calculator. This is of paramount importance for (1) the rapid and accurate assessments of the uplift and downforce loads on a PV mounting system, (2) identifying viable solutions from available mounting systems, and therefore helping reduce the cost of mounting hardware and installation. Wind loading calculations for structures are currently performed according to the American Society of Civil Engineers/ Structural Engineering Institute Standard ASCE/SEI 7; the values in this standard were calculated from simplified models that do not necessarily take into account relevant characteristics such as those from full 3D effects, end effects, turbulence generation and dissipation, as well as minor effects derived from shear forces on installation brackets and other accessories. This standard does not include provisions that address the special requirements of rooftop PV systems, and attempts to apply this standard may lead to significant design errors as wind loads are incorrectly estimated. Therefore, an accurate calculator would be of paramount importance for the preliminary assessments of the uplift and downforce loads on a PV mounting system, identifying viable solutions from available mounting systems, and therefore helping reduce the cost of the mounting system and installation. The challenge is that although a full-fledged three-dimensional computational fluid dynamics (CFD) analysis would properly and accurately capture the complete physical effects of air flow over PV systems, it would be impractical for this tool, which is intended to be a real-time web-based calculator. CFD routinely requires enormous computation times to arrive at solutions that can be deemed accurate and grid-independent even in powerful and massively parallel computer platforms.

14 SOLAR ENERGY↗

Large Eddy Simulation Modeling of Flashback and Flame Stabilization in Hydrogen-Rich Gas Turbines Using a Hierarchical Validation Approach

This project was a combined computational and experimental effort to improve predictive capability for boundary layer flashback of premixed swirl flames relevant to gas-turbine power plants operating with high-hydrogen-content fuels. During the course of this project, significant progress in modeling was made on four major fronts: 1) use of direct numerical simulation of turbulent flames to understand the coupling between the flame and the turbulent boundary layer; 2) improved modeling capability for flame propagation in stratified pre-mixtures; 3) improved portability of computer codes using the OpenFOAM platform to facilitate transfer to industry and other researchers; and 4) application of LES to flashback in swirl combustors, and a detailed assessment of its capabilities and limitations for predictive purposes. A major component of the project was an experimental program that focused on developing a rich experimental database of boundary layer flashback in swirl flames. Both methane and high-hydrogen fuels, including effects of elevated pressure (1 to 5 atm), were explored. For this project, a new model swirl combustor was developed. Kilohertz-rate stereoscopic PIV and chemiluminescence imaging were used to investigate the flame propagation dynamics. In addition to the planar measurements, a technique capable of detecting the instantaneous, time-resolved 3D flame front topography was developed and applied successfully to investigate the flow-flame interaction. The UT measurements and legacy data were used in a hierarchical validation approach where flows with increasingly complex physics were used for validation. First component models were validated with DNS and literature data in simplified configurations, and this was followed by validation with the UT 1-atm flashback cases, and then the UT high-pressure flashback cases. The new models and portable code represent a major improvement over what was available before this project was initiated.

42 ENGINEERING↗

PHENIX results on centrality dependence of yields and correlations in d plus Au collisions at root s(NN)=200 GeV

PHENIX has measured the transverse momentum (pT) spectra and two particle angular correlations for high pT particles in d+Au collisions at psNN=200 GeV using the RHIC Year-2008 run data. The azimuthal angle correlations for two particles with a large rapidity gap exhibit a ridge-like structure. Using the pi-0s reconstructed in the EMCal, we have successfully extended the pT reach of the correlation up to 8 GeV/c. We find that the azimuthal anisotropy of hadrons found at low pT persists up to 6 GeV/c with a significant centrality and pT dependence, similar to what was observed in A+A collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fluid Mechanics of an Obliquely Mounted MIV Gauge

Over the 40-yr history of shock-compression science, numerous physical phenomena have been considered for use in detecting wave profiles. Few of the devices have actually been used for a significant and persistent study.

47 OTHER INSTRUMENTATION↗

Physics based Reliability Models for Supercritical CO 2 Turbomachinery Components

CSP plants using sCO 2 power cycle can potentially achieve high thermal conversion efficiency at low capital cost due to compact turbo-machinery and other components. An sCO 2 expander and improved heat exchanger is expected to provide a major stepping stone for achieving CSP power at $0.06/kW-hr LCOE, energy conversion efficiency > 50%, and total power block cost < $1,200/kW installed. Hybrid gas bearings (HGB) and Dry gas seals (DGS) are two key components that are critical for robust and efficient operation of future sCO 2 power cycles. From a system perspective, HGBs can provide substantial benefits (~3 to 5%) to the modular CSP operation, while use of DGSs limits the leakage losses as well as the windage losses in the generator and can increase overall generator efficiency by ~ 8%. The life limiting mechanisms of HGBs and DGSs in high pressure, high temperature sCO2 environment were not well understood and hence this program aims at developing physics-based performance prediction and life prediction models for these two components. The output of the performance models serves as input to the life models. The performance models were validated using experimental data and using data from literature. DGS and HGB models have shed light on some risks associated with its operation in sCO2 expander and will guide more robust designs tailored towards sCO 2 cycle applications. A DGS was tested in 1MWe sCO2 turbo expander to study its performance and validate the performance models. To understand the effect of high pressures, high temperatures and sCO 2 environment on crack initiation, crack propagation and LCF/HCF life of these components, a novel experimental setup was developed. The LCF test data was used to calibrate and validate the life prediction models. The effect of sCO 2 corrosion on the life of Ni-based super-alloys is quantified using long duration experiments. Advanced microstructure and spectroscopic analyses conducted shed light on some key differences between IN718 and IN617 in terms of oxidation morphology, chemical species diffusion and trapping, the formation of protective corrosion resistant layers and changes in surface properties. A novel crack initiation model was developed which takes into account strain distribution, temperature, hold times and activation energy at the gas-metal interface and predicts the number of cycles to initiate the crack in the Inconel alloy. An advanced 3D fracture mechanics-based crack propagation model was also developed which takes into account the transient loading mission of the component and predict how fast the crack will grow and what will be its final shape and size. High energy X-Ray tomography measurements were conducted to study and validate 3D morphology of crack growth inside the specimen. These two models along with chemical kinetics models obtained from corrosion tests can be combined to predict the total life of the component under given long term operating conditions. A Bayesian Hybrid Modeling (BHM) probabilistic framework was developed to quantify uncertainty in physics-based performance models and in crack initiation and crack propagation models, to calibrate the models and to validate the models with statistical confidence. Advanced coupled physics multi-phase flow models were developed to study condensation, boiling and phase change behavior of CO 2 under various sCO 2 turbomachinery and heat exchanger conditions. A novel experimental setup was also designed to observe phase change of CO 2 during its transition across liquid-gas-supercritical state lines.

14 SOLAR ENERGY↗

Modeling and Simulation Capability Development Plan

This report provides an assessment of the current status of the modelling and simulation capability for the Nuclear-Renewable Hybrid Energy System (N-R HES). The current capabilities are separated in: • Creation of the system boundary condition. More in detail this refer to the capability to create either demand, renewable supply, price profiles etc.. • Creation of a numerical representation of the physical model of a N-R HES system • Derivation, from the actual performance of the physical system, the cash flow to perform the economical assessment • Optimization of N-R HES configuration and dispatch Those capabilities are the ones which enable the workflow summarize in Figure 1below.

97 MATHEMATICS AND COMPUTING↗