Help with derivation - involves cross & dot products

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
4 replies · 3K views
lets_resonate
Messages
15
Reaction score
0
Hello,

I'm reading a classical mechanics book. In it, they show a derivation of the centrifugal force equation:

[tex]\vec{ F_{\textrm{cf}}} = -m \vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right)[/tex]

I understood the derivation up to that point. However, they have a couple additional steps after that whereby they derive another equation based on this one. It is valid if the axis of rotation is chosen to lie along the z-axis:

[tex]\vec{ F_{\textrm{cf}}} = m \| \vec{ \omega } \| ^ 2 \vec \rho[/tex]

Here, [itex]\vec \rho[/itex] is the "cylindrical-radius vector to the particle from the z-axis".

The derivation only shows a couple steps:

[tex]\vec{ F_{\textrm{cf}}} = -m \vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right)[/tex]
[tex]\vec{ F_{\textrm{cf}}} = -m \left[ \vec{ \omega } \left( \vec{ \omega } \cdot \vec{ r } \right) - \vec{r} \omega^{2} \right][/tex]
[tex]\vec{ F_{\textrm{cf}}} = m \| \vec{ \omega } \| ^ 2 \left( x \mathbf{\hat{x}} + y \mathbf{\hat{y}} \right)[/tex]
[tex]\vec{ F_{\textrm{cf}}} = m \| \vec{ \omega } \| ^ 2 \vec \rho[/tex]

I'm confused about how they went from the first step to the second step. Can anyone please help?

Attempt at a solution

If the first two steps in the derivation are valid, it would imply this:

[tex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right) = \vec{ \omega } \left( \vec{ \omega } \cdot \vec{ r } \right) - \vec{r} \| \vec \omega \|^{2}[/tex]

So I decided to work with the left side of this equation and hope that it would yield the right side. I used this definition of the cross product:

[tex]\vec{ a } \times \vec{ b } = \left< a_y b_z - a_z b_y, a_z b_x - a_x b_z, a_x b_y - a_y b_x \right>[/tex]

Then I wrote out [itex]\vec{\omega}[/itex] in component form. Since it is chosen to lie along the z-axis, its x and y components are zero:

[tex]\vec{ \omega } = \left< 0, 0, \| \vec{ \omega } \| \right>[/tex]

Then I wrote out [itex]\vec r[/itex]. I simply named its components x, y, and z:

[tex]\vec{ r } = \left< x, y, z \right>[/tex]

Then [itex]\vec{ \omega } \times \vec{ r }[/itex] is:

[tex]\vec{ \omega } \times \vec{ r } = \left< -\| \vec{ \omega } \| y, \| \vec{ \omega } \| x, 0 \right>[/tex]

And [itex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right)[/itex] is:

[tex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right) = \left< -\| \vec{ \omega } \|^2 x, -\| \vec{ \omega } \|^2 y, 0 \right>[/tex]

[tex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right) = - \vec{ r } \| \vec{ \omega } \|^2[/tex]

However, this is only the second term of the equation above (the first one under my attempt). The first term is missing. What did I do wrong?
 
Last edited:
Physics news on Phys.org
lets_resonate said:
Hello,

I'm reading a classical mechanics book. In it, they show a derivation of the centrifugal force equation:

[tex]\vec{ F_{\textrm{cf}}} = -m \vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right)[/tex]

I understood the derivation up to that point. However, they have a couple additional steps after that whereby they derive another equation based on this one. It is valid if the axis of rotation is chosen to lie along the z-axis:

[tex]\vec{ F_{\textrm{cf}}} = m \| \vec{ \omega } \| ^ 2 \vec \rho[/tex]

Here, [itex]\vec \rho[/itex] is the "cylindrical-radius vector to the particle from the z-axis".

The derivation only shows a couple steps:

[tex]\vec{ F_{\textrm{cf}}} = -m \vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right)[/tex]
[tex]\vec{ F_{\textrm{cf}}} = -m \left[ \vec{ \omega } \left( \vec{ \omega } \cdot \vec{ r } \right) - \vec{r} \omega^{2} \right][/tex]
[tex]\vec{ F_{\textrm{cf}}} = m \| \vec{ \omega } \| ^ 2 \left( x \mathbf{\hat{x}} + y \mathbf{\hat{y}} \right)[/tex]
[tex]\vec{ F_{\textrm{cf}}} = m \| \vec{ \omega } \| ^ 2 \vec \rho[/tex]

I'm confused about how they went from the first step to the second step. Can anyone please help?

Attempt at a solution

If the first two steps in the derivation are valid, it would imply this:

