Chapter 27 · Magnetism

Magnetic Flux

Three equivalent forms, five worked examples, zero handwaving.

ΦB = ∫ B⃗ · dA⃗

What Is Magnetic Flux?

Magnetic flux is the answer to one question: how much magnetic field is actually passing through a surface?

Not how strong the field is in general — but how much of it is punching through a specific area you care about. Think of holding a hoop in a wind tunnel. The "wind flux" through the hoop depends on (1) how strong the wind is, (2) how big the hoop is, and (3) which way you're holding the hoop relative to the wind. Point it face-on: maximum throughput. Tilt it sideways: less gets through. Edge-on: zero. Magnetic flux works the same way.

🧲 Physical Picture

Imagine the magnetic field as a bunch of arrows flying through space. Flux counts how many of those arrows actually pierce your surface, weighted by how perpendicularly they hit it. Arrows that graze along the surface at a shallow angle barely count. Arrows that punch straight through count fully.

There are three algebraically equivalent ways to write this:

Form 1 — Angle
\[\Phi_B = \int B\cos\phi\,dA\]

Use \(\phi\), the angle between \(\textbf{B}\) and the surface normal.

Form 2 — Projection
\[\Phi_B = \int B_\perp\,dA\]

Project \(\textbf{B}\) onto the normal first, then integrate.

Form 3 — Dot Product
\[\Phi_B = \int \textbf{B}\cdot d\textbf{A}\]

Most general. The dot product does the projection automatically.

These are all the same because:

\[\textbf{B}\cdot d\textbf{A} = B\,dA\cos\phi = B_\perp\,dA\]

The dot product is defined to extract the parallel (perpendicular-to-surface) component automatically. That's literally what it does geometrically.


Every Symbol, Dissected

Symbol Name Physical Meaning
\(\Phi_B\) Magnetic flux Total magnetic "throughput" through the surface. Units: Wb (weber) = T·m².
\(B\) Field magnitude How strong the magnetic field is at a given point. Units: tesla (T).
\(\phi\) Angle to normal Angle between \(\textbf{B}\) and the outward surface normal — NOT the surface itself. This is the #1 trap.
\(dA\) Area element Infinitesimal patch of the surface. Scalar — just a tiny area magnitude.
\(d\textbf{A}\) Vector area element Same tiny patch, but now it's a vector pointing in the direction of the surface normal. Direction encodes orientation.
\(B_\perp\) Normal component The component of \(\textbf{B}\) that is perpendicular to the surface — i.e., actually piercing it. \(B_\perp = B\cos\phi\).
\(\hat{n}\) Unit normal A unit vector perpendicular to the surface. \(d\textbf{A} = \hat{n}\,dA\).
⚠ Classic Trap — Angle Convention

\(\phi\) is the angle between \(\textbf{B}\) and the surface normal. If the problem gives you the angle between \(\textbf{B}\) and the surface itself, call it \(\theta\). Then \(\phi = 90° - \theta\), and: \[\Phi_B = BA\cos(90°-\theta) = BA\sin\theta\] Swap in \(\sin\) instead of \(\cos\). Mixing these up will give you a completely wrong answer.


Five Examples, Three Forms

1
Uniform Field, Tilted Surface Form 1 — ∫ B cos φ dA

Setup: A uniform magnetic field \(B = 0.50\text{ T}\) passes through a flat rectangular surface of area \(A = 0.20\text{ m}^2\). The angle between \(\textbf{B}\) and the surface normal is \(\phi = 60°\). Find \(\Phi_B\).

🔍 Intuition Before Math

The field is tilted at 60° from the normal. \(\cos 60° = 0.5\), so only half the field's strength is punching perpendicularly through the surface. We expect an answer that's half of the "face-on" maximum.

Step 1 — Recognize Uniformity

Because both \(B\) and \(\phi\) are constant over the entire flat surface, they factor out of the integral:

\[\Phi_B = \int B\cos\phi\,dA = B\cos\phi \int dA\]

The remaining integral \(\int dA\) simply sums up all the tiny area patches, which gives the total area \(A\):

\[\int dA = A \qquad \Longrightarrow \qquad \Phi_B = BA\cos\phi\]

This is the flat-surface, uniform-field shortcut. You only pull out the integral when things are constant — if either \(B\) or \(\phi\) varied across the surface, you'd have to keep the integral.

