Green's theorem

Say we have a piecewise smooth, closed curve CC in R2\mathbb R^2 bounding a region RR in R2\mathbb R^2, oriented counter-clockwise; and a smooth vector field F=(P,Q)\mathbf F=(P, Q) also in R2\mathbb R^2.

Green’s theorem states that integrating a function along CC is equal to integrating the curl of that function over the area RR.

CFdr=Rcurl(F)dA \oint_C \mathbf F \cdot d\mathbf r = \iint_R \mathrm{curl}(\mathbf F) dA

Or

CPdx+Qdy=R(QxPy)dA \oint_C Pdx + Qdy = \iint_R (Q_x - P_y) dA

Normal form

Say we have the same vector field F\mathbf F and curve CC enclosing region RR. We can also use Green’s theorem to relate the integral across CC of Fn^\mathbf F \cdot \hat{\mathbf{n} } (the outwards-pointing vector normal to the curve) to an integral over RR.

fluxC(F)=C(Fn^)ds=Rdiv(F)dA \mathrm{flux}_C(\mathbf F) = \oint_C (\mathbf F \cdot \hat{\mathbf{n} }) ds = \iint_R \mathrm{div}(\mathbf F) dA
Proof

The proof for the normal form is given in terms of Green’s theorem.

We can write the integral as CFn^ds\oint_C \mathbf F \cdot \hat{\mathbf n} ds. Note that n^ds\hat{\mathbf n} ds has the same magnitude and is perpendicular to the differential change in position drd\mathbf r. If dr=(dx,dy)d\mathbf r = (dx,dy), then n^=(dy,dx)\hat{\mathbf n} = (dy, -dx).

The dot product is then Fn^ds=PdyQdx\mathbf F \cdot \hat{\mathbf n}ds = Pdy - Qdx.

Now let G=(Q,P)\mathbf G = (-Q, P). This means Fn^ds=Gdr\mathbf F \cdot \hat{\mathbf n}ds = \mathbf G \cdot d\mathbf r.

By Green’s theorem, CGdr=Rcurl(G)dA\oint_C \mathbf G \cdot d\mathbf r = \iint_R \mathrm{curl}(\mathbf G) dA.

curl(G)=curl(QP)=Px+Qy=div(F)\mathrm{curl}(\mathbf G) = \mathrm{curl} \begin{pmatrix}-Q \\ P\end{pmatrix} = P_x + Q_y = \mathrm{div}(\mathbf F).

Area

Green’s theorem can be used to find the area of a 2D shape. The integral we would normally use to calculate the area is just R1dA\iint_R 1\,dA. In order to use Green’s theorem, we need to find a function F=(P,Q)\mathbf F = (P, Q) such that QxPy=curl(F)=1Q_x - P_y = \mathrm{curl}(\mathbf F) = 1. By inspection we can find F=(0,x)\mathbf F = (0, x) or F=(y,0)\mathbf F = (-y, 0).

Then we can integrate CFdr\oint_C \mathbf F \cdot d\mathbf r over the boundary of the region to find its area.

Note

Remember, Green’s theorem only works if the curve is oriented counter-clockwise. If you calculate it clockwise, you’ll need to flip the sign.