I am reading this paper by Thorne, Flammang & Zytkow (1981), which discusses the dynamics of spherical accretion onto a black hole in the Schwarzschild metric ($c=G=1$ units):
$$ds^2=-(1-2M/r)dt^2+(1-2M/r)^{-1}dr^2+r^2(d\theta^2+\sin^2\theta d\phi^2)$$
For a gas with inward velocity $v$ (as measured by an observer at rest in this metric), equations 2-3 of the paper give this expression for the four-velocity: $$\mathbf{u}=y(1-2M/r)^{-1}\frac{\partial}{\partial t}-vy\frac{\partial}{\partial r}$$ $$y\equiv\mathbf{u}\cdot\frac{\partial}{\partial t}=(1-2M/r)^{1/2}(1-v^2)^{-1/2}$$ $y$ is the so-called energy parameter, and it would be a constant if the trajectory was a geodesic.
I am having trouble understanding how to obtain this expression. Here is what I tried so far:
- The velocity is $v=-\frac{dr}{dt}$ (minus sign because the gas is moving inward)
- The four-velocity components are usually written as $U^{\mu}=\frac{dx^\mu}{d\tau}$, where $\tau$ is the proper time
- By the chain rule, $\frac{dr}{d\tau}=\frac{dr}{dt}\frac{dt}{d\tau}=-v\frac{dt}{d\tau}$
- The four-velocity as a vector is $\mathbf{u}=U^\mu\partial_\mu=U^\mu\frac{\partial}{\partial x^\mu}$
So that: $$\mathbf{u}=\frac{dt}{d\tau}\frac{\partial}{\partial t}+\frac{dr}{d\tau}\frac{\partial}{\partial r} \quad\\=\frac{dt}{d\tau}\left(\frac{\partial}{\partial t}-v\frac{\partial}{\partial r}\right)$$
Further, since $ds^2=-d\tau^2$, we have: $$-1=-(1-2M/r)\left(\frac{dt}{d\tau}\right)^2+(1-2M/r)^{-1}\left(\frac{dr}{d\tau}\right)^2\\ =\left(\frac{dt}{d\tau}\right)^2\left(-(1-2M/r)+(1-2M/r)^{-1}v^2\right) $$
I don't see how this can lead to the expression for the four-velocity from the paper, so I must be making a mistake somewhere. Any help would be appreciated.