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

Time-dependent, optically thick accretion onto a black hole

A fully relativistic hydrodynamics code which incorporates diffusive radiation transport is used to study time-dependent, spherically symmetric, optically thick accretion onto a black hole. It is found that matter free-falls into the hole regardless of whether the diffusion time scale is longer or shorter than the dynamical time. Nonadiabatic heating due to magnetic field reconnection is included. The internal energy thus generated affects the flow in a purely relativistic way, again ensuring free-fall collapse of the inflowing matter. Any matter enveloping a black hole will thus be swallowed on a dynamical time scale with relatively small net release of energy. The inclusion of angular momentum will not necessarily affect this conclusion.

Gilden, D. L.↗
NASA NTRS2020

Tidal Disruptions of Main Sequence Stars - I. Observable Quantities and their Dependence on Stellar and Black Hole Mass

This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main-sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses(MM0.1510)and six black hole masses( MM105 105BH7). We present here the results most relevant to observations. The effects of internal stellar structure and relativity decouple for both the disruption cross section and the characteristic energy width of the debris. Moreover, the full disruption cross section is almost independent of Må for Må/Me3. Independent of Må, relativistic effects increase the critical pericenter distance for full disruption events by up to a factor of∼3 relative to the Newtonian prediction. The probability of a direct capture is also independent of Må; at MBH/Me;5×106 this probability is equal to the probability of a complete disruption. The breadth of the debris energy distribution ΔE can differ from the standard estimate by factors of0.35−2, depending on Må and MBH, implying a corresponding change(∝(ΔE)−3/2)in the characteristic mass-return timescale. We provide analytic forms, suitable for use in both event rate estimates and parameter inference, to describe all these trends. For partial disruptions, we find a nearly universal relation between the star’s angular momentum and the fraction of Må remaining. Within the “empty loss-cone” regime, partial disruptions must precede full disruptions. These partial disruptions can drastically affect the rate and appearance of subsequent total disruptions.

Taeho Ryu↗
NASA NTRS2020

Tidal Disruptions of Main Sequence Stars - II. Simulation Methodology and Stellar Mass Dependence of the Character of Full Tidal Disruptions

This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main-sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses (0.15 < M_/M<10)and six black hole masses( ́10<MM105 105BH7). We present here the results most relevant to observations. The effects of internal stellar structure and relativity decouple for both the disruption cross section and the characteristic energy width of the debris. Moreover, the full disruption cross section is almost independent of Må for Må/Me3. Independent of Må, relativistic effects increase the critical pericenter distance for full disruption events by up to a factor of ∼3 relative to the Newtonian prediction. The probability of a direct capture is also independent of Må; at MBH/Me5×106this probability is equal to the probability of a complete disruption. The breadth of the debris energy distribution ΔE can differ from the standard estimate by factors of 0.35−2, depending on Må and MBH, implying a corresponding change(∝(ΔE)−3/2)in the characteristic mass-return timescale. We provide analytic forms, suitable for use in both event rate estimates and parameter inference, to describe all these trends. For partial disruptions, we find a nearly universal relation between the star’s angular momentum and the fraction of Må remaining. Within the “empty loss-cone” regime, partial disruptions must precede full disruptions. These partial disruptions can drastically affect the rate and appearance of subsequent total disruptions

Taeho Ryu↗
DOE OSTI2026

Nonlinear causality and strong hyperbolicity of baryon-rich Israel-Stewart hydrodynamics

We present the first set of fully-nonlinear, necessary and sufficient conditions guaranteeing causal evolution of the initial data for the Israel-Stewart hydrodynamic equations with shear and bulk viscosity coupled to a nonzero baryon current. These constraints not only provide nonlinear causality: they also (a) guarantee the existence of a locally well-posed evolution of the initial data (they enforce strong hyperbolicity) when excluding the endpoints of the bounds, (b) arise from purely algebraic constraints that make no underlying symmetry assumptions on the degrees of freedom and (c) propagate the relevant symmetries of the degrees of freedom over the entire evolution of the problem. Our work enforces a mathematically rigorous foundation for future studies of viscous relativistic hydrodynamics with baryon-rich matter including neutron star mergers and heavy-ion collisions.

