summaryrefslogtreecommitdiff
path: root/632/CH10/EX10.2/example10_2.sce
blob: 5674056ea1838f5c0c0af406a79c826d1d1e9be1 (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
//clc()
Nflue = 100;//kmol
NCO2 = 14.84;
NCO = 1.65;
NO2 = 5.16;
NN2 = 78.35;
PCF = 85;//PERCENT CARBON IN FEED
PIF = 15;//PERCENT INERT IN FEED
//F - amount of coke charged, W - mass of coke left,W = 0.05F
NCflue = NCO2 + NCO ;
MC = 12;
mC = MC * NCflue ;
//carbon balance gives, F * PCF / 100 = W * PCF + mC 
F = mC / ( PCF / 100 - 0.05*PCF / 100);
//let A kmol air supplied, taking N2 balance,
Nair = NN2 * 100/79;
NO2supplied = Nair - NN2;
Ntheoretical = F * PCF / (100 * MC);
Pexcess = ( NO2supplied - Ntheoretical ) * 100 / ( Ntheoretical );
disp("%",Pexcess,"(a)Percentage excess air = ")
mair = Nair * 29;
m = mair / F ;//air supplied per kg of coke charged
disp("kg",m,"(b)air supplied per kg of coke charged = ")
P = 100;//kPa
T = 500;//K
V = Nflue *22.4143*101.325 * T / (F * P * 273.15);
disp("m^3",V,"(c)volume of flue gas per kg of coke = ")
W = 0.05*F;
mCr = W * PCF/100;//carbon in refuse
mir = F * (1-PCF/100);//inert in refuse
mr = mCr + mir;
C = mCr * 100 / mr;
I = mir *100/ mr;
disp("%",C,"(d)Carbon = ")
disp("%",I,"Inert = ")