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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Ultra-short pulse laser acceleration of protons to 80 MeV from cryogenic hydrogen jets tailored to near-critical density

Laser plasma-based particle accelerators attract great interest in fields where conventional accelerators reach limits based on size, cost or beam parameters. Despite the fact that particle in cell simulations have predicted several advantageous ion acceleration schemes, laser accelerators have not yet reached their full potential in producing simultaneous high-radiation doses at high particle energies. The most stringent limitation is the lack of a suitable high-repetition rate target that also provides a high degree of control of the plasma conditions required to access these advanced regimes. Here, we demonstrate that the interaction of petawatt-class laser pulses with a pre-formed micrometer-sized cryogenic hydrogen jet plasma overcomes these limitations enabling tailored density scans from the solid to the underdense regime. Our proof-of-concept experiment demonstrates that the near-critical plasma density profile produces proton energies of up to 80 MeV. Based on hydrodynamic and three-dimensional particle in cell simulations, transition between different acceleration schemes are shown, suggesting enhanced proton acceleration at the relativistic transparency front for the optimal case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The impact of r -process heating on the dynamics of neutron star merger accretion disc winds and their electromagnetic radiation

Neutron star merger accretion discs can launch neutron-rich winds of >10 -2 M⊙. This ejecta is a prime site for r-process nucleosynthesis, which will produce a range of radioactive heavy nuclei. The decay of these nuclei releases enough energy to accelerate portions of the wind by ~0.1c. Here, we investigate the effect of r-process heating on the dynamical evolution of disc winds. We extract the wind from a 3D general relativistic magnetohydrodynamic simulation of a disc from a post-merger system. This is used to create inner boundary conditions for 2D hydrodynamic simulations that continue the original 3D simulation. We perform two such simulations: one that includes the r-process heating, and another one that does not. We follow the hydrodynamic simulations until the winds reach homology (60 s). Using time-dependent multifrequency multidimensional Monte Carlo radiation transport simulations, we then calculate the kilonova light curves from the winds with and without dynamical r-process heating. We find that the r-process heating can substantially alter the velocity distribution of the wind, shifting the mass-weighted median velocity from 0.06c to 0.12c. The inclusion of the dynamical r-process heating makes the light curve brighter and bluer at $\sim 1\, \mathrm{d}$ post-merger. However, the high-velocity tail of the ejecta distribution and the early ($\lesssim 1\, \mathrm{d}$) light curves are largely unaffected.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Kaonic Hanbury-Brown–Twiss radii at 200 GeV and 5.02 TeV

Here, we use 3 + 1 dimensional quasiparticle anisotropic hydrodynamics (aHydroQP) to make predictions for kaon Hanbury-Brown–Twiss (HBT) radii in 200 GeV and 5.02 TeV heavy-ion collisions. Using previously determined aHydroQP parameters, we compute kaonic HBT radii and their ratios as a function of the mean transverse momentum of the pair k T . We first consider Au-Au collisions at 200 GeV, finding good agreement between aHydroQP predictions and experimental data up to k T ≈ 0.8 GeV. We then present predictions for kaonic HBT radii and their ratios in 5.02 TeV Pb-Pb collisions. Our aHydroQP predictions do not exhibit a clear k T scaling of the pion and kaon source radii, however, an approximate transverse mass m T scaling is observed, particularly at 200 GeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Characteristic momentum of Hydro + and a bound on the enhancement of the speed of sound near the QCD critical point

Near the critical point in the QCD phase diagram, hydrodynamics breaks down at a momentum where the frequency of the fastest hydrodynamic mode becomes comparable with the decay rate of the slowest nonhydrodynamic mode. Hydro+ was developed as a framework which extends the range of validity of hydrodynamics beyond that momentum value. This was achieved through coupling the hydrodynamic modes to the slowest nonhydrodynamic mode. In this work, analyzing the spectrum of linear perturbations in Hydro+, we find that a slow mode falls out of equilibrium if its momentum is greater than a characteristic momentum value. That characteristic momentum turns out to be set by the branch points of the dispersion relations. These branch points occur at the critical momenta of so-called spectral curves and are related to the radius of convergence of the derivative expansion. The existence of such a characteristic momentum scale suggests that a particular class of slow modes has no remarkable effect on the flow of the plasma. Based on these results, we show that there is an enhancement of the speed of sound due to the Hydro+ corrections compared to the known thermodynamic corrections near the critical point. This enhancement can not exceed a temperature-dependent upper bound which we derive near the critical point in the QCD phase diagram.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Observable signatures of initial state momentum anisotropies in nuclear collisions

