Dynamics, Lectures 11 & 12: Oxford Mathematics 1st Year Student Lecture

Overview

Oxford Mathematics' lectures 11 & 12 delve into the dynamics of central force motion, demonstrating how Kepler's laws of planetary motion arise from Newton's inverse-square law of gravity. The lectures cover the derivation of conic section trajectories, the role of energy and angular momentum, and the calculation of deflection angles for comets. Finally, they introduce the concept of systems of particles, the center of mass, and Galilean transformations, laying the groundwork for understanding rigid body dynamics.

Key takeaways

The substitution u=1/r transforms the complex radial equation of motion for a central inverse-square force into a simple linear ODE, whose solutions are conic sections (ellipses, parabolas, hyperbolas).
The effective potential energy V_eff(r) = -k/r + (mh^2)/(2mr^2) provides a powerful tool to qualitatively understand orbital behavior (bounded vs. unbounded) based on total energy.
The center of mass of a system of particles moves as if it were a single particle acted upon only by external forces, as internal forces cancel out due to Newton's third law.
Galilean transformations (translations, rotations, constant velocity boosts) are the key to understanding how different inertial frames relate, ensuring Newton's laws remain valid.
Kepler's laws can be rigorously derived from Newton's laws of motion and universal gravitation, demonstrating the predictive power of Newtonian mechanics.
The deflection angle of a comet passing near the Sun depends on its velocity, impact parameter, and the gravitational constant, calculated using the derived trajectory equation.

Chapters

0:17 Central Force Motion and Polar Coordinates

Central force motion, specifically the inverse square law (gravity), implies planar motion.
Newton's second law in polar coordinates for a central force is analyzed.
The quantity r^2 * theta_dot is shown to be constant for central forces.

0:56 The u = 1/r Substitution

The substitution u = 1/r simplifies the radial equation of motion.
The transformed equation for u(theta) becomes a linear, second-order ODE.
For f = -k/r^2, the equation simplifies to u'' + u = constant.

2:22 General Solution and Conic Sections

The general solution for u(theta) involves cosine and sine terms, representing conic sections.
The solution is written as u(theta) = (k/mh^2) * (1 + e cos(theta)).
The parameter 'e' is identified as the eccentricity of the conic section.

3:28 Geometric Interpretation of Conic Sections

Conic sections (ellipse, parabola, hyperbola) are defined by a focus at the origin and a directrix.
The geometric definition: distance to focus = e * distance to directrix.
This geometric definition yields the same polar equation r = k*e / (1 + e cos(theta)).

8:29 Relating Time to Trajectory

The conservation of angular momentum (r^2 * theta_dot = h) relates time to the polar coordinates.
Integrating dt = (r^2 / h) d(theta) allows calculation of time intervals between angular positions.
This integral is complex but shows time can be recovered from the trajectory.

15:51 Conservation of Energy in Central Force Motion

The total energy E = 1/2 * m * |r_dot|^2 + V(r) is conserved.
The radial velocity term is expressed using r_dot and theta_dot.
The potential V(r) for f = -k/r is -k/r.

19:16 Effective Potential Energy

The energy equation is rewritten using an effective potential V_eff(r).
V_eff(r) = -k/r + (mh^2) / (2mr^2).
This form separates the radial kinetic energy (r_dot^2) from the potential terms.

22:17 Qualitative Analysis of Effective Potential

The shape of V_eff(r) has a singularity at r=0 and approaches 0 as r -> infinity.
The qualitative behavior of orbits (ellipse, parabola, hyperbola) is determined by the total energy relative to V_eff(r).
Starting below the 'axis' of V_eff corresponds to bounded (elliptical) orbits.

25:04 Circular Orbits and Equilibrium

A circular orbit occurs when the particle is at a radius where V_eff'(r) = 0.
This equilibrium radius is r0 = mh^2 / k.
Starting with zero velocity at r0 results in a stable circular orbit.

