Momentum, Impulse and Collisions

Metal spheres A and B colliding with an explosion and flying sparks

  • Young et al., 20201
  • Image generated by WordPress AI
  • Graph created using Desmos

Table of Contents:

  1. Momentum and Impulse (Young et al., 2020)
    1. Newton’s Law In Terms Of Momentum
    2. The Impulse-Momentum Theorem
    3. Momentum and Kinetic Energy Compared
  2. Conservation of Momentum (Young et al., 2020)
    1. Types of Forces
  3. Momentum Conservation and Collisions (Young et al., 2020)
    1. Elastic and Inelastic Collisions
      1. Rules to Remember for Collisions
    2. Completely Inelastic Collisions
  4. Elastic Collisions (Young et al., 2020)
    1. Elastic Collisions, One Object is at Rest
    2. Elastic Collisions and Relative Velocity
  5. Center of Mass (Young et al., 2020)
    1. Center of Mass of Solids
    2. Motion of the Center of Mass
    3. External Forces and Center-of-Mass Motion
  6. Rocket Propulsion (Young et al., 2020)

Momentum and Impulse (Young et al., 2020)

Newton’s Law In Terms Of Momentum

Consider a particle with constant mass m. Newton’s second law of motion can then be written as;

ΣF=ma=mdvdt=ddt(mv)\Sigma \vec{F} = m\vec{a} = m\frac{d\vec{v}}{dt} = \frac{d}{dt}(m\vec{v})…..(1)

Note: Since the given mass m is constant, it can be moved inside the derivative.

Thus, the net force acting on the particle is equal to the rate of change of, product of mass (m) and velocity vector (v)(\vec{v}), with respect to time.

The quantity in the parenthesis in equation (1) can then be defined as momentum or linear momentum of the particle;

p=mv\vec{p} = m\vec{v}…..(2)

The magnitude of momentum is directly proportional to the particle’s mass and speed. Larger is the mass and speed of the particle, larger is it’s momentum.

Momentum is a vector quantity. It has the same direction as the particle’s velocity. It’s SI unit is kg ˙\dot{} m/s and it’s plural form is momenta.

Figure 1: A particle’s velocity and momentum share the same direction.

The momentum of a particle can be broken down into it’s components in the x, y and z directions;

px = mvx…..(3)

py = mvy…..(4)

pz = mvz…..(5)

where vx, vy and vz are the components of particle’s velocity.

Substituting equation (2) into equation (1), we get;

ΣF=dpdt\Sigma\vec{F} = \frac{d\vec{p}}{dt}…..(6)

Equation (6) says that the net external force acting on a particle it equal to the rate of change of it’s momentum with respect to time. This relationship between external forces acting on a body and it’s momentum is only applicable in inertial frames of reference since it’s derived from Newton’s laws of motion. The physical meaning of this equation is that a very large net external force corresponds to a rapid change in the particle’s momentum whereas a small net force corresponds to a gradual rate of change in particle’s momentum.

The Impulse-Momentum Theorem

Both momentum and kinetic energy depend on the particle’s mass and velocity so wherein lies the difference? Theoretically, momentum is a vector quantity whereas kinetic energy is scalar. Magnitude of momentum (p) is directly proportional to the speed (v) of the particle but the magnitude of kinetic energy (K) varies with the square of the speed (v2). In order to understand the physical difference between the two quantities, we need to define another quantity known as impulse.

Consider a particle that is under the influence of a constant net external force, ΣF\Sigma\vec{F}. If the net force acts over a time period given by Δt=t2t1\Delta t = t_2-t_1, then the impulse (J\vec{J} ) of the net external force is equal to the product of the net external force and the time period;

J=ΣFΔt\vec{J} = \Sigma \vec{F} \Delta t…..(7)

Impulse is a vector quantity and it acts in the same direction as the net external force. The SI unit of impulse is N ˙\dot{} s or kg ˙\dot{} m/s, which is the same as the momentum.

Now, going back to equation (6), if the net external force, ΣF\Sigma \vec{F} acting on a particle is constant, then the quantity dp/dtd\vec{p}/dt is also constant and thus, can be written as;

ΣF=dpdt=p2p1t2t1\Sigma\vec{F}=\frac{d\vec{p}}{dt} = \frac{\vec{p}_2-\vec{p}_1}{t_2-t_1}…..(8)

Equation (8) says that, as a consequence of constant net force, the instantaneous rate of change of momentum with respect to time is equal to the total change of momentum over the time interval Δt\Delta t divided by Δt=t2t1\Delta t = t_2-t_1.

Multiplying both sides of equation (8) with Δt\Delta t, we get;

ΣF(t2t1)=p2p1\Sigma\vec{F} (t_2-t_1)= \vec{p}_2-\vec{p}_1…..(9)

Equating equations (7) and (9), we see that;

J=p2p1\vec{J} = \vec{p}_2-\vec{p}_1…..(10)

The Impulse-Momentum Theorem states that the impulse of the net external force, ΣF\Sigma \vec{F} acting on a particle over the time interval Δt\Delta t is equal to the change in momentum of the particle during this interval.

If the force is not constant, equation (8) can integrated over the time interval Δt\Delta t;

t1t2ΣFdt=t1t2dpdtdt=p1p2dp=p2p1.....(11)\int_{t_1}^{t_2}\Sigma\vec{F}dt = \int_{t_1}^{t_2}\frac{d\vec{p}}{dt}dt=\int_{\vec{p}_1}^{\vec{p}_2}d\vec{p}=\vec{p}_2-\vec{p}_1…..(11)

Comparing equation (11) to equation (10), we see that;

J=t1t2ΣFdt.....(12)\vec{J} = \int_{t_1}^{t_2}\Sigma\vec{F}dt…..(12)

Even when net external force is not constant, equation (12) can be written as;

J=Fav(t2t1)\vec{J}=\vec{F}_{av}(t_2-t_1)…..(13)

where Fav\vec{F}_{av} is the average net external force acting on the body.

Note: If ΣF\Sigma\vec{F} is constant, Fav\vec{F}_{av} is simply equal to ΣF\Sigma\vec{F} and equation (13) becomes equivalent to equation (7).

To understand why the integral in equation (12) can be replaced with average net external force (Fav\vec{F}_{av}) acting on a body, see Figure 2.

Figure 2: The graph shows the x-component of the net external force acting on a body with respect to time. The area under the curve is equal to the impulse of the net external force.

In Figure 2, the blue curve represents a variable force that is changing with time. Since the force is not constant, the impulse (area under the curve) is given by;

Jx=t1t2ΣFxdt.....(14)J_x = \int_{t_1}^{t_2}\Sigma F_xdt…..(14)

The pink rectangle corresponds to the average net external force acting in the x-direction, and thus the impulse (area under the curve) is given as;

Jx=Favx(t2t1).....(15)J_x = F_{av-x}(t_2-t_1)…..(15)

The area under the blue curve and the area of the pink rectangle are roughly equal and thus, impulse can be simply calculated by replacing the time varying force, ΣFx\Sigma F_x with the average net external force, FavxF_{av-x}.

Aside: Since both impulse and momentum are vector quantities, they can be broken down into their component form;

