15 Vector Analysis

15.4 Flow, Flux, Green’s Theorem and the Divergence Theorem

Flow and Flux

margin:

C2

C3

C1

1

1

x

y
Figure 15.4.1: Illustrating the principles of flow and flux.

Line integrals over vector fields have the natural interpretation of computing work when F represents a force field. It is also common to use vector fields to represent velocities. In these cases, the line integral CF𝑑r is said to represent flow.

Let the vector field F=1,0 represent the velocity of water as it moves across a smooth surface, depicted in Figure 15.4.1. A line integral over C will compute “how much water is moving along the path C.”

In the figure, “all” of the water above C1 is moving along that curve, whereas “none” of the water above C2 is moving along that curve (the curve and the flow of water are at right angles to each other). Because C3 has nonzero horizontal and vertical components, “some” of the water above that curve is moving along the curve.

When C is a closed curve, we call flow circulation, represented by CF𝑑r.

The “opposite” of flow is flux, a measure of “how much water is moving across the path C.” If a curve represents a filter in flowing water, flux measures how much water will pass through the filter. Considering again Figure 15.4.1, we see that a screen along C1 will not filter any water as no water passes across that curve. Because of the nature of this field, C2 and C3 each filter the same amount of water per second.

The terms “flow” and “flux” are used apart from velocity fields, too. Flow is measured by CF𝑑r, which is the same as CFT𝑑s by Definition 15.3.1. That is, flow is a summation of the amount of F that is tangent to the curve C.

By contrast, flux is a summation of the amount of F that is orthogonal to the direction of travel. To capture this orthogonal amount of F, we use CFn𝑑s to measure flux, where n is a unit vector orthogonal to the curve C. (Later, we’ll measure flux across surfaces, too. For example, in physics it is useful to measure the amount of a magnetic field that passes through a surface.)

How is n determined? We’ll later see that if C is a closed curve, we’ll want n to point to the outside of the curve (measuring how much is “going out”). We’ll also adopt the convention that closed curves should be traversed counterclockwise.

margin:

A

-1

1

-1

1

x

y
Figure 15.4.2: Determining “counterclockwise” is not always simple without a good definition.

(If C is a complicated closed curve, it can be difficult to determine what “counterclockwise” means. Consider Figure 15.4.2. Seeing the curve as a whole, we know which way “counterclockwise” is. If we zoom in on point A, one might incorrectly choose to traverse the path in the wrong direction. So we offer this definition: a closed curve is being traversed counterclockwise if the outside is to the right of the path and the inside is to the left.)

When a curve C is traversed counterclockwise by r(t)=f(t),g(t), we rotate T clockwise 90 to obtain n:

T=f(t),g(t)r(t)n=g(t),-f(t)r(t).

Letting F=M,N, we calculate flux as:

CFn𝑑s =CFg(t),-f(t)r(t)r(t)𝑑t
=CM,Ng(t),-f(t)𝑑t
=C(Mg(t)-Nf(t))𝑑t
=CMg(t)𝑑t-CNf(t)𝑑t.
As the x and y components of r(t) are f(t) and g(t) respectively, the differentials of x and y are dx=f(t)dt and dy=g(t)dt. We can then write the above integrals as:
=CM𝑑y-CN𝑑x.
This is often written as one integral (not incorrectly, though somewhat confusingly, as this one integral has two “d ’s”):
=CM𝑑y-Ndx.

We summarize the above in the following definition.

Definition 15.4.1      Flow, Flux

Let F=M,N be a vector field with continuous components defined on a smooth curve C, parameterized by r(t)=f(t),g(t), let T be the unit tangent vector of r(t), and let n be the clockwise 90degree rotation of T.

  • The flow of F along C is

    CFT𝑑s=CF𝑑r.
  • The flux of F across C is

    CFn𝑑s=CM𝑑y-Ndx=C(Mg(t)-Nf(t))𝑑t.

This definition of flow also holds for curves in space, though it does not make sense to measure “flux across a curve” in space.

Measuring flow is essentially the same as finding work performed by a force as done in the previous examples. Therefore we practice finding only flux in the following example.

margin:

C1

C2

1

1

x

