◀ Back to Note

Tony's Intertwingled Memex • Reference Manual

Triple Integrals over a Rectangular Region

Formulas & Definitions
Core Definition

A triple integral over a rectangular box B=[a,b]×[c,d]×[r,s]⊂R3B = [a, b] \times [c, d] \times [r, s] \subset \mathbb{R}^3B=[a,b]×[c,d]×[r,s]⊂R3 is the formal three-dimensional limit of a triple Riemann sum:

∭Bf(x,y,z) dV=lim⁡l→∞m→∞n→∞∑i=1l∑j=1m∑k=1nf(xijk∗,yijk∗,zijk∗) ΔVijk\iiint_{B} f(x, y, z) \, dV = \lim_{\substack{l \to \infty \\ m \to \infty \\ n \to \infty}} \sum_{i=1}^{l} \sum_{j=1}^{m} \sum_{k=1}^{n} f(x_{ijk}^*, y_{ijk}^*, z_{ijk}^*) \, \Delta V_{ijk}∭B​f(x,y,z)dV=l→∞m→∞n→∞​lim​i=1∑l​j=1∑m​k=1∑n​f(xijk∗​,yijk∗​,zijk∗​)ΔVijk​
where ΔVijk=Δxi Δyj Δzk\Delta V_{ijk} = \Delta x_i \, \Delta y_j \, \Delta z_kΔVijk​=Δxi​Δyj​Δzk​ is the volume of sub-box BijkB_{ijk}Bijk​, and (xijk∗,yijk∗,zijk∗)(x_{ijk}^*, y_{ijk}^*, z_{ijk}^*)(xijk∗​,yijk∗​,zijk∗​) is any sample point in BijkB_{ijk}Bijk​.

Concept Explanation

A triple integral extends single and double integration into three-dimensional space, accumulating a continuous scalar property f(x,y,z)f(x, y, z)f(x,y,z) over a spatial volume. When the domain is a rectangular box (or rectangular parallelepiped) BBB, all boundary limits for xxx, yyy, and zzz are fixed constants:

a≤x≤b,c≤y≤d,r≤z≤sa \le x \le b, \quad c \le y \le d, \quad r \le z \le sa≤x≤b,c≤y≤d,r≤z≤s

The differential volume element dVdVdV represents an infinitesimal rectangular cuboid with side lengths dxdxdx, dydydy, and dzdzdz:

dV=dz dy dx=dx dy dz=…dV = dz \, dy \, dx = dx \, dy \, dz = \dotsdV=dzdydx=dxdydz=…

Fubini's Theorem for Triple Integrals

If f(x,y,z)f(x, y, z)f(x,y,z) is continuous on the rectangular box BBB, Fubini's Theorem guarantees that the triple integral is equal to the iterated integral evaluated in any of the 3!=63! = 63!=6 possible variable permutations:

∭Bf(x,y,z) dV=∫ab∫cd∫rsf(x,y,z) dz dy dx\iiint_{B} f(x, y, z) \, dV = \int_{a}^{b} \int_{c}^{d} \int_{r}^{s} f(x, y, z) \, dz \, dy \, dx∭B​f(x,y,z)dV=∫ab​∫cd​∫rs​f(x,y,z)dzdydx
Because all limits are constants, we evaluate the iterated integral sequentially from the inside layer out, treating all outer variables as constants during each integration pass.

Algebraic Separability

When the integrand factors completely into independent single-variable functions:

f(x,y,z)=g(x)⋅h(y)⋅k(z)f(x, y, z) = g(x) \cdot h(y) \cdot k(z)f(x,y,z)=g(x)⋅h(y)⋅k(z)
the triple integral over a rectangular box factors directly into the scalar product of three single-variable definite integrals:
∭Bg(x)h(y)k(z) dV=(∫abg(x) dx)(∫cdh(y) dy)(∫rsk(z) dz)\iiint_{B} g(x)h(y)k(z) \, dV = \left( \int_{a}^{b} g(x) \, dx \right) \left( \int_{c}^{d} h(y) \, dy \right) \left( \int_{r}^{s} k(z) \, dz \right)∭B​g(x)h(y)k(z)dV=(∫ab​g(x)dx)(∫cd​h(y)dy)(∫rs​k(z)dz)

Geometric and Physical Interpretations

  1. Geometric Volume: If f(x,y,z)=1f(x, y, z) = 1f(x,y,z)=1, the integral evaluates the spatial volume of the box:
    V(B)=∭B1 dV=(b−a)(d−c)(s−r)V(B) = \iiint_{B} 1 \, dV = (b - a)(d - c)(s - r)V(B)=∭B​1dV=(b−a)(d−c)(s−r)
  2. Average Value: The average value f‾\overline{f}f​ of a continuous function across the box is:
    f‾=1V(B)∭Bf(x,y,z) dV\overline{f} = \frac{1}{V(B)} \iiint_{B} f(x, y, z) \, dVf​=V(B)1​∭B​f(x,y,z)dV
  3. Total Mass and Enclosed Charge: If ρ(x,y,z)\rho(x, y, z)ρ(x,y,z) represents mass density (in kg/m3\text{kg/m}^3kg/m3) or charge density (in C/m3\text{C/m}^3C/m3), the triple integral computes total mass M=∭Bρ dVM = \iiint_B \rho \, dVM=∭B​ρdV or total charge Q=∭Bρ dVQ = \iiint_B \rho \, dVQ=∭B​ρdV.

