Search NASA⌕ Search

SEARCH · Search NASA

Results for “fluctuating hydrodynamics”

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.

268 records · Page 15
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 OSTI2023

Assessing the ultracentral flow puzzle in hydrodynamic modeling of heavy-ion collisions

An outstanding problem in heavy-ion collisions is the inability for models to accurately describe ultra-central experimental flow data, despite that being precisely the regime where a hydrodynamic description should be most applicable. Here, we reassess the status of this puzzle by computing the flow in ultra-central collisions obtained from multiple recent Bayesian models that were tuned to various observables in different collision systems at typical centralities. While central data can now be described with better accuracy than in previous calculations, tension with experimental observation remains and worsens as one goes to ultra-central collisions. Tuning the model parameters cannot remove this tension without destroying the fit at other centralities. As such, new elements are likely needed in the standard modeling of heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1993

Steady state and dynamical structure of a cosmic-ray-modified termination shock

A hydrodynamic model is developed for the structure of a cosmic-ray-modified termination shock. The model is based on the two-fluid equations of diffuse shock acceleration (Drury and Volk, 1981). Both the steady state structure of the shock and its interaction with outer heliospheric disturbances are considered. Under the assumption that the solar wind is decelerated by diffusing interstellar cosmic rates, it is shown that the natural state of the termination shock is a gradual deceleration and compression, followed by a discontinuous jump to a downstream state which is dominated by the pressure contribution of the cosmic rays. A representative model is calculated for the steady state which incorporates both interstellar cosmic ray mediation and diffusively accelerated anomalous ions through a proposed thermal leakage mechanism. The interaction of large-scale disturbances with the equilibrium termination shock model is shown to result in some unusual downstream structure, including transmitted shocks and cosmic-ray-modified contact discontinuities. The structure observed may be connected to the 2-kHz outer heliospheric radio emission (Cairns et al., 1992a, b). The time-dependent simulations also demonstrate that interaction with solar wind compressible turbulence (e.g., traveling interplanetary shocks, etc.) could induce the termination shock to continually fluctuate between cosmic-ray-dominated and gas-dynamic states. This fluctuation may represent a partial explanation of the galactic cosmic ray modulation effect and illustrates that the Pioneer and Voyager satellites will encounter an evolving shock whose structure and dynamic properties are strongly influence by the mediation of interstellar and anomalous cosmic rays.

Donohue, D. J.↗
NASA NTRS2012

X-Ray Spectra from MHD Simulations of Accreting Black Holes

We present the results of a new global radiation transport code coupled to a general relativistic magneto-hydrodynamic simulation of an accreting, nonrotating black hole. For the first time, we are able to explain from first principles in a self-consistent way the X-ray spectra observed from stellar-mass black holes, including a thermal peak, Compton reflection hump, power-law tail, and broad iron line. Varying only the mass accretion rate, we are able to reproduce the low/hard, steep power-law, and thermal-dominant states seen in most galactic black hole sources. The temperature in the corona is T(sub e) ~ 10 keV in a boundary layer near the disk and rises smoothly to T(sub e) greater than or approximately 100 keV in low-density regions far above the disk. Even as the disk's reflection edge varies from the horizon out to approximately equal to 6M as the accretion rate decreases, we find that the shape of the Fe Ka line is remarkably constant. This is because photons emitted from the plunging region are strongly beamed into the horizon and never reach the observer. We have also carried out a basic timing analysis of the spectra and find that the fractional variability increases with photon energy and viewer inclination angle, consistent with the coronal hot spot model for X-ray fluctuations.

Schnittman, Jeremy D.↗
NASA NTRS1994

A hydrodynamic approach to cosmology: The mixed dark matter cosmological scenario

