Euler's Equations and Torque-Free Precession

The Body-Fixed Perspective

Imagine a rigid body spinning in empty space with no external forces acting on it. If you observe it from a stationary point in the laboratory, the body tumbles and rotates in a complex manner. However, if you attach your coordinate system directly to the body itself, the moments of inertia become constants. The challenge lies in the fact that this body-fixed frame is rotating. When we calculate the rate of change of angular momentum in this rotating frame, we must account for the rotation of the axes themselves. This leads to a term involving the cross product of the angular velocity and the angular momentum. The resulting equations, known as Euler's equations, describe how the components of angular velocity change relative to the principal axes of the body. In the absence of torque, these equations dictate that the angular velocity vector does not generally align with the angular momentum vector, leading to a phenomenon called precession where the axis of rotation sweeps out a cone.

Euler's Equations of Motion

The fundamental law for rotational dynamics is \(\vec{\tau} = \frac{d\vec{L}}{dt}\). When expressed in the body frame, the time derivative includes a correction term for the rotating basis vectors. For a rigid body with principal moments of inertia \(I_1, I_2, I_3\), the components of the angular momentum are \(L_i = I_i \omega_i\). The general form of Euler's equations is given by:

These equations are cyclic in the indices \(1, 2, 3\). They show that the rate of change of one component of angular velocity depends on the products of the other two components. This coupling is the source of the complex tumbling motions observed in torque-free rigid bodies.

Torque-Free Precession of a Symmetric Top

Consider a symmetric top where \(I_1 = I_2 = I_\perp\) and \(I_3 = I_\parallel\). In the torque-free case (\(\vec{\tau} = 0\)), the third Euler equation simplifies significantly. Since \(I_1 = I_2\), the term \((I_1 - I_2)\omega_1 \omega_2\) vanishes, implying \(\dot{\omega}_3 = 0\). Thus, the component of angular velocity along the symmetry axis is constant. The first two equations become:

Differentiating the first equation with respect to time and substituting the second, we find that \(\omega_1\) and \(\omega_2\) oscillate harmonically. The angular velocity vector \(\vec{\omega}\) precesses about the symmetry axis (the 3-axis) with a precession frequency \(\Omega\). This frequency is given by:

This result indicates that the body-fixed angular velocity vector traces a cone around the symmetry axis. The rate of this precession depends on the difference in moments of inertia and the spin rate about the symmetry axis.

Stability of Rotation

A crucial insight from Euler's equations is the stability of rotation about different principal axes. If a rigid body is rotating about its principal axis of maximum moment of inertia, small perturbations result in bounded oscillations, meaning the rotation is stable. Similarly, rotation about the axis of minimum moment of inertia is also stable. However, rotation about the intermediate axis (where \(I_1 < I_2 < I_3\)) is unstable. Any small deviation from this axis will grow exponentially, causing the body to tumble away from the intermediate axis. This is known as the intermediate axis theorem. The stability analysis relies on the conservation of kinetic energy and angular momentum, which constrain the motion to the intersection of an ellipsoid (energy) and a sphere (angular momentum magnitude).