Triple Integrals in Cylindrical - Calculus 3

SECTION 1

Historical Context & Motivation

The challenge of computing volumes, masses, and other quantities over three-dimensional regions has driven mathematicians for centuries. In Cartesian coordinates, many regions with rotational symmetry—cylinders, cones, paraboloids—require cumbersome algebraic descriptions that make integration prohibitively difficult. The development of cylindrical coordinates arose naturally from the earlier polar coordinate system, extending a two-dimensional idea into three dimensions by appending a vertical axis. This seemingly simple extension unlocked an elegant way to evaluate triple integrals over rotationally symmetric domains, transforming intractable Cartesian integrals into straightforward computations.

1671 • Newton's Method of Fluxions

Isaac Newton develops early integration techniques in his unpublished manuscript, laying the groundwork for computing areas and volumes through infinitesimal sums.

1691 • Polar Coordinates Formalized

Jakob Bernoulli publishes a systematic treatment of polar coordinates, describing curves via a radial distance r and angle θ rather than Cartesian x and y.

1736 • Euler's Mechanica

Leonhard Euler uses extended polar-type coordinates in three dimensions to solve mechanics problems involving rotational symmetry, foreshadowing modern cylindrical coordinates.

1811 • Jacobian Determinant Introduced

Carl Gustav Jacob Jacobi formalizes the determinant that governs how volume elements transform under coordinate changes, providing the rigorous basis for the factor r in cylindrical integrals.

1900s • Modern Multivariable Calculus

Cylindrical and spherical coordinate integration become standard tools in physics and engineering curricula, essential for electromagnetic theory, fluid dynamics, and heat conduction.

The central question this topic addresses is deceptively practical: when a three-dimensional region possesses symmetry about an axis, how can we exploit that symmetry to simplify the computation of triple integrals? As we will see, the answer lies in replacing x, y, and z with the cylindrical variables r, θ, and z, along with a crucial scaling factor that accounts for how volume elements warp under the transformation.

SECTION 2

Core Principles & Definitions

Cylindrical coordinates extend the familiar polar coordinate system into three dimensions. A point in space is described by three quantities: the radial distance r from the z-axis, the azimuthal angle θ measured from the positive x-axis in the xy-plane, and the height z, which remains unchanged from Cartesian coordinates. The conversion formulas x = r cos θ, y = r sin θ, and z = z bridge the two systems, while the inverse relations r = √(x² + y²) and θ = arctan(y/x) allow translation back to Cartesian form.

Coordinate Transformation

Replace (x, y, z) with (r, θ, z) via x = r cos θ, y = r sin θ, z = z. This maps the Cartesian grid to a system built around radial symmetry.

The Jacobian Factor r

The differential volume element transforms as dV = r dr dθ dz. The factor r arises from the Jacobian determinant of the coordinate change and is essential—omitting it is the most common error students make.

Bounds of Integration

In cylindrical coordinates, limits on r, θ, and z describe the region of integration. Cylinders become constant-r surfaces, half-planes become constant-θ surfaces, and horizontal planes remain constant-z surfaces.

When to Use Cylindrical

Choose cylindrical coordinates when the integration region or the integrand involves x² + y², when boundaries are cylinders, cones, or paraboloids, or when the domain projects onto a circle or annulus in the xy-plane.

✦

KEY TAKEAWAY

Think of cylindrical coordinates like unwrapping a can label. In Cartesian coordinates, describing the can's curved surface requires the equation x² + y² = R². In cylindrical coordinates, it's simply r = R—a single constant. The integral over the can's interior becomes a straightforward sweep: r from 0 to R, θ from 0 to 2π, z from bottom to top. The factor r in dV = r dr dθ dz accounts for the fact that annular strips farther from the axis enclose more area, much like how the outer lanes of a running track are longer than the inner ones.

SECTION 3

Visualizing Cylindrical Coordinates

