Unit 5, lesson 38

Green’s Theorem

How the circulation and the outward flux of a plane field around a closed boundary equal double integrals of curl and divergence, which hypotheses this needs, and how holes and areas fit in.

  • Theory part 75
  • Practice part 76
  • Divergence, Curl, and Green’s Theorem

Walk once around the edge of a pond and add up how much the current pushes you along. That total looks like it depends only on the edge. Green’s theorem says it equals the sum of all the tiny swirls inside the pond. Small loops inside cancel along shared edges, and only the outer edge survives.

Key terms

Positively oriented piecewise smooth boundary

Definition A curve CC given by a continuous r(t)\mathbf{r}(t) for a≤t≤ba \le t \le b is a simple closed curve if r(a)=r(b)\mathbf{r}(a) = \mathbf{r}(b) and r\mathbf{r} is one-to-one on [a,b)[a, b). It is piecewise smooth if there are finitely many points a=t0<t1<⋯<tk=ba = t_0 < t_1 < \dots < t_k = b such that on each closed piece [ti−1,ti][t_{i-1}, t_i] the map r\mathbf{r} is C1C^1, with one-sided derivatives at the two ends, and

r′(t)≠0for ti−1<t<ti.\mathbf{r}'(t) \neq \mathbf{0} \quad \text{for } t_{i-1} < t < t_i.

So the velocity may vanish or jump only at the finitely many joining points tit_i. A line integral is an ordinary integral in tt, and finitely many values of tt do not affect it, so these points do no harm. This is slightly more permissive than the piecewise regular curves of the lesson on scalar line integrals, which ask for a nonzero velocity at the ends of each piece too, and every piecewise regular curve is piecewise smooth in this sense. A simple closed curve splits the plane into a bounded inside and an unbounded outside, and DD denotes the inside together with CC. The curve is positively oriented if DD lies to the left as you travel along CC. For a single closed curve this means counterclockwise.

In plain words The boundary may have finitely many corners or cusps, where the smooth pieces join. Walk along it with the region always on your left hand.

Example The boundary of the unit square [0,1]×[0,1][0, 1] \times [0, 1] is positively oriented when traversed

(0,0)→(1,0)→(1,1)→(0,1)→(0,0).\begin{gathered} (0, 0) \to (1, 0) \to (1, 1) \\ \to (0, 1) \to (0, 0). \end{gathered}

It is four straight segments joined at four corners, so it is piecewise smooth but not smooth.

A C1 field on a neighborhood of the region

Definition The theorem asks for PP and QQ to have continuous first partial derivatives on an open set UU that contains the whole closed region DD together with its boundary CC. Being defined and smooth only on the inside of DD is not enough, and a single point of DD where PP or QQ is undefined breaks the hypothesis.

In plain words The field must be smooth everywhere on the region and a little beyond its edge, with no bad points hiding inside.

Example The field ⟨x2y,ex⟩\langle x^2 y, e^{x} \rangle is C1C^1 on all of R2\mathbb{R}^2, so any region is fine. The field

F=⟨−yx2+y2, xx2+y2⟩\mathbf{F} = \left\langle \frac{-y}{x^2 + y^2},\ \frac{x}{x^2 + y^2} \right\rangle

is C1C^1 only on the plane minus the origin. It meets the hypothesis on the annulus 1≤x2+y2≤41 \le x^2 + y^2 \le 4, but not on the unit disk.

Circulation form

Definition Let CC be a positively oriented, piecewise smooth, simple closed curve in the plane, let DD be the closed region made of CC and its inside, and let PP and QQ be C1C^1 on an open set containing DD. Then

∮CP dx+Q dy=∬D(∂Q∂x−∂P∂y)dA.\begin{aligned} &\oint_C P \, dx + Q \, dy \\[0.4em] &\quad = \iint_D \left( \frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} \right) dA. \end{aligned}

The left side is the circulation of F=⟨P,Q⟩\mathbf{F} = \langle P, Q \rangle around CC, and the integrand on the right is the scalar curl.

In plain words Total circulation around the edge equals the total curl inside. Reverse the direction of CC and the left side changes sign, so the orientation is part of the theorem.

