Calculus &
Linear Algebra II

Chapter 38

38 Flux of a vector field

By the end of this section, you should be able to answer the following questions:

  • What is the flux of a constant vector field across a flat surface in 3D?
  • What is the flux of a vector field across a plane curve in 2D?

In this section we introduce the concept of flux: In three dimensions, the flux of a vector field across a given surface is defined to be the "flow rate" of the vector field through the surface.



38 Flux of a vector field

Since many vector fields involve no motion (e.g., electric fields, magnetic fields), this definition can be very difficult to comprehend at first.


38 Flux of a vector field

A nice context for working with flux in order to understand its definition is by considering the velocity vector of a fluid (so now we do have motion).

Fluid simulation by Amanda Ghassaei     🔄 (Reset/Re-start)

38 Flux of a vector field

A nice context for working with flux in order to understand its definition is by considering the velocity vector of a fluid (so now we do have motion).

38 Flux of a vector field

In three dimensions, the flux of a fluid across a surface is given in units of volume per unit time. In other words, the flux tells us how much of the fluid (volume) passes through a given surface in one second.


38 Flux of a vector field

Consider a river flowing at a constant velocity of $2m/s$ in only one direction. Now imagine placing a $3m$ square fishing net into the river so that it somehow stays perpendicular to the flow of the river. What is the flux of the water through the net?


38 Flux of a vector field

What is the flux of the water through the net?

$\text{Vol. through net surface after 1 second}$ $= 3 m\times 3 m \times 2 m \times 1 s$ $=18 m^3/s$

Then $\underbrace{\; {\large\text{Flux}}\;}_{{\tiny\text{Amount of water through net per unit of time}}}$ $= 18 ~m^3/s$



38 Flux of a vector field

Now if we rotate the net through an angle $\theta,$ what is the flux through the net?


38 Flux of a vector field

Now if we rotate the net through an angle $\theta,$ what is the flux through the net?


38 Flux of a vector field

Consider the component of the flow perpendicular to the net $\v\pd \n = \norm{\v}\norm{\n}\cos\theta,$ where $\n$ is the unit vector perpendicular to the surface.

$\Rightarrow \text{Flux}$ $=(\v \pd \n) \times (\text{area of net})$
$=\norm{\v}\norm{\n}\cos\theta \times (3\times 3)$
$=2\cos\theta \times (3\times 3) $ $= 18 \cos\theta ~m^3/s$.

Another method: In 1 second, the same volume passes through this net as a net perpendicular to the motion with dimensions $3\times 3\cos \theta .$

\[ \Rightarrow \text{Flux} = \underbrace{ 3\times 3 \cos \theta}_{\text{Perpendicular net}} \times \underbrace{2}_{\text{velocity}} = 18 \cos\theta ~m^3/s. \]


38.1 Flux in 2D

Before we look at the flux of a vector field through more general surfaces, let's look at flux in two dimensions, by considering the flow of a two dimensional fluid through a curve in the $xy$-plane. Note that in this context of a fluid in $2$D, flux has dimensions area per unit time.

To start, consider the problem of calculating the flux of a fluid with constant velocity $\v = 2~\i$ through a line segment $C$ perpendicular to the flow, where $C$ is given by \[ C = \left\{(x,y)~|~ x = 2, 2\leq y \leq 6\right\}. \]



38.1 Flux in 2D

Constant velocity $\v = 2~\i$ through a line segment $C$ perpendicular to the flow, where $C$ is given by $C = \left\{(x,y)~|~ x = 2, 2\leq y \leq 6\right\}.$


38.1 Flux in 2D

Consider the component of the flow perpendicular to the net $\v\pd \n = \norm{\v}\norm{\n}\cos\theta$, where $\n=\i$ is the unit vector perpendicular to the curve $C$.

