Unit 4, chapter 17

Build a solid from cylindrical shells

How the shell method finds the volume of a solid of revolution from strips parallel to the axis, each sweeping a thin cylindrical shell of about circumference times height times thickness, with rotation about the y-axis.

  • Chapter 17
  • Area between curves and volumes by washers and shells

Rotate the region under y=e−x2y = e^{-x^2} and above the xx-axis, for 0≤x≤10 \le x \le 1, a full turn about the yy-axis. Slices perpendicular to that axis would need xx as a function of yy and two different formulas, while a vertical strip keeps the simple height e−x2e^{-x^2}. When it turns, such a strip sweeps out a thin cylindrical shell. This page finds the volume of a shell, builds the volume integral 2π∫01xe−x2 dx2\pi\int_0^1 x e^{-x^2}\,dx from the radius, height and thickness of the shells, and evaluates it to π(1−e−1)\pi(1 - e^{-1}).

Key terms

Cylindrical shell

Definition A cylindrical shell is the solid between two circular cylinders with a common axis and a common height hh, with inner radius r1r_1 and outer radius r2r_2, where 0≤r1<r20 \le r_1 < r_2. Its volume is the difference of two cylinder volumes,

V=πr22h−πr12h=2π rˉ h Δr,\begin{aligned} V &= \pi r_2^2 h - \pi r_1^2 h \\[0.4em] &= 2\pi\,\bar r\,h\,\Delta r, \end{aligned}

where rˉ=12(r1+r2)\bar r = \tfrac12(r_1 + r_2) is the average radius and Δr=r2−r1\Delta r = r_2 - r_1 is the thickness. The estimate 2πr1h Δr2\pi r_1 h\,\Delta r, which uses the inner circumference, is too small by exactly πh(Δr)2\pi h (\Delta r)^2.

In plain words Cut the wall of a thin tube along a seam parallel to its axis and unroll it. The result is nearly a flat slab whose length is the circumference, whose height is hh and whose thickness is Δr\Delta r. The factor 2πr2\pi r is a circumference, not the area of a disk. The inner and outer circumferences differ by 2π Δr2\pi\,\Delta r, which is small for a thin shell, and measuring the circumference at the average radius makes the slab picture exact.

Example A shell with inner radius 11, outer radius 65\tfrac65 and height 33 has volume

V=π(3625−1)⋅3=33π25.\begin{aligned} V &= \pi\left(\tfrac{36}{25} - 1\right) \cdot 3 \\[0.4em] &= \frac{33\pi}{25}. \end{aligned}

The average radius 1110\tfrac{11}{10} and the thickness 15\tfrac15 give 2π⋅1110⋅3⋅15=33π252\pi \cdot \tfrac{11}{10} \cdot 3 \cdot \tfrac15 = \tfrac{33\pi}{25} as well. The estimate with the inner circumference gives 6π5=30π25\tfrac{6\pi}{5} = \tfrac{30\pi}{25}, which is too small by π⋅3⋅125=3π25\pi \cdot 3 \cdot \tfrac{1}{25} = \tfrac{3\pi}{25}.

Strip parallel to the axis

Definition In a region rotated about an axis, a thin strip parallel to the axis stays between two nearby distances from it, so it sweeps out approximately a thin cylindrical shell. A thin strip perpendicular to the axis sweeps out approximately a washer or a disk. For rotation about the yy-axis the parallel strips are vertical, their positions are values of xx and their thickness is Δx\Delta x. The shell method uses parallel strips and integrates in the direction perpendicular to the axis.

In plain words Pick the direction of the strips first, and the method follows. Strips across the axis give washers, and strips along the axis give shells. Choose the direction in which the length of the strip has one simple formula across the whole region.

Example For the region 0≤y≤e−x20 \le y \le e^{-x^2}, 0≤x≤10 \le x \le 1, rotated about the yy-axis, the vertical strip at xx has length e−x2e^{-x^2} for every xx in [0,1][0, 1]. Horizontal slices would need two formulas. Since e−x2e^{-x^2} decreases from 11 to e−1e^{-1}, a slice at a height yy with 0≤y≤e−10 \le y \le e^{-1} runs across the whole interval from x=0x = 0 to x=1x = 1, while a slice at a height between e−1e^{-1} and 11 stops at x=−ln⁡yx = \sqrt{-\ln y}. Rotated, these slices are disks, and

V=π∫01/e1 dy+π∫1/e1(−ln⁡y) dy,\begin{aligned} V &= \pi \int_0^{1/e} 1\,dy \\[0.4em] &\quad + \pi \int_{1/e}^1 (-\ln y)\,dy, \end{aligned}