Jx=t1t2ΣFxdt=Favx(t2t1)=p2xp1x=mv2xmv1x.....(16)J_x = \int_{t_1}^{t_2}\Sigma F_xdt = F_{av-x}(t_2-t_1) = p_{2x}-p_{1x} = mv_{2x}-mv_{1x}…..(16)
Jy=t1t2ΣFydt=Favy(t2t1)=p2yp1y=mv2ymv1y.....(17)J_y = \int_{t_1}^{t_2}\Sigma F_ydt = F_{av-y}(t_2-t_1) = p_{2y}-p_{1y} = mv_{2y}-mv_{1y}…..(17)

Also, note that a small force acting over a long time period produces an equivalent impulse to when a large force acts for a short duration (See Figure 3). In other words, the two forces of different magnitude will yield a similar effect on the particle.

Figure 3: Since the area under the purple and pink curve is roughly the same, the two forces produce an equal impulse.

A real-world application of impulse-momentum theorem would be airbags in a car. Initially, the driver has a momentum of p1=mv1\vec{p}_1 = m\vec{v}_1 where v1\vec{v}_1 is the velocity of the car and hence the driver. After the collision, as the car instantaneously comes to rest, driver’s momentum drops to zero, p2\vec{p}_2= 0.

According to impulse-momentum theorem;

J=p2p1.....(18)\vec{J} = \vec{p}_2-\vec{p}_1…..(18)

This means that no matter how the driver comes to rest, the impulse will remain the same.

However, instead of instantaneously coming to rest by hitting the steering wheel, airbags increase the time over which the driver comes to rest.

Since impulse is equal to;

J=FavΔt.....(19)\vec{J} = \vec{F}_{av}\Delta t…..(19)

Increasing the time interval reduces the amount of average force exerted on the driver and therefore, increases the chances of survival.

Momentum and Kinetic Energy Compared

According to the impulse-momentum theorem (J=p2p1)\vec{J}=\vec{p}_2-\vec{p}_1), changes in the momentum of a particle are caused by impulse, which depends on the time over which the net external force acts.

On the other hand, work-energy theorem (Wtot=K2K1)W_{tot}=K_2-K_1) tells us that changes in the kinetic energy are caused by work done on the particle, which depends on the distance over which the net external force acts.

Let’s say there is a particle at rest during time t = t1. Since the particle’s initial velocity, v1=0\vec{v}_1=0, it’s momentum at t1 is zero; p1=mv1=0\vec{p}_1=m\vec{v}_1=0, and it’s kinetic energy at t1 is also zero; K=12mv2=0K=\frac{1}{2}mv^2=0.

Let a constant net external force, F\vec{F} act on this particle such that at time t = t2, it has a velocity v2\vec{v}_2. During the time period Δt=t2t1\Delta t = t_2-t_1, the particle travels a distance s.

Now according to the impulse-momentum theorem, the momentum of the particle at time t2 is given by;

p2=J+p1.....(20)\vec{p}_2 = \vec{J}+\vec{p}_1…..(20)

But p1\vec{p}_1 = 0,

p2=J.....(21)\vec{p}_2=\vec{J}…..(21)

This means that momentum of a particle at any instant is equal to the impulse (J=FΔt)\vec{J} = \vec{F}\Delta t) needed to accelerate the said particle from rest to it’s speed at that instant.

On the other hand, the work-energy theorem gives the kinetic energy of the particle at time t = t2;

K2=Wtot+K1.....(22)K_2 = W_{tot}+K_1…..(22)

Since K1 = 0;

K2=Wtot.....(23)K_2 = W_{tot}…..(23)

Thus, kinetic energy of the particle is equal to the work done on the particle that accelerates it from rest.

Impulse is the product of force needed to accelerate the particle and the time period over which this acceleration occurs;

J=FΔt.....(24)\vec{J} = \vec{F} \Delta t…..(24)

whereas, work done on the particle is equal to the product of the force needed to accelerate the particle and the distance over which this acceleration occurs.

Wtot=Fs=Fs (cos0°)=Fs.....(25)W_{tot} = \vec{F} \cdot \vec{s}=Fs~(cos0°)= Fs…..(25)

A practical application of this concept would be catching either a 0.20 kg ball moving with a speed of 10 m/s or a 0.10 kg ball moving with a speed of 20 m/s.

Momentum (mv)Kinetic Energy (1/2mv2)
Ball 1(0.20kg)x(10m/s) = 2 kg\cdotm/s(1/2)(0.20kg)(10m/s)2= 10J
Ball 2(0.10kg)x(20m/s) = 2kg\cdotm/s(1/2)(0.10kg)(20m/s)2= 20J

The two balls have the same momentum, which means that the same amount of impulse in needed to bring them both to rest.

Since p1=2 kgm/s\vec{p}_1 = 2~kg \cdot m/s in both cases and p2\vec{p}_2= 0, impulse,J=p2p1=2 kgm/s\vec{J} = \vec{p}_2-\vec{p}_1 = -2~kg \cdot m/s.

However, the second ball has twice the amount of kinetic energy when compared to the first ball, which means that to catch the second ball, you need to do twice the amount of work to bring the ball to rest.

For ball 1, K1 = 10 J and K2=0; Wtot = K2 – K1 = -10 J

For ball 2, K1 = 20 J and K2=0; Wtot = K2 – K1 = -20 J

Thus, for a given amount of force that you exert on each ball to catch it, it will take you the same amount of time to stop the two balls, but the second ball will cause your hands to be pushed back twice as much with the ball covering twice as much the distance when compared to the first ball and hence, twice the amount of work will be done.


Conservation of Momentum (Young et al., 2020)

In a system with two or more objects that are interacting, the concept of momentum is very useful. To understand this, consider the following:

Types of Forces

  1. Internal Forces: The forces that particles in a system exert on each other are called internal forces.
  2. External Forces: The forces that are exerted on the system by some outside source are called external forces.

An isolated system is when there are no external forces acting on a system and the only forces present are internal to the system.

Now, consider an isolated system consisting of two particles A and B, such that they both exert a force on each other.

Figure 4: Particles A and B exert forces on each other in an isolated system.

Now, according to Newton’s third law, the forces that A and B exert on each other are equal in magnitude and opposite in direction. This means that the impulses of the forces exerted on each other and the changes in their momenta are also equal in magnitude and opposite in direction.

For particle A with net external force, FB on A\vec {F}_{B~on~A}, we can then say;

FB on A=dpAdt.....(26)\vec{F}_{B~on~A} = \frac{d\vec{p}_A}{dt}…..(26)

Similarly, for particle B;

FA on B=dpBdt.....(27)\vec{F}_{A~on~B} = \frac{d\vec{p}_B}{dt}…..(27)

Applying Newton’s third law;

FB on A=FA on B.....(28)\vec{F}_{B~on~A} = -\vec{F}_{A~on~B}…..(28)

Alternatively;

FB on A+FA on B=0…..(29)\vec{F}_{B~on~A} + \vec{F}_{A~on~B} = 0…..(29)

Substituting equations (26) and (27) into equation (29), we get;

FB on A+FA on B=dpAdt+dpBdt=d(pA+pB)dt=0…..(30)\vec{F}_{B~on~A} + \vec{F}_{A~on~B} = \frac{d\vec{p}_A}{dt} + \frac{d\vec{p}_B}{dt} = \frac{d(\vec{p}_A+\vec{p}_B)}{dt} = 0…..(30)

Since the rate of change of momentum with respect to time for particles A and B is equal and opposite, the rate of change of their sum with respect to time is zero.

We can now call P\vec{P}, the total momentum of the system;

P=pA+pB.....(31)\vec{P} = \vec{p}_A+\vec{p}_B…..(31)

Substituting equation (31) into equation (30), we get;

