## Section10.2Angular Momentum

Consider a point particle of mass $m$ at position $\vec r$ from origin $O$ is moving with speed $\vec v\text{.}$ Then, angular momentum with respect to O is given by

\begin{equation} \vec L = \vec r \times \vec p = m \vec r \times \vec v.\tag{10.2.1} \end{equation}

We can write this also using angular velocity $\vec \omega\text{.}$ Let us suppose that this coordinate system is the fixed coordinate system $O'x'y'z'$ we have discussed in the last section. When discussing angular velocity we found that

\begin{equation*} \vec v = \vec \omega \times \vec r. \end{equation*}

Therefore, angular momentum can also be written as

\begin{equation} \vec L = m \vec r \times (\vec \omega \times \vec r).\label{eq-ang-mom-vector-single-particle}\tag{10.2.2} \end{equation}

When a rigid body rotates, all of its particles rotate at the same angular velocity. But since different particles will be in different positions, their $\vec r$ will be different. Say, we represent a rigid body by a collection of point particles of masses $m_i\text{,}$ $i = 1, 2, \cdots, N\text{.}$ Then, angular momentum of the body about O will be

\begin{equation} \vec L = \sum_{i=1}^N m_i \vec r_i \times (\vec \omega \times \vec r_i).\label{eq-ang-mom-vector-multi-particle}\tag{10.2.3} \end{equation}

Expanding this in components can be done by writing $\vec r$ and $\vec \omega$ into components and carrying out the cross product.

\begin{equation*} \vec r_i = x_i\,\hat i + y_i\,\hat j + z_i\,\hat k,\ \ \ \vec \omega = \omega_x\,\hat i + \omega_y\,\hat j + \omega_z\,\hat k. \end{equation*}

Try to do this algebra to find the following expression for $x$-component $L_x\text{.}$

\begin{align*} L_x \amp = \sum_i m_i \left( y_i^2 + z_i^2\right)\omega_x - \sum_i m_i x_i y_i \omega_y - \sum_i m_i x_i z_i \omega_z. \end{align*}

The other components of angular momentum is easy to write down by analogy.

\begin{align*} L_y\amp = \sum_i m_i \left( x_i^2 + z_i^2\right)\omega_y - \sum_i m_i y_i x_i \omega_x - \sum_i m_i y_i z_i \omega_z,\\ L_z\amp = \sum_i m_i \left( x_i^2 + y_i^2\right)\omega_z - \sum_i m_i z_i x_i \omega_x - \sum_i m_i z_i y_i \omega_y. \end{align*}

Often it is helpful to drop the suffix that denumerates particles and just use the summation sign.

\begin{align} \amp L_x = \sum m \left( y^2 + z^2\right)\omega_x - \sum m x y \omega_y - \sum m x z \omega_z.\tag{10.2.4}\\ \amp L_y = \sum m \left( x^2 + z^2\right)\omega_y - \sum m y x \omega_x - \sum m y z \omega_z.\tag{10.2.5}\\ \amp L_z = \sum m \left( x^2 + y^2\right)\omega_z - \sum m z x \omega_x - \sum m z y \omega_y.\tag{10.2.6} \end{align}

Often, we write relation between components of $\vec L$ and $\vec \omega$ using a matrix multiplication notation.

\begin{equation} \begin{pmatrix} L_x\\ L_y\\ L_z \end{pmatrix} = \begin{pmatrix} \sum m \left( y^2 + z^2\right) \amp - \sum m x y \amp - \sum m x z\\ - \sum m y x \amp \sum m \left( x^2 + z^2\right) \amp - \sum m y z\\ - \sum m z x \amp - \sum m z y \amp \sum m \left( x^2 + y^2\right) \end{pmatrix} \begin{pmatrix} \omega_x\\ \omega_y\\ \omega_z \end{pmatrix} \label{eq-vec-L-I-vec-omega-tensor-notation}\tag{10.2.7} \end{equation}

We write this equation in a formal way by denoting the quantities in the matrix by bold letter ${\bf I}\text{.}$

\begin{equation} \vec L = {\bf I}\, \vec \omega.\label{eq-vec-L-I-vec-omega}\tag{10.2.8} \end{equation}

Unlike vectors $\vec L$ and $\vec \omega\text{,}$ which have three Cartesian components, ${\bf I}$ has nine components. We call ${\bf I}$ moment inertial tensor. There are other tensor type properties in physics such as electic conductivity and refractive index tensors. The components of tensor are denoted by two subscripts corresponding to the axes. Thus,

\begin{align} \amp I_{xx} = \sum m \left( y^2 + z^2\right),\ \ I_{xy} = - \sum m x y,\ \ I_{xz} = - \sum m x z\tag{10.2.9}\\ \amp I_{yx} = - \sum m y x,\ \ I_{yy} = \sum m \left( x^2 + z^2\right) ,\ \ I_{yz} = - \sum m y z\tag{10.2.10}\\ \amp I_{zx} = - \sum m z x,\ \ I_{zy} = - \sum m z y,\ \ I_{zz} = \sum m \left( x^2 + y^2\right) \tag{10.2.11} \end{align}

