Central Orbits
In mechanics, a central orbit describes the path taken by a particle that is attracted or repelled by a fixed point, known as the center of force. The force acting on the particle is always directed towards or away from this fixed center. This concept is fundamental to understanding planetary motion, the motion of satellites, and the behavior of charged particles in electric fields.
Nature of Central Force
A central force, denoted by F, can be expressed as a function of the distance 'r' from the center of force. It can be attractive (directed towards the center) or repulsive (directed away from the center). Mathematically, a central force can be written as:
F = f(r) r̂
where f(r) is a scalar function of the distance r, and r̂ is the unit vector pointing radially outwards from the center of force. If f(r) is negative, the force is attractive. If f(r) is positive, the force is repulsive.
Properties of Motion under Central Force
Motion under a central force exhibits several key properties:
- Conservation of Angular Momentum: The angular momentum of a particle moving under a central force about the center of force is constant. This means the orbit lies in a fixed plane passing through the center of force.
- Areal Velocity is Constant: The rate at which the radius vector sweeps out area is constant. This is a direct consequence of the conservation of angular momentum.
Differential Equation of Central Orbits
To derive the equation of the orbit, we use polar coordinates (r, θ). Let the position vector of the particle be r. The angular momentum L is given by L = r × p, where p = mv is the linear momentum.
Since the force is central, F is parallel to r. Therefore, the torque τ = r × F = 0. As torque is the rate of change of angular momentum (τ = dL/dt), this implies that L is constant.
The areal velocity (dA/dt) is given by (1/2) |r × v| = (1/2) |L|/m. Since L is constant, the areal velocity is also constant. Let this constant be 'h'. So, dA/dt = h/2.
In polar coordinates, dA = (1/2) r² dθ. Therefore, dA/dt = (1/2) r² (dθ/dt) = h/2. This gives us:
r² (dθ/dt) = h
Now, let's consider the acceleration in polar coordinates. The velocity vector is v = ṙ r̂ + rθ̇ θ̂, and the acceleration vector is a = (r̈ - rθ̇²) r̂ + (rθ̈ + 2ṙθ̇) θ̂.
The equation of motion is ma = F = f(r) r̂.
Equating the radial components: m(r̈ - rθ̇²) = f(r).
We want to express this equation in terms of u = 1/r and θ.
From r² (dθ/dt) = h, we have dθ/dt = h/r² = hu².
Now, let's find dr/dt:
dr/dt = (dr/dθ) (dθ/dt) = (dr/dθ) (h/r²)
Let's change variables to u = 1/r. Then r = 1/u, and dr/dθ = d(1/u)/dθ = (-1/u²) (du/dθ).
So, dr/dt = (-1/u²) (du/dθ) (hu²) = -h (du/dθ).
Now let's find the second derivative r̈:
r̈ = d/dt (dr/dt) = d/dt [-h (du/dθ)] = d/dθ [-h (du/dθ)] (dθ/dt)
r̈ = [-h (d²u/dθ²)] (h/r²) = -h²u² (d²u/dθ²)
Substitute r̈ and rθ̇² into the equation of motion:
m(-h²u² (d²u/dθ²) - (1/u)(hu²)²) = f(1/u)
m(-h²u² (d²u/dθ²) - hu²) = f(1/u)
Divide by -mh²u²:
(d²u/dθ²) + u = -f(1/u) / (mh²u²)
This is the general differential equation for a central orbit.
Kepler's Laws of Planetary Motion
Kepler's Laws are a set of rules describing the motion of planets around the Sun, which is a special case of central force motion (specifically, an inverse square force).
First Law: The Law of Ellipses
"The orbit of every planet is an ellipse with the Sun at one of the two foci."
This law is derived from the differential equation when f(r) = -GMm/r², where G is the gravitational constant and M is the mass of the Sun, m is the mass of the planet.
The equation becomes: (d²u/dθ²) + u = (GMm/r²) / (mh²u²) = GM / (h²u²)
(d²u/dθ²) + u = GM/h²
The solution to this differential equation is:
u = (GM/h²) (1 + e cos(θ - θ₀))
or
r = 1/u = [h²/(GM)] / (1 + e cos(θ - θ₀))
This is the polar equation of a conic section. For planetary orbits, e < 1, which corresponds to an ellipse. The semi-latus rectum is p = h²/GM, and the eccentricity is e.
Second Law: The Law of Areas
"A line segment joining a planet and the Sun sweeps out equal areas during equal intervals of time."
This law is a direct consequence of the conservation of angular momentum, as we derived earlier. The areal velocity dA/dt = h/2, which is constant.
Third Law: The Law of Periods
"The square of the orbital period of a planet is directly proportional to the cube of the semi-major axis of its orbit."
For an elliptical orbit with semi-major axis 'a', the period 'T' is given by:
T² ∝ a³
For an inverse square force like gravity, the constant of proportionality is 4π² / (GM), where M is the mass of the central body. So, T² = (4π² / (GM)) a³.
For a circular orbit of radius 'r', a = r, and the velocity v = √(GM/r). The period T = 2πr / v = 2πr / √(GM/r) = 2πr√(r/GM) = 2π√(r³/GM). Squaring this gives T² = (4π²/GM) r³, which matches the third law for a = r.
Types of Orbits
Depending on the energy of the particle and the nature of the force, the orbits can be elliptical, parabolic, or hyperbolic. These correspond to the conic sections.
- Elliptical Orbits (e < 1): Bound orbits, typically seen for planets and satellites under gravity. The particle has negative total energy.
- Parabolic Orbits (e = 1): Unbound orbits, where the particle has just enough energy to escape the central force. The total energy is zero.
- Hyperbolic Orbits (e > 1): Unbound orbits, where the particle has more than enough energy to escape. The total energy is positive.
Effective Potential Energy
For a particle of mass 'm' moving under a central force F = f(r) r̂, the force can be derived from a potential energy function V(r) such that f(r) = -dV/dr.
The total energy E = Kinetic Energy + Potential Energy = (1/2)m(v<0xE1><0xB5><0xA3>² + v<0xE2><0x82><0x9B>²) + V(r).
Using v<0xE1><0xB5><0xA3> = ṙ and v<0xE2><0x82><0x9B> = rθ̇, and noting that rθ̇ = h/r, we get:
E = (1/2)m(ṙ² + (h/r)²) + V(r)
We can define an 'effective potential energy' U(r) that accounts for both the actual potential energy and the angular momentum term:
U(r) = V(r) + (1/2)m(h²/r²)
Then, the energy equation becomes E = (1/2)mṙ² + U(r). This form is analogous to the one-dimensional motion of a particle under a potential U(r), where the radial velocity ṙ plays the role of velocity.
The term (1/2)m(h²/r²) is often called the "centrifugal potential energy" or "angular momentum barrier".
The motion is confined to regions where E ≥ U(r). Turning points occur when E = U(r), which means ṙ = 0.
For an attractive inverse square force, V(r) = -GMm/r.
U(r) = -GMm/r + mh²/2r².
The shape of U(r) determines the nature of the orbit.
Moment of Inertia
The moment of inertia (I) is a measure of an object's resistance to changes in its rotation. It is the rotational analog of mass. Just as mass resists linear acceleration, moment of inertia resists angular acceleration. It depends not only on the mass of the object but also on how that mass is distributed relative to the axis of rotation.
Definition
For a system of discrete particles, the moment of inertia is the sum of the product of each particle's mass and the square of its perpendicular distance from the axis of rotation.
I = Σ mᵢ rᵢ²
where mᵢ is the mass of the i-th particle and rᵢ is its perpendicular distance from the axis of rotation.
For a continuous body, the moment of inertia is calculated by integrating over the entire mass distribution:
I = ∫ r² dm
where dm is an infinitesimal mass element and r is its perpendicular distance from the axis of rotation.
Units and Dimensions
The SI unit for moment of inertia is kilogram meter squared (kg⋅m²). The dimensions are [M L²].
Factors Affecting Moment of Inertia
- Mass of the object: More mass generally means larger moment of inertia.
- Distribution of mass: Mass concentrated further from the axis of rotation increases the moment of inertia more significantly than mass concentrated closer to the axis.
- Axis of rotation: The moment of inertia is specific to a chosen axis of rotation. Changing the axis changes the moment of inertia.
Parallel Axis Theorem
The parallel axis theorem relates the moment of inertia of a body about an axis to its moment of inertia about a parallel axis passing through its center of mass.
If I<0xE1><0xB5><0xA_> is the moment of inertia about an axis passing through the center of mass, and I is the moment of inertia about a parallel axis at a distance 'd' from the center of mass axis, then:
I = I<0xE1><0xB5><0xA_> + Md²
where M is the total mass of the body.
Perpendicular Axis Theorem
This theorem applies only to planar (2D) objects. It states that the moment of inertia of a planar object about an axis perpendicular to its plane and passing through a point O is equal to the sum of the moments of inertia about two perpendicular axes lying in the plane and intersecting at O.
If the planar object lies in the xy-plane, and the axis of rotation is the z-axis, then:
I<0xE2><0x82><0x9B> = I<0xE1><0xB5><0xBD> + I<0xE1><0xB5><0xB8>
where I<0xE2><0x82><0x9B>, I<0xE1><0xB5><0xBD>, and I<0xE1><0xB5><0xB8> are the moments of inertia about the z, x, and y axes, respectively.
Moment of Inertia of Some Common Objects
It is crucial to remember the moments of inertia for standard shapes, as these are frequently used in problems.
| Object | Axis of Rotation | Moment of Inertia (I) |
|---|---|---|
| Thin Rod | Through center, perpendicular to length | (1/12)ML² |
| Thin Rod | Through one end, perpendicular to length | (1/3)ML² |
| Annulus (Ring) | Through center, in the plane of the ring | MR² |
| Annulus (Ring) | Through center, perpendicular to the plane | MR² |
| Solid Cylinder/Disk | Through center, perpendicular to the plane (axis of symmetry) | (1/2)MR² |
| Hollow Cylinder (Thin Wall) | Through center, perpendicular to the plane (axis of symmetry) | MR² |
| Solid Sphere | Through center | (2/5)MR² |
| Hollow Sphere (Thin Wall) | Through center | (2/3)MR² |
Where M is the mass and R is the relevant radius or length.
Relation between Moment of Inertia and Torque
Similar to how force causes linear acceleration (F = ma), torque causes angular acceleration (τ = Iα).
τ = Iα
where τ is the net torque acting on the object, I is the moment of inertia about the axis of rotation, and α is the angular acceleration (α = d²θ/dt²).
Angular Momentum and Kinetic Energy
The angular momentum (L) of a rigid body rotating about a fixed axis is given by:
L = Iω
where ω is the angular velocity.
The rotational kinetic energy (K<0xE1><0xB5><0xA3>) is given by:
K<0xE1><0xB5><0xA3> = (1/2)Iω² = L² / (2I)
Applications of Moment of Inertia
Moment of inertia is a crucial concept in many areas of physics and engineering:
- Rotational Dynamics: Analyzing the motion of spinning objects like wheels, gyroscopes, and motors.
- Collisions: Understanding how angular momentum is conserved in rotational collisions.
- Engineering Design: Designing rotating machinery, vehicles, and structures where rotational stability is important.
- Astronomy: Calculating the rotation of celestial bodies.
Example Calculation: Moment of Inertia of a Uniform Rod about its Center
Consider a uniform rod of mass M and length L. We want to find its moment of inertia about an axis perpendicular to the rod and passing through its center.
We use the formula I = ∫ r² dm.
Let the rod lie along the x-axis, with the center at x=0. The ends are at x = -L/2 and x = +L/2.
The linear mass density is λ = M/L.
An infinitesimal mass element dm at position x has mass dm = λ dx = (M/L) dx.
The distance of this mass element from the axis of rotation (at x=0) is r = |x|. So, r² = x².
Now, we integrate:
I = ∫<0xE1><0xB5><0xA3>₋<0xE1><0xB5><0x87>/²<0xE2><0x81><0xBB><0xE1><0xB5><0x87>/² x² dm
I = ∫<0xE1><0xB5><0xA3>₋<0xE1><0xB5><0x87>/²<0xE2><0x81><0xBB><0xE1><0xB5><0x87>/² x² (M/L) dx
I = (M/L) ∫<0xE1><0xB5><0xA3>₋<0xE1><0xB5><0x87>/²<0xE2><0x81><0xBB><0xE1><0xB5><0x87>/² x² dx
The integral of x² is x³/3.
I = (M/L) [x³/3]<0xE1><0xB5><0xA3>₋<0xE1><0xB5><0x87>/²<0xE2><0x81><0xBB><0xE1><0xB5><0x87>/²
I = (M/L) [((L/2)³/3) - ((-L/2)³/3)]
I = (M/L) [(L³/8)/3 - (-L³/8)/3]
I = (M/L) [L³/24 + L³/24]
I = (M/L) [2L³/24]
I = (M/L) [L³/12]
I = (1/12)ML²
This matches the value in the table. This step-by-step process is how moments of inertia for continuous bodies are derived.