Step 2 — Evaluate the Cosine
\[\cos 60° = \frac{1}{2} = 0.5\]

This is the geometric projection factor — it tells you that only half of the field's magnitude is contributing to the flux. The other half is sliding along the surface and contributes nothing.

Step 3 — Substitute Numbers
\[\Phi_B = (0.50\text{ T})(0.20\text{ m}^2)(0.5)\]
\[\Phi_B = (0.50)(0.20)(0.5) = 0.050\text{ T·m}^2\]
\[\boxed{\Phi_B = 5.0\times 10^{-2}\text{ Wb}}\]
Result check: The maximum possible flux here would be \(BA = (0.50)(0.20) = 0.10\text{ Wb}\) (field perpendicular to surface). We got exactly half of that — consistent with \(\cos 60° = 0.5\). ✓
2
Same Setup, Perpendicular Component View Form 2 — ∫ B⊥ dA

Setup: Same problem as Example 1 — \(B = 0.50\text{ T}\), \(A = 0.20\text{ m}^2\), \(\phi = 60°\) — but this time we use Form 2 to build different geometric intuition.

🔍 Different Framing, Same Physics

Instead of keeping the full field \(B\) and multiplying by \(\cos\phi\) at the end, Form 2 says: decompose \(\textbf{B}\) into components first, then throw away the parallel component entirely. Only the perpendicular piece — the part actually stabbing through the surface — matters.

Step 1 — Decompose the Field

The magnetic field \(\textbf{B}\) can be split into two components relative to the surface:

\[B_\perp = B\cos\phi \qquad \text{(perpendicular to surface — contributes to flux)}\] \[B_\parallel = B\sin\phi \qquad \text{(parallel to surface — contributes nothing)}\]

Plug in numbers:

\[B_\perp = (0.50)\cos 60° = (0.50)(0.5) = 0.25\text{ T}\]
Step 2 — Integrate Only B⊥

Since \(B_\perp\) is constant over the surface:

\[\Phi_B = \int B_\perp\,dA = B_\perp \int dA = B_\perp A\]
\[\Phi_B = (0.25\text{ T})(0.20\text{ m}^2)\]
\[\boxed{\Phi_B = 5.0\times 10^{-2}\text{ Wb}}\]
Why is this the same? Because \(B_\perp A = B\cos\phi \cdot A\). Forms 1 and 2 are just rearranged versions of the same algebra. Form 2 is sometimes more intuitive because you're explicitly discarding the useless parallel component upfront.
3
Component Form and the Dot Product Form 3 — ∫ B⃗ · dA⃗

Setup: The same physical scenario, but now the field is given in component form:

\[\textbf{B} = 0.433\,\hat{i} + 0.25\,\hat{j}\ \text{ T}\]

The surface lies in the \(xz\)-plane, so its outward normal points in the \(\hat{j}\) direction. Therefore:

\[d\textbf{A} = dA\,\hat{j} \qquad \Longrightarrow \qquad \textbf{A} = 0.20\,\hat{j}\ \text{ m}^2\]
🔍 Why the Dot Product Is a Boss Move

The dot product geometrically extracts only the component of one vector that lines up with another. Here it automatically selects \(B_j\) (the component aligned with the normal \(\hat{j}\)) and throws away \(B_x\) (the component tangent to the surface). No manual decomposition needed.

Step 1 — Set Up the Dot Product

For a flat surface with uniform field:

\[\Phi_B = \textbf{B}\cdot\textbf{A} = (0.433\,\hat{i} + 0.25\,\hat{j})\cdot(0.20\,\hat{j})\]
Step 2 — Distribute and Apply Orthogonality

Expand the dot product term by term:

\[\Phi_B = (0.433\,\hat{i})\cdot(0.20\,\hat{j}) + (0.25\,\hat{j})\cdot(0.20\,\hat{j})\]

Now use the fundamental dot product results for unit vectors:

\[\hat{i}\cdot\hat{j} = 0 \qquad \text{(perpendicular vectors have zero dot product)}\] \[\hat{j}\cdot\hat{j} = 1 \qquad \text{(parallel unit vectors have dot product 1)}\]

So the \(\hat{i}\) term vanishes completely — that's the part of the field running parallel to the surface, which contributes no flux:

