Circulation Form of Green’s Theorem

Learning Objectives

  • Apply the circulation form of Green’s theorem.

Extending the Fundamental Theorem of Calculus

Recall that the Fundamental Theorem of Calculus says that

[latex]\large{\displaystyle\int_a^b{F}^\prime(x)dx=F(b)-F(a)}[/latex].

As a geometric statement, this equation says that the integral over the region below the graph of [latex]F^\prime(x)[/latex] and above the line segment [latex][a,b][/latex] depends only on the value of [latex]F[/latex] at the endpoints a and b of that segment. Since the numbers [latex]a[/latex] and [latex]b[/latex] are the boundary of the line segment [latex][a,b][/latex], the theorem says we can calculate integral [latex]\displaystyle\int_a^b{F}^\prime(x)dx[/latex] based on information about the boundary of line segment [latex][a,b][/latex] (Figure 1). The same idea is true of the Fundamental Theorem for Line Integrals:

[latex]\large{\displaystyle\int_C\nabla{f}\cdot{d}{\bf{r}}=f({\bf{r}}(b))-f({\bf{r}}(a))}[/latex].

When we have a potential function (an “antiderivative”), we can calculate the line integral based solely on information about the boundary of curve [latex]C[/latex].

<img src="/apps/archive/20220422.171947/resources/1fd7145a16a0b05f17857a2d38eae6e520680739" data-media-type="image/jpeg" alt="A graph in quadrant 1 of a generic function f(x). It is an increasing concave up function for the first quarter, an increasing concave down function for the second quarter, a decreasing concave down function for the third quarter, and an increasing concave down function for the last quarter. In the second quarter, a point a is marked on the x axis, and in the third quarter, a point b is marked on the x axis. The area under the curve and between a and b is shaded. This area is labeled the integral from a to b of f(x) dx." id="2">

Figure 1. The Fundamental Theorem of Calculus says that the integral over line segment [latex][a,b][/latex] depends only on the values of the antiderivative at the endpoints of [latex][a,b][/latex].

Green’s theorem takes this idea and extends it to calculating double integrals. Green’s theorem says that we can calculate a double integral over region [latex]D[/latex] based solely on information about the boundary of [latex]D[/latex] Green’s theorem also says we can calculate a line integral over a simple closed curve [latex]C[/latex] based solely on information about the region that [latex]C[/latex] encloses. In particular, Green’s theorem connects a double integral over region [latex]D[/latex] to a line integral around the boundary of [latex]D[/latex].

Circulation Form of Green’s Theorem

The first form of Green’s theorem that we examine is the circulation form. This form of the theorem relates the vector line integral over a simple, closed plane curve [latex]C[/latex] to a double integral over the region enclosed by [latex]C[/latex]. Therefore, the circulation of a vector field along a simple closed curve can be transformed into a double integral and vice versa.

theorem: green’s theorem, circulation form


Let [latex]D[/latex] be an open, simply connected region with a boundary curve [latex]C[/latex] that is a piecewise smooth, simple closed curve oriented counterclockwise (Figure 2). Let [latex]{\bf{F}}=\langle{P},Q\rangle[/latex] be a vector field with component functions that have continuous partial derivatives on [latex]D[/latex]. Then,

[latex]\displaystyle\oint_C{\bf{F}}\cdot{d}{\bf{r}}=\displaystyle\oint_CPdx+Qdy=\displaystyle\iint_D(Q_x-P_y)dA[/latex].

<img src="/apps/archive/20220422.171947/resources/36fc30dd3c8b4bbbf13c91ec20e65b9d6aa4a4f7" data-media-type="image/jpeg" alt="A vector field in two dimensions with all of the arrows pointing up and to the right. A curve C oriented counterclockwise sections off a region D around the origin. It is a simple, closed region." id="4">

Figure 2. The circulation form of Green’s theorem relates a line integral over curve [latex]C[/latex] to a double integral over region [latex]D[/latex].

