• If you are citizen of an European Union member nation, you may not use this service unless you are at least 16 years old.

  • You already know Dokkio is an AI-powered assistant to organize & manage your digital files & messages. Very soon, Dokkio will support Outlook as well as One Drive. Check it out today!

View
 

Fick's law

Page history last edited by Joe Redish 12 years, 3 months ago

Class content > The micro to macro connectionThermodynamics and Statistical PhysicsDiffusion and random walks

 

 

Prerequisites:

 

In our discussion of diffusion, we modeled what was happening in terms of a simple mathematical model of random walkers taking a step back and forth randomly with equal probabilities. The mathematics of probabilities, finding the number of different ways that lead to the same final displacement (sort of like figuring out the odds of throwing a 7 or 11 with a pair of dice) led to a normal (Gaussian) distribution with a height that fell like one over the square root of the number of steps and a width that grew like the square root of the number of steps. We interpreted the number of steps as time.  So our analysis said that if we released a blob of chemical at a point in space that it would spread out; the value at the starting point would decrease by a factor of the square root of time, and the width of the distribution would increase like the square root of time.

 

We might or might not find this mathematical model convincing or appropriate. What we really want to know is:

 

If we have a concentration of molecules in one region, how does it spread to a region with a lower concentration?

 

And for biological mechanisms we are interested in how this works at cellular scales -- scales of microns.  Since each cubic micron typically contains trillions of molecules, we can use kinetic theory type arguments to figure out how concentrations change by diffusion.

 

A molecular model of diffusion

We'll make a more sophisticated mathematical model of molecular motions, but one simple enough to get some useful results.

 

  • We have a uniform fluid (colored black) containing a bunch of interesting molecules (colored red) that has a concentration that varies as a function of position. For simplicity, let's  assume that it varies in just one dimension, say the "x" direction.
  • Rather than working with a continuously varying concentration, n(x), let's break the x-direction up into tiny steps of width Δx and treat the concentration as if it were a constant in each bin.  (This is like the way we think about an integral.)

 

We are interested in how the concentration changes as a result of the random molecular motions. Imagine an area, A,  perpendicular to the x-axis. It lies in the y-z plane that is perpendicular to the screen at the point we have marked x = 0 and shown by a slightly thicker blue line.


 

The molecules we are interested in will cross A from both directions. What is the net number of particles will pass through the area A in a given amount of time?  That number is called the particle flux, and given the symbol J. So

 

net number of molecules crossing A in a time Δt = JAΔt.

 

(When we did this sort of thing in our study of fluid flow, we counted the volume crossing a unit area per unit time, Q, but we could get away with that because the density of our fluid was assumed to be constant. We have to be more careful here so we introduce a new symbol.)

 

Now let's simplify our model.

  • We only care about the red molecules, so let's ignore all the (black) molecules of the fluid.
  • Molecules are moving in the y and z directions, but that doesn't effect whether they cross our area, A, so let's ignore motion in those directions.
  • Molecules are also moving at different speeds, but let's simplify that by using an average speed for all our molecules, <vx>.
  • Let's choose our bin width, Δx, to equal the average distance our molecule will travel between collisions: its mean free path, λ.

 

The resulting model looks like the figure at the right. More red molecules will hit our area A from the left than from the right so there will be a net flux from negative to positive x.

We should consider a time interval that will take all the molecules in our bin to the left of A that are going to the right through A.  If they don't undergo a collision, this will make Δx equal to 

 

λ = <vxt

Notice two things:

  • The number of red molecules in each x-bin decreases as we move to the right.
  • In each x-bin there are the SAME number of molecules going in each direction!

This means:

The net flow of red molecules from left to right is a result of the change in concentration, not because the molecules are moving more preferentially in one direction than in another.

 

Finding the flux: Fick's law

Now we can work out what J must be. The number of red molecules that pass through A from the left in a time Δ is the volume of the cylinder based on A times the number density of red molecules in that volume on the left.

 

1/2  x density of red molecules in the bin to the left x volume of cylinder = 1/2 nL (Aλ )

 

where nL is the density of red molecules in the bin on the left of x = 0. The half comes from that only half of the red molecules in the bin, on the average, are moving to the right.

 

The number of red molecules that pass through A from the right in a time Δ is the volume of the cylinder based on A times the number density of red molecules in that volume on the right.

 

1/2  x density of red molecules in the bin to the right x volume of cylinder = 1/2 nR (Aλ).

 

where nR is the density of red molecules in the bin on the right of x = 0. The half comes from that only half of the red molecules in the bin, on the average, are moving to the left.

 

The resulting net flow to the right (assuming the density is greater on the left) is our J (number per unit are per unit time) times the area times the time:

 

Dividing by A and Δwe get 

Since in this model, the width of the bin divided by our time interval is the average x velocity, we get that J is proportional to the concentration gradient times a combination of constants: half the average velocity times the mean free path.

The above equation is called Fick's Law, and is simply the statement of how the changing concentration results in a flow due to diffusion.

 

Making sense of the result

In order to "find the dog" in the equation that is Fick's law, note what this analysis tells us:

  1. The flow of particles is driven by the concentration gradient -- how it changes in space.
  2. The minus sign tells us that the flow is opposite the derivative -- in the direction of decreasing concentration. So the flow is from high concentration to low.

  1. The flow results from the difference in concentration, not from the direction of the molecule's motion. In any small region, whatever the concentration, the molecules are going every which way.
  2. The diffusion coefficient, D, that tells how fast the flow happens, is proportional to the product of the average speed of the molecules times their mean free path -- the average distance they travel before their velocities get reoriented in a new direction.

 

If we have a concentration that varies in 3D, not just in x, our equation requires the gradient -- the vector derivative.  Then the result looks like this:

The "downhill" of the concentration provides the direction of the diffusive flow.

 

Ben Geller and Joe Redish 12/3/11

 

 

 

 

 

 

 

 

 

 

 

 

Comments (0)

You don't have permission to comment on this page.