1 files changed, 25 insertions, 0 deletions
diff --git a/1226/CH21/EX21.3/EX21_3.sce b/1226/CH21/EX21.3/EX21_3.sce
new file mode 100755
index 000000000..e20bc1aeb
--- /dev/null
+++ b/1226/CH21/EX21.3/EX21_3.sce
@@ -0,0 +1,25 @@
+clc;funcprot(0);//EXAMPLE 21.3
+// Initialisation of Variables
+p1=1;.........//Pressure of air while entering the turbine in bar
+t1=300;........//Temperature of air entering the turbine in K
+p2=6.2;.........//Pressure of air after compression in bar
+etac=0.88;....//Efficiency of compressor
+etat=0.9;.....//Efficiency of turbine
+far=0.017;........//Fuel air ratio
+ga=1.4;........//Ratio of specific heats for compression
+gae=1.333;........//Ratio of specific heats for expansion
+cp=1.147;.............//Specific heat at constant pressure in kJ/kgK during expansion
+cpc=1.005;.............//Specific heat at constant pressure in kJ/kgK during compression
+C=44186;.............//Calorific value of fuel in kJ/kg
+//Calculations
+t2=t1*((p2/p1)^((ga-1)/ga));...............//Ideal temperature of air after compression in K
+t21=((t2-t1)/etac)+t1;..............//Actual temperature of air after compression in K
+t3=(((C*far)/((far+1)*cpc))+t21);..............//Temperature before expansion in turbine in K
+p4=p1;p3=p2;t4=t3*((p4/p3)^((gae-1)/gae));............//Ideal temperature after expansion in turbine in K
+t41=t3-(etat*(t3-t4));.................//Actual temperature after expansion in turbine in K
+wt=(cp*(t3-t41));........//Work done by turbine  in kJ/kg of air
+wc=round(1*cpc*(t21-t1));.................//Work done by compression in kJ/kg of air
+wnet=wt-wc;..........//Net work done in kJ/kg
+qs=(far)*C;................//Heat supplied in kJ/kg of air
+etath=wnet/qs;................//Thermal efficiency
+disp(etath*100,"Thermal efficiency in %:")