In this work, we show that the correlation between the elliptic momentum anisotropy v 2 and the average transverse momentum [ p T ] at fixed multiplicity in small system nuclear collisions carries information on the origin of the observed momentum anisotropy. A calculation using a hybrid IP-Glasma+ music +UrQMD model that includes contributions from final state response to the initial geometry as well as initial state momentum anisotropies of the color glass condensate predicts a characteristic sign change of the correlator ρ ^ ( v 2 2 , [ p T ] ) as a function of charged particle multiplicity in p + Au and d + Au collisions at s = 200 GeV , and p + Pb collisions at s = 5.02 TeV . This sign change is absent in calculations without initial state momentum anisotropies. The model further predicts a qualitative difference between the centrality dependence of ρ ^ ( v 2 2 , [ p T ] ) in Au + Au collisions at s = 200 GeV and Pb + Pb collisions at s = 5.02 TeV , with only the latter showing a sign change in peripheral events. Predictions for O + O collisions at different collision energy show a similar behavior. Experimental observation of these distinct qualitative features of ρ ^ ( v 2 2 , [ p T ] ) in small and large systems would constitute strong evidence for the presence and importance of initial state momentum anisotropies predicted by the color glass condensate effective theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2010

Rapid Variability Generated at Relativistic Shocks Simulated by Particle-in- Cell Code

Plasma instabilities excited in collisionless shocks are responsible for particle acceleration. We have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic electronpositron jet propagating into an unmagnetized electron-positron plasma. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic-like shock structure. In the leading shock, electron density increases by a factor of about 3.5 in the simulation frame. Strong electromagnetic fields are generated in the trailing shock and provide an emission site. These magnetic fields contribute to the electrons transverse deflection behind the shock. We calculate the radiation from deflected electrons in the turbulent magnetic fields. Radiation from electrons near the trailing shock will be variable due to fluctuations of density and electromagnetic fields. The properties of this radiation may be important for rapid variability in relativistic jets such as AGN jets and blazars.

Nishikawal, Ken-Ichi↗
NASA NTRS2008

Modeling Flows Around Merging Black Hole Binaries

Coalescing massive black hole binaries are produced by the merger of galaxies. The final stages of the black hole coalescence produce strong gravitational radiation that can be detected by the space-borne LISA. In cases in which the black hole merger takes place in the presence of gas and magnetic fields, various types of electromagnetic signals may also be produced. Modeling such electromagnetic counterparts of the final merger requires evolving the behavior of both gas and fields in the strong-field regions around the black holes. We have taken a first step towards this problem by mapping the flow of pressureless matter in the dynamic, 3-D general relativistic spacetime around the merging black holes. We report on the results of these initial simulations and discuss their likely importance for future hydrodynamical simulations.

Centrella, Joan↗
NASA NTRS2009

Massive Black Hole Mergers: Can we see what LISA will hear?

Coalescing massive black hole binaries are formed when galaxies merge. The final stages of this coalescence produce strong gravitational wave signals that can be detected by the space-borne LISA. When the black holes merge in the presence of gas and magnetic fields, various types of electromagnetic signals may also be produced. Modeling such electromagnetic counterparts requires evolving the behavior of both gas and fields in the strong-field regions around the black holes. We have taken a first step towards this problem by mapping the flow of pressureless matter in the dynamic, 3-D general relativistic spacetime around the merging black holes. We report on the results of these initial simulations and discuss their likely importance for future hydrodynamical simulations.

Centrella, Joan↗
NASA NTRS2010

Simulation of Relativistic Shocks and Associated Self-Consistent Radiation

Plasma instabilities excited in collisionless shocks are responsible for particle acceleration. We have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized electron-positron plasma. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic-like shock structure. In the leading shock, electron density increases by a factor of about 3.5 in the simulation frame. Strong electromagnetic fields are generated in the trailing shock and provide an emission site. These magnetic fields contribute to the electrons transverse deflection behind the shock. We calculate the radiation from deflected electrons in the turbulent magnetic fields. The properties of this radiation may be important for understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets in general, and supernova remnants.

Nishikawa, K.-I.↗
NASA NTRS1993

Tidal disruptions by supermassive black holes - Hydrodynamic evolution of stars on a Schwarzschild background