\[\Phi_B = 0 + (0.25)(0.20) = 0.050\text{ T·m}^2\]
\[\boxed{\Phi_B = 5.0\times 10^{-2}\text{ Wb}}\]
Sanity check on the field components: Note that \(|\textbf{B}| = \sqrt{0.433^2 + 0.25^2} = \sqrt{0.1875 + 0.0625} = \sqrt{0.25} = 0.50\text{ T}\). The angle: \(\arctan(0.433/0.25) = \arctan(1.732) = 60°\) from the \(\hat{j}\) axis. This is exactly the same setup as Examples 1 and 2 in disguise — confirming all three forms give identical results. ✓
4
Field Parallel to Surface — Zero Flux Geometric Limit

Setup: \(B = 0.80\text{ T}\), \(A = 0.30\text{ m}^2\). The magnetic field runs parallel to the surface — it is not piercing through at all, just sliding along it. Find \(\Phi_B\).

🔍 The Zero-Flux Limit

Imagine holding a window frame and the wind blows exactly sideways — completely parallel to the glass. No air passes through the window. Same principle here: if the field runs tangent to the surface, it threads zero flux through it.

Step 1 — Identify the Angle

If \(\textbf{B}\) is parallel to the surface, it is perpendicular to the surface normal. Therefore the angle between \(\textbf{B}\) and the normal is:

\[\phi = 90°\]
Step 2 — Apply the Formula
\[\Phi_B = BA\cos\phi = BA\cos 90°\]
\[\cos 90° = 0\]
\[\Phi_B = (0.80)(0.30)(0) = 0\]
\[\boxed{\Phi_B = 0}\]
Why this matters: This is not a degenerate edge case. It's physically significant. For a closed surface (like a sphere) surrounding a magnetic dipole, every field line that enters through one side also exits through the other. The net flux is always zero — this is Gauss's Law for Magnetism: \(\oint \textbf{B}\cdot d\textbf{A} = 0\). No magnetic monopoles exist.
5
Non-Uniform Field — The Integral Earns Its Keep Full Integration

Setup: A magnetic field varies with position:

\[B(x) = 2x \quad \text{[T, with } x \text{ in meters]}\]

The field points in the \(+\hat{i}\) direction everywhere. The surface is a rectangle lying in the \(yz\)-plane, with its normal also in the \(+\hat{i}\) direction. It spans \(x = 0\) to \(x = 3.0\text{ m}\) in one direction and has height \(L = 2.0\text{ m}\) in the \(z\)-direction.

🔍 Why You Can't Factor B Out

In the previous examples, \(B\) was uniform — the same everywhere on the surface. Here \(B\) grows as you move in the \(x\)-direction, so strips of the surface near \(x = 3\text{ m}\) have much more flux passing through them than strips near \(x = 0\). You have to sum up the contribution of each infinitesimally thin strip separately.

Step 1 — Geometry: Slice the Surface

Since \(B\) depends only on \(x\), we slice the rectangle into thin vertical strips, each at position \(x\), each with width \(dx\) and height \(L = 2.0\text{ m}\). The area of one strip is:

\[dA = L\,dx = 2.0\,dx\]
Step 2 — Angle Factor

The field \(\textbf{B}\) points in \(+\hat{i}\) and the area normal also points in \(+\hat{i}\). They are parallel, so:

\[\phi = 0° \qquad \Longrightarrow \qquad \cos\phi = 1\]

The field punches straight through — no projection penalty.

Step 3 — Set Up the Integral

The flux contribution from one strip at position \(x\) is:

\[d\Phi_B = B(x)\,dA = (2x)(2.0\,dx) = 4x\,dx\]

Integrate over the full width of the surface:

\[\Phi_B = \int_0^{3.0} 4x\,dx\]
Step 4 — Evaluate the Integral

Factor out the constant \(4\):

\[\Phi_B = 4\int_0^{3.0} x\,dx\]

Use the power rule \(\int x\,dx = \frac{x^2}{2}\):

\[\Phi_B = 4\left[\frac{x^2}{2}\right]_0^{3.0}\]

Evaluate at the limits:

