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.
QCD thermalization: Ab initio approaches and interdisciplinary connections
Heavy-ion collisions at BNL's Relativistic Heavy Ion Collider and CERN's Large Hadron Collider provide strong evidence for the formation of a quark-gluon plasma, with temperatures extracted from relativistic viscous hydrodynamic simulations shown to be well above the transition temperature from hadron matter. Outstanding problems in QCD include how the strongly correlated quark gluon matter forms in a heavy-ion collision, its properties off-equilibrium, and the thermalization process in the plasma. In this work, we review the theoretical progress in this field in weak coupling QCD effective field theories and in strong coupling holographic approaches based on gauge-gravity duality. We outline the interdisciplinary connections of different stages of the thermalization process to non equilibrium dynamics in other systems across energy scales ranging from inflationary cosmology, to strong field QED, to ultracold atomic gases, with emphasis on the universal dynamics of non-thermal and hydrodynamic attractors. We survey measurements in heavy-ion collisions that are sensitive to the early non-equilibrium stages of the collision and discuss the potential for future measurements. We summarize the current state-of-the art in thermalization studies and identify promising avenues for further progress.
AthenaK: A Performance-portable Version of the Athena++ Adaptive Mesh Refinement Framework
We describe AthenaK: a new implementation of the Athena++ block-based adaptive mesh refinement framework using the Kokkos programming model. Finite volume methods for Newtonian, special relativistic, and general relativistic (GR) hydrodynamics and magnetohydrodynamics (MHD), and GR-radiation hydrodynamics and MHD, as well as a module for evolving Lagrangian tracer or charged test particles (e.g., cosmic rays) are implemented using the framework. In two companion papers, we describe (1) a new solver for the Einstein equations based on the Z4c formalism, and (2) a GRMHD solver in dynamical spacetimes also implemented using the framework, enabling new applications in numerical relativity. By adopting Kokkos, the code can be run on virtually any hardware, including CPUs, GPUs from multiple vendors, and emerging Advanced RISC Machine processors. AthenaK shows excellent performance and weak scaling, achieving over 1 billion cell updates per second for hydrodynamics in three dimensions on a single NVIDIA Grace Hopper processor. It does this with a typical parallel efficiency of 80% on 65,536 AMD GPUs on the OLCF Frontier system. Such performance portability enables AthenaK to leverage modern exascale computing systems for challenging applications in astrophysical fluid dynamics, numerical relativity, and multimessenger astrophysics.
X-rays from quasars and active galaxies
Some of the most important theoretical ideas and observations for quasars and the nuclei of active galaxies are reviewed, and areas of future research are suggested. Emphasis is on the nature of the power source, the radiation processes, and the mechanism for formation and collimation of jets. Phenomena that produce X-rays are of particular concern. Particular topics discussed are the observed and expected time variabilility, the gas supply mechanisms and luminosity evolution, thermal and nonthermal radiation processes, observed and theoretical spectra, criteria for thermalization of electrons and ions, effects of electron-positron pairs on relativistic plasmas, hydrodynamic, electrodynamic and inertial methods for producing and confining jets. The paper concludes with a list of needed observations.
Relativistic hadrons and the origin of relativistic outflows in active galactic nuclei
We examine the hydrodynamic origin of relativistic outflows in active galactic nuclei (AGN). Specifically, we propose that the presence of a population of relativistic hadrons in the AGN 'central engine' and the associated neutron production suffices to produce outflows which under rather general conditions could be relativistic. The main such condition is that the size of the neutron production region be larger than the neutron flight path tau(sub n) approximately 3 x 10(exp 13) cm. This condition guarantees that the mean energy per particle in the proton fluid, resulting from the decay of the neutrons outside their production region, be greater than the proton rest mass. The expansion of this fluid can then lead naturally to a relativistic outflow by conversion of its internal energy to directed motion. We follow the development of such flows by solving the mass, energy as well as the kinetic equation for the proton gas in steady state, taking into account the source terms due to compute accurately the adiabatic index of the expanding gas, and in conjunction with Bernoulli's equation the detailed evolution of the bulk Lorentz factor. We further examine the role of large-scale magnetic fields in confining these outflows to produce the jets observed at larger scales.
A Magnetohydrodynamic Boost for Relativistic Jets
We performed relativistic magnetohydrodynamic simulations of the hydrodynamic boosting mechanism for relativistic jets explored by Aloy & Rezzolla (2006) using the RAISHIN code. Simulation results show that the presence of a magnetic field changes the properties of the shock interface between the tenuous, overpressured jet (V^z j) flowing tangentially to a dense external medium. We find that magnetic fields can lead to more efficient acceleration of the jet, in comparison to the pure-hydrodynamic case. A "poloidal" magnetic field (B^z), tangent to the interface and parallel to the jet flow, produces both a stronger outward moving shock and a stronger inward moving rarefaction wave. This leads to a large velocity component normal to the interface in addition to acceleration tangent to the interface, and the jet is thus accelerated to larger Lorentz factors than those obtained in the pure-hydrodynamic case. Likewise, a strong "toroidal" magnetic field (B^y), tangent to the interface but perpendicular to the jet flow, also leads to stronger acceleration tangent to the shock interface relative to the pure-hydrodynamic case. Thus. the presence and relative orientation of a magnetic field in relativistic jets can significant modify the hydrodynamic boost mechanism studied by Aloy & Rezzolla (2006).
A Magnetohydrodynamic Boost for Relativistic Jets
We have performed relativistic magnetohydrodynamic simulations of the hydrodynamic boosting mechanism for relativistic jets explored by Aloy & Rezzolla (2006) using the RAISHIN code. Simulation results show that the presence of a magnetic field may change the properties of the shock interface between the tenuous, overpressured jet (V(sub j) (sup z)) flowing tangentially to a dense external medium. Magnetic fields can lead to more efficient acceleration of the jet, in comparison to the pure-hydrodynamic case. A poloidal magnetic field (B(sup z)), tangent to the interface and parallel to the jet flow, produces both a stronger outward moving shock and inward moving rarefaction wave. This leads to a large velocity component normal to the interface in addition to acceleration tangent to the interface, and the jet is thus accelerated to a larger Lorentz factors than those obtained in the pure-hydrodynamic case. In contrast, a strong toroidal magnetic field (B(sup y)), tangent to the interface but perpendicular to the jet flow, also leads to stronger acceleration tangent to the shock interface relative to the pure-hydrodynamic case, but to a lesser extent than found for the poloidal case due to the fact that the velocity component normal to the shock interface is now much smaller. Overall, the acceleration efficiency in the toroidal case is less than that of the poloidal case but both geometries still result in higher Lorentz factors than the pure-hydrodynamic case. Thus, the presence and relative orientation of a magnetic field in relativistic jets can have a significant influence on the hydrodynamic boost mechanism studied by Aloy & Rezzolla (2006).
Transport-based initial conditions for heavy-ion collisions at finite densities
Here, we employ the SMASH transport model to provide event-by-event initial conditions for the energy-momentum tensor and conserved charge currents in hydrodynamic simulations of relativistic heavy-ion collisions. We study the fluctuations and dynamical evolution of three conserved charge currents (net baryon, net electric charges, and net strangeness) with a four-dimensional lattice-QCD-based equation of state, NEOS-4D, in the hydrodynamic phase. Out-of-equilibrium corrections at the particlization are generalized to finite densities to ensure the conservation of energy, momentum, and the three types of charges. These theoretical developments are integrated within the X-SCAPE code as a unified framework for studying the nuclear matter properties in the Beam Energy Scan program.
Relativistic Guiding-Center Motion: Action Principle, Kinetic Theory and Hydrodynamics
We treat the guiding-center dynamics in a varying external Maxwell field using a relativistically covariant action principle which reproduces the known Vandervoort expression for the drift velocity and extends it to curved spacetime. We derive the corresponding kinetic theory and ideal hydrodynamic theory. In contrast to conventional five-equation hydrodynamics, the guiding-center hydrodynamics needs only three equations due to a constraint on the motion across magnetic field. Furthermore, we argue that such a hydrodynamics is applicable to strongly coupled plasmas where kinetic theory fails.
Freezing out fluctuations in Hydro + near the QCD critical point
Not Available
Perturbative approaches in relativistic kinetic theory and the emergence of first-order hydrodynamics
Hydrodynamics can be formulated in terms of a perturbative series in derivatives of the temperature, chemical potential, and flow velocity around an equilibrium state. Different formulations for this series have been proposed over the years, which consequently led to the development of various hydrodynamic theories. In this work, we discuss the relativistic generalizations of the perturbative expansions put forward by Chapman and Enskog, and Hilbert, using general matching conditions in kinetic theory. This allows us to describe, in a comprehensive way, how different out-of-equilibrium definitions for the hydrodynamic fields affect the development of the hydrodynamic perturbative series. We provide a perturbative method for systematically deriving the hydrodynamic formulation recently proposed by Bemfica, Disconzi, Noronha, and Kovtun (BDNK) from relativistic kinetic theory. The various transport coefficients that appear in BDNK (at first-order) are explicitly computed using a new formulation of the relaxation time approximation for the Boltzmann equation. Assuming Bjorken flow, we also determine the hydrodynamic attractors of BDNK theory and compare the overall hydrodynamic evolution obtained using this formulation with that generated by the Israel-Stewart equations of motion and also kinetic theory.
Dissipative magnetohydrodynamics for nonresistive relativistic plasmas: An implicit second-order flux-conservative formulation with stiff relaxation
Based on a 14-moment closure for nonresistive (general-) relativistic viscous plasmas, here we describe a new numerical scheme that is able to handle all first-order dissipative effects (heat conduction, bulk and shear viscosities), as well the anisotropies induced by the presence of magnetic fields. The latter is parametrized in terms of a thermal gyrofrequency or, equivalently, a thermal Larmor radius and allows to correctly capture the thermal Hall effect. By solving an extended Israel-Stewart-like system for the dissipative quantities that enforces algebraic constraints via stiff-relaxation, we are able to cast all first-order dissipative terms in flux-divergence form. This allows us to apply traditional high-resolution shock capturing methods to the equations, making the system suitable for the numerical study of highly turbulent flows. We present several numerical tests to assess the robustness of our numerical scheme in flat spacetime. The 14-moment closure can seamlessly interpolate between the highly collisional limit found in neutron star mergers, and the highly anisotropic limit of relativistic Braginskii magnetohydrodynamics appropriate for weakly collisional plasmas in black-hole accretion problems. We believe that this new formulation and numerical scheme will be useful for a broad class of relativistic magnetized flows.
Skewness of mean transverse momentum fluctuations in heavy-ion collisions
We propose the skewness of mean transverse momentum, $\langle{p_t}\rangle$, fluctuations as a fine probe of hydrodynamic behavior in relativistic nuclear collisions. We describe how the skewness of the $\langle{p_t}\rangle$ distribution can be analyzed experimentally, and we use hydrodynamic simulations to predict its value. We predict in particular that $\langle{p_t}\rangle$ fluctuations have positive skew, which is significantly larger than if particles were emitted independently. We elucidate the origin of this result by deriving generic formulas relating the fluctuations of $\langle{p_t}\rangle$ to the fluctuations of the early-time thermodynamic quantities. We postulate that the large positive skewness of $\langle{p_t}\rangle$ fluctuations is a generic prediction of hydrodynamic models.
QCD equation of state at finite chemical potentials for relativistic nuclear collisions
In this work, we review the equation of state of QCD matter at finite densities. We discuss the construction of the equation of state with net baryon number, electric charge, and strangeness using the results of lattice QCD simulations and hadron resonance gas models. Its application to the hydrodynamic analyses of relativistic nuclear collisions suggests that the interplay of multiple conserved charges is important in the quantitative understanding of the dense nuclear matter created at lower beam energies. Several different models of the QCD equation of state are discussed for comparison.
Smooth Particle Hydrodynamics for Neutron Star Oscillation with General Relativistic effects
Explore the source record for details and available documents.
X-Ray Spectra from MHD Simulations of Accreting Black Holes
We present new global calculations of X-ray spectra from fully relativistic magneto-hydrodynamic (MHO) simulations of black hole (BH) accretion disks. With a self consistent radiative transfer code including Compton scattering and returning radiation, we can reproduce the predominant spectral features seen in decades of X-ray observations of stellar-mass BHs: a broad thermal peak around 1 keV, power-law continuum up to >100 keV, and a relativistically broadened iron fluorescent line. By varying the mass accretion rate, different spectral states naturally emerge: thermal-dominant, steep power-law, and low/hard. In addition to the spectral features, we briefly discuss applications to X-ray timing and polarization.
Evidence of a toroidal magnetic field in the core of 3C 84
The spatial scales of relativistic radio jets, probed by relativistic magneto-hydrodynamic (RMHD) jet launching simulations and by most very long baseline interferometry (VLBI) observations differ by an order of magnitude. Bridging the gap between these RMHD simulations and VLBI observations requires selecting nearby active galactic nuclei (AGN), the parsec-scale region of which can be resolved. The radio source 3C 84 is a nearby bright AGN fulfilling the necessary requirements: it is launching a powerful, relativistic jet powered by a central supermassive black hole, while also being very bright. Using 22 GHz globe-spanning VLBI measurements of 3C 84 we studied its sub-parsec region in both total intensity and linear polarisation to explore the properties of this jet, with a linear resolution of ~0.1 parsec. We tested different simulation set-ups by altering the bulk Lorentz factor Γ of the jet, as well as the magnetic field configuration (toroidal, poloidal, helical). We confirm the persistence of a limb brightened structure, which reaches deep into the sub-parsec region. The corresponding electric vector position angles (EVPAs) follow the bulk jet flow inside but tend to be orthogonal to it near the edges. Our state-of-the-art RMHD simulations show that this geometry is consistent with a spine-sheath model, associated with a mildly relativistic flow and a toroidal magnetic field configuration.
Massive stars in quasi-static equilibrium.
Quasi-static equilibrium in massive stars during gravitational contraction including discussion of general relativistic and classical hydrodynamic results