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/CH5/EX5.1/Ex5_1.sce | 26 ++++++++++++++++++++++++++
 1 file changed, 26 insertions(+)
 create mode 100644 3511/CH5/EX5.1/Ex5_1.sce

(limited to '3511/CH5/EX5.1')

diff --git a/3511/CH5/EX5.1/Ex5_1.sce b/3511/CH5/EX5.1/Ex5_1.sce
new file mode 100644
index 000000000..9dba042b6
--- /dev/null
+++ b/3511/CH5/EX5.1/Ex5_1.sce
@@ -0,0 +1,26 @@
+clc;
+p1=1; // Pressure before compression in bar
+T1=350; // Temperature before compression in kelvin
+T3=2000; // Temperature after combustion in kelvin
+rp=1.3; // Pressure ratio
+Cp=1.005; // Specific heat at constant pressure in kJ/kg K
+r=1.4; // Specific heat ratio
+R=287; // Characteristic gas constant in J/kg K
+
+T2=T1*(rp)^((r-1)/r); // Temperature at the end of the compression
+T4=T3*(1/rp)^((r-1)/r); // Temperature after expansion
+Wc=Cp*(T2-T1); // Work done during compression
+WT=Cp*(T3-T4); // Work done during expansion
+WN=WT-Wc; // Net work done
+p2=rp*p1; // Pressure at state 2
+p3=p2; p4=p1; // Constant pressure process
+V1=R*T1/(p1*10^5); // specific Volume at state 1
+V2=R*T2/(p2*10^5); // specific Volume at state 2
+V3=R*T3/(p3*10^5); // specific Volume at state 3
+V4=R*T4/(p4*10^5); // specific Volume at state 4
+imep=WN*10^3/(V4-V2); // Mean effective pressure
+q=Cp*(T3-T2); // Heat supplied
+eff=WN/q; // Efficiency of a Joule cycle
+disp ("bar",imep*10^-5,"Mean effective pressure = ");
+disp ("%",eff*100,"Efficiency of a Joule cycle = ");
+
-- 
cgit