Volume Integrals for 3D Properties

Hello! Welcome to your next lesson in our "Mathematical Methods for Physics" module.

In the last lesson, we mastered surface integrals, which allow us to calculate the flux of a vector field through a 2D surface. This involved summing up the component of the field perpendicular to the surface. Today, we take the logical next step by extending integration to three dimensions.

Our goal is to evaluate volume integrals to determine properties of three-dimensional regions. We'll learn how to integrate a scalar function over a 3D volume, a technique essential for calculating physical quantities like the mass of a non-uniform object, the total charge in a distribution, or the moment of inertia of a solid body. This is the final type of integral we'll explore before we unite all our vector calculus concepts with the powerful integral theorems.

1. The Concept of a Volume Integral

Just as a line integral sums a function along a curve and a surface integral sums it over an area, a volume integral sums a scalar function over a three-dimensional region. We denote it as:
$\int_{V} ϕ (x) d V$
Here, V is the 3D region of integration, $ϕ (x)$ is a scalar function (or scalar field) defined at each point $x = (x, y, z)$ in the volume, and dV is an infinitesimally small volume element.

The physical meaning depends on the function $ϕ$ :

If $ϕ (x) = 1$ , the integral simply gives the total volume of the region V.
If $ϕ (x) = ρ_{m} (x)$ , the mass density, the integral gives the total mass M.
If $ϕ (x) = ρ_{q} (x)$ , the charge density, the integral gives the total charge Q.

Computationally, we evaluate a volume integral as an iterated or "triple" integral. Let's begin by seeing how this works in the familiar Cartesian coordinate system.

2 Surfaces (and Volumes)

For a formal introduction to volume integrals, please read the following section from David Tong's lecture notes on Vector Calculus. It clearly defines the concept and how it's evaluated through successive integrations.

Please read Section 2.1.3, 'Volume Integrals' (page 34). Focus on understanding the conceptual definition and the two equivalent ways of setting up the iterated integral shown in Figure 5.

2. Volume Integrals in Cartesian Coordinates

As the notes explain, we evaluate ∫_V φ dV by breaking it down into three successive one-dimensional integrals. In Cartesian coordinates, dV = dx dy dz. The primary challenge is determining the limits for each integral based on the geometry of the volume V.

The general strategy is:

Innermost Integral (e.g., dz): The limits are surfaces. You integrate from the "bottom" surface z₁ = g₁(x,y) to the "top" surface z₂ = g₂(x,y) of the volume.
Middle Integral (e.g., dy): The limits are curves. After the first integration, the problem is conceptually "flattened" into a 2D region (the projection of V onto the xy-plane). You integrate from the "lower" curve y₁ = f₁(x) to the "upper" curve y₂ = f₂(x) of this 2D projection.
Outermost Integral (e.g., dx): The limits are constants, representing the overall extent of the 2D projection along the x-axis.

Let's watch a detailed example of this process.

Introduction to triple integral finding bounds

This video provides an excellent, clear demonstration of finding the volume of a tetrahedron by setting up and solving a triple integral in Cartesian coordinates.

Watch the video from the beginning up to 20:02. The presenter will walk you through: (0:02 - 4:56) Visualizing the 3D region (a tetrahedron) from its planar equation. (5:59 - 11:42) Setting up the bounds for the innermost integral (dz), from the bottom surface (z=0) to the top surface (the plane). (11:42 - 14:09) Projecting the volume onto the xy-plane and finding the bounds for the dy and dx integrals. (14:09 - 20:02) Evaluating the complete iterated integral step-by-step. This methodical approach is fundamental to solving volume integrals.

3. Changing Coordinate Systems

For regions with spherical or cylindrical symmetry, forcing the problem into Cartesian coordinates often leads to very difficult integrals. Just as we used polar coordinates for 2D integrals, we can use spherical or cylindrical polar coordinates for 3D integrals.

When we change coordinates, the volume element dV also changes. The conversion factor is given by the absolute value of the Jacobian determinant of the coordinate transformation.

$d V = d x d y d z = ∣ J ∣ d u d v d w$
where (u, v, w) are the new coordinates and J is the Jacobian.

2 Surfaces (and Volumes)

Let's return to David Tong's notes for a concise and rigorous introduction to changing coordinates for volume integrals. This section provides the key formulas you will use constantly in physics.

Please read Section 2.1.4, 'Spherical Polar and Cylindrical Polar Coordinates' (pages 35-36). Pay close attention to: The transformation equations relating Cartesian (x,y,z) to spherical (r, θ, φ) and cylindrical (ρ, φ, z) coordinates. The resulting volume elements: dV = r² sinθ dr dθ dφ for spherical and dV = ρ dρ dφ dz for cylindrical coordinates.