The diagonal components, $I_{xx}\text{,}$ $I_{yy}\text{,}$ $I_{zz}$ are called moments of inertia about $x\text{,}$ $y\text{,}$ and $z$ axes respectively. The off-diaginal components do not have such simple interpretations.

If a body is regarded as continuous, the sum in these definitions become integrals. To go from sum over masses to integrals over space, we replace $m$ by $\rho dV\text{,}$ where $\rho$ is the density, $dV$ is a volume element. In the case of two-dimensional structure, you will have area element $dA$ and in the case of a one-dimensional structure it will be just a size $ds\text{.}$ Thus, the diagonal elment $I_{zz}$ will be

\begin{equation} I_{zz} = \int\, \rho \, \left( x^2 + y^2\right)\, dV,\tag{10.2.12} \end{equation}

and off-diagonal element $I_{zx}$ will be

\begin{equation} I_{zx} = -\int\, \rho \, z x \, dV.\tag{10.2.13} \end{equation}

The other components of ${\bf I}$ will be obtained in a similar way.

From the matrix representation of ${\bf I}\text{,}$ it is obvious that it is a symmetric matrix. We say that moment of inertia is a symmetric tensor. Since symmetric matrices can always be diagonalized, there is a choice of direction of the body-attached axes so that off-diagonal elements turn out to be zero. We call these axes principal axes. In principal axes we will have only three non-zero moment of inertia components; in that case, we just denote them by simpler symbols, $I_1\text{,}$ $I_2\text{,}$ $I_3\text{.}$

\begin{equation*} {\bf I} = = \begin{pmatrix} I_1 \amp 0 \amp 0\\ 0 \amp I_2 \amp 0\\ 0 \amp 0 \amp I_3 \end{pmatrix} \end{equation*}

Figure 10.2.4 shows a dumbbell of two point masses $m$ and length $2a$ placed on $z$-axis and rotating about the axis in the direction of unit vector $\hat u\text{.}$ Ignore the mass of the connecting rod. Find the angular momentum of the dumbbell.

Hint

First find moments of inertia using coordinates of the two masses.

$\vec L = 2 m a^2 \omega_x \hat i + 2 m a^2 \omega_y \hat j \text{.}$

Solution

The coordinates of the two point masses of the rigid body are

\begin{align*} \amp x_1 = 0,\ y_1 = 0,\ z_1 = a\\ \amp x_2 = 0,\ y_2 = 0,\ z_2 = -a \end{align*}

Now, we compute the six moments of inertia for this system.

\begin{align*} \amp I_{xx} = ma^2 + ma^2 = 2 m a^2 = I_{yy} \\ \amp I_{zz} = 0 \\ \amp I_{xy} = 0 = I_{xz} = I_{yz} \end{align*}

Let us write angular velocity vector $\hat \omega$ as

\begin{equation*} \hat \omega = \omega_x\, \hat i + \omega_y\, \hat j + \omega_z\, \hat k. \end{equation*}

Therefore, Cartesian components of angular momentum will be

\begin{align*} \amp L_x = 2 m a^2 \omega_x \\ \amp L_y = 2 m a^2 \omega_y \\ \amp L_z = 0 \end{align*}

Figure 10.2.6 shows a dumbbell of two point masses $m$ and length $2a$ that rotates about an axis at angular speed $\omega$ that is not a principal axis of the dumbbell. Pointing $z$-axis towards the axis of rotation leaves the dumbell rotating at an angle $\theta$ with respect to the axis.

Supposing that the masses are in $xz$-plane at time $t=0$ so that azimuth of one of the masses is $\phi=\omega t\text{.}$ (a) Find angular momentum vector at an arbitrary instant. (b) Show that the magnitude of angular momentum is constant in time. Ignore the mass of the connecting rod.

Hint

First find moments of inertia using coordinates of the two masses. Pick a moment to get $x$ and $y$ coordinates in terms of spherical angles.

(b) $2m\omega a^2\sin\theta\text{.}$

Solution 1 (a)

(a) Let's work out coordinates of the two masses at an instant when the azimuth angle for one of the masses is $\phi=\omega t$ as shown in Figure 10.2.7. From this figure, we can work out following coordinates of the two masses.

\begin{align*} \amp x_1 = a\,\sin\theta\cos\phi,\ y_1 = a\,\sin\theta\sin\phi,\ z_1 = -a\cos\theta\\ \amp x_2 = -x_1,\ y_2 = -y_1,\ z_2 = -z_1 \end{align*}

Now, we compute the six moments of inertia for this system.

\begin{align*} \amp I_{xx} = 2ma^2(\sin^2\theta \sin^2\phi + \cos^2\theta) \\ \amp I_{yy} = 2ma^2(\sin^2\theta \cos^2\phi + \cos^2\theta) \\ \amp I_{zz} = 2ma^2(\sin^2\theta \cos^2\phi + \sin^2\theta \sin^2\phi ) = 2ma^2\sin^2\theta \\ \amp I_{xy} = -2ma^2\sin^2\theta \sin\phi\cos\phi\\ \amp I_{xz} = 2ma^2\sin\theta\cos\theta \cos\phi \\ \amp I_{yz} = 2ma^2\sin\theta\cos\theta \sin\phi \end{align*}