y
(a)

C1

C2

1

1

x

y
(b)
Figure 15.4.3: Illustrating the curves and vector fields in Example 15.4.1. In (a) the vector field is F1, and in (b) the vector field is F2.
Example 15.4.1 Finding flux across curves in the plane

Curves C1 and C2 each start at (1,0) and end at (0,1), where C1 follows the line y=1-x and C2 follows the unit circle, as shown in Figure 15.4.3. Find the flux across both curves for the vector fields F1=y,-x+1 and F2=-x,2y-x.

SolutionWe begin by finding parameterizations of C1 and C2. As done in Example 15.3.3, parameterize C1 by creating the line that starts at (1,0) and moves in the -1,1 direction: r1(t)=1,0+t-1,1=1-t,t, for 0t1. We parameterize C2 with the familiar r2(t)=cost,sint on 0tπ/2. For reference later, we give each function and its derivative below:

r1(t)=1-t,t,r1(t)=-1,1.
r2(t)=cost,sint,r2(t)=-sint,cost.

When F=F1=y,-x+1 (as shown in Figure 15.4.3(a)), over C1 we have M=y=t and N=-x+1=-(1-t)+1=t. Using Definition 15.4.1, we compute the flux:

C1Fn𝑑s =C1(Mg(t)-Nf(t))𝑑t
=01(t(1)-t(-1))𝑑t
=012t𝑑t
=1.

Over C2, we have M=y=sint and N=-x+1=1-cost. Thus the flux across C2 is:

C1Fn𝑑s =C1(Mg(t)-Nf(t))𝑑t
=0π/2((sint)(cost)-(1-cost)(-sint))𝑑t
=0π/2sintdt
=1.

Notice how the flux was the same across both curves. This won’t hold true when we change the vector field.

When F=F2=-x,2y-x (as shown in Figure 15.4.3(b)), over C1 we have M=-x=t-1 and N=2y-x=2t-(1-t)=3t-1. Computing the flux across C1:

C1Fn𝑑s =C1(Mg(t)-Nf(t))𝑑t
=01((t-1)(1)-(3t-1)(-1))𝑑t
=01(4t-2)𝑑t
=0.

Over C2, we have M=-x=-cost and N=2y-x=2sint-cost. Thus the flux across C2 is:

C1Fn𝑑s =C1(Mg(t)-Nf(t))𝑑t
=0π/2((-cost)(cost)-(2sint-cost)(-sint))𝑑t
=0π/2(2sin2t-sintcost-cos2t)𝑑t
=π/4-1/20.285.

We analyze the results of this example below.

In Example 15.4.1, we saw that the flux across the two curves was the same when the vector field was F1=y,-x+1. This is not a coincidence. We show why they are equal in Example 15.4.6. In short, the reason is this: the divergence of F1 is 0, and when divF=0, the flux across any two paths with common beginning and ending points will be the same.

We also saw in the example that the flux across C1 was 0 when the field was F2=-x,2y-x. Flux measures “how much” of the field crosses the path from left to right (following the conventions established before). Positive flux means most of the field is crossing from left to right; negative flux means most of the field is crossing from right to left; zero flux means the same amount crosses from each side. When we consider Figure 15.4.3(b), it seems plausible that the same amount of F2 was crossing C1 from left to right as from right to left.

Green’s Theorem

There is an important connection between the circulation around a closed region R and the curl of the vector field inside of R, as well as a connection between the flux across the boundary of R and the divergence of the field inside R. These connections are described by Green’s Theorem and the Divergence Theorem, respectively. We’ll explore each in turn.

Green’s Theorem states “the counterclockwise circulation around a closed region R is equal to the sum of the curls over R.”

Theorem 15.4.1      Green’s Theorem

Let R be a closed, bounded region of the plane whose boundary C is composed of finitely many smooth curves, let r(t) be a counterclockwise parameterization of C, and let F=M,N where Nx and My are continuous over R. Then

CF𝑑r=RcurlFdA.

We’ll explore Green’s Theorem through an example.

margin:

R

-1

1

1

2

x

y
Figure 15.4.4: The vector field and planar region used in Example 15.4.2.
Example 15.4.2 Confirming Green’s Theorem

