blob: 0ad13de119d6929b4d2506562c389f4a65dda90f (
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
|
clc;
clear;
//Example 3.34
rho=960.63 //Density in [kg/m^3]
Cp=4.216*10^3 //Specific heat in [J/(kg.K)]
D=16 //Diameter in [cm]
D=D/100 //[m]
k=0.68 //Thermal conductivity in [W/m.K]
A=(%pi*(D/2)^2)
L=A/(%pi*D) //Length=A/P in [m]
Beta=0.75*10^-3 //[K^-1]
alpha=1.68*10^-7 //[m^2/s]
g=9.81 //[m/s^2]
Tw=403 //[K]
T_inf=343 //[K]
dT=Tw-T_inf //[K]
v=0.294*10^-6 //[m^2/s]
Nra=(g*Beta*(L^3)*dT)/(v*alpha)
//1.For Top surface
Nnu=0.15*(Nra)^(1.0/3.0) //Nusselt number
ht=Nnu*k/L //Heat transfer coeff for top surface in W/(m^2.K)
ht=round(ht)
//2.For bottom surface
Nnu=0.27*Nra^(1.0/4.0) //Nusselt number
hb=Nnu*k/L //[W/sq m.K]
hb=round(hb)
Q=(ht+hb)*A*dT; //[W]
printf("The rate of heat input is %f W",Q)
|
