1 files changed, 21 insertions, 0 deletions
diff --git a/3511/CH6/EX6.14/Ex6_14.sce b/3511/CH6/EX6.14/Ex6_14.sce
new file mode 100644
index 000000000..f07812768
--- /dev/null
+++ b/3511/CH6/EX6.14/Ex6_14.sce
@@ -0,0 +1,21 @@
+clc;
+T01=38+273; // Inlet temperature of compressor in kelvin
+eff_c=0.82; // Compressor efficiency
+T03=650+273; // Turbine inlet temperature in kelvin
+eff_T=0.8; // Turbine efficiency
+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
+
+t=T03/T01;
+// For maximun specific work we know that
+ropt=(sqrt (t*eff_c*eff_T))^(r/(r-1));
+T_02=T01*ropt^((r-1)/r);
+T02=T01+(T_02-T01)/eff_c;
+T_04=T03/ropt^((rg-1)/rg);
+T04=T03-eff_T*(T03-T_04);
+eff_th=((Cpg*(T03-T04))-(Cpa*(T02-T01)))/(Cpg*(T03-T02));
+
+disp (ropt,"Optimum pressure ratio = ");
+disp ("%",eff_th*100, "Overall efficiency = ");