Solving Diffusion PDE in a Hollow Cylinder

  • Context: Graduate 
  • Thread starter Thread starter wenzhe2092
  • Start date Start date
  • Tags Tags
    Diffusion Pde
Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
1 reply · 3K views
wenzhe2092
Messages
3
Reaction score
0
Dear all,

I'm trying to solve the diffusion PDE for my system, shown below:

[tex] \frac{\partial C}{\partial t} = D (\frac{\partial^2 C}{\partial r^2} + \frac{1}{r} \frac{\partial C}{\partial r})[/tex]

where C is the concentration, changing with time t and radius r. D is the diffusion coefficient.

I'm solving this using separation of variables, giving me two ODE.

[tex]T = Aexp (-\lambda^2 D t)[/tex]

where [tex]-\lambda^2[/tex] is the separation constant and A is the integration constant

[tex]R(r) = BJ_{0}(\lambda r) + CY_{0} (\lambda r)[/tex]

Principle solution given by:

[tex]C(r,t) = Aexp (-\lambda^2 D t) [BJ_{0}(\lambda r) + CY_{0} (\lambda r)][/tex]

My question is, given the boundary conditions and initial conditions of:

C(0.0135,t) = 0.433
C(0.0185,t) = 0
C(r,0) = 0.0398

How would i be able to solve this? I'm really stuck and any help would be appreciated. Note my system is a hollow cylinder.
 
Physics news on Phys.org
wenzhe2092 said:
Dear all,

Principle solution given by:

[tex]C(r,t) = Aexp (-\lambda^2 D t) [BJ_{0}(\lambda r) + CY_{0} (\lambda r)][/tex]

My question is, given the boundary conditions and initial conditions of:

C(0.0135,t) = 0.433
C(0.0185,t) = 0
C(r,0) = 0.0398

How would i be able to solve this? I'm really stuck and any help would be appreciated. Note my system is a hollow cylinder.


Look like the four unknowns A, B, C and [itex]\lambda[/itex] can just be reduced to three unknown. You only have three boundary and initial conditions.:wink: