Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “relativistic 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.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Radio jets in colliding galaxies: Testing the interaction-activity connection

The idea that galaxy interactions and merging are related to the generation of starburst and AGN activity in galactic nuclei has been the subject of intensive investigations over the past several years and is still a matter of lively debate. Peculiar morphologies, indicative of tidal interactions, have been detected in high-luminosity radio galaxies, in quasars, and in ultraluminous IRAS galaxies. In addition, low-luminosity radio and active galaxies show similar evidence for a recent merger or for nearby companions. In a recent CCD optical study of galaxies selected on the basis that they all contain well defined radio jets, it was found that almost half of the sample consists of pairs of elliptical galaxies. Many of these low-luminosity radio galaxies with companions show a well defined distorted radio jet structure at the VLA scale with an S- or C-shaped morphology. We have developed a general numerical simulation algorithm for ballistic radio jets with the intention of applying this model to the study of the bent jets seen in colliding pairs of galaxies and with the hope of testing the well documented interaction-activity connection. In our model the morphological evolution of the jets is determined by their response to the simple mechanical forces (i.e., gravity and ram pressure) imposed on them from both the host and the companion galaxies. Radiative losses, jet precession, magnetic effects, relativistic terms, and hydrodynamic instabilities have all been ignored. Starting with a previously derived collision model for the interacting pair of elliptical galaxies NGC 4782/4783, we have used our algorithm to simulate the specific two-sided jet morphology seen in the radio source 3C 278, associated with NGC 4782. This is the first time that such jet simulations have been produced for a galaxy pair whose relative orbit was determined independently from the jet modeling.

Borne, Kirk D.↗
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.↗
DOE OSTI2021

Modified equilibrium distributions for Cooper-Frye particlization

We introduce a positive-definite single-particle distribution that is suitable for describing the transition from a macroscopic hydrodynamic to a microscopic kinetic description during the late stages of heavy-ion collisions in the presence of moderately large viscous corrections. The modified equilibrium distribution function can be constructed with hydrodynamic input from either relativistic viscous fluid dynamics or anisotropic fluid dynamics. We test the modified equilibrium distribution's hydrodynamic output for a stationary hadron resonance gas subject to either shear stress, bulk pressure, or baryon diffusion current at a given freeze-out temperature and baryon chemical potential. While it does not reproduce all components of the net baryon current and energy-momentum tensor exactly, it significantly improves upon the customary linearized approximations for the nonequilibrium correction δf n which typically lead to unphysical negative distribution functions at large particle momenta. A comparison of particle spectra and p T –differential elliptic flow coefficients from the Cooper–Frye formula computed with the modified equilibrium distribution and with linearized δf n corrections is presented, for two different (2 + 1)–dimensional hypersurfaces corresponding to central and noncentral Pb + Pb collisions at the Large Hadron Collider (LHC).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Studying QGP with flow: A theory overview

I review recent developments in the phenomenological study of the quark-gluon plasma (QGP) transport properties based on a personal selection of results that were presented at Quark Matter 2019. Additionally, the constraints on the temperature dependence of QGP shear and bulk viscosity are summarized. I discuss new theory advancements towards more realistic 3D dynamical simulations of heavy-ion collisions at finite baryon density. The challenges and opportunities of applying hydrodynamics to small collision systems are highlighted.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Bayesian calibration of viscous anisotropic hydrodynamic simulations of heavy-ion collisions

Owing to large pressure gradients at early times, standard hydrodynamic model simulations of relativistic heavy-ion collisions do not become reliable until O(1) fm/c after the collision. To address this one often introduces a prehydrodynamic stage that models the early evolution microscopically, typically as a conformal, weakly interacting gas. In such an approach the transition from the prehydrodynamic to the hydrodynamic stage is discontinuous, introducing considerable theoretical model ambiguity. Alternatively, fluids with large anisotropic pressure gradients can be handled macroscopically using the recently developed viscous anisotropic hydrodynamics (VAH). In high-energy heavy-ion collisions VAH is applicable already at very early times, and at later times transitions smoothly into conventional second-order viscous hydrodynamics. Here we present a Bayesian calibration of the VAH model with experimental data for Pb-Pb collisions at the LHC at $\sqrt{^{S}NN}$ = 2.76 A TeV. We find that the VAH model has the unique capability of constraining the specific viscosities of the quark-gluon plasma at higher temperatures than other previously used models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Modeling stochastic fluctuations in relativistic kinetic theory

