General Heat Conduction Equations: Rectangular/cartesian Coordinates
General Heat Conduction Equations: Rectangular/cartesian Coordinates
Rate of heat Input + rate of Generation = rate of heat Output + rate of heat Accumulation
Rectangular/cartesian coordinates:
FKKKSA
FKKKSA
FKKKSA
T
k 2T
q
2T q
2
2
x cp x cp
x cp
conduction in 3 dimensions:
T
x
FKKKSA
2T 2T 2T q
2 2 2
x y z cp
T T T q
x y z c
2
Cartesian coordinates :
Cylindrical coordinates:
r r
T
t
T 1 2T 2T q 1 T
r 2 2 2
r r
z k t
Spherical coordinates :
1 2 T 1
T
sin
r
2
2
r r r r sin
k a bT
FKKKSA
1
2T q 1 T
r 2 sin 2 2 k t
k
= thermal diffusivity
cp
x 0 T = T1
x L T = T2
T T T q
x y z c
2
FKKKSA
T
t
d2T
0
2
dx
T
0
2
x
2
d dT
dx dx
1st integration
dT
d dx
0dx
2nd integration
FKKKSA
dT c dx
1
dT
dx
C1
Tx C1 x C 2
T C (0) C
1
1
2
C 2 T1
T2 C1 ( L) T 1
Boundary Condition
x 0 T = T1
x L T = T2
FKKKSA
C1
Tx
Chem. Eng. Dept
T2 T1
L
T2 T1
L
x T1
dT
qx Ak
dx
qx dx Ak dT
x2
T2
x1
T1
qx(x1 x2 ) Ak(T1 T2 )
(T1 T2 )
qx Ak
(x2 x1)
qx
FKKKSA
(T1 T2 )
(x2 x1 )
Ak
FKKKSA
FKKKSA
FKKKSA
FKKKSA
cp
T T T q
x y z c
2
c p
T q
x c
0
2
p
2
d T
dx
FKKKSA
Boundary condition
BC1.
BC2.
x o
x L
d T
dT
=0
dx
T = T1
dx
q dx
dT
dx
k
0
FKKKSA
0 C
dx
C1 0
d T
dx
dT
1st integration:
x C1
q
dT x C dx
k
1
Tx
Boundary condition
BC1.
BC2.
x o
x L
dT
=0
dx
T = T1
C2
FKKKSA
T1
2k
2k
2k
2
x C1 x C 2
L 0 L C2
C1 0
L T1
Temp. Profile:
Chem. Eng. Dept
Tx
2k
q
r r
T 1 2T 2 T q
1 T
r 2 2 2
r r
z k t
r r
T 1 2T 2T q 1 T
r 2 2 2
r r
z k t
0
T
r 0
r r
FKKKSA
dr
dT
r 0
dr
dr
dT
r 0
dr
First integration:
d
r
Second integration:
FKKKSA
dT
r
0dr
dr
dT
dr
C1
dT
C1
r
dr
Tr C1 ln r C 2
Tr C1 ln r C 2
Boundary condition:
BC1
r = r1
BC2
T = T1
T1 C1 ln r1 C 2
r = r2
T = T2
T 2 C1 ln r 2 C 2
T2 T1
T T
(ln r ln r ) T
ln r
ln
Temp. Profile:
Tr
r1
r1
FKKKSA
C 2 T1
r2
T2 T1
ln
r2
r1
ln r1
dT
qr Ak
dr
dT 1 T1 T2
dr r ln r1
r2
A 2rL
1 T1 T
qr (2rL)k
r ln r1
r2
FKKKSA
rearrange
T1 T2
qr
r2
ln
r1
2Lk
FKKKSA
r
r
1
r
FKKKSA
r r
T
sin
r
sin
1
1
2
sin
q
k
1 T
r
r
1
r
r r
T
sin
sin
r
1
r
FKKKSA
2 T
r
0
r
Chem. Eng. Dept
1
2
sin
0
1 T
d 2 dT
r
0
dr dr
Boundary condition:
BC1 r = r1
BC2 r = r2
T = T1
T = T2
dT
dr
Second integration:
Tr
FKKKSA
C1
r
C2
dT
r 2
0dr
dr
C1
dT
C1
r
dr
T1
Solving simultaneously gives:
C1
Temp. Profile:
FKKKSA
Tr
C1
r1
C2
T2
r
C1
r2
C2
C2
T T
C2
r1r2
r2 T2 r1 T1
r1r2
T = T1
T = T2
Tr
C1
r1
r r2 r1
T2
r2 T2 r1 T1
r2 r1
r2 r1
dT
qr Ak
dr
Fouriers
law
A 4 r
r1r2
qr 4 r k
qr
FKKKSA
4 k
T1 T2
1
r2
r1
1
r1
Chem. Eng. Dept
dT
dr
T T
1
T1 T2
qr
1 1
r1 r2
4 k
r1r2
r2 r1
T T
1
FKKKSA
FKKKSA
FKKKSA