Kinetic schemes for the relativistic hydrodynamics

Hyperbolic systems, non-relativistic and relativistic Euler equations, kinetic schemes, conservation laws, discontinuous solutions, high order accuracy

Relativistic Hydrodynamics in Isotropic and Anisotropic Cosmological Backgrounds

A seguito di recenti risultati nel campo dell’Astrofisica, con questo elaborato ci si propone di approfondire il ruolo della viscosità in una classe di modelli cosmologici. La discussione ha lo scopo di generalizzare delle tecniche applicate con successo nello studio dei fluidi ideali a sistemi dove anche la viscosità é un parametro che necessita di essere largamente preso in considerazione. Nello specifico, ci si serve di strumenti computazionali e geometrici e di teoria delle equazioni differenziali alle derivate parziali per comprendere, da un punto di vista matematico, come associare le Equazioni di Einstein ai fluidi in specifici background cosmologici. Questo elaborato parte dallo studio dei fluidi ideali in un background di tipo FLRW al fine di generalizzarlo a fluidi viscosi che scorrono in un background di tipo Bianchi I. Infine si indicano alcuni problemi ancora aperti relativi al caso dissipativo e le possibili strategie risolutive di tali questioni.

Accelerating solutions of perfect fluid hydrodynamics for initial energy density and life-time measurements in heavy ion collisions

A new class of accelerating, exact, explicit and simple solutions of relativistic hydrodynamics is presented. Since these new solutions yield a finite rapidity distribution, they lead to an advanced estimate of the initial energy density and life-time of high energy heavy ion collisions. Accelerating solutions are also given for spherical expansions in arbitrary number of spatial dimensions.

Korteveg-de Vries solitons in a cold quark-gluon plasma

The relativistic heavy ion program developed at RHIC and now at LHC motivated a deeper study of the properties of the quark-gluon plasma (QGP) and, in particular, the study of perturbations in this kind of plasma. We are interested on the time evolution of perturbations in the baryon and energy densities. If a localized pulse in baryon density could propagate throughout the QGP for long distances preserving its shape and without loosing localization, this could have interesting consequences for relativistic heavy ion physics and for astrophysics. A mathematical way to prove that this can happen is to derive (under certain conditions) from the hydrodynamical equations of the QGP a Korteveg-de Vries (KdV) equation. The solution of this equation describes the propagation of a KdV soliton. The derivation of the KdV equation depends crucially on the equation of state (EOS) of the QGP. The use of the simple MIT bag model EOS does not lead to KdV solitons. Recently we have developed an EOS for the QGP which includes both perturbative and nonperturbative corrections to the MIT one and is still simple enough to allow for analytical manipulations. With this EOS we were able to derive a KdV equation for the cold QGP.

Nonlinear waves in a quark gluon plasma

We study the propagation of perturbations in the quark gluon plasma. This subject has been addressed in other works and in most of the theoretical descriptions of this phenomenon the hydrodynamic equations have been linearized for simplicity. We propose an alternative approach, also based on hydrodynamics but taking into account the nonlinear terms of the equations. We show that these terms may lead to localized waves or even solitons. We use a simple equation of state for the QGP and expand the hydrodynamic equations around equilibrium configurations. The resulting differential equations describe the propagation of perturbations in the energy density. We solve them numerically and find that localized perturbations can propagate for long distances in the plasma. Under certain conditions our solutions mimic the propagation of Korteweg-de Vries solitons.

BREAKING WAVES IN AN IDEAL QUARK GLUON PLASMA

We study the propagation of perturbations in the energy density in a quark gluon plasma. Expanding the Euler and continuity equations of relativistic hydrodynamics around equilibrium configurations we obtain a nonlinear differential equation called the breaking wave equation. We solve it numerically and follow the time-evolution of initially localized pulses. We find that, quite unexpectedly, these pulses live for a very long time (compared to the reaction time-scales) before breaking. In practice, they mimick the Korteweg-de Vries solitons. Their existence may have some observable consequences.

Sound waves and solitons in hot and dense nuclear matter

Assuming that nuclear matter can be treated as a perfect fluid, we study the propagation of perturbations in the baryon density. The equation of state is derived from a relativistic mean field model, which is a variant of the non-linear Walecka model. The expansion of the Euler and continuity equations of relativistic hydrodynamics around equilibrium configurations leads to differential equations for the density perturbation. We solve them numerically for linear and spherical perturbations and follow the propagation of the initial pulses. For linear perturbations we find single soliton solutions and solutions with one or more solitons followed by ""radiation"". Depending on the equation of state a strong damping may occur. We consider also the evolution of perturbations in a medium without dispersive effects. In this case we observe the formation and breaking of shock waves. We study all these equations also for matter at finite temperature. Our results may be relevant for the analysis of RHIC data. They suggest that the shock waves formed in the quark gluon plasma phase may survive and propagate in the hadronic phase. (C) 2009 Elseiver. B.V. All rights reserved.

Shock wave formation in hot nuclear matter

