Diffusion equation in 1D

poiuy

The solution to the diffusion equation in 1D may be written as follows:



n'(x,t) = N/sqrt(4piDt) * exp(-x^2/4DT)

where n'(x,t) is the concentration of the particles at position x at time t, N is the total number of particles and D is the diffusion coefficient.

Write down an expression for the number of particles in a slab of thickness dx located at position x.




I assumed it would be the integral of the function between x and x+dx with respect to x.

However exp(-x^2/4Dt) can't be integrated between these values. I have a standard integral for exp(-ax^2) which is 0.5 sqrt (pi/a) but this only applies to integrating between zero and infinity.


If anybody could point me in the right direction it would be greatly appreciated, I think I am missing something obvious here and this is a really simple question.

Thanks

nrqed

Welcome to the forums!

Note that the integrale of any function f(x) between x and x +dx is simply f(x) dx!

[tex] \int_x^{x+dx} f(x') dx' \approx f(x) dx [/tex]

poiuy

Wow thanks, incredible that I could have had 14 years of education and never been taught that, thanks very much.

poiuy

Actually thinking about it, it's incredible that I couldn't work that out for myself.

tomostler13

not sure if this is in exactly the right place but its quite late and it is on 1D diffusion. I am given the task of producing a working code for a non-linear 1D diffusion equation given below, it models the diffusion of "hard" spheres in a solution, in a vessel of a particular height. Using the information below i have to construct a finite difference equation and apply neumann boundary conditions, i can obtain the finite difference equations but dont quite understand neumann boundary conditions. The equation is



i think this FDE is correct, please feel free to correct me rho_c is just a packing constant. I am stuck as to where to go to derive my boundary conditions

any help would be greatly appreciated