summaryrefslogtreecommitdiff
path: root/3821/CH14/EX14.16/Example14_16.sce
blob: 835eefb047c391c5b2ea4e57caf91701ca26cd4b (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
////Chapter No 14 Air Standard Cycles
////Example No 14.16 Page No 317
///Find constant pressure
//Input data
clc;
clear;
r=14;               //compression ratio of standard diesel cycle
P1=1;               //compression stroke in bar
T1=300;             //temperature of air in k
T3=2774;            //temperature rises in k
CP=1.005;
v1=14;
v2=1;
gamma1=1.4;
Qs=1921.43;
R=0.287*10^3;


//Calculation
T2=T1*((v1/v2)^(gamma1-1));      //Constant pressure
rho=T3/T2;                       //cut-off ratio
eta=(1-(1/r^(gamma1-1))*((1/gamma1)*(((rho^(gamma1))-1)/(rho-1))))*100; //air standard efficiency
HS=(CP*(T3-T2));                 //heat supplied
WD=(Qs*eta)*10^-2;                //Net work done
v1=(R*T1/P1)*10^-5;               //characteristics gas equation
v2=(v1/r );                       //characteristics gas equation
Sv=(v1-v2);                       //Swept volume
Pme=(WD/Sv )*10^-2;               //Mean effective pressur
Pm=((P1*r)/((r-1)*(gamma1-1)))*((gamma1*(r^(gamma1-1)))*(rho-1)-((rho^(gamma1))-1));  // mean effective pressure 


//utput
printf('constant pressure= %f K \n',T2);
printf('cut-off ratio= %f \n ',rho);
printf('air standard efficiency= %f percent \n',eta);
printf('heat supplied=  %f KJ/Kg \n',HS);
printf('Net work done= %f KJ/Kg \n',WD);
printf('characteristics gas equation= %f m^3/Kg \n',v1);
printf('characteristics gas equation=%f m^3/Kg \n ',v2);
printf('Swept volume=%f m^3/Kg \n ',Sv);
printf('Mean effective pressure=  %f bar \n',Pme);
printf('Mean effective pressure= %f bar \n ',Pm);