A geometric understanding of cylindrical coordinates is essential before writing any integrals. The diagram below illustrates how a point P in three-dimensional space is located by its radial distance r from the z-axis, the angle θ that its projection onto the xy-plane makes with the positive x-axis, and its height z above the xy-plane. Notice how the coordinate surfaces in this system are qualitatively different from those in Cartesian space: constant-r surfaces are cylinders, constant-θ surfaces are half-planes emanating from the z-axis, and constant-z surfaces are horizontal planes.

A point P in cylindrical coordinates is specified by the radial distance r (cyan) from the z-axis, the angle θ (violet) from the positive x-axis, and the height z (amber). The dashed ellipse represents a constant-r cylinder projected into the xy-plane.

In the diagram, observe that r is always non-negative and measures the perpendicular distance from the z-axis—not from the origin. The angle θ sweeps counterclockwise when viewed from above, and z behaves identically to its Cartesian counterpart. When we integrate over a region described by these coordinates, the differential volume element dV = r dr dθ dz reflects the fact that a small "box" in cylindrical coordinates is not a rectangular parallelepiped but rather a curved wedge whose cross-sectional area increases linearly with r.

SECTION 4

Mathematical Framework

The rigorous foundation for triple integrals in cylindrical coordinates rests on the change-of-variables theorem for multiple integrals. When we transform from Cartesian to cylindrical coordinates, we must multiply the integrand by the absolute value of the Jacobian determinant of the transformation, which ensures that we correctly account for how infinitesimal volume elements stretch or compress under the mapping.

COORDINATE TRANSFORMATION

x = r cos θ, y = r sin θ, z = z

where r ≥ 0, 0 ≤ θ < 2π, and z ∈ ℝ.

JACOBIAN DETERMINANT

∂(x, y, z)/∂(r, θ, z) = | cos θ −r sin θ 0 | = r | sin θ r cos θ 0 | | 0 0 1 |

The determinant evaluates to r cos²θ + r sin²θ = r. Since r ≥ 0, we have |J| = r.

TRIPLE INTEGRAL IN CYLINDRICAL COORDINATES

∭_E f(x, y, z) dV = ∫∫∫ f(r cos θ, r sin θ, z) · r dr dθ dz

The factor r is the Jacobian and must always be included. The limits of integration on r, θ, and z are chosen to describe the region E.

To set up the integral, first identify which variable is most naturally bounded by constants (often θ or z), and integrate with respect to it on the outside. Then express the remaining limits in terms of the outer variables. For many standard regions—solid cylinders, cones, paraboloids capped by planes—the bounds on θ are simply 0 to 2π, the bounds on z are functions of r (or constants), and the bounds on r are constants or functions of θ.

TYPICAL ITERATED FORM

∫₀²π ∫₀ᴿ ∫_{g₁(r)}^{g₂(r)} f(r cos θ, r sin θ, z) · r dz dr dθ

Here R is the maximum radial extent of the region, and g₁(r), g₂(r) are the lower and upper z-bounds as functions of r.

SECTION 5

The Cylindrical Volume Element in Detail

Understanding why the factor r appears in the volume element is not merely algebraic bookkeeping—it has a clear geometric interpretation. In Cartesian coordinates, the infinitesimal volume element dx dy dz is a tiny rectangular box with sides of length dx, dy, and dz. When we switch to cylindrical coordinates, the corresponding element is a small curved wedge. Its dimensions are dr in the radial direction, r dθ in the tangential direction (the arc length at distance r), and dz vertically. Multiplying these three dimensions gives the cylindrical volume element dV = r dr dθ dz. The diagram below shows this wedge-shaped element and labels its three dimensions.

The infinitesimal volume element in cylindrical coordinates is a curved wedge with radial thickness dr, tangential arc length r dθ, and vertical height dz. Their product gives dV = r dr dθ dz.

The geometric picture makes clear why forgetting the factor r is such a consequential error: without it, every annular strip at every radius would be weighted equally, drastically underestimating contributions from regions far from the axis. Points at larger r subtend the same angle dθ but sweep out a proportionally larger arc, so their volume contribution scales linearly with r.

⚠

Common Pitfall

Always include the Jacobian factor r when converting to cylindrical coordinates. The volume element is dV = r dr dθ dz, not dr dθ dz. This is analogous to the factor r that appears in polar double integrals with dA = r dr dθ.

SECTION 6

Worked Example

Let us compute the volume of the solid region E bounded below by the paraboloid z = x² + y² and above by the plane z = 4. This is a classic problem that highlights the power of cylindrical coordinates: the Cartesian description of the paraboloid involves x² + y², which simplifies to r² in cylindrical form.

Volume Under a Plane and Above a Paraboloid

1
Step 1 — Describe the Region

The paraboloid z = x² + y² opens upward and the plane z = 4 caps it from above. They intersect where x² + y² = 4, which is a circle of radius 2 in the xy-plane. In cylindrical coordinates, the paraboloid becomes z = r² and the intersection circle is r = 2.

Region: 0 ≤ θ ≤ 2π, 0 ≤ r ≤ 2, r² ≤ z ≤ 4
2
Step 2 — Write the Integral

The volume is given by ∭_E dV. Converting to cylindrical coordinates with dV = r dz dr dθ, and inserting the limits identified above:

V = ∫₀²π ∫₀² ∫_{r²}⁴ r dz dr dθ
3
Step 3 — Evaluate the Innermost Integral (dz)

Integrating with respect to z first: ∫_{r²}⁴ r dz = r · (4 − r²) = 4r − r³.

∫_{r²}⁴ r dz = 4r − r³
4
Step 4 — Evaluate the Middle Integral (dr)

Now integrate with respect to r from 0 to 2: ∫₀² (4r − r³) dr = [2r² − r⁴/4]₀² = (2 · 4 − 16/4) − 0 = 8 − 4 = 4.

∫₀² (4r − r³) dr = 4
5
Step 5 — Evaluate the Outermost Integral (dθ)

Finally, integrate with respect to θ: ∫₀²π 4 dθ = 4 · 2π = 8π.

V = 8π ≈ 25.13 cubic units

This example demonstrates the general workflow: identify the region in cylindrical coordinates, express all bounds in terms of r, θ, and z, include the Jacobian factor r, then integrate from inside out. The full rotational symmetry of the region around the z-axis allowed θ to range over [0, 2π] with no dependence on r or z, factoring out as a simple multiplicative constant.

SECTION 7

Cylindrical vs. Other Coordinate Systems

Choosing the right coordinate system is one of the most important strategic decisions when setting up a triple integral. The table below compares the three standard coordinate systems for triple integration, highlighting when each is most effective and what form the volume element takes.

Feature	Cartesian (x, y, z)	Cylindrical (r, θ, z)	Spherical (ρ, φ, θ)
Volume element	dx dy dz	r dr dθ dz	ρ² sin φ dρ dφ dθ
Best suited for	Rectangular boxes, regions with flat planar boundaries	Cylinders, cones, paraboloids, regions with axial symmetry	Spheres, hemispheres, cones with vertex at origin
Simplifies x² + y²?	No	Yes → r²	Partially → ρ² sin²φ
Simplifies x² + y² + z²?	No	Partially → r² + z²	Yes → ρ²
Jacobian complexity	1 (trivial)	r (simple)	ρ² sin φ (moderate)

Comparison of standard coordinate systems for triple integration

✦

CHOOSING THE RIGHT SYSTEM

Think of coordinate systems as specialized tools. Cartesian is your general-purpose wrench—it works on everything but excels at nothing special. Cylindrical coordinates are like a pipe wrench—specifically designed for problems with circular cross-sections and axial symmetry. Spherical coordinates are your ball-and-socket tool—ideal when the geometry radiates from a single point. Matching the tool to the geometry of the region can reduce a multi-page computation to a few lines.

