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
|
clc;
clear all;
disp("square plate")
x=0.4;//m
L=0.75;//mm
ts=28;// degree C
rhol=993.95;//kg/m^3
k=62.53*10^(-2);//W/m.C
mu=728.15*10^(-6);// kg/ms
hfg=2402*10^3;// J/kg
rhov=0.561;//kg/m^3
g=9.81;// m/s
tsat=42;// degree C
delx=(4*k*mu*(tsat-ts)*x/(g*rhol*(rhol-rhov)*hfg))^0.25;
disp("mm",delx*1000,"The film thickness at the bottom of plate =")
hx=k/delx;
disp("W/m^2.C",hx,"heat transfer coefficient =")
delL=(4*k*mu*(tsat-ts)*L/(g*rhol*(rhol-rhov)*hfg))^0.25;
disp("mm",delL*1000,"The film thickness at the bottom of plate =")
hL=4*k/(3*delL);
disp("W/m^2.C",hL,"heat transfer coefficient =")
h=1.2*hL;
disp("W/m^2.C",h,"The overall heat transfer coefficient =")
Q=h*L*L*(tsat-ts);
disp("kW",Q/1000,"rate of heat transfer per metre width, Q =")
m=Q/hfg;//kg/hr
disp("kg/hr",m*3600,"Condensate flow rate =")
hinc=h*(sin(%pi*25/180))^0.25;
disp("W/m^2.C",hinc,"heat transfer coefficient when the plate is inclined 25 degree with the horizontal")
Re=4*m/(mu*L);
disp(Re,"Re =")
|