[tex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right) = \vec{ \omega } \left( \vec{ \omega } \cdot \vec{ r } \right) - \vec{r} \| \vec \omega \|^{2}[/tex]
Yes, this is a vector identity (it's easy to prove using the Levi-Civita symbol if you have seen that)
So I decided to work with the left side of this equation and hope that it would yield the right side. I used this definition of the cross product:

[tex]\vec{ a } \times \vec{ b } = \left< a_y b_z - a_z b_y, a_z b_x - a_x b_z, a_x b_y - a_y b_x \right>[/tex]

Then I wrote out [itex]\vec{\omega}[/itex] in component form. Since it is chosen to lie along the z-axis, its x and y components are zero:

[tex]\vec{ \omega } = \left< 0, 0, \| \vec{ \omega } \| \right>[/tex]

Then I wrote out [itex]\vec r[/itex]. I simply named its components x, y, and z:

[tex]\vec{ r } = \left< x, y, z \right>[/tex]

Then [itex]\vec{ \omega } \times \vec{ r }[/itex] is:

[tex]\vec{ \omega } \times \vec{ r } = \left< -\| \vec{ \omega } \| y, \| \vec{ \omega } \| x, 0 \right>[/tex]

And [itex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right)[/itex] is:

[tex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right) = \left< -\| \vec{ \omega } \|^2 x, -\| \vec{ \omega } \|^2 y, 0 \right>[/tex]

[tex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right) = - \vec{ r } \| \vec{ \omega } \|^2[/tex]

However, this is only the second term of the equation above (the first one under my attempt). The first term is missing. What did I do wrong?
It's your very last step that you did incorrectly. The line before is NOT equal to [itex]- \vec{ r } \| \vec{ \omega } \|^2[/itex]
Do you see why?:wink:
 
Gah! I got it! But I might have one more small question:

[tex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right) = \left< -\| \vec{ \omega } \|^2 x, -\| \vec{ \omega } \|^2 y, 0 \right>[/tex]
[tex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right) = -\| \vec{ \omega } \|^2 \left< x, y, 0 \right>[/tex]
[tex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right) = -\| \vec{ \omega } \|^2 \left( \vec{r} - z \right)[/tex]
[tex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right) = \| \vec{ \omega } \|^2 z - \vec{r} \| \vec{ \omega } \|^2[/tex]
[tex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right) = \| \vec{ \omega } \| \left( \| \vec{ \omega } \| z \right) - \vec{r} \| \vec{ \omega } \|^2[/tex]
[tex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right) = \| \vec{ \omega } \| \left( \vec \omega \cdot \vec r \right) - \vec{r} \| \vec{ \omega } \|^2[/tex]

In the very last line: I have the magnitude of [itex]\vec \omega[/itex] outside the parantheses in the first term, but in the book they have the actual vector. Did I miss a small detail somewhere that led to this mistake?

====

And also, a bigger question. This classical mechanics book that I'm reading -- it seems to move along at a pretty vigorous pace. In many cases, there is very little explanation that accompanies the derivations.

Now, I'm only a pithy little high school student. I took first year physics and calculus. What other prerequisites are there for a typical classical mechanics course? Where did you guys learn your "Levi-Civita symbols" (whatever those are)? Or am I actually ready for this with my background?
 
lets_resonate said:
Gah! I got it! But I might have one more small question:

[tex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right) = \left< -\| \vec{ \omega } \|^2 x, -\| \vec{ \omega } \|^2 y, 0 \right>[/tex]
[tex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right) = -\| \vec{ \omega } \|^2 \left< x, y, 0 \right>[/tex]
[tex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right) = -\| \vec{ \omega } \|^2 \left( \vec{r} - z \right)[/tex]
It's not quite right to write this. I know what you mean, but you are not subtracting z (which is a scalar) but you are subtracting the *vector* [itex]z \vec{k}[/itex]

What you must then use is that z by itself is

[tex]z = \frac{1}{\|\vec{\omega} \|} \vec{\omega} \cdot \vec{r}[/tex]

and the unit vector [itex]\vec{k}[/itex] can be written as [itex]\frac{\vec{\omega}}{\| \vec{\omega} \|}[/itex]

Then you will get the desired result.


But notice that this is not a very satisfying approach since a lots of steps are involved and it's for a very special case (omega along z). It's much much better to do it for the general case with the general identity. But I understand the need to do a specific case the long way.

[tex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right) = \| \vec{ \omega } \|^2 z - \vec{r} \| \vec{ \omega } \|^2[/tex]
[tex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right) = \| \vec{ \omega } \| \left( \| \vec{ \omega } \| z \right) - \vec{r} \| \vec{ \omega } \|^2[/tex]
[tex]\vec{ \omega } \times \left( \vec{ \omega } \times \vec{ r } \right) = \| \vec{ \omega } \| \left( \vec \omega \cdot \vec r \right) - \vec{r} \| \vec{ \omega } \|^2[/tex]

In the very last line: I have the magnitude of [itex]\vec \omega[/itex] outside the parantheses in the first term, but in the book they have the actual vector. Did I miss a small detail somewhere that led to this mistake?

====

And also, a bigger question. This classical mechanics book that I'm reading -- it seems to move along at a pretty vigorous pace. In many cases, there is very little explanation that accompanies the derivations.

Now, I'm only a pithy little high school student. I took first year physics and calculus. What other prerequisites are there for a typical classical mechanics course? Where did you guys learn your "Levi-Civita symbols" (whatever those are)? Or am I actually ready for this with my background?

What book are you using? Is it for self-study or are you taking a class?

The key point is that in order to do more advanced mechancis, you need to build some background in maths, mostly in vector analysis and vector calculus. So At this point it might be good to pick up a math book (oriented toward physicists or engineers. I would suggest Boas or Arfken) and to just focus on these topics before going back to Mechanics.
 
Last edited: