Skip to main content

Section 7.3 Applications of unoriented line and surface integrals

While most applications of integration over curves and surfaces involve the oriented line and surface integrals that we studied throughout this course, unoriented line and surface integrals can also be useful in applications, when we want to calculate a quantity associated to a curve or a surface that should not depend on the orientation. In this section we study a few applications of unoriented line and surface integrals, such as calculating the centre of mass of a wire and a thin sheet of material.

Subsection 7.3.1 Centre of mass of a wire

In the previous sections we studied unoriented line and surface integrals of functions. Those are useful to calculate quantities associated to curves and surfaces that should not depend on a choice of orientation. The prototypical examples were the arc length of a curve and surface area of a surface, but there are many other applications.

Our first application concerns the calculation of the centre of mass of a wire in \(\mathbb{R}^n\text{.}\) Let us first recall the physical concept of centre of mass.

Suppose that there are \(k\) point particles of masses \(m_1,\ldots,m_k\) at positions \(X_1,\ldots,X_k \in \mathbb{R}\) on a line. The centre of mass of the system of particles is located at

\begin{equation*} \bar{x} = \frac{\sum_{i=1}^k m_i X_i}{m}, \end{equation*}

where \(m = \sum_{i=1}^k m_i\) is the total mass of the system. The numerator is sometimes called the “first moment of the system about the origin” and written as \(M\text{.}\)

If we are given instead a rod in \(\mathbb{R}\) between \(x=a\) and \(x=b\) with mass density \(\rho(x)\text{,}\) then to get its centre of mass we use the slicing principle. We slice the rod into small line segments. Let \(dx\) be of a typical line segment located at the point \(x\text{.}\) Its mass is \(dm = \rho(x)\ dx\text{,}\) and its first moment about the origin is \(dM = x \rho(x)\ dx\text{.}\) We then sum over slices and take the limit of an infinite number of infinitesimal slices. This turns the calculation of the total mass and first moment of the rod as a definite integral, and its centre of mass is:

\begin{equation*} \bar{x} = \frac{\int_a^b x \rho(x)\ dx}{\int_a^b \rho(x)\ dx} = \frac{1}{m} \int_a^b x \rho(x)\ dx. \end{equation*}

Now what if we want to calculate the centre of mass of a wire that is bent and twisted in \(\mathbb{R}^n\text{?}\) We apply the same idea, but thinking of the wire as a parametric curve \(\alpha: [a,b] \to \mathbb{R}^n\text{,}\) with \(\alpha(t) = (x_1(t), \ldots, x_n(t))\) and image curve \(C=\alpha([a,b])\text{.}\) Let \(\rho: C \to \mathbb{R}\) be the mass density of the wire, which we assume to be continuous. We slice the wire (the image curve \(C\)) into small curve segments. Let \(ds\) (the line element from Subsection 7.1.1) be the length of a typical curve segment located at \(\alpha(t)\text{.}\) Its mass is \(dm = \rho(\alpha(t))\ ds\text{.}\) By summing over all curve segments and taking the limit of an infinite number of segments of infinitesimal length, we calculate the total mass of the wire as an unoriented line integral:

\begin{equation*} m = \int_C \rho\ ds = \int_a^b \rho(\alpha(t)) |\mathbf{T}(t)|\ dt. \end{equation*}

To get the centre of mass, we also need to calculate the first moments of the wire. Here, we are working in \(\mathbb{R}^n\text{.}\) There are \(n\) first moments, one in each coordinate \(x_1, \ldots, x_n\text{.}\) The first moments of the curve segment are given by

\begin{equation*} x_k(t) \rho(\alpha(t))\ ds, \qquad k=1,\ldots,n. \end{equation*}

Summing over curve segments, and taking the limit as usual, we obtain that the position of the centre of mass of the wire is given by the point in \(\mathbb{R}^n\) with coordinates \(\bar{x}_1,\ldots,\bar{x}_n\text{,}\) with:

\begin{equation*} \bar{x}_k = \frac{1}{m} \int_C x_k \rho\ ds = \frac{1}{m} \int_a^b x_k(t) \rho(\alpha(t)) |\mathbf{T}(t)|\ dt, \qquad k=1,\ldots,n. \end{equation*}

In other words, to find the centre of mass of the wire, we need to evaluate the unoriented line integrals corresponding to the total mass of the wire and its \(n\) first moments.

Note that the centre of mass of the wire will not generally be located on the wire itself, since the wire is twisted and bent in the ambient space \(\mathbb{R}^n\text{;}\) this will be clear in the next example.

