Solve Fick's second law of diffusion

Click For Summary
SUMMARY

This discussion focuses on solving Fick's second law of diffusion, represented by the equation ∂c/∂t = D ∂²c/∂x². The initial conditions specified are c(x,0) = 0, c(0,t) = A, and c(∞,t) = 0, indicating no initial concentration, a constant concentration at x=0, and zero concentration as x approaches infinity. The solution derived is c(x,t) = A erfc(x/(2√(Dt))). The Laplace method, particularly the single Laplace transform with respect to time, is identified as the technique used to obtain this solution, with references to the book "Transport Phenomena" by Bird, Stewart, and Lightfoot for further insights.

PREREQUISITES
  • Understanding of partial differential equations (PDEs)
  • Familiarity with the Laplace transform method
  • Knowledge of error functions (erfc)
  • Concepts of diffusivity in physical systems
NEXT STEPS
  • Study the Laplace transform method for solving PDEs in detail
  • Explore similarity solutions in the context of transport phenomena
  • Review the error function and its applications in diffusion problems
  • Investigate the viscous flow startup problem as described in "Transport Phenomena" by Bird, Stewart, and Lightfoot
USEFUL FOR

Researchers, physicists, and engineers working in fields related to diffusion processes, particularly those interested in mathematical modeling and analysis of concentration changes over time and space.

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,592
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.
 

Similar threads

  • · Replies 8 ·
Replies
8
Views
4K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 5 ·
Replies
5
Views
4K
  • · Replies 5 ·
Replies
5
Views
2K
Replies
3
Views
2K
  • · Replies 1 ·
Replies
1
Views
538
  • · Replies 2 ·
Replies
2
Views
4K
  • · Replies 2 ·
Replies
2
Views
2K