Surface integral

Say we have a surface SS in R3\mathbb R^3 and a scalar field f:R3Rf : \mathbb R^3 \to \mathbb R. We want to integrate ff over the surface SS.

SfdA\int_S f dA

Here dAdA is the differential surface area.

Since we are integrating a 2 dimensional surface in 3 dimensional space, let’s parametrize SS with a function r(u,v):ΩR3\mathbf r(u, v) : \mathbb \Omega \to \mathbb R^3. r\mathbf r takes a point in the region Ω\Omega of the (u,v)(u, v) plane and projects it onto the 3 dimensional surface SS.

Imagine breaking up Ω\Omega into square pieces, each of some area ΔA=Δu×Δv\Delta A = \Delta u \times \Delta v. When the corresponding region is projected onto SS, the area is distorted. If we hold one variable (say, uu) constant and vary the other one by some small dvdv we will get a small vector that we will call rvdv\mathbf r_v dv. Likewise, for a small change dudu we get rudu\mathbf r_u du. These two vectors form a parallelogram in R3\mathbb R^3. The area on SS that corresponds to our differential area du×dvdu \times dv on Ω\Omega is then the area of the parallelogram formed by the two vectors rudu\mathbf r_u du and rvdv\mathbf r_v dv. Recall that the magnitude of the cross product of two vectors is the area of the parallelogram they form, so we can write this are as dA=ru×rvdudvdA = \lVert \mathbf r_u \times \mathbf r_v \rVert du dv.

Note

fxf_x is one way to write the partial derivative of the function ff with respect to xx. It means the same thing as fx\frac{\partial f}{\partial x}.

We can now write our integral in terms of Ω\Omega and r\mathbf r:

Ωf(r(u,v))ru×rvdudv\iint_\Omega f(\mathbf r(u, v)) \lVert \mathbf r_u \times \mathbf r_v \rVert du dv

Over a vector field

The integral we just discussed was over a scalar field ff, but we can just as easily compute one for F:R3R3\mathbf F: \R^3 \to \R^3. In this case, we want to integrate the dot product of F\mathbf F with the outwards pointing normal vector n^\hat{\mathbf n}.

SFn^ds\int_S \mathbf F \cdot \hat{\mathbf n} ds

Note that since n^\hat{\mathbf n} is normal to the surface, it must be perpendicular to both ru\mathbf r_u and rv\mathbf r_v. Therefore the direction of n^\hat{\mathbf n} is just the cross product ru×rv\mathbf r_u \times \mathbf r_v. Since n^\hat{\mathbf n} is a unit vector, we need to divide by the norm, finally giving us n^=ru×rvru×rv\hat{\mathbf n} = \frac{\mathbf r_u \times \mathbf r_v}{\lVert \mathbf r_u \times \mathbf r_v \rVert}.

Now we can set f=Fn^f=\mathbf F \cdot \hat{\mathbf n} and use the surface integral formula above.

ΩF(r(u,v))(ru×rv)dudv\iint_\Omega \mathbf F(\mathbf r(u,v)) \cdot (\mathbf r_u \times \mathbf r_v) du dv
Note

We almost always (unless otherwise noted) want the normal vector on a closed surface to point out. If your normal is pointing the wrong way you can either change the order of the cross product or just add a negative.