Let F=-y,x2+1 and let R be the region of the plane bounded by the triangle with vertices (-1,0), (1,0) and (0,2), shown in Figure 15.4.4. Verify Green’s Theorem; that is, find the circulation of F around the boundary of R and show that is equal to the double integral of curlF over R.

SolutionThe curve C that bounds R is composed of 3 lines. While we need to traverse the boundary of R in a counterclockwise fashion, we may start anywhere we choose. We arbitrarily choose to start at (-1,0), move to (1,0), etc., with each line parameterized by r1(t), r2(t) and r3(t), respectively.

We leave it to the reader to confirm that the following parameterizations of the three lines are accurate:

r1(t) =2t-1,0, for 0t1,with r1(t)=2,0,
r2(t) =1-t,2t, for 0t1,with r2(t)=-1,2,and
r3(t) =-t,2-2t, for 0t1,with r3(t)=-1,-2.

The circulation around C is found by summing the flow along each of the sides of the triangle. We again leave it to the reader to confirm the following computations:

C1F𝑑r1 =010,(2t-1)2+12,0𝑑t=0,
C2F𝑑r2 =01-2t,(1-t)2+1-1,2𝑑t=11/3,and
C3F𝑑r3 =012t-2,t2+1-1,-2𝑑t=-5/3.

The circulation is the sum of the flows: 2.

We confirm Green’s Theorem by computing RcurlFdA. We find curlF=2x+1. The region R is bounded by the lines y=2x+2, y=-2x+2 and y=0. Integrating with the order dxdy is most straightforward, leading to

02y/2-11-y/2(2x+1)𝑑x𝑑y=02(2-y)𝑑y=2,

which matches our previous measurement of circulation.

margin:

R

-2

2

2

2

x

y
Figure 15.4.5: The vector field and planar region used in Example 15.4.3.
Example 15.4.3 Using Green’s Theorem

Let F=sinx,cosy and let R be the region enclosed by the curve C parameterized by r(t)=2cost+110cos(10t),2sint+110sin(10t) on 0t2π, as shown in Figure 15.4.5. Find the circulation around C.

SolutionComputing the circulation directly using the line integral looks difficult, as the integrand will include terms like “sin(2cost+110cos(10t)).”

Green’s Theorem states that CF𝑑r=RcurlFdA; since curlF=0 in this example, the double integral is simply 0 and hence the circulation is 0.

Since curlF=0, we can conclude that the circulation is 0 in two ways. One method is to employ Green’s Theorem as done above. The second way is to recognize that F is a conservative field, hence there is a function z=f(x,y) wherein F=f. Let A be any point on the curve C; since C is closed, we can say that C “begins” and “ends” at A. By the Fundamental Theorem of Line Integrals, CF𝑑r=f(A)-f(A)=0.

One can use Green’s Theorem to find the area of an enclosed region by integrating along its boundary. Let C be a closed curve, enclosing the region R, parameterized by r(t)=f(t),g(t). We know the area of R is computed by the double integral R𝑑A, where the integrand is 1. By creating a field F where curlF=1, we can employ Green’s Theorem to compute the area of R as CF𝑑r.

One is free to choose any field F to use as long as curlF=1. Common choices are F=0,x, F=-y,0 and F=-y/2,x/2. We demonstrate this below.

Example 15.4.4 Using Green’s Theorem to find area

Let C be the closed curve parameterized by r(t)=t-t3,t2 on -1t1, enclosing the region R, as shown in Figure 15.4.6. Find the area of R.margin:

R

-1

1

1

x

y
Figure 15.4.6: The region R, whose area is found in Example 15.4.4.

SolutionWe can choose any field F, as long as curlF=1. We choose F=-y,0. We also confirm (left to the reader) that r(t) traverses the region R in a counterclockwise fashion. Thus

Area of R =R𝑑A
=CF𝑑r
=-11-t2,01-3t2,2t𝑑t
=-11(-t2)(1-3t2)𝑑t
=815.

The Divergence Theorem

Green’s Theorem makes a connection between the circulation around a closed region R and the sum of the curls over R. The Divergence Theorem makes a somewhat “opposite” connection: the total flux across the boundary of R is equal to the sum of the divergences over R.