FB on A+FA on B=dPdt=0…..(30)\vec{F}_{B~on~A} + \vec{F}_{A~on~B} = \frac{d\vec{P}}{dt} = 0…..(30)

Thus, even though the individual momentum of the two particles is changing, the total momentum P\vec{P} is constant, or in other words, the total momentum of an isolated system is conserved.

If external forces are present,

ΣFA+ΣFB=dPdt.....(31)\Sigma \vec{F}_A + \Sigma \vec{F}_B = \frac{d\vec{P}}{dt}…..(31)

where

ΣFA=FB on A+FExternal on A.....(32)\Sigma \vec{F}_A = \vec{F}_{B~on~A} + \vec{F}_{External~on~A}…..(32)

ΣFB=FA on B+FExternal on B.....(33)\Sigma \vec{F}_B = \vec{F}_{A~on~B} + \vec{F}_{External~on~B}…..(33)

In general, when external forces are present, the total momentum P\vec{P} is not constant.

However, consider two blocks sitting on a surface;

Figure 5: Free-body diagrams for block A and B sitting on a surface.

This system is not isolated since there are external forces (normal force and weight) exerting on the two blocks.

Equation (31) then becomes;

For block A;

FB on A+nA+wA=dpAdt.....(34)\vec{F}_{B~on~A}+\vec{n}_A+\vec{w}_A = \frac{d\vec{p}_A}{dt}…..(34)

For block B;

FA on B+nB+wB=dpBdt.....(35)\vec{F}_{A~on~B}+\vec{n}_B+\vec{w}_B = \frac{d\vec{p}_B}{dt}…..(35)

But, nA=wA\vec{n}_A = -\vec{w}_A and nB=wB\vec{n}_B = -\vec{w}_B;

So then, equations (34) and (35) reduce to equations (26) and (27), which means that once again, the total momentum of the system is constant. This means that even when external forces are present, the total momentum of the system is conserved provided that the net external forces equal zero.

Thus, the law of conservation of momentum states that, “The total momentum of the system is conserved when the vector sum of all the external forces acting on a system is equal to zero“.

The law of conservation of momentum is useful because it can be applied to a system without having much knowledge about the nature of the internal forces acting between the particles of a system. Since Newton’s laws were used to derive the principle of conservation of momentum, it is applicable only in inertial frames of reference.

The law of conservation of momentum can be generalized for a system of n number of particles with the total momentum given by;

p=p1+p2+p3+.......+pn=mv1+mv2+mv3+.......+mvn.....(36)\vec{p} = \vec{p}_1+\vec{p}_2+\vec{p}_3+…….+\vec{p}_n = m\vec{v}_1+m\vec{v}_2+m\vec{v}_3+…….+m\vec{v}_n…..(36)

Since momentum is a vector quantity, p\vec{p} can be easily calculated using it’s component form;

px=p1x+p2x+p3x+......+pnx.....(37)p_x = p_{1x}+p_{2x}+p_{3x}+……+p_{nx}…..(37)

Similarly,

py=p1y+p2y+p3y+......+pny.....(38)p_y = p_{1y}+p_{2y}+p_{3y}+……+p_{ny}…..(38)
pz=p1z+p2z+p3z+......+pnz.....(39)p_z = p_{1z}+p_{2z}+p_{3z}+……+p_{nz}…..(39)

If the net external force acting on a system is zero, then px, py and pz are all constant2.


Momentum Conservation and Collisions (Young et al., 2020)

In physics, when objects interact strongly with one another over a short period of time, it is referred to as a collision. Some examples of collisions include: billiard balls striking each other on a pool table, car accidents, neutrons bombarding with atomic nuclei in a nuclear reactor etc.

For most collisions, the forces acting between the interacting objects are much larger than external forces acting on them and therefore, the colliding objects can be treated as an isolated system. This would mean that the total momentum before and after the collision will be conserved.

Elastic and Inelastic Collisions

When the forces acting between the colliding objects are conservative, the total mechanical energy of the system is conserved and the kinetic energy of the system before and after the collision is the same. Such an interaction between two objects is called an elastic collision.

For example, collision between two billiard balls can be considered almost completely elastic.

Consider two gliders A and B moving on a frictionless air table, such that each glider is mounted with an ideal spring bumper.

Figure 6: a) Initially the gliders are moving towards each other with velocities vA1\vec{v}_{A1} and vB1\vec{v}_{B1}. b) The gliders collide with each other such that the kinetic energy is transferred to elastic potential energy that is stored in the compressed springs while the gliders momentarily come to rest. c) The gliders bounce back in opposite directions after the springs expand and the elastic potential energy is converted back to kinetic energy.

In figure 6, since the spring’s obey Hooke’s Law, they exert conservative force, F = -kx, on each other. This means that the total mechanical energy of the system remains the same before, during and after the collision;

KA1+KB1=UA2+UB2=KA3+KB3.....(40)K_{A1}+K_{B1} = U_{A2} + U_{B2} = K_{A3} + K_{B3}…..(40)
12mAvA12+12mBvB12=12kAxA22+12kBxB22=12mAvA32+12mBvB32.....(41)\frac{1}{2}m_Av_{A1}^2+\frac{1}{2}m_Bv_{B1}^2=\frac{1}{2}k_Ax_{A2}^2+\frac{1}{2}k_Bx_{B2}^2=\frac{1}{2}m_Av_{A3}^2+\frac{1}{2}m_Bv_{B3}^2…..(41)

Note: xA1,xB1,vA2,vB2,xA3,xB3x_{A1}, x_{B1}, v_{A2}, v_{B2}, x_{A3}, x_{B3} are all zero in equation (41).

Since the total mechanical energy of the system remains conserved, we can say that the collision between the two gliders A and B is elastic.

When the total kinetic energy is lost after a collision, such that the total mechanical energy is not conserved, the interaction is called an inelastic collision. A fired bullet that gets drilled into a block of wood is an example of an inelastic collision. This is because the friction of the wood that brings the bullet to rest is a non-conservative force that causes the system to lose kinetic energy after the collision. In figure 6, if the frictionless air table is replaced with a rough surface with coefficient of kinetic friction given by μk\mu_k, the collision would then be inelastic.

If during an inelastic collision, the objects stick together and start to move together as one, such an interaction between the objects is referred to as a completely inelastic collision. In figure 6, if the spring bumpers were to be replaced with velcro strips, such that the two gliders stick together after the collision and start to move as one object, the collision would then be completely inelastic (see Figure 7).

Figure 7: The installed velcro (blue strips) makes the gliders stick together after the collision and the two move together as a single object with total kinetic energy that is less than the kinetic energy of the system before the collision.

Note: It is also possible for an object to have a final kinetic energy that is greater than it’s initial value, as is the case for rifle recoil.

Rules to Remember for Collisions

  1. If the net external force acting on the system is zero, the total momentum of the system before and after the collision is same3.
  2. If a collision is elastic, the total kinetic energy of the system before and after the collision is same.

Completely Inelastic Collisions

Consider the gliders A and B as shown in figure 7. Since the two gliders stick together in a completely inelastic collision, they both travel with the same velocity post-collision;

vA2=vB2=v2.....(42)\vec{v}_{A2} = \vec{v}_{B2} = \vec{v}_2…..(42)

Applying law of conservation of momentum;

mAvA1+mBvB1=(mA+mB)v2.....(43)m_A\vec{v}_{A1}+m_B\vec{v}_{B1} = (m_A+m_B)\vec{v}_2…..(43)

