6 files changed, 251 insertions, 0 deletions
diff --git a/1187/CH13/EX13.1/1.sce b/1187/CH13/EX13.1/1.sce
new file mode 100755
index 000000000..cb461c034
--- /dev/null
+++ b/1187/CH13/EX13.1/1.sce
@@ -0,0 +1,51 @@
+clc
+
+// Maximum hydraulic efficiency occurs for minimum pressure loss, that is, when
+
+// dp1/dQ=2.38Q-1.43=0
+
+Q_opt=1.43/2.38;
+
+p1_min=1.19*Q_opt^2-1.43*Q_opt+0.47; // MPa
+
+rho=1000; // kg/m^3
+g=9.81; // m/s^2
+w=69.1; // rad/s
+P=200*10^3; // W
+Ohm_P=0.565; // rad
+d=0.5; // m
+h=0.06; // m
+
+p1=p1_min*10^6/(rho*g); // mH2O, coversion of units
+
+H=(w*P^(1/2)/(rho^(1/2)*Ohm_P))^(4/5)/g;
+
+Hydraulic_efficiency=(H-p1)/H;
+disp("Hydraulic Efficiency =")
+disp(Hydraulic_efficiency)
+
+Overall_efficiency=P/(Q_opt*rho*g*H);
+disp("Overall Efficiency =")
+disp(Overall_efficiency)
+
+H_Euler=H-p1;
+
+u1=w*0.25;
+v_w1=g*H_Euler/u1;
+A=%pi*d*h*0.95; 
+v_r=Q_opt/A;
+
+alpha1=atand(v_r/v_w1);
+disp("Outlet angles of the guide vanes =")
+disp(alpha1)
+disp("degrees")
+
+beta1=atand(v_r/(v_w1-u1));
+disp("Rotor blade angle at inlet =")
+disp(beta1)
+disp("degrees")
+u2=w*0.325/2;
+beta2=atand(v_r/u2);
+disp("Rotor blade angle at outlet =")
+disp(beta2)
+disp("degrees")
diff --git a/1187/CH13/EX13.2/2.sce b/1187/CH13/EX13.2/2.sce
new file mode 100755
index 000000000..a237355e1
--- /dev/null
+++ b/1187/CH13/EX13.2/2.sce
@@ -0,0 +1,39 @@
+clc
+
+w=6.25;
+D=0.75; // m
+gv_angle=15; // guide vane angle in degrees
+g=9.81; // m/s^2
+H=27.5; // m
+A1=0.2; // m^2
+rho=1000; // kg/m^3
+p_atm=101.3*10^3;
+p_min=35*10^3;
+
+u1=%pi*w*D;
+v1=u1*sind(105)/sind(60);
+v_r1=v1*sind(gv_angle);
+v_w1=v1*cosd(gv_angle);
+v_w2=0;
+
+n_hydraulic=u1*v_w1/g/H;
+
+n_overall=0.97*n_hydraulic;
+disp("Overall efficiency =")
+disp(n_overall)
+
+Q=A1*v_r1;
+
+P=n_overall*Q*rho*g*H;
+Ohm_P=w*2*%pi/(g*H)^(5/4)*(P/rho)^(1/2);
+
+// sigma > 0.119*(0.5)^(1.84) = 0.0331
+
+sigma=0.0331;
+
+//((p_atm-p_min)/(rho*g)-z0)/H > 0.0331
+
+z0=((p_atm-p_min)/(rho*g))-sigma*H;
+disp("Limiting value for the height of the draft tube above =")
+disp(z0)
+disp("m")
diff --git a/1187/CH13/EX13.3/3.sce b/1187/CH13/EX13.3/3.sce
new file mode 100755
index 000000000..42a72e0fe
--- /dev/null
+++ b/1187/CH13/EX13.3/3.sce
@@ -0,0 +1,42 @@
+clc
+
+// Static head upstream = -11 mm H2O = -11*1000/1.2 mm air = -9.167 m air
+
+h=9.167; // m air
+g=9.81; // m/s^2
+d1=0.75; // m, tip diameters
+d2=0.4; // m, hub diameters
+d3=0.075; // m, diameter above atmospheric pressure
+d4=0.011; // m, diameter below atmospheric pressure
+P=6500; // W
+w=25;
+rho=1000; // kg/m^3
+
+v=sqrt(2*g*h); // Velocity upstream
+Q=%pi/4*d1^2*v; // Volume flow rate
+
+H=d3+d4; // Total head rise across fans
+p=rho*g*H;
+
+n_fan=Q*p/P;
+disp("Total efficiency =")
+disp(n_fan)
+
+p_ideal=p/n_fan;
+u=%pi*w*(d1+d2)/2;
+
+v_w2_A=p_ideal/(2*1.2*u);
+
+v1=Q/(%pi/4*(d1^2-d2^2  ));
+
+beta1_A=atand(v1/u);
+
+beta2_A=atand(v1/(u-v_w2_A));
+
+beta1_B=atand(v1/(u+v_w2_A));
+
+beta2_B=atand(v1/u);
+
+printf("Inlet angles for resp. fans %f & %f \n\n", beta1_A, beta1_B)
+
+printf("Outlet angles for resp. fans %f & %f", beta2_A, beta2_B)
diff --git a/1187/CH13/EX13.5/5.sce b/1187/CH13/EX13.5/5.sce
new file mode 100755
index 000000000..60deebf50
--- /dev/null
+++ b/1187/CH13/EX13.5/5.sce
@@ -0,0 +1,39 @@
+clc
+
+Q=0.04; // m^3/s
+d=0.15; // m
+h=28; // m
+f=0.006;
+l=38; // m
+g=9.81;
+fre=50; // Hz
+n_manometer = 0.75; 
+theta=30; // degrees
+
+v=Q/(%pi/4*d^2);
+h1=(3+4*f*l/d)*v^2/2/g; // Total head loss through pipes and valves
+
+h_m=h+h1; // Manometric head
+
+// w=2*%pi*50/n; where n = number of pairs of poles.
+// Ohm_s=w*Q^(1/2)/(g*H)^(3/4) = 0.876/n rad
+
+// If n = 2, Ohm_s = 0.438 rad, which suggests pump 1 or 2, and ω = 157 rad/s. Outlet flow area = %pi*D*D/10
+
+// v_r2=0.04/(%pi*D^2/10)
+// u2= ω*D/2 = 78.54 D
+
+// v_w2= g*h_m/(n_manometer*u2) = 5.06/D; // m^2/s
+
+// tan(theta) = v_r2/(u2-v_w2)
+
+// Solving above equation, we get
+// 78.54*D^3 - 5.06*D - 0.2205 = 0;
+
+// Solving above cubic equation we get
+
+D = 0.272; // m
+disp("D = ")
+disp(D)
+disp("m")
+disp("That is near enough. So we choose pump 1")
diff --git a/1187/CH13/EX13.6/6.sce b/1187/CH13/EX13.6/6.sce
new file mode 100755
index 000000000..8ac30ca1a
--- /dev/null
+++ b/1187/CH13/EX13.6/6.sce
@@ -0,0 +1,40 @@
+clc
+
+f=0.0085;
+l=21.1; // m
+d=0.09; // m
+g=9.81; // m/s^2
+rho=1000; // kg/m^3
+
+// h1=hf=(4*f*l/d)*(16*Q^2/(2*%pi^2*d^4*g)) = (100*Q)^2
+
+disp("(a)The head loss due to pipe friction in terms of flow rate Q is given as")
+disp("(100*Q)^2")
+
+// For Pump
+Q=[0:0.006:0.042 0.052];
+H=[15 16 16.5 16.5 15.5 13.5 10.5 7 0]
+plot(Q,H,"r")
+xlabel("Q(m^3/s)")
+ylabel("H(m)")
+
+// For Pipe System
+
+// H1 = 11.5 + (100*Q)^2;
+
+Q=[0:0.01:0.06];
+plot(Q,(11.5+10000*Q^2),"b")
+
+legend("pipe system", "pump")
+
+// From the plot of the pump and pipe characteristics, the intersection is at 
+
+H=16; // m
+Q=0.021; // m^3/s
+n=0.74;
+
+P=rho*g*H*Q/n;
+
+disp("(b)Power required =")
+disp(P)
+disp("W")
diff --git a/1187/CH13/EX13.7/7.sce b/1187/CH13/EX13.7/7.sce
new file mode 100755
index 000000000..dbf1abbde
--- /dev/null
+++ b/1187/CH13/EX13.7/7.sce
@@ -0,0 +1,40 @@
+clc
+
+H=16.5; // m
+Q=0.015; // m^3/s
+n=0.63;
+H_s=11.5;
+rho=1000; // kg/m^3
+g=9.81; // m/s^2
+
+h_f=(100*Q)^2; // frictional head loss
+
+h_valve = H - H_s - h_f;
+
+P=rho*g*H*Q/n;
+disp("(i) the power consumption of the pump =")
+disp(P/1000)
+disp("kW")
+
+disp("(ii) The power dissipated in the pump =")
+P_d=P*(1-n)/1000;
+disp(P_d)
+disp("kW")
+
+disp("(iii) The power lost by pipe friction =")
+P_f=rho*g*h_f*Q;
+disp(P_f/1000)
+disp("kW")
+
+disp("(iv) The power lost in the valve =")
+P_valve=rho*g*h_valve*Q;
+disp(P_valve/1000)
+disp("kW")
+
+P_s=rho*g*H_s*Q;
+
+n_overall = P_s/P*100;
+
+disp("(b) Overall efficiency of the installation =")
+disp(n_overall)
+disp("%")