View Full Version : Analytical solution for heat equation with simple boundary conditions
uqjzhang
May18-11, 01:00 AM
I am trying to solve the following heat equation ODE:
d^2T/dr^2+1/r*dT/dr=0 (steady state) or
dT/dt=d^2T/dr^2+1/r*dT/dr (transient state)
The problem is simple: a ring with r1<r<r2, T(r1)=T1, T(r2)=T2.
I have searched the analytical solution for this kind of ODEs in polar coordinate systems, found a lot of solutions for more difficult problems :devil:
Can anybody help me with this one please?
JJacquelin
May18-11, 04:10 AM
The stady state problem :
d^2T/dr^2+1/r*dT/dr=0 is a simple linear ODE. With boundary conditions T(r1)=T1, T(r2)=T2, did you find the solution ?
The transiant problem :
dT/dt=d^2T/dr^2+1/r*dT/dr is a linear PDE. Boundary conditions T(r1)=T1, T(r2)=T2 are not enough. A unique solution requires more conditions related with time (for exemple what is the state at time 0, and/or how the T(r1), T(r2) are estabished as a function of time, or...)
Various methods are possible, with more or less difficulties depending on the initial conditions.
If the initial conditions are not stated, one can find the general solutions on the form of a series of particular solutions T(r,t)=Sum(Ak*Fk(r)*Gk(t)) where Ak are constants, Fk are Bessel functions or r and Gk is exponential functions of t. Other representations are possible.
So, first you have to clarify the initial conditions.
JJacquelin
May21-11, 01:51 AM
I found it difficult to solve the equation?
Anyone of you can help me?
Steady state solution in attachment :
JJacquelin
May21-11, 02:47 AM
I found it difficult to solve the equation?
Anyone of you can help me?
Transient solutions in attachment :
vBulletin® v3.8.7, Copyright ©2000-2012, vBulletin Solutions, Inc.