Using the information current, we develop a Lorentz-covariant framework for modeling equilibrium fluctuations in relativistic kinetic theory in the grand-canonical ensemble. The resulting stochastic theory is proven to be causal and covariantly stable, and its predictions do not depend on the choice of spacetime foliation used to define the grand-canonical probabilities. As expected, in a box containing N > 5 particles, Boltzmann’s molecular chaos postulate is broken with (almost exact) probability N -1/2 , leading to a breakdown of the Boltzmann equation in small systems. Here, we also verify that, in ultrarelativistic gases, transient hydrodynamics already accounts for at least 80% of the equilibrium fluctuations of the stress-energy tensor at a given time. Finally, we compute the correlators at nonequal times for two selected collision kernels: that of a chemically active diluted solution, and that of ultrarelativistic scalar particles self-interacting via a quartic potential. For the former, we compute the density-density correlators analytically in real space, and dehydrodynamization of the stochastic theory is proven to occur whenever the mean free path diverges at high energy.

Astronomy & Astrophysics↗
DOE OSTI2023

Colliding localized, lumpy holographic shocks with a granular nuclear structure

We apply a recent and simple technique which speeds up the calculation of localized collisions in holography to study more realistic models of the pre-hydrodynamic phase of heavy ion collisions using gauge/gravity duality. Our initial data reflects the lumpy nuclear structure of real heavy ions and our projectiles’ aspect ratio mimics the Lorentz contraction of nuclei during RHIC collisions. At the hydrodynamization time of the central region of the quark gluon plasma developed during the collision, we find that most of the system’s vorticity is located well outside the hydrodynamized part of the plasma. Only the relativistic corrections to the thermal vorticity within the hydrodynamized region are non-negligible. We compare the transverse flow shortly after the collision with previous results which did not use granular initial conditions and determine the proper energy density and fluid velocity in the hydrodynamized subregion of the plasma.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Post-Newtonian kinetic theory

A kinetic theory for relativistic gases in the presence of gravitational fields is developed in the second post-Newtonian approximation. The corresponding Boltzmann equation is determined from the evolution of the one-particle distribution function with respect to the proper time along the world line of the particle. From the knowledge of the equilibrium Maxwell–Jüttner distribution function in the second post-Newtonian approximation the components of the particle four-flow and energy–momentum tensor are obtained. The Eulerian hydrodynamic equations for the mass density, mass-energy density and momentum density in the second post-Newtonian approximation are determined from the Boltzmann equation. It is shown that the combination of the hydrodynamic equations of mass and mass-energy densities leads to the hydrodynamic equation for the internal energy density in the first post-Newtonian approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Relaxation time approximation for a multispecies relativistic gas

We generalize a recent prescription for the relaxation time approximation for the relativistic Boltzmann equation for systems with multiple particle species at finite temperature. This is performed by adding counter-terms to the traditional Anderson-Witting ansatz for each particle species. Our approach allows for the use of momentum-dependent relaxation times and the obedience of local conservation laws regardless of the definition of the local equilibrium state. As an application, we derive the first order Chapman-Enskog corrections to the equilibrium distribution and display results for the hadron-resonance gas. We also demonstrate that our collision term ansatz obeys the second law of thermodynamics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2001

Jet Stability and the Generation of Superluminal and Stationary Components

