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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Stochastic relativistic viscous hydrodynamics from the Metropolis algorithm

We propose an algorithm for simulating stochastic relativistic fluid dynamics based on Metropolis updates. Each step of the algorithm begins with an update based on ideal hydrodynamics. This is followed by proposing random (spatial) momentum transfers between fluid cells, keeping the total energy fixed. These proposals are then accepted or rejected using the change in entropy as a statistical weight. The algorithm reproduces relativistic viscous hydrodynamics in the “density frame,” which is a formulation of viscous hydrodynamics we review and clarify here. This formulation is first order in time and requires no auxiliary dynamical fields such as Π 𝜇⁢𝜈 . The only parameters are the shear and bulk viscosities and the equation of state. Here, by adopting the 3+1 split of general relativity, we extend the Metropolis algorithm to general space-time coordinates, such as Bjorken coordinates, which are commonly used to simulate heavy-ion collisions.

Hydrodynamic noise↗
DOE OSTI2023

Collective dynamics of polarized spin-half fermions in relativistic heavy-ion collisions

Standard relativistic hydrodynamics has been successful in describing the properties of the strongly interacting matter produced in the heavy-ion collision experiments. Recently, there has been a significant theoretical advancement in this field to explain spin polarization of hadrons emitted in these processes. Although current models have successfully explained some of the experimental data based on the coupling between spin polarization and vorticity of the medium, they still lack a clear understanding of the differential measurements. This is commonly interpreted as an indication that the spin needs to be treated as an independent degree of freedom whose dynamics is not entirely bound to flow circulation. In particular, if the spin is a macroscopic property of the system, in equilibrium its dynamics should follow hydrodynamic laws. Here, we develop a framework of relativistic hydrodynamics which includes spin degrees of freedom from the quantum kinetic theory for Dirac fermions and use it for modeling the dynamics of matter. Following experimental observations, we assume that the polarization effects are small and derive conservation laws for the net baryon current, the energy–momentum tensor and the spin tensor based on the de Groot–van Leeuwen–van Weert definitions of these currents. We present various properties of the spin polarization tensor and its components, analyze the propagation properties of the spin polarization components, and derive the spin-wave velocity for arbitrary statistics. We find that only the transverse spin components propagate, analogously to the electromagnetic waves. Finally, using our framework, we study the space–time evolution of the spin polarization for the systems respecting certain space–time symmetries and calculate the mean spin polarization per particle, which can be compared to the experimental data. We find that, for some observables, our spin polarization results agree qualitatively with the experimental findings and other model calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Relativistic spin hydrodynamics with torsion and linear response theory for spin relaxation

Using the second law of local thermodynamics and the first-order Palatini formalism, we formulate relativistic spin hydrodynamics for quantum field theories with Dirac fermions, such as QED and QCD, in a torsionful curved background. We work in a regime where spin density, which is assumed to relax much slower than other non-hydrodynamic modes, is treated as an independent degree of freedom in an extended hydrodynamic description. Spin hydrodynamics in our approach contains only three non-hydrodynamic modes corresponding to a spin vector, whose relaxation time is controlled by a new transport coefficient: the rotational viscosity. We study linear response theory and observe an interesting mode mixing phenomenon between the transverse shear and the spin density modes. We propose several field-theoretical ways to compute the spin relaxation time and the rotational viscosity, via the Green-Kubo formula based on retarded correlation functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Acausality-driven instabilities in relativistic viscous hydrodynamics

We investigate non-linear instabilities stemming from superluminal propagation of information in Israel-Stewart-like models of relativistic viscous fluid dynamics. In relativity, the characteristic speed of propagation of information, $w$, and the speed of the fluid, $v$, allow us to differentiate between regimes of the hydrodynamic equations that are acausal but stable ($w>1$), unstable ($v^{2} w^{2} \geq 1$), and covariantly ill-posed ($w^{2} \leq 0$). As an analytical benchmark, we present a new solution that illustrates these distinct regimes. We compare this analytical solution to the result of a numerical relativistic viscous fluid dynamics solver, and confirm that the analytical result can be recovered numerically in the stable regime, whether causal or acausal. The onset of numerical instabilities is further found to occur in the regime predicted by the analytical solution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Fluctuating relativistic dissipative hydrodynamics as a gauge theory

