Learning Objectives
- Calculate circulation and flux on more general regions.
Green’s Theorem on General Regions
Green’s theorem, as stated, applies only to regions that are simply connected—that is, Green’s theorem as stated so far cannot handle regions with holes. Here, we extend Green’s theorem so that it does work on regions with finitely many holes (Figure 1).
Before discussing extensions of Green’s theorem, we need to go over some terminology regarding the boundary of a region. Let [latex]D[/latex] be a region and let [latex]C[/latex] be a component of the boundary of [latex]D[/latex]. We say that [latex]C[/latex] is positively oriented if, as we walk along [latex]C[/latex] in the direction of orientation, region [latex]D[/latex] is always on our left. Therefore, the counterclockwise orientation of the boundary of a disk is a positive orientation, for example. Curve [latex]C[/latex] is negatively oriented if, as we walk along [latex]C[/latex] in the direction of orientation, region [latex]D[/latex] is always on our right. The clockwise orientation of the boundary of a disk is a negative orientation, for example.
Let [latex]D[/latex] be a region with finitely many holes (so that [latex]D[/latex] has finitely many boundary curves), and denote the boundary of [latex]D[/latex] by [latex]\partial{D}[/latex] (Figure 2). To extend Green’s theorem so it can handle [latex]D[/latex], we divide region [latex]D[/latex] into two regions, [latex]D_1[/latex] and [latex]D_2[/latex] (with respective boundaries [latex]\partial{D_1}[/latex] and [latex]\partial{D_2}[/latex]), in such a way that [latex]D=D_1\cup{D}_2[/latex] and neither [latex]D_1[/latex] nor [latex]D_2[/latex] has any holes (Figure 2).
Assume the boundary of [latex]D[/latex] is oriented as in the figure, with the inner holes given a negative orientation and the outer boundary given a positive orientation. The boundary of each simply connected region [latex]D_1[/latex] and [latex]D_2[/latex] is positively oriented. If [latex]{\bf{F}}[/latex] is a vector field defined on [latex]D[/latex], then Green’s theorem says that
[latex]\begin{aligned} \displaystyle\oint_{\partial{D}}{\bf{F}}\cdot{d}{\bf{r}}&=\displaystyle\oint_{\partial{D_1}}{\bf{F}}\cdot{d}{\bf{r}}+\displaystyle\oint_{\partial{D_2}}{\bf{F}}\cdot{d}{\bf{r}} \\ &=\displaystyle\iint_{D_1}Q_x-P_ydA+\displaystyle\iint_{D_2}Q_x-P_ydA \\ &=\displaystyle\iint_{D}(Q_x-P_y)dA \end{aligned}[/latex].
Therefore, Green’s theorem still works on a region with holes.
To see how this works in practice, consider annulus [latex]D[/latex] in Figure 3 and suppose that [latex]{\bf{F}}=\langle{P},Q\rangle[/latex] is a vector field defined on this annulus. Region [latex]D[/latex] has a hole, so it is not simply connected. Orient the outer circle of the annulus counterclockwise and the inner circle clockwise (Figure 3) so that, when we divide the region into [latex]D_1[/latex] and [latex]D_2[/latex], we are able to keep the region on our left as we walk along a path that traverses the boundary. Let [latex]D_1[/latex] be the upper half of the annulus and [latex]D_2[/latex] be the lower half. Neither of these regions has holes, so we have divided [latex]D[/latex] into two simply connected regions.
We label each piece of these new boundaries as [latex]P_i[/latex] for some [latex]i[/latex], as in Figure 3. If we begin at [latex]P[/latex] and travel along the oriented boundary, the first segment is [latex]P_1[/latex], then [latex]P_2[/latex], [latex]P_3[/latex], and [latex]P_4[/latex]. Now we have traversed [latex]D_1[/latex] and returned to [latex]P[/latex] Next, we start at [latex]P[/latex] again and traverse [latex]D_2[/latex]. Since the first piece of the boundary is the same as [latex]P_4[/latex] in [latex]D_1[/latex], but oriented in the opposite direction, the first piece of [latex]D_2[/latex] is [latex]-P_4[/latex]. Next, we have [latex]P_5[/latex], then [latex]-P_2[/latex], and finally [latex]P_6[/latex].
Figure 3 shows a path that traverses the boundary of [latex]D[/latex]. Notice that this path traverses the boundary of region [latex]D_1[/latex], returns to the starting point, and then traverses the boundary of region [latex]D_2[/latex]. Furthermore, as we walk along the path, the region is always on our left. Notice that this traversal of the [latex]P_i[/latex] paths covers the entire boundary of region [latex]D[/latex]. If we had only traversed one portion of the boundary of [latex]D[/latex], then we cannot apply Green’s theorem to [latex]D[/latex].
The boundary of the upper half of the annulus, therefore, is [latex]P_1\cup{P_2}\cup{P_3}\cup{P_4}[/latex] and the boundary of the lower half of the annulus is [latex]-P_4\cup{P_5}\cup-{P_2}\cup{P_6}[/latex]. Then, Green’s theorem implies
[latex]\begin{aligned} \displaystyle\int_{\partial{D}}{\bf{F}}\cdot{d}{\bf{r}}&=\displaystyle\int+{P_1}{\bf{F}}\cdot{d}{\bf{r}}+\displaystyle\int{P_2}{\bf{F}}\cdot{d}{\bf{r}}+\displaystyle\int{P_3}{\bf{F}}\cdot{d}{\bf{r}}+\displaystyle\int{P_4}{\bf{F}}\cdot{d}{\bf{r}}+\displaystyle\int{-P_4}{\bf{F}}\cdot{d}{\bf{r}}+\displaystyle\int{P_5}{\bf{F}}\cdot{d}{\bf{r}}-\displaystyle\int{P_2}{\bf{F}}\cdot{d}{\bf{r}}+\displaystyle\int{P_6}{\bf{F}}\cdot{d}{\bf{r}} \\ &=\displaystyle\int+{P_1}{\bf{F}}\cdot{d}{\bf{r}}+\displaystyle\int{P_2}{\bf{F}}\cdot{d}{\bf{r}}+\displaystyle\int{P_3}{\bf{F}}\cdot{d}{\bf{r}}+\displaystyle\int{P_4}{\bf{F}}\cdot{d}{\bf{r}}-\displaystyle\int{P_4}{\bf{F}}\cdot{d}{\bf{r}}+\displaystyle\int{P_5}{\bf{F}}\cdot{d}{\bf{r}}-\displaystyle\int{P_2}{\bf{F}}\cdot{d}{\bf{r}}+\displaystyle\int{P_6}{\bf{F}}\cdot{d}{\bf{r}} \\ &=\displaystyle\int+{P_1}{\bf{F}}\cdot{d}{\bf{r}}+\displaystyle\int{P_3}{\bf{F}}\cdot{d}{\bf{r}}+\displaystyle\int{P_5}{\bf{F}}\cdot{d}{\bf{r}}+\displaystyle\int{P_6}{\bf{F}}\cdot{d}{\bf{r}} \\ &=\displaystyle\int_{\partial{D}_1}{\bf{F}}\cdot{d}{\bf{r}}+\displaystyle\int_{\partial{D}_2}{\bf{F}}\cdot{d}{\bf{r}} \\ &=\displaystyle\iint_{D_1}(Q_x-P_y)dA+\displaystyle\iint_{D_2}(Q_x-P_y)dA \\ &=\displaystyle\iint_{D}(Q_x-P_y)dA \end{aligned}[/latex].
Therefore, we arrive at the equation found in Green’s theorem—namely,
[latex]\large{\displaystyle\oint_{\partial{D}}{\bf{F}}\cdot{d}{\bf{r}}=\displaystyle\iint_D(Q_z-P_y)dA}[/latex].
The same logic implies that the flux form of Green’s theorem can also be extended to a region with finitely many holes:
[latex]\large{\displaystyle\oint_{C}{\bf{F}}\cdot{\bf{D}}ds=\displaystyle\iint_D(P_x+Q_y)dA}[/latex].
Example: using green’s theorem on a region with holes
Calculate integral
[latex]\large{\displaystyle\oint_{\partial{D}}\left(\sin{x}-\frac{x^3}3\right)dx+\left(\frac{y^3}3+\sin{y}\right)dy}[/latex],
where [latex]D[/latex] is the annulus given by the polar inequalities [latex]1\leq{\bf{r}}\leq2,0\leq \theta\leq2\pi[/latex].
Example: using the extended form of green’s theorem
Let [latex]{\bf{F}}=\langle{P},Q\rangle=\langle\frac{y}{x^2+y^2},-\frac{y}{x^2+y^2}[/latex] and let [latex]C[/latex] be any simple closed curve in a plane oriented counterclockwise. What are the possible values of [latex]\displaystyle\oint_C{\bf{F}}\cdot{d}{\bf{r}}[/latex]?
try it
Calculate integral [latex]\displaystyle\oint_{\partial{D}}{\bf{F}}\cdot{d}{\bf{r}}[/latex], where [latex]D[/latex] is the annulus given by the polar inequalities [latex]2\leq{r}\leq5,0\leq\theta\leq2\pi[/latex], and [latex]{\bf{F}}(x,y)=\langle{x}^3,5x+e^y\sin{y}\rangle[/latex].
Watch the following video to see the worked solution to the above Try It
Activity: measuring area from a boundary: the planimeter
Imagine you are a doctor who has just received a magnetic resonance image of your patient’s brain. The brain has a tumor (Figure 5). How large is the tumor? To be precise, what is the area of the red region? The red cross-section of the tumor has an irregular shape, and therefore it is unlikely that you would be able to find a set of equations or inequalities for the region and then be able to calculate its area by conventional means. You could approximate the area by chopping the region into tiny squares (a Riemann sum approach), but this method always gives an answer with some error.
Instead of trying to measure the area of the region directly, we can use a device called a rolling planimeter to calculate the area of the region exactly, simply by measuring its boundary. In this project you investigate how a planimeter works, and you use Green’s theorem to show the device calculates area correctly.
A rolling planimeter is a device that measures the area of a planar region by tracing out the boundary of that region (Figure 6). To measure the area of a region, we simply run the tracer of the planimeter around the boundary of the region. The planimeter measures the number of turns through which the wheel rotates as we trace the boundary; the area of the shape is proportional to this number of wheel turns. We can derive the precise proportionality equation using Green’s theorem. As the tracer moves around the boundary of the region, the tracer arm rotates and the roller moves back and forth (but does not rotate).
Let [latex]C[/latex] denote the boundary of region [latex]D[/latex], the area to be calculated. As the tracer traverses curve [latex]C[/latex], assume the roller moves along the y-axis (since the roller does not rotate, one can assume it moves along a straight line). Use the coordinates [latex](x, y)[/latex] to represent points on boundary [latex]C[/latex], and coordinates [latex](0, Y)[/latex] to represent the position of the pivot. As the planimeter traces [latex]C[/latex], the pivot moves along the y-axis while the tracer arm rotates on the pivot.
Media
Watch a short animation of a planimeter in action.
Begin the analysis by considering the motion of the tracer as it moves from point [latex](x, y)[/latex] counterclockwise to point [latex](x+dx, y+dy)[/latex] that is close to [latex](x, y)[/latex] (Figure 7). The pivot also moves, from point [latex](0, Y)[/latex] to nearby point [latex]( 0, Y+dY)[/latex]. How much does the wheel turn as a result of this motion? To answer this question, break the motion into two parts. First, roll the pivot along the y-axis from [latex](0, Y)[/latex] to [latex](0, Y+dY)[/latex] without rotating the tracer arm. The tracer arm then ends up at point [latex](x,y+dY)[/latex] while maintaining a constant angle [latex]\phi[/latex] with the x-axis. Second, rotate the tracer arm by an angle [latex]d\theta[/latex] without moving the roller. Now the tracer is at point [latex](x+dx,y+dy)[/latex]. Let [latex]l[/latex] be the distance from the pivot to the wheel and let [latex]L[/latex] be the distance from the pivot to the tracer (the length of the tracer arm).
- Explain why the total distance through which the wheel rolls the small motion just described is [latex]\sin\phi{d}Y+ld\theta=\frac{x}LdY+ld\theta[/latex].
- Show that [latex]\displaystyle\oint_Cd\theta=0[/latex].
- Use step 2 to show that the total rolling distance of the wheel as the tracer traverses curve [latex]C[/latex] is
Total wheel roll [latex]=\frac1L\displaystyle\oint_CxdY[/latex].
Now that you have an equation for the total rolling distance of the wheel, connect this equation to Green’s theorem to calculate area [latex]D[/latex] enclosed by [latex]C[/latex]. - Show that [latex]x^2+(y-Y)^2=L^2[/latex].
- Assume the orientation of the planimeter is as shown in Figure 7. Explain why [latex]Y\leq{y}[/latex], and use this inequality to show there is a unique value of [latex]Y[/latex] for each point [latex](x, y)[/latex]: [latex]Y=y=\sqrt{L^2-x^2}[/latex].
- Use step 5 to show that [latex]dY=dy+\frac{x}{\sqrt{L^2-x^2}}dx[/latex].
- Use Green’s theorem to show that [latex]\displaystyle\oint_C\frac{x}{\sqrt{L^2-x^2}}dx=0[/latex].
- Use step 7 to show that the total wheel roll is
Total wheel roll [latex]=\frac1L\displaystyle\oint_CxdY[/latex].
It took a bit of work, but this equation says that the variable of integration [latex]Y[/latex] in step 3 can be replaced with [latex]y[/latex]. - Use Green’s theorem to show that the area of [latex]D[/latex] is [latex]\displaystyle\oint_Cxdy[/latex]. The logic is similar to the logic used to show that the area of [latex]D=\frac12\displaystyle\oint_C-ydx+xdy[/latex].
- Conclude that the area of [latex]D[/latex] equals the length of the tracer arm multiplied by the total rolling distance of the wheel.
You now know how a planimeter works and you have used Green’s theorem to justify that it works. To calculate the area of a planar region [latex]D[/latex], use a planimeter to trace the boundary of the region. The area of the region is the length of the tracer arm multiplied by the distance the wheel rolled.