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,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Measurements of second-harmonic Fourier coefficients from azimuthal anisotropies in p + p , p + Au , d + Au , and He 3 + Au collisions at s N N = 200 GeV

Recently, the PHENIX Collaboration has published second- and third-harmonic Fourier coefficients v 2 and v 3 for midrapidity (|η|< 0.35 ) charged hadrons in 0%–5% central p+Au, d+Au, and 3 He+Au collisions at $\sqrt{s_{NN}}$= 200 GeV, utilizing three sets of two-particle correlations for two detector combinations with different pseudorapidity acceptance [Acharya et al., Phys. Rev. C 105, 024901 (2022)]. Here, this paper extends these measurements of v 2 to all centralities in p+Au, d+Au, and 3 He+Au collisions, as well as p+p collisions, as a function of transverse momentum (p T ) and event multiplicity. The kinematic dependence of v 2 is quantified as the ratio R of v 2 between the two detector combinations as a function of event multiplicity for 0.5 < p T <1 and 2 T < 2.5 GeV/c. A multiphase-transport (AMPT) model can reproduce the observed v 2 in most-central to midcentral d Au and 3 He+Au collisions. However, the AMPT model systematically overestimates the measurements in p+p, p+Au, and peripheral d+Au and 3 He+Au collisions, indicating a higher nonflow contribution in the AMPT model than in the experimental data. The AMPT model fails to describe the observed R for 0.5 < p T < 1 GeV/c , but there is qualitative agreement with the measurements for 2 < p T < 2.5 GeV/c.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Proton and neutron density distributions at supranormal density in low- and medium-energy heavy-ion collisions. II. Central Pb + Pb collisions

This paper represents a continuation of our investigation into the limitations on total, proton, and neutron particle number densities, as well as the asymmetry of proton and neutron density distributions achievable in central heavy-ion collisions. We explore these aspects at low and medium energies within the framework of the Boltzmann-Uhlenbeck-Uehling (pBUU) transport and time-dependent Hartree-Fock (TDHF) models. Previous studies, focusing on symmetric and asymmetric collisions of Ca and Sn nuclei [1], and initial results on Pb-nuclei collisions [2], emphasized the role of the Coulomb interaction in these events. Our findings indicated that: (i) the highest total densities predicted at 𝐸 beam = 800 MeV/nucleon were on the order of ≈ 2.5⁢𝜌 0 (𝜌 0 = 0.16 fm −3 ), (ii) the proton-neutron asymmetry for maximal densities, 𝛿=(𝜌$^{max}_{n}$− 𝜌$^{max}_{p}$)/(𝜌$^{max}_{n}$+𝜌$^{max}_{p}$), did not generally exceed the asymmetry in the initial state of the collision at all beam energies and tended to decrease during the reaction, and (iii) a significant portion of this asymmetry had its microscopic origin in Coulomb forces, masking the pure nuclear contribution. These new findings, particularly relevant in the astrophysical context, are further examined in this work, focusing on the heaviest target-projectile combination 212,208 Pb accessible in an experiment. We introduce the SkT3 Skyrme force model, not previously used for the Pb system, and compare it to the SV-bas and SV-sym34 models to explore the symmetry-energy dependence of the results. Contour plots of nucleonic densities are presented, contrasting the time evolution of the density distributions in low (TDHF) and high (pBUU) models. We also present the evolution of normalized maximal proton, neutron, and total nucleon number density with increasing beam energy in the full pBUU model and the Vlasov approximation, aiming to explore the impact of correlations in the reaction. The time evolution of the proton and neutron density distributions in the plane transverse to the beam direction is illustrated at both low and high beam energy. In conclusion, our detailed examination of the Pb system in this work provides further essential evidence that the aforementioned findings (i)–(iii) are only weakly dependent on system size and a symmetry-energy model, and thus, they are of more general importance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Estimate of background baseline and upper limit on the chiral magnetic effect in isobar collisions at $\sqrt{S_{NN}}$=200 GeV at the BNL Relativistic Heavy Ion Collider

