◀ Back to Note

Tony's Intertwingled Memex • Reference Manual

Triple Integrals over a General Region

Formulas & Definitions
Core Definition

A triple integral over a general 3D bounded region E⊂R3E \subset \mathbb{R}^3E⊂R3 is defined by enclosing EEE within a larger rectangular box BBB and setting an indicator integrand F(x,y,z)=f(x,y,z)F(x, y, z) = f(x, y, z)F(x,y,z)=f(x,y,z) on EEE and 000 outside EEE:

∭Ef(x,y,z) dV=∭BF(x,y,z) dV\iiint_{E} f(x, y, z) \, dV = \iiint_{B} F(x, y, z) \, dV∭E​f(x,y,z)dV=∭B​F(x,y,z)dV
For a Type I region, the solid is bounded between continuous surfaces z=u1(x,y)z = u_1(x, y)z=u1​(x,y) and z=u2(x,y)z = u_2(x, y)z=u2​(x,y) over a planar 2D projection domain D⊂R2D \subset \mathbb{R}^2D⊂R2 in the xyxyxy-plane:
∭Ef(x,y,z) dV=∬D(∫u1(x,y)u2(x,y)f(x,y,z) dz)dA\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∭E​f(x,y,z)dV=∬D​(∫u1​(x,y)u2​(x,y)​f(x,y,z)dz)dA

Concept Explanation

In real physical and engineering problems, 3D solids are rarely simple rectangular boxes. Instead, boundaries are curved surfaces, slanted planes, and intersecting manifolds. We evaluate triple integrals over such general regions by setting functional limits on the inner variables while maintaining constant limits on the outermost variable.

Region Classifications in 3D Space

Depending on which coordinate axis connects the bottom and top bounding surfaces, 3D regions are classified into three fundamental types:

  1. Type I (zzz-simple):
    E={(x,y,z)∣(x,y)∈Dxy, u1(x,y)≤z≤u2(x,y)}E = \{(x, y, z) \mid (x, y) \in D_{xy}, \, u_1(x, y) \le z \le u_2(x, y)\}E={(x,y,z)∣(x,y)∈Dxy​,u1​(x,y)≤z≤u2​(x,y)}
    The solid is bounded below by z=u1(x,y)z = u_1(x, y)z=u1​(x,y) and above by z=u2(x,y)z = u_2(x, y)z=u2​(x,y). The 2D shadow DxyD_{xy}Dxy​ lies in the xyxyxy-plane.
  2. Type II (xxx-simple):
    E={(x,y,z)∣(y,z)∈Dyz, u1(y,z)≤x≤u2(y,z)}E = \{(x, y, z) \mid (y, z) \in D_{yz}, \, u_1(y, z) \le x \le u_2(y, z)\}E={(x,y,z)∣(y,z)∈Dyz​,u1​(y,z)≤x≤u2​(y,z)}
    The solid is bounded behind by x=u1(y,z)x = u_1(y, z)x=u1​(y,z) and in front by x=u2(y,z)x = u_2(y, z)x=u2​(y,z). The 2D shadow DyzD_{yz}Dyz​ lies in the yzyzyz-plane.
  3. Type III (yyy-simple):
    E={(x,y,z)∣(x,z)∈Dxz, u1(x,z)≤y≤u2(x,z)}E = \{(x, y, z) \mid (x, z) \in D_{xz}, \, u_1(x, z) \le y \le u_2(x, z)\}E={(x,y,z)∣(x,z)∈Dxz​,u1​(x,z)≤y≤u2​(x,z)}
    The solid is bounded on the left by y=u1(x,z)y = u_1(x, z)y=u1​(x,z) and on the right by y=u2(x,z)y = u_2(x, z)y=u2​(x,z). The 2D shadow DxzD_{xz}Dxz​ lies in the xzxzxz-plane.

The 4-Step Boundary Construction Protocol

To set up an iterated triple integral over a general region (e.g., standard order dz dy dxdz \, dy \, dxdzdydx):

  1. Inner Surface Limits (zzz): Look along the chosen inner axis (zzz). Identify the entering surface z=u1(x,y)z = u_1(x, y)z=u1​(x,y) (bottom) and exiting surface z=u2(x,y)z = u_2(x, y)z=u2​(x,y) (top).
  2. 2D Projection Shadow (DDD): Project the outermost boundary of the solid onto the remaining coordinate plane (xyxyxy-plane). This boundary is found by intersecting the upper and lower surfaces (u1(x,y)=u2(x,y)u_1(x, y) = u_2(x, y)u1​(x,y)=u2​(x,y)) or setting bounding planes to coordinate axes (z=0z = 0z=0).
  3. Middle Curve Limits (yyy): Treat the 2D projection DDD as a double integral domain. Cast an arrow along the middle axis (yyy) from the lower boundary curve y=g1(x)y = g_1(x)y=g1​(x) to the upper boundary curve y=g2(x)y = g_2(x)y=g2​(x).
  4. Outer Constant Limits (xxx): Identify the absolute minimum and maximum constant extent of the shadow along the outer axis: a≤x≤ba \le x \le ba≤x≤b.

