DOC PREVIEW
CALTECH APH 162 - Biological Physics Laboratory

This preview shows page 1-2 out of 6 pages.

Save
View full document
View full document
Premium Document
Do you want full access? Go Premium and unlock all 6 pages.
Access to all documents
Download any document
Ad free experience
View full document
Premium Document
Do you want full access? Go Premium and unlock all 6 pages.
Access to all documents
Download any document
Ad free experience
Premium Document
Do you want full access? Go Premium and unlock all 6 pages.
Access to all documents
Download any document
Ad free experience

Unformatted text preview:

Diffusion of Solid Particles Confined in a Viscous FluidAPh 162 – Biological Physics Laboratory Diffusion of Solid Particles Confined in a Viscous Fluid1 I realize that some of you have little or no experience with partial differential equations (PDEs). However we think it is important to understand the nature of the diffusion process, especially as it relates to biology, to this end I would like to go through the theory behind the experiment you are about to do. This reading is certainly of the crash-course variety, so feel free to ask Rob, Hernan, or me any questions. Solution of the 2D Diffusion Equation: The 2D diffusion equation allows us to talk about the statistical movements of randomly moving particles in two dimensions. By random, we mean that we cannot correlate the movement at one moment to movement at the next, or in other words there is no deterministic/predictive power over the exact motion of the particle. The movement of each individual particle moving in a Brownian (diffuse) way does not follow the diffusion equation. However, many identical particles each obeying the same boundary and initial conditions share some statistical properties dealing with their spatial and temporal evolution. It is those statistical properties that the diffusion equation captures. The function P(x,y,t) gives the probability of finding a perfectly average particle in the small vicinity of the point (x,y) at time t (we call this a probability distribution). Such properties are for instance the average position of a statistical particle as a function of time. The evolution of some systems does follow the diffusion equation outright: when you put a drop of dye in a beaker, there are millions of dye molecules each exhibiting Brownian motion, but as a group they exhibit the smooth, well-behaved statistical features of the diffusion equation. Consider the two-dimensional diffusion equation in Cartesian coordinates: 012=∂∂−∇tPDP → 012222=∂∂−∂∂+∂∂tPDyPxP The diffusion equation can be derived from the probabilistic nature of Brownian motion described as random walks (speak with me if you really want to see the derivation). The constant D is the diffusion coefficient whose nature we will explore in a moment, but for now we are solving a math problem. Notice we have no powers or complicated functions of any derivative of the function P. Since the derivatives of P only appear to first power, we call this a linear PDE. Also notice that the equation does not mix x’s, y’s, or time within a single term, that is there are no multiplicative, divisional, or functional terms in the equation between the independent variables, we call this feature separability of the PDE. These two properties, linearity and separability, make the problem vastly more tractable than many other PDE’s encountered in physics. I am going to attempt to walk you through a challenging but very powerful method for solving PDE’s using the Fourier Transform (FT). Many of you may know that the FT is used in signal analysis and manipulation, but it was first used by Fourier to solve this problem. Consider the following integral relations that define the 2D FT in Cartesian coordinates. We will call the function Pˆthe FT of our original function P: ∫∫∫∫∞∞−∞∞−⋅+⋅∞∞−∞∞−⋅+⋅−==yxyxykxkiykxkiyxdkdktkkPetyxPdxdytyxPetkkPyxyx),,(ˆ),,(),,(),,(ˆ)(2)(2ππ Notice the symmetry in going forward and backward in the transform. This is because switching between the normal form of the problem and what we call Fourier Space, where the problem exists after the FT, are physically identical. We will require one more fact familiar from calculus, namely integration by parts:∫∫−=bababavduuvudv Using these facts, let’s examine the spatial derivatives of the diffusion equation, where we consider the second derivative to be the function of interest. We can integrate these second derivatives by parts, identifying u and v as: )(2 ykxkiyxeu⋅+⋅−=π and xPv∂∂= dxexduykxkiyx)(2 ⋅+⋅−∂∂=π and dxxPdv22∂∂= then integration by parts tells us: ∫∫∫∫∞∞−∞∞−⋅+⋅−∞∞−⋅+⋅−∞∞−∞∞−⋅+⋅−∂∂⋅∂∂−∂∂=∂∂dxdyxPexxPedxdyxPeykxkiyxykxkiykxkiyxyxyx)(2),()(222)(2πππ ∫∫∞∞−∞∞−⋅+⋅−∞∞−⋅+⋅−∂∂⋅+∂∂= dxdyxPeikxPeykxkixyxykxkiyxyx)(2),()(22πππ Notice in the first line we effectively ‘transferred’ one derivative off P and put it on the exponential of the FT, but since the exponential is an explicit function we can just perform the derivative, giving us the constant on the right most integral of the second line. It should be physically intuitive that at an infinite distance from where we originally introduced our statistically perfect particle there should be zero chance for all time of it appearing. This leads us to the boundary condition that the function P and all of its derivatives are zero at ∞== yx , eliminating the left term of the second line, and allowing us to write: ∫∫∫∫∞∞−∞∞−⋅+⋅−∞∞−∞∞−⋅+⋅−∂∂⋅=∂∂dxdyxPeikdxdyxPeykxkixykxkiyxyx)(222)(22πππ Looking at the right hand side, it looks like we are taking the FT of a first derivative of P. We pursue this integral in the same manner: ∫∫∫∫∞∞−∞∞−⋅+⋅−∞∞−⋅+⋅−∞∞−∞∞−⋅+⋅−⋅+⋅=∂∂⋅ dxdyePikePdxdyxPeykxkixyxykxkiykxkiyxyxyx)(2),()(2)(22ππππ Recall, according to our boundary condition, P = 0 at ∞== yx , and we have: ∫∫∫∫∞∞−∞∞−⋅+⋅−∞∞−∞∞−⋅+⋅−⋅=∂∂⋅ dxdyePikdxdyxPeykxkixykxkiyxyx)(2)(22πππFinally we can say that: ()()PikdxdyxPedxdyePikdxdyxPexykxkiykxkixykxkiyxyxyxˆ22222)(2)(2222)(2πππππ=∂∂⋅=∂∂∫∫∫∫∫∫∞∞−∞∞−⋅+⋅−∞∞−∞∞−⋅+⋅−∞∞−∞∞−⋅+⋅− In general, if you continued the above procedure, you would find that: ()PikPxnxnnˆ2π=∂∂ Now keep in mind, we took the spatial FT (i.e. dealing with x and y), so the derivative in time does not change under a FT. This means we can write the 2D diffusion equation after our FT as: ()0ˆ1)(ˆ2222=∂∂++⋅tPDkkPyxπ → ()0ˆ)(2ˆ222=⋅++∂∂PkkDtPyxπ Amazing! We’ve completely eliminated our spatial dependence; this remaining equation is a simple first order ODE in time, with the


View Full Document
Download Biological Physics Laboratory
Our administrator received your request to download this document. We will send you the file to your email shortly.
Loading Unlocking...
Login

Join to view Biological Physics Laboratory and access 3M+ class-specific study document.

or
We will never post anything without your permission.
Don't have an account?
Sign Up

Join to view Biological Physics Laboratory 2 2 and access 3M+ class-specific study document.

or

By creating an account you agree to our Privacy Policy and Terms Of Use

Already a member?