We present a three-dimensional numerical study of tidal disruption of a main-sequence star by a supermassive black hole. The simulations include general relativistic effects which are important in this regime. We analyze stars in a marginally bound orbit around the black hole with pericentric separation of a few Schwarzschild radii. We show that during a close passage, as a result of relativistic effects analogous to the perihelion shift, the trajectories of the debris of the star fan out into a crescent-like shape centered on the black hole. We also discuss the increase of the central density of the star as it approaches pericentric distance, the fraction of the debris accreted by the hole, its accretion rate, the distribution of debris orbits bound to the hole, and the velocity of unbound ejected material. We compare these results with the disruption of the star by a Newtonian point mass.

Laguna, Pablo↗
DOE OSTI2026

Characterizing radial flow fluctuations in relativistic heavy-ion collisions at top RHIC and LHC energies

This study presents a systematic investigation of the transverse-momentum differential radial flow fluctuations observable 𝑣 0 ⁡ (𝑝 𝑇 ) in relativistic heavy-ion collisions at top Relativistic Heavy Ion Collider ($\sqrt{s_{NN}}$ = 200 GeV) and Large Hadron Collider ($\sqrt{s_{NN}}$ = 2.76 and 5.02 TeV) energies. Using a multistage hydrodynamic model, this study assesses the sensitivity of 𝑣 0 ⁡(𝑝 𝑇 ) to a wide range of physical effects, including bulk and shear viscosities, off-equilibrium corrections at particlization, the presence of a hadronic afterburner, and the nucleon size in the initial conditions. By employing complementary rescaling strategies, this study demonstrates how different physical effects leave distinct imprints on the shape of 𝑣 0 ⁡(𝑝 𝑇 ). A combined double-rescaling of 𝑣 0 ⁡(𝑝 𝑇 )/𝑣 0 versus 𝑝 𝑇 /⟨𝑝 𝑇 ⟩ reveals a universality across a wide range of energies and model assumptions in the low-𝑝 𝑇 regime, a robust signature of collective behavior. This allows us to disentangle the universal dynamics of the bulk medium from model-specific features that emerge at higher 𝑝 𝑇 . These results establish 𝑣 0⁡ (𝑝 𝑇 ) as a powerful and complementary observable for constraining quark-gluon plasma transport properties and initial-state granularity, offering a unique probe of the created QCD medium.

Du, Lipei [University of California, Berkeley, CA ↗
NASA NTRS1981

Compton losses, Compton rockets

The radiation force on a relativistic plasma is shown to accelerate the plasma to relativistic bulk velocities under certain conditions, thus demonstrating that relativistic bulk motion not only alleviates Compton losses, but can also result from such losses. In order to isolate Compton-rocket effects from hydrodynamical ones, it is assumed that the plasma in its instantaneous rest frame has a temporally constant density and is spatially uniform.

Cheng, A. Y. S.↗
NASA NTRS2009

Simulation of Relativistic Shocks and Associated Radiation from Turbulent Magnetic Fields

Plasma instabilities excited in collisionless shocks are responsible for particle acceleration. We have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized electron-positron plasma. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic-like shock structure. In the leading shock, electron density increases by a factor of about 3.5 in the simulation frame. Strong electromagnetic fields are generated in the trailing shock and provide an emission site. These magnetic fields contribute to the electron's transverse deflection behind the shock. The jitter'' radiation from deflected electrons in turbulent magnetic fields has different properties than synchrotron radiation, which is calculated in a uniform magnetic field. This jitter radiation may be important for understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets in general, and supernova remnants. New spectra based on simulations will be presented.

Nishikawa, K.-I.↗
DOE OSTI2026

thornado+FLASH-X: A Hybrid Discontinuous Galerkin–Implicit-explicit and Finite-volume Framework for Neutrino-radiation Hydrodynamics in Core-collapse Supernovae

We present neutrino-transport algorithms implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado) and their coupling to self-gravitating hydrodynamics within the adaptive mesh refinement–based multiphysics simulation framework FLASH-X. thornado, developed primarily for simulations of core-collapse supernovae (CCSNe), employs a spectral, six-species two-moment formulation with algebraic closure and special-relativistic observer corrections accurate to $\mathcal{O}(v/c)$, and uses discontinuous Galerkin (DG) methods for phase-space discretization combined with implicit-explicit time stepping. A key development is a nonlinear neutrino–matter coupling algorithm based on nested fixed-point iteration with Anderson acceleration, enabling fully implicit treatment of collisional processes, including energy-coupling interactions such as neutrino–electron scattering and pair production. Coupling to finite-volume (FV) hydrodynamics is achieved through a hybrid DG-FV representation of the fluid variables and operator-split evolution within FLASH-X. The implementation is verified using basic transport tests with idealized opacities and relaxation and deleptonization problems with tabulated microphysics. Spherically symmetric CCSN simulations demonstrate accuracy and robustness of the coupled scheme, including close agreement with the CCSN simulation code Chimera. An axisymmetric CCSN simulation further demonstrates the viability of DG-based neutrino transport for multidimensional supernova modeling within FLASH-X. thornado’s neutrino-transport solver is GPU-enabled using OpenMP offloading or OpenACC, and all CCSN applications included in this work use the GPU implementation. Together, these results establish a foundation for future enhancements in physics fidelity, numerical algorithms, and computational performance, for increasingly realistic large-scale CCSN simulations.

