Calculus III: Triple Integrals in Spherical Coordinates

Triple integrals are the workhorse for calculating volumes, masses, and average values in three dimensions, but their setup in rectangular coordinates can be painfully complex for regions with spherical symmetry. Spherical coordinates provide an elegant and powerful alternative, transforming intricate problems into manageable ones. This system is indispensable in engineering fields like electromagnetism, fluid dynamics, and astrophysics, where phenomena naturally radiate from a point.

The Spherical Coordinate System

The spherical coordinate system locates a point $P$ in space using three coordinates: $(\rho ,\theta ,\varphi )$. Here, $\rho$ (rho) is the distance from the origin to the point, $\theta$ (theta) is the same azimuthal angle as in cylindrical coordinates (measured from the positive x-axis), and $\varphi$ (phi) is the polar angle measured from the positive z-axis. It is crucial to note the conventions: $0\le \rho <\infty$, $0\le \theta \le 2\pi$, and $0\le \varphi \le \pi$. The relationship to rectangular coordinates $(x,y,z)$ is given by the transformation equations: $$x=\rho \sin \varphi \cos \theta ,y=\rho \sin \varphi \sin \theta ,z=\rho \cos \varphi .$$ Visualizing this, $\rho$ is the radius of a sphere, $\theta$ sweeps around the z-axis, and $\varphi$ controls the "tilt" from the North Pole. A region is best suited for spherical integration if its boundaries are spheres, cones, or planes through the origin.

The Spherical Volume Element and Jacobian

When integrating, we cannot simply use $d\rho d\theta d\varphi$. A small change in each coordinate traces out a curvy box in space, not a rectangular prism. To find the correct volume element, we calculate the Jacobian determinant of the coordinate transformation. The Jacobian factor $J$ for spherical coordinates is: $$\left|\frac{\partial (x,y,z)}{\partial (\rho ,\theta ,\varphi )}\right|={\rho }^2\sin \varphi .$$ This leads to the fundamental spherical coordinate volume element: $$dV={\rho }^2\sin \varphi d\rho d\theta d\varphi .$$ You can think of ${\rho }^2\sin \varphi d\rho d\theta d\varphi$ as the volume of an infinitesimal "spherical brick." The ${\rho }^2$ term accounts for the fact that spherical shells get larger in surface area as you move outward, and the $\sin \varphi$ term accounts for the varying circumference of horizontal slices of a sphere. Forgetting this factor is the single most common computational error.

Setting Up Integration Limits

The strategy for setting limits in spherical coordinates follows a consistent "from the inside out" approach, often described as radial, then polar, then azimuthal.

1. Radial Limit ($\rho$): For a given direction specified by fixed $\theta$ and $\varphi$, trace a ray from the origin outward. The $\rho$ limits are where the ray enters and exits the region $E$. These are often functions of $\theta$ and $\varphi$, like $\rho =0$ to $\rho =R$ for a solid sphere.
2. Polar Angle Limit ($\varphi$): After accounting for all radial distances, consider how far the direction vector can tilt from the z-axis. For a full sphere, $\varphi$ goes from $0$ (North Pole) to $\pi$ (South Pole). For a region above a cone, say $z=\sqrt{x^2+y^2}$, the cone has an angle of $\varphi =\pi /4$. The limits would be $\varphi =0$ to $\varphi =\pi /4$.
3. Azimuthal Angle Limit ($\theta$): Finally, let the entire configuration rotate around the z-axis. These are typically constants, like $0$ to $2\pi$ for full rotation, or a subset for symmetric wedges.

Example: Set up the triple integral for the volume of a sphere of radius $R$ centered at the origin.

• The ray from the origin exits the sphere at $\rho =R$. So, $\rho$ goes from $0$ to $R$.
• The angle $\varphi$ must cover the entire sphere, from $0$ to $\pi$.
• The angle $\theta$ must provide full rotation, from $0$ to $2\pi$.

The volume integral becomes: $$V={\int }_{\theta =0}^{2\pi }{\int }_{\varphi =0}^{\pi }{\int }_{\rho =0}^R{\rho }^2\sin \varphi d\rho d\varphi d\theta .$$ Evaluating this confirms the well-known formula $V=\frac{4}{3}\pi R^3$.

Converting Problems to Spherical Form

The decision to use spherical coordinates is driven by two factors: the geometry of the region and the function being integrated. A spherical region paired with a function like $f(x^2+y^2+z^2)$ is ideal. The conversion process has three key steps.

First, describe the region $E$ using spherical inequalities. For example, the solid between spheres of radii $a$ and $b$ is simply $a\le \rho \le b$. A solid cone is described by $0\le \varphi \le \alpha$ for some constant $\alpha$.

