summaryrefslogtreecommitdiff
path: root/2510/CH17/EX17.13/Ex17_13.sce
blob: 4bf8728ee06ba6ca6c078fdae9cc2ed080ca36ba (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
//Variable declaration:
w = 1                           //Width of single of fin (m)
t = 2/10**3                     //Fin base thickness (m)
l = 6/10**3                     //Fin length thickness (m)
T1 = 250                        //Surface temperature ( C)
T2 = 20                         //Ambient air temperature ( C)
h = 40                          //Surface convection coefficient (W/m^2.K)
k = 240                         //Thermal conductivity of fin (W/m.K)

//Calculation:
Ab = t*w                        //Base area of the fin (m^2)
Te = T1-T2                      //Excess temperature at the base of the fin (K)
Qw = h*Ab*Te                    //Heat transfer rate without a fin (W)
Af = 2*w*(sqrt(l**2-(t/2)**2))  //Fin surface area (m^2)
Qm = h*Af*Te                    //Maximum heat transfer rate (m^2)
Bi = h*(t/2)/k                  //Biot number
Lc = l                          //Corrected length (m)
Ap = l*t/2                      //Profile area (m^2)
A = sqrt((Lc**3*h)/k*Ap)        //Abscissa for the fin efficiency figure
//From figure 17.4:
nf = 0.99                       //Fin efficiency
Qf = nf*Qm                      //Fin heat transfer rate (W)
R = Te/Qf                       //Fin thermal resistance (K/W)
E = Qf/Qw                       //Fin effectiveness
Qm = round(Qm*10**-1)/10**-1

//Result:
printf("1. The heat transfer rate without the fin is : %.1f W .",Qw)
printf("2. The maximum heat transfer rate from the fin is : %f W .",Qm)
printf("3. The fin efficiency is : %.0f %%",nf*100)
printf("   The fin thermal resistance is : %.1f C/W .",R)
printf("   The fin effectiveness is : %.1f .",E)