summaryrefslogtreecommitdiff
path: root/839/CH6/EX6.3
diff options
context:
space:
mode:
authorpriyanka2015-06-24 15:03:17 +0530
committerpriyanka2015-06-24 15:03:17 +0530
commitb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
treeab291cffc65280e58ac82470ba63fbcca7805165 /839/CH6/EX6.3
downloadScilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip
initial commit / add all books
Diffstat (limited to '839/CH6/EX6.3')
-rwxr-xr-x839/CH6/EX6.3/Example_6_3.sce37
1 files changed, 37 insertions, 0 deletions
diff --git a/839/CH6/EX6.3/Example_6_3.sce b/839/CH6/EX6.3/Example_6_3.sce
new file mode 100755
index 000000000..c9884b1ca
--- /dev/null
+++ b/839/CH6/EX6.3/Example_6_3.sce
@@ -0,0 +1,37 @@
+//clear//
+clear;
+clc;
+
+//Example 6.3
+//Given
+pa = 2.7; //[atm]
+T = 288; //[K]
+D = 0.075; //[m]
+L = 70; //[m]
+Vbar = 60; //[m/s]
+M = 29;
+rh = D/4; //[m]
+mu = 1.74*10^-5 //[kg/m-s] Appendix 8
+rho_a = (29/22.4)*(2.7/1)*(273/288) // [kg/m^3]
+R = 82.056*10^-3;
+G = Vbar*rho_a //[kg/m^2-s]
+Nre = D*G/mu;
+kbyd = 0.00015*(0.3048/0.075);
+f = 0.0044; //[from Fig. 5.9]
+
+//Using Eq.(6.52)
+//pbar = 1.982; //[atm]
+//pb = 1.264; //[atm]
+err = 1;
+eps = 10^-3;
+pb = 1.5;
+
+while(err>eps)
+pbar = (pa+pb)/2;
+A = ((f*L/(2*rh))+log(pa/pb));
+pb_new = pa-(R*T*G^2/(pbar*29*101325))*A;
+err = pb-pb_new;
+pb = pb_new;
+end
+pb; //[atm]
+pbar = (pa+pb)/2 // [atm]