MODELOVANJE I SIMULACIJA PROCESA Vežba br. 3 ...

2X (3 − X )= K3(1 − X )873⎛⎜⎝pp0⎞⎟⎠(19)Naravno, prvo je potrebno izračunati konstantu ravnoteže na temperaturi T=873 K (idrugim temperaturama), što se dobija na osnovu zavisnosti K od temperature:d ln KdT0∆Hr= (20)2RT(U jednačini (20) je pretpostavljenjo da jeIz jednačine (20) sledi da je:∆ H = const = ∆ )0rH rKT= K298e∆HR0r⎛⎜⎝12981 ⎞− ⎟T ⎠(21)Konstanta ravnoteže na 298 K se izračunava pomoću jednačine:K298−0∆G rRT0= e(22)Standardna Gibss-ova energija i entalpija hemijske reakcije su, respektivno:0∆G r= ∑υ j∆Gf , j(23)0∆H r= ∑υ j∆Hf , j(24)Pomoću sledećeg programa u MATLABu, primenom naredbe fzero za rešavanjenelinerne jednačine, se dobijaju konverzije za različite radne uslove.% Primer 3.3 - Izracunavanje ravnotezne konverzije reakcije u gasnoj fazi% za razlicite radne uslove T i P% - nelinearna jednacina% N. Nikacevic, 2009clear allclose all% radni uslovi - promenjivo T i P - ciklusT = [873 773 773];P = [1 1 3];% temperature, K% pritisci, atm% pocetna pretpostavka konverzijeX_0 = 0.5;22