Relativistic Euler equations

This article needs additional citations for verification. Please help improve this article by adding citations to reliable sources. Unsourced material may be challenged and removed. Find sources: "Relativistic Euler equations" – news · newspapers · books · scholar · JSTOR (April 2020) (Learn how and when to remove this message)

In

${\displaystyle c=1}$ and the

Motivation

For most fluids observable on Earth, traditional fluid mechanics based on Newtonian mechanics is sufficient. However, as the fluid velocity approaches the speed of light or moves through strong gravitational fields, or the pressure approaches the energy density (${\displaystyle P\sim \rho }$), these equations are no longer valid.^[2] Such situations occur frequently in astrophysical applications. For example, gamma-ray bursts often feature speeds only ${\displaystyle 0.01\%}$ less than the speed of light,^[3] and neutron stars feature gravitational fields that are more than ${\displaystyle 10^{11}}$ times stronger than the Earth's.^[4] Under these extreme circumstances, only a relativistic treatment of fluids will suffice.

Introduction

The equations of motion are contained in the continuity equation of the stress–energy tensor ${\displaystyle T^{\mu \nu }}$:

${\displaystyle \nabla _{\mu }T^{\mu \nu }=0,}$

where ${\displaystyle \nabla _{\mu }}$ is the covariant derivative.^[5] For a perfect fluid,

${\displaystyle T^{\mu \nu }\,=(e+p)u^{\mu }u^{\nu }+pg^{\mu \nu }.}$

Here ${\displaystyle e}$ is the total mass-energy density (including both rest mass and internal energy density) of the fluid, ${\displaystyle p}$ is the

${\displaystyle u^{\mu }}$

${\displaystyle g^{\mu \nu }}$

. If ${\displaystyle n}$ is the

${\displaystyle \nabla _{\mu }(nu^{\mu })=0.}$

These equations reduce to the classical Euler equations if the fluid three-velocity is much less than the speed of light, the pressure is much less than the energy density, and the latter is dominated by the rest mass density. To close this system, an equation of state, such as an ideal gas or a Fermi gas, is also added.^[1]

Equations of Motion in Flat Space

In the case of flat space, that is ${\displaystyle \nabla _{\mu }=\partial _{\mu }}$ and using a metric signature of ${\displaystyle (-,+,+,+)}$, the equations of motion are,^[6]

${\displaystyle (e+p)u^{\mu }\partial _{\mu }u^{\nu }=-\partial ^{\nu }p-u^{\nu }u^{\mu }\partial _{\mu }p}$

Where ${\displaystyle e=\gamma \rho c^{2}+\rho \epsilon }$ is the energy density of the system, with ${\displaystyle p}$ being the pressure, and ${\displaystyle u^{\mu }=\gamma (1,{\frac {\mathbf {v} }{c}})}$ being the four-velocity of the system.

Expanding out the sums and equations, we have, (using ${\displaystyle {\frac {d}{dt}}}$ as the material derivative)

${\displaystyle (e+p){\frac {\gamma }{c}}{\frac {du^{\mu }}{dt}}=-\partial ^{\mu }p-{\frac {\gamma }{c}}{\frac {dp}{dt}}u^{\mu }}$

Then, picking ${\displaystyle u^{\nu }=u^{i}={\frac {\gamma }{c}}v_{i}}$ to observe the behavior of the velocity itself, we see that the equations of motion become

${\displaystyle (e+p){\frac {\gamma }{c^{2}}}{\frac {d}{dt}}(\gamma v_{i})=-\partial _{i}p-{\frac {\gamma ^{2}}{c^{2}}}{\frac {dp}{dt}}v_{i}}$

Note that taking the non-relativistic limit, we have ${\displaystyle {\frac {1}{c^{2}}}(e+p)=\gamma \rho +{\frac {1}{c^{2}}}\rho \epsilon +{\frac {1}{c^{2}}}p\approx \rho }$. This says that the energy of the fluid is dominated by its

In this limit, we have ${\displaystyle \gamma \rightarrow 1}$ and ${\displaystyle c\rightarrow \infty }$, and can see that we return the Euler Equation of ${\displaystyle \rho {\frac {dv_{i}}{dt}}=-\partial _{i}p}$.