Example Take F=⟨xy,x2⟩\mathbf{F} = \langle xy, x^2 \rangle on the unit square. Here Qx−Py=2x−x=xQ_x - P_y = 2x - x = x, so

∬Dx dA=∫01 ⁣ ⁣∫01x dx dy=12.\iint_D x \, dA = \int_0^1 \!\! \int_0^1 x \, dx \, dy = \frac{1}{2}.

Directly, the bottom edge has y=0y = 0 and gives 00. The right edge has x=1x = 1, dx=0dx = 0 and gives ∫011 dy=1\int_0^1 1 \, dy = 1. The top edge has y=1y = 1 and runs from x=1x = 1 to x=0x = 0, giving

∫10x dx=−12.\int_1^0 x \, dx = -\frac{1}{2}.

The left edge has x=0x = 0 and gives 00. The total is 12\frac{1}{2}, as promised.

Flux form

Definition Under the same hypotheses, let n\mathbf{n} be the outward unit normal to CC, which for positive orientation is the unit tangent turned clockwise by a right angle. Then

∮CF⋅n ds=∮CP dy−Q dx=∬D(Px+Qy)dA.\begin{aligned} \oint_C \mathbf{F} \cdot \mathbf{n} \, ds &= \oint_C P \, dy - Q \, dx \\[0.4em] &= \iint_D \left( P_x + Q_y \right) dA. \end{aligned}

The left side is the outward flux of F\mathbf{F} across CC, and the integrand on the right is the divergence.

In plain words Net flow out through the edge equals the total source strength inside. Circulation and flux are different measurements of the same field, so keep them apart.

Example Let F=⟨x,y⟩\mathbf{F} = \langle x, y \rangle and let DD be the disk of radius 22. The divergence is 22, so the right side is 2⋅4π=8π2 \cdot 4\pi = 8\pi. Directly, on the circle of radius 22 the outward normal is ⟨x,y⟩/2\langle x, y \rangle / 2, so F⋅n=2\mathbf{F} \cdot \mathbf{n} = 2, and the circle has length 4π4\pi. The flux is 2⋅4π=8π2 \cdot 4\pi = 8\pi.

Simple regions and regions with holes

Definition The proof first handles simple regions, which are both of type I, lying between two graphs y=g1(x)y = g_1(x) and y=g2(x)y = g_2(x), and of type II, lying between x=h1(y)x = h_1(y) and x=h2(y)x = h_2(y). It extends to any region that can be cut into finitely many simple regions, because the line integrals along the cuts are traversed once in each direction and cancel.

A region with holes is handled the same way. Let C0C_0 be a piecewise smooth simple closed curve, and let C1,…,CmC_1, \dots, C_m be piecewise smooth simple closed curves inside C0C_0 that do not meet C0C_0 or each other, none of them inside another. Let DD be the closed region inside C0C_0 and outside every CjC_j. Its boundary is the union of all these curves, and the positive orientation keeps DD on the left, which means C0C_0 counterclockwise and each inner curve CjC_j clockwise. If PP and QQ are C1C^1 on an open set containing DD, then both forms of the theorem hold with the line integral taken as the sum over all boundary curves. With every curve written counterclockwise, the circulation form reads

∮C0F⋅dr−∑j=1m∮CjF⋅dr=∬D(Qx−Py)dA.\begin{aligned} &\oint_{C_0} \mathbf{F} \cdot d\mathbf{r} - \sum_{j=1}^{m} \oint_{C_j} \mathbf{F} \cdot d\mathbf{r} \\[0.4em] &\quad = \iint_D \left( Q_x - P_y \right) dA. \end{aligned}

The field may be undefined inside the holes, since the holes are not part of DD.

In plain words Holes are allowed as long as the field is smooth on what remains, including the inner edges. Traverse the inner curves clockwise and add them in.

Example Let F=⟨−y,x⟩/(x2+y2)\mathbf{F} = \langle -y, x \rangle / (x^2 + y^2) and let DD be the annulus 1≤x2+y2≤41 \le x^2 + y^2 \le 4. On the plane minus the origin Qx−Py=0Q_x - P_y = 0. Let C1C_1 and C2C_2 be the circles of radius 11 and 22 about the origin, both counterclockwise. Green’s theorem gives

