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

Model study of the energy dependence of the correlation between anisotropic flow and the mean transverse momentum in Au + Au collisions

A hybrid model that employs the hadron-string transport model UrQMD and the (3+1)-dimensional relativistic viscous hydrodynamic code vHLLE, is used to investigate the beam energy dependence of the correlation coefficient ρ(v$^{2}_{2}$,[p T ]) between the average transverse momentum [p T ] of hadrons emitted in an event and the square of the anisotropic flow coefficient v$^{2}_{2}$. For Au+Au collisions, the model predicts characteristic patterns for the energy and event-shape dependence of the variances for [p T ] and v$^{2}_{n}$ [Var([p T ] and Var(v$^{2}_{2}$)], and the covariance of v$^{2}_{n}$ and [p T ][cov(v$^{2}_{2}$,[p T ])], consistent with the attenuation effects of the specific shear viscosity η/s. Here in contrast, ρ(v$^{2}_{2}$,[p T ]) is predicted to be insensitive to the beam energy but sensitive to the initial-state geometry of the collisions. These observations suggest that a precise set of measurements for Var([p T ]), Var(v$^{2}_{2}$), cov(v$^{2}_{2}$,[p T ]) and ρ(v$^{2}_{2}$,[p T ]) as a function of beam energy and event shape, could serve to constrain better the eccentricities and their fluctuations over the entire span of beam energies, and consequently aid precision extraction of the temperature and baryon chemical-potential dependence of η/s from the wealth of Au+Au data obtained in the RHIC beam energy scan.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

A boundary integral formalism for stochastic ray tracing in billiards

Determining the flow of rays or non-interacting particles driven by a force or velocity field is fundamental to modelling many physical processes. These include particle flows arising in fluid mechanics and ray flows arising in the geometrical optics limit of linear wave equations. In many practical applications, the driving field is not known exactly and the dynamics are determined only up to a degree of uncertainty. This paper presents a boundary integral framework for propagating flows including uncertainties, which is shown to systematically interpolate between a deterministic and a completely random description of the trajectory propagation. A simple but efficient discretisation approach is applied to model uncertain billiard dynamics in an integrable rectangular domain.

42 ENGINEERING↗
DOE OSTI2022

Ratios of collective flow observables in high-energy isobar collisions are insensitive to final-state interactions

The ratios of bulk observables, such as harmonic flow v 2 and v 3 , between high-energy 96 Ru + 96 Ru and 96 Zr + 96 Zr collisions were recently argued to be a clean probe of the nuclear structure differences between 96 Ru and 96 Zr . Using a transport model simulation of isobar collisions, we quantify this claim from the dependence of the ratios v 2,Ru /v 2, Zr and v 3,Ru /v 3,Zr on various final-state effects, such as the shear viscosity, hadronization, and hadronic cascade. Although the v 2 and v 3 change by more than 50% when varying the final-state effects, the ratios are unchanged. In addition, these ratios are independent of the transverse momentum p T and hadron species, despite of up to a factor of 2 change in v n . Further, the ratio of mean transverse momentum $\langle$p T $\rangle$ is found to be controlled by the nuclear skin and nuclear radius but is only slightly impacted by the final-state effects. Therefore, these isobar ratios serve as a clean probe of the initial condition of the quark-gluon plasma, which, in turn, is controlled by the collective structure of the colliding nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Experimental Investigation of Uranium and Iron Condensation from High-Temperature Plasma Conditions

We used a plasma flow reactor (PFR) to generate synthetic fallout nanoparticles from vapor-phase condensation of two different input concentrations of uranium and iron analytes (U/Fe = 1:1 and 1:2). Synthetic fallout from complex chemical matrices (e.g., mixtures of U and Fe) has not been generated in this setup before and allows for the observation of relative condensation and fractionation of nuclear debris. Experiments were conducted under two different temperature histories with variations in particle flow patterns along the PFR. Transmission electron microscopy (TEM) observation and analysis of the nanoparticles revealed variations in speciation of uranium oxides (UO 2 and α-UO 3 ) depending on the competition between flow mixing and oxygen sequestration by iron. A ternary metal oxide, UFeO 4 , was observed in addition to iron oxides (e.g., FeO and Fe 3 O 4 ), which suggests that fallout models should account for chemical speciation of ternary metal oxides (i.e., UFeO 4 ) and their relative condensation behaviors in addition to those of singular metal oxides (e.g., FeO and UO 2 ). X-ray Energy Dispersive Spectroscopy (EDS) elemental maps showed that some particles had Fe-rich cores surrounded by U-rich regions. This suggests that either condensed U oxides coagulate onto molten Fe oxides or that the increase in iron analyte concentration might drive the system toward a higher degree of supersaturation, leading to earlier formation of iron oxide particles and providing an energetically favored pathway for nucleation of uranium oxides around iron oxide particles.

and nuclear chemistry↗
DOE OSTI2023

Strange hadron collectivity in pPb and PbPb collisions

The collective behavior of K$_{S}^{0}$ and $Λ/ \overline{Λ}$ strange hadrons is studied by measuring the elliptic azimuthal anisotropy ($\mathcal{v}$ 2 ) using the scalar-product and multiparticle correlation methods. Proton-lead (pPb) collisions at a nucleon-nucleon center-of-mass energy $\sqrt{s_{NN}}$ = 8.16 TeV and lead-lead (PbPb) collisions at $\sqrt{s_{NN}}$ = 5.02 TeV collected by the CMS experiment at the LHC are investigated. Nonflow effects in the pPb collisions are studied by using a subevent cumulant analysis and by excluding events where a jet with transverse momentum greater than 20 GeV is present. The strange hadron $\mathcal{v}$ 2 values extracted in pPb collisions via the four- and six-particle correlation method are found to be nearly identical, suggesting the collective behavior. Comparisons of the pPb and PbPb results for both strange hadrons and charged particles illustrate how event-by-event flow fluctuations depend on the system size.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2007

Runtime and Pressurization Analyses of Propellant Tanks

Multi-element unstructured CFD has been utilized at NASA SSC to carry out analyses of propellant tank systems in different modes of operation. The three regimes of interest at SSC include (a) tank chill down (b) tank pressurization and (c) runtime propellant draw-down and purge. While tank chill down is an important event that is best addressed with long time-scale heat transfer calculations, CFD can play a critical role in the tank pressurization and runtime modes of operation. In these situations, problems with contamination of the propellant by inclusion of the pressurant gas from the ullage causes a deterioration of the quality of the propellant delivered to the test article. CFD can be used to help quantify the mixing and propellant degradation. During tank pressurization under some circumstances, rapid mixing of relatively warm pressurant gas with cryogenic propellant can lead to rapid densification of the gas and loss of pressure in the tank. This phenomenon can cause serious problems during testing because of the resulting decrease in propellant flow rate. With proper physical models implemented, CFD can model the coupling between the propellant and pressurant including heat transfer and phase change effects and accurately capture the complex physics in the evolving flowfields. This holds the promise of allowing the specification of operational conditions and procedures that could minimize the undesirable mixing and heat transfer inherent in propellant tank operation. It should be noted that traditional CFD modeling is inadequate for such simulations because the fluids in the tank are in a range of different sub-critical and supercritical states and elaborate phase change and mixing rules have to be developed to accurately model the interaction between the ullage gas and the propellant. We show a typical run-time simulation of a spherical propellant tank, containing RP-1 in this case, being pressurized with room-temperature nitrogen at 540 R. Nitrogen, shown in blue on the right-hand side of the figures, enters the tank from the diffuser at the top of the figures and impinges on the RP-1, shown in red, while the propellant is being continuously drained at the rate of 1050 lbs/sec through a pipe at the bottom of the tank. The sequence of frames in Figure 1 shows the resultant velocity fields and mixing between nitrogen and RP-1 in a cross-section of the tank at different times. A vortex is seen to form in the incoming nitrogen stream that tends to entrain propellant, mixing it with the pressurant gas. The RP-1 mass fraction contours in Figure 1 are also indicative of the level of mixing and contamination of the propellant. The simulation is used to track the propagation of the pure propellant front as it is drawn toward the exit with the evolution of the mixing processes in the tank. The CFD simulation modeled a total of 10 seconds of run time. As is seen from Figure 1d, after 5.65 seconds the propellant front is nearing the drain pipe, especially near the center of the tank. Behind this pure propellant front is a mixed fluid of compromised quality that would require the test to end when it reaches the exit pipe. Such unsteady simulations provide an estimate of the time that a high-quality propellant supply to the test article can be guaranteed at the modeled mass flow rate. In the final paper, we will discuss simulations of the LOX and propellant tanks at NASA SSC being pressurized by an inert ullage. Detailed comparisons will be made between the CFD simulations and lower order models as well as with test data. Conditions leading to cryo collapse in the tank will also be identified.

Field, Robert E.↗
DOE OSTI2025

Characterizing IHE Response to Multiple Shock Loading

The response of high explosives to shock loading is traditionally measured with a steady loading pressure. In many accident scenarios involving fragment impact, however, a loading duration that is shorter than the build up to detonation may occur. Fragments passing through multiple materials before reaching a high explosive charge may produce loading that is comprised of more than one shock wave. Additionally, the build up to detonation in high explosive corner turning loads the explosive a short duration pressure pulse, since rarefactions can often rapidly overtake the reactive wave. For these reasons, we have studied the response of the insensitive high explosive (IHE) materials PBX 9502 and LX-17 to complex loadings of varied intensity and duration. We refer to a single loading of limited duration as a “thin pulse”, whereas more complex scenarios were studied with an impactor that produces a double shock in the explosive. The following report presents experimental data and analyses of thin pulse shock initiation and double shock experiments designed to guide development of models of Insensitive High Explosives (IHEs) under controlled one-dimensional conditions relevant to accident scenarios and corner turning. Thin pulse shock initiation data on PBX 9502 and LX-17 were obtained under varied pulse duration, pressed density, and temperature conditions in order to probe various parameters essential for the development of a physics-based Cheetah reactive flow hotspot model. In situ pressure gauges provide insight into the degree of reaction in the explosive that are not obtainable with optical PDV measurements or distance measurements such as run to detonation. Double shock data was obtained to inform a Composition Aware Cheetah model which can be applied to any TATB-based IHEs. This model supports efforts to find a new IHE formulation and potentially incorporate new binders into IHE formulations. Simulations of each experiment are included to demonstrate the utility of these focused experiments to developing models of HE behavior. One-dimensional gas gun experiments are essential for characterizing shocked HE behavior and informing HE models.

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

Collision-System and Beam-Energy Dependence of Anisotropic Flow Fluctuations

Here, elliptic flow measurements from two-, four- and six-particle correlations are used to investigate flow fluctuations in collisions of U+U at $\sqrt{s_{NN}}$ = 193 GeV, Cu+Au at $\sqrt{s_{NN}}$ = 200 GeV and Au+Au spanning the range $\sqrt{s_{NN}}$ = 11.5 - 200 GeV. The measurements show a strong dependence of the flow fluctuations on collision centrality, a modest dependence on system size, and very little if any, dependence on particle species and beam energy. The results, when compared to similar LHC measurements, viscous hydrodynamic calculations, and Glauber model eccentricities, indicate that initial-state-driven fluctuations predominate the flow fluctuations generated in the collisions studied.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Measurement of kaon directed flow in Au+Au collisions in the high baryon density region

Rapidity-odd directed flow v 1 measurements are presented for $K^±$ and $K^0_S$ in Au + Au collisions for $\sqrt{s_{NN}}$ from 3.0 to 3.9 GeV with the STAR experiment. For comparison, v 1 of π ± , protons, and Λ from the same collisions are also discussed. The mid-rapidity v 1 slope dv 1 /dy| y=0 for protons and Λ is positive in these collisions. On the other hand, v 1 slope of kaons exhibits a strong dependence: negative at p T < 0.6 GeV/c and positive at higher p T . A similar p T dependence is also evident for the v 1 slope of charged pions. Compared to the spectator-removed calculations in Au+Au collisions at $\sqrt{s_{NN}}$ = 3.0–3.9 GeV, the JAM model demonstrates a pronounced shift of the v 1 slopes of mesons towards the negative direction. It suggests that the shadowing effect of the spectators plays an important role in the observed kaon anti-flow at low p T in the high baryon density region of non-central collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

On velocity and reactive scalar spectra in turbulent premixed flames

Abstract Kinetic energy and reactive scalar spectra in turbulent premixed flames are studied from compressible three-dimensional direct numerical simulations (DNS) of a temporally evolving rectangular slot-jet premixed flame, a statistically one-dimensional configuration. The flames correlate to a lean premixed hydrogen–air mixture at an equivalence ratio of 0.7, preheated to 700 K and at 1 atm, and three DNS are considered with a fixed jet Reynolds number of 10 000 and a jet Damköhler number varying between 0.13 and 0.54. For the study of spectra, motivated by the need to account for density change, which can be locally strong in premixed flames, a new density-weighted definition for two-point velocity/scalar correlations is proposed. The density-weighted two-point correlation tensor retains the essential properties of its constant-density (incompressible) counterpart and recovers the density-weighted Reynolds stress tensor in the limit of zero separation. The density weighting also allows the derivation of balance equations for velocity and scalar spectrum functions in the wavenumber space that illuminate physics unique to combusting flows. Pressure–dilatation correlation is a source of kinetic energy at high wavenumbers and, analogously, reaction rate–scalar fluctuation correlation is a high-wavenumber source of scalar energy. These results are verified by the spectra constructed from the DNS data. The kinetic energy spectra show a distinct inertial range with a $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}-5/3$ scaling followed by a ‘diffusive–reactive’ range at higher wavenumbers. The exponential drop-off in this range shows a distinct inflection in the vicinity of the wavenumber corresponding to a laminar flame thickness, $\delta _L$ , and this is attributed to the contribution from the pressure–dilatation term in the energy balance in wavenumber space. Similarly, a clear spike in spectra of major reactant species (hydrogen) arising from the reaction-rate term is observed at wavenumbers close to $\delta _L$ . It appears that in the inertial range classical scaling laws for the spectra involving the Kolmogorov scale are applicable, but in the high-wavenumber range where chemical reactions have a strong signature the laminar flame thickness produces a better collapse. It is indicated that a full scaling should perhaps involve the Kolmogorov scale, laminar flame thickness, Damköhler number and Karlovitz number.