Putting these together yields the complete iterated integral:

∭Ef(x,y,z) dV=∫ab∫g1(x)g2(x)∫u1(x,y)u2(x,y)f(x,y,z) dz dy dx\iiint_{E} f(x, y, z) \, dV = \int_{a}^{b} \int_{g_1(x)}^{g_2(x)} \int_{u_1(x, y)}^{u_2(x, y)} f(x, y, z) \, dz \, dy \, dx∭E​f(x,y,z)dV=∫ab​∫g1​(x)g2​(x)​∫u1​(x,y)u2​(x,y)​f(x,y,z)dzdydx

Physical Quantities: Volume, Mass, and Center of Mass

  • Spatial Volume: V(E)=∭E1 dVV(E) = \iiint_E 1 \, dVV(E)=∭E​1dV.
  • Total Mass: If ρ(x,y,z)\rho(x, y, z)ρ(x,y,z) is mass density, then M=∭Eρ(x,y,z) dVM = \iiint_E \rho(x, y, z) \, dVM=∭E​ρ(x,y,z)dV.
  • First Moments:
    Myz=∭Exρ dV,Mxz=∭Eyρ dV,Mxy=∭Ezρ dVM_{yz} = \iiint_E x \rho \, dV, \quad M_{xz} = \iiint_E y \rho \, dV, \quad M_{xy} = \iiint_E z \rho \, dVMyz​=∭E​xρdV,Mxz​=∭E​yρdV,Mxy​=∭E​zρdV
  • Center of Mass Coordinates:
    xˉ=MyzM,yˉ=MxzM,zˉ=MxyM\bar{x} = \frac{M_{yz}}{M}, \quad \bar{y} = \frac{M_{xz}}{M}, \quad \bar{z} = \frac{M_{xy}}{M}xˉ=MMyz​​,yˉ​=MMxz​​,zˉ=MMxy​​

Engineering Applications

  • Aerospace & Mechanical Engineering: Calculating the total fuel mass and dynamic center of gravity r⃗CG=⟨xˉ,yˉ,zˉ⟩\vec{r}_{\text{CG}} = \langle \bar{x}, \bar{y}, \bar{z} \ranglerCG​=⟨xˉ,yˉ​,zˉ⟩ for non-rectangular wingtip fuel cells, fuselage tanks, and liquid propellant slosh baffles.
  • Civil & Structural Engineering: Computing total concrete volume, dead load, and foundation bearing moments for asymmetric bridge abutments, gravity retaining walls, and dam spillways.
  • Ocean & Naval Architecture: Determining hydrostatic displacement volume and center of buoyancy for ship hulls and submarine ballast tanks bounded by parabolic waterlines and inclined transverse bulkheads.

Example Problems

Example 1: Integrating over a First-Octant Tetrahedron

Evaluate ∭Ez dV\iiint_{E} z \, dV∭E​zdV, where EEE is the solid tetrahedron in the first octant bounded by the planes x=0x = 0x=0, y=0y = 0y=0, z=0z = 0z=0, and 2x+y+z=22x + y + z = 22x+y+z=2.

Solution

Step 1: Determine inner zzz-limits:
The solid lies above z=0z = 0z=0 and below the inclined plane 2x+y+z=2  ⟹  z=2−2x−y2x + y + z = 2 \implies z = 2 - 2x - y2x+y+z=2⟹z=2−2x−y.
Thus, 0≤z≤2−2x−y0 \le z \le 2 - 2x - y0≤z≤2−2x−y.

Step 2: Determine 2D projection DDD onto the xyxyxy-plane:
Setting z=0z = 0z=0 in the plane equation gives 2x+y=2  ⟹  y=2−2x2x + y = 2 \implies y = 2 - 2x2x+y=2⟹y=2−2x. In the first octant, x≥0,y≥0x \ge 0, y \ge 0x≥0,y≥0.
Middle limits: 0≤y≤2−2x0 \le y \le 2 - 2x0≤y≤2−2x.
Outer limits: Setting y=0y = 0y=0 yields 2x=2  ⟹  0≤x≤12x = 2 \implies 0 \le x \le 12x=2⟹0≤x≤1.

Step 3: Set up and evaluate the iterated integral:

I=∫01∫02−2x∫02−2x−yz dz dy dxI = \int_{0}^{1} \int_{0}^{2-2x} \int_{0}^{2-2x-y} z \, dz \, dy \, dxI=∫01​∫02−2x​∫02−2x−y​zdzdydx

Inner zzz-integral:

∫02−2x−yz dz=[z22]02−2x−y=12(2−2x−y)2\int_{0}^{2-2x-y} z \, dz = \left[ \frac{z^2}{2} \right]_{0}^{2-2x-y} = \frac{1}{2}(2 - 2x - y)^2∫02−2x−y​zdz=[2z2​]02−2x−y​=21​(2−2x−y)2

Middle yyy-integral:
Use uuu-substitution: let u=2−2x−y  ⟹  dy=−duu = 2 - 2x - y \implies dy = -duu=2−2x−y⟹dy=−du.
When y=0,u=2−2xy = 0, u = 2 - 2xy=0,u=2−2x; when y=2−2x,u=0y = 2 - 2x, u = 0y=2−2x,u=0.