$\Rightarrow \text{Flux}$ $=(\v \pd \n) \times (\text{length of }C)$
$=\norm{\v}\norm{\n}\cos\theta \times 4$
$=2\cos(0) \times 4 $ $= 8 ~m^2/s$.

If $\n = -\i$, we have that $\theta = \pi$.

\[ \Rightarrow \text{Flux} = (\v \pd \n) \times (\text{length of }C) = -8 ~m^2/s. \]



38.1 Flux in 2D

Now consider calculating the flux of the velocity vector $\v(x,y)$ in the $x$-$y$ plane through a curve $C$. We first divide $C$ up into arcs of length $\Delta S$, and approximate $\v$ as constant over each arc.



38.1 Flux in 2D

This constant vector over each arc shall be evaluated at a representative point in each arc, say $P^* = (x^*, y^*)$.

We also approximate the arc as a straight line, so that

$\displaystyle \Delta S \approx \sqrt{(\Delta x)^2+ (\Delta y)^2} \approx \norm{\r'(t)}\Delta t$     (see Def. in 35.1)


38.1 Flux in 2D

The component of $\v$ which is perpendicular to $C$ (over $\Delta S$) is $\approx \v(P^*)\pd \n (P^*)$.

flux through one arc $\approx \v(P^*)\pd \n (P^*)~\Delta S$.

$\Rightarrow $ total flux through $C$ $\approx \displaystyle \sum_{i}\v(P_i^*)\pd \n (P_i^*)~\Delta S_i$.

If we take the limit as $\Delta S$, we obtain an exact expression for the flux over the entire curve C as a line integral: \[ \text{Flux}= \int_C \v\pd \n ~dS \] where $\n$ is a unit vector normal to $C$.



38.1 Flux in 2D

If we take the limit as $\Delta S$, we obtain an exact expression for the flux over the entire curve C as a line integral: \[ \text{Flux}= \int_C \v\pd \n ~dS \] where $\n$ is a unit vector normal to $C$.

We use this expression as a definition of flux of any two dimensional vector field $\v$ across a plane curve $C$. Note then that \[ \text{dimensions of flux (in 2D)} = (\text{dimensions of }\v) \times (\text{distance}). \]



38.1.1 Evaluating flux in 2D

To evaluate the $2$D flux integral \[ \int_C \v\pd \n ~dS, \] we need to parameterise $C$ and express all quantities in the line integral in terms of this parameterisation (similar to the computation of line integral in the definition of work done).




38.1.1 Evaluating flux in 2D
  1. Parameterise $C$ by finding a $\r(t) = x(t)~\i + y(t)~\j$ with $t \in[a, b]$ that describes $C$.
  2. Write $\v(x,y)$ restricted to $C$ as $\v(\r(t))= \v(x(t),y(t))$.
  3. Compute a unit tangent vector to $C$ by using $\T(x,y) = \dfrac{\r'(t)}{\norm{\r'(t)}}$, where $\r'(t) = x'(t)~\i + y'(t)~\j$ is a tangent vector to $C$.




