1 files changed, 52 insertions, 0 deletions
diff --git a/32/CH19/EX19.01/19_01.sce b/32/CH19/EX19.01/19_01.sce
new file mode 100755
index 000000000..eed30b1e3
--- /dev/null
+++ b/32/CH19/EX19.01/19_01.sce
@@ -0,0 +1,52 @@
+//pathname=get_absolute_file_path('19.01.sce')
+//filename=pathname+filesep()+'19.01-data.sci'
+//exec(filename)
+//Specific heat of gases(in kJ/kg.K):
+Cpg=1.13 //kJ/kg.K
+Cpa=1.005 //kJ/kg.K
+rg=1.33
+ra=1.4
+C=41.84*10^3 //kJ/kg of fuel
+//Temperatures(in K):
+T1=272
+T3=1000
+//Compression efficiency:
+nc=0.84
+p3=3
+p2=3
+p1=0.5
+p5=0.4
+//Turbine efficiency:
+nt=0.82
+//Nozzle efficiency:
+nn=0.92
+//Speed(in m/s):
+Ca=200
+//Temperature at state 2(in K):
+T2=T1*(p2/p1)^((ra-1)/ra)
+//Temperature at state 2'(in K):
+T2a=T1+(T2-T1)/nc
+//Compressive work(in kW):
+Wc=Cpa*(T2a-T1)
+printf("\n RESULT \n")
+printf("\nPower required for compressor = %f kW/kg",Wc)
+//Air fuel ratio:
+r=(C-Cpg*T3)/(Cpg*T3-Cpa*T2a)
+printf("\nAir fuel ratio = %f",r)
+//Temperature at state 4'(in K):
+T4a=T3-Cpa/Cpg*(T2a-T1)/(1+r)
+T4a=810.46
+T4=T3-(T3-T4a)/nt
+//Pressure of gas leaving turbine(in bar):
+p4=p3*(T4/T3)^(rg/(rg-1))
+printf("\nPressure of gas leaving turbine = %f bar",p4)
+//Temperature at state 5(in K):
+T5=T4a*(p5/p4)^((rg-1)/rg)
+//Temperature at state 5'(in K):
+T5a=T4a-nn*(T4a-T5)
+//Exit jet velocity(in m/s):
+C5a=sqrt(2*Cpg*(T4a-T5a)*10^3)
+Ce=C5a
+//Thrust per kg of air per second:
+T=(1+1/r)*Ce-Ca
+printf("\nThrust = %f N/kg/s",T)
+\ No newline at end of file