1 files changed, 32 insertions, 0 deletions
diff --git a/965/CH9/EX9.17/17.sci b/965/CH9/EX9.17/17.sci
new file mode 100644
index 000000000..c6112b5a0
--- /dev/null
+++ b/965/CH9/EX9.17/17.sci
@@ -0,0 +1,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 =")
+