Friday 10 July 2020

homework and exercises - Is this expression for the kinetic energy of a spinning disk revolving about a second axis correct?


My question is motivated from a question from another user. You can see the configuration of the rotating system here: https://physics.stackexchange.com/q/143377/. enter image description here



I am not interested in all the complicated arguments of his question, but only of the expression for the total kinetic energy. My answer was that the rotational KE can be expressed as the addition of the KE of the center of mass plus the KE relative to the center of mass, which results in this expression:


$$E_k=\frac{1}{4}mr^2\omega_2^2+\frac{1}{2}md^2\omega_1^2$$


(Notice that this result is independent of the sign of $\omega_2$).


But the original OP claims that the right expression is $$E_k=\frac{1}{2}md^2\omega_1^2+\frac{1}{2}mr^2(\omega_1-\omega_2)^2,$$ based on answers obtained on other forums (which I checked) and even the moderators in those forums seem to agree with it. The OP itself does not know enough physics to come up with its own answer, but still does not believe mine for the reasons given above.


So, my question is: I am missing something pretty obvious here? which of the expressions is the correct one (if any?) Thanks!



Answer



The answer depends on what the symbols mean. The question does not make it clear how the symbols are defined. The most confusing quantity is $\omega_2$. How is this defined? Is it the angular velocity of the disc relative to the fixed lab axes or relative to the axle about which it is rotating (where this axle itself will be rotating at $\omega_1$)? Also what is the sign convention for $\omega_1$? The problem states that $\omega_1$ is rotating clockwise, so is a positive $\omega_1$ supposed to mean clockwise rotation, or counterclockwise rotation? We will see that different answers to these questions give different expressions for the kinetic energy---one gives your answer and the other gives his answer. Thus I think the ultimate reason for disagreement is confusion about what the symbols mean.


Let us first solve the problem using one choice of meaning for the symbols, and obtain the expression for the kinetic energy, then we will see how the expression changes when we use different meanings for the symbols. I will use "your" definition, where positive $\omega_1$ and $\omega_2$ both indicate rotations in a counterclockwise direction, and the meaning of $\omega_2$ is the angular velocity of the disc about its center of mass relative to the fixed lab frame.


Then we know the center of mass of the disc has a speed $\omega_1 d$, and the disc rotates at an angular velocity $\omega_2$ with respect to the lab axes. Then we can obtain the kinetic energy using the following result: the kinetic energy of a rigid object is the sum of a translational piece $T_{\textrm{trans}}$ given by $T_{\textrm{trans}}=\dfrac{1}{2} m v_{cm}^2$ and a rotation piece $T_{\textrm{rot}}$ given by $T_{\textrm{rot}}=\dfrac{1}{2}I \omega^2$ where $m$ is the mass of the object, $v_{cm}$ is its center of mass speed, $\omega$ is its angular velocity about its center of mass, and $I$ is its moment of inertia about its center of mass (this is shown in the appendix). Using $v_{cm} = \omega_1 d$, $\omega = \omega_2$, and $I=\frac{1}{2}mr^2$, we find the total kinetic energy $T$ is given by $T=\dfrac{1}{2}m \omega_1^2 d^2 + \dfrac{1}{4} m \omega_2^2 r^2$. This is the expression you have.


Now let's see what happens if we instead express the kinetic energy not in terms of the absolute disc angular velocity $\omega_2$, but the angular velocity $\omega'_2$ of the disc relative to the axle. Then the angular velocity of the disc relative to the lab frame is the angular velocity of the disc relative to the axle plus the angular velocity of the axle relative to the disc. That is, $\omega_2 = \omega'_2 + \omega_1$. Plugging this in to our equation for kinetic energy, we get $T=\dfrac{1}{2}m \omega_1^2 d^2 + \dfrac{1}{4} m (\omega'_2 + \omega_1)^2 r^2$. Now suppose we want to express this kinetic energy in terms of $\omega'_1$, where $\omega'_1$ is positive for clockwise rotations and negative for counterclockwise rotations; that is, $\omega'_1= -\omega_1$. Then the expression for the kinetic energy becomes $T=\dfrac{1}{2}m \omega_1^2 d^2 + \dfrac{1}{4} m (\omega'_2-\omega'_1)^2 r^2$. This is now his expression. So it seems the only difference was what you meant by your variables




$\newcommand{\r}{\mathbf{r}}$ $\newcommand{\v}{\mathbf{v}}$ $\newcommand{\w}{\boldsymbol\omega}$ Here I will explain the only physics in the problem, which is that the kinetic energy of a rigid body can be decomposed into a part that gives the translation energy of the center of mass and a part for rotation about the center of mass. At some point in time this object is described by a spatial density profile $\rho(\r)$, and a spatial velocity profile $\mathbf{v}(\r)$. Since the body is rigid, the velocity profile $\mathbf{v}(\r)$ must have the form $\mathbf{v}(\r) = \mathbf{v}_0 + \w \times \r$, where $\w$ is the angular velocity. The kinetic energy is given by $T=\dfrac{1}{2}\int \rho(\r) v^2(\r) d\r$.


