Calculus III: Triple Integrals

Moving from single to double integrals lets you calculate areas and accumulate quantities over two-dimensional regions. Triple integrals extend this powerful idea into the third dimension, allowing you to model and compute physical properties of three-dimensional solids. Whether you're determining the mass of an irregularly shaped object with varying density, finding its balance point, or analyzing its resistance to rotation, triple integrals provide the essential mathematical framework. Mastering them is crucial for advanced engineering fields, including fluid dynamics, thermodynamics, and structural mechanics, where thinking and calculating in three dimensions is the norm.

Definition and Geometric Interpretation

At its core, a triple integral is a natural extension of the Riemann sum concept to three dimensions. Given a scalar function $f(x,y,z)$ defined over a solid region $E$ in space, the triple integral represents the limit of a sum of products $f(x_i,y_j,z_k)\Delta V_{ijk}$, where $\Delta V$ is the volume of a small rectangular box within the region.

$${\iiint }_{E}f(x,y,z)dV=\underset{\max \Delta V_{ijk}\to 0}{\lim }\sum _{i=1}^m\sum _{j=1}^n\sum _{k=1}^{p}f(x_i^{\ast },y_j^{\ast },z_k^{\ast })\Delta V_{ijk}$$

Geometrically, if the function is simply $f(x,y,z)=1$, then the triple integral calculates the volume of the solid $E$:

$$V={\iiint }_{E}1dV$$

This is the 3D analogue to using a single integral for area under a curve. For any other function, the integral accumulates the quantity described by $f$ throughout the volume—imagine adding up tiny contributions of density, charge, or temperature from every infinitesimal piece of the solid.

Evaluating Iterated Triple Integrals

To compute a triple integral, we typically evaluate it as an iterated triple integral using Fubini's Theorem, which allows us to integrate in any order (given the function is continuous and the bounds are properly chosen). The notation $dV$ can be $dzdydx$, $dydxdz$, etc., depending on the chosen order of integration.

The most common setup for a region $E$ is when it is "Type 1": bounded between two surfaces. For instance, if $E$ lies between graphs $z=u_1(x,y)$ and $z=u_2(x,y)$ over a 2D region $D$ in the $xy$-plane, the integral becomes:

$${\iiint }_{E}f(x,y,z)dV={\iint }_D\left[{\int }_{u_1(x,y)}^{u_2(x,y)}f(x,y,z)dz\right]dA$$

You then evaluate the remaining double integral over $D$ using bounds in $x$ and $y$. The key skill is correctly identifying these bounds by projecting the 3D solid onto a coordinate plane and describing the projection with inequalities.

Example: Volume of a Tetrahedron Find the volume of the tetrahedron bounded by the coordinate planes and the plane $x+2y+3z=6$.

1. The solid is bounded below by $z=0$ and above by the plane $z=(6-x-2y)/3$.
2. The projection $D$ in the $xy$-plane is the triangle bounded by $x=0$, $y=0$, and the line $x+2y=6$ (from setting $z=0$).
3. Describe $D$: For a given $x$ from 0 to 6, $y$ ranges from 0 to $(6-x)/2$.
4. Set up and evaluate the iterated integral:

$$V={\int }_{x=0}^6{\int }_{y=0}^{(6-x)/2}{\int }_{z=0}^{(6-x-2y)/3}1dzdydx$$ Solving the innermost integral gives $\frac{6-x-2y}{3}$. Proceeding through the $y$ and $x$ integrations yields a final volume of 6 cubic units.

Applications: Mass, Center of Mass, and Moments

Triple integrals truly shine in their physical applications. The most direct is calculating the mass of a solid with non-uniform variable density $\rho (x,y,z)$ (mass per unit volume):

$$m={\iiint }_E\rho (x,y,z)dV$$

Once mass is known, you can find its balance point, the center of mass $(\bar{x},\bar{y},\bar{z})$. The coordinates are weighted averages of position. For example, the first moment about the $yz$-plane is $M_{yz}={\iiint }_{E}x\rho (x,y,z)dV$, and $\bar{x}=M_{yz}/m$. The formulas are:

$$\bar{x}=\frac{1}{m}{\iiint }_{E}x\rho dV,\bar{y}=\frac{1}{m}{\iiint }_{E}y\rho dV,\bar{z}=\frac{1}{m}{\iiint }_{E}z\rho dV$$

Related to this are the moments of inertia, which measure a solid's resistance to rotational acceleration about an axis. For rotation about the $z$-axis, the moment of inertia is the integral of the square of the distance from the axis ($x^2+y^2$) times the density:

$$I_z={\iiint }_E(x^2+y^2)\rho (x,y,z)dV$$

Similarly, $I_x={\iiint }_E(y^2+z^2)\rho dV$ and $I_y={\iiint }_E(x^2+z^2)\rho dV$. These are critical values in mechanical design for predicting rotational dynamics.

