1 files changed, 21 insertions, 0 deletions
diff --git a/1052/CH16/EX16.4/164.sce b/1052/CH16/EX16.4/164.sce
new file mode 100755
index 000000000..67e214f25
--- /dev/null
+++ b/1052/CH16/EX16.4/164.sce
@@ -0,0 +1,21 @@
+clc;
+//Example 16.4
+//page no 187
+printf("Example 16.4 page no 187\n\n");
+//if flow for both phases is laminar then cal pressure drop total
+//a. apply eq. 16.22b to obtain Y_G
+X=1.66
+Y_G=10-10.405*X+8.6786*X^2-0.9167*X^3
+printf("\n Y_G=%f ",Y_G);
+//pressure drop from eq 16.1
+P_drop_g=2.71
+P_drop=Y_G*P_drop_g
+printf("\n pressure drop P_drop=%f psf/100 ft",P_drop);
+//b. apply eq 16.23b to generate Y_l
+Y_l=6.4699*X^-0.556
+printf("\n Y_l =%f ",Y_l);
+//pressure drop from eq. 16.2
+P_drop_l=7.50
+P_drop_b=Y_l*P_drop_l
+printf("\n pressure drop P_drop_b=%f psf/100 ft",P_drop_b);
+