A wire is bent into the semi-circle \(x^2+y^2=4\text{,}\) \(y \geq 0\text{.}\) Its mass density is given by the function

\begin{equation*} \rho(x,y) = 4-y \end{equation*}

on the wire. Find the centre of mass of the wire.

Let us see what we expect first. The wire is heavier near its base (\(y=0\)) than at the top (\(y=2\)). The mass density is however symmetric about the \(y\)-axis (it does not vary in \(x\)). We thus expect the centre of mass to be located at a point \((0, \bar{y})\text{,}\) with \(\bar{y} \lt 2\text{,}\) since it should be below the top of the wire.

To calculate the centre of mass, we will need to calulate unoriented line integrals along the curve, so we first parametrize the curve as \(\alpha: [0,\pi] \to \mathbb{R}^2\) with

\begin{equation*} \alpha(\theta) = (2 \cos(\theta), 2 \sin(\theta)). \end{equation*}

The tangent vector is

\begin{equation*} \mathbf{T}_\theta = (- 2\sin(\theta), 2\cos(\theta) ), \end{equation*}

whose norm is

\begin{equation*} |\mathbf{T}_\theta| = \sqrt{4 \sin^2(\theta) + 4 \cos^2(\theta)} = 2. \end{equation*}

The line element is then

\begin{equation*} ds = |\mathbf{T}_\theta| d\theta = 2 d\theta. \end{equation*}

We first calculate the total mass of the wire. It is given by the unoriented line integral:

\begin{align*} m =\amp \int_C \rho\ ds\\ =\amp \int_0^\pi (4-2\sin(\theta)) 2 d\theta\\ =\amp 4( 2 \pi - 2)\\ =\amp 8 (\pi-1). \end{align*}

We then calculate the coordinates \(\bar{x}\) and \(\bar{y}\) of the centre of mass.

\begin{align*} \bar{x} =\amp \frac{1}{m} \int_C x \rho\ ds\\ =\amp \frac{1}{8(\pi-1)} \int_0^\pi (2 \cos(\theta) ) (4 - 2 \sin(\theta) ) 2 \ d\theta\\ =\amp 0, \end{align*}

since the trigonometric integral vanishes. As for \(\bar{y}\text{,}\) we get:

\begin{align*} \bar{y} =\amp \frac{1}{m} \int_C y \rho \ ds\\ =\amp \frac{1}{8(\pi-1)} \int_0^\pi (2 \sin(\theta)) (4-2\sin(\theta)) 2\ d\theta\\ =\amp \frac{8}{8(\pi-1)}\left(4 - \frac{\pi}{2} \right)\\ =\amp \frac{8 - \pi}{2(\pi-1)} \end{align*}

Therefore, the centre of mass of the wire is located at the point

\begin{equation*} (\bar{x}, \bar{y}) = \left(0, \frac{8 - \pi}{2(\pi-1)} \right). \end{equation*}

To make sure that this is consistent with our expectation, we can find a numerical value for the location of the centre of mass. We get:

\begin{equation*} (\bar{x}, \bar{y}) \simeq (0, 1.134). \end{equation*}

This is consistent with our expectation. Phew!

Subsection 7.3.2 Centre of mass of a thin sheet

We can do a very similar calculation to obtain the centre of mass of a thin sheet of material (such as aluminium foil, or paper) in \(\mathbb{R}^3\text{.}\) Suppose that the sheet of material takes the shape of a parametric surface \(\alpha: D \to \mathbb{R}^3\text{,}\) with \(\alpha(u,v)= (x(u,v), y(u,v), z(u,v))\) and image surface \(S = \alpha(D)\text{.}\) Suppose that \(\rho: S \to \mathbb{R}\) is the mass density of the sheet (mass per unit area). Using the same slicing approach, but with small pieces of surface of area \(dS\text{,}\) we obtain that the total mass of the sheet can be written as the unoriented surface integral

\begin{equation*} m = \iint_S \rho\ dS = \iint_D \rho(\alpha(u,v)) |\mathbf{T}_u \times \mathbf{T}_v|\ dA. \end{equation*}

We calculate the centre of mass as before, by calculating the first moments in the three coordinates \(x,y,z\text{.}\) Using the slicing process, those become unoriented surface integrals. The result is that position \((\bar{x},\bar{y},\bar{z})\) of the centre of mass is given by