\[\Phi_B = 4\left(\frac{3.0^2}{2} - \frac{0^2}{2}\right) = 4\left(\frac{9.0}{2}\right) = 4(4.5)\]
\[\boxed{\Phi_B = 18\text{ Wb}}\]
Interpretation: The field at \(x = 0\) is \(B = 0\text{ T}\) — no flux there. At \(x = 3\text{ m}\), \(B = 6\text{ T}\) — a strong field. The integral captures this continuously growing contribution across the entire surface. If you had naively used the "average" field \(B_\text{avg} = \frac{0+6}{2} = 3\text{ T}\) times the area \(A = (3)(2) = 6\text{ m}^2\), you'd get \(18\text{ Wb}\) — which actually works here because \(B(x) = 2x\) is linear. For nonlinear fields, only the integral gives the right answer.
6
3D Vector Field, Tilted Plane — Building the Normal Vector Form + Explicit n̂

Setup: A uniform magnetic field is given in full 3D component form:

\[\textbf{B} = 3\,\hat{i} - 2\,\hat{j} + 5\,\hat{k} \quad \text{[T]}\]

The surface is a flat rectangle of area \(A = 0.40\text{ m}^2\) lying in a tilted plane. The plane's outward unit normal is:

\[\hat{n} = \frac{1}{\sqrt{2}}\,\hat{j} + \frac{1}{\sqrt{2}}\,\hat{k}\]