Here we argue that different formulations of hydrodynamics are related to uncertainties in the definitions of local thermodynamic and hydrodynamic variables. We show that this ambiguity can be resolved by viewing different formulations of hydrodynamics as particular “gauge choices” which lead to the same physical behavior of the system. Using the example of bulk viscosity, we show that Bemfica–Disconzi–Noronha–Kovtun (BDNK) and Israel–Stewart hydrodynamics are particular “gauge choices” of this type, related by a well-defined transformation of thermodynamic and hydrodynamic variables. We argue that this gauge ambiguity is necessary to ascertain the causality of stochastic hydrodynamic evolution and conjecture that it could explain the applicability of hydrodynamics outside its expected regime of validity since far from equilibrium and close to equilibrium may be related through transformations of this type.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Machine Learning-Driven Conservative-to-Primitive Conversion in Hybrid Piecewise Polytropic and Tabulated Equations of State

We present a novel machine learning (ML)-based method to accelerate conservative-to-primitive inversion, focusing on hybrid piecewise polytropic and tabulated equations of state. Traditional root-finding techniques are computationally expensive, particularly for large-scale relativistic hydrodynamics simulations. To address this, we employ feedforward neural networks (NNC2PS and NNC2PL), trained in PyTorch (2.0+) and optimized for GPU inference using NVIDIA TensorRT (8.4.1), achieving significant speedups with minimal accuracy loss. The NNC2PS model achieves 𝐿 1 and 𝐿 ∞ errors of 4.54 × 10 −7 and 3.44 × 10−6, respectively, while the NNC2PL model exhibits even lower error values. TensorRT optimization with mixed-precision deployment substantially accelerates performance compared to traditional root-finding methods. Specifically, the mixed-precision TensorRT engine for NNC2PS achieves inference speeds approximately 400 times faster than a traditional single-threaded CPU implementation for a dataset size of 1,000,000 points. Ideal parallelization across an entire compute node in the Delta supercomputer (dual AMD 64-core 2.45 GHz Milan processors and 8 NVIDIA A100 GPUs with 40 GB HBM2 RAM and NVLink) predicts a 25-fold speedup for TensorRT over an optimally parallelized numerical method when processing 8 million data points. Moreover, the ML method exhibits sub-linear scaling with increasing dataset sizes. We release the scientific software developed, enabling further validation and extension of our findings. By exploiting the underlying symmetries within the equation of state, these findings highlight the potential of ML, combined with GPU optimization and model quantization, to accelerate conservative-to-primitive inversion in relativistic hydrodynamics simulations.

conservative-to-primitive conversion↗
NASA NTRS1988

Cosmic-ray hydrodynamics at relativistic shocks

A test particle description of first-order Fermi acceleration of cosmic rays at plane-parallel relativistic shocks is presented in which the pitch angle distribution is obtained explicitly and used to determine the pressure anisotropy or variable Eddington factor as a function of position. The Eddington factor is in turn used to find the hydrodynamical cosmic-ray sound speed, taking the anisotropy into account. For upstream flows speeds above c/sq rt 3, it is shown that the cosmic-ray sound speed exceeds the shock speed in a region just upstream of the shock and decreases montonically with distance further upstream, reaching the value c/sq rt 3 far from the shock. Thus, the shock succeeds in maintaining an anisotropy in the cosmic-ray distribution which is sufficient to keep the flow subsonic with respect to the cosmic-ray sound speed both downstream and in a finite region upstream. The astrophysical implications of these findings are discussed.

Kirk, J. G.↗
NASA NTRS1993

Hydrodynamics of Relativistic Fireballs

Many models of gamma-ray bursts involve a fireball, which is an optically thick concentration of radiation energy with a high ratio of energy density to rest mass. We examine analytically and numerically the evolution of a relativistic fireball. We show that, after an early rearrangement phase, most of the matter and energy in the fireball is concentrated within a narrow shell. The shell propagates at nearly the speed of light, with a frozen radial profile, and according to a simple set of scaling laws. The spectrum of the escaping radiation is harder at early times and softer later on. Depending on the initial energy-to-mass ratio, the final outcome of a fireball is either photons with roughly the initial temperature or ultrarelativistic baryons. In the latter case, the energy could be converted back to gamma-rays via interaction with surrounding material.

Piran, Tsvi↗
DOE OSTI2020

Far From Equilibrium Hydrodynamics and the Beam Energy Scan

The existence of hydrodynamic attractors in rapidly expanding relativistic systems has shed light on the success of relativistic hydrodynamics in describing heavy-ion collisions at zero chemical potential. As the search for the QCD critical point continues, it is important to investigate how out of equilibrium effects influence the trajectories on the QCD phase diagram. In this proceedings, we study a Bjorken expanding hydrodynamic system based on DMNR equations of motion with initial out of equilibrium effects and finite chemical potential in a system with a critical point. We find that the initial conditions are not unique for a specific freeze-out point, but rather the system can evolve to the same final state freeze-out point with a wide range of initial baryon chemical potential, µ B . For the same initial energy density and baryon density, depending on how far out of equilibrium the system begins, the initial µ B can vary by Δµ B ~ 350 MeV. Furthermore, our results indicate that knowledge of the out-of-equilibrium effects in the initial state provide vital information that influences the search for the QCD critical point.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Stochastic fluctuations in relativistic fluids: Causality, stability, and the information current