\begin{align*} \bar{x} =\amp \frac{1}{m} \iint_S x \rho\ dS = \frac{1}{m} \iint_D x(u,v) \rho(\alpha(u,v)) |\mathbf{T}_u \times \mathbf{T}_v|\ dA,\\ \bar{y} =\amp \frac{1}{m} \iint_S y \rho\ dS = \frac{1}{m} \iint_D y(u,v) \rho(\alpha(u,v)) |\mathbf{T}_u \times \mathbf{T}_v|\ dA,\\ \bar{z} =\amp \frac{1}{m} \iint_S z \rho\ dS = \frac{1}{m} \iint_D z(u,v) \rho(\alpha(u,v)) |\mathbf{T}_u \times \mathbf{T}_v|\ dA. \end{align*}

Note that, as was the case for the wire, the centre of mass of the sheet is not generally expected to lie on the sheet itself, since the sheet can be bent and twisted in the ambient space \(\mathbb{R}^3\text{.}\)

Suppose the a sheet of paper is bent in the shape of the cylinder \(x^2+y^2=9\text{,}\) with \(z \in [0,1]\text{.}\) Suppose that its mass density is given by the function

\begin{equation*} \rho(x,y,z) = z+1. \end{equation*}

Find the centre of mass of the cylinder.

First, let us see what we expect. The cylinder is heavier at the top than at the bottom. However, its mass density has circular symmetry above the \(z\)-axis, as it only depends on the height \(z\) on the cylinder. Thus we expect the centre of mass to be in the middle of the cylinder, i.e. with coordinates \((0,0,\bar{z})\text{.}\) Furthermore, we expect its \(z\)-coordinate to be a little bit higher than half-way up the cylinder, since the cylinder is heavier at the top than at the bottom. So we expect \(0.5 \lt \bar{z} \lt 1\text{.}\)

To calculate the required unoriented surface integrals, we first parametrize the surface as \(\alpha:D \to \mathbb{R}^3\) with

\begin{equation*} D = \{(u,\theta) \in \mathbb{R}^2\ |\ u \in [0,1], \theta \in [0,2\pi] \} \end{equation*}

and

\begin{equation*} \alpha(u,\theta) = (3 \cos(\theta), 3 \sin(\theta), u). \end{equation*}

The tangent vectors are

\begin{equation*} \mathbf{T}_u = (0,0,1), \qquad \mathbf{T}_\theta = (- 3 \sin(\theta), 3 \cos(\theta), 0), \end{equation*}

and the cross product is

\begin{equation*} \mathbf{T}_u \times \mathbf{T}_\theta = ( - 3 \cos(\theta), - 3 \sin(\theta), 0). \end{equation*}

Its norm is

\begin{equation*} |\mathbf{T}_u \times \mathbf{T}_\theta| = \sqrt{9 \cos^2(\theta) + 9 \sin^2(\theta)} = 3, \end{equation*}

and thus the surface element is

\begin{equation*} dS = 3 du d\theta. \end{equation*}

We calculate the total mass of the sheet of paper. We get:

\begin{align*} m =\amp \iint_S \rho\ dS\\ =\amp \int_0^1 \int_0^{2 \pi} (u+1) 3\ d\theta du\\ =\amp 6 \pi \int_0^1 (u+1)\ du\\ =\amp 9 \pi. \end{align*}

As for the coordinates of the centre of mass, we get:

\begin{align*} \bar{x} =\amp \frac{1}{m}\iint_S x \rho\ dS\\ =\amp\frac{1}{9 \pi} \int_0^1 \int_0^{2 \pi} (3 \cos(\theta) ) (u+1) 3\ d \theta du \\ =\amp 0, \end{align*}

since the integral over \(\theta\) is zero. Similarly,

\begin{align*} \bar{y} =\amp \frac{1}{m}\iint_S y \rho\ dS\\ =\amp \frac{1}{9 \pi}\int_0^1 \int_0^{2 \pi} (3 \sin(\theta) ) (u+1) 3\ d \theta du \\ =\amp 0. \end{align*}

As for \(\bar{z}\text{,}\) we get:

\begin{align*} \bar{z} =\amp \frac{1}{m}\iint_S z \rho\ dS\\ =\amp \frac{1}{9 \pi} \int_0^1 \int_0^{2\pi} (u) (u+1) 3\ d \theta du\\ =\amp \frac{2}{3} \int_0^1 (u^2+u)\ du\\ =\amp \frac{5}{9}. \end{align*}

As a result, the centre of mass of the sheet of paper is located at the point

\begin{equation*} (\bar{x}, \bar{y}, \bar{z}) = \left(0,0,\frac{5}{9} \right). \end{equation*}

