Solving Diffusion PDE in a Hollow Cylinder

  • Context: Graduate 
  • Thread starter Thread starter wenzhe2092
  • Start date Start date
  • Tags Tags
    Diffusion Pde
Click For Summary
SUMMARY

The discussion focuses on solving the diffusion partial differential equation (PDE) in a hollow cylinder, represented by the equation \(\frac{\partial C}{\partial t} = D \left(\frac{\partial^2 C}{\partial r^2} + \frac{1}{r} \frac{\partial C}{\partial r}\right)\). The user employs the method of separation of variables, yielding the time-dependent solution \(T = A \exp(-\lambda^2 D t)\) and the radial component \(R(r) = B J_{0}(\lambda r) + C Y_{0}(\lambda r)\). Given the boundary conditions \(C(0.0135,t) = 0.433\), \(C(0.0185,t) = 0\), and \(C(r,0) = 0.0398\), the user seeks guidance on solving for the unknowns A, B, C, and \(\lambda\), noting that only three unknowns can be determined from the provided conditions.

PREREQUISITES
  • Understanding of partial differential equations (PDEs)
  • Familiarity with the method of separation of variables
  • Knowledge of Bessel functions, specifically \(J_{0}\) and \(Y_{0}\)
  • Basic concepts of boundary and initial value problems
NEXT STEPS
  • Study the application of boundary conditions in solving PDEs
  • Learn about Bessel function properties and their applications in cylindrical coordinates
  • Explore numerical methods for solving PDEs, such as finite difference methods
  • Investigate the uniqueness and existence theorems for solutions to PDEs
USEFUL FOR

Mathematicians, physicists, and engineers working on diffusion processes in cylindrical geometries, as well as students studying advanced topics in partial differential equations.

wenzhe2092
Messages
3
Reaction score
0
Dear all,

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

<br /> \frac{\partial C}{\partial t} = D (\frac{\partial^2 C}{\partial r^2} + \frac{1}{r} \frac{\partial C}{\partial r})<br />

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.

T = Aexp (-\lambda^2 D t)

where -\lambda^2 is the separation constant and A is the integration constant

R(r) = BJ_{0}(\lambda r) + CY_{0} (\lambda r)

Principle solution given by:

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

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:

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

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 \lambda can just be reduced to three unknown. You only have three boundary and initial conditions.:wink:
 

Similar threads

  • · Replies 3 ·
Replies
3
Views
4K
  • · Replies 36 ·
2
Replies
36
Views
6K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 4 ·
Replies
4
Views
4K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 5 ·
Replies
5
Views
3K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 8 ·
Replies
8
Views
4K
  • · Replies 5 ·
Replies
5
Views
2K