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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2024

Dynamical Aspects of the Quark Gluon Plasma

Shortly after the Big Bang the entire universe was filled with a nearly perfect fluid known as the Quark Gluon Plasma. Relativistic heavy ion collisions can now reproduce this fluid in the laboratory where the Quark Gluon Plasma exhibits a rapid but smooth cross-over phase transition into hadrons at vanishing net-baryon densities. Recent experiments plan to explore finite baryon densities where a critical point is expected. If found, this would mark the first discovery of a critical point in a relativistic system described by a fundamental theory of nature, which would have far-reaching consequences for high-energy nuclear physics and nuclear astrophysics (such as in neutron star mergers). Characteristic temperatures of equilibrium (e.g. the inflection point of the entropy density) and transport coefficients (e.g. minimum of the shear viscosity over entropy density) vary widely at a cross-over phase transition but converge at a critical point, and extracting the behavior of these characteristic temperatures is a major focal point of this research. Specifically, the interplay between strange and light hadrons is exploited to study the flavor hierarchy in the cross-over region. To investigate this, a viscous relativistic hydrodynamics framework with two conserved charges is being developed into a new open-source code along with initial conditions that contain baryon number and strangeness. Flow observables sensitive to the equation of state and transport coefficients are calculated across beam energies. New techniques are being developed to study this flavor hierarchy from first principles and to extract the characteristic temperatures from experimental data. Through this new dynamical framework, this project provides essential guidance to the Beam Energy Scan II runs at Relativistic Heavy-Ion Collider and the future Facility for Antiproton and Ion Research facility in the search for the Quantum Chromodynamic critical point and subsequent investigation of the baryon-rich Quark Gluon Plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2022

Λ spin polarization in event-by-event relativistic heavy-ion collisions

We present a systematic study of $\Lambda$ hyperon's polarization observables using event-by-event (3+1)D relativistic hydrodynamics. The effects of initial hot spot size and QGP's specific shear viscosity on the polarization observables are quantified. We examine the effects of the two formulations of the thermal shear tensor on the polarization observables using the same hydrodynamic background. With event-by-event simulations, we make predictions for the Fourier coefficients of $\Lambda$'s longitudinal polarization $P^z$ with respect to the event planes of different orders of anisotropic flow. We propose new correlations among the Fourier coefficients of $P^z$ and charged hadron anisotropic flow coefficients to further test the mapping from fluid velocity gradients to hyperon's polarization. Finally, we present a system size scan with Au+Au, Ru+Ru, and O+O collisions at $\sqrt{s_\mathrm{NN}} = 200$ GeV to study the system size dependence of polarization observables at the Relativistic Heavy-ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

3D structure of anisotropic flow in small collision systems at energies available at the BNL Relativistic Heavy Ion Collider

Here, we present (3 + 1)-dimensional [(3 + 1)D] dynamical simulations of asymmetric nuclear collisions at the BNL Relativistic Heavy Ion Collider (RHIC). Employing a dynamical initial state model coupled to (3 + 1)D viscous relativistic hydrodynamics, we explore the rapidity dependence of anisotropic flow in the RHIC small system scan at 200 GeV center-of-mass energy. We calibrate parameters to describe central 3 He + Au collisions and make extrapolations to d + Au and p + Au collisions. Our calculations demonstrate that approximately 50% of the v 3 (p T ) difference between the measurements by the STAR and PHENIX Collaborations can be explained by the use of reference flow vectors from different rapidity regions. This emphasizes the importance of longitudinal flow decorrelation for anisotropic flow measurements in asymmetric nuclear collisions, and the need for (3 + 1)D simulations. We also present results for the beam energy dependence of particle spectra and anisotropic flow in d + Au collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 OSTI2022

Relativistic wind farm effect: Possibly turbulent flow of a charged, massless relativistic fluid in graphene

