Solve Fick's second law of diffusion

abstracted6
Messages
37
Reaction score
0
I'm curious how to solve Fick's second law of diffusion \frac{∂c}{∂t}=D \frac{∂^2c}{∂x^2}For conditions:c(x,0)=0c(0,t)=Ac(\infty,t)=0Physically this means:
-c(x,t) is the concentration at point x at time t.
-Initially there is no concentration of diffusing species.
-At x=0 for all t the is a constant concentration "a".
-As x goes to infinity for all time, the concentration is 0.
-D is the diffusivity, assume it is a constant.

The solution is:c(x,t)=A erfc(\frac{x}{2\sqrt{Dt}})
What method was used to arrive at that solution?
 
Physics news on Phys.org
Hi !

May be this formula was obtained thanks to the Laplace method for PDE resolution.
Normally we would have to use the double Laplace transform (relatively to x AND t), which would be rather arduous.
But the PDE and boundary conditions are simple enough to use the usual single Laplace transform (relatively to t only).
 

Attachments

  • Laplace.JPG
    Laplace.JPG
    38.9 KB · Views: 6,573
Possibly try a similarity solution?
 
Transport Phenomena by Bird, Stewart, and Lightfoot show how to solve this (singular perturbation boundary layer problem) using similarity solutions. Look for the analogous viscous flow startup problem.
 
Back
Top