1 files changed, 34 insertions, 0 deletions
diff --git a/3821/CH7/EX7.21/Example7_21.sce b/3821/CH7/EX7.21/Example7_21.sce
new file mode 100644
index 000000000..29010c20c
--- /dev/null
+++ b/3821/CH7/EX7.21/Example7_21.sce
@@ -0,0 +1,34 @@
+///Chapter No 7 Fluid Mechanics
+///Example 7.21 Page No:131-132
+///Find Bernoulli's equation
+///Input data
+clc;
+clear;
+L1=300;                  //Length of pipe in m
+D11=0.9;                 //Diameter at higher end in m
+D12=0.6;                 //Diameter at lower end in m
+S=0.85;                 //Specific gravity 
+Q1=0.08;                 //Flow in l/s
+P1=40*10^3;            //Pressure at higher end in kPa
+pi1=3.14;                //pi constant 
+rho=1000;               //Rho constant
+g1=9.81;                 //Gravity constant
+slop=1/50;              //1 in 50
+
+
+//Calculation
+//Datum line is passing through the center of the low end,therefore
+Z2=0; 
+Z1=slop*L1;
+//Q=A1*V1=A2*V2        Continuity eqation
+V11=Q1/((pi1/4)*(D11^2)); //Frome continuity eqation, discharge
+V12=Q1/((pi1/4)*(D12^2));     //Frome continuity eqation, discharge
+///Bernoulli's equation 
+P2=(((((P1/(rho*S*g1))+(V11^2/(2*g1))+Z1)-(V12^2/(2*g1))+Z2))*(S*rho*g1))*10^-3;
+
+
+///Output
+printf('Z1= %f m \n ',Z1);
+printf('continuity eqation discharge V11= %f m^3 \n ',V11);
+printf('continuity eqation, discharge V12= %f m^3 \n',V12);
+printf('bernoullis equation= %f Kpa \n ',P2);