From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 2063/CH1/EX1.24/1_24.sce | 25 +++++++++++++++++++++++++ 1 file changed, 25 insertions(+) create mode 100755 2063/CH1/EX1.24/1_24.sce (limited to '2063/CH1/EX1.24/1_24.sce') 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) -- cgit