Engineering Applications

  • Computer & Mechanical Engineering: Evaluating total volumetric heat generation Q=∭q(x,y,z) dVQ = \iiint q(x, y, z) \, dVQ=∭q(x,y,z)dV and spatial temperature profiles across microchip heat sinks and thermoelectric cooling blocks.
  • Electrical & Semiconductor Engineering: Computing total enclosed electrostatic charge in 3D doped semiconductor substrates and MOS transistor channels governed by Gauss's Law: ∭(∇⋅E⃗) dV=Qencε0\iiint (\nabla \cdot \vec{E}) \, dV = \frac{Q_{\text{enc}}}{\varepsilon_0}∭(∇⋅E)dV=ε0​Qenc​​.
  • Materials Science & Solid Mechanics: Determining total mass and gravitational centers of uniform composite structural blocks, acoustic dampening panels, and additive manufacturing test coupons.

Example Problems

Example 1: Iterated Integration over a 3D Box

Evaluate the triple integral:

I=∭B12xyz2 dVwhere B=[0,1]×[1,2]×[0,2]I = \iiint_{B} 12xyz^2 \, dV \quad \text{where } B = [0, 1] \times [1, 2] \times [0, 2]I=∭B​12xyz2dVwhere B=[0,1]×[1,2]×[0,2]

Solution

Apply Fubini's Theorem for Triple Integrals:

I=∫01∫12∫0212xyz2 dz dy dxI = \int_{0}^{1} \int_{1}^{2} \int_{0}^{2} 12xyz^2 \, dz \, dy \, dxI=∫01​∫12​∫02​12xyz2dzdydx

Step 1: Integrate innermost layer with respect to zzz (holding x,yx, yx,y constant):

∫0212xyz2 dz=[12xyz33]02=[4xyz3]02=4xy(2)3−0=32xy\int_{0}^{2} 12xyz^2 \, dz = \left[ 12xy \frac{z^3}{3} \right]_{0}^{2} = \left[ 4xyz^3 \right]_{0}^{2} = 4xy(2)^3 - 0 = 32xy∫02​12xyz2dz=[12xy3z3​]02​=[4xyz3]02​=4xy(2)3−0=32xy

Step 2: Integrate middle layer with respect to yyy (holding xxx constant):

∫1232xy dy=[16xy2]12=16x(2)2−16x(1)2=64x−16x=48x\int_{1}^{2} 32xy \, dy = \left[ 16xy^2 \right]_{1}^{2} = 16x(2)^2 - 16x(1)^2 = 64x - 16x = 48x∫12​32xydy=[16xy2]12​=16x(2)2−16x(1)2=64x−16x=48x

Step 3: Integrate outermost layer with respect to xxx:

∫0148x dx=[24x2]01=24(1)2−0=24\int_{0}^{1} 48x \, dx = \left[ 24x^2 \right]_{0}^{1} = 24(1)^2 - 0 = 24∫01​48xdx=[24x2]01​=24(1)2−0=24

Thus, I=24I = 24I=24.

Example 2: Average Value of a Spatial Function

Find the average value f‾\overline{f}f​ of the scalar field f(x,y,z)=x+2y−zf(x, y, z) = x + 2y - zf(x,y,z)=x+2y−z over the box B=[0,2]×[0,3]×[1,2]B = [0, 2] \times [0, 3] \times [1, 2]B=[0,2]×[0,3]×[1,2].

Solution

Step 1: Compute box volume V(B)V(B)V(B):

V(B)=(2−0)(3−0)(2−1)=2×3×1=6V(B) = (2 - 0)(3 - 0)(2 - 1) = 2 \times 3 \times 1 = 6V(B)=(2−0)(3−0)(2−1)=2×3×1=6

Step 2: Formulate and evaluate the triple integral:

I=∫02∫03∫12(x+2y−z) dz dy dxI = \int_{0}^{2} \int_{0}^{3} \int_{1}^{2} (x + 2y - z) \, dz \, dy \, dxI=∫02​∫03​∫12​(x+2y−z)dzdydx

Inner zzz-integral:

∫12(x+2y−z) dz=[(x+2y)z−z22]12=(2(x+2y)−2)−(1(x+2y)−12)=x+2y−32\int_{1}^{2} (x + 2y - z) \, dz = \left[ (x + 2y)z - \frac{z^2}{2} \right]_{1}^{2} = \left(2(x + 2y) - 2\right) - \left(1(x + 2y) - \frac{1}{2}\right) = x + 2y - \frac{3}{2}∫12​(x+2y−z)dz=[(x+2y)z−2z2​]12​=(2(x+2y)−2)−(1(x+2y)−21​)=x+2y−23​

Middle yyy-integral:

∫03(x+2y−32) dy=[xy+y2−32y]03=3x+9−92=3x+92\int_{0}^{3} \left(x + 2y - \frac{3}{2}\right) \, dy = \left[ xy + y^2 - \frac{3}{2}y \right]_{0}^{3} = 3x + 9 - \frac{9}{2} = 3x + \frac{9}{2}∫03​(x+2y−23​)dy=[xy+y2−23​y]03​=3x+9−29​=3x+29​

Outer xxx-integral:

∫02(3x+92) dx=[32x2+92x]02=32(4)+92(2)=6+9=15\int_{0}^{2} \left(3x + \frac{9}{2}\right) \, dx = \left[ \frac{3}{2}x^2 + \frac{9}{2}x \right]_{0}^{2} = \frac{3}{2}(4) + \frac{9}{2}(2) = 6 + 9 = 15∫02​(3x+29​)dx=[23​x2+29​x]02​=23​(4)+29​(2)=6+9=15

Step 3: Calculate the average value:

f‾=IV(B)=156=2.5\overline{f} = \frac{I}{V(B)} = \frac{15}{6} = 2.5f​=V(B)I​=615​=2.5


Connections & References

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