Search NASA⌕ Search

SEARCH · Search NASA

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Effective Action for Relativistic Hydrodynamics from the Crooks Fluctuation Theorem

A new effective theory framework for fluctuating hydrodynamics in the relativistic regime is derived using standard thermodynamical principles and general properties of nonequilibrium stochastic dynamics. For the first time, we establish clear and concise conditions for ensuring that the resulting effective theories are causal, stable, and well-posed within general relativity. Here, these properties are independent of spacetime foliation and are valid in the full nonlinear regime. Out-of-equilibrium fluctuations are constrained by a relativistically covariant version of the Crooks fluctuation theorem, which determines how the entropy production is distributed even when the system is driven by an external force. This leads to an emerging Z 2 symmetry responsible for imposing fluctuation-dissipation relations for n-point correlation functions, which matches the standard constraints for the Schwinger-Keldysh effective action.

Fluctuation theorems↗
NASA NTRS2003

A Numerical Gamma-Ray Burst Simulation Using Three-Dimensional Relativistic Hydrodynamics: The Transition from Spherical to Jet-like Expansion

Utilizing 3D relativistic hydrodynamical calculations, we have examined the evolution of an expanding relativistic blob of gas intended to be representative of a jet associated with ejecta from an extremely energetic event such as a hypernova, that produces a gamma-ray burst (Aloy et al. 2000; Tan, Matzner, & McKee 2001; MacFadyen, Woosley, & Heger 2001, Zhang, Woosley, & Heger 2003, Zhang, Woosley, & MacFadyen 2003). Since these are the first such calculations applied to the blob during the time in which the afterglow radiation is produced, we have purposely kept them simple in an effort to concentrate on the most fundamental aspects of the physics. We restrict our attention to the transition from spherical to jetlike expansion that occurs during the time that the Lorentz factor becomes less than the reciprocal of the jet spreading angle. We have not yet attached specific numbers to our results. From the SRHD equations, one sees that the relevant quantities are the ratios of pressure to density, and of distance to time. If we specify either one of these two sets of numbers, the other one is also determined.

Cannizzo, John K.↗
DOE OSTI2026

A new effective theory for stochastic relativistic hydrodynamics

Thermal fluctuations are a fundamental feature of dissipative systems that are essential for understanding physics near the expected critical point of QCD and in small systems. When such fluctuations are modeled naively in relativistic systems, strange features can appear such as negative self-correlation functions. We construct an effective theory for nonlinear stochastic relativistic hydrodynamics that ensure a well-posed mathematical formulation. Using Crooks fluctuation theorem, we derive a symmetry of the effective action that incorporates fluctuations through a suitable free energy functional. For divergence type theories, the action can then be fully specified using a single vector generating current. The equations of motion obtained using this procedure are guaranteed to be flux conservative and symmetric hyperbolic when the dynamics is causal. This ensures that these equations are well-posed (for suitable initial data) and are in a form that can easily be simulated, including with Metropolis techniques.

Mullins, Nicki [University of Illinois at Urbana-C↗
DOE OSTI2025

General relativistic hydrodynamic simulations of binary strange star mergers

We perform fully general-relativistic simulations of binary strange star mergers considering two different approaches for thermal effects. The first uses a cold equation of state (EOS) derived from a modified version of the MIT bag model which is then supplemented by a Γ-law correction. The second approach employs a microphysical description of the finite-temperature effects. We describe results obtained with the two treatments, highlighting the influence of thermal effects. We find that the postmerger dynamics differs significantly in the two cases, leading to quantitative differences in the postmerger gravitational-wave spectrum and ejecta mass. The peak frequency of the postmerger gravitational-wave emission is consistent with the established quasi-universal relations for binary neutron star mergers and as a result, our simulations cannot distinguish between mergers of neutron stars and those of strange stars. Our models with realistic treatment of finite-temperature effects produce a significant amount of ejecta ≳0.02 M ⊙ ​. Here, the resulting flux of strangelets near the Earth, computed assuming that all neutron star mergers are in fact strange-stars mergers and that the binary considered here is representative, is in tension with experimental upper limits. As such, our results tentatively disfavor a scenario in which strange-quark matter is the lowest energy state of matter.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2002

Computational Relativistic Astrophysics Using the Flow Field-Dependent Variation Theory

We present our method for solving general relativistic nonideal hydrodynamics. Relativistic effects become pronounced in such cases as jet formation from black hole magnetized accretion disks which may lead to the study of gamma-ray bursts. Nonideal flows are present where radiation, magnetic forces, viscosities, and turbulence play an important role. Our concern in this paper is to reexamine existing numerical simulation tools as to the accuracy and efficiency of computations and introduce a new approach known as the flow field-dependent variation (FDV) method. The main feature of the FDV method consists of accommodating discontinuities of shock waves and high gradients of flow variables such as occur in turbulence and unstable motions. In this paper, the physics involved in the solution of relativistic hydrodynamics and solution strategies of the FDV theory are elaborated. The general relativistic astrophysical flow and shock solver (GRAFSS) is introduced, and some simple example problems for computational relativistic astrophysics (CRA) are demonstrated.

Richardson, G. A.↗
DOE OSTI2026

Radial oscillations of viscous neutron stars: Zero diffusion case

In this article, the spectrum of radial oscillations of neutron stars is systematically studied within two frameworks of viscous relativistic hydrodynamics: the relativistic Navier-Stokes and Israel-Stewart theories. A correspondence is established between the discrete stellar eigenmodes and the continuous dispersion relation of perturbations around a homogeneous fluid, providing a basis for interpreting our numerical results. We analyze the Newtonian limit and assess the impact of relativistic corrections, such as the gravitational redshifting of microscopic relaxation timescales. We show that bulk viscosity can significantly affect the behavior of both hydrodynamic and nonhydrodynamic fundamental modes, and that, depending on the magnitude of the viscous effects, it is the nonhydrodynamic mode that becomes unstable beyond the turning point in a sequence of equilibrium configurations. These results provide a useful step toward systematic studies of neutron star quasinormal modes in the presence of viscosity.

79 ASTRONOMY AND ASTROPHYSICS↗
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 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 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↗
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 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↗