For the search of the chiral magnetic effect (CME), STAR previously presented the results from isobar collisions ($^{96}_{44}$Ru + $^{96}_{44}$Ru, $^{96}_{40}$Zr + $^{96}_{40}$Zr) obtained through a blind analysis. The ratio of results in Ru+Ru to Zr+Zr collisions for the CME-sensitive charge-dependent azimuthal correlator (Δ⁢𝛾), normalized by elliptic anisotropy (𝑣2), was observed to be close to but systematically larger than the inverse multiplicity ratio. The background baseline for the isobar ratio, 𝑌= (Δ⁢𝛾/𝑣2) Ru /(Δ⁢𝛾/𝑣2) Zr , is naively expected to be (1/𝑁) Ru /(1/𝑁) Zr ; however, genuine two- and three-particle correlations are expected to alter it. We estimate the contributions to 𝑌 from those correlations, utilizing both the isobar data and hijing simulations. After including those contributions, we arrive at a final background baseline for 𝑌, which is consistent with the isobar data. Here, we extract an upper limit for the CME fraction in the Δ⁢𝛾 measurement of approximately 10% at a 95% confidence level on in isobar collisions at $\sqrt{S_{NN}}$=200 GeV, with an expected 15% difference in their squared magnetic fields.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Model emulation and closure tests for (3+1)D relativistic heavy-ion collisions

In nuclear and particle physics, reconciling sophisticated simulations with experimental data is vital for understanding complex systems like the Quark Gluon Plasma (QGP) generated in heavy ion collisions. However, computational demands pose challenges, motivating using Gaussian Process emulators for efficient parameter extraction via Bayesian calibration. We conduct a comparative analysis of Gaussian Process emulators in heavy-ion physics to identify the most adept emulator for parameter extraction with minimal uncertainty. Furthermore, our study contributes to advancing computational techniques in heavy-ion physics, enhancing our ability to interpret experimental data and understand QGP properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Bayesian analysis of (3 +1)⁢D relativistic nuclear dynamics with the RHIC beam energy scan data

This work presents a Bayesian inference study for relativistic heavy-ion collisions in the beam energy scan program at the BNL Relativistic Heavy-Ion Collider. The theoretical model simulates event-by-event (3+1)-dimensional [(3+1)⁢D] collision dynamics using hydrodynamics and hadronic transport theory. We analyze the model's 20-dimensional posterior distributions obtained using three model emulators with different accuracy and demonstrate the essential role of training an accurate model emulator in the Bayesian analysis. Our analysis provides robust constraints on the quark-gluon plasma's transport properties and various aspects of (3+1)⁢D relativistic nuclear dynamics. By running full model simulations with 100 parameter sets sampled from the posterior distribution, we make predictions for p T -differential observables and estimate their systematic theory uncertainty. Here, a sensitivity analysis is performed to elucidate how individual experimental observables respond to different model parameters, providing useful physics insights into the phenomenological model for heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Polarization of massive fermions in a vortical fluid

For this study, fermions become polarized in a vortical fluid due to spin-vorticity coupling. Such a polarization can be calculated from the Wigner function in a quantum kinetic approach. By extending previous results for chiral fermions, we derive the Wigner function for massive fermions up to next-to-leading order in spatial gradient expansion. The polarization density of fermions can be calculated from the axial vector component of the Wigner function and is found to be proportional to the local vorticity ω. The polarizations per particle for fermions and antifermions decrease with the chemical potential and increase with energy (mass). Both quantities approach the asymptotic value ℏω/4 in the large energy (mass) limit. The polarization per particle for fermions is always smaller than that for antifermions, whose ratio of fermions to antifermions also decreases with the chemical potential. The polarization per particle on the Cooper-Frye freeze-out hypersurface can also be formulated and is consistent with the previous result of Becattini et al.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Exploring new small system geometries in heavy ion collisions

Relativistic heavy ion collisions produce nuclei-sized droplets of quark-gluon plasma whose expansion is well described by viscous hydrodynamic calculations. Over the past half decade, this formalism was also found to apply to smaller droplets closer to the size of individual nucleons, as produced in p + p and p + A collisions. The hydrodynamic paradigm was further tested with a variety of collision species, including p + Au,d + Au, and 3 He + Au producing droplets with different geometries. Nevertheless, questions remain regarding the importance of pre-hydrodynamic evolution and the exact medium properties during the hydrodynamic evolution phase, as well as the applicability of alternative theories that argue the agreement with hydrodynamics is accidental. In this study we explore options for new collision geometries including p + O and O + O proposed for running at the Large Hadron Collider, as well as 4 He + Au,C + Au,O + Au, and 7,9 Be + Au at the Relativistic Heavy Ion Collider.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Assessing saturation physics explanations of collectivity in small collision systems with the ip-jazma model

