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 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 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 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 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↗
NASA NTRS2001

Computational Relativistic Astrophysics Using the Flowfield-Dependent Variation Theory

Theoretical models, observations and measurements have preoccupied astrophysicists for many centuries. Only in recent years, has the theory of relativity as applied to astrophysical flows met the challenges of how the governing equations can be solved numerically with accuracy and efficiency. Even without the effects of relativity, the physics of magnetohydrodynamic flow instability, turbulence, radiation, and enhanced transport in accretion disks has not been completely resolved. Relativistic effects become pronounced in such cases as jet formation from black hole magnetized accretion disks and also in the study of Gamma-Ray bursts (GRB). Thus, 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 flowfield-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 OSTI2024

Causal and stable first-order chiral hydrodynamics

We derive the set of inequalities that is necessary and sufficient for nonlinear causality and linear stability of first-order relativistic hydrodynamics with either a U ( 1 ) V conserved current or a U ( 1 ) A current with a chiral anomaly or both. Our results apply to generic hydrodynamic frames in which no relations among the transport parameters are imposed. Furthermore, our analysis yields, to the best of our knowledge, the first theory of viscous chiral hydrodynamics proven to be causal and stable. We find that causality demands the absence of vorticity-induced heat flux, forcing a departure from the thermodynamic frame in the chiral case. The inequalities for causality and stability define a hypervolume in the space of transport parameters, wherein each point corresponds to a consistent formulation. Notably, causality is determined by just three combinations of transport parameters. We present our results in a form amenable to numerical hydrodynamic simulations. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2025

Energy Spectrum and Mass Composition of Ultra-high-energy Cosmic Rays Originating from Relativistic Jets of Nearby Radio Galaxies

Relativistic jets of radio galaxies (RGs) are possible sources of ultra-high-energy cosmic rays (UHECRs). Recent studies combining relativistic hydrodynamic simulations with Monte Carlo particle transport have demonstrated that UHECRs can be accelerated to energies beyond 10 20 eV through shocks, turbulence, and relativistic shear in jet-induced flows of Fanaroff–Riley type RGs. The resulting time-asymptotic UHECR spectrum is well modeled by a double power law with an “extended” exponential cutoff, primarily shaped by relativistic shear acceleration. In this study, we adopt this novel source spectrum and simulate the propagation of UHECRs from nearby RGs using the CRPropa code. We focus on Virgo A (Vir A), Centaurus A (Cen A), Fornax A (For A), and Cygnus A (Cyg A), expected to be the most prominent UHECR sources among RGs. We then analyze the energy spectrum and mass composition of UHECRs arriving at Earth. We find that, due to the extended high-energy tail in the source spectrum, UHECRs from Vir A, which has a higher Lorentz factor, exhibit a higher flux at the highest energies and a lighter mass composition at Earth compared to those from Cen A and For A with lower Lorentz factors. Despite Cyg A having an even higher Lorentz factor, the large distance limits its contribution. With a small number of nearby prominent RGs, our findings suggest that if RGs are the major sources of UHECRs, the energy spectrum and mass composition of observed UHECRs would exhibit hemispheric differences between the Northern and Southern skies at the highest energies.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2020

Tidal Disruptions of Main-sequence Stars. IV. Relativistic Effects and Dependence on Black Hole Mass

Using a suite of fully relativistic hydrodynamic simulations applied to main-sequence stars with realistic internal density profiles, we examine full and partial tidal disruptions across a wide range of black hole mass( MM105 105BH7)and stellar mass (MM0.33) as larger MBH leads to stronger relativistic effects. For fixed Må, as MBH increases, the ratio of the maximum pericenter distance yielding full disruptions (t) to its Newtonian prediction rises rapidly, becoming triple the Newtonian value for= ́MM510BH7, while the ratio of the energy width of the stellar debris for full disruptions to the Newtonian prediction decreases steeply, resulting in a factor of 2 correction at= ́MM510BH7. We provide approximate formulae that express the relativistic corrections of both t and the energy width relative to their Newtonian approximate estimates. For partial disruptions, we find that the fractional remnant mass for a given ratio of the pericenter to t is higher for larger MBH. These results have several implications. As MBH increases above~M107, the cross section for complete disruptions is suppressed by competition with direct capture. However, the cross-section ratio for partial to complete disruptions depends only weakly on MBH. The relativistic correction to the debris energy width delays the time of peak mass-return rate and diminishes the magnitude of the peak return rate. ForMM10BH7, the MBH-dependence of the full disruption cross section and the peak mass-return rate and time is influenced more by relativistic effects than by Newtonian dynamics.

Taeho Ryu↗
DOE OSTI2024

Four-dimensional QCD equation of state with multiple chemical potentials

Here, we construct a four-dimensional version of the equation of state (EoS) model neos, neos-4d, as a function of the temperature and chemical potentials of baryon, electric charge, and strangeness for the hot and dense quantum chromodynamics (QCD) matter created in relativistic nuclear collisions. This EoS enables multiple conserved charge current evolution in a relativistic fluid. Input from lattice QCD simulations and a hadron resonance gas model is considered for constructing the equation of state. We investigate its applicability to the relativistic hydrodynamic description of nuclear collisions and present a method for efficient numerical implementation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2020

Tidal Disruptions of Main-sequence Stars. IV. Relativistic Effects and Dependence on Black Hole Mass

Using a suite of fully relativistic hydrodynamic simulations applied to main-sequence stars with realistic internal density profiles, we examine full and partial tidal disruptions across a wide range of black hole mass (10(exp 5)≤M(BH)/M(ʘ)≤5 x 10(exp 7) and stellar mass (0.3≤M(*)/M(ʘ)≤3) as larger M(BH) leads to stronger relativistic effects. For fixed M( ), as M(BH) increases, the ratio of the maximum pericenter distance yielding full disruptions ( ) to its Newtonian prediction rises rapidly, becoming triple the Newtonian value for M(BH)=5 x 10(exp 7)M(ʘ), while the ratio of the energy width of the stellar debris for full disruptions to the Newtonian prediction decreases steeply, resulting in a factor of 2 correction at M(BH)=5 x 10(exp 7)M(ʘ). We provide approximate formulae that express the relativistic corrections of both R(t) and the energy width relative to their Newtonian approximate estimates. For partial disruptions, we find that the fractional remnant mass for a given ratio of the pericenter to R(t) is higher for larger M(BH). These results have several implications. As M(BH) increases above ~10(exp 7)M(ʘ), the cross section for complete disruptions is suppressed by competition with direct capture. However, the cross-section ratio for partial to complete disruptions depends only weakly on M(BH). The relativistic correction to the debris energy width delays the time of peak mass-return rate and diminishes the magnitude of the peak return rate. For M(BH)≳10(exp 7)M(ʘ), the M(BH)-dependence of the full disruption cross section and the peak mass-return rate and time is influenced more by relativistic effects than by Newtonian dynamics.

Taeho Ryu↗