Core Definition
For a continuous function of two variables f(x,y) defined over a closed rectangular domain R=[a,b]Γ[c,d]={(x,y)βR2β£aβ€xβ€b,cβ€yβ€d}, the double integral β¬Rβf(x,y)dA evaluates the limit of Riemann sums over planar grid partitions:
β¬Rβf(x,y)dA=Ξxβ0Ξyβ0βlimβi=1βmβj=1βnβf(xijββ,yijββ)ΞxΞy
By
Fubini's Theorem, the double integral equals the iterated single integrals evaluated in either order:
β¬Rβf(x,y)dA=β«cdββ«abβf(x,y)dxdy=β«abββ«cdβf(x,y)dydx
Concept Explanation
When integrating over a rectangle in the xy-plane, all four boundaries are fixed geometric constants (x=a,x=b,y=c,y=d). This geometric symmetry provides unique analytical simplifications:
1. Fubini's Theorem (Iterated Integrals)
If f(x,y) is continuous on the rectangle R=[a,b]Γ[c,d], the double integral can be evaluated as two successive single-variable definite integrals:
- Integrating dydx (Vertical Slices First): Hold x constant and evaluate the inner integral β«cdβf(x,y)dy to produce a single-variable function A(x) representing the cross-sectional area of the solid at position x. Then evaluate β«abβA(x)dx.
- Integrating dxdy (Horizontal Slices First): Hold y constant and evaluate the inner integral β«abβf(x,y)dx to produce A(y), then evaluate β«cdβA(y)dy.
Because the boundaries are constants, reversing the order of integration simply requires exchanging the differential order (dxdyβdydx) and swapping the corresponding constant limits of integration without altering boundary equations.
2. Product Separability Property
If the multivariable integrand factors cleanly into the product of a single-variable function of x and a single-variable function of y, such that f(x,y)=g(x)β
h(y), then over a rectangular domain R=[a,b]Γ[c,d]:
β¬Rβg(x)h(y)dA=(β«abβg(x)dx)(β«cdβh(y)dy)
This decouples the multivariable calculation into the product of two independent single-variable integrals, significantly accelerating evaluation.
3. Spatial Average Value Formula
In single-variable calculus, the average value of f(x) over [a,b] is favgβ=bβa1ββ«abβf(x)dx. For a two-variable function f(x,y) over a region R, the spatial average value is defined by dividing the total accumulated double integral by the planar surface area of the domain:
favgβ=Area(R)1ββ¬Rβf(x,y)dA
For a rectangular domain R=[a,b]Γ[c,d], the area is simply Area(R)=(bβa)(dβc).
Engineering Applications
- Optoelectronics & Display Engineering: Evaluating average surface irradiance Eavgβ and total emitted radiant luminous flux across flat-panel LED displays and solar photovoltaic panels.
- Electrical & Computer Engineering: Calculating total electrostatic charge accumulation Q=β¬RβΟ(x,y)dA across rectangular microchip dies and capacitive sensor pads.
- Mechanical & Thermal Engineering: Determining total heat transfer rate q=β¬Rβh(x,y)(TsββTββ)dA across rectangular radiator cooling fins and heat sinks.
- Civil & Structural Engineering: Computing total contact load and average soil bearing pressure beneath rectangular concrete foundation footings.
Example Problems
Example 1: Basic Iterated Integration & Spatial Average
Evaluate β¬Rβ(3x2β2y)dA over the rectangular region R=[0,2]Γ[1,3], and determine the spatial average value favgβ.
Solution
Step 1: State Governing Formula
State Fubini's Theorem for Rectangular Regions:
β¬Rβf(x,y)dA=β«cdββ«abβf(x,y)dxdy
Here
a=0,b=2 and
c=1,d=3. Set up the iterated integral:
I=β«13ββ«02β(3x2β2y)dxdy
Step 2: Evaluate the Inner Integral (dx)
Treat y as a constant:
β«02β(3x2β2y)dx=[x3β2yx]x=0x=2β=((2)3β2y(2))β0=8β4y
Step 3: Evaluate the Outer Integral (dy)
β«13β(8β4y)dy=[8yβ2y2]13β=(8(3)β2(3)2)β(8(1)β2(1)2)
=(24β18)β(8β2)=6β6=0
Step 4: Compute Spatial Average Value
State the Spatial Average Value Formula:
favgβ=Area(R)1ββ¬Rβf(x,y)dA
The area of the rectangle is
Area(R)=(2β0)(3β1)=2Γ2=4.
favgβ=41β(0)=0
Example 2: Strategic Order Selection via Fubini's Theorem
Evaluate the double integral:
β¬RβxexydAwhereΒ R=[0,1]Γ[0,2]
Solution
Step 1: Analyze Integration Order Strategy
If we integrate with respect to x first, the factor x in xexy forces integration by parts. However, if we integrate with respect to y first, x is treated as a constant multiplier that perfectly matches the internal derivative of the exponent u=xyβΉdu=xdy. Therefore, we choose order dydx.
Step 2: Set up the Iterated Integral
State Fubini's Theorem for Rectangular Regions:
β¬RβxexydA=β«01ββ«02βxexydydx
Step 3: Evaluate Inner Integral (dy)
Treating x as constant:
β«02βxexydy=[exy]y=0y=2β=e2xβe0=e2xβ1
Step 4: Evaluate Outer Integral (dx)
β«01β(e2xβ1)dx=[21βe2xβx]01β=(21βe2β1)β(21βe0β0)
=21βe2β1β21β=21βe2β23β
Connections & References