Transcript Hydrodynamics and Flow

Intensive Lecture
YITP, December 11th, 2008
Relativistic Ideal and
Viscous Hydrodynamics
Tetsufumi Hirano
Department of Physics
The University of Tokyo
TH, N. van der Kolk, A. Bilandzic, arXiv:0808.2684[nucl-th];
to be published in “Springer Lecture Note in Physics”.
Plan of this Lecture
1st Day
• Hydrodynamics in Heavy Ion Collisions
• Collective flow
• Dynamical Modeling of heavy ion
collisions (seminar)
2nd Day
• Formalism of relativistic ideal/viscous
hydrodynamics
• Bjorken’s scaling solution with viscosity
• Effect of viscosity on particle spectra
(discussion)
PART 3
Formalism of
relativistic
ideal/viscous
hydrodynamics
Relativistic Hydrodynamics
Equations of motion in relativistic hydrodynamics
Energy-momentum conservation
Energy-Momentum tensor
Current conservation
The i-th conserved current
In H.I.C., Nim = NBm (net baryon current)
Tensor/Vector Decomposition
Tensor decomposition with a given time-like
and normalized four-vector um
where,
“Projection” Tensor/Vector
•um is local four flow velocity. More precise
meaning will be given later.
•um is perpendicular to Dmn.
•Local rest frame (LRF):
•Naively speaking, um (Dmn) picks up time(space-)like component(s).
Intuitive Picture of Projection
time like
flow vector
field
Decomposition of Tmu
:Energy density
:Energy (Heat) current
:Shear stress tensor
:(Hydrostatic+bulk) pressure
P = Ps + P
<…>: Symmetric, traceless and transverse to um
& un
Decomposition of Nm
:charge density
:charge current
Q. Count the number of unknowns
in the above decomposition and confirm
that it is 10(Tmn)+4k(Nim).
Here k is the number of independent currents.
Note: If you consider um as independent variables,
you need additional constraint for them.
If you also consider Ps as an independent
variable, you need the equation of state Ps=Ps(e,n).
Ideal and Dissipative Parts
Energy Momentum tensor
Charge current
Ideal part
Dissipative part
Meaning of um
um is four-velocity of “flow”. What kind of flow?
Two major definitions of flow are
1. Flow of energy (Landau)
2. Flow of conserved charge (Eckart)
Meaning of um (contd.)
Landau
(Wm=0, uLmVm=0)
Vm
uLm
Eckart
(Vm=0,uEmWm=0)
Wm
uEm
Just a choice of local reference frame.
Landau frame might be relevant in H.I.C.
Relation btw.
Landau and Eckart
Relation btw.
Landau and Eckart (contd.)
Entropy Conservation
in Ideal Hydrodynamics
Neglect “dissipative part” of energy momentum
tensor to obtain “ideal hydrodynamics”.
Therefore,
Q. Derive the above equation.
Entropy Current
Assumption (1st order theory):
Non-equilibrium entropy current vector has
linear dissipative term(s) constructed from
(Vm, P, pmn, (um)).
(Practical) Assumption:
•Landau frame (omitting subscript “L”).
•No charge in the system.
Thus, a = 0 since Nm = 0, Wm = 0 since considering
the Landau frame, and g = 0 since um Sm should
be maximum in equilibrium (stability condition).
The 2nd Law of Thermodynamics
and Constitutive Equations
The 2nd thermodynamic law tells us
Q. Check the above calculation.
Constitutive Equations (contd.)
Newton
Stokes
Thermodynamic
force
Transport
coefficient
tensor
shear
scalar
bulk
“Current”
Equation of Motion
: Lagrange (substantial)
derivative
: Expansion scalar (Divergence)
Equation of Motion (contd.)
Q. Derive the above equations of motion
from energy-momentum conservation.
Note: We have not used the constitutive
equations to obtain the equations of motion.
Intuitive Interpretation of EoM
Change of volume
•Dilution
•Compression
Work done by pressure
Production of entropy
Conserved Current Case
Lessons from (Non-Relativistic)
Navier-Stokes Equation
Assuming incompressible fluids such that
, Navier-Stokes eq. becomes
Source of flow
(pressure gradient)
Diffusion of flow
(Kinematic viscosity, h/r,
plays a role of diffusion
constant.)
Final flow velocity comes from interplay between
these two effects.
Generation of Flow
P
Pressure gradient
Expand
Source of flow
Expand
x
 Flow phenomena
are important in H.I.C
to understand EOS
Diffusion of Flow
Heat equation
(k: heat conductivity
~diffusion constant)
For illustrative purpose, one discretizes the
equation in (2+1)D space:
Diffusion ~ Averaging
~ Smoothing
R.H.S. of descretized heat/diffusion eq.
y
y
j
subtract
i
j
i
x
x
Suppose Ti,j is larger (smaller) than an average value
around the site, R.H.S. becomes negative (positive).
2nd derivative w.r.t. coordinates  Smoothing
Shear Viscosity Reduces
Flow Difference
Shear flow
(gradient of flow) Smoothing of flow
Next time step
Microscopic interpretation can be made.
Net momentum flow in space-like direction.
 Towards entropy maximum state.