Equation (43) can be used to determine the common velocity after the collision for gliders A and B.

As an example, if an object A that is travelling with x-component of velocity given by vA1xv_{A1x}, collides inelastically with an object B which is at rest, equation (43) reduces to;

v2x=mAmA+mBvA1x.....(44)v_{2x} = \frac{m_A}{m_A+m_B}v_{A1x}…..(44)

If we compare the total kinetic energy before and after the collision;

K1=12mAvA1x2.....(45)K_1 = \frac{1}{2}m_Av_{A1x}^2…..(45)
K2=12(mA+mB)(mAmA+mB)2vA1x2.....(46)K_2 = \frac{1}{2}(m_A+m_B)\left(\frac{m_A}{m_A+m_B}\right)^2v_{A1x}^2…..(46)

Taking ratio;

K2K1=mAmA+mB.....(47)\frac{K_2}{K_1}=\frac{m_A}{m_A+m_B}…..(47)

Since mA+mB>mA,K2<K1m_A+m_B>m_A, K_2<K_1 and it can be confirmed that the kinetic energy post collision is less than the total initial kinetic energy of the system, as is the case in an inelastic collision.

Elastic Collisions (Young et al., 2020)

The forces between two colliding objects during an elastic collision are conservative. When two billiard balls strike, each ball deforms by a very small amount resulting in momentary conversion of some kinetic energy into elastic potential energy, which is then converted back to kinetic energy after the collision without any losses to the total mechanical energy of the system.

Consider two objects A and B undergoing an elastic collision in one dimensional space, say the x-direction.

Since the kinetic energy of the system before and after the collision is the same, we get;

12mAvA1x2+12mBvB1x2=12mAvA2x2+12mBvB2x2.....(48)\frac{1}{2}m_Av_{A1x}^2+\frac{1}{2}m_Bv_{B1x}^2 = \frac{1}{2}m_Av_{A2x}^2+\frac{1}{2}m_Bv_{B2x}^2…..(48)

And the conservation of momentum gives;

mAvA1x+mBvB1x=mAvA2x+mBvB2x.....(49)m_Av_{A1x}+m_Bv_{B1x} = m_Av_{A2x}+m_Bv_{B2x}…..(49)

If the masses (mA,mB)m_A, m_B) and initial velocities of the two objects (vA1x,vB1x)v_{A1x}, v_{B1x}) are known, equations (48) and (49) can be used to determine the final velocities of object A and B (vA2x,vB2x)v_{A2x}, v_{B2x}).

Elastic Collisions, One Object is at Rest

If object B is initially at rest, then equation (48) reduces to;

12mAvA1x2=12mAvA2x2+12mBvB2x2.....(48)\frac{1}{2}m_Av_{A1x}^2 = \frac{1}{2}m_Av_{A2x}^2+\frac{1}{2}m_Bv_{B2x}^2…..(48)

And equation (49) becomes;

mAvA1x=mAvA2x+mBvB2x.....(49)m_Av_{A1x} = m_Av_{A2x}+m_Bv_{B2x}…..(49)

Rearranging equations (48);

mAvA1x2=mAvA2x2+mBvB2x2.....(50)m_Av_{A1x}^2 = m_Av_{A2x}^2+m_Bv_{B2x}^2…..(50)
mAvA1x2mAvA2x2=mBvB2x2.....(51)m_Av_{A1x}^2 -m_Av_{A2x}^2 = m_Bv_{B2x}^2…..(51)
mA(vA1x2vA2x2)=mBvB2x2.....(52)m_A(v_{A1x}^2-v_{A2x}^2)=m_Bv_{B2x}^2…..(52)

Using the identity, a2b2=(a+b)(ab)a^2-b^2 = (a+b)(a-b),

mA(vA1x+vA2x)(vA1xvA2x)=mBvB2x2.....(53)m_A(v_{A1x}+v_{A2x})(v_{A1x}-v_{A2x})=m_Bv_{B2x}^2…..(53)

Similarly, rearranging equation (49);

mAvA1xmAvA2x=mBvB2x.....(54)m_Av_{A1x}-m_Av_{A2x}=m_Bv_{B2x}…..(54)
mA(vA1xvA2x)=mBvB2x.....(55)m_A(v_{A1x}-v_{A2x})=m_Bv_{B2x}…..(55)

Dividing equation (53) by equation (55), we get;

vA1x+vA2x=vB2x.....(56)v_{A1x}+v_{A2x}=v_{B2x}…..(56)

Now substituting vB2xv_{B2x} in equation (53) with the result from equation (56), we get;

mA(vA1x+vA2x)(vA1xvA2x)=mB(vA1x+vA2x)2.....(57)m_A(v_{A1x}+v_{A2x})(v_{A1x}-v_{A2x})=m_B(v_{A1x}+v_{A2x})^2…..(57)
mA(vA1xvA2x)=mB(vA1x+vA2x).....(58)m_A(v_{A1x}-v_{A2x})=m_B(v_{A1x}+v_{A2x})…..(58)

Opening the brackets and rearranging equation (58);

mAvA1xmAvA2x=mBvA1x+mBvA2x.....(59)m_Av_{A1x}-m_Av_{A2x}=m_Bv_{A1x}+m_Bv_{A2x}…..(59)
mAvA1xmBvA1x=mBvA2x+mAvA2x.....(60)m_Av_{A1x}-m_Bv_{A1x}=m_Bv_{A2x}+m_Av_{A2x}…..(60)
vA1x(mAmB)=vA2x(mA+mB).....(61)v_{A1x}(m_A-m_B)=v_{A2x}(m_A+m_B)…..(61)

Substituting the value of vA2xv_{A2x} from equation (62) into equation (56);

vB2x=vA1x+(mAmB)(mA+mB)vA1x.....(63)v_{B2x}=v_{A1x}+\frac{(m_A-m_B)}{(m_A+m_B)}v_{A1x}…..(63)

Equations (62) and (63) give a general result to calculate the final velocities of two objects A and B undergoing an elastic collision, given that the object B is initially at rest.

Let’s look at special cases:

Case 1: When mA>>>mBm_A>>>m_B,

Then mA+mBmAm_A+m_B \approx m_A and equation (62) becomes;

vA2x=(mAmB)mAvA1x=(1mBmA)vA1x.....(64)v_{A2x}=\frac{(m_A-m_B)}{m_A}v_{A1x} = \left(1-\frac{m_B}{m_A}\right)v_{A1x}…..(64)

Since mA>>>mBm_A>>>m_B, mA/mB<<1;m_A/m_B<<1;

vA2xvA1x.....(65)v_{A2x} \approx v_{A1x}…..(65)

Equation (63) reduces to;

vB2x2vA1x.....(66)v_{B2x} \approx 2v_{A1x}…..(66)

An example would be if a bowling ball moving with velocity vA1x\vec{v}_{A1x} were to collide with a stationary ping-pong ball. The bowling ball would keep moving with approximately the same velocity as if the collision never occurred. However, the ping-pong ball will move with a speed that is twice that of the bowling ball and in the same direction as the bowling ball (see Figure 8).

Figure 8: A bowling ball and a ping pong before and after an elastic collision given that the ping pong ball is initially at rest.

You might think that if the bowling ball is moving with approximately the same velocity as before, where is the ping pong ball getting the kinetic energy from that is causing it to move with twice the speed?

