Math Core

Lesson 5.5 · Vector Calculus

Curl and divergence

A vector field can do two basic things near a point: swirl around it or spread out from it. The curl measures the swirl and the divergence measures the spreading. Both are computed from partial derivatives using the "del" operator ∇\nabla, and they are the key ingredients in the two big theorems still to come, Stokes' theorem and the divergence theorem.

The del operator

Write ∇\nabla as a vector of differentiation operators:

∇=⟨∂∂x, ∂∂y, ∂∂z⟩.\nabla = \left\langle \frac{\partial}{\partial x},\ \frac{\partial}{\partial y},\ \frac{\partial}{\partial z}\right\rangle.

Applied to a scalar function, ∇f\nabla f is the gradient. Applied to a vector field F=⟨P,Q,R⟩\mathbf{F} = \langle P, Q, R\rangle with a cross product or a dot product, it gives the curl and the divergence.

Curl

Definition

Curl

The curl of F=⟨P,Q,R⟩\mathbf{F} = \langle P, Q, R\rangle is the vector field

curl⁡F=∇×F=∣ijk∂∂x∂∂y∂∂zPQR∣=⟨Ry−Qz, Pz−Rx, Qx−Py⟩.\operatorname{curl}\mathbf{F} = \nabla\times\mathbf{F} = \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k} \\ \dfrac{\partial}{\partial x} & \dfrac{\partial}{\partial y} & \dfrac{\partial}{\partial z} \\ P & Q & R \end{vmatrix} = \left\langle R_y - Q_z,\ P_z - R_x,\ Q_x - P_y\right\rangle.

Meaning. Imagine F\mathbf{F} is the velocity of a fluid and you drop a tiny paddle wheel into it. The wheel spins fastest when its axle points along curl⁡F\operatorname{curl}\mathbf{F}, and the length of the curl is proportional to how fast it spins (twice the angular speed, in fact). If curl⁡F=0\operatorname{curl}\mathbf{F} = \mathbf{0}, the field is irrotational. Notice that the third component, Qx−PyQ_x - P_y, is exactly the integrand in Green's theorem.

Worked example: Computing curl and divergence

Let F=⟨xy, yz2, x2z⟩\mathbf{F} = \langle xy,\ yz^2,\ x^2z\rangle. Then P=xyP = xy, Q=yz2Q = yz^2, R=x2zR = x^2z, and

curl⁡F=⟨0−2yz, 0−2xz, 0−x⟩=⟨−2yz, −2xz, −x⟩.\operatorname{curl}\mathbf{F} = \left\langle 0 - 2yz,\ 0 - 2xz,\ 0 - x\right\rangle = \langle -2yz,\ -2xz,\ -x\rangle.

The divergence (defined below) is ∂∂x(xy)+∂∂y(yz2)+∂∂z(x2z)=y+z2+x2\dfrac{\partial}{\partial x}(xy) + \dfrac{\partial}{\partial y}(yz^2) + \dfrac{\partial}{\partial z}(x^2z) = y + z^2 + x^2.

Curl and conservative fields

Two facts connect curl to the previous lessons.

Curl detects conservative fields

  1. For any ff with continuous second partial derivatives, curl⁡(∇f)=0\operatorname{curl}(\nabla f) = \mathbf{0}. So every conservative field is irrotational.
  2. Conversely, if F\mathbf{F} has continuous partial derivatives on all of R3\mathbb{R}^3 (or any simply connected region) and curl⁡F=0\operatorname{curl}\mathbf{F} = \mathbf{0}, then F\mathbf{F} is conservative.

Fact 1 is Clairaut's Theorem: for instance, the first component of curl⁡(∇f)\operatorname{curl}(\nabla f) is fzy−fyz=0f_{zy} - f_{yz} = 0. Fact 2 is the three-dimensional version of the Py=QxP_y = Q_x test.

Worked example: Is it conservative?

Let F=⟨y2z3, 2xyz3, 3xy2z2⟩\mathbf{F} = \langle y^2z^3,\ 2xyz^3,\ 3xy^2z^2\rangle. Then

Ry−Qz=6xyz2−6xyz2=0,Pz−Rx=3y2z2−3y2z2=0,Qx−Py=2yz3−2yz3=0.\begin{aligned} R_y - Q_z &= 6xyz^2 - 6xyz^2 = 0, \\ P_z - R_x &= 3y^2z^2 - 3y^2z^2 = 0, \\ Q_x - P_y &= 2yz^3 - 2yz^3 = 0. \end{aligned}

So curl⁡F=0\operatorname{curl}\mathbf{F} = \mathbf{0} on R3\mathbb{R}^3 and F\mathbf{F} is conservative. Integrating as in the previous lesson gives the potential f=xy2z3f = xy^2z^3.

Divergence

Definition

Divergence

The divergence of F=⟨P,Q,R⟩\mathbf{F} = \langle P, Q, R\rangle is the scalar function

div⁡F=∇⋅F=∂P∂x+∂Q∂y+∂R∂z.\operatorname{div}\mathbf{F} = \nabla\cdot\mathbf{F} = \frac{\partial P}{\partial x} + \frac{\partial Q}{\partial y} + \frac{\partial R}{\partial z}.

In the plane, div⁡⟨P,Q⟩=Px+Qy\operatorname{div}\langle P, Q\rangle = P_x + Q_y.

Meaning. If F\mathbf{F} is a fluid velocity, div⁡F(P)\operatorname{div}\mathbf{F}(P) is the net rate at which fluid flows out of a tiny box around PP, per unit volume. Positive divergence means a source (fluid spreading out), negative means a sink, and div⁡F=0\operatorname{div}\mathbf{F} = 0 everywhere means the fluid is incompressible. For example, ⟨x,y,z⟩\langle x, y, z\rangle has divergence 33: it pushes outward everywhere.

