Calculating Flux with Surface Integrals

Hello! Welcome to your next lesson in the "Mathematical Methods for Physics" module.

In our previous lesson, we learned how to integrate a vector field along a one-dimensional path to calculate quantities like the work done by a force, W = ∫_C **F** ⋅ d**r**. We saw that this involved summing up the component of the field that is tangent to the path.

Today, we extend this idea of integrating fields into a higher dimension. Instead of a 1D path, we will consider a 2D surface living in 3D space. Our goal is to evaluate surface integrals of vector fields, a procedure that allows us to answer a fundamental physical question: "How much of a vector field passes through a given surface?" This quantity is known as flux.

The concept of flux is ubiquitous in physics, from the flow of a fluid through a membrane to the electric field passing through a closed surface in Gauss's Law. Mastering its calculation is essential for your journey into electromagnetism, fluid dynamics, and continuum mechanics.

1. The Concept of Flux

Imagine a river flowing, represented by a velocity vector field v. If you dip a net into the river, the amount of water passing through the net per second depends on three things: the speed of the water, the size of the net, and the orientation of the net relative to the flow. If the net is face-on to the current, you get maximum flow. If it's parallel to the current, nothing passes through it.

This "amount of flow through a surface" is what we call flux.

To quantify this, consider a tiny patch of the surface with area dS. It has a unit normal vector **n** that is perpendicular to the surface at that point. The component of the vector field **F** that is perpendicular to the surface is given by the dot product **F** ⋅ **n**. The flux through this tiny patch is then dΦ = (**F** ⋅ **n**) dS.

To find the total flux, we sum up these contributions over the entire surface S:

Flux = Φ = \iint_{S} (F \cdot n) d S

The following videos introduce this core idea.

Flux of a Vector Field Across a Surface // Vector Calculus

First, let's watch a short segment from Dr. Trefor Bazett's 'Flux of a Vector Field Across a Surface' video. It provides an excellent intuitive picture of flux and the role of the normal vector.

Watch the first four minutes of the video, covering: (0:00 - 1:49) The analogy between flux across a 2D curve and flux across a 3D surface. (1:49 - 4:00) The definition of the flux integral, ∬ F⋅n dS, and the importance of choosing an orientation for the surface (i.e., choosing which way the normal vectors point).

Orientation

As the video highlighted, a crucial first step is to define the orientation of the surface by choosing a continuous field of unit normal vectors.

For an open surface (like a sheet), you might choose "upward" or "downward" normals.
For a closed surface that encloses a volume (like a sphere), there is a standard convention:
- Positive orientation: The normal vectors point outward from the enclosed volume.
- Negative orientation: The normal vectors point inward.

Unless specified otherwise, "flux" usually implies flux in the positive (outward) direction for a closed surface.

Calculus III - Surface Integrals of Vector Fields

For a formal definition of oriented and closed surfaces, please read the first section of this document from Paul's Online Math Notes.

Read the introductory section up to the heading 'In order to work with surface integrals...'. This will formally define an oriented surface, a closed surface, and the convention for positive/negative orientation.

2. The Computational Method

The formula ∬_S (**F** ⋅ **n**) dS is conceptually clear but computationally difficult. As with line integrals, the key is to use parameterization to convert the surface integral into a standard double integral over a flat region in a parameter plane (e.g., the uv-plane or xy-plane).

Let's define a vector surface element d**S** = **n** dS. The flux integral is then written more compactly as ∬_S **F** ⋅ d**S**. Our goal is to find a practical expression for d**S**.

For Parametric Surfaces: `r(u,v)`

Suppose our surface S is parameterized by **r**(u,v) = x(u,v) **î** + y(u,v) **ĵ** + z(u,v) **k̂**.

A normal vector to the surface is given by the cross product of the partial derivatives: **r**_u × **r**_v.
The unit normal is **n** = (**r**_u × **r**_v) / ||**r**_u × **r**_v||.
The scalar surface area element is dS = ||**r**_u × **r**_v|| du dv.

Substituting these into the flux definition yields a wonderful simplification:

d S = n d S = \frac{r _{u} \times r _{v}}{∣∣ r _{u} \times r _{v} ∣∣} \cdot ∣∣ r_{u} \times r_{v} ∣∣ d u d v = (r_{u} \times r_{v}) d u d v

The magnitude terms cancel out! This gives us our main computational formula for flux over a parametric surface defined on a domain D in the uv-plane:

Flux = \iint_{D} F (r (u, v)) \cdot (r_{u} \times r_{v}) d u d v

You must always check if **r**_u × **r**_v points in the desired direction for the orientation. If it points the wrong way, you simply use its negative: - (**r**_u × **r**_v).

Surface Integrals with Vector Fields for Flux Calculation — This image provides a great visual summary. The flux integral `∬_S **F**⋅**n** dS` is computed by converting it to a double integral over the parameter domain `D` using the cross product of the partial derivatives of the parameterization `r(u,v)`.

For Explicit Surfaces: `z = g(x,y)`

This is a special case where we can use x and y as the parameters. The parameterization is **r**(x,y) = <x, y, g(x,y)>.
The cross product simplifies to:

r_{x} \times r_{y} = ⟨ - g_{x}, - g_{y}, 1 ⟩

This vector has a positive z-component, so it corresponds to an upward orientation. For a downward orientation, you would take its negative. The flux integral becomes:

Flux = \iint_{D} F (x, y, g (x, y)) \cdot ⟨ - g_{x}, - g_{y}, 1 ⟩ d x d y

where D is the projection of the surface S onto the xy-plane.

