Physical applications

$\newenvironment {prompt}{}{} \newcommand {\todo }[0]{} \newcommand {\oiint }[0]{{\large \bigcirc }\kern -1.56em\iint } \newcommand {\mooculus }[0]{\textsf {\textbf {MOOC}\textnormal {\textsf {ULUS}}}} \newcommand {\npnoround }[0]{\nprounddigits {-1}} \newcommand {\npnoroundexp }[0]{\nproundexpdigits {-1}} \newcommand {\npunitcommand }[1]{\ensuremath {\mathrm {#1}}} \newcommand {\RR }[0]{\mathbb R} \newcommand {\R }[0]{\mathbb R} \newcommand {\N }[0]{\mathbb N} \newcommand {\Z }[0]{\mathbb Z} \newcommand {\sagemath }[0]{\textsf {SageMath}} \newcommand {\d }[0]{\mathop {}\!d} \newcommand {\l }[0]{\ell } \newcommand {\ddx }[0]{\frac {d}{\d x}} \newcommand {\zeroOverZero }[0]{\ensuremath {\boldsymbol {\tfrac {0}{0}}}} \newcommand {\inftyOverInfty }[0]{\ensuremath {\boldsymbol {\tfrac {\infty }{\infty }}}} \newcommand {\zeroOverInfty }[0]{\ensuremath {\boldsymbol {\tfrac {0}{\infty }}}} \newcommand {\zeroTimesInfty }[0]{\ensuremath {\small \boldsymbol {0\cdot \infty }}} \newcommand {\inftyMinusInfty }[0]{\ensuremath {\small \boldsymbol {\infty -\infty }}} \newcommand {\oneToInfty }[0]{\ensuremath {\boldsymbol {1^\infty }}} \newcommand {\zeroToZero }[0]{\ensuremath {\boldsymbol {0^0}}} \newcommand {\inftyToZero }[0]{\ensuremath {\boldsymbol {\infty ^0}}} \newcommand {\numOverZero }[0]{\ensuremath {\boldsymbol {\tfrac {\#}{0}}}} \newcommand {\dfn }[0]{\textbf } \newcommand {\unit }[0]{\mathop {}\!\mathrm } \newcommand {\eval }[1]{\bigg [ #1 \bigg ]} \newcommand {\seq }[1]{\left ( #1 \right )} \newcommand {\epsilon }[0]{\varepsilon } \newcommand {\phi }[0]{\varphi } \newcommand {\iff }[0]{\Leftrightarrow } \DeclareMathOperator {\arccot }{arccot} \DeclareMathOperator {\arcsec }{arcsec} \DeclareMathOperator {\arccsc }{arccsc} \DeclareMathOperator {\si }{Si} \DeclareMathOperator {\scal }{scal} \DeclareMathOperator {\sign }{sign} \newcommand {\arrowvec }[1]{{\overset {\rightharpoonup }{#1}}} \newcommand {\vec }[1]{{\overset {\boldsymbol {\rightharpoonup }}{\mathbf {#1}}}} \DeclareMathOperator {\proj }{\vec {proj}} \newcommand {\veci }[0]{{\boldsymbol {\hat {\imath }}}} \newcommand {\vecj }[0]{{\boldsymbol {\hat {\jmath }}}} \newcommand {\veck }[0]{{\boldsymbol {\hat {k}}}} \newcommand {\vecl }[0]{\vec {\boldsymbol {\l }}} \newcommand {\uvec }[1]{\mathbf {\hat {#1}}} \newcommand {\utan }[0]{\mathbf {\hat {t}}} \newcommand {\unormal }[0]{\mathbf {\hat {n}}} \newcommand {\ubinormal }[0]{\mathbf {\hat {b}}} \newcommand {\dotp }[0]{\bullet } \newcommand {\cross }[0]{\boldsymbol \times } \newcommand {\grad }[0]{\boldsymbol \nabla } \newcommand {\divergence }[0]{\grad \dotp } \newcommand {\curl }[0]{\grad \cross } \newcommand {\lto }[0]{\mathop {\longrightarrow \,}\limits } \newcommand {\bar }[0]{\overline } \newcommand {\surfaceColor }[0]{violet} \newcommand {\surfaceColorTwo }[0]{redyellow} \newcommand {\sliceColor }[0]{greenyellow} \newcommand {\vector }[1]{\left \langle #1\right \rangle } \newcommand {\sectionOutcomes }[0]{} \newcommand {\HyperFirstAtBeginDocument }[0]{\AtBeginDocument } \newcommand {\epstopdfsetup }[0]{\setkeys {ETE}} changed the theorem header of amsthm failed\MessageBreak }}} \newcommand {\GetTitleStringSetup }[0]{\setkeys {gettitlestring}} \newcommand {\Sectionformat }[2]{#1}$

We apply the procedure of “Slice, Approximate, Integrate” to model physical situations.

Thus far, we have studied several geometric applications of the procedure of “Slice, Approximate, Integrate”. Indeed, it can be used to find lengths, areas, volumes. This procedure is not limited to model only geometric situations. Many problems from other STEM fields requires the same technique. This section examines many of these situations. First, we make a more generalized observation of the philosophy behind the “Slice, Approximate, Integrate” procedure.

A broader perspective

Let’s take a step back and try to think about each of the situations we’ve applied the “Slice, Approximate, Integrate” procedure to model. A fundamental aspect to develop each formula lies in knowing how to approximate each slice. For each quantity of interest - length, area, or volume- we used an object whose length, area, or volume we could calculate. For instance, for the area between curves, we know how to find the area of rectangles, so we approximate that each slice is a rectangle. For solids of revolution, we can calculate the volume of washers or shell, so we approximate that each slice is a washer or a shell.

In physics, we take measurable quantities from the real world, and attempt to find meaningful relationships between them. These relationships are often expressed using formulas, but these formulas come with inherent assumptions. By utilizing the “Slice, Approximate, Integrate” procedure, we are able to generalize the scenarios in which we can compute physical quantities of interest.

Mass of a wire with variable density

Given a physical object, the density is a measure of how the mass of the object is distributed. In the case where the density of the object is constant, we have:

3D: $\left <\textrm {density}\right > = \frac {\left <\textrm {mass}\right >}{\left <\textrm {volume}\right >}$

2D: $\left <\textrm {density}\right > = \frac {\left <\textrm {mass}\right >}{\left <\textrm {aera}\right >}$

1D: $\left <\textrm {density}\right > = \frac {\left <\textrm {mass}\right >}{\left <\textrm {length}\right >}$

We can and often do approximate physical objects, like wires, as one dimensional. In this case, the linear density is a measure of its mass per unit length.

If a wire has linear density $5 \unit {g}/\unit {m}$ , how many grams will $4 \unit {m}$ of this wire be? $\mathrm {Mass} = \answer {20} \unit {g}$

In this case, we need just multiply the density by the length.

Sometimes the linear density of an object can vary from one part of the object to another. In this case density will be a non-constant function, which we will represent by the Greek letter, rho, $\rho$ .

The kilogram ( $kg$ ) and meter ( $m$ ) are the standard SI units of mass and length. The standard unit for linear density is thus $kg/m$ .

Suppose that a wire that extends from $x=0$ to $x=2$ has density profile $\rho (x) = \frac {1}{(1+2x)^2}$ . How do we calculate its mass?

The length of the wire is 3, so the mass is $\rho (3) \cdot 3 = \frac {7}{49}$ units. We cannot use the result $\left <\textrm {density}\right > = \frac {\left <\textrm {mass}\right >}{\left <\textrm {length}\right >}$ because the density is not constant!

What should we do? Let’s try the “Slice, Approximate, Integrate” procedure.

We cut the wire into numerous pieces. One such piece of width $Delta x$ is shown.

Step 2: Approximate The only instance in which we know how to calculate the mass of a 1D object is when the density is constant. This is a rather poor approximation for the entire rod, but if each slice is very small, it is a fairly good approximation.

For example, consider the slice that extends from $x=1$ to $x=1.001$ . Using the formula $\rho = \frac {1}{(1+2x)^2}$ , we find to 5 decimal places: $\rho (1) = \answer [given, tolerance=.0005]{.11111}$ and $\rho (1.001) = \answer [given,tolerance=.0005]{ .11096}$ .

Since we know how to calculate the mass of an object with constant density, we thus approximate that:

the density along each slice is constant, but each slice may have a different density than the others the density along each slice is constant, and each slice has the same density as all of the other slices the density of the whole rod is constant.

Since the density along is slice is constant, we find that the mass $\Delta m$ of a single slice is given by:

$\Delta m = \rho (x^*) \Delta x$ where $x^*$ is the $x$ -value along the slice where the density is approximated.

Since there are $n$ slices used, let $x^*_k$ denote the $x$ -value used to determine the density of the $k$ -th slice, and $\Delta x_k$ be the width of the $k$ -slice. The total approximate mass $m$ of the rod is thus:

$m = \sum _{k=0}^n \rho (x_k^*) \Delta x_k$

The usual procedure converts this approximate mass into the exact mass; indeed the definite integral will perform the simultaneous limiting process of shrinking the width of each slice while adding up the contributions from all of them. Indeed, since the mass of a single slice was found to be $\Delta m = \rho \Delta x$ , the total exact mass is found the same way as before:

$m = \int _{x=0}^{x=2} \rho \d x$

In this case:

$\begin{align*} m &= \int_{x=0}^{x=2} \frac{1}{(1+2x)^2} \d x \\ &= \eval{\answer[given]{-\frac{1}{2}(1+2x)^{-1}} }_0^2 \\ &= \answer[given]{\frac{2}{5}} \end{align*}$

There was nothing particularly special about this example. In fact:

Suppose a thin wire has a continuous density profile given by $\rho (x)$ on $x=a$ to $x=b$ . The mass of the wire from $x=a$ to $x=b$ is given by:

$m = \int _{x=a}^{x=b} \rho \d x$

A thin wire that extends from $x=0$ to $x=4$ has density profile given by:

$\rho (x) = \left \{ \begin {array}{cl} 4x, & 0 \leq x\leq 2 \\ 8 ,& 2 < x \leq 4 \end {array} \right .$

Without performing any calculations, which half of the rod has more mass?

the left half the right half they both have equal mass

Note that the density of the left half of the rod is always less than 8 since $4x \leq 8$ for $0\leq x \leq 2$ . Thus, the left half should have less mass than the right side.

We now find the total mass of the rod.

It is still true that:

$m = \int _{x=0}^{x=4} \rho \d x$ but since $\rho$ is now a piecewise function, we will need to compute the integral in pieces. Write with me:

$\begin{align*} m = \int_{x=0}^{x=4} \rho \d x & = \int_{x=0}^{x=2} \answer[given]{4x} \d x + \int_{x=\answer{2}}^{x=\answer{4}} \answer[given]{8} \d x \\ &= \eval{\answer[given]{2x^2}}_0^2+\eval{\answer[given]{8x}}_2^4 \\ &= \answer[given]{24} \textrm{ units} \end{align*}$

Let’s now find the exact value of $a$ so that the mass of the wire from $x=0$ to $x=a$ is exactly half of the total mass of the wire.

Since we have reasoned that the right half of the rod has more mass than the left half, we expect:

$a<2$ $a=2$ a¿2

Since we expect that $a>2$ , we can find $a$ by setting:

$\int _{x=a}^{x=4} \rho \d x = \frac {1}{2} \int _{x=0}^{x=4} \rho \d x \\$

Since $a>2$ , we can use the bottom piece for $\rho$ in the first integral, and since we found that the total mass of the wire is $24$ units, we have:

$\begin{align*} \int_{x=a}^{x=4} 8 \d x &= \frac{1}{2} \cdot 24 \eval{8x}_a^4 =12 \\ \answer[given]{32-8a} = 12 \\ a &= \answer[given]{\frac{5}{2}} \end{align*}$

This value of $a$ is called the center of mass of the wire.

Note that we could have found it by setting $\int _{x=0}^{x=a} \rho \d x$ equal to 12, but this would have required two integrals since we expected that $a>2$ . Nevertheless, its not bad practice to show that the part of the rod from $x=0$ to $x=a=5/2$ is exactly half the mass of the rod:

$\int _{x=0}^{x=5/2} \rho \d x = \int _{x=0}^{x=2} 4x \d x + \int _{x=2}^{x=5/2} \rho \d x =\eval {2x^2}_0^2+\eval {8x}_0^{5/2} = 12$

Work

One of the most important concepts in physics is that of work, which measures the change in energy that occurs when a force moves an object over a certain displacement. For those familiar with physics, the Work-Energy Theorem is a powerful tool for studying situations in Newtonian mechanics (such as a box sliding down an incline plane).

For a constant force $F$ acting on a particle over a displacement on $d$ in the direction of displacement, work is given by the formula:

$\left <\textrm {Work} \right >= \left <\textrm {Force} \right > \cdot \left <\textrm {displacement} \right >$

By denoting work by $W$ , force by $F$ , and displacement by $d$ , we can write:

$W = F \cdot d$

The Newton ( $N$ ) and Joule ( $J$ ) are the standard SI units of force and work/energy. For those curious, the Newton is defined to be the amount of force required to give a mass of 1 $kg$ and acceleration of $1m/s^2$ . The Joule is then defined to be the amount of work done by a force of 1 $N$ while displacing an object 1 $m$ in the direction of the force.

While some background knowledge with physics certainly gives important context to these examples, it is not necessary to perform the required mathematics. In fact, we can think of the formula $W=F\cdot d$ as “When the quantity $F$ is constant over a displacement $d$ , $W$ is some quantity calculated by the formula $W=F \cdot d$ ”.

Note that the formula $W=F \cdot d$ comes with the requirement that the force be constant over the motion! There are many examples of important forces in physics that are not constant:

The force required to stretch a spring $x$ units from its equilibrium position is given by $F(x) = k x$ , where $k$ is a constant of proportionality.
The attractive gravitational force between two objects of masses $M$ and $m$ that are separated by a distance of $r$ units is $F(r) = \frac {GMm}{r^2}$ , where $G$ is a constant.
The electrostatic force between two particles with charges $Q$ and $q$ that are separated by a distance of $r$ units is $F(r) = \frac {kQq}{r^2}$ , where $k$ is a constant.

In what follows, we will look at two scenarios in which the formula $W=F \cdot d$ does not immediately apply. They result when either the force required to move an object is not constant, or when different parts of the object we move must be moved different distances. They can be summarized by:

Work done by a non-constant force on a particle

Suppose that a non-constant force $F(x)$ acts on an object from $x=a$ to $x=b$ . To find how much work does the force do on the object, what should we do? Let’s try the “Slice, Approximate, Integrate” procedure.

We divide the displacement into numerous pieces. One such piece of width $\Delta x$ is shown.

The only instance in which we know how to calculate the work is when force is constant. If the force is continuous and each slice is very small, it is a good approximation. Indeed, since we know how to calculate $W$ when $F$ is constant by the formula $W=F\cdot d$ , we approximate that the force along each piece:

is constant, but the force along neighboring pieces may be different. is constant, and the force along all neighboring pieces must be the same as the force on the indicated one. the force along the entire displacement is constant.

Since the force along is slice is constant, we find that the work $\Delta W$ done by the force along to displace the object $\Delta x$ along the single slice is given by:

$\Delta W = F(x^*) \Delta x$ where $x^*$ is the $x$ -value along the slice where the force is approximated.

Since there are $n$ slices used, let $x^*_k$ denote the $x$ -value used to determine the force along the $k$ -th slice, and $\Delta x_k$ be the width of the $k$ -slice. The total approximate work $W$ done by the force is thus:

$W = \sum _{k=0}^n F(x_k^*) \Delta x_k$

The usual procedure converts this approximate work into the exact work; indeed the definite integral will perform the simultaneous limiting process of shrinking the width of each slice while adding up the contributions from all of the, Indeed, since the work done by the force along a single slice was found to be $\Delta W =F \Delta x$ , the total exact mass is found by:

$W = \int _{x=a}^{x=b} F \d x$

Let’s try some examples:

A common example of a variable force is the force required to stretch a spring.

The force required to compress or stretch a spring $x$ units from its natural length (its unstretched length) is proportional to $x$ . That is, the force is given by: $F(x) = k\cdot x$ for some constant $k$ , known as the spring constant.

The SI units for the spring constant is $N/m$ (force/distance).

Units are important! Make sure that before starting a problem, all quantities are expressed in terms of the same units. If they are not, you will have to convert them!

If a a force of $1\unit {N}$ stretches a given spring $2\unit {cm}$ , then how far will a force of $5\unit {N}$ stretch the spring? The spring will stretch $\answer {10}\unit {cm}$ .

What is the spring constant in $N/m$ in this case?

Convert the distances to meters.

$k= \answer {50}\unit {N}/\unit {m}$ .

Suppose that $40N$ of force is required to stretch a spring $4m$ from its equilibrium position. How much work is done by the force?

Since the force is not constant, we want to use:

$W=\int _{x=a}^{x=b} F \d x = \int _0^4 kx \d x$

We do not yet know the value of the spring constant, but since $40 N$ is required to displace the spring $4m$ , we find by using Hooke’s Law:

$\begin{align*} F &= kx \\ 40 &= k(4) \\ k &= \answer{10} N/m \end{align*}$

Thus,

$\begin{align*} W= \int_0^4 kx \d x &= \int_0^4 10x \d x \\ &= \eval{\answer[given]{5x^2}}_0^4 \\ &= \answer[given]{80} J \end{align*}$

Work done by a constant force on a collection of particles

A cylindrical storage tank with a radius of $10 \unit {m}$ and a height of $30\unit {m}$ is filled with water, which has density $\rho =1000 \unit {kg}/\unit {m}^3$ . Compute the amount of work performed by pumping the water up to a point $5$ meters above the top of the tank. You may use the fact that $g=9.8m/s^2$ and that the force exerted by gravity on an object of mass $m$ is given by $F=mg$ .

Note that we need to overcome the force exerted by gravity to lim the liquid, which in this case is a constant $62.4 \unit {kg}/\unit {m}^3$ .

Find the weight of the liquid of the tank. Then, use $W=F \cdot d$ , where $F$ will be the weight of the liquid in the tank and $d$ will be $5$ meters. We cannot use $W=F \cdot d$ because different parts of the water have to be moved different distances!

Note that if we were liming the entire tank $35$ meters, there would be no need for calculus since each particle of water would be moved the same distance. Here, different parts of the water must be moved different distances, Indeed, the water at the top of the tank has a much shorter distance to travel than the water at the bottom of the tank.

A key observation is that each particle of water at the same height must be moved the same distance, so we will slice the tank into vertical pieces. Setting $y=0$ at the base of the tank, and examine a piece at height $y$ :

We know how to compute the work done by a constant force if we displace an object a certain distance. Here, the force is constant, but each particle in the slice must be moved a different distance. If the slice is thin enough, however, we can approximate that each particle within it must be moved the same distance. That is, we can treat the slice as a single object and find the work $\Delta W$ required to move it using $\Delta W = \Delta F \cdot d$

Note here that the slice is at a height $y$ and needs to be limed to a height of $35$ , so $d=\answer [given]{35-y}$ .

The amount of force $\Delta F$ is small because the slice width $\Delta y$ is small. We must express $\Delta F$ in terms of $y$ since we sliced with respect to $y$ . First, note that gravity exerts a force of $\Delta F = (\Delta m) g$ on the slice. Since the density of the liquid is constant, we have $\Delta m = \rho \Delta V$ . Note here that the cross-sections of the tank are circular, so $\Delta V = A \Delta y$ . We note that the cross-sectional area is constant no matter at what height $y$ the slice is drawn, so $A = \pi r^2 = 100 \pi$ .

Putting it all together:

$\Delta W = \Delta F \cdot d = (\Delta m) g (35-y) = (\Delta V) \rho g (35-y) = A (\Delta y) \rho g (35-y)$

Rearranging this gives the approximate work required to move a single slice of liquid:

$\Delta W = \rho g A (35-y) \Delta y = 98000 \pi (35-y) \Delta y$

We can now find the exact work in the usual way. We need to add up the contributions of all of the slices of liquid, which begin at $y=\answer [given]{0}$ and extend to $y=\answer {30}$ . Thus, the exact work is given by:

$W = \int _{y=0}^{y= \answer [given]{30}} 98000\pi (35-y) \d y$

Evaluating this integral is easier if we factor out the $98000\pi$ first:

$\begin{align*} W = \int_{y=0}^{y= 30} 98000\pi(35-y) \d y &= 98000\pi \int_{y=0}^{y= 30} (35-y) \d y\\ &= 98000\pi \eval{35y-\frac{1}{2}y^2}_0^{30} &= \answer[tolerance=100, given]{ 5.88 \times 10^7 \pi} J \end{align*}$

Final thoughts

The technique of “Slice, Approximate, Integrate” can be used to solve physical problems as well as geometric ones. These are not the only examples of physical problems that can be modeled and solved by using this technique. The exercises will explore other problems, and you will run into many more in other STEM courses.

“Mathematics is the abstract key which turns the lock of the physical universe” - John Polkinghome