42 ENGINEERING↗
DOE OSTI2021

Thermal conductivity modeling of monodispersed microspheres using discrete element method

Particle beds are widely used in various systems and processes, such as particle heat exchangers, granular flow reactors, and additive manufacturing. Accurate modeling of the thermal conductivity of particle beds and understanding of their heat transfer mechanisms are important. However, previous models were based on simple cubic packing of particles, which could not accurately represent the actual heat transfer processes under certain conditions. Here, we examine the effect of the packing structure on the thermal conductivity of particle beds. We use monodispersed silica microspheres with average particle sizes ranging from 23 to 330 μm as a model material. We employ a transient hot-wire technique to measure the thermal conductivity of the particle beds with packing density of 43%–57% within a temperature range of room temperature to 500 °C and under N2 gaseous pressures of 20–760 Torr. We then use a discrete element method (DEM) to obtain the realistic packing structure of the particles, which is then fed into a finite-element model (FEM) to calculate the thermal conductivity, with the consideration of solid conduction, gas conduction, and radiation heat transfer. Our results show that the thermal conductivity model based on the more realistic random packing structure derived from the DEM shows better agreement with the experimental data compared to that based on the simple cubic-packing structure. The combined DEM and FEM methodology can serve as a useful tool to predict the effective thermal conductivity of particle beds and to quantify different heat transfer mechanisms under various conditions.