A second identity pairs with curl⁡(∇f)=0\operatorname{curl}(\nabla f) = \mathbf{0}:

div⁡(curl⁡F)=0\operatorname{div}(\operatorname{curl}\mathbf{F}) = 0

for any F\mathbf{F} with continuous second partials. (Expand it and the mixed partials cancel in pairs.) Also, div⁡(∇f)=fxx+fyy+fzz\operatorname{div}(\nabla f) = f_{xx} + f_{yy} + f_{zz} is called the Laplacian of ff, written ∇2f\nabla^2 f.

Worked example: A field that is not a curl

Is there a vector field G\mathbf{G} with curl⁡G=⟨x, y, z⟩\operatorname{curl}\mathbf{G} = \langle x,\ y,\ z\rangle?

If there were, then div⁡(curl⁡G)\operatorname{div}(\operatorname{curl}\mathbf{G}) would be 00. But div⁡⟨x,y,z⟩=1+1+1=3≠0\operatorname{div}\langle x, y, z\rangle = 1 + 1 + 1 = 3 \ne 0. So no such G\mathbf{G} exists.

Common mistake

Curl is a vector and divergence is a scalar. The expression curl⁡(div⁡F)\operatorname{curl}(\operatorname{div}\mathbf{F}) is meaningless, because you can't take the curl of a scalar. Also keep the curl's middle component straight: it is Pz−RxP_z - R_x, not Rx−PzR_x - P_z. The cofactor expansion's minus sign flips it.

Green's theorem in vector form

Curl and divergence let you write Green's theorem in two ways for a plane field F=⟨P,Q⟩\mathbf{F} = \langle P, Q\rangle and a positively oriented curve CC around DD.

Two forms of Green's theorem

Circulation form: ∮CF⋅T ds=∬D(curl⁡F)⋅k dA=∬D(Qx−Py) dA\displaystyle\oint_C \mathbf{F}\cdot\mathbf{T}\,ds = \iint_D (\operatorname{curl}\mathbf{F})\cdot\mathbf{k}\,dA = \iint_D (Q_x - P_y)\,dA.

Flux form: ∮CF⋅n ds=∬Ddiv⁡F dA=∬D(Px+Qy) dA\displaystyle\oint_C \mathbf{F}\cdot\mathbf{n}\,ds = \iint_D \operatorname{div}\mathbf{F}\,dA = \iint_D (P_x + Q_y)\,dA,

where n\mathbf{n} is the outward unit normal to CC.

The flux form says the net outflow across the boundary equals the total of all the sources inside. Stokes' theorem and the divergence theorem are these two statements lifted into three dimensions.

Worked example: Outward flux across a circle

Find the outward flux of F=⟨x,y⟩\mathbf{F} = \langle x, y\rangle across the unit circle.

div⁡F=1+1=2\operatorname{div}\mathbf{F} = 1 + 1 = 2, so the flux is ∬D2 dA=2π\iint_D 2\,dA = 2\pi. Check directly: on the unit circle the outward normal is n=⟨x,y⟩\mathbf{n} = \langle x, y\rangle, so F⋅n=x2+y2=1\mathbf{F}\cdot\mathbf{n} = x^2 + y^2 = 1, and the flux is the length 2π2\pi.

Practice

Practice 1

Find div⁡F\operatorname{div}\mathbf{F} for F=⟨x2y, yz, xz2⟩\mathbf{F} = \langle x^2y,\ yz,\ xz^2\rangle.

Enter an expression, e.g. 3x^2 - 2x + 1

Practice 2

Find curl⁡F\operatorname{curl}\mathbf{F} for F=⟨y, z, x⟩\mathbf{F} = \langle y,\ z,\ x\rangle.

Enter a point like (2, -3)

Practice 3

Let F=⟨x2z, yz2, xy⟩\mathbf{F} = \langle x^2z,\ yz^2,\ xy\rangle. Find curl⁡F\operatorname{curl}\mathbf{F} at the point (1,2,1)(1, 2, 1).

Enter a point like (2, -3)

Practice 4

Find div⁡F\operatorname{div}\mathbf{F} for F=⟨exsin⁡y, excos⁡y, z⟩\mathbf{F} = \langle e^x\sin y,\ e^x\cos y,\ z\rangle.

Enter a number. Fractions like 3/4 and sqrt(2) are OK.

Practice 5

Assume ff and F\mathbf{F} have continuous second partial derivatives. Which expression is always equal to 00?

Practice 6

Find the constant aa that makes F=⟨axy+z, x2, x⟩\mathbf{F} = \langle axy + z,\ x^2,\ x\rangle conservative on R3\mathbb{R}^3.

Enter a number. Fractions like 3/4 and sqrt(2) are OK.

Practice 7

F=⟨2xyz, x2z, x2y+1⟩\mathbf{F} = \langle 2xyz,\ x^2z,\ x^2y + 1\rangle has zero curl. Evaluate ∫CF⋅dr\displaystyle\int_C \mathbf{F}\cdot d\mathbf{r} along any path from (0,0,0)(0, 0, 0) to (1,2,3)(1, 2, 3).

Enter a number. Fractions like 3/4 and sqrt(2) are OK.

Practice 8

Use the flux form of Green's theorem to find the outward flux of F=⟨x3, y3⟩\mathbf{F} = \langle x^3,\ y^3\rangle across the unit circle.

Enter a number. Fractions like 3/4 and sqrt(2) are OK.