Theorem 15.4.2      The Divergence Theorem (in the plane)

Let R be a closed, bounded region of the plane whose boundary C is composed of finitely many smooth curves, let r(t) be a counterclockwise parameterization of C, and let F=M,N where Mx and Ny are continuous over R. Then

CFn𝑑s=RdivFdA.
margin:

R

-2

2

-2

2

x

y
Figure 15.4.7: The region R used in Example 15.4.5.
Example 15.4.5 Confirming the Divergence Theorem

Let F=x-y,x+y, let C be the circle of radius 2 centered at the origin and define R to be the interior of that circle, as shown in Figure 15.4.7. Verify the Divergence Theorem; that is, find the flux across C and show it is equal to the double integral of divF over R.

SolutionWe parameterize the circle in the usual way, with
r(t)=2cost,2sint, 0t2π. The flux across C is

CFn𝑑s =C(Mg(t)-Nf(t))𝑑t
=02π((2cost-2sint)(2cost)-(2cost+2sint)(-2sint))𝑑t
=02π4𝑑t=8π.

We compute the divergence of F as divF=Mx+Ny=2. Since the divergence is constant, we can compute the following double integral easily:

RdivFdA=R2𝑑A=2R𝑑A=2(area of R)=8π,

which matches our previous result.

Example 15.4.6 Flux when divF=0

Let F be any field where divF=0, and let C1 and C2 be any two nonintersecting paths, except that each begin at point A and end at point B (see Figure 15.4.8). Show why the flux across C1 and C2 is the same.

Solutionmargin:

C1

B

A

C2
Figure 15.4.8: As used in Example 15.4.6, the vector field has a divergence of 0 and the two paths only intersect at their initial and terminal points.
By referencing Figure 15.4.8, we see we can make a closed path C that combines C1 with C2, where C2 is traversed with its opposite orientation. We label the enclosed region R. Since divF=0, the Divergence Theorem states that

CFn𝑑s=RdivFdA=R0𝑑A=0.

Using the properties and notation given in Theorem 15.3.1, consider:

0 =CFn𝑑s
=C1Fn𝑑s+C2*Fn𝑑s
(where C2* is the path C2 traversed with opposite orientation)
=C1Fn𝑑s-C2Fn𝑑s.
C2Fn𝑑s =C1Fn𝑑s.

Thus the flux across each path is equal.

In this section, we have investigated flow and flux, quantities that measure interactions between a vector field and a planar curve. We can also measure flow along spatial curves, though as mentioned before, it does not make sense to measure flux across spatial curves.

It does, however, make sense to measure the amount of a vector field that passes across a surface in space — i.e, the flux across a surface. We will study this, though in the next section we first learn about a more powerful way to describe surfaces than using functions of the form z=f(x,y).

Exercises 15.4

 

Terms and Concepts

  1. 1.

    Let F be a vector field and let C be a curve. Flow is a measure of the amount of F going                           C; flux is a measure of the amount of F going                           C.

  2. 2.

    What is circulation?

  3. 3.

    Green’s Theorem states, informally, that the circulation around a closed curve that bounds a region R is equal to the sum of                           across R.

  4. 4.

    The Divergence Theorem states, informally, that the outward flux across a closed curve that bounds a region R is equal to the sum of                           across R.

  5. 5.

    Let F be a vector field and let C1 and C2 be any nonintersecting paths except that each starts at point A and ends at point B. If             =0, then C1FT𝑑s=C2FT𝑑s.

  6. 6.

    Let F be a vector field and let C1 and C2 be any nonintersecting paths except that each starts at point A and ends at point B. If             =0, then C1Fn𝑑s=C2Fn𝑑s.

Problems