Experimental measurements in relativistic collisions of small systems from $p + p$ to $p/d/ ^3He + A$ at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) reveal particle emission patterns that are strikingly similar to those observed in $A + A$ collisions of large nuclei. One explanation of these patterns is the formation of small droplets of quark-gluon plasma (QGP) followed by hydrodynamic evolution. A geometry engineering program was proposed to further investigate these emission patterns, and the experimental data from that program in $p + Au, d + Au, ^3 He + Au$ collisions for elliptic and triangular anisotropy coefficients $v_2$ and $v_3$ follow the pattern predicted by hydrodynamic calculations [C. Aidala et al. (PHENIX Collaboration), Nat. Phys. 15, 214 (2019)]. One alternative approach, referred to as initial-state correlations, suggests that for small systems the patterns observed in the final-state hadrons are encoded at the earliest moments of the collision and therefore require no final-state parton scattering or hydrodynamic evolution. Recently, new calculations using only initial-state correlations, in the dilute-dense approximation of gluon saturation physics, reported striking agreement with the $v_2$ patterns observed in $p/d/ ^3He + Au$ data at RHIC [M. Mace, V. V. Skokov, P. Tribedy, and R. Venugopalan, Phys. Rev. Lett. 121, 052301 (2018)]. The results reported by Mace, Skokov, Tribedy and Venugopalan (MSTV) are counterintuitive and thus we aim here to reproduce some of the basic features of these calculations. In this first investigation, we provide a description of our publicly available model, ip-jazma, and investigate its implications for saturation scales, multiplicity distributions, and eccentricities, reserving for later work the analysis of momentum spectra and azimuthal anisotropies. We find that our implementation of the saturation physics model reproduces the results of the MSTV calculation of the multiplicity distribution in $d + Au$ collisions at RHIC. However, additional aspects of studies, together with existing data, call into question some of the essential elements reported by MSTV. Resolution of these issues will require further developments of ip-jazma, in order to determine if it can replicate the qualitative agreement with the $v_2$ reported by MSTV. Both the work reported here and future studies will establish which features in the experimental data are uniquely attributable to the color glass condensate description.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Rapidity dependence of initial state geometry and momentum correlations in p + Pb collisions

Event geometry and initial state correlations have been invoked as possible explanations of long range azimuthal correlations observed in high multiplicity p+p and p+Pb collisions. We study the rapidity dependence of initial state momentum correlations and event-by-event geometry in √s = 5.02 TeV p+Pb collisions within the 3+1D IP-Glasma model [1], where the longitudinal structure is governed by JIMWLK rapidity evolution of the incoming nuclear gluon distributions. Furthermore, we find that the event geometry is correlated across large rapidity intervals whereas initial state momentum correlations are relatively short range in rapidity. Based on our results, we discuss implications for the relevance of both effects in explaining the origin of collective phenomena in small systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Self-consistent feedback mechanism for the sudden viscous dissipation of finite-Mach-number compressing turbulence

Previous work [Davidovits and Fisch, Phys. Rev. Lett. 116, 105004 (2016)] demonstrated that the compression of a turbulent field can lead to a sudden viscous dissipation of turbulent kinetic energy (TKE), and that paper suggested this mechanism could potentially be used to design new fast-ignition schemes for inertial confinement fusion (ICF). We expand on previous work here by accounting for finite Mach numbers, rather than relying on a zero-Mach-limit assumption as previously done. The finite-Mach-number formulation is necessary to capture a self-consistent feedback mechanism in which dissipated TKE increases the temperature of the system, which in turn modifies the viscosity and thus the TKE dissipation itself. Direct numerical simulations with a tenth-order accurate Padé scheme were carried out to analyze this self-consistent feedback loop for compressing turbulence. Results show that, for finite Mach numbers, the sudden viscous dissipation of TKE still occurs, for both the solenoidal and dilatational turbulent fields. As the domain is compressed, oscillations in dilatational TKE are encountered due to the highly oscillatory nature of the pressure dilatation. An analysis of the source terms for the internal energy shows that the mechanical-work term dominates the viscous turbulent dissipation. As a result, the effect of the suddenly dissipated TKE on temperature is minimal for the Mach numbers tested. Moreover, an analytical expression is derived that confirms the dissipated TKE does not significantly alter the temperature evolution for low Mach numbers, regardless of compression speed. The self-consistent feedback mechanism is thus quite weak for subsonic turbulence, which could limit its applicability for ICF.