We present a numerical simulation of the response of an expanding relativistic jet to the ejection of a superluminal component. The simulation has been performed with a relativistic time-dependent hydrodynamical code from which simulated radio maps are computed by integrating the transfer equations for synchrotron radiation. The interaction of the superluminal component with the underlying jet results in the formation of multiple conical shocks behind the main perturbation. These trailing components can be easily distinguished because they appear to be released from the primary superluminal component instead of being ejected from the core. Their oblique nature should also result in distinct polarization properties. Those appearing closer to the core show small apparent motions and a very slow secular decrease in brightness and could be identified as stationary components. Those appearing farther downstream are weaker and can reach superluminal apparent motions. The existence of these trailing components indicates that not all observed components necessarily represent major perturbations at the jet inlet; rather, multiple emission components can be generated by a single disturbance in the jet. While the superluminal component associated with the primary perturbation exhibits a rather stable pattern speed, trailing components have velocities that increase with distance from the core but move at less than the jet speed. The trailing components exhibit motion and structure consistent with the triggering of pinch modes by the superluminal component. The increase in velocity of the trailing components is an indirect consequence of the acceleration of the expanding fluid, which is assumed to be relativistically hot; if observed, such accelerations would therefore favor an electron-positron (as opposed to proton rest mass) dominated jet.

Agudo, Ivan↗
DOE OSTI2022

Pseudogauge dependence of the spin polarization and of the axial vortical effect

The mean spin polarization vector of spin 1/2 particles in a relativistic fluid at local thermal equilibrium (LTE) in different pseudogauges (PGs), i.e., with different choices for the decomposition of orbital and spin angular momentum, is obtained. The spin polarization obtained in the canonical PG differs from the one obtained in the Belinfante PG. It is found that this difference can not be written by replacing the thermal vorticity with the spin potential in the usual polarization formula. Other PG choices affect the contribution of thermal shear to the spin polarization. So, the choice of a PG affects the predictions for the polarization measured in heavy-ion collisions. In general, it is shown that the Wigner function of a noninteracting Dirac field at LTE is affected by the PG transformations. Explicit expressions of the mean axial current are also calculated and it is found that even the axial vortical effect conductivity depends on the PG.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Global polarization of ΛΛ¯ hyperons in Pb-Pb collisions at √ s NN = 2.76 and 5.02 TeV

The global polarization of the Λ and Λ ¯ hyperons is measured for Pb-Pb collisions at s N N = 2.76 and 5.02 TeV recorded with the ALICE at the Large Hadron Collider (LHC). The results are reported differentially as a function of collision centrality and hyperon's transverse momentum ( p T ) for the range of centrality 5–50%, 0.5 < p T < 5 GeV/ c , and rapidity | y | < 0.5 . The hyperon global polarization averaged for Pb-Pb collisions at s N N = 2.76 and 5.02 TeV is found to be consistent with zero, ( P H ) ( % ) ≈ 0.01 ± 0.06 (stat.) ± 0.03 (syst.) in the collision centrality range 15–50%, where the largest signal is expected. The results are compatible with expectations based on an extrapolation from measurements at lower collision energies at the Relativistic Heavy Ion Collider, hydrodynamical model calculations, and empirical estimates based on collision energy dependence of directed flow, all of which predict the global polarization values at LHC energies of the order of 0.01 % .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Kinematic and dynamical origins of mean-p T fluctuations in heavy-ion collisions

Event-by-event fluctuations of the mean transverse momentum (mean-p T ) provide a sensitive probe of collective dynamics beyond single-particle spectra and anisotropic flow. We present a systematic study of mean-p T fluctuation observables using a Bayesian-calibrated multistage hydrodynamic framework, including quantitative comparisons to RHIC measurements and model-based investigations of beam-energy and kinematic-acceptance effects. The experimental definitions employed by the STAR and ALICE Collaborations are implemented explicitly and found to yield consistent results within controlled limits. We study the centrality and beam-energy dependence of the observable, its sensitivity to key soft-sector ingredients, and the impact of the kinematic p T acceptance. By introducing scaled-p T cuts, we demonstrate that a part of the apparent energy dependence arises from kinematic projection effects, while the remaining trends reflect genuine collective dynamics. Our results establish mean-p T fluctuations as a nontrivial and independent validation of calibrated hydrodynamic descriptions of the quark–gluon plasma.

