summaryrefslogtreecommitdiff
path: root/629/CH13
diff options
context:
space:
mode:
Diffstat (limited to '629/CH13')
-rw-r--r--629/CH13/EX13.1/ex13_1.txt3
-rw-r--r--629/CH13/EX13.1/example13_1.sce22
-rw-r--r--629/CH13/EX13.10/ex13_10.txt3
-rw-r--r--629/CH13/EX13.10/example13_10.sce20
-rw-r--r--629/CH13/EX13.3/example13_3.sce28
-rw-r--r--629/CH13/EX13.4/ex13_4.txt3
-rw-r--r--629/CH13/EX13.4/example13_4.sce16
-rw-r--r--629/CH13/EX13.5/ex13_5.txt2
-rw-r--r--629/CH13/EX13.5/example13_5.sce19
-rw-r--r--629/CH13/EX13.6/ex13_6.txt2
-rw-r--r--629/CH13/EX13.6/example13_6.sce16
-rw-r--r--629/CH13/EX13.7/ex13_7.txt3
-rw-r--r--629/CH13/EX13.7/example13_7.sce12
-rw-r--r--629/CH13/EX13.8/ex13_8.txt2
-rw-r--r--629/CH13/EX13.8/example13_8.sce8
-rw-r--r--629/CH13/EX13.9/ex13_9.txt3
-rw-r--r--629/CH13/EX13.9/example13_9.sce17
17 files changed, 179 insertions, 0 deletions
diff --git a/629/CH13/EX13.1/ex13_1.txt b/629/CH13/EX13.1/ex13_1.txt
new file mode 100644
index 000000000..98ed2bc53
--- /dev/null
+++ b/629/CH13/EX13.1/ex13_1.txt
@@ -0,0 +1,3 @@
+
+The volume rate of flow = 29.7 m^3/s.
+ \ No newline at end of file
diff --git a/629/CH13/EX13.1/example13_1.sce b/629/CH13/EX13.1/example13_1.sce
new file mode 100644
index 000000000..b23f5d559
--- /dev/null
+++ b/629/CH13/EX13.1/example13_1.sce
@@ -0,0 +1,22 @@
+clear
+clc
+//Example 13.1 DISCHARGE FROM VELOCITY DATA
+r=[0 5 10 15 20 25 30 35 40 45 47.5 50]; //[cm]
+V=[50 49.5 49 48 46.5 45 43 40.5 37.5 34 25 0]; //velocity[m/s]
+
+//dA=2*pi*r*dr, Q=V*dA
+del_A(1)=0; //dA
+q(1)=0;
+for i=2:1:11
+ del_r(i)=(r(i+1)-r(i-1))/2;//dr
+ del_A(i)=2*%pi*r(i)*del_r(i)*10^-4; //[m^2]
+ q(i)=V(i)*del_A(i); //[m^3/s]
+end
+//for i=12,
+ del_r(12)=(r(12)-r(11)/2);
+ del_A(12)=2*%pi*r(12)*del_r(12)*10^-4; //[m^2]
+ q(12)=V(12)*del_A(12); //[m^3/s]
+
+//Discharge
+Q=sum(q) //m^3/s
+printf("\nThe volume rate of flow = %.1f m^3/s.\n",Q) \ No newline at end of file
diff --git a/629/CH13/EX13.10/ex13_10.txt b/629/CH13/EX13.10/ex13_10.txt
new file mode 100644
index 000000000..360821622
--- /dev/null
+++ b/629/CH13/EX13.10/ex13_10.txt
@@ -0,0 +1,3 @@
+
+The uncertainty of the calculated discharge = 0.0046 m^3/s.
+ \ No newline at end of file
diff --git a/629/CH13/EX13.10/example13_10.sce b/629/CH13/EX13.10/example13_10.sce
new file mode 100644
index 000000000..71455724e
--- /dev/null
+++ b/629/CH13/EX13.10/example13_10.sce
@@ -0,0 +1,20 @@
+clear
+clc
+//Example 13.10 UNCERTAINTY ESTIMATE FOR AN ORIFICE METER
+K=0.66;
+d=0.15; //[m]
+A=%pi*d^2/4 //area[m^2]
+D=0.24; //[m]
+g=9.81; //[m/s^2]
+h=0.25; //[m]
+S=13.6; //specific gravity of Hg
+del_h=h*(S-1) //piezometric head[m]
+Q=K*A*sqrt(2*g*del_h)
+//From equation of uncertainty, (UQ/Q)^2=(UK/K)^2+(2Ud/d)^2+(Uh/2h)^2
+//from example 13.2
+UK=0.03;
+Ud=0.00015;//[m]
+Uh=0.01;//[m]
+//Uncertainty in Q
+UQ=Q*sqrt((UK/K)^2+(2*Ud/d)^2+(Uh/(2*h))^2)
+printf("\nThe uncertainty of the calculated discharge = %.4f m^3/s.\n",UQ) \ No newline at end of file
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
diff --git a/629/CH13/EX13.4/ex13_4.txt b/629/CH13/EX13.4/ex13_4.txt
new file mode 100644
index 000000000..8568fc440
--- /dev/null
+++ b/629/CH13/EX13.4/ex13_4.txt
@@ -0,0 +1,3 @@
+
+The deflection on the manometer = 1.1 ft.
+ \ No newline at end of file
diff --git a/629/CH13/EX13.4/example13_4.sce b/629/CH13/EX13.4/example13_4.sce
new file mode 100644
index 000000000..61ee6759b
--- /dev/null
+++ b/629/CH13/EX13.4/example13_4.sce
@@ -0,0 +1,16 @@
+clear
+clc
+//Example 13.4 MANOMETER DEFLECTION FOR AN ORIFICE METER
+d=8/12; //[ft]
+Ao=%pi*d^2/4 //area [ft^2]
+v=1.22*10^-5; //[ft^2/s]
+Q=2; //flow rate [ft^3/s]
+g=32.2; //[ft/s^2]
+Re=4*Q/(%pi*d*v) //Reynolds number
+//interpolate for d/D=8/12 from fig. 13.14
+K=0.68
+del_h=Q^2/(2*g*K^2*Ao^2) //change in piezometric head [ft]
+//del_h=del_l(g_w-g_air)/g_w
+//g_w>>g_a
+del_l=del_h //manometer deflection [ft]
+printf("\nThe deflection on the manometer = %.1f ft.\n",del_l) \ No newline at end of file
diff --git a/629/CH13/EX13.5/ex13_5.txt b/629/CH13/EX13.5/ex13_5.txt
new file mode 100644
index 000000000..732f57a6b
--- /dev/null
+++ b/629/CH13/EX13.5/ex13_5.txt
@@ -0,0 +1,2 @@
+
+The mass flow rate of natural gas = 0.302 kg/s. \ No newline at end of file
diff --git a/629/CH13/EX13.5/example13_5.sce b/629/CH13/EX13.5/example13_5.sce
new file mode 100644
index 000000000..59606cf2c
--- /dev/null
+++ b/629/CH13/EX13.5/example13_5.sce
@@ -0,0 +1,19 @@
+clear
+clc
+//Example 13.5 MASS FLOW RATE OF NATURAL GAS
+d0=0.07; //[m]
+d1=0.1; //[m]
+A0=%pi*d0^2/4 //area[m^2]
+A1=%pi*d1^2/4 //area[m^2]
+rho=0.678; //density of methane[kg/m^3]
+k=1.31;
+v=1.59*10^-5; //[m^2/s]
+p1=101000; //pressure[Pa]
+delp=10000; //difference in pressure[Pa]
+p2=p1-delp //[Pa]
+ReK=(d0/v)*sqrt(2*delp/rho) //ReK=Re/K
+//From fig 13.14,
+K=0.7; //flow coefficient
+Y=1-((1/k)*(1-p2/p1)*(0.41+0.35*(A0/A1)^2)) //Compressibility factor
+m=Y*A0*K*sqrt(2*rho*(p1-p2)) //mass flow rate[kg/s]
+printf("\nThe mass flow rate of natural gas = %.3f kg/s.\n",m) \ No newline at end of file
diff --git a/629/CH13/EX13.6/ex13_6.txt b/629/CH13/EX13.6/ex13_6.txt
new file mode 100644
index 000000000..747059965
--- /dev/null
+++ b/629/CH13/EX13.6/ex13_6.txt
@@ -0,0 +1,2 @@
+
+The discharge of water at 10°C is 0.268 m^3/s. \ No newline at end of file
diff --git a/629/CH13/EX13.6/example13_6.sce b/629/CH13/EX13.6/example13_6.sce
new file mode 100644
index 000000000..aa5acb49b
--- /dev/null
+++ b/629/CH13/EX13.6/example13_6.sce
@@ -0,0 +1,16 @@
+clear
+clc
+//Example 13.6 FLOW RATE USING A VENTURI METER
+del_p=35000; //pressure difference [Pa]
+Gamma=9810; //specific weight[N/m^3]
+del_h=del_p/Gamma //change in piezometric head [m]
+d=0.2; //[m]
+v=1.31*10^-6; //[m^2/s]
+g=9.81; //[m/s^2]
+ReK=(d/v)*sqrt(2*g*del_h) //(ReK=Re/K)
+//Interpolating from fig.13.14
+K=1.02; //flow coefficient
+A=%pi*d^2/4 //area [m^2]
+//Discharge
+Q=1.02*A*sqrt(2*g*del_h) //[m^3/s]
+printf("\nThe discharge of water at 10°C is %.3f m^3/s.\n",Q) \ No newline at end of file
diff --git a/629/CH13/EX13.7/ex13_7.txt b/629/CH13/EX13.7/ex13_7.txt
new file mode 100644
index 000000000..b5e8e9b96
--- /dev/null
+++ b/629/CH13/EX13.7/ex13_7.txt
@@ -0,0 +1,3 @@
+
+The discharge of water over the weir = 0.23 m^3/s.
+ \ No newline at end of file
diff --git a/629/CH13/EX13.7/example13_7.sce b/629/CH13/EX13.7/example13_7.sce
new file mode 100644
index 000000000..8d600774b
--- /dev/null
+++ b/629/CH13/EX13.7/example13_7.sce
@@ -0,0 +1,12 @@
+clear
+clc
+//Example 13.7 FLOW RATE FOR A RECTANGULAR WEIR
+H=0.21; //head on weir[m]
+P=0.6; //height of weir[m]
+L=1.3; //width[m]
+g=9.81; //acceleration due to gravity[m/s^2]
+//Flow coefficient
+K=0.4+0.05*(H/P)
+//Discharge
+Q=K*L*sqrt(2*g*H^3) //[m^3/s]
+printf("\nThe discharge of water over the weir = %.2f m^3/s.\n",Q) \ No newline at end of file
diff --git a/629/CH13/EX13.8/ex13_8.txt b/629/CH13/EX13.8/ex13_8.txt
new file mode 100644
index 000000000..3a0a2fe2e
--- /dev/null
+++ b/629/CH13/EX13.8/ex13_8.txt
@@ -0,0 +1,2 @@
+
+ The flow of water over the weir = 0.096 m^3/s. \ No newline at end of file
diff --git a/629/CH13/EX13.8/example13_8.sce b/629/CH13/EX13.8/example13_8.sce
new file mode 100644
index 000000000..86af6e81b
--- /dev/null
+++ b/629/CH13/EX13.8/example13_8.sce
@@ -0,0 +1,8 @@
+clear
+clc
+//Example 13.8 FLOW RATE FOR A TRIANGULAR WEIR
+H=0.43; //head on weir[m]
+g=9.81; //[m/s^2]
+//Discharge
+Q=0.179*sqrt(2*g*(H^5)) //[m^3/s]
+printf("\n The flow of water over the weir = %.3f m^3/s.\n",Q) \ No newline at end of file
diff --git a/629/CH13/EX13.9/ex13_9.txt b/629/CH13/EX13.9/ex13_9.txt
new file mode 100644
index 000000000..a5d3e2441
--- /dev/null
+++ b/629/CH13/EX13.9/ex13_9.txt
@@ -0,0 +1,3 @@
+
+The mass flow rate of air flowing through a venturi meter = 0.0264 kg/s.
+ \ No newline at end of file
diff --git a/629/CH13/EX13.9/example13_9.sce b/629/CH13/EX13.9/example13_9.sce
new file mode 100644
index 000000000..d169454c0
--- /dev/null
+++ b/629/CH13/EX13.9/example13_9.sce
@@ -0,0 +1,17 @@
+clear
+clc
+//Example 13.9 COMPRESSIBLE FLOW
+k=1.4;
+p1=150000; //upstream pressure[Pa]
+p2=100000; //throat pressure[Pa]
+T1=300; //temperature[K]
+R=287; //[J/Kg.K]
+//Ideal gas law
+rho1=p1/(R*T1) //[Kg/m^3]
+D1=0.03; //[m]
+D2=0.01; //[m]
+A2=%pi*D2^2/4 //area[m^2]
+Cd=1;
+//Mass flow rate
+m=Cd*A2*((p2/p1)^(1/k))*{([2*k/(k-1)]*p1*rho1*[1-(p2/p1)^((k-1)/k)])/(1-(p2/p1)^(2/k)*(D2/D1)^4)}^(1/2) //[Kg/s]
+printf("\nThe mass flow rate of air flowing through a venturi meter = %.4f kg/s.\n",m) \ No newline at end of file
generated by cgit (git 2.34.1) at 2024-12-19 12:43:29 +0000