42 ENGINEERING↗
DOE OSTI2021

Suppression of electroconvective and morphological instabilities by an imposed cross flow of the electrolyte

Electroconvective (EC) instability and its influence on surface morphological perturbations (Morph) are important in many applications, including electrodialysis, batteries, and fuel cells. In this work, we study the effects of a two-dimensional channel flow on the EC and Morph instabilities using two approaches. In the bulk analysis, we derive the asymptotic solutions for small and large wave numbers by neglecting the space charge layer and imposing a second kind electroosmosis slip velocity boundary condition on the electroneutral bulk region. In the full analysis, the instability of the entire region of the liquid electrolyte is analyzed using the ultraspherical spectral method. Both studies show that the flow significantly affects the EC instability. The imposed flow distorts the concentration field, causing a sheltering effect which hinders the ion transport from low- to high-concentration regions, and therefore suppresses the EC instability below a certain wave number. In combination with the viscous stabilization of the high wave number modes by the space charge layer, a sufficiently strong imposed flow can fully suppress the EC instability. The increment in the critical voltage for the instability onset is roughly proportional to the square root of the product of imposed velocity and double layer thickness. The imposed flow has a smaller effect on the Morph instability, except that it may remove the morphological modes resulting from the EC instability and thereby change the wave number of the most unstable mode.

42 ENGINEERING↗
DOE OSTI2019

Azimuthal Harmonics in Small and Large Collision Systems at RHIC Top Energies

The first (v$even\atop{1}$), second (v 2 ) and third (v 3 ) harmonic coefficients of the azimuthal particle distribution at mid-rapidity, are extracted for charged hadrons and studied as a function of transverse momentum (p T ) and mean charged particle multiplicity density $\langle$N ch $\rangle$ in U+U ($\sqrt{s_ {NN}}$ = 193 GeV), Au+Au, Cu+Au, Cu+Cu, d+Au and p+Au collisions at $\sqrt{s_ {NN}}$ = 200 GeV with the STAR Detector. For the same $\langle$N ch $\rangle$, the v$even\atop{1}$ and v 3 coefficients are observed to be independent of collision system, while v 2 exhibits such a scaling only when normalized by the initial-state eccentricity (ε 2 ). The data also show that ln(v 2 /ε 2 ) scales linearly with $\langle$N ch $\rangle$ -1/3 . Finally, these measurements provide insight into initial-geometry fluctuations and the role of viscous hydrodynamic attenuation on v n from small to large collision systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Observable signatures of initial state momentum anisotropies in nuclear collisions

In this work, we show that the correlation between the elliptic momentum anisotropy v 2 and the average transverse momentum [ p T ] at fixed multiplicity in small system nuclear collisions carries information on the origin of the observed momentum anisotropy. A calculation using a hybrid IP-Glasma+ music +UrQMD model that includes contributions from final state response to the initial geometry as well as initial state momentum anisotropies of the color glass condensate predicts a characteristic sign change of the correlator ρ ^ ( v 2 2 , [ p T ] ) as a function of charged particle multiplicity in p + Au and d + Au collisions at s = 200 GeV , and p + Pb collisions at s = 5.02 TeV . This sign change is absent in calculations without initial state momentum anisotropies. The model further predicts a qualitative difference between the centrality dependence of ρ ^ ( v 2 2 , [ p T ] ) in Au + Au collisions at s = 200 GeV and Pb + Pb collisions at s = 5.02 TeV , with only the latter showing a sign change in peripheral events. Predictions for O + O collisions at different collision energy show a similar behavior. Experimental observation of these distinct qualitative features of ρ ^ ( v 2 2 , [ p T ] ) in small and large systems would constitute strong evidence for the presence and importance of initial state momentum anisotropies predicted by the color glass condensate effective theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