We develop a general formalism for introducing stochastic fluctuations around thermodynamic equilibrium which takes into account, for the first time, recent developments in the causality and stability properties of relativistic hydrodynamic theories. The method is valid for any covariantly stable theory of relativistic viscous fluid dynamics derived from a covariant maximum entropy principle. We illustrate the formalism with some applications, showing how it could be used to consistently introduce fluctuations in a model of relativistic heat diffusion and in conformally invariant Israel-Stewart theory in a general hydrodynamic frame. Furthermore, the latter example is used to study the hydrodynamic frame dependence of the symmetric two-point function of fluctuations of the energy-momentum tensor.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1997

Gamma-Ray Bursts and Relativistic MHD around Coalescing Neutron Stars

In this proposal we performed the first relativistic hydrodynamic simulations of the binary neutron stars near coalescence. Calculations were performed in three spatial dimensions and were supplemented with magnetohydrodynamic simulations of the magnetic reconnection as the neutron star fluid responds to the relativistic forces of the binary. The ultimate goal of the proposed work was to develop a theoretical model for the temporal and spectral evolution of cosmological gamma-ray bursts produced during the magnetic reconnection. This is the first time such hydrodynamic calculations were performed and some surprising results were obtained. The neutrons stars were observed to compress and heat in the simulations. They could even collapse to black holes many orbits before coalescence. The question then arose as to whether this strange phenomenon could power a gamma-ray burst. The bulk of computing effort in this proposal was therefore devoted to an analysis of the implications of this possible new gamma-ray burst paradigm. The initial focus of the work performed under this proposal was therefore to utilize (3+1) dimensional and spherical numerical general relativistic hydrodynamics to study the origin, evolution, and parametric sensitivity of the collapse instability. We derived modified conditions of hydrostatic equilibrium for stars in the curved space of quasi-static orbits. We followed the magnetic field evolution for approximately one orbit (10 msec). The magnetic field energy was observed to exponentiate with an e-folding time of about 1 msec. Thus, the field should build up very quickly to a magnitude such that reconnection and back reaction of the fluid inhibits further growth.

Mathews, Grant J.↗
DOE OSTI2020

Effective viscosities in a hydrodynamically expanding boost-invariant QCD plasma

Background: The near-equilibrium properties of a QCD plasma can be encoded into transport coefficients such as bulk and shear viscosity. In QCD, the ratio of these transport coefficients to entropy density, ζ/$\textit{s}$ and η/$\textit{s}$, depends nontrivially on the plasma's temperature. This is unlike in conformal systems where they take constant values such as η/$\textit{s}$ = 1/(4π). Purpose: In this work, we show that in a (0 + 1)D boost-invariant fluid with no transverse expansion, a temperature-dependent ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ can be described by an equivalent effective viscosity $\langleζ/s\rangle_{\text{eff}}$ or $\langleη/s\rangle_{\text{eff}}$. This effective viscosity combines the actual temperature-dependent ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ with the temperature profile of the fluid. Further, we extend the concept of effective viscosity in systems with transverse expansion and discuss how effective viscosities can be used to identify families of ζ/$\textit{s(T)}$ and η /$\textit{s(T)}$ that lead to similar hydrodynamic evolution. Methods: The Navier-Stokes relativistic hydrodynamic equations are used to provide a first definition of effective viscosity, in (0 + 1)D and (1 + 1)D. In the (0 + 1)D case, the analysis is extended to Israel-Stewart-type second-order hydrodynamics to clarify the effect of higher-order hydrodynamics corrections on the effective viscosity. Results: In a boost-invariant fluid with no transverse expansion [(0 + 1)D], the effective viscosity is expressed as a simple integral of ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ over temperature, with a weight determined by the speed of sound of the fluid. The result is general for any equation of state with a moderate temperature dependence of the speed of sound, including the QCD equation of state. This definition of effective viscosity can be used to identify infinite families of ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ that produce essentially indistinguishable temperature profiles. In a boost-invariant cylindrical system [(1 + 1)D], a similar definition of effective viscosity is obtained in terms of characteristic trajectories in time and transverse direction. This leads to an infinite number of constraints on an infinite functional space for ζ/$\textit{s(T)}$ and η/$\textit{s(T)}$. Realistic examples are presented by using a finite number of constraints on a finite functional space. Conclusions: The definition of effective viscosity in a (0 + 1)D system clarifies how infinite families of ζ/$\textit{s(T)}$ and η/$\textit{s(T)}$ can result in nearly identical hydrodynamic temperature profiles. By extending the study to a boost-invariant cylindrical [(1 + 1)D] fluid, we identify an approximate but more general definition of effective viscosity that highlights the potential and limits of the concept of effective viscosity in fluids with limited symmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Anomaly-induced edge currents in hydrodynamics with parity anomaly

