Physics Mechanics 100% Free Open Access
Chapter 5 โ€ข Theory & Derivations

Conservation of Linear Momentum

Center of mass for discrete and continuous systems, Tsiolkovsky rocket equation, variable-mass systems, and 1D/2D collisions in Laboratory and Center-of-Mass frames.

ยง5.1 Center of Mass of Discrete Systems and Continuous Bodies

For an extended body or a collection of $N$ interacting particles of masses $m_i$ at positions $\mathbf{r}_i$:

1. Discrete Particle Systems

The center of mass position vector $\mathbf{R}_{\text{cm}}$ is the mass-weighted average position: $$\mathbf{R}_{\text{cm}} = \frac{1}{M} \sum_{i=1}^N m_i \mathbf{r}_i, \quad M = \sum_{i=1}^N m_i$$

2. Continuous Mass Distributions

For a continuous body with mass density $\rho(\mathbf{r})$: $$\mathbf{R}_{\text{cm}} = \frac{1}{M} \iiint_V \mathbf{r} \, \rho(\mathbf{r}) \, dV, \quad M = \iiint_V \rho(\mathbf{r}) \, dV$$

3. Analytical Center of Mass for Standard Geometries

  • Uniform Semicircular Wire of Radius $R$: Let wire lie in $xy$-plane ($y \ge 0$). Linear density $\lambda = M/(\pi R)$. $dm = \lambda R d\theta$: $$y_{\text{cm}} = \frac{1}{M} \int_0^\pi (R \sin \theta) (\lambda R d\theta) = \frac{\lambda R^2}{M} [-\cos\theta]_0^\pi = \frac{2 R}{\pi} \approx 0.637 R$$
  • Uniform Semicircular Disc of Radius $R$: Surface density $\sigma = \frac{2M}{\pi R^2}$. Slice into concentric rings of radius $r$: $$y_{\text{cm}} = \frac{1}{M} \int_0^R \left(\frac{2r}{\pi}\right) (\pi r \sigma dr) = \frac{2\sigma}{M} \int_0^R r^2 dr = \frac{4 R}{3\pi} \approx 0.424 R$$
  • Uniform Solid Hemisphere of Radius $R$: Slice into horizontal discs of height $z$ ($z \in [0, R]$): $$z_{\text{cm}} = \frac{3 R}{8} = 0.375 R$$

ยง5.2 Multi-Particle Dynamics and Linear Momentum Conservation

The total linear momentum $\mathbf{P}$ of an $N$-particle system is: $$\mathbf{P} = \sum_{i=1}^N \mathbf{p}_i = \sum_{i=1}^N m_i \mathbf{v}_i = M \mathbf{v}_{\text{cm}}$$

1. Equation of Motion for the Center of Mass

Differentiating total momentum with respect to time: $$\frac{d\mathbf{P}}{dt} = M \mathbf{a}_{\text{cm}} = \sum_{i=1}^N \mathbf{F}_i^{\text{ext}} + \sum_{i=1}^N \sum_{j \neq i} \mathbf{F}_{ij}^{\text{int}}$$ By Newton's third law in strong form, mutual internal forces between particles cancel identically in pairs: $$\mathbf{F}_{ij}^{\text{int}} + \mathbf{F}_{ji}^{\text{int}} = \mathbf{0}$$ Therefore: $$\frac{d\mathbf{P}}{dt} = M \mathbf{a}_{\text{cm}} = \mathbf{F}_{\text{net}}^{\text{ext}}$$ Law of Conservation of Linear Momentum: If the net external force on a system vanishes ($\mathbf{F}_{\text{net}}^{\text{ext}} = \mathbf{0}$): $$\mathbf{P} = M \mathbf{v}_{\text{cm}} = \text{constant}$$ The center of mass moves with constant velocity regardless of complex internal explosions or collisions!

ยง5.3 Variable-Mass Dynamics and the Tsiolkovsky Rocket Equation

