diff options
| author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
|---|---|---|
| committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
| commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
| tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /629/CH13/EX13.3/example13_3.sce | |
| parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
| download | Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2 Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip | |
initial commit / add all books
Diffstat (limited to '629/CH13/EX13.3/example13_3.sce')
| -rw-r--r-- | 629/CH13/EX13.3/example13_3.sce | 28 |
1 files changed, 28 insertions, 0 deletions
diff --git a/629/CH13/EX13.3/example13_3.sce b/629/CH13/EX13.3/example13_3.sce new file mode 100644 index 000000000..3a23e8846 --- /dev/null +++ b/629/CH13/EX13.3/example13_3.sce @@ -0,0 +1,28 @@ +clear
+clc
+//Example 13.3 ANALYSIS OF AN ORIFICE METER
+g=9.81; //[m/s^2]
+l=0.25; //deflection[m]
+S=13.6; //specific gravity of Hg
+h=l*(S-1) //piezometric head[m]
+d=0.15; //[m]
+D=0.24; //[m]
+Ao=%pi*d^2/4 //[m^2]
+A1=%pi*D^2/4 //[m^2]
+v=10^-6; //[m^2/s]
+//Flow coefficient
+ReK=(d/v)*sqrt(2*g*h) //ReK=Re/K
+//From fig.13.14, for Re/K, d/D=0.625
+K=0.66;
+Q=K*Ao*sqrt(2*g*h) //[m^3/s]
+printf("\nThe discharge in the system = %.3f m^3/s.\n",Q)
+Cv=0.98;
+//K=Cv*Cc/(1-(Cc*Ao/A1)^2)
+P=[K*(Ao/A1)^2 Cv -K] //polynomial in Cc
+q=roots(P) //roots of P
+Cc=q(2) //positive root
+V1=Q/A1 //[m/s]
+V2=Q/(Cc*Ao) //[m/s]
+//Head loss
+hL=(V2-V1)^2/(2*g) //[m]
+printf("\nThe head loss produced by the orifice = %.2f m.\n",hL)
\ No newline at end of file |