1 files changed, 64 insertions, 0 deletions
diff --git a/905/CH4/EX4.6/4_6.sce b/905/CH4/EX4.6/4_6.sce
new file mode 100755
index 000000000..8d5a2eb31
--- /dev/null
+++ b/905/CH4/EX4.6/4_6.sce
@@ -0,0 +1,64 @@
+clear;
+clc;
+
+// Illustration 4.6
+// Page: 255
+
+printf('Illustration 4.6 -  Page: 255\n\n');
+
+// solution
+//*****Data*****//
+Ff = 0.9; // [foaming factor]
+sigma = 70; // [liquid surface tension, dyn/cm]
+Do = 5; // [mm]
+//From Example 4.4
+// X = 0.016;
+p = 15 // [pitch, mm]
+// From equ 4.35
+// Ah/Aa = A
+A = 0.907*(Do/p)^2; // [ratio of vapor hole area to tray active area]
+
+// Assume 
+t = 0.5; // [m]
+// From equ 4.32
+alpha = 0.0744*t+0.01173;
+beeta = 0.0304*t+0.015;
+
+// Since X<0.1, therefore
+X = 0.1;
+// From equ 4.31
+Cf = alpha*log10(1/X) + beeta; 
+// Since Ah/Aa > 0.1, therefore
+Fha = 1;
+Fst = (sigma/20)^0.2; // [surface tension factor]
+// From equ 4.30
+C = Fst*Ff*Fha*Cf;
+
+// From Example 4.4
+rowg = 1.923; // [kg/cubic m]
+rowl = 986; // [kg/cubic m]
+Qg = 1.145; // [cubic m/s]
+// From equation 4.29
+vgf = C*(sqrt((rowl-rowg)/rowg)); // [m/s]
+// Since X<0.1
+// Equ 4.34 recommends Ad/At = B = 0.1
+B = 0.1;
+// For an 80% approach to flooding, equation 4.33 yields
+f = 0.8; 
+D = sqrt((4*Qg)/(f*vgf*%pi*(1-B))); // [m]
+// At this point, the assumed value of tray spacing ( t = 0.5 m) must be   // checked against the recommended values of Table 4.3. Since the calculated
+// value of D < 1.0 m, t = 0.5 m is the recommended tray spacing, and no
+// further iteration is needed.
+
+deff('[y] = f14(Q)','y = B-((Q-sin(Q))/(2*%pi))');
+Q = fsolve(1.5,f14);
+Lw = D*sin(Q/2); // [m]
+rw = D/2*cos(Q/2); // [m]
+
+At = %pi/4*D^2; // [total cross sectional area, square m]
+Ad = B*At; // [Downcomer area, square m]
+Aa = At-2*Ad; // [ Active area over the tray, square m]
+Ah = 0.101*Aa; // [Total hole area, square m]
+
+printf('Summarizing, the details of the sieve-tray design are as follows:\n\n');
+printf(" Diameter = %f m\n Tray spacing = %f m\n Total cross-sectional area  = %f square m\n Downcomer area = %f square m\n Active area over the tray = %f square m\n Weir length = %f m\n Distance from tray center to weir = %f m\n Total hole area = %f square m\n Hole arrangement: 5 mm diameter on an equilateral-triangular pitch 15 mm between hole centers, punched in stainless steel sheet metal 2 mm thick\n\n",D,t,At,Ad,Aa,Lw,rw,Ah); 
+\ No newline at end of file