1 files changed, 28 insertions, 0 deletions
diff --git a/3511/CH9/EX9.4/Ex9_4.sce b/3511/CH9/EX9.4/Ex9_4.sce
new file mode 100644
index 000000000..e37bf1934
--- /dev/null
+++ b/3511/CH9/EX9.4/Ex9_4.sce
@@ -0,0 +1,28 @@
+clc;
+rp=4; // Total head pressure ratio
+eff_O=0.85; // Overall total head isentropic efficiency
+T01=290; // Total head inlet temperature in kelvin
+alpha_1=10; // Inlet air angle in degree
+alpha_2=45; // Outlet air angle in degree
+u=220; // Blade velocity in m/s
+wf=0.86; // Wok done factor
+R=284.6; // Characteristic gas constant in kJ/kg K
+Cp=1005; // Specific heat in J/kg K
+r=1.4; // Specific heat ratio
+
+eff_P=1/(log10(((rp^((r-1)/r)-1)/eff_O)+1)/(log10(rp)*((r-1)/r)));; 
+// From velocity triangle
+ca=u/(tand(alpha_1)+tand(alpha_2)); // Axial velocity
+Del_Tstage=wf*u*ca*(tand(alpha_2)-tand(alpha_1))/Cp; // Stage temperature rise
+T02=T01*(rp)^((r-1)/(r*eff_P)); // Total head temperature 
+T02_T01=T02-T01; // Total temperature rise
+n=T02_T01/Del_Tstage; // Total number of stages
+// from velocty traingles
+w1=ca/cosd(alpha_2);
+c1=ca/cosd(alpha_1);
+T1=T01-c1^2/(2*Cp); // Static temperature
+M=w1/sqrt(r*R*T1); // Mach number at inlet
+
+disp (eff_P*100,"Polytropic efficiency of the compressor = ");
+disp (n,"Total number of stages = ");
+disp (M,"Mach number at inlet = ");