Dif Finitas Numerical Fin.
Dif Finitas Numerical Fin.
McKelliget, 2002
= hP ( T T )
Divide the fin x axis into N points, at which the temperature is to be calculated. Since
this is a one-dimensional problem, and we have very powerful computers, separate these
node points by equal intervals, .
xi-1 xi xi+1
1
N
x i+1/2
= ( hP )i (Ti T )
i
dT
dT
kA
dx i +1 / 2
dx i 1 / 2
= ( hP )i (T i T )
( kA)i +1/2
( Ti+1 Ti )
( Ti Ti1 ) =
( kA )i 1 / 2
( hP )i (Ti T )
Rearranging gives
1/1
J. McKelliget, 2002
This equation can be applied to nodes i = 2 to N-1 and constitutes a finite difference
representation of the fin equation.
If necessary the thermal conductivity can be interpolated using the arithmetic mean,
ki +1 / 2 = 12 ( ki + ki +1 )
At the base of the fin (i=1) the first boundary condition is used
T1 = Tb
At the end of the fin the boundary condition
k
dT
= h ( T ( L) T )
dx x = L
( TN TN 1 ) =
( h ) N (T N T )
I=N
k N TN 1 + ( k N + hN ) TN = hN T
This approach is called theTaylor Series Finite Difference technique
2/2
J. McKelliget, 2002
Qconv
Qx,i-1/2
Qx,i+1/2
3/3