summaryrefslogtreecommitdiff
path: root/965/CH9/EX9.18/18.sci
blob: 2a9409f4551bda86336d5068317ec3fcb470a518 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
clc;
clear all;
disp("heat transfer rate")
L=3.2;//m
d=0.006;//m diameter
ts=54;// degree C
tsat=100;//degree C
rhol=973.7;//kg/m^3
mu=365*10^(-6);// Ns/m^2
k=0.668;// W/m.C
rhov=0.596;// kg/m^3
hfg=2257*10^3;// J/kg
g=9.81;//m/s
disp("h=0.0077*(rhol*(rhol-rhov)*k^2*g/(mu^2))^0.333*Re^0.4")
disp("Eliminating h from euqation we get the condition that the flow will be turbulent if ")
disp("0.00296*((rhol*(rhol-rhov)*k^3*g*(tsat-ts)^3*L^3/(mu^5*hfg^3))^(5/9)>1800")
x=0.00296*(rhol*(rhol-rhov)*k^3*g*(tsat-ts)^3*L^3/(mu^5*hfg^3))^(5/9);
if(x>1800)
Re=x
h=0.0077*(rhol*(rhol-rhov)*k^2*g/(mu^2))^0.333*Re^0.4
Q=h*L*1*(tsat-ts);
disp("kW/m",Q/1000,"hear transfer rate per unit width =")