Setting Up Bounds for Complex Regions

The most challenging aspect of triple integrals is often not the integration itself, but correctly describing the region $E$ with limits of integration. A systematic approach is essential.

1. Visualize and Sketch: Understand the bounding surfaces (planes, cylinders, spheres, etc.).
2. Choose an Integration Order: Select an order ($dzdydx$, etc.) that makes the bounds as simple as possible. If the top and bottom surfaces are functions of $x$ and $y$, integrate with respect to $z$ first.
3. Project onto a Coordinate Plane: After "integrating out" the first variable, you are left with the projection of $E$ onto a plane. For a $dz$-first order, this is the $xy$-plane. Sketch this 2D region $D$.
4. Describe the 2D Projection ($D$): Find the bounds for the next two variables by describing $D$ using the remaining coordinates, just as you would for a double integral.

For regions with cylindrical or spherical symmetry, switching to cylindrical coordinates $(r,\theta ,z)$ or spherical coordinates $(\rho ,\theta ,\varphi )$ can dramatically simplify both the integrand and the bounds. Remember that the volume differential changes: $dV=rdrd\theta dz$ in cylindrical coordinates and $dV={\rho }^2\sin \varphi d\rho d\theta d\varphi$ in spherical coordinates.

Example: Cylindrical Region Bound Consider the solid bounded by the paraboloid $z=x^2+y^2$ and the plane $z=4$. Using cylindrical coordinates is ideal because the surfaces become $z=r^2$ and $z=4$. The projection $D$ is the circle $r^2\le 4$, so $0\le r\le 2$ and $0\le \theta \le 2\pi$. The $z$-bounds for a fixed $(r,\theta )$ are from $z=r^2$ to $z=4$. An integral over this region becomes: $${\int }_{\theta =0}^{2\pi }{\int }_{r=0}^2{\int }_{z=r^2}^{4}f(r\cos \theta ,r\sin \theta ,z)rdzdrd\theta$$

Common Pitfalls

1. Incorrect Order of Integration: The most frequent error is setting bounds that do not correspond to the chosen order. Remember that the limits for the inner integral can depend on the outer variables, but the limits for an outer integral must be constants. Always check that your bounds describe every point in the region $E$ exactly once.

• Correction: After writing bounds, test them with a point inside the region. Verify it satisfies all inequalities. Also, check a point on the boundary to ensure your equations match the surfaces.

1. Forgetting the Jacobian: When changing from rectangular coordinates $(x,y,z)$ to cylindrical or spherical coordinates, you must multiply the integrand by the appropriate factor ($r$ or ${\rho }^2\sin \varphi$). Omitting this is equivalent to using a wrong volume element.

• Correction: Treat $dV=dxdydz$ as a variable to be transformed. The Jacobian determinant of the coordinate transformation provides the necessary scaling factor, which is $r$ for cylindrical and ${\rho }^2\sin \varphi$ for spherical coordinates.

1. Misidentifying the Projection Region: When integrating with respect to $z$ first, a common mistake is to use the intersection curve of the top and bottom surfaces as the boundary for the 2D projection $D$. This is only correct if the surfaces meet. Often, the projection is determined by where the side of the solid intersects the plane $z=0$ (or another constant).

• Correction: To find the projection $D$, set the equations for the top and bottom surfaces equal to each other. This gives the curve *within the $xy$-plane* that is the "shadow" of where the surfaces meet. The full region $D$ is the area enclosed by this curve.

1. Mixing Density and Volume: When finding mass, center of mass, or moments, students sometimes forget to include the density function $\rho (x,y,z)$, especially when it is variable. Using $f(x,y,z)=1$ calculates volume, not mass.

• Correction: Always pause and state what physical quantity the integrand represents. For mass, the integrand is "little bit of density times little bit of volume," so it is $\rho dV$. For a $z$-coordinate of the center of mass, it's "little bit of $z$-moment" or $z\cdot \rho dV$.

Summary

• Triple integrals ${\iiint }_{E}fdV$ generalize integration to three dimensions, accumulating a quantity described by $f(x,y,z)$ over a solid region $E$.
• Evaluation is done via iterated integrals, where the crucial step is correctly establishing the bounds of integration by projecting the 3D solid onto a coordinate plane.
• A fundamental application is computing the volume of a solid when $f=1$, and the mass when $f$ is a density function $\rho (x,y,z)$.
• From mass, you can calculate the center of mass (the weighted average position) and moments of inertia (which quantify resistance to rotation), both essential for engineering analysis.
• Mastering the setup requires careful visualization and a systematic approach to describing 3D regions with inequalities. For symmetric regions, cylindrical or spherical coordinates simplify the process, but you must remember to include the correct Jacobian factor ($r$ or ${\rho }^2\sin \varphi$) in the volume element $dV$.