Cordeiro, Ian [University of Illinois Urbana-Champ↗
DOE OSTI2025

Stochastic fluctuations and the relaxation time in transient relativistic fluids

We argue that the ratio between the shear viscosity and the shear relaxation time, η/τ π , should be defined as a thermodynamic quantity obtained from the equal-time symmetric correlator of the shear-stress tensor. In kinetic theory, we show that this ratio does not depend on the type of interaction. Similarly, an exact expression for this ratio is obtained for holographic gauge theories. We also determine how stochastic fluctuations change η/τ π in transient relativistic hydrodynamics and show that thermal fluctuations do not spoil causality and stability.

Denicol, Gabriel S. [Universidade Federal Fluminen↗
DOE OSTI2024

Initial data for first-order causal viscous conformal fluids in general relativity

We solve the Einstein constraint equations for a first-order causal viscous relativistic hydrodynamic theory in the case of a conformal fluid. For such a theory, a direct application of the conformal method does not lead to a decoupling of the equations, even for constant-mean curvature initial data. We combine the conformal method applied to a background perfect fluid theory with a perturbative argument in order to obtain the result.

Disconzi, Marcelo (ORCID:0000000234497778)↗
DOE OSTI2025

Collision-Energy Dependence in Heavy-Ion Collisions from Nonlinear QCD Evolution

We explore the effects of including the energy dependence determined from evolution equations within the color glass condensate framework on observables in ultrarelativistic heavy-ion collisions. This amounts to integrating the JIMWLK evolution equations into the impact parameter-dependent glasma model, which is then coupled to viscous relativistic hydrodynamics. This methodology allows for a systematic representation of nuclei at specific Bjorken-x values, which are probed at different center-of-mass energies of the collision and rapidities of final state particles. Comparing the methodology to the conventional impact parameter-dependent glasma model, we find significant effects on multiplicity distributions and particle spectra, especially in smaller collision systems at the highest center-of-mass energies. Our results highlight the importance of incorporating nonlinear QCD evolution in the description of heavy-ion collisions at varying center-of-mass energies, as the precise extraction of transport coefficients will be affected. This Letter establishes a robust framework for understanding the quark gluon plasma and nuclear structure at high energy, integrating small-x physics into the initial conditions of heavy-ion collisions.

Mäntysaari, Heikki↗
DOE OSTI2024

Covariant formulation of spinodal decomposition in rapidly expanding quark gluon plasma

Quantum chromodynamics (QCD) is expected to have a first order phase transition between the confined hadron gas and the deconfined quark gluon plasma at high baryon densities. This will result in phase boundary effects in the metastable and unstable regions. It is important to include these effects in phenomenological models of heavy ion collisions to identify experimental signatures of a phase transition. This requires building intuition on phase separation in rapidly expanding fluids. In this work we present the covariant equations of relativistic hydrodynamics with a phase boundary, provide prescriptions to extend the equation of state to metastable and unstable regions, and show the effects of spinodal separation in a Bjorken flow. Published by the American Physical Society 2024

Kapusta, Joseph I. (ORCID:0000000259429835)↗
DOE OSTI2026

Impact of QCD energy evolution on observables in heavy-ion collisions

Here, we study how the inclusion of energy dependence as dictated by quantum chromodynamic (QCD) small-x evolution equations affects key observables in ultra-relativistic heavy-ion collisions. Specifically, we incorporate JIMWLK evolution into the IP-Glasma framework, which serves as the initial condition for a simulation pipeline that includes viscous relativistic hydrodynamics and a hadronic afterburner. This approach enables a consistent modeling of highly energetic nuclei across varying Bjorken-x values, which are relevant for different collision energies and rapidity regions. In comparison to the standard IP-Glasma setup without small-x evolution, we observe pronounced changes in particle multiplicities and spectral distributions, especially in smaller systems and at the highest available energies. We further explore effects on anisotropic flow observables and correlations between mean transverse momentum and elliptic flow. Our findings underscore the critical role of nonlinear QCD evolution in accurately modeling the early stages of heavy-ion collisions, as well as its implications for extracting transport properties of the quark-gluon plasma.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Thermalization from quantum entanglement: Jet simulations in the massive Schwinger model

We investigate the emergence of thermalization in a quantum-field-theoretic model mimicking the production of jets in QCD: the massive lattice Schwinger model coupled to external sources. Specifically, we compute the expectation values of local operators as functions of time and compare them to their thermal counterparts, quantify the overlap between the evolving density matrix and the thermal one, and compare the dynamics of the energy-momentum tensor to predictions from relativistic hydrodynamics. Through these studies, we find that the system approaches thermalization at late times and elucidate the mechanisms by which quantum entanglement drives thermalization in closed field-theoretic systems. Our results show how thermodynamic behavior emerges in real time from unitary quantum dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

r -process Nucleosynthesis from Hyperaccreting Neutron Stars in Common Envelopes

We investigate nuclear reactions and feedback in hyperaccreting neutron star environments, considering accretion rates in the range 0.3–3 × 10 4 M ⊙ yr −1 , typical of short-period compact-object binaries in common envelopes. Our models account for weak reactions, neutrino energy loss, nuclear energy release, pair production, degenerate equations of state, and general relativistic hydrodynamics. Depending on the accretion rates, these systems can develop both proton- and neutron-rich atmospheres with strong convective instabilities linking the neutrino sphere to the outgoing accretion shock inside the radiation trapping zone. Convection drives nucleons through multiple heating and cooling cycles, with photodisintegration dominating during the heating phase and heavy element synthesis during the cooling phase, ejecting material with abundances that depend on the accretion rate and depth of the final decompression trajectory. The turbulent nature of convective currents is conducive to creating a wide range of nuclear products through a variety of effects, including nuclear statistical equilibrium freeze-out and the r-, p-, and γ-processes. We also observe a novel multistep process in reheated trajectories, consisting of proton-capture and photodissociation reactions operating on r-process seeds, producing overall neutron-deficient isotopes. A significant amount of infalling gas experiences high entropy and short (millisecond) freeze-out timescales capable of making r-process elements with high overabundances through a disequilibrium effect between neutrons and α-particles that does not require an excess of neutrons.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1985

The cosmic-ray shock structure problem for relativistic shocks

The time asymptotic behaviour of a relativistic (parallel) shock wave significantly modified by the diffusive acceleration of cosmic-rays is investigated by means of relativistic hydrodynamical equations for both the cosmic-rays and thermal gas. The form of the shock structure equation and the dispersion relation for both long and short wavelength waves in the system are obtained. The dependence of the shock acceleration efficiency on the upstream fluid spped, long wavelength Mach number and the ratio N = P sub co/cP sub co+P sub go)(Psub co and P sub go are the upstream cosmic-ray and thermal gas pressures respectively) are studied.

Webb, G. M.↗