There are different derivations for the variable-mass system motion equation, depending on whether the mass is entering or leaving a body (in other words, whether the moving body's mass is increasing or decreasing, respectively). To simplify calculations, all bodies are considered as particles. It is also assumed that the mass is unable to apply external forces on the body outside of accretion/ablation events.

Mass accretion edit

The following derivation is for a body that is gaining mass (accretion). A body of time-varying mass m moves at a velocity v at an initial time t. In the same instant, a particle of mass dm moves with velocity u with respect to ground. The initial momentum can be written as^[5]

${\displaystyle \mathbf {p} _{\mathrm {1} }=m\mathbf {v} +\mathbf {u} \mathrm {d} m}$ 

Now at a time t + dt, let both the main body and the particle accrete into a body of velocity v + dv. Thus the new momentum of the system can be written as

${\displaystyle \mathbf {p} _{\mathrm {2} }=(m+\mathrm {d} m)(\mathbf {v} +\mathrm {d} \mathbf {v} )=m\mathbf {v} +m\mathrm {d} \mathbf {v} +\mathbf {v} \mathrm {d} m+\mathrm {d} m\mathrm {d} \mathbf {v} }$ 

Since dmdv is the product of two small values, it can be ignored, meaning during dt the momentum of the system varies for

${\displaystyle \mathrm {d} \mathbf {p} =\mathbf {p} _{\mathrm {2} }-\mathbf {p} _{\mathrm {1} }=(m\mathbf {v} +m\mathrm {d} \mathbf {v} +\mathbf {v} \mathrm {d} m)-(m\mathbf {v} +\mathbf {u} \mathrm {d} m)=m\mathrm {d} \mathbf {v} -(\mathbf {u} -\mathbf {v} )\mathrm {d} m}$ 

Therefore, by Newton's second law

${\displaystyle \mathbf {F} _{\mathrm {ext} }={\frac {\mathrm {d} \mathbf {p} }{\mathrm {d} t}}={\frac {m\mathrm {d} \mathbf {v} -(\mathbf {u} -\mathbf {v} )\mathrm {d} m}{\mathrm {d} t}}=m{\frac {\mathrm {d} \mathbf {v} }{\mathrm {d} t}}-(\mathbf {u} -\mathbf {v} ){\frac {\mathrm {d} m}{\mathrm {d} t}}}$ 

Noting that u - v is the velocity of dm relative to m, symbolized as v_rel, this final equation can be arranged as^[6]

${\displaystyle \mathbf {F} _{\mathrm {ext} }+\mathbf {v} _{\mathrm {rel} }{\frac {\mathrm {d} m}{\mathrm {d} t}}=m{\mathrm {d} \mathbf {v} \over \mathrm {d} t}}$ 

Mass ablation/ejection edit

In a system where mass is being ejected or ablated from a main body, the derivation is slightly different. At time t, let a mass m travel at a velocity v, meaning the initial momentum of the system is

${\displaystyle \mathbf {p} _{\mathrm {1} }=m\mathbf {v} }$ 

Assuming u to be the velocity of the ablated mass dm with respect to the ground, at a time t + dt the momentum of the system becomes

${\displaystyle \mathbf {p} _{\mathrm {2} }=(m-\mathrm {d} m)(\mathbf {v} +\mathrm {d} \mathbf {v} )+\mathbf {u} \mathrm {d} m=m\mathbf {v} +m\mathrm {d} \mathbf {v} -\mathbf {v} \mathrm {d} m-\mathrm {d} m\mathrm {d} \mathbf {v} +\mathbf {u} \mathrm {d} m}$ 

where u is the velocity of the ejected mass with respect to ground, and is negative because the ablated mass moves in opposite direction to the mass. Thus during dt the momentum of the system varies for

${\displaystyle \mathrm {d} \mathbf {p} =\mathbf {p} _{\mathrm {2} }-\mathbf {p} _{\mathrm {1} }=(m\mathbf {v} +m\mathrm {d} \mathbf {v} -\mathrm {d} \mathbf {m} \mathrm {d} \mathbf {v} -\mathbf {v} \mathrm {d} m+\mathbf {u} \mathrm {d} m)-(m\mathbf {v} )=m\mathrm {d} \mathbf {v} +[\mathbf {u} -(\mathbf {v} +\mathrm {d} \mathbf {v} )]\mathrm {d} m}$ 

Relative velocity v_rel of the ablated mass with respect to the mass m is written as

${\displaystyle \mathbf {v} _{\mathrm {rel} }=\mathbf {u} -(\mathbf {v} +\mathrm {d} \mathbf {v} )}$ 

Therefore, change in momentum can be written as

${\displaystyle \mathrm {d} \mathbf {p} =m\mathrm {d} \mathbf {v} +\mathbf {v} _{\mathrm {rel} }\mathrm {d} m}$ 

Therefore, by Newton's second law

${\displaystyle \mathbf {F} _{\mathrm {ext} }={\frac {\mathrm {d} \mathbf {p} }{\mathrm {d} t}}={\frac {m\mathrm {d} \mathbf {v} +\mathbf {v} _{\mathrm {rel} }\mathrm {d} m}{\mathrm {d} t}}=m{\frac {\mathrm {d} \mathbf {v} }{\mathrm {d} t}}+\mathbf {v} _{\mathrm {rel} }{\frac {\mathrm {d} m}{\mathrm {d} t}}}$ 

Therefore, the final equation can be arranged as

${\displaystyle \mathbf {F} _{\mathrm {ext} }-\mathbf {v} _{\mathrm {rel} }{\frac {\mathrm {d} m}{\mathrm {d} t}}=m{\mathrm {d} \mathbf {v} \over \mathrm {d} t}}$ 

By the definition of acceleration, a = dv/dt, so the variable-mass system motion equation can be written as

${\displaystyle \mathbf {F} _{\mathrm {ext} }+\mathbf {v} _{\mathrm {rel} }{\frac {\mathrm {d} m}{\mathrm {d} t}}=m\mathbf {a} }$ 

In bodies that are not treated as particles a must be replaced by a_cm, the acceleration of the center of mass of the system, meaning

${\displaystyle \mathbf {F} _{\mathrm {ext} }+\mathbf {v} _{\mathrm {rel} }{\frac {\mathrm {d} m}{\mathrm {d} t}}=m\mathbf {a} _{\mathrm {cm} }}$ 

Often the force due to thrust is defined as ${\displaystyle \mathbf {F} _{\mathrm {thrust} }=\mathbf {v} _{\mathrm {rel} }{\frac {\mathrm {d} m}{\mathrm {d} t}}}$  so that

${\displaystyle \mathbf {F} _{\mathrm {ext} }+\mathbf {F} _{\mathrm {thrust} }=m\mathbf {a} _{\mathrm {cm} }}$ 

This form shows that a body can have acceleration due to thrust even if no external forces act on it (F_ext = 0). Note finally that if one lets F_net be the sum of F_ext and F_thrust then the equation regains the usual form of Newton's second law:

${\displaystyle \mathbf {F} _{\mathrm {net} }=m\mathbf {a} _{\mathrm {cm} }}$ 

Ideal rocket equation edit

The ideal rocket equation, or the Tsiolkovsky rocket equation, can be used to study the motion of vehicles that behave like a rocket (where a body accelerates itself by ejecting part of its mass, a propellant, with high speed). It can be derived from the general equation of motion for variable-mass systems as follows: when no external forces act on a body (F_ext = 0) the variable-mass system motion equation reduces to^[2]

${\displaystyle \mathbf {v} _{\mathrm {rel} }{\frac {\mathrm {d} m}{\mathrm {d} t}}=m{\frac {\mathrm {d} \mathbf {v} }{\mathrm {d} t}}}$ 

If the velocity of the ejected propellant, v_rel, is assumed have the opposite direction as the rocket's acceleration, dv/dt, the scalar equivalent of this equation can be written as

${\displaystyle -v_{\mathrm {rel} }{\frac {\mathrm {d} m}{\mathrm {d} t}}=m{\mathrm {d} v \over \mathrm {d} t}}$ 

from which dt can be canceled out to give

${\displaystyle -v_{\mathrm {rel} }\mathrm {d} m=m\mathrm {d} v\,}$ 

Integration by separation of variables gives

${\displaystyle -v_{\mathrm {rel} }\int _{m_{0}}^{m_{1}}{\frac {\mathrm {d} m}{m}}=\int _{v_{0}}^{v_{1}}\mathrm {d} v}$ 

${\displaystyle v_{\mathrm {rel} }\ln {\frac {m_{0}}{m_{1}}}=v_{1}-v_{0}}$ 

By rearranging and letting Δv = v₁ - v₀, one arrives at the standard form of the ideal rocket equation:

${\displaystyle \Delta v=v_{\mathrm {rel} }\ln {\frac {m_{0}}{m_{1}}}}$ 

where m₀ is the initial total mass, including propellant, m₁ is the final total mass, v_rel is the effective exhaust velocity (often denoted as v_e), and Δv is the maximum change of speed of the vehicle (when no external forces are acting).