Since angular velocity is pointed towards $z$-axis, we have

\begin{equation*} \omega_x=0,\ \ \omega_y=0,\ \ \omega_z = \omega. \end{equation*}

Therefore, after using $\phi=\omega t\text{,}$ Cartesian components of angular momentum will be

\begin{align*} \amp L_x = I_{xz} \omega_z = 2m\omega a^2\sin\theta\cos\theta \cos\omega t\\ \amp L_y = I_{yz} \omega_z = 2m\omega a^2\sin\theta\cos\theta \sin\omega t\\ \amp L_z = I_{zz}\omega_z = 2ma^2\sin^2\theta \end{align*}
Solution 2 (b)

(b) From the components of angular momentum vector we get the magnitude by

\begin{equation*} L = \sqrt{L_x^2 + L_y^2 + L_z^2}. \end{equation*}

Let's work out each square and then add them.

\begin{align*} \amp L_x^2 = (2m\omega a^2)^2 \sin^2\theta\cos^2\theta \cos^2\phi\\ \amp L_y^2 = (2m\omega a^2)^2 \sin^2\theta\cos^2\theta \sin^2\phi \\ \amp L_z^2 = (2m\omega a^2)^2 \sin^4\theta \end{align*}

Adding the first two we get

\begin{equation*} L_x^2 + L_y^2 = (2m\omega a^2)^2 \sin^2\theta\cos^2\theta. \end{equation*}

Adding this to $L_z^2$ involves following algebra step.

\begin{equation*} \sin^2\theta\cos^2\theta + \sin^4\theta = \sin^2\theta. \end{equation*}

Therefore, magnitude of angular momentum will be

\begin{equation*} L = 2m\omega a^2\sin\theta. \end{equation*}

Since we do not have $\phi$ in the magnitude, it is constant in time.

(Indian JEE 2016) Two thin circular discs of mass $m$ and $4m\text{,}$ having radii of $a$ and $2a\text{,}$ respectively, are rigidly fixed by a massless, rigid rod of length = $\sqrt{24}\,a$ through their centres. This assembly is laid on a firm and flat surface, and set rolling without slipping on the surface so that the angular speed about the axis of the rod is $\omega\text{.}$ The angular momentum of the entire assembly about the point 'O' is $\vec L$ (see Figure 10.2.9). Which of the following statement(s) is(are) true?

1. The centre of mass of the assembly rotates about the $z$-axis with an angular speed of $\omega/5\text{.}$
2. The magnitude of angular momentum of centre of mass of the assembly about the point O is $81 ma^2\omega\text{.}$
3. The magnitude of angular momentum of the assembly about its centre of mass is $17 ma^2\omega/2\text{.}$
4. The magnitude of the $z$-component of $L$ is $55 ma^2\omega\text{.}$
Hint

Solution 1 (A)

(A) Center of mass of the assymbly will rotate at same angular speed as the rolling of the smaller disk.

Figure 10.2.10 shows bottom point of smaller disk goes around $z$-axis in circle C1 of radius $\sqrt{l^2 + a^2} = 5a\text{.}$ That means, in the time a point on the disk makes once around the $z$-axis, the disk rotates $5$ times.

Hence, angular frequency will be one-fifth. Therefore, centre of mass of the assembly rotates about the $z$-axis with an angular speed of $\omega/5\text{.}$

Solution 2 (B)

(B) In Figure 10.2.11 I have replaced the disks by point masses at their centers. The CM will be between these two points. Temoprarily you can place origin at $m$ and use an axis towards the other mass. This gives

\begin{equation*} c = \frac{4ml}{m+4m} =\frac{4}{5}\,l. \end{equation*}

Therefore, the CM moves in a circle of radius

\begin{equation*} R = (l+c)\,\cos\theta = \frac{9l}{5}\times\frac{l}{5a} = 8.64a. \end{equation*}

It is rotating about $z$-axis at $\omega/5$ with moment of inertia $(m+4m)R^2\text{.}$ Therefore,

\begin{equation*} L_{\text{cm},z} = 5mR^2\frac{\omega}{5} = 74.6\,ma^2\omega. \end{equation*}

Therefore, (B) is incorrect.

Solution 3 (C)

(C) The two-disk structure rotates about axis through the CM at angular speed $\omega\text{.}$ The moment of inertia of the structure will be

\begin{equation*} I = \frac{1}{2}ma^2 + \frac{1}{2}(4m)(2a)^2 = \frac{17}{2}ma^2. \end{equation*}

(D) We can get $z$-component of net angular momentum $\vec L$ by adding components of the angular momentum about the axis through the CM and angular momentum of the CM about $z$-axis. We have already worked these out above.