Register to reply

Diffusion equation with an external force

by tehrv
Tags: diffusion, equation, external, force
Share this thread:
tehrv
#1
May26-13, 11:16 PM
P: 11
I thought I would do a little physics for fun, since it has been over 20 years since I last did any. I picked up Sethna's Stat Mech book. He gives a derivation of the diffusion equation that goes as follows:

A particle makes random steps, ie x(t+Δt)=x(t) + l(t). The steps l(t) are given according to a distritubion χ whose moments are:

[itex]
\int \chi (z) dz = 1
[/itex]


[itex]
\int z\chi (z) dz = 0
[/itex]


[itex]
\int z^2\chi (z) dz = a^2
[/itex]

For a particle to go from x' at t to x at t+Δt the step, l(t), must be x-x'. This happens with probability χ(x-x') times the probability density ρ(x', t) that it started at x'.

[itex]
\rho (x, t+\Delta t) = \int \rho (x', t)\chi (x-x') dx'
[/itex]

[itex]
\rho (x, t+\Delta t) = \int \rho (x-z, t)\chi (z) dz
[/itex]

Assume ρ is broad and do a Taylor series.

[itex]
\rho (x, t+\Delta t) ≈ \int \left[\rho (x, t)-z\frac{\partial \rho}{\partial x}+z^2/2\frac{\partial^2 \rho}{\partial x^2} \right]\chi (z) dz
[/itex]

Using the moments, we get:

[itex]
\rho (x, t+\Delta t) ≈ \rho (x, t)+1/2\frac{\partial^2 \rho}{\partial x^2}a^2
[/itex]

Now assume ρ is slow and so changes little during this time step:

[itex]
\frac{\partial \rho}{\partial t}\Delta t ≈ \rho (x, t+\Delta t) - \rho (x, t)
[/itex]

So we get:

[itex]
\frac{\partial \rho}{\partial t} = D\frac{\partial^2 \rho}{\partial x^2}
[/itex]

where

[itex]
D = \frac{a^2}{2\Delta t}
[/itex]

------------------------------------------------------------------------------------

So far so good. But he says if we add an external force F whereby x(t+Δt)=x(t) + FγΔt + l(t), then we get

[itex]
\frac{\partial \rho}{\partial t} = -γF\frac{\partial \rho}{\partial x}+D\frac{\partial^2 \rho}{\partial x^2}
[/itex]

But I am having trouble getting that result. It seems that the second derivative term is more like
[itex]
D+F^2γ^2Δt/2.
[/itex]

Seeing as this is my first Taylor series expansion in over 20 years, I was wondering if one of you could help me out and tell me if I have made a mistake. It shouldn't take you very long.

Thanks!
Phys.Org News Partner Physics news on Phys.org
Physicists discuss quantum pigeonhole principle
First in-situ images of void collapse in explosives
The first supercomputer simulations of 'spin?orbit' forces between neutrons and protons in an atomic nucleus

Register to reply

Related Discussions
Diffusion coefficient in diffusion equation and random walk ? Atomic, Solid State, Comp. Physics 3
Diffusion equation and neutron diffusion theory Nuclear Engineering 15
Replace external force with an internal force Classical Physics 2
No resultant force means no external force? Classical Physics 12
1D wave equation with dirac delta function as an external force. Differential Equations 7