The following content is
provided under a Creative Commons license.
Your support will help MIT OpenCourseWare continue to offer
high quality educational resources for free.
To make a donation or to view additional materials from
hundreds of MIT courses, visit MIT OpenCourseWare at
ocw.mit.edu. Today I am going to tell you
about flux of a vector field for a curve.
In case you have seen flux in physics,
probably you have seen flux in space,
and we are going to come to that in a couple of weeks,
but for now we are still doing everything in the plane.
So bear with me if you have seen a more complicated version
of flux. We are going to do the easy one
first. What is flux?
Well, flux is actually another kind of line integral.
Let's say that I have a plane curve and a vector field in the
plane. Then the flux of F across a
curve C is, by definition, a line integral,
but I will use notation F dot n ds.
I have to explain to you what it means, but let me first box
that because that is the important formula to remember.
That is the definition. What does that mean?
First, mostly I have to tell you what this little n is.
The notation suggests it is a normal vector,
so what does that mean? I have a curve in the plane and
I have a vector field. Let's see.
The vector field will be yellow today.
And I will want to integrate along the curve the dot product
of F with the normal vector to the curve, a unit normal vector
to the curve. That means a vector that is at
every point of the curve perpendicular to the curve and
has length one. N everywhere will be the unit
normal vector to the curve C pointing 90 degrees clockwise
from T. What does that mean?
That means I have two normal vectors, one that is pointing
this way, one that is pointing that way.
I have to choose a convention. And the convention is that the
normal vector that I take goes to the right of the curve as I
am traveling along the curve. You mentioned that you were
walking along this curve, then you look to your right,
that is that direction. What we will do is just,
at every point along the curve, the dot product between the
vector field and the normal vector.
And we will sum that along the various pieces of the curve.
What this notation means is that if we actually break C into
small pieces of length delta(s) then the flux will be the limit,
as the pieces become smaller and smaller,
of the sum of F dot n delta S. I take each small piece of my
curve, I do the dot product between F and n and I multiply
by the length of a piece. And then I add these together.
That is what the line integral means.
Of course that is, again, not how I will compute
it. Just to compare this with work,
conceptually it is similar to the line integral we did for
work except the line integral for work -- Work is the line
integral of F dot dr, which is also the line integral
of F dot T ds. That is how we reformulated it.
That means we take our curve and we figure out at each point
how big the tangent component -- I guess I should probably take
the same vector field as before. Let's see.
My field was pointing more like that way.
What I do at any point is project F to the tangent
direction, I figure out how much F is going along my curve and
then I sum these things together.
I am actually summing -- -- the tangential component of my field
F. Roughly-speaking the work
measures, you know, when I move along my curve,
how much I am going with or against F.
Flux, on the other hand, measures, when I go along the
curve, roughly how much the field is going to across the
curve. Counting positively what goes
to the right, negatively what goes to the
left. Flux is integral F dot n ds,
and that one corresponds to summing the normal component of
a vector field. But apart from that
conceptually it is the same kind of thing.
Just the physical interpretations will be very
different, but for a mathematician these
are two line integrals that you set up and compute in pretty
much the same way. Let's see.
I should probably tell you what it means.
Why do we make this definition? What does it correspond to?
Well, the interpretation for work made a lot of sense when F
was representing a force. The line integral was actually
the work done by the force. The interpretation for flux
makes more sense if you think of F as a velocity field.
What is the interpretation? Let's say that for F is a
velocity field. That means I am thinking of
some fluid that is moving, maybe water or something,
and it is moving at a certain speed.
And my vector field represents how things are moving at every
point of the plane. I claim that flux measures how
much fluid passes through -- -- the curve C per unit time.
If you imagine that maybe you have a river and you are somehow
building a damn here, a damn with holes in it so that
the water still passes through, then this measures how much
water passes through your membrane per unit time.
Let's try to figure out why this is true.
Why does this make sense? Let's look at what happens on a
small portion of our curve C. I am zooming in on my curve C.
I guess I need to zoom further. That is a little piece of my
curve, of length delta S, and there is a fluid flow.
On my picture things are flowing to the right.
Here I am drawing a constant vector field because if you zoom
in enough then your vectors will pretty much be the same
everywhere. If you enlarge the picture
enough then things will be pretty much a uniform flow.
Now, how much stuff goes through this little piece of
curve per unit time? Well, what happens over time is
the fluid is moving while my curve is staying the same place
so it corresponds to something like this.
I claim that what goes through C in unit time is actually going
to be a parallelogram. Here is a better picture.
I claim that what will be going through C is this shaded
parallelogram to the left of C. Let's see.
If I move for unit time it works.
That is the stuff that goes through my curve,
for a small portion of curve in unit time.
And, of course, I would need to add all of
these together to get the entire curve.
Let's try to understand how big this parallelogram is.
To know how big this parallelogram is I would like to
use base times height or something like that.
And maybe I want to actually flip my picture so that the base
and the height make more sense to me.
Let me actually turn it this way.
And, in case you have trouble reading the rotated picture,
let me redo it on the board. What passes through a portion
of C in unit time is the contents of a parallelogram
whose base is on C. So it has length delta s.
That is a piece of C. And the other side is going to
be given by my velocity vector F.
And to find the height of this thing, I need to know what
actually the normal component of this vector is.
If I call n the unit normal vector to the curve then the
area is base times height. The base is delta S and the
height is the normal component of F, so it is F dot n.
And so you see that when you sum these things together you
get, what I said, flux.
Now, if you are worried about the fact that actually -- If
your unit time is too long then of course things might start
changing as it flows. You have to take the time unit
and the length unit that are sufficiently small so that
really this approximation where C is a straight line and where
flow is at constant speed are valid.
You want to take maybe a segment here that is a few
micrometers. And the time unit might be a
few nanoseconds or whatever, and then it is a good
approximation. What I mean by per unit time
is, well, actually, that works, but you want to
think of a really, really small time.
And then the amount of matter that passes in that really,
really small time is the flux times the amount of time.
Let's be a tiny bit more careful.
And what I am saying is the amount of stuff that passes
through C depends actually on whether n is going this way or
the opposite way. Actually,
what is implicit in this explanation is that I am
counting positively all the stuff that flows across C in the
direction of n and negatively what flows in the opposite
direction. What flows to the right of C,
well, across C from left to right is counted positively.
While what flows right to left is counted negatively.
So, in fact, it is the net flow through C
per unit time. Any questions about the
definition or the interpretation or things like that?
Yes? Well, you can have both not in
the same small segment. But it could be that,
well, imagine that my vector field accidentally goes in the
opposite direction then this part of the curve,
while things are flowing to the left,
contributes negatively to flux. And here maybe the field is
tangent so the normal component becomes zero.
And then it becomes positive and this part of the curve
contributes positively. For example,
if you imagine that you have a round tank in which the fluid is
rotating and you put your dam just on a diameter across then
things are going one way on one side,
the other way on the other side,
and actually it just evens out. We don't have complete
information. It is just the total net flux.
OK. If there are no other questions
then I guess we will need to figure out how to compute this
guy and how to actually do this line integral.
Well, let's start with a couple of easy examples.
Let's say that C is a circle of radius (a) centered at the
origin going counterclockwise. And let's say that our vector
field is xi yj. What does that look like?
Remember, xi plus yj is a vector field that is pointing
radially away from the origin. Because at every point it is
equal to the vector from the origin to that point.
Now, if we have a circle and let's say we are going
counterclockwise. Actually, I have a nicer
picture. Let me do it here.
That is my curve and my vector field.
And the normal vector, see, when you go counterclockwise in
a closed curve, this convention that a normal
vector points to the right of curve makes it point out.
The usual convention, when you take flux for a closed
curve, is that you are counting the flux going out of the region
enclosed by the curve. And, of course,
if you went clockwise it would be the other way around.
You choose to do it the way you want, but the most common one is
to count flux going out of the region.
Let's see what happens. Well, if I am anywhere on my
circle, see, the normal vector is sticking straight out of the
circle. That is a property of the
circle that the radial direction is perpendicular to the circle.
Actually, let me complete this picture.
If I take a point on the circle, I have my normal vector
that is pointing straight out so it is parallel to F.
Along C we know that F is parallel to n,
so F dot n will be equal to the magnitude of F times,
well, the magnitude of n, but that is one.
Let me put it anywhere, but that is the unit normal
vector. Now, what is the magnitude of
this vector field if I am at a point x, y?
Well, it is square root of x squared plus y squared,
which is the same as the distance from the origin.
So if this distance, if this radius is a then the
magnitude of F will just be a. In fact, F dot n is constant,
always equal to a. So the line integral will be
pretty easy because all I have to do is the integral of F dot n
ds becomes the integral of a ds. (a) is a constant so I can take
it out. And integral ds is just a
length of C which is 2pi a, so I will get 2pi a squared.
And that is positive, as we expected,
because stuff is flowing out of the circle.
Any questions about that? No.
OK. Just out of curiosity,
let's say that we had taken our other favorite vector field.
Let's say that we had the same C, but now the vector field
. Remember, that one goes
counterclockwise around the origin.
If you remember what we did several times,
well, along the circle that vector field now is tangent to
the circle. If it is tangent to the circle
it doesn't have any normal component.
The normal component is zero. Things are not flowing into the
circle or out of it. They are just flowing along the
circle around and around so the flux will be zero.
F now is tangent to C. F dot n is zero and,
therefore, the flux will be zero.
These are examples where you can compute things
geometrically. And I would say,
generally speaking, with flux, well,
if it is a very complicated field then you cannot.
But, if a field is fairly simple,
you should be able to get some general feeling for whether your
answer should be positive, negative or zero just by
thinking about which way is my flow going.
Is it going across the curve one way or the other way?
Still no questions about these examples?
The next thing we need to know is how we will actually compute
these things because here, yeah, it works pretty well,
but what if you don't have a simple geometric interpretation.
What if I give you a really complicated curve and then you
have trouble finding the normal vector?
It is going to be annoying to set up things this way.
Actually, there is a better way to do it in coordinates.
Just as we do work, when we compute this line
integral, usually we don't do it geometrically like this.
Most of the time we just integrate M dx plus N dy in
coordinates. That is a similar way to do it
because it is, again, a line integral so it
should work the same way. Let's try to figure that out. How do we do the calculation in
coordinates, or I should say using components?
That is the general method of calculation when we don't have
something geometric to do. Remember, when we were doing
things for work we said this vector dr, or if you prefer T
ds, we said just becomes symbolically dx and dy.
When you do the line integral of F dot dr you get line
integral of n dx plus n dy. Now let's think for a second
about how we would express n ds. Well, what is n ds compared to
T ds? Well, M is just T rotated by 90
degrees, so n ds is T ds rotated by 90 degrees.
That might sound a little bit outrageous because these are
really symbolic notations but it works.
I am not going to spend too much time trying to convince you
carefully. But if you go back to where we
wrote this and how we tried to justify this and you work your
way through it, you will see that n ds can be
analyzed the same way. N is T rotated 90 degrees
clockwise. That tells us that n ds is --
How do we rotate a vector by 90 degrees?
Well, we swept the two components and we put a minus
sign. You have dy and dx.
And you have to be careful where to put the minus sign.
Well, if you are doing it clockwise, it is in front of dx. Well, actually,
let me just convince you quickly.
Let's say we have a small piece of C.
If we do T delta S, that is also vector delta r.
That is going to be just the vector that goes along the curve
given by this. Its components will be indeed
the change in x, delta x, and the change in y,
delta y. And now, if I want to get n
delta S, well, I claim now that it is
perfectly valid and rigorous to just rotate that by 90 degrees.
If I want to rotate this by 90 degrees clockwise then the x
component will become the same as the old y component.
And the y component will be minus delta x.
Then you take the limit when the segment becomes shorter and
shorter, and that is how you can justify this.
That is the key to computing things in practice.
It means, actually, you already know how to compute
line integrals for flux. Let me just write it explicitly.
Let's say that our vector field has two components.
And let me just confuse you a little bit and not call them M
and N for this time just to stress the fact that we are
doing a different line integral. Let me call them P and Q for
now. Then the line integral of F dot
n ds will be the line integral of
dot product . That will be the integral of -
Q dx P dy. Well, I am running out of space
here. It is integral along C of
negative Q dx plus P dy. And from that point onwards you
just do it the usual way. Remember, here you have two
variables x and y but you are integrating along a curve.
If you are integrating along a curve x and y are related.
They depend on each other or maybe on some other parameter
like T or theta or whatever. You express everything in terms
of a single variable and then you do a usual single integral.
Any questions about that? I see a lot of confused faces
so maybe I shouldn't have called my component P and Q. If you prefer,
if you are really sentimentally attached to M and N then this
new line integral becomes the integral of - N dx M dy.
If a problem tells you compute flux instead of saying compute
work, the only thing you change is
instead of doing M dx plus N dy you do minus N dx plus M dy.
And I am sorry to say that I don't have any good way of
helping you remember which one of the two gets the minus sign,
so you just have to remember this formula by heart.
That is the only way I know. Well, you can try to go through
this argument again, but it is really best if you
just remember that formula. I am not going to do an example
because we already know how to do line integrals.
Hopefully you will get to see one at least in recitation on
Monday. That is all pretty good.
Let me tell you now what if I have to compute flux along a
closed curve and I don't want to compute it?
Well, remember in the case of work we had Green's theorem.
We saw yesterday Green's theorem.
Let's us replace a line integral along a closed curve by
a double integral. Well, here it is the same.
We have a line integral along a curve.
If it is a closed curve, we should be able to replace it
by a double integral. There is a version of Green's
theorem for flux. And you will see it is not
scarier than the other one. It is perhaps less scarier or
perhaps just as scary or just not as scary,
depending on how you feel about it, but it works pretty much the
same way. What does Green's theorem for
flux say? It says if C is a curve that
encloses a region R counterclockwise and if I have a
vector field that is defined everywhere,
not just on C but also inside, so also on R.
Well, maybe I should give names to the components.
If you will forgive me for a second, I will still use P and Q
for now. You will see why.
It is defined and differentiable in R.
Then I can actually -- -- replace the line integral for
flux by a double integral over R of some function.
And that function is called the divergence of F dA.
This is the divergence of F. And I have to define for you
what this guy is. The divergence of a vector
field with components P and Q is just P sub x Q sub y.
This one is actually easier to remember than curl because you
just take the x component, take its partial with respect
to x, take the y component,
take its partial with respect to y and add them together.
No signs. No switching things around.
This one is pretty straightforward.
The picture again is if I have my curve C going
counterclockwise around a region R and I want to find the flux of
some vector field F that is everywhere in here.
Maybe some parts of C will contribute positively and some
parts will contribute negatively.
Just to reiterate what I said, positively here means,
because we are going counterclockwise,
the normal vector points out of the region.
This guy here is the flux out of R through C.
That is the formula. Any questions about what the
statement says or how to use it concretely?
No. OK.
It is pretty similar to Green's theorem for work.
Actually, I should say -- This is called Green's theorem in
normal form also. Not that the other one is
abnormal, but just that the old one for work was,
you could say, in tangential form.
That just means, well, Green's theorem,
as seen yesterday was for the line integral F dot T ds,
integrating the tangent component of F.
The one today is for integrating the normal component
of F. OK. Let's prove this.
Good news. It is much easier to prove than
the one we did yesterday because we are just going to show that
it is the same thing just using different notations. How do I prove it?
Well, maybe actually it would help if first,
before proving it, I actually rewrite what it
means in components. We said the line integral of F
dot n ds is actually the line integral of - Q dx P dy.
And we want to show that this is equal to the double integral
of P sub x Q sub y dA. This is really one of the
features of Green's theorem. No matter which form it is,
it relates a line integral to a double integral.
Let's just try to see if we can reduce it to the one we had
yesterday. Let me forget what these things
mean physically and just focus on the math.
On the math it is a line integral of something dx plus
something dy. Let's call this guy M and let's
call this guy N. Let M equal negative Q and N
equal P. Then this guy here becomes
integral of M dx plus N dy. And I know from yesterday what
this is equal to, namely using the tangential
form of Green's theorem. Green for work.
This is the double integral of curl of this guy.
That is Nx minus My dA. But now let's think about what
this is in terms of M and N. Well, we said that M is
negative Q so this is negative My.
And we said P is the same as N, so this is Nx.
Just by renaming the components, I go from one form
to the other one. So it is really the same
theorem. That's why it is also called
Green's theorem. But the way we think about it
when we use it is different, because one of them computes
the work done by a force along a closed curve,
the other one computes the flux maybe of a velocity field out of
region. Questions?
Yes? That is correct.
If you are trying to compute a line integral for flux,
wait, where did I put it? A line integral for flux just
becomes this. And once you are here you know
how to compute that kind of thing.
The double integral side does not even have any kind of
renaming to do. You know how to compute a
double integral of a function. This is just a particular kind
of function that you get out of a vector field,
but it is like any function. The way you would evaluate
these double integrals is just the usual way.
Namely, you have a function of x and y, you have a region and
you set up the bounds for the isolated integral.
The way you would evaluate the double integrals is really the
usual way, by slicing the region and
setting up the bounds for iterated integrals in dx,
dy or dydx or maybe rd, rd theta or whatever you want.
In fact, in terms of computing integrals, we just have two sets
of skills. One is setting up and
evaluating double integrals. The other one is setting up and
evaluating line integrals. And whether these line
integrals or double integrals are representing work,
flux, integral of a curve, whatever,
the way that we actually compute them is the same.
Let's do an example. Oh, first. Sorry.
This renaming here, see, that is why actually I call my
components P and Q because the argument would have gotten very
messy if I had told you now I call M ,N and I call N minus M
and so on. But, now that we are through
with this, if you still like M and N
better, you know, what this says -- The
formulation of Green's theorem in this language is just
integral of minus N dx plus M dy is the double integral over R of
Mx plus Ny dA. Now let's do an example.
Let's look at this picture again, the flux of xi plus yj
out of the circle of radius A. We did the calculation directly
using geometry, and it wasn't all that bad.
But let's see what Green's theorem does for us here. Example.
Let's take the same example as last time.
F equals xi yj. C equals circle of radius a
counterclockwise. How do we set up Green's
theorem. Well, let's first figure out
the divergence of F. The divergence of this field,
I take the x component, which is x, and I take its
partial respect to x. And then I do the same with the
y component, and I will get one plus one equals two.
So, the divergence of this field is two.
Now, Green's theorem tells us that the flux out of this region
is going to be the double integral of 2 dA.
What is R now? Well, R is the region enclosed
by C. So if C is the circle,
R is the disk of radius A. Of course, we can compute it,
but we don't have to because double integral of 2dA is just
twice the double integral of dA so it is twice the area of R.
And we know the area of a circle of radius A.
That is piA2. So, it is 2piA2.
That is the same answer that we got directly,
which is good news. Now we can even do better.
Let's say that my circle is not at the origin.
Let's say that it is out here. Well, then it becomes harder to
calculate the flux directly. And it is harder even to guess
exactly what will happen because on this side here the vector
field will go into the region so the contribution to flux will be
negative here. Here it will be positive
because it is going out of the region.
There are positive and negative terms.
Well, it looks like positive should win because here the
vector field is much larger than over there.
But, short of computing it, we won't actually know what it
is. If you want to do it by direct
calculation then you have to parametize this circle and
figure out what the line integral will be.
But if you use Green's theorem, well, we never used the fact
that it is the circle of radius A at the origin.
It is true actually for any closed curve that the flux out
of it is going to be twice the area of the region inside.
It still will be 2piA2 even if my circle is anywhere else in
the plane. If I had asked you a trick
question where do you want to place this circle so that that
the flux is the largest? Well, the answer is it doesn't
matter. Now, let's just finish quickly
by answering a question that some of you,
I am sure, must have, which is what does divergence
mean and what does it measure? I mean, we said for curl,
curl measures how much things are rotating somehow.
What does divergence mean? Well, the answer is divergence
measures how much things are diverging.
Let's be more explicit. Interpretation of divergence.
You can think of it, you know, what do I want to say
first? If you take a vector field that
is a constant vector field where everything just translates then
there is no divergence involved because the derivatives will be
zero. If you take the guy that
rotates things around you will also compute and find zero for
divergence. This is not sensitive to
translation motions where everything moves together or to
rotation motions, but instead it is sensitive to
explaining motions. A possible answer is that it
measures how much the flow is expanding areas.
If you imagine this flow that we have here on the picture,
things are moving away from the origin and they fill out the
plane. If we mention this fluid
flowing out there, it is occupying more and more
space. And so that is what it means to
have positive divergence. If you took the opposite vector
field that contracts everything to the origin that will have
negative divergence. That is a good way to think
about it if you are thinking of a gas maybe that can expand to
fill out more volume. If you thinking of water,
well, water doesn't really shrink or expand.
The fact that it is taking more and more space actually means
that there is more and more water.
The other way to think about it is divergence is the source
rate, it is the amount of fluid that
is being inserted into the system,
that is being pumped into the system per unit time per unit
area. What div F equals two here
means is that here you actually have matter being created or
being pumped into the system so that you have more and more
water filling more and more space as it flows.
But, actually, divergence is not two just at
the origin. It is two everywhere.
So, in fact, to have this you need to have a
system of pumps that actually is in something water absolutely
everywhere uniformly. That is the only way to do this.
I mean if you imagine that you just have one spring at the
origin then, sure, water will flow out,
but as you go further and further away it will do so more
and more slowly. Well, here it is flowing away
faster and faster. And that means everywhere you
are still pumping more water into it.
So, that is what divergence measures.