Assuming that nuclear matter can be treated as a perfect fluid, we study the propagation of perturbations in the baryon density at high temperature. The equation of state is derived from the non-linear Walecka model. The expansion of the Euler and continuity equations of relativistic hydrodynamics around equilibrium configurations lead to the breaking wave equation for the density perturbation. We solve it numerically for this perturbation and follow the propagation of the initial pulses.

On the radial expansion of tubular structures in a quark-gluon plasma

We study the radial expansion of cylindrical tubes in a hot QGP. These tubes are treated as perturbations in the energy density of the system which is formed in heavy ion collisions at RHIC and LHC. We start from the equations of relativistic hydrodynamics in two spatial dimensions and cylindrical symmetry and perform an expansion of these equations in a small parameter, conserving the nonlinearity of the hydrodynamical formalism. We consider both ideal and viscous fluids and the latter are studied with a relativistic Navier-Stokes equation. We use the equation of state of the MIT bag model. In the case of ideal fluids we obtain a breaking wave equation for the energy density fluctuation, which is then solved numerically. We also show that, under certain assumptions, perturbations in a relativistic viscous fluid are governed by the Burgers equation. We estimate the typical expansion time of the tubes. (C) 2012 Elsevier B.V. All rights reserved.

FLUCTUATIONS AND INITIAL STATE GRANULARITY IN HEAVY ION COLLISIONS AND THEIR EFFECTS ON OBSERVABLES FROM HYDRODYNAMICS

A comparison is made between results obtained using smooth initial conditions and event-by-event initial conditions in the hydrodynamical description of relativistic nuclear collisions. Some new results on directed flow are also included.

Topology Studies of Hydrodynamics Using Two-Particle Correlation Analysis

The effects of fluctuating initial conditions are studied in the context of relativistic heavy ion collisions where a rapidly evolving system is formed. Two-particle correlation analysis is applied to events generated with the NEXSPHERIO hydrodynamic code, starting with fluctuating nonsmooth initial conditions (IC). The results show that the nonsmoothness in the IC survives the hydroevolution and can be seen as topological features of the angular correlation function of the particles emerging from the evolving system. A long range correlation is observed in the longitudinal direction and in the azimuthal direction a double peak structure is observed in the opposite direction to the trigger particle. This analysis provides clear evidence that these are signatures of the combined effect of tubular structures present in the IC and the proceeding collective dynamics of the hot and dense medium.

KINETICS VERSUS HYDRODYNAMICS: GENERALIZATION OF LANDAU/COOPER-FRYE PRESCRIPTION FOR FREEZE-OUT

The problem of spectra formation in hydrodynamic approach to A + A collisions is considered within the Boltzmann equations. It is shown analytically and illustrated by numerical calculations that the particle momentum spectra can be presented in the Cooper-R-ye form despite freeze-out is not sharp and has the finite temporal width. The latter is equal to the inverse of the particle collision rate at points (t(sigma) (r, p), r) of the maximal emission at a fixed momentum p. The set of these points forms the hypersurfaces t(sigma)(r,p) which strongly depend on the values of p and typically do not enclose completely the initially dense matter. This is an important difference from the standard Cooper-Frye prescription (CFp), with a common freeze-out hypersurface for all p, that affects significantly the predicted spectra. Also, the well known problem of CFp as for negative contributions to the spectra from non-space-like parts of the freeze-out hypersurface is naturally eliminated in this improved prescription.

Directed flow at midrapidity in event-by-event hydrodynamics

Fluctuations in the initial geometry of a nucleus-nucleus collision have been recently shown to result in a new type of directed flow (v(1)) that, unlike the usual directed flow, is also present at midrapidity. We compute this new v(1) versus transverse momentum and centrality for Au-Au collisions at RHIC using the hydrodynamic code NeXSPheRIO. We find that the event plane of v(1) is correlated with the angle of the initial dipole of the distribution, as predicted, though with a large dispersion. It is uncorrelated with the reaction plane. Our results are in excellent agreement with results inferred from STAR correlation data.

Hydro-kinetic approach to relativistic heavy ion collisions

We develop a combined hydro-kinetic approach which incorporates a hydrodynamical expansion of the systems formed in A + A collisions and their dynamical decoupling described by escape probabilities. The method corresponds to a generalized relaxation time (tau(rel)) approximation for the Boltzmann equation applied to inhomogeneous expanding systems; at small tau(rel) it also allows one to catch the viscous effects in hadronic component-hadron-resonance gas. We demonstrate how the approximation of sudden freeze-out can be obtained within this dynamical picture of continuous emission and find that hypersurfaces, corresponding to a sharp freeze-out limit, are momentum dependent. The pion m(T) spectra are computed in the developed hydro-kinetic model, and compared with those obtained from ideal hydrodynamics with the Cooper-Frye isothermal prescription. Our results indicate that there does not exist a universal freeze-out temperature for pions with different momenta, and support an earlier decoupling of higher p(T) particles. By performing numerical simulations for various initial conditions and equations of state we identify several characteristic features of the bulk QCD matter evolution preferred in view of the current analysis of heavy ion collisions at RHIC energies.

Shear viscosity from thermal fluctuations in relativistic conformal fluid dynamics

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)