Necessity of Relaxation Time
Non-relativistic case (Cattaneo(1948))
Balance eq.:
Constitutive eq.:
Fourier’s law
t : “relaxation time”
Parabolic equation (heat equation)
ACAUSAL!
Finite t
Hyperbolic equation (telegraph equation)
Heat Kernel
x
perturbation
on top of
background
x
Heat
transportation
causality
Instability
• The 1st order equation is not only
acausal but also unstable under
small perturbation on a moving
back-ground. W.A.Hiscock and L.Lindblom,
PRD31,725(1985).
• For particle frame with new EoM, see
K.Tsumura and T.Kunihiro, PLB668, 425(2008).
• For a possible relation btw. stability and
causality, see G.S.Danicol et al., J.Phys.G35,
115102(2008).
Entropy Current (2nd)
Assumption (2nd order theory):
Non-equilibrium entropy current vector has
linear + quadratic dissipative term(s)
constructed from (Vm, P, pmn, (um)).
Stability condition O.K.
The 2nd Law of Thermodynamics:
2nd order case
Sometimes
omitted,
but needed.
 Generalization of thermodynamic force!?
Same equation, but different definition of p and P.
Summary:
Constitutive Equations
w: vorticity
•Relaxation terms appear (tp and tP are relaxation
time).
•No longer algebraic equations! Dissipative currents
become dynamical quantities like thermodynamic
variables.
•Employed in recent viscous fluid simulations.
(Sometimes the last term is neglected.)
Plan of this Lecture
1st Day
• Hydrodynamics in Heavy Ion Collisions
• Collective flow
• Dynamical Modeling of heavy ion
collisions (seminar)
2nd Day
• Formalism of relativistic ideal/viscous
hydrodynamics
• Bjorken’s scaling solution with viscosity
• Effect of viscosity on particle spectra
(discussion)
PART 4
Bjorken’s Scaling Solution
with Viscosity
“Bjorken” Coordinate
t
Boost  parallel shift
Boost invariant
 Independent of hs
z
0
Bjorken’s Scaling Solution
Assuming boost invariance for thermodynamic
variables P=P(t) and 1D Hubble-like flow
Hydrodynamic equation for perfect fluids with a
simple EoS,
Conserved and Non-Conserved
Quantity in Scaling Solution
expansion
pdV work
Bjorken’s Equation
in the 1st Order Theory
(Bjorken’s solution)
= (1D Hubble flow)
Q. Derive the above equation.
Viscous Correction
Correction from shear viscosity
(in compressible fluids)
Correction from
bulk viscosity
 If these corrections vanish, the above equation
reduces to the famous Bjorken equation.
Expansion scalar = theta = 1/tau in scaling solution
Recent Topics on Transport
Coefficients
Need microscopic theory (e.g., Boltzmann eq.)
to obtain transport coefficients.
•
is obtained from
Super Yang-Mills theory.
Kovtun, Son, Starinet,…
•
is obtained from lattice.
Nakamura, Sakai,…
•Bulk viscosity has a prominent peak around Tc.
Kharzeev, Tuchin, Karsch, Meyer…
Bjorken’s Equation
in the 2nd Order Theory
where
New terms appear in the 2nd order theory.
 Coupled differential equations
Sometimes, the last terms are neglected.
Importance of these terms
 see Natsuume and Okamura, 0712.2917[hep-th].