26:50 Geostationary Orbit Calculation

A geostationary orbit requires the satellite's period to match Earth's rotation (2 pi per day).
This implies a specific theta_dot value.
Solving for r0 using the angular momentum (h = r^2 * theta_dot) and period yields the required orbital radius of ~42,000 km.

33:21 Energy Expression and Eccentricity

The total energy E can be expressed in terms of the eccentricity 'e'.
E = (k^2 / 2mh^2) * (e^2 - 1).
E < 0 for bounded (elliptical) orbits, and E >= 0 for unbounded (parabolic/hyperbolic) orbits.

35:56 Comet Deflection Example: Setup

A comet of mass 'm' approaches from a large distance with speed 'v'.
The problem asks for the deflection angle caused by the Sun's gravity.
A key parameter 'p' is defined as the closest approach distance if the Sun were absent.

40:38 Translating Initial Conditions

The 'infinity' initial condition (large distance, speed v) needs to be translated into the polar coordinate system.
The initial velocity vector is decomposed into radial (r_dot) and tangential (r*theta_dot) components.
Using geometry, r_dot = -v*cos(alpha) and r*theta_dot = v*sin(alpha), where alpha is related to p and r.

50:30 Calculating Angular Momentum (h)

Angular momentum h = r^2 * theta_dot.
Substituting r*theta_dot = v*sin(alpha) and r from the geometry gives h = R * v * sin(alpha).
Using sin(alpha) = p/R, the angular momentum simplifies to h = vp.

55:40 Initial Conditions for u(theta)

As r -> infinity, alpha -> 0, and u = 1/r -> 0. So, u(0) = 0.
The radial velocity r_dot -> -v as r -> infinity.
Using the chain rule, du/d(theta) = (-1/h) * dr/dt, leading to du/d(theta)|_0 = v/h = 1/p.

1:03:33 Solving for u(theta) with Initial Conditions

The general solution u(theta) = A*cos(theta) + B*sin(theta) + C is used.
The conditions u(0)=0 and du/d(theta)|_0 = 1/p determine the constants A, B, and C.
The specific solution for the comet's trajectory is u(theta) = (1/p)*sin(theta) + (k/mh^2)*(1-cos(theta)).

1:09:06 Calculating the Angle of Deflection

The angle of deflection is related to the asymptotic directions of the hyperbola.
The asymptotic direction theta_star is found where u(theta_star) = 0.
Solving for theta_star using half-angle formulas leads to tan(theta_star/2) = -mpv^2/k.

1:23:26 Final Deflection Angle Formula

The angle of deflection delta is theta_star - pi.
The final formula for deflection is delta = pi - 2 * arctan(mpv^2/k).
This shows the deflection depends on the impact parameter 'p' and the comet's velocity 'v'.

1:30:46 Kepler's Laws of Planetary Motion

K1: Planetary paths are ellipses with the Sun at one focus.
K2: A line joining the Sun and planet sweeps out equal areas in equal times.
K3: The square of the orbital period is proportional to the cube of the semi-major axis.

1:38:31 Derivation of Kepler's Second Law

The area swept out in polar coordinates is given by integral(1/2 * r^2 d(theta)).
Area as a function of time is integral(1/2 * r^2 * (d(theta)/dt) dt).
Since r^2 * theta_dot = h (constant), the rate of area swept is constant (h/2), proving Kepler's Second Law.

1:44:01 Derivation of Kepler's Third Law

The period T = 2 * Area / (dA/dt) = 2 * (pi*a*b) / (h/2) = 4*pi*a*b / h.
Relating semi-major axis 'a' and semi-minor axis 'b' to r0 and eccentricity 'e'.
This leads to T^2 proportional to a^3, with the constant of proportionality involving m, h, and k.

Summary, takeaways, and chapters were generated by AI from the video's transcript and may contain errors. The video belongs to its creator, Oxford Mathematics.