Dear all,(adsbygoogle = window.adsbygoogle || []).push({});

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 seperation 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 Forums - The Fusion of Science and Community**

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Diffusion PDE solve

Loading...

Similar Threads - Diffusion solve | Date |
---|---|

A Pseudospectral method using FFT to solve diff equation | Sep 20, 2016 |

Numerically solving system of four PDEs | Jan 21, 2015 |

Solve Fick's second law of diffusion | Aug 5, 2012 |

Solving Diffusion equation with Convection | Apr 27, 2011 |

Is there a general method for solving advection-diffusion equations? | Apr 22, 2011 |

**Physics Forums - The Fusion of Science and Community**