∫2−2x012u2(−du)=∫02−2x12u2 du=[u36]02−2x=(2−2x)36=8(1−x)36=43(1−x)3\int_{2-2x}^{0} \frac{1}{2}u^2 (-du) = \int_{0}^{2-2x} \frac{1}{2}u^2 \, du = \left[ \frac{u^3}{6} \right]_{0}^{2-2x} = \frac{(2 - 2x)^3}{6} = \frac{8(1 - x)^3}{6} = \frac{4}{3}(1 - x)^3∫2−2x0​21​u2(−du)=∫02−2x​21​u2du=[6u3​]02−2x​=6(2−2x)3​=68(1−x)3​=34​(1−x)3

Outer xxx-integral:

∫0143(1−x)3 dx=[−13⋅43(1−x)4]01=[−49(1−x)4]01=0−(−49)=49\int_{0}^{1} \frac{4}{3}(1 - x)^3 \, dx = \left[ -\frac{1}{3} \cdot \frac{4}{3}(1 - x)^4 \right]_{0}^{1} = \left[ -\frac{4}{9}(1 - x)^4 \right]_{0}^{1} = 0 - \left(-\frac{4}{9}\right) = \frac{4}{9}∫01​34​(1−x)3dx=[−31​⋅34​(1−x)4]01​=[−94​(1−x)4]01​=0−(−94​)=94​

Thus, ∭Ez dV=49\iiint_E z \, dV = \frac{4}{9}∭E​zdV=94​.

Example 2: Volume Bounded by a Parabolic Cylinder

Find the volume of the solid EEE bounded by the parabolic cylinder y=x2y = x^2y=x2 and the planes z=0z = 0z=0 and y+z=4y + z = 4y+z=4.

Solution

Step 1: Inner zzz-limits:
The solid extends from the floor z=0z = 0z=0 up to the ceiling y+z=4  ⟹  z=4−yy + z = 4 \implies z = 4 - yy+z=4⟹z=4−y.
Thus, 0≤z≤4−y0 \le z \le 4 - y0≤z≤4−y.

Step 2: 2D projection DDD onto the xyxyxy-plane:
Setting z=0z = 0z=0 in y+z=4y + z = 4y+z=4 gives y=4y = 4y=4. The lateral boundary is y=x2y = x^2y=x2.
Middle limits (yyy): for each xxx, yyy goes from the parabola up to the line: x2≤y≤4x^2 \le y \le 4x2≤y≤4.
Outer limits (xxx): the intersection points are x2=4  ⟹  −2≤x≤2x^2 = 4 \implies -2 \le x \le 2x2=4⟹−2≤x≤2.

Step 3: Evaluate volume integral:

V=∫−22∫x24∫04−y1 dz dy dx=∫−22∫x24(4−y) dy dxV = \int_{-2}^{2} \int_{x^2}^{4} \int_{0}^{4-y} 1 \, dz \, dy \, dx = \int_{-2}^{2} \int_{x^2}^{4} (4 - y) \, dy \, dxV=∫−22​∫x24​∫04−y​1dzdydx=∫−22​∫x24​(4−y)dydx

Middle yyy-integral:

∫x24(4−y) dy=[4y−y22]x24=(16−8)−(4x2−x42)=8−4x2+x42\int_{x^2}^{4} (4 - y) \, dy = \left[ 4y - \frac{y^2}{2} \right]_{x^2}^{4} = \left(16 - 8\right) - \left(4x^2 - \frac{x^4}{2}\right) = 8 - 4x^2 + \frac{x^4}{2}∫x24​(4−y)dy=[4y−2y2​]x24​=(16−8)−(4x2−2x4​)=8−4x2+2x4​

Outer xxx-integral (exploiting symmetry of even integrand):

V=2∫02(8−4x2+x42) dx=2[8x−43x3+x510]02V = 2 \int_{0}^{2} \left(8 - 4x^2 + \frac{x^4}{2}\right) \, dx = 2 \left[ 8x - \frac{4}{3}x^3 + \frac{x^5}{10} \right]_{0}^{2}V=2∫02​(8−4x2+2x4​)dx=2[8x−34​x3+10x5​]02​
V=2(16−323+3210)=2(16−323+165)=2(240−160+4815)=2(12815)=25615V = 2 \left( 16 - \frac{32}{3} + \frac{32}{10} \right) = 2 \left( 16 - \frac{32}{3} + \frac{16}{5} \right) = 2 \left( \frac{240 - 160 + 48}{15} \right) = 2\left(\frac{128}{15}\right) = \frac{256}{15}V=2(16−332​+1032​)=2(16−332​+516​)=2(15240−160+48​)=2(15128​)=15256​


Connections & References

  • Parent Concepts: Multiple Integrals, Double Integrals over a General Region
  • Sub-concepts & Extensions:
    • Triple Integrals over a Rectangular Region
    • Changing the Order of Integration in Triple Integrals
    • Change of Variables in Multiple Integrals
  • Course Links:
    • Calc-3 Session 16 - Triple Integrals