Derivation of the Equations of Motion

In order to determine the equations of motion, we take advantage of the following spatial projection tensor condition:

${\displaystyle \partial _{\mu }T^{\mu \nu }+u_{\alpha }u^{\nu }\partial _{\mu }T^{\mu \alpha }=0}$

We prove this by looking at ${\displaystyle \partial _{\mu }T^{\mu \nu }+u_{\alpha }u^{\nu }\partial _{\mu }T^{\mu \alpha }}$ and then multiplying each side by ${\displaystyle u_{\nu }}$. Upon doing this, and noting that ${\displaystyle u^{\mu }u_{\mu }=-1}$, we have ${\displaystyle u_{\nu }\partial _{\mu }T^{\mu \nu }-u_{\alpha }\partial _{\mu }T^{\mu \alpha }}$. Relabeling the indices ${\displaystyle \alpha }$ as ${\displaystyle \nu }$ shows that the two completely cancel. This cancellation is the expected result of contracting a temporal tensor with a spatial tensor.

Now, when we note that

${\displaystyle T^{\mu \nu }=wu^{\mu }u^{\nu }+pg^{\mu \nu }}$

where we have implicitly defined that ${\displaystyle w\equiv e+p}$, we can calculate that

${\displaystyle {\begin{aligned}\partial _{\mu }T^{\mu \nu }&=(\partial _{\mu }w)u^{\mu }u^{\nu }+w(\partial _{\mu }u^{\mu })u^{\nu }+wu^{\mu }\partial _{\mu }u^{\nu }+\partial ^{\nu }p\\\partial _{\mu }T^{\mu \alpha }&=(\partial _{\mu }w)u^{\mu }u^{\alpha }+w(\partial _{\mu }u^{\mu })u^{\alpha }+wu^{\mu }\partial _{\mu }u^{\alpha }+\partial ^{\alpha }p\end{aligned}}}$

and thus

${\displaystyle u^{\nu }u_{\alpha }\partial _{\mu }T^{\mu \alpha }=(\partial _{\mu }w)u^{\mu }u^{\nu }u^{\alpha }u_{\alpha }+w(\partial _{\mu }u^{\mu })u^{\nu }u^{\alpha }u_{\alpha }+wu^{\mu }u^{\nu }u_{\alpha }\partial _{\mu }u^{\alpha }+u^{\nu }u_{\alpha }\partial ^{\alpha }p}$

Then, let's note the fact that ${\displaystyle u^{\alpha }u_{\alpha }=-1}$ and ${\displaystyle u^{\alpha }\partial _{\nu }u_{\alpha }=0}$. Note that the second identity follows from the first. Under these simplifications, we find that

${\displaystyle u^{\nu }u_{\alpha }\partial _{\mu }T^{\mu \alpha }=-(\partial _{\mu }w)u^{\mu }u^{\nu }-w(\partial _{\mu }u^{\mu })u^{\nu }+u^{\nu }u^{\alpha }\partial _{\alpha }p}$

and thus by ${\displaystyle \partial _{\mu }T^{\mu \nu }+u_{\alpha }u^{\nu }\partial _{\mu }T^{\mu \alpha }=0}$, we have

${\displaystyle (\partial _{\mu }w)u^{\mu }u^{\nu }+w(\partial _{\mu }u^{\mu })u^{\nu }+wu^{\mu }\partial _{\mu }u^{\nu }+\partial ^{\nu }p-(\partial _{\mu }w)u^{\mu }u^{\nu }-w(\partial _{\mu }u^{\mu })u^{\nu }+u^{\nu }u^{\alpha }\partial _{\alpha }p=0}$

We have two cancellations, and are thus left with

${\displaystyle (e+p)u^{\mu }\partial _{\mu }u^{\nu }=-\partial ^{\nu }p-u^{\nu }u^{\alpha }\partial _{\alpha }p}$

See also

• Relativistic heat conduction
• Equation of state (cosmology)

References

This relativity-related article is a stub. You can help Wikipedia by expanding it.

• v
• t
• e