1 files changed, 44 insertions, 0 deletions
diff --git a/3511/CH6/EX6.4/Ex6_4.sce b/3511/CH6/EX6.4/Ex6_4.sce
new file mode 100644
index 000000000..1b4731b6d
--- /dev/null
+++ b/3511/CH6/EX6.4/Ex6_4.sce
@@ -0,0 +1,44 @@
+clc;
+eff_c=0.8; // Isentropic efficiency of compression each stage
+eff_CT=0.88; // Isentropic efficiency of compressor turbine
+eff_PT=0.88; // Isentropic efficiency of power turbine
+eff_trans=0.98; // Turbine to compressor transmission efficiency
+rp=3; // Pressure ratio in each stage of compression
+T08=297; // Temperature after intercooler in kelvin
+ma=15; // Air mass flow in kg/s
+eff_reg=0.8; // Regenerator effectiveness
+del_P=0.1; // Regenerator gas side pressure loss in bar
+T01=327; // Ambient temperature in kelvin
+p01=1; // Ambient pressure in bar
+T03=1000; // Maximum cycle temperature in kelvin
+CV=43.1; // Calorific value in MJ/kg
+Cpa=1.005;// Specific heat of air  at constant pressure in kJ/kg K
+Cpg=1.147;// Specific heat of fuel at constant pressure in kJ/kg K
+rg=1.33;// Specific heat ratio of fuel
+r=1.4; // Specific heat ratio of air
+p03=rp^2; // Pressre at state 3 in bar
+T_07=T01*(rp)^((r-1)/r);
+T07=T01+(T_07-T01)/eff_c;
+WLPC=ma*Cpa*(T07-T01); //  Work of low pressue compressor
+T_02=T08*(rp)^((r-1)/r);
+T02=T08+(T_02-T08)/eff_c;
+WHPC=ma*Cpa*(T02-T08);
+WC=WLPC+WHPC; // Compressor work
+WCa=WC/eff_trans; // Actual compressor work
+// Neglecting effect of mf
+T09=T03-(WCa/(ma*Cpg));
+T_09=T03-(T03-T09)/eff_PT;
+p09=p03/(T03/T_09)^(rg/(rg-1));
+p04=p01+del_P;
+T_04=T09*(p04/p09)^((rg-1)/rg);
+T04=T09-eff_PT*(T09-T_04);
+WTP=ma*Cpg*(T09-T04); // Work output of power turbine
+T05=T02+eff_reg*(T04-T02);
+mf=(ma*Cpg*(T03-T05))/(CV*10^3);
+sfc=mf*3600/(WTP);//Specifc fuel consumption
+eff_th=WTP/(mf*CV*10^3); // Thermal efficiency
+
+
+disp ("kW   (roundoff error)",WTP,"Work output of power turbine = ");
+disp ("kg/kW h",sfc,"Specifc fuel consumption = ");
+disp ("%",eff_th*100,"Thermal efficiency = ");