Here in this paper, we discuss relativistic hydrodynamics for a massless Dirac fermion in (2 + 1) dimensions, which has the parity anomaly—a global ’t Hooft anomaly between U(1) and parity symmetries. We investigate how hydrodynamics implements the party anomaly, particularly focusing on the transport phenomena at the boundary. Based on the parity anomaly matching and the second law of local thermodynamics, we find U(1) and entropy currents localized at the boundary, as well as the bulk anomalous current with vanishing divergence. These edge currents are similar to the ( 1 + 1)-dimensional chiral transports, but the coefficients are given by half of theirs. We also generalize our discussion to more general anomalies among multiple U(1) symmetries and single $\mathbb{Z}_2$ symmetry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Causal hydrodynamic fluctuations in non-static and inhomogeneous backgrounds

Highlights: •The fluctuation–dissipation relation is modified by inhomogeneous backgrounds. •The constitutive equations are restricted by causality and fluctuation properties. •The structure of constitutive equation gives a natural cutoff scale of the theory. •The memory function with tensor structures in inhomogeneous backgrounds is obtained. -- Abstract: To integrate hydrodynamic fluctuations, namely thermal fluctuations of hydrodynamics, into dynamical models of high-energy nuclear collisions based on relativistic hydrodynamics, the property of the hydrodynamic fluctuations given by the fluctuation–dissipation relation should be carefully investigated. The fluctuation–dissipation relation for causal dissipative hydrodynamics with the finite relaxation time is naturally given in the integral form of the constitutive equation by the linear-response theory. While, the differential form of the constitutive equation is commonly used in analytic investigations and dynamical calculations for practical reasons. We give the fluctuation–dissipation relation for the general linear-response differential form and discuss the restrictions to the structure of the differential form, which comes from the causality and the positive semi-definiteness of the noise autocorrelation, and also the relation of those restrictions to the cutoff scale of the hydrodynamic fluctuations. We also give the fluctuation–dissipation relation for the integral form in non-static and inhomogeneous background by introducing new tensors, the pathline projectors. We find new modification terms to the fluctuation–dissipation relation for the differential form in non-static and inhomogeneous background which are particularly important in dynamical models to describe rapidly expanding systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

thornado-hydro: A Discontinuous Galerkin Method for Supernova Hydrodynamics with Nuclear Equations of State

This paper describes algorithms for non-relativistic hydrodynamics in the toolkit for high-order neutrino radiation hydrodynamics (thornado), which is being developed for multiphysics simulations of core-collapse supernovae (CCSNe) and related problems with Runge–Kutta discontinuous Galerkin (RKDG) methods. More specifically, thornado employs a spectral type nodal collocation approximation, and we have extended limiters — a slope limiter to prevent non-physical oscillations and a bound-enforcing limiter to prevent non-physical states — from the standard RKDG framework to be able to accommodate a tabulated nuclear equation of state (EoS). To demonstrate the efficacy of the algorithms with a nuclear EoS, we first present numerical results from basic test problems in idealized settings in one and two spatial dimensions, employing Cartesian, spherical-polar, and cylindrical coordinates. Then, we apply the RKDG method to the problem of adiabatic collapse, shock formation, and shock propagation in spherical symmetry, initiated with a 15 M ⊙ progenitor. Herein, we find that the extended limiters improve the fidelity and robustness of the RKDG method in idealized settings. The bound-enforcing limiter improves robustness of the RKDG method in the adiabatic collapse application, while we find that slope limiting in characteristic fields is vulnerable to structures in the EoS — more specifically, in the phase transition from nuclei and nucleons to bulk nuclear matter. The success of these applications marks an important step toward applying RKDG methods to more realistic CCSN simulations with thornado in the future.

Nuclear astrophysics (1129)↗
DOE OSTI2024

Unveiling Baryon Charge Carriers through Charge Stopping in Isobar Collisions

Utilizing a comprehensive (3 + 1)⁢D relativistic hydrodynamic framework with multiple conserved charge currents and charge-dependent lattice-QCD-based equation of state, we study the baryon and electric charge number deposition at midrapidity in isobar Ru + Ru and Zr + Zr collisions at the center of mass energy $\sqrt{s}_{NN}$ = 200 GeV. Comparing our predictions with upcoming experimental data from the Relativistic Heavy Ion Collider will shed light on the existence of baryon junctions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