blob: 39c871a649207d797443c1f2ad1cad52f179e190 (
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
|
clc
clear
//Input data
n=6;//Number of cylinder
bp=130;//Brake power in kW
N=1800;//The speed of the engine in rpm
CV=42000;//The calorific value of the fuel in kJ/kg
C=86;//The composition of carbon in the fuel in percent
H=13;//The composition of Hydrogen in the fuel in percent
NC=1;//The non combustibles present in the fuel in percent
na=85;//The absolute volumetric efficiency in percent
ni=38;//The indicated thermal efficiency in percent
nm=80;//The mechanical efficiency in percent
ac=110;//The excess consumption of air in percent
sb=1.2;//The stroke to the bore ratio
da=1.3;//The density of air in kg/m^3
pi=3.141;//Mathematical constant of pi
//Calculations
saf=(((C/100)*(32/12))+((H/100)*(8/1)))*(1/0.23);//The stoichiometric air fuel ratio
aaf=saf*(1+1.1);//The actual air fuel ratio
Ma=(0.23*32)+(0.77*28);//The molecular weight of air in kg/kmol
a=(C/100)/12;//For carbon balance
b=(H/100)/2;//For hydrogen balace
x=aaf/Ma;//Number of kmol of air per kg of fuel
c=(0.21*x)-a-(b/2);//For oxygen balance
d1=0.79*x;//For nitrogen balance
ip=bp/(nm/100);//The indicated power in kW
mf=ip/[(ni/100)*CV];//The mass flow rate of fuel in kg/s
ma=mf*aaf;//The mass flow rate of air in kg/s
Va=ma/da;//Actual volume flow rate in m^3/s
Vs=Va/(na/100);//The swept volume per second in m^3/s
d=[[Vs*(4/pi)*(1/1.2)*((2*60)/N)*(1/n)]^(1/3)]*1000;//The diameter of the bore in mm
L=1.2*d;//The length of the stroke in mm
T=a+c+d1;//The total composition in kmol
CO2=(a/T)*100;//The volume of CO2 in %
O2=(c/T)*100;//The volume of O2 in %
N2=(d1/T)*100;//The volume of N2 in %
//Output
printf(' The volumetric composition of dry exhaust gas : \n 1) CO2 = %3.5f kmol and volume = %3.2f percent \n 2) O2 = %3.5f kmol and volume = %3.2f percent \n 3) N2 = %3.5f kmol and volume = %3.2f percent \n The bore of the engine = %3.0f mm \n The stroke of the engine = %3.1f mm ',a,CO2,c,O2,d1,N2,d,L)
|