Since \(\frac{5}{9} \simeq 0.556,\text{,}\) this is consistent with our expectation that the centre of mass should be on the \(z\)-axis, a little bit higher than half-way up the cylinder, which is at \(z=0.5\text{.}\)

Exercises 7.3.3 Exercises

1.

Find the centre of mass of a wire that is bent in the shape of a circle of radius \(R\) centered at the origin in \(\mathbb{R}^2\text{,}\) with mass density \(\rho(x,y) = 2R - y\text{.}\) Is your result consistent with your expectations?

Solution.

Let us first see what we expect. As the wire is bent in a circle centered at the origin, if the mass density was constant (or symmetric about the \(x\)- and \(y\)-axes), the centre of mass would be at the origin. However, the mass density is not constant here: it decreases linearly with \(y\text{.}\) It is constant in the \(x\)-direction, so we expect that centre of mass to lie on the \(y\)-axis. As the mass decreases as \(y\) increases, we expect the centre of mass to be at a position \((0,\bar{y})\text{,}\) with \(-R \lt \bar{y} \lt 0\text{.}\)

We parametrize the circle of radius \(R\) as \(\alpha:[0,2\pi] \to \mathbb{R}^2\) with \(\alpha(\theta) = (R \cos(\theta), R\sin(\theta) )\text{.}\) The tangent vector is

\begin{equation*} \mathbf{T}(\theta) = (- R \sin(\theta), R \cos(\theta) ). \end{equation*}

Its norm is

\begin{equation*} |\mathbf{T}(\theta)| = \sqrt{R^2 \sin^2(\theta) + R^2 \cos^2(\theta)} = R. \end{equation*}

The total mass of the wire is

\begin{align*} m = \int_C \rho\ ds =\amp \int_0^{2\pi} (2R - R \sin(\theta) ) R d\theta\\ =\amp 4 \pi R^2. \end{align*}

The \(\bar{x}\)-coordinate of the centre of mass is

\begin{align*} \bar{x} =\amp \frac{1}{m} \int_C x \rho\ ds\\ =\amp \frac{1}{4 \pi R^2} \int_0^{2\pi} R \cos(\theta) (2 R - R \sin(\theta) ) R \ d\theta\\ =\amp 0, \end{align*}

as expected. As for the \(\bar{y}\)-coordinate, we get:

\begin{align*} \bar{y} =\amp \frac{1}{m} \int_C y \rho\ ds\\ =\amp \frac{1}{4 \pi R^2} \int_0^{2 \pi} R \sin(\theta) (2 R - R \sin(\theta) ) R\ d \theta\\ =\amp \frac{1}{4 \pi R^2} (-R^3 \pi)\\ =\amp - \frac{R}{4}. \end{align*}

Therefore, the position of the centre of mass is

\begin{equation*} (\bar{x},\bar{y}) = \left(0, - \frac{R}{4} \right). \end{equation*}

As expected, it is on the \(y\)-axis, and the \(y\)-coordinate is between \(-R\) and \(0\text{.}\)

2.

Given a wire with density \(\rho(x,y)\) that lies on a curve \(C \subset \mathbb{R}^2\text{,}\) its moments of inertia about the \(x\)- and the \(y\)-axes are defined by

\begin{equation*} I_x = \int_C y^2 \rho\ ds, \qquad I_y = \int_C x^2 \rho\ ds. \end{equation*}

Find the moments of inertia of a wire that lies along the line \(2x+y = 5\) between \(x=0\) and \(x=1\text{,}\) with density \(\rho(x,y) = x\text{.}\)

Solution.

We parametrize the line as \(\alpha: [0,1] \to \mathbb{R}^2\) with

\begin{equation*} \alpha(t) = (t,5-2 t). \end{equation*}

The tangent vector and its norm are:

\begin{equation*} \mathbf{T}(t) = (1,-2), \qquad |\mathbf{T}(t)| = \sqrt{1+2^2} = \sqrt{5}. \end{equation*}

We can then evaluate the moments of inertia. First,

\begin{align*} I_x =\amp \int_C y^2 \rho \ ds\\ =\amp \sqrt{5} \int_0^1 (5-2t)^2 t \ dt\\ =\amp \sqrt{5} \int_0^1 (25 t - 20 t^2 + 4 t^3)\ dt\\ =\amp \sqrt{5} \left( \frac{25}{2} - \frac{20}{3} + 1 \right)\\ =\amp \frac{41 \sqrt{5}}{6}. \end{align*}

Second,

\begin{align*} I_y =\amp \int_C x^2\rho \ ds\\ =\amp \sqrt{5} \int_0^1 t^3\ dt\\ =\amp \frac{\sqrt{5}}{4}. \end{align*}

