1 files changed, 34 insertions, 0 deletions
diff --git a/3511/CH8/EX8.7/Ex8_7.sce b/3511/CH8/EX8.7/Ex8_7.sce
new file mode 100644
index 000000000..c48b8ea8a
--- /dev/null
+++ b/3511/CH8/EX8.7/Ex8_7.sce
@@ -0,0 +1,34 @@
+clc;
+N=10000; // Speed in rpm
+Q=600; // Flow rate m^2/min
+rp=4; // Pressure ratio 
+eff_c=0.82; // Compressor efficiency
+T01=293; // Inlet temperature in kelvin
+p01=1.0; // Inlet pressure in bar
+Cp=1.005;// Specific heat   at constant pressure in kJ/kg K
+Cv=0.717;// Specific heat   at constant volume in kJ/kg K
+r=1.4; // Specific heat ratio 
+R=287; // Characteristic gas constant in J/kg K
+ca=60; // Axial velocity in m/s
+D2_D1=2 ;// Diameter ratio
+
+T_03=T01*rp^((r-1)/r);
+T03=T01+(T_03-T01)/eff_c;
+u2=sqrt (Cp*10^3*(T03-T01));
+Wc=u2^2; // Work of compression
+D2=(u2*60/(3.14*N));
+D1=D2/D2_D1;
+T1=T01-(ca^2/(2-Cp*10^3));
+p1=p01*(T1/T01)^(r/(r-1));
+row1=p1*10^5/(R*T1);
+Wroot=(Q/60)*(1/(ca*3.14*D1));
+u1=3.14*N*D1/60;
+alpha_root=atand (ca/u1);
+alpha_tip= atand (ca/u2);
+
+disp ("(i).Power input ");
+disp ("kW/kg/s",Wc/1000,"Wc = ");
+disp ("(ii).Impeller Diameters");
+disp ("m",D2,"D2 = ","m",D1,"D1 = ");
+disp ("(iii).Impeller and diffuser blade angles at inlet");
+disp ("degree",alpha_tip,"alpha_tip = ","degree",alpha_root,"alpha_root = ");