blob: 27b9c3e48442f0e7074fca9fb02c299b2e79eb57 (
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
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
|
clear ;
clc;
// Example 13.7
printf('Example 13.7\n\n');
//Page no. 416
// Solution fig E13.7
G = 100 ;// Basis: Pyrolysis Gas-[lb mol]
ub_CO = 10/100 ;// fraction of CO left unburnt
ex_air = 40/100 ;;// fraction of excess air
m_vol = 359 ;// molar volume of gas at std. cond.-[cubic feet]
Ts = 492 ;// Standard temperature -[degree Rankine]
Ps = 29.92 ;//Standard pressure -[in. Hg]
// Analysis of entering gas of entering gas
Tf1 = 90 ;// Temperature of gas-[degree F]
T_gas = Tf1 + 460 ;//Temperature of gas-[degree Rankine]
P_gas = 35 ;//Pressure-[in. Hg]
CO2 = 6.4/100 ;// mol fraction of CO2
O2 = 0.1/100 ;// mol fraction of O2
CO = 39/100 ;// mol fraction of CO
H2 = 51.8/100 ;// mol fraction of H2
CH4 = 0.6/100 ;// mol fraction of CH4
N2 = 2.1/100 ;// mol fraction of N2
// Analysis of entering air
Tf2 = 70 ;// Temperature of air -[degree F]
T_air = Tf2 + 460 ;//Temperature of air-[degree Rankine]
P_air = 29.4 ;//Pressure of air [in. Hg]
f_N2 = 79/100 ;// mol fraction of N2
f_O2 = 21/100 ;// mol fraction of O2
// Get O2 required for combustion of CO,H2 & CH4 according to the following equation
// CO + 1/2O2-->CO2
//H2 + 1/2O2-->H20
//CH4 + 2O2--> CO2 + 2H2O
O2r_O2 = O2 * G ;// O2 required by O2-[lb mol]
O2r_CO = CO * G/2 ;// O2 required by CO-[lb mol]
O2r_H2 = H2 * G/2 ;// O2 required by H2-[lb mol]
O2r_CH4 = G * CH4 * 2 ;// O2 required by CH4-[lb mol]
O2r_total = O2r_O2 + O2r_CO + O2r_H2 + O2r_CH4 ;// Total O2 required-[lb mol]
ex_O2 = ex_air * O2r_total ;// Excess O2-[lb mol]
total_O2 = ex_O2 + O2r_total ;// Total amt of O2 in air-[lb mol]
total_N2 = total_O2 * (f_N2/f_O2);// Total amt of in air-[lb mol]
air = total_O2 + total_N2 ;// Total air entering -[lb mol]
// Product analysis
P_CO = ub_CO * CO * G ;//Unburnt CO in P-[lb mol]
//Element balance of 2N
P_N2 = N2 * G + total_N2 ;// N2 in P-[lb mol]
//Element balance of C
P_CO2 = (CO2 + CO + CH4) * G - 1 * P_CO;//CO2 in P-[lb mol]
// Element balance of 2H
P_H2O = (H2 + 2 * CH4) * G ;// H2 in P-[lb mol]
// Element balance of 2O
P_O2 = (CO2 + O2 + 0.5 * CO) * G + total_O2 -P_CO2-0.5 * (P_H2O + P_CO);// O2 in P-[lb mol]
P = P_CO + P_N2 + P_CO2 + P_H2O + P_O2 ;// Product-[lb mol]
Tf3 = 400 ;// Temperature of product-[degree F]
T_prod = Tf3 + 460 ;//Temperature of product-[degree Rankine]
P_prod = 35 ;// Pressure of product -[in.Hg]
V_gas = (G * m_vol * T_gas * Ps)/(Ts * P_gas);
V_air = (air * m_vol * T_air * Ps)/(Ts * P_air);
V_prod = (P * m_vol * T_prod * Ps)/(Ts * P_prod);
air_ft3 = V_air/V_gas ;//Air supplied per ft^3 of gas entered-[cubic feet]
P_ft3 = V_prod/V_gas ;//Product gas produced per ft^3 of gas entered-[cubic feet]
printf(' Air supplied per ft^3 of gas entered %.2f cubic feet.\n ',air_ft3);
printf(' Product gas produced per ft^3 of gas entered %.2f cubic feet.\n',P_ft3);
|
