summaryrefslogtreecommitdiff
path: root/551/CH13/EX13.22/22.sce
diff options
context:
space:
mode:
Diffstat (limited to '551/CH13/EX13.22/22.sce')
-rwxr-xr-x551/CH13/EX13.22/22.sce46
1 files changed, 46 insertions, 0 deletions
diff --git a/551/CH13/EX13.22/22.sce b/551/CH13/EX13.22/22.sce
new file mode 100755
index 000000000..40d145d54
--- /dev/null
+++ b/551/CH13/EX13.22/22.sce
@@ -0,0 +1,46 @@
+clc
+r1=15.3; //V1/V2
+r2=7.5; //V4/V3
+p1=1; //bar
+T1=300; //K
+n_mech=0.8;
+C=42000; //kJ/kg
+y=1.4;
+R=287;
+cp=1.005;
+cv=0.718;
+V2=1; ////Assuming V2=1 m^3
+
+T2=T1*r1^(y-1);
+p2=p1*r1^y;
+T3=r1/r2*T2;
+m=p2*10^5*V2/R/T2;
+T4=T3/r2^(y-1);
+
+Q_added=m*cp*(T3-T2);
+Q_rejected=m*cv*(T4-T1);
+W=Q_added-Q_rejected;
+
+pm=W/(r1-1)/V2/100;
+disp("Mean effective pressure =")
+disp(pm)
+disp("bar")
+
+ratio=p2/pm;
+disp("Ratio of maximum pressure to mean effective pressure =")
+disp(ratio)
+
+n_cycle=W/Q_added;
+disp("Cycle efficiency =")
+disp(n_cycle)
+
+n_thI=0.5;
+n_cycle1=n_thI*n_cycle;
+
+n_thB=n_mech*n_cycle1;
+
+BP=1;
+mf=BP/C/n_thB*3600;
+disp("Fuel consumption per kWh =")
+disp(mf)
+disp("kg/kWh") \ No newline at end of file