Second, convert the integrand $f(x,y,z)$ using the transformation equations: $$f(x,y,z)=f(\rho \sin \varphi \cos \theta ,\rho \sin \varphi \sin \theta ,\rho \cos \varphi ).$$ A common simplification occurs with distance terms: $x^2+y^2+z^2={\rho }^2$.

Third, write the full triple integral with the correct volume element: $${\iiint }_{E}f(x,y,z)dV={\int }_{\theta =\alpha }^{\beta }{\int }_{\varphi =g_1(\theta )}^{g_2(\theta )}{\int }_{\rho =h_1(\theta ,\varphi )}^{h_2(\theta ,\varphi )}f(\rho ,\theta ,\varphi ){\rho }^2\sin \varphi d\rho d\varphi d\theta .$$ The order of integration can be changed if it simplifies the limits, but the ${\rho }^2\sin \varphi$ factor must always precede $d\rho d\varphi d\theta$.

Applications: Gravitational and Electric Potential

The true power of spherical integration shines in physics applications. A classic problem is finding the gravitational potential or electric potential due to a solid object with a given density or charge distribution. For a point mass or charge $dm$ located at $(\rho ,\theta ,\varphi )$, its contribution to the potential at the origin is proportional to $dm/\rho$, where $\rho$ is the distance to the origin.

Consider a solid ball of radius $R$ with constant density $\delta$. To find the gravitational potential at the origin (inside the ball), we integrate the contribution from all mass elements. Here, $dm=\delta dV=\delta ({\rho }^2\sin \varphi d\rho d\theta d\varphi )$. The potential $\Phi$ is: $$\Phi \propto {\iiint }_E\frac{\delta }{\rho }dV=\delta {\int }_0^{2\pi }{\int }_0^{\pi }{\int }_0^R\frac{1}{\rho }\cdot {\rho }^2\sin \varphi d\rho d\varphi d\theta .$$ Notice the beautiful cancellation: ${\rho }^2/\rho =\rho$. The integral simplifies to: $$\Phi \propto \delta {\int }_0^{2\pi }d\theta {\int }_0^{\pi }\sin \varphi d\varphi {\int }_0^R\rho d\rho .$$ This evaluates cleanly to a finite number, demonstrating the utility of spherical coordinates. For electric potential calculations with a non-uniform charge density $\sigma (\rho )$, the setup is identical, with $\sigma (\rho )$ replacing the constant $\delta$.

Common Pitfalls

1. Incorrect Jacobian/Volume Element: Using $dV=d\rho d\theta d\varphi$ is a critical error. Always remember: $dV={\rho }^2\sin \varphi d\rho d\theta d\varphi$. A quick dimensional check can help: ${\rho }^2$ has units of length${}^2$, so the product has units of volume.
2. Misidentifying the Angle $\varphi$: Confusing the polar angle $\varphi$ (measured from the z-axis) with an angle in the xy-plane is common. Remember, for a cone $z=\sqrt{x^2+y^2}$, the relationship is $z=\rho \cos \varphi$ and $\sqrt{x^2+y^2}=\rho \sin \varphi$. Setting them equal gives $\cos \varphi =\sin \varphi$, so $\varphi =\pi /4$.
3. Wrong Integration Order and Limits: The most logical order is typically $d\rho$ then $d\varphi$ then $d\theta$. When finding limits, sketch the region. The $\rho$ limit often depends on the surface (e.g., $\rho =\sec \varphi$ for a plane $z=1$), not just a constant. Always ensure your limits describe the entire region.
4. Forgetting the $\sin \varphi$ in the Integral: Even if you remember ${\rho }^2$, the $\sin \varphi$ factor can be missed, especially when $\varphi$ limits are constants. This factor is essential and comes directly from the Jacobian. When you integrate $d\varphi$, the integral of $\sin \varphi$ from $0$ to $\pi$ is 2, which is why the surface area of a sphere is $4\pi R^2$ and not $2{\pi }^2{R}^2$.

Summary

• Spherical coordinates $(\rho ,\theta ,\varphi )$ are optimal for integrating over regions with spherical symmetry, such as spheres, balls, and cones.
• The correct volume element is $dV={\rho }^2\sin \varphi d\rho d\theta d\varphi$, derived from the absolute value of the Jacobian determinant for the coordinate transformation.
• Setting limits follows a radial-first strategy: $\rho$ from origin to surface, $\varphi$ from axis to bounding cone, and $\theta$ for full or partial rotation.
• Converting a problem involves translating both the region's boundaries and the integrand function into the spherical coordinate language.
• Key applications in engineering and physics include calculating gravitational potential and electric potential for symmetric mass and charge distributions, where the math simplifies elegantly.