This normal is a unit vector (we'll verify that) pointing at 45° between the \(+y\) and \(+z\) axes. Find \(\Phi_B\).

🔍 What the Normal Encodes

The area vector \(d\textbf{A} = \hat{n}\,dA\) is the normal dressed up with an area magnitude. The direction of \(\hat{n}\) tells the dot product which way "through" the surface means. A normal of \(\frac{1}{\sqrt{2}}\hat{j} + \frac{1}{\sqrt{2}}\hat{k}\) means the surface is tilted — it faces equally toward \(+y\) and \(+z\). So both the \(\hat{j}\) and \(\hat{k}\) components of \(\textbf{B}\) will contribute to flux, but the \(\hat{i}\) component (parallel to this surface) contributes nothing.

Step 1 — Verify \(\hat{n}\) Is a Unit Vector

Always check that \(|\hat{n}| = 1\) before using it. If it's not normalized, the area element is wrong.

\[|\hat{n}| = \sqrt{\left(\frac{1}{\sqrt{2}}\right)^2 + \left(\frac{1}{\sqrt{2}}\right)^2} = \sqrt{\frac{1}{2} + \frac{1}{2}} = \sqrt{1} = 1 \quad \checkmark\]
Step 2 — Write the Area Vector

For a flat surface with uniform field, the entire area contributes uniformly, so:

\[d\textbf{A} = \hat{n}\,dA \qquad \Longrightarrow \qquad \textbf{A} = \hat{n}\,A = 0.40\left(\frac{1}{\sqrt{2}}\,\hat{j} + \frac{1}{\sqrt{2}}\,\hat{k}\right)\text{ m}^2\]
\[\textbf{A} = \frac{0.40}{\sqrt{2}}\,\hat{j} + \frac{0.40}{\sqrt{2}}\,\hat{k}\]
Step 3 — Compute the Dot Product

Now execute \(\Phi_B = \textbf{B}\cdot\textbf{A}\), expanding component by component:

\[\Phi_B = (3\,\hat{i} - 2\,\hat{j} + 5\,\hat{k})\cdot\left(\frac{0.40}{\sqrt{2}}\,\hat{j} + \frac{0.40}{\sqrt{2}}\,\hat{k}\right)\]

Distribute. Only matching unit vectors survive (\(\hat{i}\cdot\hat{j}=0\), \(\hat{i}\cdot\hat{k}=0\), \(\hat{j}\cdot\hat{k}=0\), \(\hat{j}\cdot\hat{j}=1\), \(\hat{k}\cdot\hat{k}=1\)):

\[\Phi_B = (3)\left(\frac{0.40}{\sqrt{2}}\right)\underbrace{\hat{i}\cdot\hat{j}}_{0} + (-2)\left(\frac{0.40}{\sqrt{2}}\right)\underbrace{\hat{j}\cdot\hat{j}}_{1} + (5)\left(\frac{0.40}{\sqrt{2}}\right)\underbrace{\hat{k}\cdot\hat{k}}_{1}\]
\[\Phi_B = 0 + \frac{-2 \times 0.40}{\sqrt{2}} + \frac{5 \times 0.40}{\sqrt{2}}\]
\[\Phi_B = \frac{-0.80 + 2.00}{\sqrt{2}} = \frac{1.20}{\sqrt{2}}\]
\[\Phi_B = \frac{1.20}{1.4142} \approx 0.849\text{ T·m}^2\]
\[\boxed{\Phi_B \approx 0.849\text{ Wb}}\]
What just happened geometrically: The \(\hat{i}\) component of \(\textbf{B}\) (\(= 3\text{ T}\)) ran completely parallel to the tilted surface — zero contribution. The \(-\hat{j}\) component tried to push flux in the \(-y\) direction through a surface that faces \(+y\) — negative contribution. The \(+\hat{k}\) component pushed flux in a direction the surface faces — positive contribution. The net is positive because the \(\hat{k}\) term won. The \(\frac{1}{\sqrt{2}}\) factors appear because the surface is tilted at 45°, so neither \(\hat{j}\) nor \(\hat{k}\) gets full credit.
7
Position-Dependent Vector Field Over a Flat Surface ∫∫ B⃗(x,y) · dA⃗ — Double Integral

Setup: A magnetic field depends on position in the \(xy\)-plane:

\[\textbf{B}(x,y) = xy\,\hat{i} + 3y^2\,\hat{j} + 2\,\hat{k} \quad \text{[T]}\]

The surface is a flat rectangle lying in the \(xy\)-plane, with its outward normal pointing in the \(+\hat{k}\) direction. It spans \(0 \le x \le 2\text{ m}\) and \(0 \le y \le 1\text{ m}\). Find \(\Phi_B\).

🔍 Three Field Components, One Survives

The surface lies in the \(xy\)-plane, so its normal is \(\hat{k}\). The dot product \(\textbf{B}\cdot\hat{k}\) will kill the \(\hat{i}\) and \(\hat{j}\) parts of \(\textbf{B}\) immediately — only the \(\hat{k}\) component of the field threads through the surface. The \(xy\) and \(3y^2\) terms run tangent to the surface and vanish. Only the constant \(2\,\hat{k}\) punches through. That's the whole game before we even integrate.

Step 1 — Identify \(d\textbf{A}\)

The surface lies in the \(xy\)-plane with normal \(\hat{k}\). The area element is:

\[d\textbf{A} = \hat{k}\,dA = \hat{k}\,dx\,dy\]

Here we write \(dA = dx\,dy\) because we're integrating over a rectangle — tiny rectangular patches of width \(dx\) and height \(dy\).

Step 2 — Compute \(\textbf{B}\cdot d\textbf{A}\)
\[\textbf{B}\cdot d\textbf{A} = (xy\,\hat{i} + 3y^2\,\hat{j} + 2\,\hat{k})\cdot(\hat{k}\,dx\,dy)\]

Distribute the dot product through each term:

\[= xy(\hat{i}\cdot\hat{k})\,dx\,dy + 3y^2(\hat{j}\cdot\hat{k})\,dx\,dy + 2(\hat{k}\cdot\hat{k})\,dx\,dy\]
\[= xy(0)\,dx\,dy + 3y^2(0)\,dx\,dy + 2(1)\,dx\,dy\]
\[\textbf{B}\cdot d\textbf{A} = 2\,dx\,dy\]

The first two terms vanish. The integrand collapses to just the constant \(2\).

Step 3 — Set Up and Evaluate the Double Integral
\[\Phi_B = \int_0^1\int_0^2 2\,dx\,dy\]

Inner integral (over \(x\), holding \(y\) fixed):

\[\int_0^2 2\,dx = \left[2x\right]_0^2 = 4\]

Outer integral (over \(y\)):

\[\Phi_B = \int_0^1 4\,dy = \left[4y\right]_0^1 = 4\text{ Wb}\]
\[\boxed{\Phi_B = 4\text{ Wb}}\]
Why did the spatial variation in \(\textbf{B}\) not matter here? The \(xy\) and \(3y^2\) terms varied over the surface, but they were aligned with \(\hat{i}\) and \(\hat{j}\) — running tangent to the surface, contributing zero flux regardless of how they varied. Only the \(\hat{k}\) component mattered, and it was the constant \(2\text{ T}\). So the flux is just \(2\text{ T} \times (2\times 1)\text{ m}^2 = 4\text{ Wb}\). The double integral confirmed this — and when fields aren't constant in the surviving component, you'll really need it.
8
Curved Surface — Normal from the Cross Product Parameterized Surface + dA⃗ = (∂r/∂u × ∂r/∂v) du dv

Setup: A uniform magnetic field points in the \(+z\)-direction:

\[\textbf{B} = B_0\,\hat{k} \quad \text{[T, constant]}\]

The surface is the curved wall of a cylinder of radius \(R\) and height \(h\), oriented with its axis along \(+z\). Find the flux through this curved lateral surface.

This one requires building \(d\textbf{A}\) from scratch using a parameterization.

🔍 The Answer Before the Math

The curved wall of a cylinder has its outward normal pointing radially outward — in the \(\hat{r}\) direction in cylindrical coordinates, which has no \(z\)-component. The field \(\textbf{B} = B_0\hat{k}\) points purely in \(z\). These two vectors are perpendicular everywhere on the curved surface. Dot product of perpendicular vectors is zero. So \(\Phi_B = 0\) — the field slides along the surface without piercing through it. The math should confirm this. Let's grind it out anyway because the machinery is the lesson.

Step 1 — Parameterize the Surface

A cylinder of radius \(R\) can be parameterized using two parameters: \(\phi \in [0, 2\pi)\) (angle around the axis) and \(z \in [0, h]\) (height up the axis). The position vector of any point on the curved wall is:

\[\textbf{r}(\phi, z) = R\cos\phi\,\hat{i} + R\sin\phi\,\hat{j} + z\,\hat{k}\]

This traces out every point on the cylinder as \(\phi\) and \(z\) vary over their ranges.

Step 2 — Compute the Partial Derivatives

To build \(d\textbf{A}\), we need the tangent vectors to the surface in the \(\phi\) and \(z\) directions.

Partial derivative with respect to \(\phi\) (tangent in the "around" direction):

\[\frac{\partial \textbf{r}}{\partial \phi} = -R\sin\phi\,\hat{i} + R\cos\phi\,\hat{j} + 0\,\hat{k}\]

Partial derivative with respect to \(z\) (tangent in the "up" direction):

\[\frac{\partial \textbf{r}}{\partial z} = 0\,\hat{i} + 0\,\hat{j} + 1\,\hat{k}\]

These two tangent vectors lie in the surface. Their cross product gives a vector perpendicular to the surface — the outward normal — and its magnitude gives the area scaling factor.

Step 3 — Compute the Cross Product \(\frac{\partial\textbf{r}}{\partial\phi}\times\frac{\partial\textbf{r}}{\partial z}\)

Evaluate the \(3\times 3\) determinant:

\[\frac{\partial\textbf{r}}{\partial\phi}\times\frac{\partial\textbf{r}}{\partial z} = \begin{vmatrix}\hat{i} & \hat{j} & \hat{k} \\ -R\sin\phi & R\cos\phi & 0 \\ 0 & 0 & 1\end{vmatrix}\]

Expand along the top row:

\[= \hat{i}\begin{vmatrix}R\cos\phi & 0 \\ 0 & 1\end{vmatrix} - \hat{j}\begin{vmatrix}-R\sin\phi & 0 \\ 0 & 1\end{vmatrix} + \hat{k}\begin{vmatrix}-R\sin\phi & R\cos\phi \\ 0 & 0\end{vmatrix}\]
\[= \hat{i}(R\cos\phi\cdot 1 - 0\cdot 0) - \hat{j}((-R\sin\phi)\cdot 1 - 0\cdot 0) + \hat{k}((-R\sin\phi)\cdot 0 - R\cos\phi\cdot 0)\]
\[= R\cos\phi\,\hat{i} + R\sin\phi\,\hat{j} + 0\,\hat{k}\]

This is the outward normal vector field of the cylinder! It points in the radial direction \(\hat{r}\) and has magnitude \(R\):

\[\left|\frac{\partial\textbf{r}}{\partial\phi}\times\frac{\partial\textbf{r}}{\partial z}\right| = \sqrt{(R\cos\phi)^2 + (R\sin\phi)^2} = R\sqrt{\cos^2\phi+\sin^2\phi} = R\]

So the vector area element is:

\[d\textbf{A} = \left(\frac{\partial\textbf{r}}{\partial\phi}\times\frac{\partial\textbf{r}}{\partial z}\right)d\phi\,dz = (R\cos\phi\,\hat{i} + R\sin\phi\,\hat{j})\,d\phi\,dz\]

The magnitude \(R\,d\phi\,dz\) is the area of a tiny rectangular patch on the cylinder (arc length \(R\,d\phi\) times height \(dz\)). Makes sense.

Step 4 — Compute \(\textbf{B}\cdot d\textbf{A}\)
\[\textbf{B}\cdot d\textbf{A} = (B_0\,\hat{k})\cdot(R\cos\phi\,\hat{i} + R\sin\phi\,\hat{j})\,d\phi\,dz\]
\[= B_0\,R\cos\phi\underbrace{(\hat{k}\cdot\hat{i})}_{0}\,d\phi\,dz + B_0\,R\sin\phi\underbrace{(\hat{k}\cdot\hat{j})}_{0}\,d\phi\,dz = 0\]

Every term vanishes. The integrand is identically zero everywhere on the surface.

Step 5 — Integrate (trivially)
\[\Phi_B = \int_0^h\int_0^{2\pi} 0\,d\phi\,dz = 0\]
\[\boxed{\Phi_B = 0}\]
Geometric interpretation: The cross product \(\frac{\partial\textbf{r}}{\partial\phi}\times\frac{\partial\textbf{r}}{\partial z} = R\cos\phi\,\hat{i} + R\sin\phi\,\hat{j}\) is exactly the outward radial vector \(R\hat{r}\). It has no \(\hat{k}\) component — the cylinder's wall never faces up or down, only outward. So \(\textbf{B} = B_0\hat{k}\) is always perpendicular to the outward normal of the curved wall, everywhere. The dot product is zero at every single point. This is why the lateral surface catches zero flux from a vertical field — verified rigorously through the full parameterization machinery.
Why this machinery matters for EE: This cross-product method for building \(d\textbf{A}\) is the general tool for any surface you can parameterize. It's how you handle solenoids, toroids, and any curved geometry where the normal isn't obvious. The same framework appears in Faraday's law — computing the EMF induced in a loop requires the flux integral through an arbitrary surface bounded by that loop.

Why the Three Forms Are Identical

This isn't magic — it's just vector decomposition. Start from the dot product definition:

\[\textbf{B}\cdot d\textbf{A} = |\textbf{B}||d\textbf{A}|\cos\phi = B\,dA\cos\phi\]

And by definition of the perpendicular component:

\[B_\perp \equiv B\cos\phi \qquad \Longrightarrow \qquad B_\perp\,dA = B\cos\phi\,dA\]

Putting it together:

\[\underbrace{\textbf{B}\cdot d\textbf{A}}_{\text{Form 3}} = \underbrace{B\cos\phi\,dA}_{\text{Form 1}} = \underbrace{B_\perp\,dA}_{\text{Form 2}}\]

They are three different ways of thinking about the same geometric operation: extracting the component of \(\textbf{B}\) that is perpendicular to the surface.


When to Use Each Form

Form 1 — B cos φ dA
\[\Phi_B = \int B\cos\phi\,dA\]

You know the angle \(\phi\) between \(\textbf{B}\) and the normal. Clearest for angle problems.

Form 2 — B⊥ dA
\[\Phi_B = \int B_\perp\,dA\]

You want to decompose and discard the parallel component explicitly. Geometrically intuitive.

Form 3 — B⃗ · dA⃗
\[\Phi_B = \int \textbf{B}\cdot d\textbf{A}\]

Field in component form, or surface normal is explicit. Most general flat or tilted surface form.

Form 4 — Parameterized
\[\Phi_B = \iint \textbf{B}\cdot\!\left(\tfrac{\partial\textbf{r}}{\partial u}\times\tfrac{\partial\textbf{r}}{\partial v}\right)du\,dv\]

Curved surfaces. Build \(d\textbf{A}\) from tangent vector cross product. Most general form of all.

⚠ Final Angle Warning — Burned Into Your Brain

\(\phi\) is measured from the surface normal, not from the surface plane. If you're given the complement angle \(\theta\) (angle from the surface), then \(\cos\phi = \sin\theta\) and your formula becomes \(\Phi_B = BA\sin\theta\). Every semester, this distinction takes people out. Don't be that person.