Why only p00-pzz?
In EoM of energy density,
appears in spite of constitutive equations.
According to the Bjorken solution,
Relaxation Equation?
Digression:
Full 2nd order equation?
Beyond I-S equation, see R.Baier et al., JHEP
0804,100 (2008); Tsumura-Kunihiro?; D. Rischke,
talk at SQM 2008. According to Rischke’s talk,
constitutive equations with vanishing heat flow
are
Digression (contd.): Bjorken’s Equation
in the “full” 2nd order theory
See also, R.Fries et al.,PRC78,034913(2008).
Note that the equation for shear is valid only for conformal EOS and
that no 2nd and 3er terms for bulk.
Model EoS (crossover)
Crossover EoS:
Tc = 0.17GeV
D = Tc/50
dH = 3, dQ = 37
Energy-Momentum Tensor
at t0 in Comoving Frame
In what follows, bulk viscosity is neglected.
Numerical Results
(Temperature)
T0 = 0.22 GeV
t0 = 1 fm/c
h/s = 1/4p
tp = 3h/4p
Same initial
condition
(Energy momentum
tensor is isotropic)
Numerical code (C++) is available upon request.
Numerical Results
(Temperature)
T0 = 0.22 GeV
t0 = 1 fm/c
h/s = 1/4p
tp = 3h/4p
Same initial
condition
(Energy momentum
tensor is anisotropic)
Numerical code (C++) is available upon request.
Numerical Results
(Temperature)
Numerical code (C++) is available upon request.
Numerical Results (Entropy)
T0 = 0.22 GeV
t0 = 1 fm/c
h/s = 1/4p
tp = 3h/4p
Same initial
condition
(Energy momentum
tensor is isotropic)
Numerical code (C++) is available upon request.
Numerical Results (Entropy)
T0 = 0.22 GeV
t0 = 1 fm/c
h/s = 1/4p
tp = 3h/4p
Same initial
condition
(Energy momentum
tensor is anisotropic)
Numerical code (C++) is available upon request.
Numerical Results (Entropy)
Numerical code (C++) is available upon request.
Numerical Results
(Shear Viscosity)
Numerical code (C++) is available upon request.
Numerical Results
(Initial Condition Dependence
in the 2nd order theory)
Numerical code (C++) is available upon request.
Numerical Results
(Relaxation Time dependence)
Saturated values
non-trivial
Relaxation time larger
Maximum p is smaller
Relaxation time smaller
Suddenly relaxes to
1st order theory
Remarks
• Sometimes results from ideal hydro are
compared with the ones from 1st order
theory. But initial conditions must be
different.
• Be careful what is attributed for the
difference between two results.
• Sensitive to initial conditions and new
parameters (relaxation time for stress
tensor)
Plan of this Lecture
1st Day
• Hydrodynamics in Heavy Ion Collisions
• Collective flow
• Dynamical Modeling of heavy ion
collisions (seminar)
2nd Day
• Formalism of relativistic ideal/viscous
hydrodynamics
• Bjorken’s scaling solution with viscosity
• Effect of viscosity on particle spectra
(discussion)
PART 5
Effect of Viscosity on
Particle Spectra
Particle Spectra
in Hydrodynamic Model
• How to compare with experimental data
(particle spectra)?
• Free particles (l/L>>1) eventually
stream to detectors.
• Need prescription to convert
hydrodynamic (thermodynamic) fields
(l/L<<1) into particle picture.
• Need kinetic (or microscopic)
interpretation of hydrodynamic behavior.
Microscopic Interpretation
Single particle phase space density in local
thermal equilibrium:
Kinetic definition of current and energy momentum
tensor are
Matter in (Kinetic) Equilibrium
Kinetically equilibrated
matter at rest
Kinetically equilibrated
matter at finite velocity
um
py
py
px
Isotropic distribution
px
Lorentz-boosted distribution
Cooper-Frye Formula
•No dynamics of evaporation.
•Just counting the net number of particles
(out-going particles) - (in-coming particles)
through hypersurface S
•Negative contribution can appear at some
space-like hyper surface elements.
1st Moment
um is normalized, so we can always choose amn such that
1st Moment (contd.)
Vanishing for n= i due to odd function in integrant.
Q. Go through all steps in the above derivation.
2nd Moment
where,
Deviation from Equilibrium
Distribution
Neglecting anti-particles,
1
4
9
Grad’s 14 moments
EoM for epsilons can be also obtained from BE.
Taylor Expansion around
Equilibrium Distribution
Taylor Expansion around
Equilibrium Distribution (contd.)
14 Conditions
“Landau” conditions (2)
Viscosities (12)
Epsilons can be
expressed by
dissipative
currents
Relation btw. Coefficients and
Dissipative Currents
Finally,
For details, see Sec.6 and Appendix C in
Israel-Stewart paper. (Ann.Phys.118,341(1979))
Remarks
• Boltzmann equation gives a microscopic
interpretation of hydrodynamics.
• However, hydrodynamics may be
applied for gas/liquid where the
Boltzmann equation can not be applied.
• Deviation from equilibrium can be
incorporated into phase space
distribution.
References
Far from Complete List…
• General
–
–
–
–
L.D.Landau, E.M.Lifshitz, Fluid Mechanics, Section 133-136
L.P.Csernai, Introduction to Relativistic Heavy Ion Collisions
D.H.Rischke, nucl-th/9809044.
J.-Y.Ollitrault, 0708.2433[nucl-th].
• Viscous hydro + transport coefficient
–
–
–
–
–
–
–
–
–
C.Eckart, Phys.Rev.15,919(1940).
M.Namiki, C.Iso,Prog.Theor.Phys.18,591(1957)
C.Iso, K.Mori,M.Namiki, Prog.Theor.Phys.22,403(1959)
I.Mueller, Z. Phys. 198, 329 (1967)
W.Israel, Ann.Phys.100,310(1976)
W.Israel. J.M.Stewart, Ann.Phys.118,341(1979)
A.Hosoya, K.Kajantie, Nucl.Phys.B250, 666(1985)
P.Danielewicz,M.Gyulassy, Phys.Rev.D31,53(1985)
I.Muller, Liv.Rev.Rel 1999-1.
References for
Astrophysics/Cosmology
• N.Andersson, G.L.Comer, Liv.Rev.Rel.,2007-1
• R.Maartens, astro-ph/9609119.
Disclaimer:
I’m not familiar with these kinds of review…