Interestingly, we simply have a different use of terminology here. In aircraft dynamics, for example, if a small perturbation isn't damped, then the aircraft isn't considered "stable."
If there's dissipation, small oscillations will damp out so long as they don't grow due to some dynamical instability. So the two uses are consistent.
So I woke up this morning troubled by the fact that in my view of rotating rigid bodies, the motions wouldn't necessarily be reversible. If one magically and perfectly reversed the rotation, the rigid body wouldn't eventually work its way back to the condition it was in when the perturbation was applied.
Precisely. Damping is an irreversible process that involves an increase in entropy, so it can only happen when the "rigid" body is coupled to many other degrees of freedom (such as molecules of fuel).
But I was definitely confused on rigid body rotation. Now I'm not entirely sure that I'm unconfused on the rotation of a body with significant quantities of viscous fluids, but . . . that may be a bit ambitious for a Saturday morning.
There's a beautiful geometric way to visualize the motion which this thread reminded me of. Attach a coordinate frame to the rigid body, with the origin at the center of mass and the axes aligned with the principal axes of the body. Now use that coordinate system. It's non-inertial, which means the angular momentum vector in those coordinates doesn't remain constant. But because these coordinates are related to inertial coordinates by a rotation, the
length of the angular momentum vector is constant. So is the total energy, of course. Therefore
[latex]${\rm total \, angular \, momentum}=\vec L^2 = L_x^2+L_y^2+L_z^2 ={\rm constant}$[/latex]
[latex]${\rm kinetic \, energy}=T={1 \over 2}\left( L_x^2/A+L_y^2/B+L_z^2/C \right) = {\rm constant}$[/latex]
The first equation is a sphere, the second is an ellipsoid. Therefore for some given angular momentum and energy, the motion will be along a line of intersection between those two shapes. That line of intersection is some closed curve on the sphere. It can be a small circle around the min axis or a small circle around the max axis. But it can't be a small circle around the mid axis, as you can see if you visualize the intersection - it must be a very large circle that only occasionally passes close to it.
So you can see the stability purely from that, which is pretty nice. No need to write any differential equations at all.