1 files changed, 34 insertions, 0 deletions
diff --git a/3821/CH7/EX7.18/Example7_18.sce b/3821/CH7/EX7.18/Example7_18.sce
new file mode 100644
index 000000000..c16fd79c9
--- /dev/null
+++ b/3821/CH7/EX7.18/Example7_18.sce
@@ -0,0 +1,34 @@
+///Chapter No 7 Fluid Mechanics
+///Example 18 Page No:129
+//// Find rate of water flow l/min
+//Input data
+clc;
+clear;
+L1=200;               //length of pipe in m
+D11=1;                //Diameter at high end in m
+D12=0.4;              //Diameter at low end in m
+P1=50*10^3;         //Pressure at high end  in kPa
+Q1=4000;              //Rate of water flow l/min
+S=1;                 //Slope of pipe 1 in 100
+Z2=0;                //Datum line is passing through the center of the low end,therefore
+pi1=3.14;
+rho=1000;
+g1=9.81;
+
+
+///Calculation
+Q1=(4000*10^-3)/60;   //rate of water flow l/min in m**3/s
+Z1=1/100*L1;           //slope of pipe 1 in 100 is in m
+//Q=A1*V11=A2V2        //continuity eqation ,discharge
+V11=Q1/((pi1/4)*(D11^2));//in m^3
+V12=Q1/((pi1/4)*(D12**2));//in m^3
+//bernoulli's equation 
+P2=(((((P1/(rho*g1))+(V11^2/(2*g1))+Z1)-(V12^2/(2*g1))-Z2))*(rho*g1))*10^-3; 
+
+
+///Output
+printf('rate of water flow=%f m^3/s \n ',Q1);
+printf('slope of pipe=%f m \n',Z1);
+printf('continuity eqation ,discharge= %f m^3 \n ",V11);
+printf('continuity eqation ,discharge= %f m^3 \n ",V12);
+printf('bernoullis equation for discharge= %f kpa \n ',P2);