You need to remember that vA1xv_{A1x} is only approximately equal to vA2xv_{A2x}. The ping pong has so little mass compared to the bowling ball that it requires very little kinetic energy to gain twice the speed of the bowling ball.

Average mass of a bowling ball = 6 kg

Average mass of a ping pong ball = 0.0027 kg

Let’s say the bowling ball is initially moving with a speed of 6.7 m/s (vA1x=6.7 m/s)v_{A1x} = 6.7~m/s).

Then the speed attained by the ping pong ball after the collision, vB2x=13.4 m/sv_{B2x} = 13.4~m/s.

KEbowling ball=12(6kg)(6.7m/s)2=20.1 JKE_{bowling~ball} = \frac{1}{2}(6kg)(6.7m/s)^2 = 20.1~J
KEping pong ball=12(0.0027kg)(13.4m/s)2=0.0181 JKE_{ping~pong~ball} = \frac{1}{2}(0.0027kg)(13.4m/s)^2 = 0.0181~J

Clearly, the kinetic energy of the ping pong ball is insignificant compared to the kinetic energy of the bowling ball, even though the ping pong ball has twice the speed of the bowling ball.

Case 2: When mB>>>mAm_B>>>m_A,

Then mA+mBmBm_A+m_B \approx m_B and equation (62) becomes;

vA2x=(mAmB)mBvA1x=(mAmB1)vA1x.....(67)v_{A2x}=\frac{(m_A-m_B)}{m_B}v_{A1x} = \left(\frac{m_A}{m_B}-1\right)v_{A1x}…..(67)

Since mB>>>mAm_B>>>m_A, mA/mB<<1;m_A/m_B<<1;

vA2xvA1x.....(68)v_{A2x} \approx -v_{A1x}…..(68)

Equation (63) reduces to;

vB2x=2mAmBvA1x.....(69)v_{B2x}=\frac{2m_A}{m_B}v_{A1x}…..(69)

As shown earlier, mA/mB<<1;m_A/m_B<<1;

vB2x<<vA1x.....(70)v_{B2x} << v_{A1x}…..(70)

Now if the ping pong ball is moving towards the bowling ball with initial velocity vA1x\vec{v}_{A1x} such that the bowling ball is at rest, then the ping pong ball will strike the bowling ball and reverse its direction without much change to it’s magnitude of the velocity vector. On the other hand, the bowling ball will move in the same direction as vA1x\vec{v}_{A1x} but it’s magnitude of velocity will be much smaller compared to the ping pong ball (see Figure 9).

Figure 9: A ping pong ball and a bowling ball before and after an elastic collision given that the bowling ball is initially at rest.

Case 3: When mA=mBm_A = m_B,

Plugging this condition into equations (62) and (63), we get;

vA2x=0…..(71)v_{A2x} = 0…..(71)
vB2x=vA1x.....(72)v_{B2x} = v_{A1x}…..(72)

In case of two billiard balls on a pool table, when one ball moving with velocity vA1x\vec{v}_{A1x} strikes another ball that is initially at rest, the moving ball loses all of it’s momentum and kinetic energy, and passes it on to the stationary ball. Therefore, in the case when balls have equal mass, the first ball comes to a rest after the elastic collision while the second ball starts moving with the same velocity as the initial velocity of the first ball (see figure 10).

Figure 10: Two billiard balls of equal mass before and after an elastic collision, given that one ball is initially at rest.

Elastic Collisions and Relative Velocity

For two objects A and B that undergo an elastic collision, equation (56) can be written as;

vA1x=vB2xvA2x.....(73)v_{A1x} = v_{B2x}-v_{A2x}…..(73)

Now, vB2xvA2xv_{B2x}-v_{A2x} is the relative velocity of B with respect to A after the collision, which according to equation (73) is equal to vA1xv_{A1x}.

Before the collision,

vB1xvA1xv_{B1x} – v_{A1x} is the relative velocity of B with respect to A.

Since vB1x=0v_{B1x} = 0, we can say that the relative velocity of B with respect to A after the collision is equal to the negative of the relative velocity of B with respect to A before the collision.

Thus, the relative velocity before and after the collision has same magnitude but opposite direction. The direction of relative velocity before and after the collision is opposite to one another because before the collision, A and B are moving towards each other but after the collision, A and B start to drift apart.

Let’s say we were to observe the same collision from a different coordinate system that is moving with constant velocity relative to the first (stationary) coordinate system, then the velocities of A and B will be different (B will no longer have zero initial velocity with respect to the moving coordinate system), but their relative velocities will have the same magnitude and direction.

We can therefore, deduce a general statement for the relative velocities of A and B undergoing an elastic collision in one-dimensional space (regardless of whether one object is at rest or not):

For two objects undergoing an elastic collision in straight line motion, their relative velocities before and after collision have same magnitude but opposite direction.

vB2xvA2x=(vB1xvA1x).....(74)v_{B2x}-v_{A2x} = -\left(v_{B1x} – v_{A1x}\right)…..(74)

A similar relationship between relative velocities also exists in 3-D space,

The relative velocities of two objects undergoing an elastic collision have the same magnitude in three dimensional space before and after the collision.

If the above condition is met, we can say that the kinetic energy of the system before and after the collision will be conserved.


Center of Mass (Young et al., 2020)

A helpful way of using the principle of conservation of momentum is by describing it in terms of center of mass.

Consider a system of n particles with masses given as m1,m2,m3,.....,mnm_1, m_2, m_3,…..,m_n such that the location of m1m_1 is given by coordinates (x1,y1)(x_1, y_1), location of m2m_2 is given by coordinates (x2,y2)(x_2, y_2) and so on.

The center of mass of the system, (xcm,ycm)(x_{cm}, y_{cm}) is then defined as;

xcm=m1x1+m2x2+m3x3+.....+mnxnm1+m2+m3+.....+mn=Σi=1nmixiΣi=1nmi.....(75)x_{cm}=\frac{m_1x_1+m_2x_2+m_3x_3+…..+m_nx_n}{m_1+m_2+m_3+…..+m_n}=\frac{\Sigma_{i=1}^nm_ix_i}{\Sigma_{i=1}^n m_i}…..(75)
ycm=m1y1+m2y2+m3y3+.....+mnynm1+m2+m3+.....+mn=Σi=1nmiyiΣi=1nmi.....(76)y_{cm}=\frac{m_1y_1+m_2y_2+m_3y_3+…..+m_ny_n}{m_1+m_2+m_3+…..+m_n}=\frac{\Sigma_{i=1}^nm_iy_i}{\Sigma_{i=1}^n m_i}…..(76)

In vector form, the center of mass, (rcm)(\vec{r}_{cm}) can be written as;

rcm=m1r1+m2r2+m3r3+.....+mnrnm1+m2+m3+.....+mn=Σi=1nmiriΣi=1nmi.....(77)\vec{r}_{cm} = \frac{m_1\vec{r}_1+m_2\vec{r}_2+m_3\vec{r}_3+…..+m_n\vec{r}_n}{m_1+m_2+m_3+…..+m_n}=\frac{\Sigma_{i=1}^nm_i\vec{r}_i}{\Sigma_{i=1}^nm_i}…..(77)

Thus, center of mass is a “mass-weighted average position of the particles”.

Center of Mass of Solids

In order to determine the center of mass of solids that have a continuous distribution of particles, the sum in equation (75) and (76) needs to be switched to an integral. Solving these integrals can often get quite complex.

