1 files changed, 25 insertions, 0 deletions
diff --git a/1226/CH3/EX3.38/EX3_38.sce b/1226/CH3/EX3.38/EX3_38.sce
new file mode 100755
index 000000000..88c24902a
--- /dev/null
+++ b/1226/CH3/EX3.38/EX3_38.sce
@@ -0,0 +1,25 @@
+clc;funcprot(0);//EXAMPLE 3.38
+// Initialisation of Variables
+t1=300;.................//Inlet temperature in K
+p1=1;....................//Inlet pressure in bar
+ma=1;....................//Mass of air in kg
+rp=6.25;.............//Pressure ratio
+t3=1073;...........//Maximum temperature in K
+etac=0.8;............//Efficiency of compressor
+etat=0.8;.............//Efficiency of turbine
+ga=1.4;.................//Ratio of specific heats
+cp=1.005;.............//Specific heat at constant pressure in kJ/kgK
+//Calculations
+t2=t1*(rp^((ga-1)/ga));...........//Ideal Temperature of air while leaviing the compressor in K
+t21=((t2-t1)/etac)+t1;............//Actual Temperature of air while leaviing the compressor in K
+Wcomp=ma*cp*(t21-t1);.............//Compressor work in kJ/kg
+t4=t3/(rp^((ga-1)/ga));........//Ideal temperature of air while leaving the turbine in K
+t41=t3-(etat*(t3-t4));..........//Actual temperature of air while leaving the turbine in K
+Wtur=ma*cp*(t3-t41);..............//Turbine work in kJ/kg
+Wnet=Wtur-Wcomp;.................//Net work produced in kJ/kg
+Qs=ma*cp*(t3-t21);.................//Heat supplied in kJ/kg
+disp(Wcomp,"Compressor work in kJ/kg:")
+disp(Wtur,"Turbine work in kJ/kg:")
+disp(Qs,"Heat supplied in kJ/kg:")
+disp((Wnet/Qs)*100,"Cycle efficiency in %:")
+disp(t41,"Actual exhaust temperature of turbine in K")