Diffusion equation

cosmos 6th March 2017 at 10:55am
Diffusion

Pt=D2P\frac{\partial P}{\partial t} = D \nabla^2 P,

where D is the diffusion coefficient, which when derived from a random walk is

D=Δx22dtD = \frac{\langle \Delta x^2 \rangle}{2dt}

where dd is the dimension of space. The 22 comes from the fact that walker can jump in any of two directions, per dimension. Δx2\langle \Delta x^2 \rangle and tt represent the expected distance squared, and the time step in the random walk, respectively. See for example these notes, for derivation.

See also a simple kinetic derivation of diffusion coefficient (in the context of solid state diffusion), see page 7 Also see Alex's notes on kinetic theory.

https://en.wikipedia.org/wiki/Diffusion

Solutions of the diffusion equation

Can show uniqueness.

Solutions to diffusion equation, using Fourier transformation, and using Green functions. Can also derive from Fokker-Planck equation.

Solutions to diffusion equation for free, absorbing, and reflecting boundary conditions.

Solution by Separation of variables

Assume P=X(x)T(t)P=X(x)T(t), and get an equation

TT=XX\frac{T'}{T}=\frac{X''}{X}

which implies both sides are equal to the same constant, which can be shown to need to be negative, so: TT=α2\frac{T'}{T}=-\alpha^2 for TT and XX=α2\frac{X''}{X}=-\alpha^2

This has solutions

T(t)=Ceα2DtT(t)=Ce^{-\alpha^2 Dt} X(x)=Asinαx+BcosαxX(x)=A\sin{\alpha x} + B\cos{\alpha x}

and can then apply boundary conditions.