Notice that Green’s theorem can be used only for a two-dimensional vector field [latex]{\bf{F}}[/latex]. If [latex]{\bf{F}}[/latex] is a three-dimensional field, then Green’s theorem does not apply. Since

[latex]\large{\displaystyle\int_CPdx+Qdy=\displaystyle\int_C{\bf{F}}\cdot{\bf{T}}ds}[/latex],

this version of Green’s theorem is sometimes referred to as the tangential form of Green’s theorem.

The proof of Green’s theorem is rather technical, and beyond the scope of this text. Here we examine a proof of the theorem in the special case that [latex]D[/latex] is a rectangle. For now, notice that we can quickly confirm that the theorem is true for the special case in which [latex]{\bf{F}}\langle{P},Q\rangle[/latex] is conservative. In this case,

[latex]\large{\displaystyle\oint_CPdx+Qdy=0}[/latex]

because the circulation is zero in conservative vector fields. By Cross-Partial Property of Conservative Fields Theorem, [latex]{\bf{F}}[/latex] satisfies the cross-partial condition, so [latex]P_y=Q_x[/latex]. Therefore,

[latex]\large{\displaystyle\iint_D(Q_x-P_y)dA=\displaystyle\iint_D0dA=0=\displaystyle\oint_CPdx+Qdy}[/latex],

which confirms Green’s theorem in the case of conservative vector fields.

Proof

Let’s now prove that the circulation form of Green’s theorem is true when the region [latex]D[/latex] is a rectangle. Let [latex]D[/latex] be the rectangle [latex][a,b]\times[c,d][/latex] oriented counterclockwise. Then, the boundary [latex]C[/latex] of [latex]D[/latex] consists of four piecewise smooth pieces [latex]C_1[/latex], [latex]C_2[/latex], [latex]C_3[/latex], and [latex]C_4[/latex] (Figure 3). We parameterize each side of [latex]D[/latex] as follows:

[latex]\begin{aligned} C_1:{\bf{r}}_1(t)&=\langle{t},c\rangle,a\leq{t}\leq{b} \\ C_2:{\bf{r}}_2(t)&=\langle{b},t\rangle,c\leq{t}\leq{d} \\ -C_3:{\bf{r}}_3(t)&=\langle{t},d\rangle,a\leq{t}\leq{b} \\ -C_4:{\bf{r}}_4(t)&=\langle{a},t\rangle,c\leq{t}\leq{d} \end{aligned}[/latex].

<img src="/apps/archive/20220422.171947/resources/0177a6156983d7186594e82b72926d649a2e18f1" data-media-type="image/jpeg" alt="A diagram in quadrant 1. Rectangle D is oriented counterclockwise. Points a and b are on the x axis, and points c and d are on the y axis with b > a and d > c. The sides of the rectangle are side c1 with endpoints at (a,c) and (b,c), side c2 with endpoints at (b,c) and (b,d), side c3 with endpoints at (b,d) and (a,d), and side c4 with endpoints at (a,d) and (a,c)." id="5">

Figure 3. Rectangle [latex]D[/latex] is oriented counterclockwise.

Then,