SECTION 8

Connections to Advanced Theory

Triple integrals in cylindrical coordinates serve as a gateway to several important areas of advanced mathematics and physics. The Jacobian determinant that gives rise to the factor r is a special case of the general theory of differential forms and the pullback of a volume form under a smooth mapping. In courses on differential geometry or real analysis, you will see that dV = r dr ∧ dθ ∧ dz is the pullback of the Cartesian 3-form dx ∧ dy ∧ dz under the cylindrical coordinate map, with the wedge product ensuring proper orientation.

This Lesson	Advanced Extension
Jacobian determinant \|J\| = r	General change-of-variables via differential forms and pullbacks on manifolds
Triple integral of density → mass	Moments of inertia, center of mass in engineering mechanics and rotational dynamics
Volume beneath z = f(r, θ)	Solutions to Laplace's and Poisson's equations in cylindrical symmetry (electrostatics, heat flow)
Limits 0 ≤ θ ≤ 2π for full rotation	Fourier series in θ; separation of variables in cylindrical PDEs (Bessel functions)

From triple integrals in cylindrical coordinates to advanced topics

In physics, cylindrical coordinates are indispensable for problems involving long wires, pipes, solenoids, and rotating shafts—any geometry where translational symmetry along one axis combines with rotational symmetry about that axis. The triple integrals you set up in this lesson are the same integrals that appear in computing electric fields via Gauss's law, gravitational potentials of cylindrical mass distributions, and the temperature profiles of cylindrical conductors. Mastering the setup and evaluation of these integrals now provides a foundation that will carry through electromagnetism, thermodynamics, and advanced engineering courses.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

Explain why the volume element in cylindrical coordinates is dV = r dr dθ dz rather than simply dr dθ dz. What geometric feature does the factor r account for, and what would go wrong if it were omitted?

Show answer

The factor r arises because the tangential dimension of an infinitesimal volume element at distance r from the z-axis is r dθ (an arc length), not simply dθ. Without r, the element would be dr × dθ × dz, which has units of length × radians × length—not a proper volume. Geometrically, omitting r would treat thin annular strips near the axis the same as strips far from it, even though the latter enclose far more volume. This is the Jacobian determinant of the coordinate transformation from Cartesian to cylindrical, and it corrects for the non-uniform scaling of area as one moves away from the axis.

Question 2: BASIC CALCULATION

Evaluate the triple integral ∭_E z dV where E is the solid cylinder x² + y² ≤ 9, 0 ≤ z ≤ 2.

Show answer

Convert to cylindrical: x² + y² ≤ 9 becomes r ≤ 3, so the limits are 0 ≤ θ ≤ 2π, 0 ≤ r ≤ 3, 0 ≤ z ≤ 2. The integral becomes ∫₀²π ∫₀³ ∫₀² z · r dz dr dθ. Evaluating inner to outer: ∫₀² z dz = [z²/2]₀² = 2; then ∫₀³ r · 2 dr = 2[r²/2]₀³ = 9; then ∫₀²π 9 dθ = 18π. The final answer is 18π.

Question 3: INTERMEDIATE

Find the volume of the solid bounded above by the sphere x² + y² + z² = 16 and below by the cone z = √(x² + y²). Set up and evaluate the integral in cylindrical coordinates.

Show answer

