Оптимизация химико-технологических процессов органического синтеза. Бочкарев В.В - 47 стр.

UptoLike

ВУЗ: 

ТПУ | Томск

Составители: 

Рубрика: 

47
Расчет констант равновесия можно проводить по приближенным
уравнениям, справедливым для интервала температур 800 - 1200 К:
lg K = 7.54 - 7120 / T
p,1
0
,
lg K = - 5.60 + 4600 / T
p,2
0
,
lg K = 9.28 - 9000 / T
p,3
0
.
Решение системы трех нелинейных уравнений можно провести каким-
либо из описанных в литературе способами [10]. Наиболее часто для этой
цели используют метод Ньютона-Рафсона. Ниже приведена программа
решения данного примера.
Алгоритм работы программы:
Вводятся исходные данные – температура, давление, требуемое
мольное соотношение H
2
/CO на выходе из реактора.
Проводится расчет констант равновесия.
При определенном соотношении СН
4
/Н
2
О в исходной смеси
проводится решение системы нелинейных уравнений (процедура
noLineUraw).
Рассчитывается мольное соотношение H
2
/CO.
Определение мольного соотношении СН
4
/Н
2
О в исходной смеси
проводится методом “золотого сечения” (процедура GoldMin) путем
минимизации функции (X
H2/CO,треб.
- X
H2/CO,расч.
)
2
.
program OptRatio;
Uses Crt;
type massiv=array[1..20] of real;
massiv2=array[1..20,1..2] of real;
matrix=array[1..20,1..20] of real;
var n : integer;
f,x : massiv;
intX : massiv2;
KP1,KP2,KP3,p,beta,alfamax,e,RatioOpt,Ratio : real;
xco,xh2,xch4,xco2,xh2o,sx,T,t0,alfa,Fgold : real;
{***********************************************************}
procedure systemUraw;
begin
xco:=x[1]+2*x[3]; xh2:=x[1]-2*x[2]; xch4:=x[2]+alfa;
xco2:=beta-x[3]; xh2o:=1-x[1]; sx:=1+beta+alfa+x[1]-x[2]+x[3];
f[1]:=xco*xh2*p-KP1*(xh2o*sx);
f[2]:=xch4*sx-KP2*(xh2*xh2*p);
f[3]:=xco*xco*p-KP3*(xco2*sx);
end;
{***********************************************************}
Nitro PDF Trial
www.nitropdf.com

Страницы