How Does Angular Momentum Conservation Affect Asteroid Collision Dynamics?

AI Thread Summary
Angular momentum conservation plays a crucial role in asteroid collision dynamics, particularly when a minor asteroid collides with a major one and becomes embedded. Before the collision, both linear and angular momentum are conserved, allowing for the calculation of post-collision velocities and angular velocities. After the collision, the combined system moves as a single rigid body, with the minor asteroid contributing to the overall moment of inertia. The discussion highlights the importance of correctly determining the linear velocity of the center of mass and the angular velocity, especially in the context of centripetal forces acting on the embedded asteroid. Overall, the conversation emphasizes the complexities of applying conservation laws in inelastic collisions involving extended bodies.
  • #251
erobz said:
##\dots## but I believe it was also stated by @kuruman that angular momentum is only conserved about the systems center of mass.
I never stated or implied that. In post #178 I conserved angular momentum about the center of the planet, not about the CM, and showed that you end up with the same equation for ##\omega.##
 
Physics news on Phys.org
  • #252
kuruman said:
I never stated or implied that. In post #178 I conserved angular momentum about the center of the planet, not about the CM, and showed that you end up with the same equation for ##\omega.##
Ah, this is a big problem. I agreed with @erobz, but at this point, I can't really understand why my friend is wrong, mistaking ##d## for ##R##. Why? Sorry to persist with this kind of soap opera....
 
  • #253
kuruman said:
I never stated or implied that. In post #178 I conserved angular momentum about the center of the planet, not about the CM, and showed that you end up with the same equation for ##\omega.##
I don't mean to misquote you...I could have sworn you said that talking about angular momentum not being conserved for a differential mass ##dm## away from the center of mass axis. Maybe I'm imaging it.