[latex]\begin{aligned} \displaystyle\int_C{\bf{F}}\bullet{d}{\bf{r}}&=\displaystyle\int_{C_1}{\bf{F}}\bullet{d}{\bf{r}}+\displaystyle\int_{C_2}{\bf{F}}\bullet{d}{\bf{r}}+\displaystyle\int_{C_3}{\bf{F}}\bullet{d}{\bf{r}}+\displaystyle\int_{C_4}{\bf{F}}\bullet{d}{\bf{r}} \\ &=\displaystyle\int_{C_1}{\bf{F}}\bullet{d}{\bf{r}}+\displaystyle\int_{C_2}{\bf{F}}\bullet{d}{\bf{r}}-\displaystyle\int_{-C_3}{\bf{F}}\bullet{d}{\bf{r}}-\displaystyle\int_{-C_4}{\bf{F}}\bullet{d}{\bf{r}} \\ &=\displaystyle\int_a^b{\bf{F}}({\bf{r}}_1(t))\bullet{\bf{r}}_1(t)dt+\displaystyle\int_c^d{\bf{F}}({\bf{r}}_2(t))\bullet{\bf{r}}_2(t)dt-\displaystyle\int_a^b{\bf{F}}({\bf{r}}_3(t))\bullet{\bf{r}}_3(t)dt-\displaystyle\int_c^d{\bf{F}}({\bf{r}}_4(t))\bullet{\bf{r}}_4(t)dt \\ &=\displaystyle\int_a^bP(t,c)dt+\displaystyle\int_c^dQ(b,t)dt-\displaystyle\int_a^bP(t,d)dt-\displaystyle\int_c^dQ(a,t)dt \\ &=\displaystyle\int_a^b(P(t,c)-P(t,d))dt+\displaystyle\int_c^d(Q(b,t)-Q(a,t))dt \\ &=-\displaystyle\int_a^b(P(t,d)-P(t,c))dt+\displaystyle\int_c^d(Q(b,t)-Q(a,t))dt \end{aligned}[/latex].

By the Fundamental Theorem of Calculus,

[latex]\large{P(t,d)-P(t,c)=\displaystyle\int_c^d\frac{\partial}{\partial{y}}P(t,y)dy\text{ and }Q(b,t)-Q(a,t)=\displaystyle\int_a^b\frac{\partial}{\partial{x}}Q(x,t)dx}[/latex].

Therefore,

