diff options
Diffstat (limited to '2063/CH1/EX1.24/1_24.sce')
-rwxr-xr-x | 2063/CH1/EX1.24/1_24.sce | 25 |
1 files changed, 25 insertions, 0 deletions
diff --git a/2063/CH1/EX1.24/1_24.sce b/2063/CH1/EX1.24/1_24.sce new file mode 100755 index 000000000..ab67c95da --- /dev/null +++ b/2063/CH1/EX1.24/1_24.sce @@ -0,0 +1,25 @@ +clc
+clear
+//Input data
+P1=1;//pressure at the beginning of compression stroke in bar
+T1=298;//Temperature at the beginning of compression stroke in K
+P3=38;//Pressure at the end of constant volume heat addition in bar
+T4=1573;//Temperature at the end of constant volume heat addition in K
+r=9.5;//Compression ratio
+Cp=1.004;//Specific heat of air at constant pressure
+Cv=0.717;//Specific heat of air at constant volume
+
+//Calculations
+r1=Cp/Cv;//Isentropic ratio
+T2=T1*r^(r1-1);//Temperature at the end of compression stroke in K
+P2=P1*r^r1;//Pressure at the end of compression stroke in bar
+T3=T2*(P3/P2);//Temperature at the end of constant volume heat addition in K
+rho=T4/T3;//Cut off ratio
+T5=T4*(rho/r)^(r1-1);//Temperature at the end of expansion stroke in K
+Qs=Cv*(T3-T2)+Cp*(T4-T3);//Heat supplied per kg in kJ
+Qr=Cv*(T5-T1);//Heat rejected per kg in kJ
+W=Qs-Qr;//Work done per kg of air in kJ
+n=(W/Qs)*100;//Efficiency of the air standard dual cycle in percent
+
+//Output
+printf('(a)The work done per kg of air is %3.1f kJ\n (b)Cycle efficiency is %3.2f percent',W,n)
|