ВУЗ:
Составители:
42
R1out:=K1(Tx)*(P*(H2-2*Xt-Yt)*exp(0.25*Ln((CO-Xt+Yt)/(Met+Xt)))-exp(0.25*Ln((Met+Xt)/(CO-
Xt+Yt)))/exp(0.5*Ln(Kp1(Tx))));
R2out:=K2(Tx)*P*((CO2-Yt)*(H2-2*Xt-Yt)-(CO-Xt+Yt)*(H2O+Yt)/Kp2(Tx))/(0.8*(H2O+Yt)+(CO2-Yt));
IF n = 1 THEN
begin
tau1:=60*1000*dX/R1out;
tau2:=60*1000*ABS(Yt-Y0)/R2out;
end
ELSE
begin
tau1:=60*1000*dX/((R1in+R1out)/2);
tau2:=60*1000*ABS(Yt-Y0)/((R2in+R2out)/2);
end;
IF ABS(tau1-tau2)>epsilon THEN GoldZero(0, Xt, 0.000001, xG, Fg);
tau0:=tau0+tau1; {Время контакта}
tau[i,n]:=tau0;
W[i,n]:=U*tau0/3600; {Объем катализатора}
X[i,n]:=Xt; {текущая конверсия СО (выход метанола в 1-й реакции)}
Y[i,n]:=Yt; {текущая конверсия СО
2
(выход воды во 2-й реакции)}
T[i,n]:=Tx; {температура в слоеdx}
Y0:=Yt;
if (U*(Met+Xt)*Ro/F)>=Gm then Exit;
{Проверка первого условия}
until Tx>=Tout;
end;
begin
ClrScr;
{Подготовка файла вывода}
Writeln('Введите имя файла вывода');
Readln(FileName);
Assign(Dev, FileName);
ReWrite(Dev);
{Задание исходных значений}
epsilon:=0.000001;
dX:=0.0001;
Xt:=0.0;
Yt:=0.0;
n:=0;
CO:=CO0;
CO2:=CO20;
H2:=H20;
Inr:=Inr0;
U:=282000.0;
Tx:=Tn;
Nitro PDF Trial
www.nitropdf.com
Страницы
- « первая
- ‹ предыдущая
- …
- 41
- 42
- 43
- 44
- 45
- …
- следующая ›
- последняя »
