summaryrefslogtreecommitdiff
path: root/1073/CH3/EX3.33/3_33.sce
blob: 9735551d9202ea11f1020159c4767fecfdf727af (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
clc;
clear;
//Example 3.33
k=0.03406   //[W/(m/K)]
Beta=2.47*10^-3 //K^-1
Npr=0.687   //Prandtl number
v=26.54*10^-6   //m^2/s
g=9.81  //[m/s^2]
Tw=523  //[K]
T_inf=288   //[K]
dT=Tw-T_inf //[K]
D=0.3048    //[m]
Ngr=(g*Beta*dT*(D^3))/(v^2)     //Grashof number
Nra=Ngr*Npr 
//For Nra less than 10^9,we have for horizontal cylinder
Nnu=0.53*(Nra^(1.0/4.0))    //Nusselt number
h=Nnu*k/D   //[W/sq m.K]
Q_by_l=h*%pi*D*dT;  //W/m
printf("Heat loss of heat transfer per meter lengh is %f W/m",Q_by_l);