Most popular derivation Consider the following system: In the following derivation, "the rocket" is taken to mean "the rocket and all of its unexpended propellant".
Newton's second law of motion relates external forces (\vec{F}_i) to the change in linear momentum of the whole system (including rocket and exhaust) as follows: \sum_i \vec{F}_i = \lim_{\Delta t \to 0} \frac{\vec{P}_{\Delta t} - \vec{P}_0}{\Delta t} where \vec{P}_0 is the momentum of the rocket at time t = 0: \vec{P}_0 = m \vec{V} and \vec{P}_{\Delta t} is the momentum of the rocket and exhausted mass at time t = \Delta t: \vec{P}_{\Delta t} = \left(m - \Delta m \right) \left(\vec{V} + \Delta \vec{V} \right) + \Delta m \vec{V}_\text{e} and where, with respect to the observer: • \vec{V} is the velocity of the rocket at time t = 0 • \vec{V} + \Delta \vec{V} is the velocity of the rocket at time t = \Delta t • \vec{V}_\text{e} is the velocity of the mass added to the exhaust (and lost by the rocket) during time \Delta t • m is the mass of the rocket at time t = 0 • \left( m - \Delta m \right) is the mass of the rocket at time t = \Delta t The velocity of the exhaust \vec{V}_\text{e} in the observer frame is related to the velocity of the exhaust in the rocket frame v_\text{e} by: \vec {v}_\text{e} = \vec{V}_\text{e} - \vec{V} thus, \vec {V}_\text{e} = \vec{V} + \vec{v}_\text{e} Solving this yields: \vec{P}_{\Delta t} - \vec{P}_0 = m\Delta \vec{V} + \vec{v}_\text{e} \Delta m - \Delta m \Delta \vec{V} If \vec{V} and \vec{v}_\text{e} are opposite, \vec{F}_\text{i} have the same direction as \vec{V}, \Delta m \Delta \vec{V} are negligible (since dm \, d\vec{v} \to 0), and using dm = -\Delta m (since ejecting a positive \Delta m results in a decrease in rocket mass in time), \sum_i F_i = m \frac{dV}{dt} + v_\text{e} \frac{dm}{dt} If there are no external forces then \sum_i F_i = 0 (
conservation of linear momentum) and -m\frac{dV}{dt} = v_\text{e}\frac{dm}{dt} Assuming that v_\text{e} is constant (known as
Tsiolkovsky's hypothesis \Delta v = v_\text{eff} \sum ^{j=N}_{j=1} \frac{\phi/N}{\sqrt{(1-j\phi/N)(1-j\phi/N+\phi/N)}} Notice that for large N the last term in the denominator \phi/N\ll 1 and can be neglected to give \Delta v \approx v_\text{eff} \sum^{j=N}_{j=1}\frac{\phi/N}{1-j\phi/N} = v_\text{eff} \sum ^{j=N}_{j=1} \frac{\Delta x}{1-x_j} where \Delta x = \frac{\phi}{N} and x_j = \frac{j\phi}{N} . As N\rightarrow \infty this
Riemann sum becomes the definite integral \lim_{N\to\infty}\Delta v = v_\text{eff} \int_{0}^{\phi} \frac{dx}{1-x} = v_\text{eff}\ln \frac{1}{1-\phi} = v_\text{eff} \ln \frac{m_0}{m_f} , since the final remaining mass of the rocket is m_f = m_0(1-\phi).
Special relativity If
special relativity is taken into account, the following equation can be derived for a
relativistic rocket, with \Delta v again standing for the rocket's final velocity (after expelling all its reaction mass and being reduced to a rest mass of m_1) in the
inertial frame of reference where the rocket started at rest (with the rest mass including fuel being m_0 initially), and c standing for the
speed of light in vacuum: \frac{m_0}{m_1} = \left[\frac{1 + {\frac{\Delta v}{c}}}{1 - {\frac{\Delta v}{c}}}\right]^{\frac{c}{2v_\text{e}}} Writing \frac{m_0}{m_1} as R allows this equation to be rearranged as \frac{\Delta v}{c} = \frac{R^{\frac{2v_\text{e}}{c}} - 1}{R^{\frac{2v_\text{e}}{c}} + 1} Then, using the
identity R^{\frac{2v_\text{e}}{c}} = \exp \left[ \frac{2v_\text{e}}{c} \ln R \right] (here "exp" denotes the
exponential function;
see also Natural logarithm as well as the "power" identity at
logarithmic identities) and the identity \tanh x = \frac{e^{2x} - 1} {e^{2x} + 1} (
see Hyperbolic function), this is equivalent to \Delta v = c \tanh\left(\frac {v_\text{e}}{c} \ln \frac{m_0}{m_1} \right) ==Terms of the equation==