which needs the inverse function and an integration by parts. The vertical strips avoid both.

Shell radius and height

Definition Let a region g(x)≤y≤f(x)g(x) \le y \le f(x), a≤x≤ba \le x \le b, be rotated about a vertical axis x=kx = k that does not pass through the interior of the region. For the vertical strip at xx, the shell radius is the distance ∣x−k∣\lvert x - k \rvert from the strip to the axis, and the shell height is the length of the strip, top minus bottom, h(x)=f(x)−g(x)h(x) = f(x) - g(x). The thickness is the width Δx\Delta x of the strip, and the shell has volume approximately

2π∣x−k∣ h(x) Δx.2\pi \lvert x - k \rvert\,h(x)\,\Delta x.

For the yy-axis, k=0k = 0, and the radius is xx when x≥0x \ge 0.

In plain words Three separate measurements go into each shell. The radius says how far the strip is from the axis, the height says how long the strip is, and the thickness says how wide it is. The radius comes from the position of the strip, and the height comes from the curves.

Example For the region under y=e−x2y = e^{-x^2} on [0,1][0, 1] rotated about the yy-axis, the strip at xx has radius xx, height e−x2e^{-x^2} and thickness Δx\Delta x, so its shell has volume approximately 2πxe−x2 Δx2\pi x e^{-x^2}\,\Delta x. At x=12x = \tfrac12 the radius is 12\tfrac12 and the height is e−1/4≈0.779e^{-1/4} \approx 0.779. Rotating the same region about the line x=−1x = -1 would keep the height e−x2e^{-x^2} but change the radius to x+1x + 1.

Rotation about the y-axis

Definition Rotating a region of the half-plane x≥0x \ge 0 a full turn about the yy-axis moves each point (x,y)(x, y) of the region around a circle of radius xx that is centered on the axis at height yy and lies in the plane perpendicular to the axis at that height. The shell formula with radius xx assumes the region lies in x≥0x \ge 0. If a region crosses the axis, its parts on the two sides sweep overlapping solids, and the solid of revolution is the one swept by the right part together with the mirror image of the left part.

In plain words The yy-axis is vertical, so strips parallel to it are vertical and the shell radius is simply the xx-coordinate of the strip. Changing the axis changes the solid, even when the region stays the same.

Example The region under y=e−x2y = e^{-x^2} on [0,1][0, 1] gives two different solids. About the yy-axis it gives the solid of this page. About the xx-axis its slices are disks of radius e−x2e^{-x^2}, and the volume

π∫01e−2x2 dx≈1.879\pi\int_0^1 e^{-2x^2}\,dx \approx 1.879

has an integrand with no elementary antiderivative, so it is found numerically. The shell integral about the yy-axis carries an extra factor xx, which lets a substitution finish it. The region under the same curve on [−1,1][-1, 1] sweeps exactly the same solid about the yy-axis as the region on [0,1][0, 1], because its left half is the mirror image of its right half.

Shell method

Definition Let ff and gg be continuous on [a,b][a, b], where 0≤a<b0 \le a < b, with g(x)≤f(x)g(x) \le f(x) for every xx in [a,b][a, b], and write h(x)=f(x)−g(x)h(x) = f(x) - g(x) for the shell height. Rotating the region of points with a≤x≤ba \le x \le b and g(x)≤y≤f(x)g(x) \le y \le f(x) about the yy-axis gives a solid of volume

V=2π∫abx h(x) dx.V = 2\pi \int_a^b x\,h(x)\,dx.

It is the limit of the sums of shell volumes 2πxˉi h(xˉi) Δxi2\pi \bar x_i\,h(\bar x_i)\,\Delta x_i, with xˉi\bar x_i the midpoint of the iith strip, as the widest strip shrinks to width 00. For a vertical axis x=kx = k with the region on one side, the radius xx becomes ∣x−k∣\lvert x - k \rvert. For a horizontal axis the strips are horizontal and the integral is in yy.

In plain words Multiply circumference by height by thickness for each strip, and let the integral add the shells from the innermost to the outermost. The method needs the region on one side of the axis, and it is the natural choice when the strip lengths have one formula in the variable that measures distance from the axis.

Example For the region 0≤y≤e−x20 \le y \le e^{-x^2}, 0≤x≤10 \le x \le 1, rotated about the yy-axis, the radius is xx and the height is e−x2e^{-x^2}, so

V=2π∫01xe−x2 dx.V = 2\pi \int_0^1 x e^{-x^2}\,dx.