However, there are three general statements that can be made regarding the center of mass of solids:

  1. In case of homogeneous objects4 such as a sugar cube, a billiard ball or a can of fruit juice, the center of mass lies at the same location as the geometric center. For the homogeneous cube, sphere and cylinder shown below, the black dot corresponds to the geometric center which is also the center of mass of the object.
  2. If an object has an axis of symmetry5, the center of mass will always be located along this axis, as is the case for a wheel shown below:
  3. Center of mass (black dot) does not have to lie inside the object, as is the case for a donut shown below:

Motion of the Center of Mass

The significance of center of mass can be understood by examining what happens to it when the particles in the system are moving.

If we take the time derivative of equations (75) and (76), we get the x and y components of the velocity vector of the center of mass.

The x-component of the velocity vector can be calculated as follows;

vcmx=dxcmdt=ddt(m1x1+m2x2+m3x3+.....+mnxnm1+m2+m3+.....+mn)...v_{cm-x} = \frac{dx_{cm}}{dt} = \frac{d}{dt}\left(\frac{m_1x_1+m_2x_2+m_3x_3+…..+m_nx_n}{m_1+m_2+m_3+…..+m_n}\right)…
...=(1m1+m2+m3+.....+mn)ddt(m1x1+m2x2+m3x3+.....+mnxn)...…=\left(\frac{1}{m_1+m_2+m_3+…..+m_n}\right)\frac{d}{dt}\left(m_1x_1+m_2x_2+m_3x_3+…..+m_nx_n\right)…
...=1m1+m2+m3+.....+mn[m1(dx1dt)+m2(dx2dt)+m3(dx3dt).....+mn(dxndt)]...(78)…=\frac{1}{m_1+m_2+m_3+…..+m_n}\left[m_1\left(\frac{dx_1}{dt}\right)+m_2\left(\frac{dx_2}{dt}\right)+m_3\left(\frac{dx_3}{dt}\right)…..+m_n\left(\frac{dx_n}{dt}\right)\right]…(78)

Now dx1dt=vx1\frac{dx_1}{dt} = v_{x1}, dx2dt=vx2\frac{dx_2}{dt}=v_{x2} , dx3dt=vx3\frac{dx_3}{dt}=v_{x3},…….dxndt=vxn\frac{dx_n}{dt}=v_{xn};

vcmx=m1vx1+m2vx2+m3vx3+.....+mnvxnm1+m2+m3+.....+mn.....(79)v_{cm-x} = \frac{m_1v_{x1}+m_2v_{x2}+m_3v_{x3}+…..+m_nv_{xn}}{m_1+m_2+m_3+…..+m_n}…..(79)

Similarly, the y-component of the velocity vector is given by;

vcmy=m1vy1+m2vy2+m3vy3+.....+mnvynm1+m2+m3+.....+mn.....(80)v_{cm-y} = \frac{m_1v_{y1}+m_2v_{y2}+m_3v_{y3}+…..+m_nv_{yn}}{m_1+m_2+m_3+…..+m_n}…..(80)

The vector equation for the velocity of the center of mass is then given by;

vcm=m1v1+m2v2+m3v3+.....+mnvnm1+m2+m3+.....+mn.....(81)\vec{v}_{cm} = \frac{m_1\vec{v}_{1}+m_2\vec{v}_{2}+m_3\vec{v}_{3}+…..+m_n\vec{v}_{n}}{m_1+m_2+m_3+…..+m_n}…..(81)

Let M=m1+m2+m3+......+mnM = m_1+m_2+m_3+……+m_n be the total mass of the system, then equation (81) can be written as;

Mvcm=m1v1+m2v2+m3v3+.....+mnvn.....(82)M\vec{v}_{cm}=m_1\vec{v}_{1}+m_2\vec{v}_{2}+m_3\vec{v}_{3}+…..+m_n\vec{v}_{n}…..(82)

The sum of the momenta of each particle, m1v1+m2v2+m3v3+.....+mnvnm_1\vec{v}_{1}+m_2\vec{v}_{2}+m_3\vec{v}_{3}+…..+m_n\vec{v}_{n} equals the total momentum of the system, P\vec{P}.

We can now see that the total momentum of the system simply equals the product of total mass (M) and velocity of the center of mass (vcm\vec{v}_{cm});

P=Mvcm.....(83)\vec{P}=M\vec{v}_{cm}…..(83)

When a cricket batsman strikes the ball with his bat, he is striking a collection of particles that make up the ball. The impulse delivered to the ball causes each of these particles of mass m1,m2,m3.....m_1, m_2, m_3….. to move with a momentum given by m1v1,m2v2,m3v3.....m_1\vec{v}_1, m_2\vec{v}_2, m_3\vec{v}_3…..The total momentum of the system, m1v1+m2v2+m3v3m_1\vec{v}_1+m_2\vec{v}_2+m_3\vec{v}_3…..is the same as if the batsman were to hit a single particle of mass M=m1+m2+m3+....M=m_1+m_2+m_3+…. that moves with a velocity vcm\vec{v}_{cm} after the strike.

Thus, the concept of center of mass justifies the practice of treating a macroscopic object as a point particle.

The principle of conservation of momentum states that if the net external forces acting on a system equal zero, then the total momentum of the system is conserved, or in other words, P\vec{P} is constant.

Since vcm=P/M\vec{v}_{cm} = \vec{P}/M, we can say that if the net external forces equal zero, the velocity of the center of mass is also constant.

External Forces and Center-of-Mass Motion

When the net external forces acting on a system do not equal zero, the total momentum of the system is no longer constant. This means that the velocity of the center of mass is also not constant.

If we take the time derivative of equation (81), we get the acceleration of the center of mass;

acm=dvcmdt=ddt[m1v1+m2v2+m3v3+.....+mnvnm1+m2+m3+.....+mn]...\vec{a}_{cm} = \frac{d\vec{v}_{cm}}{dt}=\frac{d}{dt}\left[\frac{m_1\vec{v}_{1}+m_2\vec{v}_{2}+m_3\vec{v}_{3}+…..+m_n\vec{v}_{n}}{m_1+m_2+m_3+…..+m_n}\right]…
...1m1+m2+m3+...+mn[m1ddtv1+m2ddtv2+m3ddtv3+....+mnddtvn]...…\frac{1}{m_1+m_2+m_3+…+m_n}\left[m_1\frac{d}{dt}\vec{v}_1+m_2\frac{d}{dt}\vec{v}_2+m_3\frac{d}{dt}\vec{v}_3+….+m_n\frac{d}{dt}\vec{v}_n\right]…
=m1a1+m2a2+m3a3+....+mnanm1+m2+m3+....+mn.....(84)=\frac{m_1\vec{a}_1+m_2\vec{a}_2+m_3\vec{a}_3+….+m_n\vec{a}_n}{m_1+m_2+m_3+….+m_n}…..(84)

If total mass of the system, M=m1+m2+m3+......mnM=m_1+m_2+m_3+……m_n;

Macm=m1a1+m2a2+m3a3+....+mnan.....(85)M\vec{a}_{cm} = m_1\vec{a}_1+m_2\vec{a}_2+m_3\vec{a}_3+….+m_n\vec{a}_n…..(85)

According to Newton’s second law, m1a1m_1\vec{a}_1 is the vector sum of all the forces on particle 1, m2a2m_2\vec{a}_2 is the vector sum of all the forces on particle 2 and so on. This means that the vector sum of all the forces acting on all the particles in the system equals m1a1+m2a2+m3a3+.....+mnanm_1\vec{a}_1+m_2\vec{a}_2+m_3\vec{a}_3+…..+m_n\vec{a}_n;

