From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001
From: prashantsinalkar
Date: Tue, 10 Oct 2017 12:27:19 +0530
Subject: initial commit / add all books

---
 3511/CH6/EX6.15/Ex6_15.sce | 28 ++++++++++++++++++++++++++++
 1 file changed, 28 insertions(+)
 create mode 100644 3511/CH6/EX6.15/Ex6_15.sce

(limited to '3511/CH6/EX6.15/Ex6_15.sce')

diff --git a/3511/CH6/EX6.15/Ex6_15.sce b/3511/CH6/EX6.15/Ex6_15.sce
new file mode 100644
index 000000000..08a920784
--- /dev/null
+++ b/3511/CH6/EX6.15/Ex6_15.sce
@@ -0,0 +1,28 @@
+clc;
+p01=1; // Stagnation pressure at entry in bar
+pa=0.93; // Static pressure at entry in bar
+T1=10+273;// Static temperature in entry in kelvin
+p02=6; // Pressure at state 2 in bar
+T02=230+273; // Temperature at state 2 in kelvin
+P=5100; // Turbine output power in kW
+A=0.1; // Compressor entry area in m^2
+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
+R=287; // Characteristic constant in J/kg K
+T04=460+273; // Exhaust pipe temperature in kelvin
+
+M=sqrt (((p01/pa)^((r-1)/r)-1)/((r-1)/2));
+T01=T1*(1+(r-1)/2*M^2);
+T_02=T01*(p02/p01)^((r-1)/r);
+eff_c=(T_02-T01)/(T02-T01);
+row_s=(pa*10^5)/(R*T1);
+a=sqrt (r*R*T1);
+V=M*a;
+m=row_s*A*V;
+T03=(P/(m*Cpg))+T04;
+
+disp ("%",eff_c*100,"Compressor efficiency = ");
+disp ("kg/s",m,"Mass flow rate = ");
+disp ("K   (roundoff error)",T03,"Turbine inlet stagnation temperature = ");
-- 
cgit