Solve Fick's second law of diffusion

Click For Summary

Discussion Overview

The discussion revolves around solving Fick's second law of diffusion, specifically the equation \(\frac{∂c}{∂t}=D \frac{∂^2c}{∂x^2}\) under given boundary conditions. Participants explore various methods for finding a solution, including mathematical techniques and references to literature.

Discussion Character

  • Technical explanation, Mathematical reasoning, Debate/contested

Main Points Raised

  • One participant presents the equation and boundary conditions, proposing a solution in the form of \(c(x,t)=A \, \text{erfc}\left(\frac{x}{2\sqrt{Dt}}\right)\).
  • Another participant suggests that the solution may have been obtained using the Laplace method for partial differential equations, noting the potential complexity of a double Laplace transform.
  • A different participant proposes considering a similarity solution as an alternative approach.
  • Another participant references the book "Transport Phenomena" by Bird, Stewart, and Lightfoot, indicating that it discusses solving similar problems using similarity solutions and suggests looking at analogous viscous flow startup problems.

Areas of Agreement / Disagreement

Participants present multiple methods and approaches for solving the problem, indicating that there is no consensus on a single solution method. The discussion remains open with various competing views on how to approach the problem.

Contextual Notes

Participants express different mathematical techniques and references, but there are no settled assumptions or resolutions regarding the best method to solve the equation.

Who May Find This Useful

This discussion may be useful for individuals interested in mathematical methods for solving partial differential equations, particularly in the context of diffusion processes in physics and engineering.

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,603
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
3K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 5 ·
Replies
5
Views
4K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 0 ·
Replies
0
Views
3K
Replies
3
Views
2K
  • · Replies 1 ·
Replies
1
Views
863
  • · Replies 2 ·
Replies
2
Views
4K