initial commit / add all books

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

+