EDIT: ( yeah, scrolling though the saga, I don't see it ...very sorry for the misquote)
 
Last edited:
  • #254
kuruman said:
$$\begin{align} & I_{cm}=\frac{2}{5}MR^2+\frac{Mm}{M+m}R^2 \nonumber \\
& I_{cm}=\frac{2}{5}MR^2\left(1+\frac{5}{2}\frac{m}{M+m}\right) \nonumber \\
\end{align}$$When ##m<<M~## you can approximate in the denominator ##M+m\approx M##. Hence $$I_{cm}\approx\frac{2}{5}MR^2\left(1+\frac{5}{2}\frac{m}{M}\right).$$
Ok, thank you, but...
Hak said:
Moreover, you had said that the same result as ##\omega## would be arrived at without preliminary simplifications. By using ##\omega = \frac{L}{I_{cm}}##, with ##I_{cm}## simplified moment of inertia, you arrive at the result ##\omega = \frac{5mMv}{(2M + 5m) R}##, which does not correspond to the correct one. Try redoing the calculations, most likely I am wrong, however, I would like to understand it more. Thank you very much.
What about it, @kuruman?
 
Last edited:
  • #255
kuruman said:
I never stated or implied that. In post #178 I conserved angular momentum about the center of the planet, not about the CM, and showed that you end up with the same equation for ##\omega.##
But where then is that to be applied? If that ##\omega ## was to be applied about ##O## then ##F = m \omega^2 R ## This is the crux of the issue with the @Hak friends solution to the problem.
 
  • #256
erobz said:
But where then is that to be applied? If that ##\omega ## was to be applied about ##O## then ##F = m \omega^2 R ##
This is exactly what I wanted to say... It's really a dilemma...
 
  • #257
erobz said:
But where then is that to be applied? If that ##\omega ## was to be applied about ##O## then ##F = m \omega^2 R ## Becuase this is the crux of the issue with the @Hak friends solution to the problem.
The algebra for finding ##\omega## in post #206 is quite correct and involves no approximation. But as you note, for finding the force, we must decide what radius to use.
Note that whatever force is exerted on m, an equal and opposite force is exerted on M. I think that makes it clear that the correct radius to use is the distance to the common mass centre.
 
  • #258
haruspex said:
The algebra for finding ##\omega## in post #206 is quite correct and involves no approximation. But as you note, for finding the force, we must decide what radius to use.
Note that whatever force is exerted on m, an equal and opposite force is exerted on M. I think that makes it clear that the correct radius to use is the distance to the common mass centre.
Thank you, but what is the connection between Newton's Third Law and the choice of the most appropriate radius to use? Thank you if you would like to answer.
 
  • #259
Hak said:
Thank you, but what is the connection between Newton's Third Law and the choice of the most appropriate radius to use? Thank you if you would like to answer.
I don't know the details of Newtons Third Law argument, but for me it is clear that masses are orbiting their common center of mass, just like we expect of any other composite body freely rotating.
 
  • #260
erobz said:
I don't know the details of Newtons Third Law argument, but for me it is clear that masses are orbiting their common center of mass, just like we expect of any other composite body freely rotating.
Yes, it is correct to me. However, I don't understand what this has to do with the Newton's Third Law argument...
 
  • #261
erobz said:
I don't know the details of Newtons Third Law argument, but for me it is clear that masses are orbiting their common center of mass, just like we expect of any other composite body freely rotating.
I view this as a single rigid body rotating in space with angular velocity ##\vec {\omega}.## It follows that mass element ##dm## at any point P on the body at position vector ##\vec {r}'_p## relative to the CM will have velocity (also relative to the CM) ##\vec v'_p=\vec {\omega}\times \vec {r}'_p##. Now if the CM is moving relative to a point in an inertial frame with velocity ##\vec V_{cm}##, the velocity of P relative to the inertial frame will be $$\vec v_p=\vec V_{cm}+\vec v'_p=\vec V_{cm}+\vec {\omega}\times \vec {r}'_p.$$The centripetal force on mass ##m## is ##\vec F_c=-(dm) \omega^2 \vec{r}'_p.##
 
  • #262
kuruman said:
I view this as a single rigid body rotating in space with angular velocity ##\vec {\omega}.## It follows that mass element ##dm## at any point P on the body at position vector ##\vec {r}'_p## relative to the CM will have velocity (also relative to the CM) ##\vec v'_p=\vec {\omega}\times \vec {r}'_p##. Now if the CM is moving relative to a point in an inertial frame with velocity ##\vec V_{cm}##, the velocity of P relative to the inertial frame will be $$\vec v_p=\vec V_{cm}+\vec v'_p=\vec V_{cm}+\vec {\omega}\times \vec {r}'_p.$$The centripetal force on mass ##m## is ##\vec F_c=-(dm) \omega^2 \vec{r}'_p.##
This would imply that since the minor asteroid is a homogeneous body, it has mass ##dm = m##? And also that ##r'_p = d##?
 
  • #263
Hak said:
This would imply that since the minor asteroid is a homogeneous body, it has mass ##dm = m##? And also that ##r'_p = d##?
No, ##dm## is a very small mass at position ##\vec {r}'## relative to the center of the planet. Imagine chopping up the composite body into many little pieces. ##dm## is the mass of the little piece that you will find at position ##\vec {r}'##; it could be small piece of the asteroid or a small piece of the planet. The only thing you can say is that $$\int dm=M+m.$$
 
  • #264
kuruman said:
No, ##dm## is a very small mass at position ##\vec {r}'## relative to the center of the planet. Imagine chopping up the composite body into many little pieces. ##dm## is the mass of the little piece that you will find at position ##\vec {r}'##; it could be small piece of the asteroid or a small piece of the planet. The only thing you can say is that $$\int dm=M+m.$$
OK, so yours is just a general overview of the physical situation. The result of ##F## you wrote does not lead to the form we found earlier, right?
 
  • #265
Hak said:
OK, so yours is just a general overview of the physical situation. The result of ##F## you wrote does not lead to the form we found earlier, right?
At this point, with 264 posts and counting, I have no clue what form you are referring to. Please post it here so that I can see what you mean.
 
  • #266
kuruman said:
At this point, with 264 posts and counting, I have no clue what form you are referring to. Please post it here so that I can see what you mean.
##F = m \omega^2 d##
 
  • #267
Hak said:
Thank you, but what is the connection between Newton's Third Law and the choice of the most appropriate radius to use? Thank you if you would like to answer.
There's a symmetry here. Each body is pulling the other towards the common mass centre. If that is displaced d from M's centre then ##m(R-d)=Md##, hence##m(R-d)\omega^2=Md\omega^2##.
 
  • #268
haruspex said:
There's a symmetry here. Each body is pulling the other towards the common mass centre. If that is displaced d from M's centre then ##m(R-d)=Md##, hence##m(R-d)\omega^2=Md\omega^2##.
OK, thank you very much. I cannot understand, however, how it can be deduced from this symmetric equation that ##d## is the radius to be chosen for the rotation. Sorry if I don't fully understand...
 
  • #269
haruspex said:
There's a symmetry here. Each body is pulling the other towards the common mass centre. If that is displaced d from M's centre then ##m(R-d)=Md##, hence##m(R-d)\omega^2=Md\omega^2##.
Isn't it simpler to note that any point on the composite rigid object rotates about the CM with constant angular velocity ##\vec{\omega}##? Then that point has centripetal acceleration ##\vec a_c=-\omega^2 \vec {r}'## where ##\vec {r}'## is the position vector relative to the CM. Thus the net force on mass element ##dm## is ##\vec {F}_c=-(dm)\omega^2 \vec {r}'.##
 
  • #270
Hak said:
##F = m \omega^2 d##
Yes, I see now what you are asking. The general expression for the force at any point is ##\vec {F}_c=-(dm)\omega^2 \vec {r}'## where ##\vec {r}'## is the location of the point relative to the CM. The asteroid is at location ##\vec {r}'=\vec d## and has mass ##dm = m##. So the force on the asteroid is $$\vec {F}_c=-m\omega^2 \vec {d}.$$
 
  • #271
kuruman said:
Yes, I see now what you are asking. The general expression for the force at any point is ##\vec {F}_c=-(dm)\omega^2 \vec {r}'## where ##\vec {r}'## is the location of the point relative to the CM. The asteroid is at location ##\vec {r}'=\vec d## and has mass ##dm = m##. So the force on the asteroid is $$\vec {F}_c=-m\omega^2 \vec {d}.$$
Is it correct?
 
  • #272
Hak said:
Ah, this is a big problem. I agreed with @erobz, but at this point, I can't really understand why my friend is wrong, mistaking ##d## for ##R##. Why? Sorry to persist with this kind of soap opera....
It seems to me that acting as an information conduit between your friend and PF is confusing you because you don't seem to have the expertise to diagnose where and why your friend is wrong. Your friend needs to speak for him/herself and post here. Then we will be able to diagnose your friend's problem. Maybe your friend is stringing you along trying to see how many posts this thread will add up to before it is closed down.
 
  • #273
kuruman said:
It seems to me that acting as an information conduit between your friend and PF is confusing you because you don't seem to have the expertise to diagnose where and why your friend is wrong. Your friend needs to speak for him/herself and post here. Then we will be able to diagnose your friend's problem. Maybe your friend is stringing you along trying to see how many posts this thread will add up to before it is closed down.
Yes, you are completely right. Forget my friend's reasoning. One issue, however, remained unresolved: that of the simplified moment of inertia...
 
  • #274
kuruman said:
Isn't it simpler to note that any point on the composite rigid object rotates about the CM with constant angular velocity ##\vec{\omega}##? Then that point has centripetal acceleration ##\vec a_c=-\omega^2 \vec {r}'## where ##\vec {r}'## is the position vector relative to the CM. Thus the net force on mass element ##dm## is ##\vec {F}_c=-(dm)\omega^2 \vec {r}'.##
Yes, I was just providing confirmation from other evidence, since @Hak and @erobz seemed unsure.
 
  • #275
Hak said:
Yes, you are completely right. Forget my friend's reasoning. One issue, however, remained unresolved: that of the simplified moment of inertia...
Please remind us, what issue is this?
 
  • #276
Hak said:
Thank you, that's very comforting to me. Could I know what method you used to arrive at this result? How did you get there through the expansions?

Moreover, you had said that the same result as ##\omega## would be arrived at without preliminary simplifications. By using ##\omega = \frac{L}{I_{cm}}##, with ##I_{cm}## simplified moment of inertia, you arrive at the result ##\omega = \frac{5mMv}{(2M + 5m)(M+m) R}##, which does not correspond to the correct one. Try redoing the calculations, most likely I am wrong, however, I would like to understand it more. Thank you very much.
The second paragraph. Thanks.
 
Last edited:
  • #277
kuruman said:
Funny thing, I did too. My corrected exact value for ##\omega## when ##M=9m## is $$\omega =\frac{1}{5}\frac{v}{R}.$$It turns out that no approximation is needed if one uses the simplified moment of inertia in post#128.

(Original expression edited to add missing factor of 5 in the numerator.)
This is the post I am referring to.
 
  • #278
Hak said:
The second paragraph. Thanks.
Your expression in the second paragraph is dimensionally incorrect. The angular velocity has dimensions of [v/R]. Your expression has dimensions of [Mv/r].
 
  • #279
kuruman said:
Your expression in the second paragraph is dimensionally incorrect. The angular velocity has dimensions of [v/R]. Your expression has dimensions of [Mv/r].
I edited my message, correcting my result. I had mistakenly copied from the paper to the screen. My expression is ##\omega = \frac{5mMv}{(2M + 5m)(M+m) R}##
 
  • #280
but I believe it was also stated by @kuruman that angular momentum is only conserved about the systems center of mass.
kuruman said:
I never stated or implied that. In post #178 I conserved angular momentum about the center of the planet, not about the CM, and showed that you end up with the same equation for ##\omega.##
What was your objection to the statement; that you never said it, and/or that it is not true? Because in the end, we can take angular momentum about any axis ( as you showed previously in post #178), but it is ultimately conserved about the systems center of mass in either case. i.e. ##\omega ## is to be applied about the systems center of mass, not about the arbitrary axis?
 
Last edited:
  • #281
erobz said:
What was your objection to the statement; that you never said it, and/or that it is not true? Because in the end, we can take angular momentum about any axis ( as you showed previously in post #178), but it is ultimately conserved about the systems center of mass in either case. i.e. ##\omega ## is to be applied about the systems center of mass, not about the arbitrary axis?
My objection was to both. I never said it because in this case angular momentum is conserved about any point in space since there are no torques acting on the two-mass system about any point in space. I think I'll write something up to clarify these issues, because the idea is quite simple.
 
  • Like
Likes erobz and hutchphd
  • #282
Consider an irregular rigid object of mass ##M## simultaneously rotating and translating in space in the absence of external forces and torques. This means that both its linear and angular momentum are conserved. We will calculate its angular momentum about an arbitrary point A (see diagram below) given that the velocity of the center of mass is ##\vec V_{cm}## and its rotational angular velocity is ##\vec \omega.##

AngularMomentum.png


First we point out that the axis of rotation must pass through the CM. That's because linear momentum is conserved which means that ##\vec V_{cm}= \rm{const.}## If the axis of rotation did not pass through the CM, the CM would rotate about that hypothetical axis and its velocity would change direction which means that it cannot be constant. Thus, the linear velocity relative to the CM of an arbitrary point at position ##\vec {r}'## from the CM is ##\vec {v}'=\vec{\omega}\times \vec {r}'.##

Now consider mass element ##dm## at an arbitrary point P the position vector of which relative to A is ##\vec r##. Let ##\vec {r}' ## be the position vector of ##dm## relative to the CM and ##\vec R## be the position of the CM relative to A. Finally, let ##\vec v_P## be the instantaneous velocity of point P relative to point A. Addition of velocities requires that the velocity of point P relative to A is ##\vec v_P=\vec V_{cm}+\vec{\omega}\times \vec {r}'.##

The angular momentum contribution of ##dm## about point A is $$d\vec L=(dm)\vec r\times \vec v_P=(dm)(\vec R+\vec {r}')\times(\vec V_{cm}+\vec{\omega}\times \vec {r}')$$ We multiply out the cross product on the RHS and consider each of the four resulting terms separately.

1. ##~(dm) \vec R \times \vec V_{cm}##
The contribution from the first term to the total angular momentum is $$\vec L_1=\int (dm) \vec R \times \vec V_{cm}=\left( \int dm \right) \vec R \times \vec V_{cm}=M\vec R \times \vec V_{cm}.$$ This is often called the angular momentum of the center of mass.

2. ##~(dm)\vec R\times (\vec{\omega}\times \vec {r}')##
Using the triple-product rule,
$$\vec R\times (\vec{\omega}\times \vec {r}')=\vec{\omega}(\vec R\cdot \vec {r}') -\vec {r}'(\vec R \cdot \vec{\omega}).$$The second term on the RHS is zero because ##\vec R## is perpendicular to ##\vec {\omega}##. To find the contribution of the first term, we integrate $$\vec L_2=\int (dm)\vec{\omega}(\vec R\cdot \vec {r}')=\vec{\omega}\left(\vec R\cdot \int (dm)\vec {r}'\right)=0$$ Note that by definition of the CM coordinates, ##\int (dm)\vec {r}'=0.##

3. ##~(dm)\vec {r}'\times \vec V_{cm}##
This term also vanishes upon integration for the same reason when we take constant ##\vec R## out of the integral. $$\vec L_3=\int (dm)\vec {r}'\times \vec V_{cm}=\left(\int (dm)\vec {r}'\right)\times \vec V_{cm}=0.$$ 4. ##~(dm)\vec {r}'\times (\vec{\omega}\times \vec {r}')##
Triple-product rule once more
$$\vec {r}'\times (\vec{\omega}\times \vec {r}')=\vec{\omega}(\vec {r}'\cdot \vec {r}')-\vec {r}'(\vec{\omega}\cdot \vec {r}').$$The second term on the RHS vanishes because the vectors are orthogonal. Integrating the first term, $$\vec L_4=\int \vec{\omega}(\vec {r}'\cdot \vec {r}')=\vec{\omega}\int (dm){r'}^2=I_{cm}~\vec{\omega}.$$This is called the angular momentum about the center of mass. Thus, the total angular momentum about point A of the translating and rotating rigid body is $$\begin{align} & \vec L=\vec L_1+\vec L_4=M\vec R \times \vec V_{cm}+I_{cm}~\vec{\omega}.\end{align}$$ The first term in this expression, angular momentum of the center of mass, depends on the choice of reference point A. The second term, angular momentum about the center of mass, does not.

The bottom line is that one can always write the angular momentum of a translating and rotating rigid body in the form of Equation (1).

Note that a simplification occurs if point A is chosen anywhere on the straight line path of the CM. Then ##\vec {R}## and ##\vec V_{cm}## are parallel and the angular momentum of the CM vanishes. This is why in the asteroid problem and all problems like it, it is convenient to choose the CM as a point of reference for the angular momentum and immediately write the angular momentum conservation equation in the form $$m~v~d=I_{cm}~\omega$$ where ##d## is the vertical distance between the path of the asteroid and the CM of the combined object. This equation says that the initial orbital angular momentum about the CM of the asteroid before the collision is converted to spin angular momentum of the combined object after the collision. Furthermore, as we have seen, choosing point A away from the path of the CM adds equal amounts of (orbital) angular momentum of the CM on both sides of the momentum conservation equation which cancel out. (See post #178.)
 
  • #283
kuruman said:
Consider an irregular rigid object of mass ##M## simultaneously rotating and translating in space in the absence of external forces and torques. This means that both its linear and angular momentum are conserved. We will calculate its angular momentum about an arbitrary point A (see diagram below) given that the velocity of the center of mass is ##\vec V_{cm}## and its rotational angular velocity is ##\vec \omega.##

View attachment 332757

First we point out that the axis of rotation must pass through the CM. That's because linear momentum is conserved which means that ##\vec V_{cm}= \rm{const.}## If the axis of rotation did not pass through the CM, the CM would rotate about that hypothetical axis and its velocity would change direction which means that it cannot be constant. Thus, the linear velocity relative to the CM of an arbitrary point at position ##\vec {r}'## from the CM is ##\vec {v}'=\vec{\omega}\times \vec {r}'.##

Now consider mass element ##dm## at an arbitrary point P the position vector of which relative to A is ##\vec r##. Let ##\vec {r}' ## be the position vector of ##dm## relative to the CM and ##\vec R## be the position of the CM relative to A. Finally, let ##\vec v_P## be the instantaneous velocity of point P relative to point A. Addition of velocities requires that the velocity of point P relative to A is ##\vec v_P=\vec V_{cm}+\vec{\omega}\times \vec {r}'.##

The angular momentum contribution of ##dm## about point A is $$d\vec L=(dm)\vec r\times \vec v_P=(dm)(\vec R+\vec {r}')\times(\vec V_{cm}+\vec{\omega}\times \vec {r}')$$ We multiply out the cross product on the RHS and consider each of the four resulting terms separately.

1. ##~(dm) \vec R \times \vec V_{cm}##
The contribution from the first term to the total angular momentum is $$\vec L_1=\int (dm) \vec R \times \vec V_{cm}=\left( \int dm \right) \vec R \times \vec V_{cm}=M\vec R \times \vec V_{cm}.$$ This is often called the angular momentum of the center of mass.

2. ##~(dm)\vec R\times (\vec{\omega}\times \vec {r}')##
Using the triple-product rule,
$$\vec R\times (\vec{\omega}\times \vec {r}')=\vec{\omega}(\vec R\cdot \vec {r}') -\vec {r}'(\vec R \cdot \vec{\omega}).$$The second term on the RHS is zero because ##\vec R## is perpendicular to ##\vec {\omega}##. To find the contribution of the first term, we integrate $$\vec L_2=\int (dm)\vec{\omega}(\vec R\cdot \vec {r}')=\vec{\omega}\left(\vec R\cdot \int (dm)\vec {r}'\right)=0$$ Note that by definition of the CM coordinates, ##\int (dm)\vec {r}'=0.##

3. ##~(dm)\vec {r}'\times \vec V_{cm}##
This term also vanishes upon integration for the same reason when we take constant ##\vec R## out of the integral. $$\vec L_3=\int (dm)\vec {r}'\times \vec V_{cm}=\left(\int (dm)\vec {r}'\right)\times \vec V_{cm}=0.$$ 4. ##~(dm)\vec {r}'\times (\vec{\omega}\times \vec {r}')##
Triple-product rule once more
$$\vec {r}'\times (\vec{\omega}\times \vec {r}')=\vec{\omega}(\vec {r}'\cdot \vec {r}')-\vec {r}'(\vec{\omega}\cdot \vec {r}').$$The second term on the RHS vanishes because the vectors are orthogonal. Integrating the first term, $$\vec L_4=\int \vec{\omega}(\vec {r}'\cdot \vec {r}')=\vec{\omega}\int (dm){r'}^2=I_{cm}~\vec{\omega}.$$This is called the angular momentum about the center of mass. Thus, the total angular momentum about point A of the translating and rotating rigid body is $$\begin{align} & \vec L=\vec L_1+\vec L_4=M\vec R \times \vec V_{cm}+I_{cm}~\vec{\omega}.\end{align}$$ The first term in this expression, angular momentum of the center of mass, depends on the choice of reference point A. The second term, angular momentum about the center of mass, does not.

The bottom line is that one can always write the angular momentum of a translating and rotating rigid body in the form of Equation (1).

Note that a simplification occurs if point A is chosen anywhere on the straight line path of the CM. Then ##\vec {R}## and ##\vec V_{cm}## are parallel and the angular momentum of the CM vanishes. This is why in the asteroid problem and all problems like it, it is convenient to choose the CM as a point of reference for the angular momentum and immediately write the angular momentum conservation equation in the form $$m~v~d=I_{cm}~\omega$$ where ##d## is the vertical distance between the path of the asteroid and the CM of the combined object. This equation says that the initial orbital angular momentum about the CM of the asteroid before the collision is converted to spin angular momentum of the combined object after the collision. Furthermore, as we have seen, choosing point A away from the path of the CM adds equal amounts of (orbital) angular momentum of the CM on both sides of the momentum conservation equation which cancel out. (See post #178.)
Thank you very much for your interesting contribution!
 
  • #284
kuruman said:
Consider an irregular rigid object of mass ##M## simultaneously rotating and translating in space in the absence of external forces and torques. This means that both its linear and angular momentum are conserved. We will calculate its angular momentum about an arbitrary point A (see diagram below) given that the velocity of the center of mass is ##\vec V_{cm}## and its rotational angular velocity is ##\vec \omega.##

View attachment 332757

First we point out that the axis of rotation must pass through the CM. That's because linear momentum is conserved which means that ##\vec V_{cm}= \rm{const.}## If the axis of rotation did not pass through the CM, the CM would rotate about that hypothetical axis and its velocity would change direction which means that it cannot be constant. Thus, the linear velocity relative to the CM of an arbitrary point at position ##\vec {r}'## from the CM is ##\vec {v}'=\vec{\omega}\times \vec {r}'.##

Now consider mass element ##dm## at an arbitrary point P the position vector of which relative to A is ##\vec r##. Let ##\vec {r}' ## be the position vector of ##dm## relative to the CM and ##\vec R## be the position of the CM relative to A. Finally, let ##\vec v_P## be the instantaneous velocity of point P relative to point A. Addition of velocities requires that the velocity of point P relative to A is ##\vec v_P=\vec V_{cm}+\vec{\omega}\times \vec {r}'.##

The angular momentum contribution of ##dm## about point A is $$d\vec L=(dm)\vec r\times \vec v_P=(dm)(\vec R+\vec {r}')\times(\vec V_{cm}+\vec{\omega}\times \vec {r}')$$ We multiply out the cross product on the RHS and consider each of the four resulting terms separately.

1. ##~(dm) \vec R \times \vec V_{cm}##
The contribution from the first term to the total angular momentum is $$\vec L_1=\int (dm) \vec R \times \vec V_{cm}=\left( \int dm \right) \vec R \times \vec V_{cm}=M\vec R \times \vec V_{cm}.$$ This is often called the angular momentum of the center of mass.

2. ##~(dm)\vec R\times (\vec{\omega}\times \vec {r}')##
Using the triple-product rule,
$$\vec R\times (\vec{\omega}\times \vec {r}')=\vec{\omega}(\vec R\cdot \vec {r}') -\vec {r}'(\vec R \cdot \vec{\omega}).$$The second term on the RHS is zero because ##\vec R## is perpendicular to ##\vec {\omega}##. To find the contribution of the first term, we integrate $$\vec L_2=\int (dm)\vec{\omega}(\vec R\cdot \vec {r}')=\vec{\omega}\left(\vec R\cdot \int (dm)\vec {r}'\right)=0$$ Note that by definition of the CM coordinates, ##\int (dm)\vec {r}'=0.##

3. ##~(dm)\vec {r}'\times \vec V_{cm}##
This term also vanishes upon integration for the same reason when we take constant ##\vec R## out of the integral. $$\vec L_3=\int (dm)\vec {r}'\times \vec V_{cm}=\left(\int (dm)\vec {r}'\right)\times \vec V_{cm}=0.$$ 4. ##~(dm)\vec {r}'\times (\vec{\omega}\times \vec {r}')##
Triple-product rule once more
$$\vec {r}'\times (\vec{\omega}\times \vec {r}')=\vec{\omega}(\vec {r}'\cdot \vec {r}')-\vec {r}'(\vec{\omega}\cdot \vec {r}').$$The second term on the RHS vanishes because the vectors are orthogonal. Integrating the first term, $$\vec L_4=\int \vec{\omega}(\vec {r}'\cdot \vec {r}')=\vec{\omega}\int (dm){r'}^2=I_{cm}~\vec{\omega}.$$This is called the angular momentum about the center of mass. Thus, the total angular momentum about point A of the translating and rotating rigid body is $$\begin{align} & \vec L=\vec L_1+\vec L_4=M\vec R \times \vec V_{cm}+I_{cm}~\vec{\omega}.\end{align}$$ The first term in this expression, angular momentum of the center of mass, depends on the choice of reference point A. The second term, angular momentum about the center of mass, does not.

The bottom line is that one can always write the angular momentum of a translating and rotating rigid body in the form of Equation (1).

Note that a simplification occurs if point A is chosen anywhere on the straight line path of the CM. Then ##\vec {R}## and ##\vec V_{cm}## are parallel and the angular momentum of the CM vanishes. This is why in the asteroid problem and all problems like it, it is convenient to choose the CM as a point of reference for the angular momentum and immediately write the angular momentum conservation equation in the form $$m~v~d=I_{cm}~\omega$$ where ##d## is the vertical distance between the path of the asteroid and the CM of the combined object. This equation says that the initial orbital angular momentum about the CM of the asteroid before the collision is converted to spin angular momentum of the combined object after the collision. Furthermore, as we have seen, choosing point A away from the path of the CM adds equal amounts of (orbital) angular momentum of the CM on both sides of the momentum conservation equation which cancel out. (See post #178.)
@kuruman I got it. My only doubt is with respect to the conservation of linear momentum. In post #178, you say that:
kuruman said:
Note that after the collision, the CM is the only point that moves in a straight line with velocity ##V_{cm}## while all other points rotate about it with angular velocity ##\omega.## This means that the linear momentum of any mass element ##dm## away from the CM is not conserved.
In post #214, @haruspex says:
haruspex said:
It can be written as equations in more than one way, but what the author is saying is that the direction of motion of m is initially unchanged, and M's motion can be thought of as the sum of a rotation about its centre and linear motion of its centre (which is always true).
If M's rotation rate is ##\omega## and initial post-collision mass centre velocity is v' then m's velocity is ##v'+R\omega## and we have ##Mv'+m(v'+R\omega)=mv##.
As can be seen, in @haruspex's message the linear momentum is not conserved with respect to the CoM. Instead, you say that linear momentum is only ever conserved with respect to the CoM. How is this discrepancy possible. Where am I wrong? What wrong assumptions am I making? Thank you very much.
 
  • #285
Hak said:
Thank you, that's very comforting to me. Could I know what method you used to arrive at this result? How did you get there through the expansions?

Moreover, you had said that the same result as ##\omega## would be arrived at without preliminary simplifications. By using ##\omega = \frac{L}{I_{cm}}##, with ##I_{cm}## simplified moment of inertia, you arrive at the result ##\omega = \frac{5mMv}{(2M + 5m)(M+m)R}##, which does not correspond to the correct one. Try redoing the calculations, most likely I am wrong, however, I would like to understand it more. Thank you very much.
However, I have been trying since yesterday, with the simplified moment of inertia calculated by you (even with subsequent corrections) I always arrive at the same result, the one quoted above. This result does not correspond to the original one, whereas you said we would get the same value of ##\omega## with this moment of inertia. I can't see where the error lies and trying still doesn't make sense. Could you tell me which way is correct for what you are saying to start making sense? Thank you very much. @kuruman
 
  • #286
Hak said:
in @haruspex's message the linear momentum is not conserved with respect to the CoM.
I don’t know how you deduce that from what I wrote, and I don't know what you mean by "with respect to the CoM" in the context of linear momentum. Conservation of linear momentum is not in respect of any particular point.
 
  • #287
haruspex said:
I don’t know how you deduce that from what I wrote.
I deduced it from the fact that that equation relates to the conservation of linear momentum with respect to the centre of the major asteroid, not with respect to the centre of mass, as @erobz and I have said many times...
 
  • #288
Hak said:
I deduced it from the fact that that equation relates to the conservation of linear momentum with respect to the centre of the major asteroid, not with respect to the centre of mass, as @erobz and I have said many times...
Please see the text I added to post #286.
 
  • #289
kuruman said:
the linear momentum of any mass element ##dm## away from the CM is not conserved.
So what does this assertion mean? How can it be justified by that equation, which considers ##v'## as the velocity of the centre of the major asteroid, and not the centre of mass? Thank you.
 
  • #290
Hak said:
So what does this assertion mean? How can it be justified by that equation, which considers ##v'## as the velocity of the centre of the major asteroid, and not the centre of mass? Thank you.
@kuruman refers to the momentum conservation of an individual element mass dm. The equations I posted are for momentum conservation of the whole system.
One mass element may gain momentum in a given direction while another loses it, but the total sum is conserved.
 
  • #291
haruspex said:
@kuruman refers to the momentum conservation of an individual element mass dm. The equations I posted are for momentum conservation of the whole system.
One mass element may gain momentum in a given direction while another loses it, but the total sum is conserved.
OK, everything is now clear. Thank you very much.
 
  • #292
Hak said:
This result does not correspond to the original one, whereas you said we would get the same value of ω with this moment of inertia.
Are you saying that when you use the approximate expression $$I_{cm}\approx\frac{2}{5}MR^2\left(1+\frac{5}{2}\frac{m}{M}\right)$$ instead of the exact expression $$I_{cm}=\frac{2}{5}MR^2\left(1+\frac{5}{2}\frac{m}{M+m}\right)$$ in ##\omega = \dfrac{mvd}{I_{cm}}## you get expressions that look different? If so, that is is to be expected because the expressions for ##I_{cm}## are different.

The real question is whether that difference affects ##\omega## significantly. The quickest way to find out is to plug in numbers and see for yourself how close the approximate expression is to the exact expression. Make sure that ##m<<M##, say ##m=10^{-6}M##.

On edit: Corrected typo. See #295.
 
Last edited:
  • #293
kuruman said:
Are you saying that when you use the approximate expression $$I_{cm}\approx\frac{2}{5}MR^2\left(1+\frac{5}{2}\frac{m}{M}\right)$$ instead of the exact expression $$I_{cm}=\frac{2}{5}MR^2\left(1+\frac{5}{2}\frac{m}{M+m}\right)$$ in ##\omega = \dfrac{mv}{I_{cm}}## you get expressions that look different? If so, that is is to be expected because the expressions for ##I_{cm}## are different.

The real question is whether that difference affects ##\omega## significantly. The quickest way to find out is to plug in numbers and see for yourself how close the approximate expression is to the exact expression. Make sure that ##m<<M##, say ##m=10^{-6}M##.
@kuruman Why ##\omega = \frac {mv}{I_{cm}}##? Shouldn't the expression in the numerator be the angular momentum before the collision, thus equal to ##mvd##? I don't understand.
Besides, it was you who said we would get the same expression for ##\omega## with the simplified moment of inertia without approximations, so I didn't understand (and the maths don't add up). What did you mean? I am very uncertain and I don't seem to have understood what you meant (and mean).
 
  • #294
kuruman said:
Funny thing, I did too. My corrected exact value for ##\omega## when ##M=9m## is $$\omega =\frac{1}{5}\frac{v}{R}.$$It turns out that no approximation is needed if one uses the simplified moment of inertia in post#128.

(Original expression edited to add missing factor of 5 in the numerator.)
This is the post I am referring to. What did you mean?
 
  • #295
Hak said:
This is the post I am referring to. What did you mean?
I meant that when I first posted my answer I had ##\omega =\dfrac{1}{25}\dfrac{v}{R}## which I then corrected. It's entirely possible that the correction is wrong but his point is totally irrelevant to the discussion.

Hak said:
@kuruman Why ##\omega = \frac {mv}{I_{cm}}##? Shouldn't the expression in the numerator be the angular momentum before the collision, thus equal to ##mvd##? I don't understand.
Yes, I meant to say ##\omega = \dfrac {mvd}{I_{cm}}.## It was a typo. Thanks for pointing it out.
 
  • #296
Hak said:
Thank you, that's very comforting to me. Could I know what method you used to arrive at this result? How did you get there through the expansions?

Moreover, you had said that the same result as ##\omega## would be arrived at without preliminary simplifications. By using ##\omega = \frac{L}{I_{cm}}##, with ##I_{cm}## simplified moment of inertia, you arrive at the result ##\omega = \frac{5mMv}{(2M + 5m)(M+m)R}##, which does not correspond to the correct one. Try redoing the calculations, most likely I am wrong, however, I would like to understand it more. Thank you very much.
Thanks. I have noticed, however, that by placing ##m \ll M##, one obtains: ##\omega = \frac{5mv}{2MR}##, which corresponds to the value found previously for ##m \ll M##. So, can we say that the simplified moment of inertia is an alternative way to find the simplified angular velocity?
 
  • #297
kuruman said:
Consider an irregular rigid object of mass ##M## simultaneously rotating and translating in space in the absence of external forces and torques. This means that both its linear and angular momentum are conserved. We will calculate its angular momentum about an arbitrary point A (see diagram below) given that the velocity of the center of mass is ##\vec V_{cm}## and its rotational angular velocity is ##\vec \omega.##

View attachment 332757

First we point out that the axis of rotation must pass through the CM. That's because linear momentum is conserved which means that ##\vec V_{cm}= \rm{const.}## If the axis of rotation did not pass through the CM, the CM would rotate about that hypothetical axis and its velocity would change direction which means that it cannot be constant. Thus, the linear velocity relative to the CM of an arbitrary point at position ##\vec {r}'## from the CM is ##\vec {v}'=\vec{\omega}\times \vec {r}'.##

Now consider mass element ##dm## at an arbitrary point P the position vector of which relative to A is ##\vec r##. Let ##\vec {r}' ## be the position vector of ##dm## relative to the CM and ##\vec R## be the position of the CM relative to A. Finally, let ##\vec v_P## be the instantaneous velocity of point P relative to point A. Addition of velocities requires that the velocity of point P relative to A is ##\vec v_P=\vec V_{cm}+\vec{\omega}\times \vec {r}'.##

The angular momentum contribution of ##dm## about point A is $$d\vec L=(dm)\vec r\times \vec v_P=(dm)(\vec R+\vec {r}')\times(\vec V_{cm}+\vec{\omega}\times \vec {r}')$$ We multiply out the cross product on the RHS and consider each of the four resulting terms separately.

1. ##~(dm) \vec R \times \vec V_{cm}##
The contribution from the first term to the total angular momentum is $$\vec L_1=\int (dm) \vec R \times \vec V_{cm}=\left( \int dm \right) \vec R \times \vec V_{cm}=M\vec R \times \vec V_{cm}.$$ This is often called the angular momentum of the center of mass.

2. ##~(dm)\vec R\times (\vec{\omega}\times \vec {r}')##
Using the triple-product rule,
$$\vec R\times (\vec{\omega}\times \vec {r}')=\vec{\omega}(\vec R\cdot \vec {r}') -\vec {r}'(\vec R \cdot \vec{\omega}).$$The second term on the RHS is zero because ##\vec R## is perpendicular to ##\vec {\omega}##. To find the contribution of the first term, we integrate $$\vec L_2=\int (dm)\vec{\omega}(\vec R\cdot \vec {r}')=\vec{\omega}\left(\vec R\cdot \int (dm)\vec {r}'\right)=0$$ Note that by definition of the CM coordinates, ##\int (dm)\vec {r}'=0.##

3. ##~(dm)\vec {r}'\times \vec V_{cm}##
This term also vanishes upon integration for the same reason when we take constant ##\vec R## out of the integral. $$\vec L_3=\int (dm)\vec {r}'\times \vec V_{cm}=\left(\int (dm)\vec {r}'\right)\times \vec V_{cm}=0.$$ 4. ##~(dm)\vec {r}'\times (\vec{\omega}\times \vec {r}')##
Triple-product rule once more
$$\vec {r}'\times (\vec{\omega}\times \vec {r}')=\vec{\omega}(\vec {r}'\cdot \vec {r}')-\vec {r}'(\vec{\omega}\cdot \vec {r}').$$The second term on the RHS vanishes because the vectors are orthogonal. Integrating the first term, $$\vec L_4=\int \vec{\omega}(\vec {r}'\cdot \vec {r}')=\vec{\omega}\int (dm){r'}^2=I_{cm}~\vec{\omega}.$$This is called the angular momentum about the center of mass. Thus, the total angular momentum about point A of the translating and rotating rigid body is $$\begin{align} & \vec L=\vec L_1+\vec L_4=M\vec R \times \vec V_{cm}+I_{cm}~\vec{\omega}.\end{align}$$ The first term in this expression, angular momentum of the center of mass, depends on the choice of reference point A. The second term, angular momentum about the center of mass, does not.

The bottom line is that one can always write the angular momentum of a translating and rotating rigid body in the form of Equation (1).

Note that a simplification occurs if point A is chosen anywhere on the straight line path of the CM. Then ##\vec {R}## and ##\vec V_{cm}## are parallel and the angular momentum of the CM vanishes. This is why in the asteroid problem and all problems like it, it is convenient to choose the CM as a point of reference for the angular momentum and immediately write the angular momentum conservation equation in the form $$m~v~d=I_{cm}~\omega$$ where ##d## is the vertical distance between the path of the asteroid and the CM of the combined object. This equation says that the initial orbital angular momentum about the CM of the asteroid before the collision is converted to spin angular momentum of the combined object after the collision. Furthermore, as we have seen, choosing point A away from the path of the CM adds equal amounts of (orbital) angular momentum of the CM on both sides of the momentum conservation equation which cancel out. (See post #178.)
@kuruman Also, I would like to ask you how you drew the diagram above, and how you obtained the depiction of the irregular rigid body. Thank you very much, if you feel like telling me.
 
  • #298
Hak said:
@kuruman Also, I would like to ask you how you drew the diagram above, and how you obtained the depiction of the irregular rigid body. Thank you very much, if you feel like telling me.
Like @erobz I use Powerpoint to make figures. There is a "Freeform" feature to draw closed 2d shapes freehand. In other versions of Powerpoint the icon shown below might look different.

Freeform.png
 
  • #299
kuruman said:
Like @erobz I use Powerpoint to make figures. There is a "Freeform" feature to draw closed 2d shapes freehand. In other versions of Powerpoint the icon shown below might look different.

View attachment 332784
You can also use the curve function ( which I believe is just to the left ). It will make smooth/continuous curve, and closes automatically in the vicinity of the start point.
 
  • #300
erobz said:
You can also use the curve function ( which I believe is just to the left ). It will make smooth/continuous curve, and closes automatically in the vicinity of the start point.
Actually both functions on either side of the one with the arrow draw closed surfaces.
 
Back
Top