Newton's second law $\mathbf{F} = \frac{d\mathbf{p}}{dt}$ must be applied rigorously to open systems where mass continuously enters or leaves the control volume.

1. Derivation of the Rocket Equation of Motion

Consider a rocket of mass $m(t)$ moving with velocity $\mathbf{v}(t)$. During time $dt$, propellant mass $(-dm > 0)$ is ejected with exhaust velocity $\mathbf{u}_{\text{ex}}$ relative to the rocket. Linear momentum at time $t$: $p(t) = m v$. Linear momentum at time $t + dt$: $$p(t+dt) = (m + dm)(v + dv) + (-dm)(v - u_{\text{ex}}) = mv + m\,dv + u_{\text{ex}}\,dm$$ Change in momentum: $$dp = p(t+dt) - p(t) = m\,dv + u_{\text{ex}}\,dm$$ Applying external force $F_{\text{ext}}$: $$F_{\text{ext}} = \frac{dp}{dt} = m \frac{dv}{dt} + u_{\text{ex}} \frac{dm}{dt}$$ Rearranging gives the rocket dynamical equation: $$m \frac{dv}{dt} = -u_{\text{ex}} \frac{dm}{dt} + F_{\text{ext}} = T_{\text{thrust}} + F_{\text{ext}}$$ where $T_{\text{thrust}} = u_{\text{ex}}|\dot{m}|$ is the thrust force.

2. Free Space Tsiolkovsky Equation

In deep space with zero external forces ($F_{\text{ext}} = 0$): $$m dv = -u_{\text{ex}} dm \implies dv = -u_{\text{ex}} \frac{dm}{m}$$ Integrating from initial state $(m_0, v_0)$ to final burnout $(m_f, v_f)$: $$\Delta v = v_f - v_0 = u_{\text{ex}} \ln\left( \frac{m_0}{m_f} \right)$$

3. Vertical Ascent Under Gravity

For vertical launch against uniform gravity $g$: $$v(t) = -gt + u_{\text{ex}} \ln\left( \frac{m_0}{m_0 - \alpha t} \right)$$

ยง5.4 Collision Phenomena in Laboratory and Center-of-Mass Frames

Collisions are intense brief interactions where internal contact impulses vastly exceed external forces.

1. Classification by Kinetic Energy and Coefficient of Restitution

The coefficient of restitution $e$ along the line of impact is: $$e = -\frac{v_{2f} - v_{1f}}{v_{2i} - v_{1i}} = \frac{\text{relative separation speed}}{\text{relative approach speed}}$$
  • Elastic ($e = 1$): Total mechanical kinetic energy is conserved: $K_f = K_i$.
  • Inelastic ($0 < e < 1$): Energy is partially dissipated into deformation and heat: $K_f < K_i$.
  • Completely Inelastic ($e = 0$): Bodies coalesce and move together with velocity $\mathbf{v}_f = \mathbf{v}_{\text{cm}}$.

2. Laboratory vs Center-of-Mass (CM) Frame

In the CM frame, total momentum is identically zero: $m_1 \mathbf{u}_1 + m_2 \mathbf{u}_2 = \mathbf{0}$. In an elastic collision in the CM frame, the speeds of the particles are completely unchanged ($u_1' = u_1, u_2' = u_2$); the collision merely rotates the relative velocity vector by scattering angle $\theta^*$. The laboratory scattering angle $\theta$ relates to CM angle $\theta^*$ by: $$\tan \theta = \frac{\sin \theta^*}{\cos \theta^* + \frac{m_1}{m_2}}$$

๐Ÿ“ Chapter Worked Examples & Exercises

Complete derivations & analytical proofs
Medium Example 5.1: Two-Stage Rocket Burnout Velocity Optimization

A rocket has initial total mass $M_0$ and structural payload fraction $f_s = 0.10$ for each stage. The effective exhaust velocity is $u_{\text{ex}} = 3000\text{ m/s}$. Compare the final burnout velocity $\Delta v$ achieved by: (a) a single-stage rocket consuming $80\%$ of its mass in fuel, versus (b) a two-stage rocket with equal mass ratio per stage.