In Exercises 7–12, a vector field F and a curve C are given. Evaluate CFn𝑑s, the flux of F over C.

  1. 7.

    F=x+y,x-y; C is the curve with initial and terminal points (3,-2) and (3,2), respectively, parametrized by r(t)=3t2,2t on -1t1.

  2. 8.

    F=x+y,x-y; C is the curve with initial and terminal points (3,-2) and (3,2), respectively, parametrized by r(t)=3,t on -2t2.

  3. 9.

    F=x2,y+1; C is line segment from (0,0) to (2,4).

  4. 10.

    F=x2,y+1; C is the portion of the parabola y=x2 from (0,0) to (2,4).

  5. 11.

    F=y,0; C is the line segment from (0,0) to (0,1).

  6. 12.

    F=y,0; C is the line segment from (0,0) to (1,1).

In Exercises 13–26, a vector field F and a closed curve C, enclosing a region R, are given. Verify Green’s Theorem by evaluating CF𝑑r and RcurlFdA, showing they are equal.

  1. 13.

    F=x-y,x+y; C is the closed curve composed of the parabola y=x2 on 0x2 followed by the line segment from (2,4) to (0,0).

  2. 14.

    F=-y,x; C is the unit circle.

  3. 15.

    F=0,x2; C the triangle with corners at (0,0), (2,0) and (1,1).

  4. 16.

    F=x+y,2x; C the curve that starts at (0,1), follows the parabola y=(x-1)2 to (3,4), then follows a line back to (0,1).

  5. 17.

    F=x2-y2,2xy; C is the boundary of R={(x,y):0x1,2x2y2x}

  6. 18.

    F=x2y,2xy; C is the boundary of R={(x,y):0x1,x2yx}

  7. 19.

    F=2y,-3x; C is the circle x2+y2=1

  8. 20.

    F=(ex2+y2),(ey2+x2); C is the boundary of the triangle with vertices (0,0), (4,0) and (0,4)

  9. 21.

    F=xy,y; C is the square with vertices (0,0), (1,0), (1,1), (0,1).

  10. 22.

    F=y2,x2; C is the triangle with vertices (0,0), (1,0), (1,1).

  11. 23.

    F=ysinx,xy; C is the boundary of the region bordered by cosx and sinx for π4x5π4.

  12. 24.

    F=y3,siny-x3; where C is the boundary of the annular region between x2+y2=1/4 and x2+y2=1, traversed so that the region is always on the left.

  13. 25.

    F=exsiny,y3+excosy; C is the boundary of the rectangle with vertices (1,-1), (1,1), (-1,1) and (-1,-1), traversed counterclockwise.

  14. 26.

    F=-yx2+y2,xx2+y2; C is the boundary of the annulus R={(x,y):1/4x2+y21} traversed so that R is always on the left.

  15. 27.

    For F=y2,x2, evaluate Cf𝑑r where C is the boundary of the half disk of radius 1 centered at the origin, in the first and second quadrants, traversed clockwise.

  16. 28.

    Repeat the previous problem, but use the unit half-disk (centered at the origin) in the fourth and first quadrants.

In Exercises 29–32, a closed curve C enclosing a region R is given. Find the area of R by computing CF𝑑r for an appropriate choice of vector field F.

  1. 29.

    C is the ellipse parametrized by r(t)=4cost,3sint on 0t2π.

  2. 30.

    C is the curve parametrized by r(t)=cost,sin(2t) on -π/2tπ/2.

  3. 31.

    C is the curve parametrized by r(t)=3t2-2t-t3,2(t-1)2 on 0t2.

  4. 32.

    C is the curve parametrized by r(t)=2cost+110cos(10t),2sint+110sin(10t) on 0t2π.

In Exercises 33–36, a vector field F and a closed curve C, enclosing a region R, are given. Verify the Divergence Theorem by evaluating CFn𝑑s and RdivFdA, showing they are equal.

  1. 33.

    F=x-y,x+y; C is the closed curve composed of the parabola y=x2 on 0x2 followed by the line segment from (2,4) to (0,0).

  2. 34.

    F=-y,x; C is the unit circle.

  3. 35.

    F=0,y2; C the triangle with corners at (0,0), (2,0) and (1,1).

  4. 36.

    F=x2/2,y2/2; C the curve that starts at (0,1), follows the parabola y=(x-1)2 to (3,4), then follows a line back to (0,1).

  5. 37.

    Show that for any constants a, b and any closed simple curve C, Ca𝑑x+bdy=0.

Modern Campus CMS