1 files changed, 39 insertions, 0 deletions
diff --git a/2510/CH12/EX12.5/Ex12_5.sce b/2510/CH12/EX12.5/Ex12_5.sce
new file mode 100755
index 000000000..b35245a04
--- /dev/null
+++ b/2510/CH12/EX12.5/Ex12_5.sce
@@ -0,0 +1,39 @@
+//Variable declaration:
+T1 = 99.0                           //Mean film temperature (°C)
+T2 = 98.0                           //Plate surface temperature (°C)
+g = 9.807                           //Gravitational acceleration (m/s^2)
+//From Appendix:
+T3 = 100.0                          //Saturation temperatre (°C)
+h_vap1 = 970.3                      //Latent heat of steam in Btu/lb (Btu/lb)
+h_vap2 = 2.255*10**6                //Latent heat of steam in J/kg (J/kg)
+p_v = 0.577                         //Density of steam (kg/m^3)
+p_l = 960.0                         //Density of liquid water condensate (kg/m^3)
+mu_l = 2.82*10**-4                  //Absolute viscosity of liquid water condensate (kg/m.s)
+k = 0.68                            //Thermal conductivity of water (W/m.K)
+//From table 12.2
+Z = 0.4                             //Height of rectangular plate (m)
+Pw = 0.2                            //Wetted perimeter of rectangular plate (m)
+syms h                              //Average heat transfer coefficient (W/m^2.K)
+
+//Calculation:
+A = Z*Pw                            //Heat transfer area of plate (m^2)
+R = A/Pw                            //Ratio A/Pw (m)
+v_l = mu_l/p_l                      //Kinematic viscosity of liquid water condensate (m^2/s)
+Co1 = (h/k)*(v_l**2/g/(1-p_v/p_l))**(1/3)    //Condensation number (in terms of the average heat transfer coefficient)
+Re = 4*h*Z*(T3-T2)/(mu_l*h_vap2)    //Reynolds number in terms of the average heat transfer coefficient
+//From equation 12.14:
+CO1 = 0.0077*Re**Z                  //Co in terms of Reynolds number for flow type 1
+x1 = solve(h,Co1-CO1)               //Solving heat transfer coefficient (W/m^2.K)
+h1 =x1(2);                          //Average heat transfer coefficient for flow type 1 (W/m^2.K)
+Re1 = subst(h1,h,Re)                //Reynolds number for flow type 1
+CO2 = 1.874*Re**(-1/3)              //Co in terms of Reynolds number for flow tupe 2
+x2 = solve(Co1-CO2,h)               //Solving average heat transfer coefficient for flow type 2 (W/m^2.K)
+h2 = x2(1);                         //Average heat transfer coefficient for flow type 2 (W/m^2.K)
+Re2 = subst(h2,h,Re)                 //Reynolds number for flow type 2
+h2 = round(h2*10**-1)/10**-1
+
+//Result:
+printf("The type of condensation flow type 2 is laminar.")
+disp("And the condensation heat transfer coefficient is : ")
+disp(h2)
+disp("W/m^2.K.")