Endeve, Eirik [Oak Ridge National Laboratory (ORNL↗
DOE OSTI2026

Quark flavor equilibration of the quark-gluon plasma

The early stage of a heavy-ion collision is marked by rapid entropy production and the transition from a gluon saturated initial condition to a plasma of quarks and gluons that evolves hydrodynamically. However, during the early times of the hydrodynamic evolution, the chemical composition of the QCD medium is still largely unknown. We present a study of quark chemical equilibration in the (Q)GP using a novel model of viscous hydrodynamic evolution in partial chemical equilibrium. Motivated by the success of gluon saturated initial condition models, we initialize the QCD medium as a completely gluon dominated state. Local quark production during the hydrodynamic phase is then simulated through the evolution of time-dependent fugacities for each independent quark flavor, with the timescales set as free parameters to compare different rates of equilibration. We present the results of complete heavy-ion collision simulations using this partial chemical equilibrium model, and show the effects on hadronic and electromagnetic observables. In particular, we show that the development of flow is sensitive to the equilibration timescale, providing an empirical way to probe the chemical equilibration of the QCD medium.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1994

Analysis of the early spectra and light curve of SN 1987A

Numerical modeling of supernova spectra, light curves, and hydrodynamics requires physical inputs, numerical techniques, approximations, and assumptions which must be thoroughly understood in order to study the details of supernova explosions. Here, we discuss some of these in the context of the early evolution of supernova 1987A. Gray radiation-hydrodynamics is used to calculate the bolometric light curve and the hydrodynamic evolution of the supernova. Synthetic spectra are then obtained for the resulting density and velocity structure. The spectrum calculations are performed using a special-relativistic treatment of the radiative transfer equation in the comoving frame, line blanketing by about 10(exp 5) spectral lines, and departures from local thermodynamic equilibrium (LTE) for H I, He I, Mg II, and Ca II. We find that we are able to simultaneously fit the early light curve and spectra reasonably well, using a progenitor model from Arnett (1991a), without fine-tuning the free parameters. Temperature structures and radiative equilibrium, non-LTE effects, homologous expansion, and mean opacities are discussed.

Hauschildt, Peter H.↗
DOE OSTI2024

Linear Boltzmann transport for jet propagation in the quark-gluon plasma: Inelastic processes and jet modification

A linear Boltzmann transport (LBT) Monte Carlo model has been developed to describe jet propagation and interaction with the quark-gluon plasma (QGP) in relativistic heavy-ion collisions. A complete set of elastic-scattering processes and medium-induced gluon emissions based on the higher-twist formalism are incorporated for both jet shower and medium recoil partons. It has been employed to describe experimental data on large transverse momentum hadron and jet spectra, correlation and jet substructures in high-energy heavy-ion collisions. Here, we document in detail the structure of the model and validation of the Monte Carlo implementations of the physics processes in LBT, in particular, the inelastic process of medium-induced gluon radiation. We carry out a comprehensive examination of the jet-medium interaction as implemented in LBT through energy loss and momentum broadening of a single hard parton, the energy and transverse momentum transfer from leading partons to medium-induced gluons and jet-induced medium excitation, and medium modification of reconstructed jets in a static and uniform medium. With realistic and event-by-event hydrodynamic medium in heavy-ion collisions, we compute and compare with experimental data on the jet cone-size dependence of the single inclusive jet suppression at both the BNL Relativistic Heavy-Ion Collider (RHIC) and the CERN Large Hadron Collider (LHC), the dijet asymmetry at the LHC and γ-jet correlation at RHIC. Effects of medium-induced gluon emissions and jet-induced medium excitation on jet observables are systematically examined. Rescatterings of the radiated gluons and recoil partons with the QGP are found essential to account for the enhancement of soft particle yield toward the edge of the jet cone.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