Second Fick's law with nonconstant diffusion coefficient

  • Thread starter mabiondi
  • Start date
  • #1
2
0
Hello everybody.
I should solve a modified version of second Fick's law in nondimensional spherical coordinates; t is the time and rho [0,1] the nondimensional radius. In this equation the diffusion coefficient is vraiable with t and rho.

The equation is the following:

dC/dt = (1/rho^2)*d/d(rho)(rho^2*D(rho,t)*dC/d(rho))

Initial conditions:
C(0,rho) = a*rho^4 + b*rho^2 -(a + b);

Boundary conditions:
dC(t,0)/d(rho) = 0
C(t,1) = 0

Is it possible to solve such an equation by Laplace transforms?
Alternatively, is there a user friendly tool for numerical solution of a PDE like this?

Thanks in advance, cheers.
Marco
 

Answers and Replies

  • #2
20,844
4,543
I don't think an analytic solution exists. It would have to be solved numerically. I would just convert it to a set of ODE's using the method of lines, and then solve the equations numerically in FORTRAN using an automatic (stiff) integration package.

Chet
 
  • Like
Likes mabiondi
  • #3
2
0
Thanks a lot for the suggestion, I really appreciate it. I'll try to develop a routine using MOL.
Cheers,
Marco
 

Related Threads on Second Fick's law with nonconstant diffusion coefficient

Replies
1
Views
2K
  • Last Post
Replies
0
Views
4K
  • Last Post
Replies
1
Views
5K
  • Last Post
Replies
6
Views
21K
  • Last Post
Replies
2
Views
2K
  • Last Post
Replies
1
Views
2K
Replies
4
Views
1K
  • Last Post
Replies
1
Views
1K
Replies
7
Views
2K
Replies
1
Views
4K
Top