initial commit / add all books

1 files changed, 35 insertions, 0 deletions

diff --git a/49/CH7/EX7.3/ex3.sce b/49/CH7/EX7.3/ex3.sce
new file mode 100755
index 000000000..f7f3dbdeb
--- /dev/null
+++ b/49/CH7/EX7.3/ex3.sce
@@ -0,0 +1,35 @@
+//CHAPTER 7_ Flow Measurement

+//Caption : Gross volume flow rate(venturi)

+// Example 3// Page 438

+dp=0.02     //('entering the diameter of the line in which water is flowing=:')

+dt=0.01    //('entering the diameter of venturi=:')

+B=0.5;   // given

+// The discharge coefficients remains in the flat portion of the curve for reynolds numbers 10^4 to 10^6 Cd=0.95

+u=8.6*10^-4    //('entering the viscosity=:')

+Cd=0.95;

+Rn_min=10^4;

+disp("Minimum flow rate is given by:")

+disp("mdot_min=%pi*dp*u*Rn_min/4")

+mdot_min=%pi*dp*u*Rn_min/4

+g=9.81;

+printf('Minimum flow rate at 25 deg cent is %1.3f kg/s\n',mdot_min)

+pf=1000   // density of water

+At=78.53*10^-6    //('entering the throat area=:')

+pm=13.6      //('entering the density of manometer  fluid=:')

+

+//h is the height of mercury column due to flow

+disp("To calculate the mercury reading corresponding to minimum flow, using-")

+disp("h_min=((mdot_min*sqrt(1-B^4))/((sqrt(2*g*(pm-pf/pf))*pf*At*Cd)))^2;")

+h_min=((mdot_min*sqrt(1-B^4))/((sqrt(2*g*(pm-pf/pf))*pf*At*Cd)))^2;

+//in mm

+H_min=h_min*1000

+printf('So the pressure reading observed for the given flow ratre is %1.1f  mm of Hg\n',H_min)

+h_max=.25    //('entering the  value of h maximum=:')

+m_max=(pf*At*Cd*sqrt(2*g*(pm-pf/pf))*sqrt(h_max))/sqrt(1-B^4);

+printf('The maximum flow rate is %1.1f kg/s\n',m_max)

+

+

+

+

+

+