ΣF=m1a1+m2a2+m3a3+....+mnan.....(86)\Sigma\vec{F} = m_1\vec{a}_1+m_2\vec{a}_2+m_3\vec{a}_3+….+m_n\vec{a}_n…..(86)

Breaking the net external force acting on all particles, ΣF\Sigma\vec{F} into it’s external (ΣFext\Sigma\vec{F}_{ext}) and internal (ΣFint\Sigma\vec{F}_{int}) components;

ΣF=ΣFext+ΣFint=m1a1+m2a2+m3a3+....+mnan.....(87)\Sigma\vec{F} = \Sigma\vec{F}_{ext}+\Sigma\vec{F}_{int}= m_1\vec{a}_1+m_2\vec{a}_2+m_3\vec{a}_3+….+m_n\vec{a}_n…..(87)

Combining equations (85) and (87), we see that;

ΣF=ΣFext+Fint=Macm.....(88)\Sigma\vec{F} = \Sigma\vec{F}_{ext}+\vec{F}_{int}=M\vec{a}_{cm}…..(88)

In accordance with Newton’s third law, the internal forces exist in pairs for a system of n particles. For example one such pair between particle 1 and 2 is given by;

F1 on 2=F2 on 1.....(89)\vec{F}_{1~on~2}=-\vec{F}_{2~on~1}…..(89)

Since the two forces in a pair are equal in magnitude but opposite in direction, the sum of all the internal forces (ΣFint\Sigma\vec{F}_{int}) equals zero;

Thus, equation (90) reveals that when external forces act on a group of particles in a system, the center of mass moves in such a manner as if all of the mass was condensed at the location of the center of mass and it was being acted on by a net external force that equals the vector sum of all the external forces acting on the particles6.

~o~

An example of motion of center of mass:

Consider a cannon shell travelling only under the influence of gravity such that air resistance is negligible. Since gravity is the only force acting on the cannon shell, it will travel in a parabolic trajectory. Now let’s say the shell breaks off mid-flight into two pieces of equal mass. After the shell breaks, the pieces fly off and continue to travel in their own parabolic trajectories. However, the center of mass of the two pieces still travels along the original path as though the total mass of the system was still concentrated at the point of the center of mass (see figure 11).

Figure 11: At t = 0s, a cannon shell is launched at an angle. Partway through it’s parabolic trajectory, the shell breaks into two equal pieces at t = 3s. After t = 3s, the two pieces follow their own trajectories but the center of mass (black dot) still travels along the original trajectory.

What happens when one piece of the cannon shell hits the ground?

As long as the fragments were in the air, the center of mass acted like a free-falling object and followed a parabolic trajectory. According to equation (90);

ΣFext=Mg=Macm....(91)\Sigma\vec{F}_{ext} = M\vec{g} = M\vec{a}_{cm}….(91)

When one fragment of the cannon shell hits the ground, normal force acts on it. As soon as the normal force enters the equation, the net external force acting on the system changes;

ΣFext=Mg+n=Macm....(92)\Sigma\vec{F}_{ext} = M\vec{g}+\vec{n} = M\vec{a}_{cm}….(92)

Thus, the acceleration of the center of mass is no longer equal to g\vec{g} and the center of mass will deviate from it’s parabolic trajectory.

~o~

A practical application of motion of center of mass is seen in astrophysics.

Consider the motion of earth and moon as a system. It is incorrect to say that the moon orbits around the earth. In reality, both the earth and the moon orbit their common center of mass (see figure 12).

Figure 12: The earth and the moon orbit their common center of mass.

This technique has been used to identify planets orbiting distant stars. When the planets are so faint that they cannot be detected with a telescope, the motion of the star is carefully studied. If the star seems to be wobbling around a point, it is due to the presence of a planet/planets (see figure 13).

Figure 13: The star and planet obiting their common center of mass causing the star to wobble, which can be detected from earth.

~o~

The motion of a system of particles can also be described in the following fashion by rewriting equation (85) as follows:

Macm=Mdvcmdt.....(93)M\vec{a}_{cm}=M\frac{d\vec{v}_{cm}}{dt}…..(93)

Since mass M is constant, it can be moved inside the derivative;

Macm=d(Mvcm)dt.....(94)M\vec{a}_{cm}=\frac{d(M\vec{v}_{cm})}{dt}…..(94)

From equation (83), we see that MvcmM\vec{v}_{cm} just equals the total momentum of the system (P\vec{P});

Macm=dPdt.....(95)M\vec{a}_{cm}=\frac{d\vec{P}}{dt}…..(95)

Substituting equation (95) into equation (90), we get;

ΣFext=dPdt.....(96)\Sigma\vec{F}_{ext}=\frac{d\vec{P}}{dt}…..(96)

Note: Equation (96) looks similar to equation (6) with the difference being that equation (96) is applicable for when dealing with a system of particles such as a macroscopic object whereas equation (4) is to be used for when there is only a single particle.

Equation (96) tells us that the total momentum of the system can only be changed by external forces. If there are any internal forces present, they can change the individual momenta of the particles but have no effect on the total momentum of the system. Lastly, if the net external force acting on the system is zero, then according to equation (96), the total momentum of the system (P\vec{P}) is constant and according to equation (95), the acceleration of the center of mass (acm\vec{a}_{cm}) is zero or in other words, the velocity of the center of mass (vcm\vec{v}_{cm}) is constant.


Rocket Propulsion (Young et al., 2020)

When mass of an object varies with time, the concept of momentum proves useful. Newton’s second law (ΣF=ma\Sigma\vec{F}=m\vec{a}) is not applicable in this case because mass (m) is not constant.

When fuel is expelled, it is done so by exerting a force on it in the backwards direction. Now according to Newton’s third law, the fuel exerts an equal and opposite force on the rocket that propels it in the forward direction. When the fuel is ejected at such large quantities, the mass of the rocket decreases with time. However, the mass of the system (rocket plus fuel) remains constant.

Consider a rocket of initial mass, m in outer space such that no external forces act on it (force of gravity and air resistance are absent). Let x-axis be along the direction in which the rocket is moving.

Figure 14: At time t, the rocket with mass m is travelling with velocity v. After a short while, at time t + dt, the rocket is now travelling with velocity v + dv and the velocity of the burned fuel is equal to v – vex. Burned fuel of mass of -dm is ejected from the rocket and the rocket’s new mass is given by m + dm. (Move the slider left and right to compare the two images)

As can be seen in figure 14, at time t, the rocket is moving along the positive x-direction with velocity v (relative to our coordinate system) and has mass m. Thus, the x-component of it’s momentum at time t is given by;

P1=mv.....(97)P_1 = mv…..(97)

In an infinitesimal interval dt, the rocket burns fuel of mass dm and gains a speed of dv. (Note: dm is inherently a negative quantity which makes sense because as the rocket loses fuel, m+dm < m).

During this interval, the burned fuel of mass -dm is ejected from the rocket. (Note: Since dm is inherently negative, -dm is then a positive quantity). Let the speed at which the burned fuel is ejected7 relative to the rocket be vex. This means that the magnitude of the velocity of the burned fuel relative to the rocket is vex. Since the burned fuel is ejected in a direction opposite to the motion of the rocket, the velocity of the burned fuel relative to the rocket is given by -vex.

Now relative to our coordinate system, the x-component of velocity of the burned fuel, vfuel is given by;

