blob: 2d4f8566e467ca84a143f9776f5057f10bef0b26 (
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
|
clc;funcprot(0)//EXAMPLE 7.4
//Initializing the variables
n=1;...................//No of cylinders
k=0.5;
Vs=7000;............//displacement volume in cm^3
Vssi= Vs*(10^(-6));........//displacement volume in m^3
OP=14.7;...................//Power developed in kW
N=450;..................//Engine rpm
sfc=0.272;................//Specific fuel consumption in kg/kWh
//Fuel used is C7H16
mC=12;.............//mass of carbon in amu
mH=1;.................//mass of hydrogen in amu
mO=16;.................//mass of oxygen in amu
pi=1.013 * (10^5);................//initial pressure in pascal
T=30;...................//initial temperature in degree celsius
Tk=30+273;................//initial temperature in degree kelvin
R=287;..................//Gas constant for air in J/kg.K
//calculations
disp("The given fuel C7H16 requires 11(O2) for complete combustion")
mf=(7*mC)+(16*mH);
disp (mf,"Hence, Mass of fuel is :")
MO=11* 2 * mO;
disp (MO,"Mass of Oxygen required is :")
ma = MO/0.23;.......//mass of air
//Air contains 23% of oxygen by weight
afr = ma/mf;...............//air fuel ratio is the ratio of mass of air to mass of fuel
disp(afr,"The air fuel ratio is :")
MF = sfc * OP;...........//actual fuel consumed in kg/h
MA = afr*MF;
AAS = MA * (1+0.3);....................//actual air supplied in kg/h
M = AAS + MF;................//mass of charge in kg/h
VCS = ((M/60)*R*Tk)/pi;.............//Volume of charge sucked in m^3/min
DVM = Vssi * (N/2);..............//Displacement volume/min
eta = VCS/DVM;
disp (eta*100,"The volumetric efficiency of the engine is (%) :")
|