Step 1: Compute single-stage velocity
$$\text{Initial mass } m_0 = M_0, \quad \text{Final mass } m_f = 0.20 M_0$$ $$\Delta v_1 = u_{\text{ex}} \ln\left(\frac{M_0}{0.20 M_0}\right) = 3000 \ln(5) \approx 3000(1.6094) = 4828 \text{ m/s}$$

Single stage burnout reaches $4828\text{ m/s}$.

Step 2: Compute two-stage rocket velocity
$$\text{For two equal stages with stage mass ratio } \frac{m_{0,i}}{m_{f,i}} = 3.162$$ $$\Delta v_2 = 2 \times u_{\text{ex}} \ln(3.162) = 2 \times 3000 \times 1.151 \approx 6907 \text{ m/s}$$

Staging sheds empty structural mass, delivering over $2000\text{ m/s}$ higher final velocity for the same fuel mass.

Hard Example 5.2: Derivation of Center of Mass for a Solid Uniform Hemisphere

Derive by analytical volume integration the center of mass position $z_{\text{cm}}$ of a solid homogeneous hemisphere of radius $R$ and uniform mass density $\rho_0$, bounded by $z \ge 0$ and $x^2 + y^2 + z^2 \le R^2$.

Step 1: Slice the hemisphere into thin horizontal circular discs
$$\text{At height } z \in [0, R], \text{ radius of disc is } r(z) = \sqrt{R^2 - z^2}$$ $$dV = \pi r(z)^2 dz = \pi (R^2 - z^2) dz$$ $$M = \rho_0 \int_0^R \pi (R^2 - z^2) dz = \rho_0 \pi \left[ R^2 z - \frac{z^3}{3} \right]_0^R = \frac{2}{3}\pi \rho_0 R^3$$

This confirms the total volume of the hemisphere is $\frac{2}{3}\pi R^3$.

Step 2: Evaluate first moment of mass integral
$$\int z \, dm = \rho_0 \pi \int_0^R z(R^2 - z^2) dz = \rho_0 \pi \int_0^R (R^2 z - z^3) dz = \rho_0 \pi \left[ \frac{R^2 z^2}{2} - \frac{z^4}{4} \right]_0^R = \frac{1}{4}\pi \rho_0 R^4$$ $$z_{\text{cm}} = \frac{\int z \, dm}{M} = \frac{\frac{1}{4}\pi \rho_0 R^4}{\frac{2}{3}\pi \rho_0 R^3} = \frac{3}{8} R$$

The center of mass of a solid hemisphere lies exactly $3/8 R$ above the planar base.

Medium Example 5.3: Ballistic Pendulum Velocity Determination

A bullet of mass $m = 10\text{ g}$ is fired horizontally with unknown velocity $v_0$ into a wooden block of mass $M = 1.99\text{ kg}$ suspended by a light vertical cord of length $L = 2.0\text{ m}$. The bullet embeds completely into the block, and the combination swings upward through a vertical height $h = 10.0\text{ cm}$. Determine the bullet launch speed $v_0$.

Step 1: Apply momentum conservation during impact
$$m v_0 = (m + M) V \implies V = \frac{m}{m + M} v_0$$

Collision is completely inelastic, conserving linear momentum during the brief impact time.

Step 2: Apply energy conservation during subsequent swing
$$\frac{1}{2}(m + M) V^2 = (m + M) g h \implies V = \sqrt{2 g h} = \sqrt{2(9.80)(0.10)} = \sqrt{1.96} = 1.40 \text{ m/s}$$ $$v_0 = \left(\frac{m + M}{m}\right) V = \left(\frac{2.00}{0.010}\right) (1.40) = 200 \times 1.40 = 280 \text{ m/s}$$

The bullet launch velocity was exactly $280\text{ m/s}$.