LOG#028. Rockets and relativity.

The second post in this special thread of 3 devoted to Neil Armstrong memory has to do with rocketry.

Firstly, for completion, we are going to study the motion of a rocket in “vacuum” according to classical physics. Then, we will deduce the relatistic rocket equation and its main properties.

CLASSICAL NON-RELATIVISTIC ROCKETS

The fundamental law of Dynamics, following Sir Isaac Newton, reads:

\mathbf{F}=\dfrac{d\mathbf{p}}{dt}

Suppose a rocket with initial mass M_i and initial velocity u_i=0. It ejects mass of propellant “gas” with “gas speed” (particles of gas have a relative velocity or speed with respect to the rest observer when the rocket move at speed \mathbf{v}) equals to u_0 (note that the relative speed will be u_{rel} and the propellant mass is latex m_0$. Generally, this speed is also called “exhaust velocity” by engineers. The motion of a variable mass or rocket is given by the so-called Metcherski’s equation:

\boxed{M\dfrac{d\mathbf{v}}{dt}=-\mathbf{u_0}\dfrac{dM}{dt}+\mathbf{F}}

where -\mathbf{u_0}=\mathbf{v_{gas}}-\mathbf{v}. The Metcherski’s equation can be derived as follows: the rocket changes its mass and velocity so M'=M+dM and V'=V+dV, so the change in momentum is equal to M'V'=(M+dM)(V+dV), plus an additional term v_{gas}dv_{gas} and -mV. Therefore, the total change in momentum:

dP=Fdt=(M+dM)(V+dV)+v_{gas}dv_{gas}-mV

Neglecting second order differentials, and setting the conservation of mass (we are in the non-relativistic case)

dM+dm_{gas}=0

we recover

MdV=v_{rel}dM+Fdt

that represents (with the care of sign in relative speed) the Metcherski equation we have written above.

Generally speaking, the “force” due to the change in “mass” is called thrust.  With no external force, from the remaining equation of the thrust and velocity, and it can be easily integrated

u_f=-u_0\int_{M_i}^{M_f}dM

and thus we get the Tsiolkowski’s rocket equation:

\boxed{\mathbf{u_f}=\mathbf{u_0}\ln \dfrac{M_i}{M_f}}

Engineers use to speak about the so-called mass ratio R=\dfrac{M_f}{M_i}, although sometimes the reciprocal definition is also used for such a ratio so be aware, and in terms of this the Tsiokolski’s equation reads:

\boxed{\mathbf{u_f}=\mathbf{u_0}\ln \dfrac{1}{R}}

We can invert this equation as well, in order to get

\boxed{R=\dfrac{M_f}{M_i}=\exp\left(-\dfrac{u_f}{u_0}\right)}

Example: Calculate the fraction of mass of a one-stage rocket to reach the Earth’s orbit. Typical values for u_f=8km/s and u_0=4km/s show that the mass ratio is equal to R=0.14. Then, only the 14\% of the initial mass reaches the orbit, and the remaining mass is fuel.

Multistate rockets offer a good example of how engineer minds work. They have discovered that a multistage rocket is more effective than the one-stage rocket in terms of maximum attainable speed and mass ratios. The final n-stage lauch system for rocketry states that the final velocity is the sum of the different gains in the velocity after the n-th stage, so we can obtain

\displaystyle{u_f=\sum_{i=1}^{n}u_i^f=u_1^f+\cdots+u_n^f}

After the n-th step, the change in velocity reads

u_i^f=c_i\ln \dfrac{1}{R_i}

where the i-th mass ratios are definen recursively as the final mass in the n-th step and the initial mass in that step, so we have

\displaystyle{u_f=\sum_i c_i\ln \dfrac{1}{R_i}}

and we define the total mass ratio:

\displaystyle{R_T=\prod_i R_i}

If the average effective rocket exhaust velocity is the same in every step/stage, e.g. c_i=c, we get

\displaystyle{u_f=c\ln \left( \prod_{i=1}^{n} R_i^{-1}\right)}

or

\displaystyle{u_f=c \ln \left[ \left(\dfrac{M_0}{M_f}\right)_1\left(\dfrac{M_0}{M_f}\right)_2\cdots \left(\dfrac{M_0}{M_f}\right)_n\right]=c\ln \left[\left(\dfrac{M_0}{M_f}\right)_T\right]}

The influence of the number of steps, for a given exhaust velocity, in the final attainable velocity can be observed in the next plots:

 RELATIVISTIC ROCKETS

We proceed now to the relativistic generalization of the previous rocketry. An observer in the laboratory frame observes that total momentum is conserved, of course, and so:

M'du'=-u'_0dM'

where du' is the velocity increase in the rocket with a rest mass M’ in the instantaneous reference frame of the moving rocket S’. It is NOT equal to its velocity increase measured in the unprimed reference frame, du. Due to the addition theorem of velocities in SR, we have

u+du=\dfrac{u+du'}{1+\dfrac{udu'}{c^2}}