At low Reynolds numbers, the wind flow in the wake of a single wind turbine is generally not turbulent. However, turbines in wind farms affect each other's wakes so that a turbulent flow can arise. In the present work, an analog of this effect for the massless charge carrier flow around obstacles in graphene is outlined. We use a relativistic hydrodynamic simulation to analyze the flow in a sample containing impurities. Depending on the density of impurities in the sample, we indeed find evidence for potentially turbulent flow and discuss experimental consequences.

Watson, Mark (ORCID:0000000258959949)↗
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 OSTI2021

Tidal disruption discs formed and fed by stream–stream and stream–disc interactions in global GRHD simulations

When a star passes close to a supermassive black hole (BH), the BH’s tidal forces rip it apart into a thin stream, leading to a tidal disruption event (TDE). In this work, we study the post-disruption phase of TDEs in general relativistic hydrodynamics (GRHD) using our GPU-accelerated code h-amr. We carry out the first grid-based simulation of a deep-penetration TDE (β = 7) with realistic system parameters: a black hole-to-star mass ratio of 10 6 , a parabolic stellar trajectory, and a non-zero BH spin. We also carry out a simulation of a tilted TDE whose stellar orbit is inclined relative to the BH midplane. We show that for our aligned TDE, an accretion disc forms due to the dissipation of orbital energy with ~20 percent of the infalling material reaching the BH. The dissipation is initially dominated by violent self-intersections and later by stream–disc interactions near the pericentre. The self-intersections completely disrupt the incoming stream, resulting in five distinct self-intersection events separated by approximately 12 h and a flaring in the accretion rate. We also find that the disc is eccentric with mean eccentricity e ≈ 0.88. For our tilted TDE, we find only partial self-intersections due to nodal precession near pericentre. Although these partial intersections eject gas out of the orbital plane, an accretion disc still forms with a similar accreted fraction of the material to the aligned case. These results have important implications for disc formation in realistic tidal disruptions. For instance, the periodicity in accretion rate induced by the complete stream disruption may explain the flaring events from Swift J1644+57.

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

Minimal one-dimensional model of bad metal behavior from fast particle-hole scattering

A strongly interacting plasma of linearly dispersing electron and hole excitations in two spatial dimensions (2D), also known as a Dirac fluid, can be captured by relativistic hydrodynamics and shares many universal features with other quantum critical systems. Here, we propose a one-dimensional (1D) model to capture key aspects of the 2D Dirac fluid while including lattice effects and being amenable to nonperturbative computation. When interactions are added to the Dirac-like 1D dispersion without opening a gap, we show that this kind of irrelevant interaction is able to preserve Fermi-liquid-like quasiparticle features while relaxing a zero-momentum charge current via collisions between particle-hole excitations, leading to resistivity that is linear in temperature via a mechanism previously discussed for large-diameter metallic carbon nanotubes. We further provide a microscopic lattice model and obtain numerical results via density-matrix renormalization group simulations, which support the above physical picture. The limits on such fast relaxation at strong coupling are of considerable interest because of the ubiquity of bad metals in experiments.

1-dimensional systems↗
DOE OSTI2023

Utilizing the maximum likelihood estimator for flow analysis

We explore the possibility of evaluating flow harmonics by employing the maximum likelihood estimator (MLE). For a given finite multiplicity, the MLE simultaneously furnishes estimations for all the parameters of the underlying distribution function while efficiently suppressing the variance of measures. Also, the method provides a means to assess a specific class of mixed harmonics, which is not straightforwardly feasible by the approaches primarily based on particle correlations. The results are analyzed using the Wald, likelihood ratio, and score tests of hypotheses. Besides, the resultant flow harmonics obtained using MLE are compared with those derived using particle correlations and event plane methods. Here, the dependencies of extracted flow harmonics on the multiplicity of individual events and the total number of events are analyzed. It is shown that the proposed approach works efficiently to deal with the deficiency in detector acceptability. Moreover, we elaborate on a fictitious scenario where the event plane is not a well-defined quantity in the distribution function. For the latter case, the MLE is shown to largely perform better than the two-particle correlation estimator. In this regard, one concludes that the MLE furnishes a meaningful alternative to the existing approaches for flow analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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↗