vfuel=v+(vex)=vvex.....(98)v_{fuel} = v+(-v_{ex})=v-v_{ex}…..(98)

To see how the above expression was derived, review the relative velocity discussion here;

Let’s denote our coordinate system by the letter C, the rocket by the letter R and fuel by the letter F, then the relative velocity of fuel with respect to our coordinate system is given as;

vF/C=vF/R+vR/C.....(99)v_{F/C}=v_{F/R}+v_{R/C}…..(99)

As discussed above, vF/R=vexv_{F/R} = -v_{ex} and vR/C=vv_{R/C} = v,

vfuel=vex+v.....(100)v_{fuel} = -v_{ex}+v…..(100)

or

vfuel=vvex.....(101)v_{fuel} = v-v_{ex}…..(101)

which is in agreement with equation (98).

And consequently, the x-component of the momentum of the ejected fuel is equal to;

pburned fuel=dm(vfuel)=dm(vvex).....(102)p_{burned~fuel} = -dm(v_{fuel}) = -dm(v-v_{ex})…..(102)

Now at time t+dt, the momentum of the rocket is given by (see figure 14);

procket+unburned fuel=(m+dm)(v+dv).....(103)p_{rocket+unburned~fuel}=(m+dm)(v+dv)…..(103)

Thus, the total momentum of the system (rocket + ejected fuel) is given by;

P2=pburned fuel+procket+unburned fuelP_2 = p_{burned~fuel}+p_{rocket+unburned~fuel}
P2=(m+dm)(v+dv)dm(vvex).....(104)P_2=(m+dm)(v+dv)-dm(v-v_{ex})…..(104)

Since, we assumed that there are no external forces acting on the system (rocket +fuel), then according to the law of conservation of momentum;

P1=P2.....(105)P_1 = P_2…..(105)
mv=(m+dm)(v+dv)dm(vvex).....(106)mv = (m+dm)(v+dv)-dm(v-v_{ex})…..(106)

Opening the brackets;

mv=mv+m dv+v dm+dm dvv dm+vex dm.....(107)mv = mv+m~dv+v~dm+dm~dv-v~dm+v_{ex}~dm…..(107)

Cancelling terms and rearranging;

m dv=dm dvvex dm.....(108)m~dv=-dm~dv-v_{ex}~dm…..(108)

Since dm and dv are infinitesimal changes in the mass and velocity of the rocket, their product dmdv is an even smaller quantity and can be neglected;

mdv=vexdm.....(109)mdv=-v_{ex}dm…..(109)

Dividing both sides by dt;

mdvdt=vexdmdt.....(110)m\frac{dv}{dt}=-v_{ex}\frac{dm}{dt}…..(110)

Now dv/dt equals the acceleration of the rocket;

ma=vexdmdt.....(111)ma = -v_{ex}\frac{dm}{dt}…..(111)

The left side of the equation (111) simply equals the net external force acting on the rocket alone (remember net external force on the rocket and fuel as a system is zero). This net external force on the rocket is referred to as thrust and we can denote it by the letter F;

F=vexdmdt.....(112)F = -v_{ex}\frac{dm}{dt}…..(112)

F is a positive quantity because -dm/dt is positive (as discussed before, dm is inherently negative because the rocket loses mass with time). A positive net external force makes sense because the rocket accelerates in the positive x-direction. Thrust is directly proportional to the exhaust speed, vexv_{ex} which means that faster is the fuel ejected relative to the rocket, larger is the force exerted on the rocket by the burned fuel. Additionally, thrust is also directly proportional to the rate of change of rocket mass with respect to time. This makes sense because as more fuel it burnt per unit of time, larger is the force on the rocket. Thus, an effective rocket burns fuel rapidly and ejects it at a high relative speed.

Using equation (111), the x-component of acceleration is given by;

a=dvdt=vexmdmdt.....(113)a =\frac{dv}{dt}= -\frac{v_{ex}}{m}\frac{dm}{dt}…..(113)

Again, the acceleration of the rocket, a is positive because both exhaust speed vexv_{ex} and -dm/dt are positive. Now if both vexv_{ex} and dm/dt are constant, the acceleration keeps increasing until the rocket runs out of fuel.

~o~

A common misunderstanding

It is common to think that a rocket has to push off the ground to ascend. If that were the case, rockets could not use fuel to accelerate in outer space because there is nothing to push against. However, this is not the case because a rocket is able to lift off by utilizing the law of conservation of momentum. This is why rocket propulsion is even more effective in outer space because there is likely no gravity or definitely no air resistance.

~o~

Now, if the exhaust speed, vexv_{ex} is constant, then equation (113) can be integrated with respect to time in order to find an expression for the velocity of the rocket at any given time;

dv=vexmdm.....(114)dv = -\frac{v_{ex}}{m}dm…..(114)

Changing the integration variables from v to v’ and m to m’ in equation (114) so that we can use v and m as upper limits of integration;

dv=vexmdm.....(115)dv’ = -\frac{v_{ex}}{m’}dm’…..(115)

Let the rocket have an initial mass of m0 and an initial velocity of v0. Integrating equation (115);

v0vdv=m0mvexmdm.....(116)\int_{v_0}^{v}dv’ = -\int_{m_0}^m\frac{v_{ex}}{m’}dm’…..(116)

Since the exhaust speed vexv_{ex} is constant;

v0vdv=vexm0m1mdm.....(117)\int_{v_0}^{v}dv’ = -v_{ex}\int_{m_0}^m\frac{1}{m’}dm’…..(117)
vv0=vex [ln (m)ln (m0)].....(118)v-v_0=-v_{ex}~[ln~(m)-ln~(m_0)]…..(118)
vv0=vex [ln (m0)ln (m)].....(119)v-v_0=v_{ex}~[ln~(m_0)-ln~(m)]…..(119)

Using the logarithmic identity: ln(a) – ln(b) = ln(a/b);

For the rocket to gain speed that is larger than the exhaust speed, vex, the term;

ln(m0m)>1…..(121)ln\left(\frac{m_0}{m}\right)>1…..(121)

The graph of ln(x) is shown below:

Figure 15: Graph of ln x with respect to x.

Clearly, for ln(m0/m) to be greater than 1, m0/m needs to be greater than 2.72 which equals the constant e (Euler’s number).

To maximize the speed attained by the rocket, m0/m needs to be as large as possible. This is done by making the initial mass of the rocket almost all fuel resulting in m0>>>m (m is the mass of the rocket after fuel has been burnt).


  1. Young, H.D. et al. (2020) Sears and Zemansky’s university physics: With modern physics. 15th edn. Boston: Pearson. ↩︎
  2. Law of conservation of momentum is even applicable in cases where the internal forces are not conservative, contrary to the law of conservation of mechanical energy. ↩︎
  3. Momentum conservation can even be used when external forces are present because oftentimes the internal forces between the colliding objects are much larger than the external forces acting on them. ↩︎
  4. Objects made of particles that have uniform composition and same physical properties. ↩︎
  5. An axis of symmetry is an imaginary line, a rotation around which would not affect the appearance of the object. ↩︎
  6. This result is the reason why in mechanics, we can safely treat a macroscopic object as a point particle. Additionally, it also explains why the motion of a macroscopic object can only be affected by external forces (internal forces cancel each other out). This is the reason why you cannot lift yourself up by pulling on your belt, the belt will exert a force that is equal and opposite to the force applied by your hands. ↩︎
  7. The speed at which a gas is ejected is referred to as the exhaust speed. ↩︎