[latex]\begin{aligned} &-\displaystyle\int_a^b(P(t,d)-P(t,c)dt+\displaystyle\int_c^d(Q(b,t)-Q(a,t))dt \\ &=-\displaystyle\int_a^b\displaystyle\int_c^d\frac{\partial}{\partial{y}}P(t,y)dydt+\displaystyle\int_c^d\displaystyle\int_a^b\frac{\partial}{\partial{x}}Q(x,t)dxdt \end{aligned}[/latex].

But,

[latex]\begin{aligned} -\displaystyle\int_a^b\displaystyle\int_c^d\frac{\partial}{\partial{y}}P(t,y)dydt+\displaystyle\int_c^d\displaystyle\int_a^b\frac{\partial}{\partial{x}}Q(x,t)dxdt&=-\displaystyle\int_a^b\displaystyle\int_c^d\frac{\partial}{\partial{y}}P(x,y)dydx+\displaystyle\int_c^d\displaystyle\int_a^b\frac{\partial}{\partial{x}}Q(x,t)dxdy \\ &=\displaystyle\int_a^b\displaystyle\int_c^d(Q_x-P_y)dydx \\ &=\displaystyle\int \ \displaystyle\int_D(Q_x-P_y)dA \end{aligned}[/latex].

Therefore, [latex]\displaystyle\int_C{\bf{F}}\bullet{d}{\bf{r}}=\int\int_D(Q_{x}-P_{y})dA[/latex] and we have proved Green’s theorem in the case of a rectangle.

To prove Green’s theorem over a general region [latex]D[/latex], we can decompose [latex]D[/latex] into many tiny rectangles and use the proof that the theorem works over rectangles. The details are technical, however, and beyond the scope of this text.

[latex]_\blacksquare[/latex]

Example: applying Green’s Theorem over a rectangle

Calculate the line integral

[latex]\displaystyle\oint_Cx^2ydx+(y-3)dy[/latex]

where [latex]C[/latex] is a rectangle with vertices [latex](1, 1)[/latex], [latex](4, 1)[/latex], [latex](4, 5)[/latex], and [latex](1, 5)[/latex] oriented counterclockwise.

Analysis

If we were to evaluate this line integral without using Green’s theorem, we would need to parameterize each side of the rectangle, break the line integral into four separate line integrals, and use the methods from Line Integrals to evaluate each integral. Furthermore, since the vector field here is not conservative, we cannot apply the Fundamental Theorem for Line Integrals. Green’s theorem makes the calculation much simpler.

 

Example: Applying Green’s Theorem to Calculate work

Calculate the work done on a particle by force field

[latex]{\bf{F}}(x,y)=\langle{y}+\sin{x},e^y-x\rangle[/latex]

as the particle traverses circle [latex]x^2+y^2=4[/latex] exactly once in the counterclockwise direction, starting and ending at point [latex](2, 0)[/latex].

try it

Use Green’s theorem to calculate line integral

[latex]\displaystyle\oint_C\sin(x^2)dx+(3x-y)dy[/latex],

where [latex]C[/latex] is a right triangle with vertices [latex](-1, 2)[/latex], [latex](4, 2)[/latex], and [latex](4, 5)[/latex] oriented counterclockwise.

Watch the following video to see the worked solution to the above Try It

You can view the transcript for “CP 6.34” here (opens in new window).
In the preceding two examples, the double integral in Green’s theorem was easier to calculate than the line integral, so we used the theorem to calculate the line integral. In the next example, the double integral is more difficult to calculate than the line integral, so we use Green’s theorem to translate a double integral into a line integral.

Example: applying green’s theorem over an ellipse

Calculate the area enclosed by ellipse [latex]\frac{x^2}{a^2}+\frac{y^2}{b^2}=1[/latex] (Figure 6).

<img src="/apps/archive/20220422.171947/resources/d7d97e8443810811ba4304de3fe2d855f57a3cae" data-media-type="image/jpeg" alt="A horizontal ellipse graphed in two dimensions. It has vertices at (-a, 0), (0, -b), (a, 0), and (0, b), where the absolute value of a is between 2.5 and 5 and the absolute value of b is between 0 and 2.5." id="13">

Figure 6. Ellipse [latex]\frac{x^2}{a^2}+\frac{y^2}{b^2}=1[/latex] is denoted by [latex]C[/latex].

In Example “Applying Green’s Theorem over an Ellipse”, we used vector field [latex]{\bf{F}}(x,y)=\langle{P},Q\rangle=\left\langle-\frac{y}2,\frac{x}2\right\rangle[/latex] to find the area of any ellipse. The logic of the previous example can be extended to derive a formula for the area of any region [latex]D[/latex]. Let [latex]D[/latex] be any region with a boundary that is a simple closed curve [latex]C[/latex] oriented counterclockwise. If [latex]{\bf{F}}(x,y)=\langle{P},Q\rangle=\left\langle-\frac{y}2,\frac{x}2\right\rangle[/latex]. Therefore, by the same logic as in Example “Applying Green’s Theorem over an Ellipse”,

area of [latex]D=\displaystyle\iint_DdA=\frac12\displaystyle\oint_C-ydx+xdy[/latex].

It’s worth noting that if [latex]{\bf{F}}=\langle{P},Q\rangle[/latex] is any vector field with [latex]Q_x-P_y=1[/latex], then the logic of the previous paragraph works. So, the equation above is not the only equation that uses a vector field’s mixed partials to get the area of a region.

try it

Find the area of the region enclosed by the curve with parameterization [latex]{\bf{r}}(t)=\langle\sin{t}\cos{t},\sin{t}\rangle,\text{ }0\leq{t}\leq\pi[/latex].

An image of a curve in quadrants 1 and 2. The curve begins at the origin, curves up and to the right until about (.5, .8), curves to the left nearly horizontally, goes through (0,1), continues until about (-1, .7), and then curves down and to the right until it hits the origin again.

Figure 7. This region is enclosed by the curve with parameterization [latex]{\bf{r}}(t)=\langle\sin{t}\cos{t},\sin{t}\rangle,\text{ }0\leq{t}\leq\pi[/latex].