summaryrefslogtreecommitdiff
path: root/632/CH12/EX12.15/example12_15.sce
blob: 6e023294a5bc8f79c81223a0155289f5930a8ea4 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
//clc()
//CO(g) + H2O(g) = CO2(g) + H2(g)            H298 = -41.190
T1 = 298;//K
Pconv = 75;//%
T2 = 800;//K
H298 = -41.190;//kJ
Hco = 30.35;//J/molK
Hco2 = 45.64;//J/molK
Hwater = 36;//J/molK
Hh2 = 29.3;//J/molK
Nco = 1;//mol
Nh2o = 1;//mol
Ncofinal = Nco * (1 - Pconv/100);
Nwaterf = Ncofinal;
Nco2final = Nco - Ncofinal;
Nh2final = Nco2final;
H2 = (Nco2final * Hco2 + Nh2final * Hh2 + Ncofinal * Hco + Nwaterf * Hwater) * (T2 - T1);
Hr1 = H298 * (Nco - Ncofinal);
Hr2 = Hr1 * 1000 + H2;
mh2 = Nh2final * 2.016 * 10^-3;//kg
//therfore for 1000k H2,
Hr = Hr2 * 1000 / (mh2 * 1000);//kJ
disp("kJ",Hr,"Amount of heat change for 1000kg of hydrogen produced = ")