1 files changed, 19 insertions, 0 deletions
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