To see how to break this up into center of mass quantities, we will need to define the center of mass quantites. Let us define the total mass $m$ of the object to be $m=\int \rho(\r)d\r$, and let us define the center of mass position $\r_{cm}$ by $\r_{cm} = \dfrac{1}{m} \int \r \rho(\r) d \r$. Let us define the center of mass velocity $\v_{cm}$ to be the derivate of the center of mass position, which is given by the expression $\v_{cm} = \dfrac{1}{m} \int \v(\r) \rho(\r) d \r$.


First we will prove an intuitive result, which is that $\v(\r_{cm}) = \v_{cm}$. This is intuitive because the body is rigid so the center of mass has to move with the object. To prove this just see that $\begin{equation} \begin{aligned} \v_{cm} &= \dfrac{1}{m} \int \v(\r) \rho(\r) d \r \\ &= \dfrac{1}{m} \int (\mathbf{v}_0 + \w \times \r) \rho(\r) d \r \\ &= \dfrac{1}{m} \int \mathbf{v}_0 \rho(\r) d \r + \dfrac{1}{m} \int \w \times \r\rho(\r) d \r \\ &= \mathbf{v}_0 \dfrac{1}{m} \int \rho(\r) d \r + \w \times \dfrac{1}{m} \int \r\rho(\r) d \r \\ &= \mathbf{v}_0 + \w \times \r_{cm} \\ &= \v(\r_{cm}) \end{aligned} \end{equation}$


Now the first we want to do is re-express $\v$ in terms of $\r_{cm}$ and $\v_{cm}$. We have that $\begin{equation} \begin{aligned} \mathbf{v}(\r) &= \mathbf{v}_0 + \w \times \r \\ &= \mathbf{v}_0 + \w \times \r - \v(\r_{cm}) + \v(\r_{cm}) \\ &= \mathbf{v}_0 + \w \times \r - (\mathbf{v}_0 + \w \times \r_{cm}) + \v(\r_{cm}) \\ &= \w \times (\r - \r_{cm}) + \v(\r_{cm}) \\ &= \w \times (\r - \r_{cm}) + \v_{cm} \end{aligned} \end{equation}$


Now let us plug this form of $\v(\r)$ into our formula for the kinetic energy $T$: $\begin{equation} \begin{aligned} 2T&= \int \rho(\r) v^2(\r) d\r \\ &= \int \rho(\r) (\w \times (\r - \r_{cm}) + \v_{cm})^2 d\r \\ &= \int \rho(\r) ((\w \times (\r - \r_{cm}))^2 + 2 (\w \times (\r - \r_{cm}))\cdot\v_{cm} + v^2_{cm}) d\r \\ &= \int \rho(\r) (\w \times (\r - \r_{cm}))^2 d\r + \int \rho(\r) 2 (\w \times (\r - \r_{cm}))\cdot\v_{cm} d\r + \int \rho(\r) v^2_{cm} d\r \end{aligned} \end{equation}$


Let's look at this last line term by term. Let's start with the first term. $\begin{equation} \begin{aligned} \int \rho(\r) (\w \times (\r - \r_{cm}))^2 d\r &= \omega^2 \int \rho(\r) (\hat{\w} \times (\r - \r_{cm}))^2 d\r \end{aligned} \end{equation}$ Now the integral on the right had side is the weighted sum of the square of the component perpendicular to $\w$ of displacement from the center of mass. That is, it is the scalar moment of inertia for rotations about the axis parallel to $\w$ passing through the center of mass. Denoting this moment of inertia $I$, we find the first term is $I \omega^2$.


Moving on to the second term for the kinetic energy, we see $\begin{equation} \begin{aligned} \int \rho(\r) 2 (\w \times (\r - \r_{cm}))\cdot\v_{cm} d\r &= 2 \v_{cm} \cdot \left( \w \times \int \rho(\r) (\r - \r_{cm}) d\r \right) \\ &= 2 \v_{cm} \cdot \left(\w \times \left( \int \rho(\r) \r d\r - \int \rho(\r)\r_{cm} d\r \right)\right) \\ &=2 \v_{cm} \cdot \left(\w \times (m \r_{cm} - m \r_{cm}) \right) \\ &=0 \end{aligned} \end{equation}$


Finally let's look at the third term. The third term is


$\begin{equation} \begin{aligned} \int \rho(\r) v^2_{cm} d\r &= v^2_{cm} \int \rho(\r) d\r \\ &= m v^2_{cm} \end{aligned} \end{equation}$



Putting the three terms back together, we find that $2T = I \omega^2 + m v_{cm}^2$, or $T=\dfrac{1}{2}I \omega^2 +\dfrac{1}{2} m v_{cm}^2$ where $I$ is the moment of inertia about the center of mass. This is what we needed to show.


No comments:

Post a Comment

Understanding Stagnation point in pitot fluid

What is stagnation point in fluid mechanics. At the open end of the pitot tube the velocity of the fluid becomes zero.But that should result...