3.

Consider a sheet of aluminium foil that is bent in the shape of a paraboloid \(z = x^2 + y^2\) with \(z \in [0,1]\text{.}\) Suppose that it mass density is \(\rho = k\) where \(k\) is a constant. Find the total mass of the sheet and its centre of mass. Does it agree with your expectation?

Solution.

As the mass density is constant, and the paraboloid has circular symmetry about the \(z\)-axis, we expect the centre of mass to lie on the \(z\)-axis, somewhere between \(z=0\) and \(z=1\text{.}\)

We parametrize the paraboloid as \(\alpha: D \to \mathbb{R}^3\) with

\begin{equation*} D = \{(r,\theta) \in \mathbb{R}^2\ |\ r \in [0,1], \theta \in [0,2\pi] \} \end{equation*}

and

\begin{equation*} \alpha(r,\theta) = (r \cos(\theta), r \sin(\theta), r^2). \end{equation*}

The tangent vectors are

\begin{equation*} \mathbf{T}_r = (\cos(\theta), \sin(\theta), 2 r), \qquad \mathbf{T}_\theta = (- r\sin(\theta), r\cos(\theta),0). \end{equation*}

The cross product is

\begin{equation*} \mathbf{T}_r \times \mathbf{T}_\theta = \left( - 2 r^2 \cos(\theta), - 2 r^2 \sin(\theta), r \right). \end{equation*}

Its norm is

\begin{equation*} |\mathbf{T}_r \times \mathbf{T}_\theta| = \sqrt{4 r^4 \cos^2(\theta) + 4 r^4 \sin^2(\theta) + r^2} = r \sqrt{4 r^2 + 1}. \end{equation*}

We first calculate the total mass. We get:

\begin{align*} m =\amp \iint_S \rho\ dS\\ =\amp k \int_0^1 \int_0^{2 \pi} r \sqrt{4 r^2 + 1}\ d\theta r\\ =\amp 2 \pi k \int_0^1 r \sqrt{4 r^2 + 1}\ dr\\ =\amp \frac{\pi k}{4} \int_1^{5} \sqrt{u}\ du\\ =\amp \frac{\pi k}{6} \left( 5^{3/2} - 1 \right). \end{align*}

In the process of evaluating the integral we did the substitution \(u = 4 r^2 + 1\text{.}\)

Next we calculate the coordinates of the centre of mass. First,

\begin{align*} \bar{x} =\amp \frac{1}{m} \iint_S x \rho\ dS\\ =\amp \frac{k}{m} \int_0^1 \int_0^{2\pi} r^2 \cos(\theta) \sqrt{4 r^2 + 1}\ d\theta r\\ =\amp 0, \end{align*}

as expected. Next,

\begin{align*} \bar{y} =\amp \frac{1}{m} \iint_S y \rho\ dS\\ =\amp \frac{k}{m} \int_0^1 \int_0^{2\pi} r^2 \sin(\theta) \sqrt{4 r^2 + 1}\ d\theta r\\ =\amp 0, \end{align*}

as expected. Finally,

\begin{align*} \bar{z} =\amp \frac{1}{m} \iint_S z \rho\ dS\\ =\amp \frac{6 k}{\pi k \left( 5^{3/2} - 1 \right)} \int_0^1 \int_0^{2\pi} r^3 \sqrt{4 r^2 + 1}\ d\theta r\\ =\amp \frac{12}{ \left( 5^{3/2} - 1 \right)} \int_0^1 r^3 \sqrt{4 r^2 + 1}\ dr\\ =\amp \frac{3}{8 \left( 5^{3/2} - 1 \right)} \int_1^5 (u-1) \sqrt{u}\ du\\ =\amp \frac{3}{8 \left( 5^{3/2} - 1 \right)}\left(\frac{2}{5} 5^{5/2} - \frac{2}{5} - \frac{2}{3} 5^{3/2} + \frac{2}{3} \right)\\ =\amp \frac{1}{10} \frac{5^{5/2}+1}{5^{3/2}-1}. \end{align*}

As before, to evaluate the integral we did the substitution \(u = 4 r^2+1\text{.}\)

Therefore, the centre of mass is located at

\begin{equation*} (\bar{x}, \bar{y}, \bar{z} ) = \left(0,0, \frac{1}{10} \frac{5^{5/2}+1}{5^{3/2}-1} \right) \simeq (0,0,0.559). \end{equation*}

As expected, it is located on the \(z\)-axis, somewhere between \(z=0\) and \(z=1\text{.}\)