38.1.1 Evaluating flux in 2D
  1. Be careful of the direction of $\n$. By the definition of vector cross product, and since $\k$ is a unit vector normal to the $x$-$y$ plane, we can take $\n = \T \times \k$.
  2. Write $dS = \norm{\r'(t)}dt$.
  3. Evaluate the $2$D flux integral as a definite integral in terms of the parameter $t$:
    $\displaystyle \int_C \v \pd \n ~dS= \int_a^b \left[ \v\left(\r(t)\right)\pd \dfrac{\r'(t)\times \k}{\norm{\r'(t)}}\right]\norm{\r'(t)}~dt$
    $\displaystyle \quad\;\,= \int_a^b \v\left(\r(t)\right)\pd \left( \r'(t) \times \k\right) dt$.

38.1.1 Evaluating flux in 2D
  1. Parameterise $C$ by finding a $\r(t) = x(t)~\i + y(t)~\j$ with $t \in[a, b]$ that describes $C$.
  2. Write $\v(x,y)$ restricted to $C$ as $\v(\r(t))= \v(x(t),y(t))$.
  3. Compute a unit tangent vector to $C$ by using $\T(x,y) = \dfrac{\r'(t)}{\norm{\r'(t)}}$, where $\r'(t) = x'(t)~\i + y'(t)~\j$ is a tangent vector to $C$.
  4. Be careful of the direction of $\n$. By the definition of vector cross product, and since $\k$ is a unit vector normal to the $x$-$y$ plane, we can take $\n = \T \times \k$.
  5. Write $dS = \norm{\r'(t)}dt$.
  6. Evaluate the $2$D flux integral as a definite integral in terms of the parameter $t$:
    $\displaystyle \int_C \v \pd \n ~dS= \int_a^b \left[ \v\left(\r(t)\right)\pd \dfrac{\r'(t)\times \k}{\norm{\r'(t)}}\right]\norm{\r'(t)}~dt$
    $\displaystyle \quad\;\,= \int_a^b \v\left(\r(t)\right)\pd \left( \r'(t) \times \k\right) dt$.

38.1.1 Evaluating flux in 2D

Note: We may take $\n = \k\times \T$ in step 4 above for an open curve $C$. Namely, one is free to choose the direction of $\n$ (if it is not specified in a question) as long as one is consistent along the entire $C$.


38.1.2 Calculate the flux of $\v = -y~\i+x~\j$ (in the positive $x$ direction) across the line $x = 2$ (for $ 2 \leq y \leq 6$).

38.1.2 Calculate the flux of $\v = -y~\i+x~\j$ (in the positive $x$ direction) across the line $x = 2$ (for $ 2 \leq y \leq 6$).

Flux $\displaystyle =\int_C\v \pd \n ~dS=\;?\;\;$ Parametrise the curve $C$: \[ \r(t) = 2~\i + t~\j,\;\;2\leq t\leq 6. \]

So $\r'(t) = \j,$ $dS = \norm{\r'(t)}~dt=dt$ and

$\v\pd \n $ $=\left(-y~\i + x~\j\right)\pd \i$ $=\left(-t~\i+2~\j\right)\pd \i$ $= -t$.
$\Rightarrow \text{Flux} $ $=\displaystyle \int_C \v\pd \n~dS$ $=\displaystyle \int_2^6 \v\pd \n ~\norm{\r'(t)}~dt$
$=\displaystyle \int_2^6 (-t)~dt$ $=-16$.

38.2 Outward flux across a closed curve
in the plane

Let $C$ be a piecewise-smooth, simple closed curve. Let $v_1(x, y)$, $v_2(x, y)$ be continuous in the region bounded by $C$. (Note that these are some of the conditions of Green's theorem!)

The net outward flux of $\v = v_1~\i + v_2~\j$ across $C$ is given by \[ \text{Net outward flux} = \oint_C \v \pd \n ~dS, \] where $\n$ is a unit vector normal to $C$, directed outward from the region bounded by $C$.



38.2.1 Calculate the flux of $\,\v = xy~\i+xy~\j\,$ across the curve

from $(2,0)$ to $(-2,0)$ via the semicircle of radius $2$ centred at the origin (for $y \geq 0$) followed by the straight line from $(-2,0)$ to $(2,0).$


38.2.1 Calculate the flux of $\,\v = xy~\i+xy~\j\,$ across the curve

from $(2,0)$ to $(-2,0)$ via the semicircle of radius $2$ centred at the origin (for $y \geq 0$) followed by the straight line from $(-2,0)$ to $(2,0)$.

Net Outward Flux $\displaystyle =\oint_C\v \pd \n ~dS=\;?\;\;$ Where $ C = C_1 + C_2. $ So

$\displaystyle \oint_C\v \pd \n ~dS$ $\displaystyle =\int_{C_1}\v \pd \n_1 ~dS+\int_{C_2}\v \pd \n_2 ~dS$.






38.2.1 Calculate the flux of $\,\v = xy~\i+xy~\j\,$ across the curve

from $(2,0)$ to $(-2,0)$ via the semicircle of radius $2$ centred at the origin (for $y \geq 0$) followed by the straight line from $(-2,0)$ to $(2,0)$.

Parametrisation of $C_1$: \[ \r (t) = 2 \cos t~\i+ 2 \sin t~\j, \;\; 0\leq t \leq \pi. \]

So $\r'(t) = -2 \sin t~\i+ 2 \cos t~\j,$ $\norm{\r'(t)} = 2$ and $$dS = \norm{\r'(t)} ~dt = 2~dt.$$

$\Rightarrow \;\displaystyle \n_1$ $\displaystyle = \,\T \times \k$ $\displaystyle = \frac{\r'(t)}{\norm{\r'(t)}} \times \k$ $\displaystyle = \cos t ~\i + \sin t~\j$

$$\Rightarrow \;\v(\r(t)) = 4 \cos t \sin t~\i + 4 \cos t \sin t~\j.$$

38.2.1 Calculate the flux of $\,\v = xy~\i+xy~\j\,$ across the curve

from $(2,0)$ to $(-2,0)$ via the semicircle of radius $2$ centred at the origin (for $y \geq 0$) followed by the straight line from $(-2,0)$ to $(2,0)$.

$$\Rightarrow \v \pd \n_1 = 4 \cos^2 t \sin t + 4 \cos t \sin^2 t, \quad dS = 2 ~dt$$

$ \displaystyle \int_{C_1} \v\pd \n_1~dS $ $\displaystyle = \int_0^{\pi} \big( 4 \cos^2 t \sin t + 4 \cos t \sin^2 t \big)2~dt$
$\displaystyle = \bigg. \left(2 \sin(t) - \frac{2}{3} \sin(3 t)\right.\bigg.$
$ \displaystyle \qquad \bigg. \left. - 2 \cos(t) - \frac{2}{3} \cos(3 t)\right)\bigg|_0^{\pi}$
$\displaystyle =\frac{16}{3}$.

38.2.1 Calculate the flux of $\,\v = xy~\i+xy~\j\,$ across the curve

from $(2,0)$ to $(-2,0)$ via the semicircle of radius $2$ centred at the origin (for $y \geq 0$) followed by the straight line from $(-2,0)$ to $(2,0)$.

Parametrisation of $C_2$: \[ \r (t) = t~\i+ 0~\j, \;\; -2\leq t \leq 2. \]

So $\r'(t) = \i$, $\norm{\r'(t)} = 1$ and $dS = \norm{\r'(t)} ~dt = dt$.

$\Rightarrow \displaystyle \n_2 $ $\displaystyle = \T \times \k$ $\displaystyle = \frac{\r'(t)}{\norm{\r'(t)}} \times \k$ $\displaystyle = \i \times \k$ $\displaystyle = - \j$

$$\Rightarrow \v(\r(t)) = 0~\i + 0~\j.$$



38.2.1 Calculate the flux of $\,\v = xy~\i+xy~\j\,$ across the curve

from $(2,0)$ to $(-2,0)$ via the semicircle of radius $2$ centred at the origin (for $y \geq 0$) followed by the straight line from $(-2,0)$ to $(2,0)$.

$$\Rightarrow \;\v \pd \n_2 = \left( 0~\i + 0~\j\right)\pd (-\j)=0, \quad dS = dt$$

$\Rightarrow \; \displaystyle \int_{C_2} \v\pd \n_2~dS $ $ = 0$.

$$\Rightarrow \text{Net Outward Flux} = \int_{C_1} \v\pd \n_2+ \int_{C_2} \v\pd \n_2 = \frac{16}{3}.$$





Line integrals of vector fields: Flux  


Credits