Расчеты реакционной аппаратуры химических производств. Иванов Г.Н - 38 стр.

UptoLike

37
repeat
i:=0;
n:=n+1 {Номер полки};
F:=1-2*Xt {Подсчет количества молей};
U1:=U*F*(Tx-Tin)/(Tin-Tp) {Подсчет количества поддува};
P1:=U1/U;
{Расчет мольного содержания компонентов реакционной массы}
CO2:=((CO2-Yt)/F+P1*CO20)/(1+P1);
CO:=((CO-Xt+Yt)/F+P1*CO0)/(1+P1);
H2:=((H2-2*Xt-Yt)/F+P1*H20)/(1+P1);
Inr:=(Inr/F+P1*Inr0)/(1+P1);
H2O:=((H2O+Yt)/F+P1*H2O0)/(1+P1);
Met:=((Met+Xt)/F+P1*Met0)/(1+P1);
U:=U+U1 {Расчет объемной скорости};
If n<>1 Then {Если полка не первая}
Begin {Расчет скорости процесса на входе в полку n}
R1in:=K1(Tin)*(P*H2*exp(0.25*Ln(CO/Met))-exp(0.25*Ln(Met/CO))/exp(0.5*Ln(Kp1(Tin))));
R2in:=K2(Tin)*P*(CO2*H2-CO*H2O/Kp2(Tin))/(0.8*H2O+CO2);
end;
{Вычисление параметров реакции в слое dx }
repeat
i:=i+1;
Xt:=Xt+dX;
Yt:=Yt+dX;
F:=1-2*Xt;
Label1:
{Вычисление текущей температуры в слоеdx}
Tx:=Tin+2835*Xt+1296*Yt;
{Вычисление скорости реакции на выходе из слоя dx}
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;
{Подбор значения шага по Y}
Nitro PDF Trial
www.nitropdf.com