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
43
44
45
46
47
48
49
50
51
|
clc;funcprot(0);//EXAMPLE 17.42
// Initialisation of Variables
Cpw=4.18;..............//Specific heat of water in kJ/kgK
n=1;................//No of cylinders
N=350;.......//Engine rpm
pmi=2.8;..........//Mean effective pressure in bar
bl=590;..........//Brake load in N
mf=4.3;............//Fuel consumption in kg
mw=500;..............//Mass of cooling water
tw1=25;...............//Water inlet temperature in C
tw2=50;................//Water outlet temperature in C
ma=33;..................//Mass of air used per kg of fuel in kg
tr=25;.................//Room temperature in C
tg=400;.................//Exhaust temperature in C
D=0.22;.................//Engine bore in m
L=0.28;.................//Engine stroke in m
Db=1;......................//Brake drum diameter in m
C=43900;...................//Calorirfic value of fuel in kJ/kg
Cps=2.09;..................//Specific heat of steamm in exhaust in kJ/kgK
Cpg=1.0;...................//Specific heat of dry exhaust gases in kJ/kgK
k=1;....................//Two stroke engiine
perh=15;...................//Percentage of hydrogen
//Calculations
IP=(n*pmi*N*D*D*L*k*10*(%pi/4))/6;...................//Indicated power in kW
disp(IP,"Indicated power in kW:")
BP=(bl*%pi*Db*N)/(60*1000);......................//Brake power in kW
disp(BP,"Brake power in kW:")
//Heat supplied
hf=(mf/60)*C;................//heat supplied by fuel
hip=IP*60;...........//Heat equivalent of BP in kJ/min
hcw=(mw/60)*Cpw*(tw2-tw1);..........//Heat carried away by cooling water
mg=(mf+(mf*ma))/60;....................//Mass of exhaust gases in kg/min
mst=9*(perh/100)*(mf/60);..................//Mass of steam formed
mdg=mg-mst;..............................//Mass of dry exhaust gases per min
hg=(mdg)*Cpg*(tg-tr);..........//Heat carried by exhaust gasses
hst=mst*(417.5+2257.9+(Cps*(400-99.6)));....................//Heat carried by exhaust steam, the obtained values are from steam tables at NTP
mg=mf+(ma*mf);....................//Mass of exhaust gases in kg/min
ha=round(hf)-round(hip+hg+hst+hcw);............//Unaccounted heat
pf=100;pip=(hip/hf)*100;pcw=(hcw/hf)*100;pg=(hg/hf)*100;pa=(ha/hf)*100;pst=(hst/hf)*100;
printf("\n\n")
printf("HEAT BALANCE TABLE\n")
printf("_______________________________________________________________________\n")
printf("Item kJ Percent\n")
printf("_______________________________________________________________________\n")
printf("Heat supplied by fuel %d %f\n",hf,pf)
printf("Heat absorbed in IP %d %f\n",hip,pip)
printf("Heat taken away by cooling water %d %f\n",hcw,pcw)
printf("Heat carried away by dry exhaust gases %d %f\n",hg,pg)
printf("Heat carried away by steam in exhaust gases %d %f\n",hst,pst)
printf("Unaccounted heat %d %f\n",ha,pa)
printf("_____________________________________________________________________")
|
