Past exam of the mathematics course of the University of Cambridge / 2016 / ia / Paper 4 / 12B / Solution

Use an Earth-centred rotating reference frame, with Earth's constant angular velocity vector $\mathbf{\omega }$. The vectors $\mathbf{r}$, $\mathbf{v}$, and $\mathbf{a}$ are the particle's position, velocity, and acceleration measured in that rotating frame. The vector $\mathbf{g}$ is the gravitational acceleration due to the actual gravitational force, expressed in the same axes; it excludes the separate centrifugal acceleration. The remaining term $-2\mathbf{\omega }\times \mathbf{v}$ is the Coriolis acceleration. A constant rotation axis and an unaccelerated Earth-centred origin account for the absence of an Euler term or a translational-origin term.

For the size estimates, Earth's angular speed isWith drop height $H=20\mathrm{m}$, unperturbed free fall gives the fall timeand final speed $gT=20{\mathrm{m}\mathrm{s}}^{-1}$. Thus the gravitational term has size $10{\mathrm{m}\mathrm{s}}^{-2}$, the Coriolis acceleration is at most about $2\omega gT\approx 2.8\times 10^{-3}{\mathrm{m}\mathrm{s}}^{-2}$, and the centrifugal acceleration is at most ${\omega }^2{R}_0\approx 2.9\times 10^{-2}{\mathrm{m}\mathrm{s}}^{-2}$. These are characteristic sizes $10^1$, $10^{-3}$, and $10^{-2}$ in acceleration units; geometric factors can reduce them, and the Coriolis term is zero at release.

For a suitable approximation to the rotational deflection of a falling body, take $\mathbf{g}$ constant over the short fall, start at rest relative to Earth, and use the zeroth-order velocity ${\mathbf{v}}^{(0)}(t)=\mathbf{g}t$ inside the small Coriolis term. In the centrifugal term replace $\mathbf{r}$ by ${\mathbf{R}}_0$, because $H/R_0\ll 1$. Integrating the resulting acceleration twice, with initial position ${\mathbf{r}}_{\mathrm{i}}$, givesSince ${\mathbf{r}}_{\mathrm{i}}+\mathbf{g}{T}^2/2={\mathbf{R}}_0$, evaluating at the unperturbed fall time yields the printed first-iterate estimate

The approximation keeps the leading Coriolis contribution and the leading centrifugal contribution enhanced by Earth's large radius. Feeding the Coriolis velocity correction back into the equation produces terms of size ${\omega }^{2}g{T}^4$, smaller than the retained centrifugal displacement by a factor of order $H/R_0$. Replacing $\mathbf{r}$ by ${\mathbf{R}}_0$ in that term has the same relative error. Here $\omega T\approx 1.4\times 10^{-4}$ and $H/R_0\approx 3.3\times 10^{-6}$, so both approximations are well controlled. Strictly, the displayed vector is the position at time $T$: its small upward component is removed by the slightly longer actual fall time. That arrival-time correction changes the retained horizontal displacements only at higher order. Thus the horizontal components of the printed estimate give the landing point on the locally flat surface; its vertical component should not be interpreted as a landing above the ground.

Let ${\hat{\mathbf{e}}}_E,{\hat{\mathbf{e}}}_N,{\hat{\mathbf{e}}}_U$ point east, north, and vertically upwards. They form a right-handed basis, with ${\hat{\mathbf{e}}}_N\times {\hat{\mathbf{e}}}_U={\hat{\mathbf{e}}}_E$. At latitude $\lambda$, in the spherical-Earth approximation,The required cross products areAt $\lambda =45^{\circ }$, the northerly and easterly displacements are thereforeThe negative northerly component means a southward displacement. With the supplied values, these are approximately $2.9\mathrm{cm}$ south and $1.3\mathrm{mm}$ east. The eastward effect comes from Coriolis acceleration, while the leading southward effect here comes from centrifugal acceleration.