In cylindrical coordinates the cone becomes z = r and the sphere becomes r² + z² = 16, so z = √(16 − r²). They intersect when r = √(16 − r²), giving r² = 8, so r = 2√2. The region has full rotational symmetry: 0 ≤ θ ≤ 2π, 0 ≤ r ≤ 2√2, and for each r the z-limits run from the cone z = r up to the sphere z = √(16 − r²). The volume integral is V = ∫₀²π ∫₀^{2√2} ∫ᵣ^{√(16−r²)} r dz dr dθ. The inner integral gives r(√(16 − r²) − r). Using the substitution u = 16 − r² for the first term and direct integration for the second, one obtains V = 2π[−(1/3)

(16 − r²)^{3/2}

− r³/3]₀^{2√2} = 2π[(−(1/3)

(8)^{3/2}

− (2√2)³/3) − (−64/3)] = 2π · 64(2 − √2)/6 = 64π(2 − √2)/3.

Question 4: APPLIED

A cylindrical tank of radius 3 m and height 5 m is filled with a fluid whose density varies as δ(r, θ, z) = 1000 − 50z kg/m³ (density decreases with height due to thermal stratification). Find the total mass of the fluid.

Show answer

Mass is M = ∭_E δ dV. In cylindrical coordinates: M = ∫₀²π ∫₀³ ∫₀⁵ (1000 − 50z) r dz dr dθ. Innermost: ∫₀⁵ (1000 − 50z) dz = [1000z − 25z²]₀⁵ = 5000 − 625 = 4375. Middle: ∫₀³ 4375 · r dr = 4375[r²/2]₀³ = 4375 · 9/2 = 19687.5. Outer: ∫₀²π 19687.5 dθ = 39375π ≈ 123,700 kg. The total mass is 39375π kg, or approximately 1.237 × 10⁵ kg. This result makes physical sense: it is less than the mass of 1000 kg/m³ × π(9)(5) = 45000π ≈ 141,372 kg that would result from uniform density 1000 throughout, reflecting the reduced density at higher elevations.

Question 5: CRITICAL THINKING

Consider the integral ∫₋₂² ∫₋√(4−x²)^{√(4−x²)} ∫_{x²+y²}^{4} (x² + y²) dz dy dx. Without evaluating in Cartesian coordinates, convert this integral to cylindrical coordinates, explain why each limit transforms the way it does, and then evaluate it.

Show answer

The Cartesian limits describe: x from −2 to 2, y from −√(4 − x²) to √(4 − x²) (a disk of radius 2 in the xy-plane), and z from x² + y² to 4 (between paraboloid z = x² + y² and plane z = 4). In cylindrical form, the disk becomes 0 ≤ r ≤ 2, 0 ≤ θ ≤ 2π, and z = x² + y² becomes z = r², so the z-limits are r² ≤ z ≤ 4. The integrand x² + y² = r². Including the Jacobian: ∫₀²π ∫₀² ∫_{r²}⁴ r² · r dz dr dθ = ∫₀²π ∫₀² ∫_{r²}⁴ r³ dz dr dθ. Inner: r³(4 − r²) = 4r³ − r⁵. Middle: ∫₀² (4r³ − r⁵) dr = [r⁴ − r⁶/6]₀² = 16 − 64/6 = 16 − 32/3 = 16/3. Outer: (16/3)(2π) = 32π/3. The conversion is advantageous because the circular domain and the integrand x² + y² both simplify dramatically, reducing a difficult iterated integral in three Cartesian variables to a routine computation.

SUMMARY

Summary

Triple integrals in cylindrical coordinates provide a powerful tool for evaluating volume, mass, and other quantities over regions with axial symmetry. The system replaces x and y with the polar variables r and θ while retaining z, yielding the conversion formulas x = r cos θ, y = r sin θ, z = z. The critical step in any conversion is including the Jacobian factor r in the volume element dV = r dr dθ dz, which accounts for the geometric stretching of area elements at larger radii.

Setting up the integral requires expressing the bounds of integration for r, θ, and z to describe the region, then integrating from the inside out. Cylindrical coordinates are especially effective for regions bounded by cylinders, cones, and paraboloids, where the expression x² + y² simplifies to r². Mastery of this technique prepares you for spherical coordinates, applications in physics and engineering, and the general theory of coordinate transformations via the Jacobian determinant.