∮C2F⋅dr−∮C1F⋅dr=0.\oint_{C_2} \mathbf{F} \cdot d\mathbf{r} - \oint_{C_1} \mathbf{F} \cdot d\mathbf{r} = 0.

The minus sign comes from traversing the inner circle clockwise. Each counterclockwise circle of radius RR gives

∫02π(sin⁡2t+cos⁡2t)dt=2π,\int_0^{2\pi} \left( \sin^2 t + \cos^2 t \right) dt = 2\pi,

and indeed 2π−2π=02\pi - 2\pi = 0. The same argument, with a small circle about the origin as the hole, shows that every positively oriented, piecewise smooth, simple closed curve with the origin inside gives circulation 2π2\pi.

Area formula

Definition Let DD satisfy the hypotheses of Green’s theorem with positively oriented boundary CC. Choosing fields with Qx−Py=1Q_x - P_y = 1 gives

area⁡(D)=∮Cx dy=−∮Cy dx=12∮Cx dy−y dx.\begin{aligned} \operatorname{area}(D) &= \oint_C x \, dy = -\oint_C y \, dx \\[0.4em] &= \frac{1}{2} \oint_C x \, dy - y \, dx. \end{aligned}

In plain words You can find an area by walking around its edge. This is how a planimeter works.

Example The ellipse x=3cos⁡tx = 3\cos t, y=2sin⁡ty = 2\sin t for 0≤t≤2π0 \le t \le 2\pi is traversed counterclockwise. The third formula gives

12∫02π(6cos⁡2t+6sin⁡2t)dt=6π,\frac{1}{2} \int_0^{2\pi} \left( 6\cos^2 t + 6\sin^2 t \right) dt = 6\pi,

which matches the formula πab\pi a b with a=3a = 3 and b=2b = 2.

Check yourself

  1. Use Green’s theorem to find the circulation of F=⟨y2,3xy⟩\mathbf{F} = \langle y^2, 3xy \rangle around the triangle with vertices (0,0)(0, 0), (1,0)(1, 0), (0,1)(0, 1), oriented counterclockwise.
Show the answer

The field is C1C^1 everywhere and the boundary is positively oriented and piecewise smooth. The scalar curl is 3y−2y=y3y - 2y = y, so

∫01 ⁣ ⁣∫01−xy dy dx=16.\int_0^1 \!\! \int_0^{1 - x} y \, dy \, dx = \frac{1}{6}.
  1. Find the area enclosed by the astroid x=cos⁡3tx = \cos^3 t, y=sin⁡3ty = \sin^3 t, 0≤t≤2π0 \le t \le 2\pi.
Show the answer

The astroid is a simple closed curve traversed counterclockwise. It is piecewise smooth with joining points t=0,π2,π,3π2,2πt = 0, \frac{\pi}{2}, \pi, \frac{3\pi}{2}, 2\pi. On each quarter the velocity ⟨−3cos⁡2tsin⁡t,3sin⁡2tcos⁡t⟩\langle -3\cos^2 t \sin t, 3\sin^2 t \cos t \rangle is continuous and nonzero inside, and it vanishes at the joining points, which are the four cusps. This is exactly what the definition allows. Here xy′−yx′=3sin⁡2tcos⁡2tx y' - y x' = 3\sin^2 t \cos^2 t, so the area is

12∫02π3sin⁡2tcos⁡2t dt=3π8.\frac{1}{2} \int_0^{2\pi} 3 \sin^2 t \cos^2 t \, dt = \frac{3\pi}{8}.
  1. For F=⟨−y,x⟩/(x2+y2)\mathbf{F} = \langle -y, x \rangle / (x^2 + y^2) the scalar curl is 00, yet the circulation around the unit circle is 2π2\pi. Which hypothesis of Green’s theorem fails on the unit disk?
Show the answer

The field must be C1C^1 on an open set containing the whole disk. It is undefined at the origin, which lies inside, so the theorem does not apply. On regions that exclude the origin, such as an annulus, it applies and gives consistent results.

How this lesson connects

It also uses Vector line integrals, divergence and planar curl and Double integrals, boundary orientation and regular planar regions.

Resources for this lesson