# Derivation of the centrifugal and coriolis force

I was wondering how easily these two pseudo-forces can be derived mathematically in order to exhibit a clear physical meaning.

How would you proceed?

Okay, here is my (hopefully rigorous) demonstration of the origin of these forces here, from first principles. I've tried to be pretty clear what's happening with the maths. Bear with me, it's a bit lengthy!

### Angular velocity vector

Let us start with the principal equation defining angular velocity in three dimensions,

$$\dot{\mathbf{r}} = \mathbf{\omega} \times \mathbf{r}\; .$$

(This can be derived roughly by considering a centripetal force acting on a particle. Note that this equation applies symmetrically in inertial and rotating reference frames.)

Notice that we can in fact generalise this statement in terms of $r$ for an arbitrary vector $a$ that is known to be fixed in the rotating body.

### Transformation between inertial and rotating frames

Now consider a vector $a$, which we can write in Cartesian coordinates (fixed within the body) as

$$\mathbf{a} = a_x \mathbf{\hat i} + a_y \mathbf{\hat j} + a_z \mathbf{\hat k}\; .$$

In Newtonian mechanics, scalar quantities must be invariant for any given choice of frame, so we can say

$$\left.\frac{\mathrm da_x}{\mathrm dt}\right|_I = \left.\frac{\mathrm da_x}{\mathrm dt}\right|_R$$

where $I$ indicates the value is for the inertial frame, and $R$ that the value is for the rotating frame. Equivalent statements apply for $a_y$ and $a_z$, of course. Hence, any transformation of $a$ between frames must be due to changes in the unit vectors of the basis.

Now by the product rule,

\begin{align}\left.\frac{\mathrm d\mathbf{a}}{\mathrm dt}\right|_I &= \frac{\mathrm d}{\mathrm dt} \left( a_x \mathbf{\hat i} + a_y \mathbf{\hat j} + a_z \mathbf{\hat k} \right) \\& = \left( \frac{\mathrm da_x}{\mathrm dt} \mathbf{\hat i} + \frac{\mathrm da_y}{\mathrm dt} \mathbf{\hat j} + \frac{\mathrm da_z}{\mathrm dt} \mathbf{\hat k} \right) + \left( a_x \frac{\mathrm d\mathbf{\hat i}}{\mathrm dt} + a_y \frac{\mathrm d\mathbf{\hat j}}{\mathrm dt} + a_z \frac{\mathrm d\mathbf{\hat k}}{\mathrm dt} \right) .\end{align}

Using the previous equation for angular velocity, we then have

\begin{align}\left.\frac{\mathrm d\mathbf{a}}{\mathrm dt}\right|_I &= \left( \frac{\mathrm da_x}{\mathrm dt} \mathbf{\hat i} + \frac{\mathrm da_y}{\mathrm dt} \mathbf{\hat j} + \frac{\mathrm da_z}{\mathrm dt} \mathbf{\hat k} \right) + \left( a_x \mathbf{\omega} \times \mathbf{\hat i} + a_y \mathbf{\omega} \times \mathbf{\hat j} + a_z \mathbf{\omega} \times \mathbf{\hat k} \right) \\&= \left.\frac{\mathrm d\mathrm{a}}{\mathrm dt}\right|_R + \mathbf{\omega} \times \mathbf{a} \;.\end{align}

Now consider a position vector on the surface of a rotating body. We can write

$$\mathbf{v}_I = \left.\frac{\mathrm d\mathbf{r}}{\mathrm dt}\right|_I = \left.\frac{\mathrm d\mathbf{r}}{\mathbf dt}\right|_R + \mathbf{\omega} \times \mathbf{r} ,$$

and similarly for $\mathbf{a} = \mathbf{v}_I$,

\begin{align}\left.\frac{\mathrm d^2\mathbf{r}}{\mathrm dt^2}\right|_I &= \left( \left.\frac{\mathrm d}{\mathrm dt}\right|_R + \mathrm{\omega} \times \right)^2 \mathbf{r} \\&= \left.\frac{\mathrm d^2\mathbf{r}}{\mathrm dt^2}\right|_R + 2\mathbf{\omega} \times \left.\frac{\mathrm d\mathbf{r}}{\mathrm dt}\right|_R + \mathbf{\omega} \times (\mathbf{\omega} \times \mathbf{r}) \;.\end{align}

### Forces on body in rotating frame

Now consider a force acting on an object at position $\mathbf{r}$ (for example, gravity). Newton's third law states

$$\mathbf{F} = m \left.\frac{\mathrm d^2\mathbf{r}}{\mathrm dt^2}\right|_I .$$

And so substituting this into the previous equation for $\left.\frac{\mathrm d^2\mathbf{r}}{\mathrm dt^2}\right|_I$ and rearranging we get

\begin{align}\mathbf{F}_\textrm{net} &= m \left.\frac{\mathrm d^2\mathbf{r}}{\mathrm dt^2}\right|_R \\&= \mathbf{F} - 2m \mathbf{\omega} \times \mathbf{v}_R - m \mathbf{\omega} \times (\mathbf{\omega} \times \mathbf{r})\\ &= \mathbf{F} - 2m \mathbf{\omega} \times \mathbf{v}_R + m \mathbf{\omega}^2 \mathbf{r}.\end{align}

And here we have it. The second term on the right is the Coriolis force, and the third term is the centrifugal force (clearly pointing away from the centre of rotation). Any interpretation of the Coriolis and centrifugal forces then follow naturally from this single important equation.

• +1 nice derivation ;-) You put a plus sign instead of an equals in your formula for $\vec{v}_I = d\vec{r}/dt$, and I think you're missing an exponent of 2 on the operator $\left.d/dt\right|_R+\vec\omega\times$ in the next equation. Nov 9, 2010 at 2:26
• @Noldorin: In the expression for $\vec{a} = a_x \vec{i} + a_y \vec{j} + a_z \vec{k}$ you meant $\vec{a} = a_x \hat{i} + a_y \hat{j} + a_z \hat{k}$, right? Nov 9, 2010 at 6:48
• @Robert: What's the difference ? Nov 9, 2010 at 11:29
• @Robert: Not really. In Britain, we represent the unit vectors for the Cartesian system simply as i, j, and k vectors, no hats. (Or equivalently; x, y, and z vector, with hats.) Nov 9, 2010 at 13:46
• @Noldorin: Interesting, I didn't know the unit vector notation was country-dependent. (I've always preferred the $\hat{x},\hat{y},\hat{z}$ way myself, though... i, j, and k seem like an odd choice of letters) Nov 9, 2010 at 17:05