A Quick Note on Notation: The physics community almost universally uses θ for the polar angle (measured from the z-axis, 0 ≤ θ ≤ π) and φ for the azimuthal angle (measured in the xy-plane, 0 ≤ φ < 2π). Some math texts swap these. We will stick to the physics convention used in the notes above.

Evaluating Triple Integrals Using Spherical Coordinates for Volume — This image shows a region bounded by a hemisphere and a cone. Calculating the volume of such a region is much simpler in spherical coordinates, where the boundaries are often described by constant values of `r`, `θ`, or `φ`.

4. Applications and Worked Examples

Let's apply this to a foundational problem: calculating the volume of a sphere. This demonstrates the power of choosing the right coordinate system.

Finding Volume of a Sphere using Triple Integrals in Spherical Coordinates

The following video derives the spherical volume element from first principles and then uses it to find the volume of a sphere. This provides excellent intuition for where the r² sin(angle) term comes from.

Please watch the video from 6:32 to 22:06. (6:32 - 12:36) This part provides a beautiful geometric derivation of the volume element dV in spherical coordinates. (12:36 - 22:06) Here, the presenter sets up and evaluates the triple integral for the volume of a sphere. Important: The presenter uses φ for the polar angle and θ for the azimuthal angle. This is the reverse of our convention. Simply swap them in your mind as you watch: his φ is our θ, and his θ is our φ.

Calculating Physical Properties

We can integrate any scalar function, not just 1. A common application in mechanics is finding the mass or center of mass of an object with non-uniform density ρ(r).

Total Mass: $M = \int_{V} ρ (r) d V$
Center of Mass: $R_{C M} = \frac{1}{M} \int_{V} r ρ (r) d V$

Note that the center of mass calculation is a vector integral. In practice, you compute each component separately:
$X_{C M} = \frac{1}{M} \int_{V} x ρ (r) d V, Y_{C M} = \frac{1}{M} \int_{V} y ρ (r) d V, Z_{C M} = \frac{1}{M} \int_{V} z ρ (r) d V$

Test your understanding!

Consider a solid sphere of radius R centered at the origin. Its mass density is not uniform, but varies with the distance from the center as ρ(r) = ρ₀ (r/R), where ρ₀ is a constant. Set up the integral for the total mass M of the sphere. You don't need to solve it, just provide the complete triple integral with the correct integrand and bounds.

Show answer

This problem has spherical symmetry, so we must use spherical coordinates.

Integrand: The density is ρ(r) = ρ₀ (r/R).
Volume Element: dV = r² sinθ dr dθ dφ.
Bounds: For a full sphere, r goes from 0 to R, θ (polar angle) goes from 0 to π, and φ (azimuthal angle) goes from 0 to 2π.

The integral for the total mass is:
$M = \int_{V} ρ d V = \int_{0}^{2 π} \int_{0}^{π} \int_{0}^{R} (ρ_{0} \frac{r}{R}) (r^{2} sin θ) d r d θ d ϕ$
Combining terms inside the integral gives:
$M = \frac{ρ _{0}}{R} \int_{0}^{2 π} d ϕ \int_{0}^{π} sin θ d θ \int_{0}^{R} r^{3} d r$

Conclusion

In this lesson, we completed our journey through the fundamental types of integrals in vector calculus by adding volume integrals to our toolkit. You now have the ability to integrate scalar functions over lines, surfaces, and volumes.

Key Takeaways:

Volume integrals, ∫_V φ dV, are used to sum a scalar quantity φ over a 3D region V.
They are computed as iterated triple integrals. The main challenge is determining the integration bounds from the region's geometry.
For regions with symmetry, changing to cylindrical or spherical coordinates is crucial. This requires transforming the volume element: dV = ρ dρ dφ dz (cylindrical) or dV = r² sinθ dr dθ dφ (spherical).
Volume integrals are used to calculate essential physical properties like volume (φ=1), mass (φ=ρ_m), total charge (φ=ρ_q), and center of mass.

Preview of the Next Lesson:
We have now studied all the key differential operators (grad, div, curl) and all the key integral types (line, surface, volume). The grand finale of this module is to connect them. In the next lesson, we will learn the Divergence Theorem. This remarkable theorem states that the flux of a vector field through a closed surface is equal to the volume integral of the divergence of that field within the volume.
$\oint_{S} F \cdot d S = \int_{V} (\nabla \cdot F) d V$
You now have the skills to compute both sides of this equation. The theorem itself will provide a profound physical insight and a powerful computational shortcut, forming the mathematical heart of Gauss's Law in electrostatics and the continuity equation in fluid dynamics.

Recommended resources