3. Worked Examples

Let's see these methods in action. The best way to learn this is to work through examples with standard geometries.

Lec 27: Vector fields in 3D; surface integrals & flux | MIT 18.02 Multivariable Calculus, Fall 2007

The following video from MIT OpenCourseWare is an excellent, rigorous demonstration of calculating flux. We will focus on an example that requires the full machinery of parameterization and integration.

Please watch the section from 23:28 to 33:52. This part calculates the flux of the vector field F = zk̂ through a sphere of radius 'a'. Pay close attention to: How the sphere is parameterized using spherical coordinates (with angles φ and θ as parameters). How the surface element dS is derived for a sphere. How the vector field **F and the dot product **F**⋅**n are expressed in terms of the parameters. The final setup and evaluation of the double integral.

Test your understanding!

In the video example you just watched (**F** = z**k̂** through a sphere), the professor calculated the flux through the entire sphere. What would the flux be if you only considered the upper hemisphere (z ≥ 0)?

Hint: The parameterization and integrand are the same, but what are the new bounds for the integration?

Show answer

The setup is (1/a) ∫∫ z² dS. In spherical parameters, this became:
$Φ = a^{3} \int_{0}^{2 π} \int_{0}^{π} cos^{2} ϕ sin ϕ d ϕ d θ$
For the upper hemisphere, the angle φ (the polar angle from the z-axis) only goes from 0 (North Pole) to π/2 (equator). The bounds for θ remain 0 to 2π.

Φ_{upper} = a^{3} \int_{0}^{2 π} [\int_{0}^{π /2} cos^{2} ϕ sin ϕ d ϕ] d θ

The inner integral (with u = cos(φ), du = -sin(φ)dφ) evaluates to:

\int_{1}^{0} - u^{2} d u = \int_{0}^{1} u^{2} d u = [\frac{u ^{3}}{3}]_{0}^{1} = \frac{1}{3}

So, the total flux is:

Φ_{upper} = a^{3} \int_{0}^{2 π} \frac{1}{3} d θ = a^{3} \cdot \frac{1}{3} \cdot 2 π = \frac{2 π a ^{3}}{3}

A More Challenging Problem: A Closed Surface

Now let's tackle a problem at the level you're aiming for. It involves a closed surface composed of two different parts, requiring us to evaluate the flux through each and add them together, paying careful attention to the orientation of each part.

Calculus III - Surface Integrals of Vector Fields

The resource from Paul's Online Math Notes contains a superb example that matches your desire for rigorous, multi-step problems. We will work through 'Example 1'.

Please read and follow the solution for Example 1. The problem is to find the flux of F = yĵ - zk̂ through the closed surface S formed by the paraboloid y = x² + z² (for 0 ≤ y ≤ 1) and the disk x² + z² ≤ 1 at y=1. The orientation is positive (outward). Focus on how the problem is broken down: Surface S₁ (the paraboloid): Notice it's an explicit surface of the form y = g(x,z). Pay attention to how the normal vector is derived and why a negative sign is chosen to ensure it points 'outward' (in the negative y-direction for the paraboloid). Surface S₂ (the disk): This is a simple plane. See how the normal vector is just ĵ and the integral simplifies. The final result is the sum of the flux through S₁ and S₂.

This example is an excellent template for how to approach flux integrals over composite surfaces, a common scenario in physics problems (e.g., calculating flux through a cylinder with its top and bottom caps).

4. Physical Interpretation and Application

As mentioned, flux has direct physical meaning. The MIT video you watched briefly explains its connection to fluid flow. Let's formalize that.

If **v** is the velocity field of a fluid, then ∬_S **v** ⋅ d**S** represents the volume of fluid flowing through the surface S per unit time (volume flow rate).
If the field is **F** = ρ**v** where ρ is the fluid's density, then ∬_S **F** ⋅ d**S** is the mass of fluid flowing through S per unit time (mass flow rate).
In electromagnetism, if **E** is the electric field, then Φ_E = ∬_S **E** ⋅ d**S** is the electric flux. Gauss's Law states that the net electric flux through any closed surface is proportional to the total electric charge enclosed by that surface. This is one of the four fundamental Maxwell's Equations, and you now have the mathematical tool to compute one side of it!

Conclusion

Today we have generalized integration from 1D paths to 2D surfaces. This has unlocked the powerful concept of flux, a cornerstone of vector calculus with deep roots in physical law.

Key Takeaways:

Flux measures the net flow of a vector field **F** through an oriented surface S. It is defined as Φ = ∬_S **F** ⋅ **n** dS.
Orientation (the choice of normal vector direction) is critical and must be established before computing. For closed surfaces, "positive" orientation means outward-pointing normals.
We compute flux by converting the surface integral into a regular double integral using a parameterization **r**(u,v). The key computational formula is Φ = ∬_D **F** ⋅ (**r**_u × **r**_v) du dv.
For surfaces given explicitly as z=g(x,y), this simplifies to Φ = ∬_D **F** ⋅ <-g_x, -g_y, 1> dx dy (for upward orientation).

Preview of the Next Lesson:
We have now covered differentiation (grad, div, curl) and integration (line and surface integrals) of fields. The next step is to integrate over a 3D volume. After that, we will learn the two grand theorems of vector calculus: the Divergence Theorem, which relates the flux through a closed surface to the integral of the divergence over the volume it encloses, and Stokes' Theorem, which relates the work done around a closed loop to the flux of the curl through the surface bounded by that loop. You will see how all these concepts beautifully connect.

Recommended resources