The substitution u=x2u = x^2 is continuously differentiable, has du=2x dxdu = 2x\,dx and sends the bounds 00 and 11 to 00 and 11, so

V=π∫01e−u du=π[−e−u]01=π(1−e−1)≈1.986.\begin{aligned} V &= \pi \int_0^1 e^{-u}\,du \\[0.4em] &= \pi\bigl[-e^{-u}\bigr]_0^1 \\[0.4em] &= \pi(1 - e^{-1}) \approx 1.986. \end{aligned}

The volume is in cubic units. Two cylinders of radius 11 bound it. Since e−1≤e−x2≤1e^{-1} \le e^{-x^2} \le 1 on [0,1][0, 1], the solid contains the cylinder of height e−1e^{-1} and fits inside the cylinder of height 11, so πe−1≤V≤π\pi e^{-1} \le V \le \pi, roughly 1.156≤V≤3.1421.156 \le V \le 3.142. The disk integral from the card on parallel strips agrees, since ∫1/e1(−ln⁡y) dy=1−2e−1\int_{1/e}^1 (-\ln y)\,dy = 1 - 2e^{-1} and πe−1+π(1−2e−1)=π(1−e−1)\pi e^{-1} + \pi(1 - 2e^{-1}) = \pi(1 - e^{-1}).

Check yourself

  1. Rotate the region under y=1−xy = 1 - x for 0≤x≤10 \le x \le 1 about the yy-axis. Find the volume with shells, and compare it with the volume of a cone.
Show the answer

The strip at xx has radius xx and height 1−x1 - x, so

V=2π∫01x(1−x) dx=2π(12−13)=π3.\begin{aligned} V &= 2\pi\int_0^1 x(1 - x)\,dx \\[0.4em] &= 2\pi\left(\frac12 - \frac13\right) = \frac{\pi}{3}. \end{aligned}

The region is the triangle with vertices (0,0)(0, 0), (1,0)(1, 0) and (0,1)(0, 1), so the solid is a cone with base radius 11 and height 11, and the cone formula 13πr2h\tfrac13 \pi r^2 h gives the same π3\tfrac{\pi}{3}.

  1. Rotate the region under y=sin⁡(x2)y = \sin(x^2) for 0≤x≤π0 \le x \le \sqrt{\pi} about the yy-axis. First check that the shell height is nonnegative, then find the volume.
Show the answer

For 0≤x≤π0 \le x \le \sqrt{\pi} the inner value x2x^2 lies in [0,π][0, \pi], where the sine is nonnegative, so the height sin⁡(x2)\sin(x^2) is nonnegative. With radius xx,

V=2π∫0πxsin⁡(x2) dx.V = 2\pi\int_0^{\sqrt{\pi}} x \sin(x^2)\,dx.

With u=x2u = x^2, du=2x dxdu = 2x\,dx and the bounds 00 and π\pi,

V=π∫0πsin⁡u du=π[−cos⁡u]0π=2π.\begin{aligned} V &= \pi\int_0^{\pi} \sin u\,du \\[0.4em] &= \pi\bigl[-\cos u\bigr]_0^{\pi} = 2\pi. \end{aligned}

The solid fits inside a cylinder of radius π\sqrt{\pi} and height 11, whose volume is π2≈9.87\pi^2 \approx 9.87, and 2π≈6.282\pi \approx 6.28 is smaller, as it must be.

  1. Rotate the region under y=4−x2y = 4 - x^2 for −2≤x≤2-2 \le x \le 2 about the yy-axis. A classmate writes
V=2π∫−22x(4−x2) dx=0.V = 2\pi\int_{-2}^{2} x(4 - x^2)\,dx = 0.

What went wrong, and what is the volume?

Show the answer

The region crosses the axis. On the left half the factor xx is negative, so it is not a radius there, and the integrand is odd, so the two halves cancel. The left half is the mirror image of the right half and sweeps the same solid, so only [0,2][0, 2] counts, and

V=2π∫02x(4−x2) dx=2π(8−4)=8π.\begin{aligned} V &= 2\pi\int_0^2 x(4 - x^2)\,dx \\[0.4em] &= 2\pi(8 - 4) = 8\pi. \end{aligned}

Replacing xx by ∣x∣\lvert x \rvert over [−2,2][-2, 2] would give 16π16\pi, which counts the same solid twice. Washers agree with 8π8\pi, since the horizontal slices are disks of squared radius 4−y4 - y and π∫04(4−y) dy=8π\pi\int_0^4 (4 - y)\,dy = 8\pi.