## 12 The Newtonian Limit and the Euler Equations

One reason relativistic fluids are needed is that they can be used to model neutron stars. However, even though neutron stars are clearly general relativistic objects, one often starts with good old Newtonian physics when considering new applications. The main reason is that effects (such as acoustic modes of oscillation) which are primarily due to the fluids that make up the star can often be understood qualitatively from Newtonian calculations. There are also certain regimes in a neutron star where the Newtonian limit is sufficient quantitatively (such as the outer layers).

There has been much progress recently in the analysis of Newtonian multiple fluid systems. Prix [91] has developed an action-based formalism, analogous to what has been used here. Carter and Chamel [222324] have done the same except that they use a fully spacetime covariant formalism. We will be somewhat less ambitious (for example, as in [5]) by extracting the Newtonian equations as the non-relativistic limit of the fully relativistic equations. Given the results of Prix as well as of Carter and Chamel we can think of this exercise as a consistency check of our equations.

The Newtonian limit consists of writing the general relativistic field equations to order where is the speed of light. The Newtonian equations are obtained, strictly speaking, in the limit where . To order the metric becomes

where the are Cartesian coordinates and is the metric of Equation (22). The gravitational potential, denoted , is assumed to be small in the sense that
We will be somewhat lax by simply writing the general relativistic equations to , and thereby avoiding some of the additional complexities that occur because the spacetime metric clearly becomes singular as . For the reader interested in a rigorous construction of the Newtonian equations, we recommend the series of papers by Carter and Chamel [222324].

If is the proper time as measured along a fluid element worldline, then the curve it traces out can be written

Recall that its four-velocity is given by (cf. Equation (6))
but because the speed of light has been restored, . The proper time is now given by
where is the Newtonian three-velocity of the fluid. It is considered to be small in the sense that
Hence, to the correct order the four-velocity components are
where . Note that the particle number current becomes
The factor of appears because we define , regardless of whether or the meter/second value.

In order to write the single fluid Euler equations, keeping terms to the required order, it is necessary to explicitly break up the master function into its mass, kinetic, and “potential” energy parts, i.e. to write as

where is the particle mass. The internal energy is small compared to the mass, meaning
With this choice, a variation of that leaves the metric fixed yields
where
The coefficient reduces to
In terms of these variables, the pressure is seen to be
The Newtonian limit of the general relativistic fluid equations, in a general coordinate basis, is then
and
where is the Euclidean space covariant derivative. Of course, the gravitational potential is determined by a Poisson equation with as source.

A Newtonian two-fluid system can be obtained in a similar fashion. As discussed in Section 10, the main difference is that we need two sets of worldlines, describable, say, by curves where is the proper time along a constituent’s worldline. Of course, entrainment also comes into play. Its presence implies that the relative flow of the fluids is required to specify the local thermodynamic state of the system, and that the momentum of a given fluid is not simply proportional to that fluid’s flux. This is the situation for superfluid  [94110], where the entropy can flow independently of the superfluid Helium atoms. Superfluid can also be included in the mixture, in which case there will be a relative flow of the isotope with respect to , and relative flows of each with respect to the entropy [113].

Let us consider a two-fluid model like a mixture of and , or neutrons and protons in a neutron star. We will denote the two fluids as fluids and . The magnitude squared of the difference of three-velocities

is denoted so that the equation of state takes the form . Hence,
where
The coefficient reflects the effect of entrainment on the equation of state. Similarly, entrainment causes the fluid momenta to be modified to

The number density of each fluid obeys a continuity equation:

Each fluid is also seen to satisfy an Euler-type equation, which ensures the conservation of total momentum. This equation can be written
where
i.e.  is the relevant chemical potential per unit mass, and the entrainment is included via the coefficients
For a detailed discussion of these equations, see [916].