We compute the evolution of spatially flat, mixed cold and hot dark matter models containing both baryonic matter and two kinds of dark matter. Hydrodynamics is treated with a highly developed Eulerian hydrodynamic code (see Cen 1992). A standard particle-mesh (PM) code is also used in parallel to calculate the motion of the dark matter components. We adopt the following parameters: h equivalent to (sub 0)/100 km/s Mpc(exp -1) = 0.5, OMEGA(sub C) = 0.3, and OMEGA(sub B) = 0.06, with amplitude of the perturbation spectrum fixed by the Cosmic Background Explorer Satellite (COBE) Dark Matter Radiation (DMR) measurements (Smoot et al. 1992) being sigma (sub 8) = 0.67. Four different boxes are simulated with box sizes of L = (64, 16, 4, 1) h(exp -1) Mpc, respectively, the two small boxes providing good resolution but little valid information due to absence of large-scale power. We use 128(exp 3) approximate 10(exp 6.3) baryonic cells, 128(exp .3) cold dark matter particles, and 2 x 128(exp 3) hot dark matter particles. In addition to the dark matter we follow separately six baryonic species (H, H(+), He, He(+), He(++), e(-)) with allowance for both (nonequilibrium) collisional and radiative ionization in every cell. The background radiation field is also followed in detail with allowance made for both continuum and line processes, to allow nonequilibrium heating and cooling processes to be followed in detail. The mean final Zeldovich-Sunyaev y parameter is estimated to be y Bar = (5.4 + or - 2.7) x 10(exp -7) below currently attainable observations, with a rms fluctuation of approximately delta bar y = (0.6 + or - 3.0) x 10(exp -7) on arcminute scales. The rate of galaxy formation peaks at an even later epoch (z approximate 0.3) than in the standard (OMEGA = 1, sigma sub 8 = 0.67) cold dark matter (CDM) model (z approximate 0.5) and, at a redshift of z = 4, is nearly a factor of 100 lower than for the CDM model with the same value of sigma sub 8. With regard to mass function, the smallest objects are stabilized against collapse by thermal energy: the mass-weighted mass spectrum has a broad peak in the vicinity of M(sub B) = 10(exp 9.5) solar mass with a reasonable fit to the Schechter luminosity function if the ratio of baryon mass to blue light is approximately 4. In addition, one very large PM simulation was made in a box with size (320 h(exp - 1) Mpc) containing 3 x 200(exp 3) = 10(exp 7.4) particles. Utilizing this simulation we find that the model yields a cluster mass function which is about a factor of 4 higher than observed, but a cluster-cluster correlation length marginally lower than observed, but that both are closer to observations than in the (COBE) normalized CDM model. The one-dimensional pairwise velocity dispersion is 605 + or - 8 km/s at 1/h separation, lower than that of the DCM model normalized to COBE, but still significant higher than observations (Davis & Peebles 1983). A plausible velocity bias b(sub v) = 0.8 + or - 0.1 on this scale will reduce but not remove the discrepancy. The velocity auto-correlat ion function has a coherence length of 40/h Mpc, which is somewhat lower than the observed counterpart. In all these respects the model would be improved by decreasing the cold fraction of the dark OMEGA(sub CDM)/ (OMEGA(sub CDM) + OMEGA(sub HDB). But formation of galaxies and clusters of galaxies is much later in this model than in COBE-normalized CDM, perhaps too late. To improve on these constraints a larger ratio of OMEGA(sub CDM)/ (OMEGA(sub CDM) + OMEGA(sub HDM)) is required than the value of 0.67 adopted here. It does not seem possible to find a value for this ratio which would satisfy all tests. Overall, the model is similar both on large and intermediate scales to the standard CDM model normalized to the same value of sigma(sub B), but the problem with regard to late formation of galaxies is more severe in this model than in that CDM model. Adding hot dark matter, significantly improves the ability of the COBE-normalized CDM scenario to fit existing observations, but the model is in fact not as good as the CDM model with the same sigma(sub 8) and is still probably unsatisfactory with regard to several critical tests.

Cen, Renyue↗
DOE OSTI2022

Energy conserving and well-balanced discontinuous Galerkin methods for the Euler–Poisson equations in spherical symmetry

ABSTRACT This paper presents high-order Runge–Kutta (RK) discontinuous Galerkin methods for the Euler–Poisson equations in spherical symmetry. The scheme can preserve a general polytropic equilibrium state and achieve total energy conservation up to machine precision with carefully designed spatial and temporal discretizations. To achieve the well-balanced property, the numerical solutions are decomposed into equilibrium and fluctuation components that are treated differently in the source term approximation. One non-trivial challenge encountered in the procedure is the complexity of the equilibrium state, which is governed by the Lane–Emden equation. For total energy conservation, we present second- and third-order RK time discretization, where different source term approximations are introduced in each stage of the RK method to ensure the conservation of total energy. A carefully designed slope limiter for spherical symmetry is also introduced to eliminate oscillations near discontinuities while maintaining the well-balanced and total-energy-conserving properties. Extensive numerical examples – including a toy model of stellar core collapse with a phenomenological equation of state that results in core bounce and shock formation – are provided to demonstrate the desired properties of the proposed methods, including the well-balanced property, high-order accuracy, shock-capturing capability, and total energy conservation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Underground Hydrogen Storage: Transforming Subsurface Science into Sustainable Energy Solutions

As the global economy moves toward net-zero carbon emissions, large-scale energy storage becomes essential to tackle the seasonal nature of renewable sources. Underground hydrogen storage (UHS) offers a feasible solution by allowing surplus renewable energy to be transformed into hydrogen and stored in deep geological formations such as aquifers, salt caverns, or depleted reservoirs, making it available for use on demand. This study thoroughly evaluates UHS concepts, procedures, and challenges. This paper analyzes the most recent breakthroughs in UHS technology and identifies special conditions needed for its successful application, including site selection guidelines, technical and geological factors, and the significance of storage characteristics. The integrity of wells and caprock, which is important for safe and efficient storage, can be affected by the operating dynamics of the hydrogen cycle, notably the fluctuations in pressure and stress within storage formations. To evaluate its potential for broader adoption, we also examined economic elements such as cost-effectiveness and the technical practicality of large-scale storage. We also reviewed current UHS efforts and identified key knowledge gaps, primarily in the areas of hydrogen–rock interactions, geochemistry, gas migration control, microbial activities, and geomechanical stability. Resolving these technological challenges, regulatory frameworks, and environmental sustainability are essential to UHS’s long-term and extensive integration into the energy industry. This article provides a roadmap for UHS research and development, emphasizing the need for further research to fully realize the technology’s promise as a pillar of the hydrogen economy.

25 ENERGY STORAGE↗
DOE OSTI2025

Theory of x-ray photon correlation spectroscopy for multiscale flows

Complex multiscale flows associated with instabilities and turbulence are commonly induced under high-energy density (HED) conditions, but accurate measurement of their transport properties has been challenging. X-ray photon correlation spectroscopy (XPCS) with coherent x-ray sources can, in principle, probe material dynamics to infer transport properties using time autocorrelation of density fluctuations. Here we develop a theoretical framework for utilizing XPCS to study material diffusivity in multiscale flows. We extend single-scale shear flow theories to broadband flows using a multiscale analysis that captures shear and diffusion dynamics. Our theory is validated with simulated XPCS for Brownian particles advected in multiscale flows. We demonstrate the versatility of the method over several orders of magnitude in timescale using sequential-pulse XPCS, single-pulse x-ray speckle visibility spectroscopy (XSVS), and double-pulse XSVS.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

A versatile pressure-cell design for studying ultrafast molecular-dynamics in supercritical fluids using coherent multi-pulse x-ray scattering

Supercritical fluids (SCFs) can be found in a variety of environmental and industrial processes. They exhibit an anomalous thermodynamic behavior, which originates from their fluctuating heterogeneous micro-structure. Characterizing the dynamics of these fluids at high temperature and high pressure with nanometer spatial and picosecond temporal resolution has been very challenging. The advent of hard x-ray free electron lasers has enabled the development of novel multi-pulse ultrafast x-ray scattering techniques, such as x-ray photon correlation spectroscopy (XPCS) and x-ray pump x-ray probe (XPXP). These techniques offer new opportunities for resolving the ultrafast microscopic behavior in SCFs at unprecedented spatiotemporal resolution, unraveling the dynamics of their micro-structure. However, harnessing these capabilities requires a bespoke high-pressure and high-temperature sample system that is optimized to maximize signal intensity and address instrument-specific challenges, such as drift in beamline components, x-ray scattering background, and multi-x-ray-beam overlap. We present a pressure cell compatible with a wide range of SCFs with built-in optical access for XPCS and XPXP and discuss critical aspects of the pressure cell design, with a particular focus on the design optimization for XPCS.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Centrality and transverse-momentum dependence of higher-order flow harmonics of identified hadrons in Au+Au collisions at s N N = 200 GeV

Here we present high-precision measurements of elliptic, triangular, and quadrangular flow v 2 , v 3 , and v 4 , respectively, at midrapidity for identified hadrons π, p, K, φ, K s , Λ as a function of centrality and transverse momentum in Au+Au collisions at the center-of-mass energy √ s NN = 200 GeV. We observe similar v n trends between light and strange mesons which indicates that the heavier strange quarks flow as strongly as the lighter up and down quarks. The number-of-constituent-quark scaling for v 2 , v 3 , and v 4 is found to hold within statistical uncertainty for 0–10%, 10–40%, and 40–80% collision centrality intervals. The results are compared to several viscous hydrodynamic calculations with varying initial conditions, and could serve as an additional constraint to the development of hydrodynamic models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Study of the Bose-Einstein correlations of same-sign pions in proton-lead collisions

Correlations of same-sign charged pions are analysed using proton-lead collision data collected by the LHCb experiment at a nucleon-nucleon centre-of-mass energy of 5.02 TeV, corresponding to an integrated luminosity of 1.06 nb -1 . Bose-Einstein correlations are observed in the form of an enhancement of pair production for same-sign charged pions with a small four-momentum difference. The dependence of the correlation radius and the intercept parameter on the reconstructed charged-particle multiplicity is investigated. The measured correlation radii scale linearly with the cube root of the reconstructed charged-particle multiplicity, being compatible with predictions of hydrodynamic models on the collision system evolution.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Measurements of long-range two-particle correlation over a wide pseudorapidity range in p–Pb collisions at $ \sqrt{s_{\textrm{NN}}}$ = 5.02 TeV

Correlations in azimuthal angle extending over a long range in pseudorapidity between particles, usually called the “ridge” phenomenon, were discovered in heavy-ion collisions, and later found in pp and p–Pb collisions. In large systems, they are thought to arise from the expansion (collective flow) of the produced particles. Extending these measurements over a wider range in pseudorapidity and final-state particle multiplicity is important to understand better the origin of these long-range correlations in small collision systems. In this Letter, measurements of the long-range correlations in p–Pb collisions at $\sqrt{^SNN}$ = 5.02 TeV are extended to a pseudorapidity gap of Δη ~ 8 between particles using the ALICE forward multiplicity detectors. After suppressing non-flow correlations, e.g., from jet and resonance decays, the ridge structure is observed to persist up to a very large gap of Δη ~ 8 for the first time in p–Pb collisions. This shows that the collective flow-like correlations extend over an extensive pseudorapidity range also in small collision systems such as p–Pb collisions. The pseudorapidity dependence of the second-order anisotropic flow coefficient, v 2 (η), is extracted from the long-range correlations. The v 2 (η) results are presented for a wide pseudorapidity range of –3.1 < η < 4.8 in various centrality classes in p–Pb collisions. To gain a comprehensive understanding of the source of anisotropic flow in small collision systems, the v 2 (η) measurements are compared with hydrodynamic and transport model calculations. The comparison suggests that the final-state interactions play a dominant role in developing the anisotropic flow in small collision systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Shape of atomic nuclei in heavy ion collisions

In the hydrodynamic model description of heavy ion collisions, the final-state anisotropic flows v n are linearly related to the strengths of the multipole shape of the distribution of nucleons in the transverse plane, ε n : v n ∝ε n . The ε n , for n=1, 2, 3, 4, are sensitive to the shapes of the colliding ions, characterized by the quadrupole β 2 , octupole β 3 , and hexadecapole β 4 deformations. This sensitivity is investigated analytically and also in a Monte Carlo Glauber model. One observes a robust linear relation, $\langleε^{2}_{n}\rangle$ = $a^{'}_{n} + b^{'}_{n}β^{2}_{n}$, for events in a fixed centrality. The $\langleε^{2}_{1}\rangle$ has a contribution from β 3 and β 4 , and $\langleε^{2}_{3}\rangle$ from β 4 . In ultracentral collisions, there are little cross contributions between β 2 and ε 3 and between β 3 and ε 2 , but clear cross contributions are present in noncentral collisions. Additionally, $\langleε^{2}_{n}\rangle$ are insensitive to nonaxial shape parameters such as the triaxiality. This is good news because the measurements of v 2 , v 3 , and v 4 can be used to constrain simultaneously the β 2 , β 3 , and β 4 values. This is best done by comparing two colliding ions with similar mass numbers and therefore nearly identical $a^{'}_{n}$, to obtain a simple equation that relates the β n of the two species. Finally, this opens up the possibility to map the shape of the atomic nuclei at a timescale (<10 –24 s) much shorter than probed by low-energy nuclear structure physics (<10 –21 s), which ultimately may provide information complementary to that obtained in the nuclear structure experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Global Λ-hyperon polarization in Au + Au collisions at $\sqrt{S_{NN}}$ = 3 GeV

Global hyperon polarization, $\overline{P}$ H , in Au + Au collisions over a large range of collision energy, $\sqrt{S_{NN}}$, was recently measured and successfully reproduced by hydrodynamic and transport models with intense fluid vorticity of the quark-gluon plasma. While na¨ıve extrapolation of data trends suggests a large $\overline{P}$ H as the collision energy is reduced, the behavior of $\overline{P}$ H at small $\sqrt{S_{NN}}$ < 7.7 GeV is unknown. Operating the STAR experiment in fixed-target mode, we measured the polarization of Λ hyperons along the direction of global angular momentum in Au + Au collisions at $\sqrt{S_{NN}}$ = 3 GeV. As a result, the observation of substantial polarization of 4.91 ± 0.81 (stat.) ± 0.15 (syst.) % in these collisions may require a reexamination of the viscosity of any fluid created in the collision, of the thermalization timescale of rotational modes, and of hadronic mechanisms to produce global polarization.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Two-particle Bose-Einstein correlations and their Lévy parameters in PbPb collisions at s N N = 5.02 TeV

Two-particle Bose–Einstein momentum correlation functions are studied for charged-hadron pairs in lead-lead collisions at a center-of-mass energy per nucleon pair of s NN = 5.02 TeV . The data sample, containing 4.27 × 10 9 minimum bias events corresponding to an integrated luminosity of 0.607 nb − 1 , was collected by the CMS experiment in 2018. The experimental results are discussed in terms of a Lévy-type source distribution. The parameters of this distribution are extracted as functions of particle pair average transverse mass and collision centrality. These parameters include the Lévy index or shape parameter α , the Lévy scale parameter R , and the correlation strength parameter λ . The source shape, characterized by α , is found to be neither Cauchy nor Gaussian, implying the need for a full Lévy analysis. Similarly to what was previously found for systems characterized by Gaussian source radii, a hydrodynamical scaling is observed for the Lévy R parameter. The λ parameter is studied in terms of the core-halo model. ©2024 CERN, for the CMS Collaboration 2024 CERN

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

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↗