where u is the instantenous velocity of the rocket with respect to the laboratory frame S. We can perform a Taylor expansion of the denominator in the last equation, in order to obtain:

u+du=(u+du')\left(1-\dfrac{udu'}{c^2}\right)

and then

u+du=u+du'\left(1-\dfrac{u^2}{c^2}\right)

and finally, we get

du'=\dfrac{du}{1-\dfrac{u^2}{c^2}}=\gamma^2_u du

Plugging this equation into the above equation for mass (momentum), and integrating

\displaystyle{\int_{0}^{u_f}\dfrac{du}{1-\dfrac{u^2}{c^2}}=-u'_0\int_{M'_0}^{M'_f}dM}

we deduce that the relativistic version of the Tsiolkovski’s rocket equation, the so-called relativistic rocket equation, can be written as:

\dfrac{c}{2}\ln \dfrac{1+\dfrac{u_f}{c}}{1-\dfrac{u_f}{c}}=u'_0\ln\dfrac{M'_i}{M'_f}

We can suppress the primes if we remember that every data is in the S’-frame (instantaneously), and rewrite the whole equation in the more familiar way:

\boxed{u_f=c\dfrac{1-\left(\dfrac{M_f}{M_0}\right)^{\frac{2u_0}{c}}}{1+\left(\dfrac{M_f}{M_0}\right)^{\frac{2u_0}{c}}}=c\dfrac{1-R^{\frac{2u_0}{c}}}{1+R^{\frac{2u_0}{c}}}}

where the mass ratio is defined as before R=\dfrac{M_f}{M_i}. Now, comparing the above equation with the rapidity/maximum velocity in the uniformly accelerated motion:

u_f=c\tanh \left(\dfrac{g\tau}{c}\right)

we get that relativistic rocket equation can be also written in the next manner:

u_f=c\tanh \left[ -\dfrac{u_0}{c}\ln \left(\dfrac{1}{R}\right)\right]

or equivalently

u_f=c\tanh \left[ \dfrac{u_0}{c}\ln R\right]

since we have in this case

\dfrac{g\tau}{c}=-\dfrac{u_0}{c}\ln \left(\dfrac{1}{R}\right)=\dfrac{u_0}{c}\ln R

and thus

R^{\frac{u_0}{c}}=\left(\dfrac{M_f}{M_i}\right)^{\frac{u_0}{c}}=\exp \left(-\dfrac{g\tau}{c}\right)

If the propellant particles move at speed of light, e.g., they are “photons” or ultra-relativistic particles that move close to the speed of light we have the celebrated “photon rocket”. In that case, setting u_0=c, we would obtain that:

\boxed{u_f=c\dfrac{1-\left(\dfrac{M_f}{M_0}\right)^{2}}{1+\left(\dfrac{M_f}{M_0}\right)^{2}}=c\dfrac{1-R^{2}}{1+R^{2}}=c\tanh \ln R}

and where for the photon rocket (or the ultra-relativistic rocket) we have as well

\dfrac{g\tau}{c}=-\ln \left(\dfrac{1}{R}\right)=\ln R

Final remark: Instead of the mass ratio, sometimes is more useful to study the ratio fuel mass/payload. In that case, we set M_f=m and M_0=m+M, where M is the fuel mass and m is the payload. So, we would write

R=\dfrac{m}{m+M}

so then the ratio fuel mass/payload will be

\dfrac{M}{m}=R^{-1}-1=\exp \left(\dfrac{g\tau}{c}\right)-1

We are ready to study the interstellar trip with our current knowledge of Special Relativity and Rocketry. We will study the problem in the next and final post of this fascinating thread. Stay tuned!

Advertisements


Leave a Reply

Please log in using one of these methods to post your comment:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s