Event-by-event correlations↗
DOE OSTI2024

Secular Outflows from Long-Lived Neutron Star Merger Remnants

We study mass ejection from a binary neutron star merger producing a long-lived massive neutron star remnant with general-relativistic neutrino-radiation hydrodynamics simulations. In addition to outflows generated by shocks and tidal torques during and shortly after the merger, we observe the appearance of a wind driven by spiral density waves in the disk. This spiral-wave-driven outflow is predominantly located close to the disk orbital plane and have a broad distribution of electron fractions. At higher latitudes, a high electron-fraction wind is driven by neutrino radiation. The combined nucleosynthesis yields from all the ejecta components is in good agreement with Solar abundance measurements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Long-lived neutron-star remnants from asymmetric binary neutron star mergers: element formation, kilonova signals and gravitational waves

We present 3D general-relativistic neutrino-radiation hydrodynamics simulations of two asymmetric binary neutron star mergers producing long-lived neutron stars remnants and spanning a fraction of their cooling time scale. The mergers are characterized by significant tidal disruption with neutron rich material forming a massive disc around the remnant. The latter develops one-armed dynamics that is imprinted in the emitted kilo-Hertz gravitational waves. Angular momentum transport to the disc is initially driven by spiral-density waves and enhanced by turbulent viscosity and neutrino heating on longer timescales. The mass outflows are composed by neutron-rich dynamical ejecta of mass ∼10 −3 –10 −2 M ⊙ followed by a persistent spiral-wave/neutrino-driven wind of ≳ 10 −2 M ⊙ with material spanning a wide range of electron fractions, ∼0.1–0.55. Dynamical ejecta (winds) have fast velocity tails up to ∼0.8 (∼0.4) c. The outflows are further evolved to days timescale using 2D ray-by-ray radiation-hydrodynamics simulations that include an online nuclear network. We find complete r-process yields and identify the production of 56 Ni and the subsequent decay chain to 56 Co and 56 Fe. Synthetic kilonova light curves predict an extended (near-) infrared peak a few days postmerger originating from r-process in the neutron-rich/high-opacity ejecta and UV/optical peaks at a few hours (ten minutes) postmerger originating from weak r-process (free-neutron decay) in the faster ejecta components. Additionally, the fast tail of tidal origin generates kilonova afterglows potentially detectable in radio and X band on a few to ten years time scale. Quantitative effects originating from the tidal disruption merger dynamics are reflected in the multimessenger emissions.

abundances↗
DOE OSTI2020

Improved single-particle phase-space distributions for viscous fluid dynamic models of relativistic heavy ion collisions

Comparison of hydrodynamic calculations with experimental data inevitably requires a model for converting the fluid to particles. For shear viscous fluids, several different conversion models have been proposed, such as the quadratic Grad corrections, the Strickland-Romatschke (SR) ansatz, self-consistent shear corrections from linearized kinetic theory, or corrections from the relaxation time approach. We show, via comparison to nonlinear 2 → 2 kinetic theory evolution, that even for low specific shear viscosities 0.03 ≲ $\textit{η/s}$ ≲ 0.2 , these four models have varying accuracy and still significant errors in reproducing particle phase space densities from hydrodynamic fields. Moreover, we demonstrate that the overall reconstruction error of additive shear viscous $f = f_{eq} + δf$ models dramatically reduces by up to one order of magnitude, if one ensures through exponentiation that $\textit{f}$ is always positive. Further, we also illustrate how even more accurate viscous $\textit{δf}$ models can be constructed by incorporating the first time derivatives of hydrodynamic fields to account for the past evolution of the system. Such time derivatives are readily available in hydrodynamic simulations, though usually not included in the output. Though our comparisons are limited to a one-component massless system undergoing a (0 + 1)-dimensional (0 + 1D) longitudinal boost-invariant expansion, we expect that the improved models will be useful in 2 + 1D and 3 + 1D hydrodynamic studies as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