 Torque-free tumbling

A useful readout of a tumble is the alignment of the intermediate principal axis ${\hat{e}}_2$ with the conserved angular momentum direction $\hat{L}$. In the principal body frame, write $L_i=I_i{\omega }_i$ and $L=|\mathbf{L}|$. ThenA change from $c_2\simeq +1$ to $c_2\simeq -1$ means that the directed body axis has turned nearly $180^{\circ }$ relative to $\mathbf{L}$. It does not mean that $\mathbf{L}$ reversed in the fixed space frame. The unsigned line through that principal axis hides this distinction, so a marked end matters in an actual demonstration.

For the normalized intermediate-axis model with $\mathbf{\omega }(0)=(0.01,1,0)$, the rotational kinetic energy and squared angular momentum are $E=1.00005$ and $L^2=4.0001$. Hence $L^2-2E{I}_2=-0.0001$: this initial condition lies close to, but off, the intermediate-axis separatrix. An independent calculation evolves ${\dot{\mathbf{L}}}_{\mathrm{body}}={\mathbf{L}}_{\mathrm{body}}\times I^{-1}{\mathbf{L}}_{\mathrm{body}}$ and a rotation matrix taking body components to space components. The intermediate principal axis crosses $c_2=0$ at normalized times approximately $10.3773$ and $31.1319$, alternating long nearly aligned intervals with quicker reversals.

Both calculations use step size $0.005$ without rescaling either invariant. Over $0\le t\le 50$, the tumbling calculation changes rotational kinetic energy and squared angular momentum by less than $4\times 10^{-13}$ relative to their initial values. Reconstructing $\mathbf{L}$ in the fixed space frame gives relative error below $3\times 10^{-12}$. Halving the step changes the body components and rotation matrix entries by less than $3\times 10^{-10}$. These checks distinguish physical tumbling from numerical drift. On the exact intermediate-axis separatrix, approach to an intermediate-axis rotation takes infinite time; the recurrent off-separatrix example here is different.