14 SOLAR ENERGY↗
DOE OSTI2019

Nanoparticle diffusion in sheared cellular blood flow

Using a multiscale blood flow solver, the complete diffusion tensor of nanoparticles (NPs) in sheared cellular blood flow is calculated over a wide range of shear rate and haematocrit. In the short-time regime, NPs exhibit anomalous dispersive behaviors under high shear and high haematocrit due to the transient elongation and alignment of the red blood cells (RBCs). In the long-time regime, the NP diffusion tensor features high anisotropy. Particularly, there exists a critical shear rate ( ${\sim}100~\text{s}^{-1}$ ) around which the shear-rate dependence of the diffusivity tensor changes from linear to nonlinear scale. Above the critical shear rate, the cross-stream diffusivity terms vary sublinearly with shear rate, while the longitudinal term varies superlinearly. The dependence on haematocrit is linear in general except at high shear rates, where a sublinear scale is found for the vorticity term and a quadratic scale for the longitudinal term. Through analysis of the suspension microstructure and numerical experiments, the nonlinear haemorheological dependence of the NP diffusion tensor is attributed to the streamwise elongation and cross-stream contraction of RBCs under high shear, quantified by a capillary number. The RBC size is shown to be the characteristic length scale affecting the RBC-enhanced shear-induced diffusion (RESID), while the NP submicrometre size exhibits negligible influence on the RESID. Based on the observed scaling behaviours, empirical correlations are proposed to bridge the NP diffusion tensor to specific shear rate and haematocrit. As a result, the characterized NP diffusion tensor provides a constitutive relation that can lead to more effective continuum models to tackle large-scale NP biotransport applications.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Reaction plane correlated triangular flow in Au + Au collisions at $\sqrt{s_{NN}}$ = 3 GeV

Here, we measure triangular flow relative to the reaction plane at 3 GeV center-of-mass energy in Au + Au collisions at the BNL Relativistic Heavy Ion Collider. A significant v 3 signal for protons is observed, which increases for higher rapidity, higher transverse momentum, and more peripheral collisions. The triangular flow is essentially rapidity-odd with a slope at midrapidity, dv 3 /dy| (y=0) , opposite in sign compared to the slope for directed flow. No significant v 3 signal is observed for charged pions and kaons. Comparisons with models suggest that a mean field potential is required to describe these results, and that the triangular shape of the participant nucleons is the result of stopping and nuclear geometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