1 files changed, 30 insertions, 0 deletions
diff --git a/3407/CH8/EX8.3/Ex8_3.sce b/3407/CH8/EX8.3/Ex8_3.sce
new file mode 100644
index 000000000..e2f102cc9
--- /dev/null
+++ b/3407/CH8/EX8.3/Ex8_3.sce
@@ -0,0 +1,30 @@
+clear;
+clc;
+funcprot(0);
+
+//given data
+Z = 12;//number of vanes
+delW = 230;//in kW
+T01 = 1050;//stagnation temperature in K
+mdot = 1;//flow rate in kg/s
+eff_ts = 0.81;//total-to-static efficiency
+Cp = 1.1502;//in kJ/(kg.K)
+gamma = 1.333;
+R = 287;//gas constant
+
+//Calculations
+S = delW/(Cp*T01);
+alpha2 = (180/%pi)*acos(sqrt(1/Z));
+beta2 = 2*(90-alpha2);
+p3_p01 = (1-(S/eff_ts))^(gamma/(1-gamma));
+M02 = sqrt((S/(gamma-1))*((2*cos(beta2*%pi/180))/(1+cos(beta2*%pi/180))));
+M2 = sqrt((M02^2)/(1-0.5*(gamma-1)*(M02^2)));
+U2 = sqrt((gamma*R*T01)*(1/cos(beta2*%pi/180))*(S/(gamma-1)));
+
+//Results
+printf('(i) The absolut and relative flow angles:\n alpha2 = %.2f deg\n beta2 = %.2f deg',alpha2,beta2);
+printf('\n (ii) The overall pressure ratio = %.3f',p3_p01);
+printf('\n (iii) The rotor rip speed = %.1f m/s\n The inlet absolute Mach number = %.3f',U2,M2);
+
+
+//there are small errors in the answers given in textbook