Learning Objectives
- Apply the flux form of Green’s theorem.
The circulation form of Green’s theorem relates a double integral over region [latex]D[/latex] to line integral [latex]\displaystyle\oint_C{\bf{F}}\cdot{\bf{T}}ds[/latex], where [latex]C[/latex] is the boundary of [latex]D[/latex]. The flux form of Green’s theorem relates a double integral over region [latex]D[/latex] to the flux across boundary [latex]C[/latex]. The flux of a fluid across a curve can be difficult to calculate using the flux line integral. This form of Green’s theorem allows us to translate a difficult flux integral into a double integral that is often easier to calculate.
theorem: Green’s theorem, flux 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 that is oriented counterclockwise (Figure 1). Let [latex]{\bf{F}}=\langle{P},Q\rangle[/latex] be a vector field with component functions that have continuous partial derivatives on an open region containing [latex]D[/latex]. Then,
[latex]\displaystyle\oint{\bf{F}}\cdot{\bf{N}}ds=\displaystyle\iint_DP_x+Q_ydA[/latex].
Figure 1. The flux form of Green’s theorem relates a double integral over region [latex]D[/latex] to the flux across curve [latex]C[/latex].
Because this form of Green’s theorem contains unit normal vector [latex]\bf{N}[/latex], it is sometimes referred to as the normal form of Green’s theorem.
Proof
Recall that [latex]\displaystyle\oint{\bf{F}}\cdot{\bf{N}}ds=\displaystyle\oint_C-Qdx+Pdy[/latex]. Let [latex]M=-Q[/latex] and [latex]N=P[/latex]. By the circulation form of Green’s theorem,
[latex]\begin{aligned} \displaystyle\oint_C-Qdx+Pdy&=\displaystyle\oint_CMdx+nDy \\ &=\displaystyle\iint_DN_x-M_ydA \\ &=\displaystyle\iint_DP_x-(-Q)_ydA \\ &=\displaystyle\iint_DP_x+Q_ydA \end{aligned}[/latex].
[latex]_\blacksquare[/latex]
Example: Applying Green’s theorem for flux across a circle
Let [latex]C[/latex] be a circle of radius [latex]r[/latex] centered at the origin (Figure 2) and let [latex]{\bf{F}}(x,y)=\langle{x},y\rangle[/latex]. Calculate the flux across [latex]C[/latex].
Figure 2. Curve [latex]C[/latex] is a circle of radius [latex]r[/latex] centered at the origin.
Show Solution
Let [latex]D[/latex] be the disk enclosed by [latex]C[/latex]. The flux across [latex]C[/latex] is [latex]\displaystyle\oint_C{\bf{F}}\cdot{\bf{N}}ds[/latex]. We could evaluate this integral using tools we have learned, but Green’s theorem makes the calculation much more simple. Let [latex]P(x, y)=x[/latex] and [latex]Q(x, y)=y[/latex] so that [latex]{\bf{F}}=\langle{P},Q\rangle[/latex]. Note that [latex]P_x=1=Q_y[/latex], and therefore [latex]P_x+Q_y=2[/latex]. By Green’s theorem,
[latex]\displaystyle\int_C{\bf{F}}\bullet{\bf{N}}=\displaystyle\int \ \displaystyle\int_D2dA=2\displaystyle\int \ \displaystyle\int_DdA[/latex].
Since [latex]\displaystyle\int_DdA[/latex] is the area of the circle, [latex]\displaystyle\int_DdA=\pi{r}^2[/latex]. Therefore, the flux across [latex]C[/latex] is [latex]2\pi{r}[/latex].
Example: applying green’s theorem for flux across a triangle
Let [latex]S[/latex] be the triangle with vertices [latex](0, 0)[/latex], [latex](1, 0)[/latex], and [latex](0, 3)[/latex] oriented clockwise (Figure 3). Calculate the flux of [latex]{\bf{F}}(x,y)=\langle{P}(x,y),Q(x,y)\rangle=\langle{x}^2+e^y,x+y\rangle[/latex] across [latex]S[/latex].
Figure 3. Curve [latex]S[/latex] is a triangle with vertices [latex](0, 0)[/latex], [latex](1, 0)[/latex], and [latex](0, 3)[/latex] oriented clockwise.
Show Solution
To calculate the flux without Green’s theorem, we would need to break the flux integral into three line integrals, one integral for each side of the triangle. Using Green’s theorem to translate the flux line integral into a single double integral is much more simple.
Let [latex]D[/latex] be the region enclosed by [latex]S[/latex]. Note that [latex]P_x=2x[/latex] and [latex]Q_y=1[/latex] therefore, [latex]P_x+Q_y=2x+1[/latex]. Green’s theorem applies only to simple closed curves oriented counterclockwise, but we can still apply the theorem because [latex]\displaystyle\oint_C{\bf{F}}\cdot{\bf{N}}ds=-\displaystyle\oint_{-S}{\bf{F}}\cdot{\bf{N}}ds[/latex] and [latex]-S[/latex] is oriented counterclockwise. By Green’s theorem, the flux is
[latex]\begin{aligned} \displaystyle\oint_C{\bf{F}}\cdot{\bf{N}}ds&=\displaystyle\oint_{-S}{\bf{F}}\cdot{\bf{N}}ds \\ &=-\displaystyle\iint_D(P_x+Q_y)dA \\ &=-\displaystyle\iint_D(2x+1)dA \end{aligned}[/latex]
Notice that the top edge of the triangle is the line [latex]y=-3x=3[/latex] Therefore, in the iterated double integral, the y-values run from [latex]y=0[/latex] to [latex]y=-3x+3[/latex], and we have
[latex]\begin{aligned} -\displaystyle\iint_D(2x+1)dA&=-\displaystyle\int_0^1\displaystyle\int_0^{-3x+3}(2x+1)dydx \\ &=--\displaystyle\int_0^1(2x+1)(-3x+3)dx=-\displaystyle\int_0^1(-6x^2+3x+3)dx \\ &=-\left[-2x^3+\frac{3x^2}2+3x\right]_0^1=-\frac52 \end{aligned}[/latex].
try it
Calculate the flux of [latex]{\bf{F}}(x,y)=\langle{x}^3,y^3\rangle[/latex] across a unit circle oriented counterclockwise.
Show Solution
[latex]\frac{3\pi}2[/latex].
Watch the following video to see the worked solution to the above Try It
Example: applying green’s theorem for water flow across a rectangle
Water flows from a spring located at the origin. The velocity of the water is modeled by vector field [latex]{\bf{v}}(x,y)=\langle5x+y,x+3y\rangle[/latex] m/sec. Find the amount of water per second that flows across the rectangle with vertices [latex](-1, -2)[/latex], [latex](1, -2)[/latex], [latex](1, 3)[/latex], and [latex](-1, 3)[/latex], oriented counterclockwise (Figure 4).
Figure 4. Water flows across the rectangle with vertices [latex](-1, -2)[/latex], [latex](1, -2)[/latex], [latex](1, 3)[/latex], and [latex](-1, 3)[/latex], oriented counterclockwise.
Show Solution
Let [latex]C[/latex] epresent the given rectangle and let [latex]D[/latex] be the rectangular region enclosed by [latex]C[/latex]. To find the amount of water flowing across [latex]C[/latex], we calculate flux [latex]\displaystyle\int_C{\bf{v}}\cdot{\bf{N}}ds[/latex].
Let [latex]P(x, y)=5x+y[/latex] and [latex]Q(x, y)=x+3y[/latex] so that [latex]{\bf{v}}=(P,Q)[/latex]. Then, [latex]P_x=5[/latex] and [latex]Q_y=3[/latex]. By Green’s theorem,
[latex]\begin{aligned} \displaystyle\int_C{\bf{v}}\cdot{\bf{N}}ds&=\displaystyle\iint_D(P_x+Q_y)dA \\ &=\displaystyle\iint_D8dA \\ &=8(\text{area of }D)=80 \end{aligned}[/latex].
Therefore, the water flux is [latex]80[/latex] m2/sec.
Recall that if vector field [latex]{\bf{F}}[/latex] is conservative, then [latex]{\bf{F}}[/latex] does no work around closed curves—that is, the circulation of [latex]{\bf{F}}[/latex] around a closed curve is zero. In fact, if the domain of [latex]{\bf{F}}[/latex] is simply connected, then [latex]{\bf{F}}[/latex] is conservative if and only if the circulation of [latex]{\bf{F}}[/latex] around any closed curve is zero. If we replace “circulation of [latex]{\bf{F}}[/latex]” with “flux of [latex]{\bf{F}}[/latex],” then we get a definition of a source-free vector field. The following statements are all equivalent ways of defining a source-free field [latex]{\bf{F}}=\langle{P},Q\rangle[/latex] on a simply connected domain (note the similarities with properties of conservative vector fields):
- The flux [latex]\displaystyle\oint_C{\bf{F}}\cdot{\bf{N}}ds[/latex] across any closed curve [latex]C[/latex] is zero.
- If [latex]C_1[/latex] and [latex]C_2[/latex] are curves in the domain of [latex]{\bf{F}}[/latex] with the same starting points and endpoints, then [latex]\displaystyle\int_{C_1}{\bf{F}}\cdot{\bf{N}}ds=\displaystyle\int_{C_2}{\bf{F}}\cdot{\bf{N}}ds[/latex]. In other words, flux is independent of path.
- There is a stream function [latex]g(x, y)[/latex] for [latex]{\bf{F}}[/latex]. A stream function for [latex]{\bf{F}}=\langle{P},Q\rangle[/latex] is a function [latex]g[/latex] such that [latex]P=g_y[/latex] and [latex]Q=-g_x[/latex]. Geometrically, [latex]{\bf{F}}(a,b)[/latex] is tangential to the level curve of [latex]g[/latex] at [latex](a, b)[/latex]. Since the gradient of [latex]g[/latex] is perpendicular to the level curve of [latex]g[/latex] at [latex](a, b)[/latex], stream function [latex]g[/latex] has the property [latex]{\bf{F}}(a,b)\bullet\nabla{g}(a,b)=0[/latex] for any point [latex](a, b)[/latex] in the domain of [latex]g[/latex]. (Stream functions play the same role for source-free fields that potential functions play for conservative fields.)
- [latex]P_x+Q_y=0[/latex]
Example: finding a stream function
Verify that rotation vector field [latex]{\bf{F}}(x,y)=\langle{y},-x\rangle[/latex] is source free, and find a stream function for [latex]{\bf{F}}[/latex].
Show Solution
Note that the domain of [latex]{\bf{F}}[/latex] is all of [latex]\mathbb{R}^2[/latex], which is simply connected. Therefore, to show that [latex]{\bf{F}}[/latex] is source free, we can show any of items 1 through 4 from the previous list to be true. In this example, we show that item 4 is true. Let [latex]P(x, y)=y[/latex] and [latex]Q(x, y)=-x[/latex]. Then [latex]P_x+0=Q_y[/latex], and therefore [latex]P_x+Q_y=0[/latex]. Thus, [latex]{\bf{F}}[/latex] is source free.
To find a stream function for [latex]{\bf{F}}[/latex], proceed in the same manner as finding a potential function for a conservative field. Let [latex]g[/latex] be a stream function for [latex]{\bf{F}}[/latex]. Then [latex]g_y=y[/latex], which implies that
[latex]\large{g(x,y)=\frac{y^2}2+h(x)}[/latex].
Since [latex]-g_x=Q=-x[/latex] we have [latex]h'(x)=x[/latex]. Therefore,
[latex]\large{h(x)=\frac{x^2}2+C}[/latex].
Letting [latex]C=0[/latex] gives stream function
[latex]\large{g(x,y)=\frac{x^2}2+\frac{y^2}2}[/latex].
To confirm that [latex]g[/latex] is a stream function for [latex]{\bf{F}}[/latex], note that [latex]g_y=y=P[/latex] and [latex]-g_x=-x=Q[/latex].
Notice that source-free rotation vector field [latex]{\bf{F}}(x,y)=\langle{y},-x\rangle[/latex] is perpendicular to conservative radial vector field [latex]\nabla{g}=\langle{x},y\rangle[/latex] (Figure 5).
Figure 5. (a) In this image, we see the three-level curves of [latex]g[/latex] and vector field [latex]{\bf{F}}[/latex]. Note that the [latex]{\bf{F}}[/latex] vectors on a given level curve are tangent to the level curve. (b) In this image, we see the three-level curves of [latex]g[/latex] and vector field [latex]\nabla{g}[/latex]. The gradient vectors are perpendicular to the corresponding level curve. Therefore, [latex]{\bf{F}}(a,b)\bullet\nabla{g}(a,b)=0[/latex] for any point in the domain of [latex]g[/latex].
try it
Find a stream function for vector field [latex]{\bf{F}}(x,y)=\langle{x}\sin{y},\cos{y}\rangle[/latex].
Show Solution
[latex]g(x,y)=-x\cos{y}[/latex].
Vector fields that are both conservative and source free are important vector fields. One important feature of conservative and source-free vector fields on a simply connected domain is that any potential function [latex]f[/latex] of such a field satisfies Laplace’s equation [latex]f_{xx}+f{yy}=0[/latex]. Laplace’s equation is foundational in the field of partial differential equations because it models such phenomena as gravitational and magnetic potentials in space, and the velocity potential of an ideal fluid. A function that satisfies Laplace’s equation is called a harmonic function. Therefore any potential function of a conservative and source-free vector field is harmonic.
To see that any potential function of a conservative and source-free vector field on a simply connected domain is harmonic, let [latex]f[/latex] be such a potential function of vector field [latex]{\bf{F}}=\langle{P},Q\rangle[/latex]. Then, [latex]f_x=P[/latex] and [latex]f_x=Q[/latex] because [latex]\nabla{f}={\bf{F}}[/latex]. Therefore, [latex]f_{xx}=P_x[/latex] and [latex]f_{yy}=[/latex]. Since [latex]{\bf{F}}[/latex] is source free, [latex]f_{xx}+f_{yy}=P_x+Q_y=0[/latex], and we have that [latex]f[/latex] is harmonic.
Example: satisfying laplace’s equation
For vector field [latex]{\bf{F}}(x,y)=\langle{e}^x\sin{y},e^x\cos{y}\rangle[/latex], verify that the field is both conservative and source free, find a potential function for [latex]{\bf{F}}[/latex], and verify that the potential function is harmonic.
Show Solution
Let [latex]P(x,y)={e}^x\sin{y}[/latex] and [latex]Q(x,y)=e^x\cos{y}[/latex]. Notice that the domain of [latex]{\bf{F}}[/latex] is all of two-space, which is simply connected. Therefore, we can check the cross-partials of [latex]{\bf{F}}[/latex] to determine whether [latex]{\bf{F}}[/latex] is conservative. Note that [latex]P_y=e^x\cos{y}[/latex], so [latex]{\bf{F}}[/latex] is conservative. Since [latex]P_x={e}^x\sin{y}[/latex] and [latex]Q_y={e}^x\sin{y}, P_x+Q_y=0[/latex] and the field is source free.
To find a potential function for [latex]{\bf{F}}[/latex], let [latex]f[/latex] be a potential function. Then, [latex]\nabla{f}={\bf{F}}[/latex], so [latex]f_x={e}^x\sin{y}[/latex]. Integrating this equation with respect to [latex]x[/latex] gives [latex]f(x,y)={e}^x\sin{y}+h(y)[/latex]. Since [latex]f_y=e^x\cos{y}[/latex], differentiating [latex]f[/latex] with respect to [latex]y[/latex] gives [latex]e^x\cos{y}=e^x\cos{y}+h^\prime(y)[/latex]. Therefore, we can take [latex]h(y)=0[/latex], and [latex]f(x,y)={e}^x\sin{y}[/latex] is a potential function for [latex]f[/latex].
To verify that [latex]f[/latex] is a harmonic function, note that [latex]f_{xx}=\frac{\partial}{\partial{x}}({e}^x\sin{y})={e}^x\sin{y}[/latex] and [latex]f_{yy}=\frac{\partial}{\partial{x}}(e^x\cos{y})=-{e}^x\sin{y}[/latex]. Therefore, [latex]f_{xx}+f_{yy}=0[/latex], and [latex]f[/latex] satisfies Laplace’s equation.
try it
Is the function [latex]f(x,y)=e^{x+5y}[/latex] harmonic?