1 files changed, 31 insertions, 0 deletions
diff --git a/3511/CH9/EX9.5/Ex9_5.sce b/3511/CH9/EX9.5/Ex9_5.sce
new file mode 100644
index 000000000..b44bb889c
--- /dev/null
+++ b/3511/CH9/EX9.5/Ex9_5.sce
@@ -0,0 +1,31 @@
+clc;
+Q=1000; // Flow rate of free air in m^3/min
+P1=0.98; // Inlet pressure in bar
+T1=15+273; // Inlet temperature in kelvin
+Dm=0.6; // Mean diameter in m
+h=6.75; // blade length in cm
+CL=0.6; CD=0.05; // At zero angle of incidence
+Cp=1.005; // Specific heat in kJ/kg K
+r=1.4; // Specific heat ratio
+R=287; // Characteristic gas constant in J/kg K
+k=1-0.1; //Blade occupys 10% of axial area
+N=6000; // speed in rpm
+Ac=19.25*10^-4; // Projected area in m^2
+n=50;
+eff_C=1; // Efficiency of compressor
+
+row=(P1*10^5)/(R*T1); // Density
+A=k*3.14*Dm*h*10^-2; // Area of axial
+ca=Q/(60*A); // Axial velocity
+u=3.14*Dm*N/60; // Blade velocity
+beta_1=atand(u/ca); // Blade angle at inlet
+w=sqrt (ca^2+u^2); // From velocity triangle
+L=CL*row*w^2*Ac/2;
+D=CD*row*w^2*Ac/2;
+P=(L*cosd(beta_1)+D*sind (beta_1))*u*n*10^-3; // Power input / stage
+m=Q*row/60;// mass flow rate
+rp=((P*eff_C/(m*Cp*T1))+1)^(r/(r-1)); // pressure ratio
+P2=rp*P1; // Pressure
+
+disp ("kW   (Roundoff error )",P,"Power input/stage = ");
+disp ("bar",P2,"Pressure at outlet = ");