1 files changed, 29 insertions, 0 deletions
diff --git a/3685/CH13/EX13.10/Ex13_10.sce b/3685/CH13/EX13.10/Ex13_10.sce
new file mode 100644
index 000000000..af530b5b0
--- /dev/null
+++ b/3685/CH13/EX13.10/Ex13_10.sce
@@ -0,0 +1,29 @@
+clc
+//Given that
+v = 300 // Aircraft velocity in m/s
+p1 = 0.35 // Pressure in bar
+t1 = -40 // Temperature in degree centigrade
+rp = 10 // The pressure ratio of compressor 
+t4 = 1100 // Temperature of gases at turbine intlet in degree centigrade
+ma = 50 // Mass flow rate of air at the inlet of compressor in kg/s
+cp = 1.005 // Heat capacity of air at constant pressure in kJ/kg-K
+gama=1.4 // Ratio of heat capacities for air
+printf("\n Example 13.10 \n")
+T1 = t1+273
+T4 = t4+273
+T2 = T1 + (v^2)/(2*cp)*(10^-3)
+p2 = p1*(100)*((T2/T1)^(gama/(gama-1)))
+p3 = rp*p2
+p4 =p3
+T3 = T2*((p3/p2)^((gama-1)/gama))
+T5 = T4-T3+T2
+p5 = ((T5/T4)^(gama/(gama-1)))*(p4)
+p6 = p1*100
+T6 = T5*((p6/p5)^((gama-1)/gama))
+V6 = (2*cp*(T5-T6)*1000)^(1/2)
+Wp = ma*(V6-v)*v*(10^-6)
+Q1 = ma*cp*(T4-T3)*(10^-3)
+np = Wp/Q1
+printf("\n The temperature of the gases at the turbine exit is %f K,\n The pressure of the gases at the turbine exit is %f kN/m^2,\n The velocity of gases at the nozzle exit is %f m/sec,\n The propulsive efficiency of the cycle is %f percent",T5,p5,V6,np*100)
+//The answers vary due to round off error
+