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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

40
{Вычисление скорости реакции на выходе из слоя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));
{Вычисление времени контакта в слоеdx}
IF n = 1 THEN
begin
tau1:=60*1000*dX/R1out;
FGold := R2out * tau1 / (60*1000) + Y0 - Yt;
tau2:=60*1000*ABS(Yt-Y0)/R2out;
end
ELSE
begin
tau1:=60*1000*dX/((R1in+R1out)/2);
FGold := (R2in + R2out) / 2 * tau1 /(60 * 1000) + Y0 - Yt;
tau2:=60*1000*ABS(Yt-Y0)/((R2in+R2out)/2);
end;
end;
Procedure GoldZero (A, B, E:real; var x,F:real);
{"Метод золотого сечения". ПроцедураGoldZero производит поиск интервала[X1,X2],
в котором лежит решение уравненияF(x)=0. F(x) определяется в процедуре Fх.
А,В - интервал возможного варьирования переменной х; А - нижний предел; В -
верхний предел; Е - относительная погрешность вычислений (относительно интер-
вала А - В) так, что (Х2-Х1)<=Е*(В-А)}
begin
t1:=0.3819660113;
t2:=1 - t1;
x0:=A;
x1:=A + t1*(B - A);
x2:=A + t2*(B - A);
x3:=B;
Fx(F0,x0);
Fx(F1,x1);
Fx(F2,x2);
Fx(F3,x3);
i:=x3 - x1;
cicle:=0;
While (i>E*(B-A)) do
begin
cicle:=cicle+1;
if ((F0 < 0) and (F3 < 0)) or ((F0 > 0) and (F3 > 0)) then
begin
if cicle>7 then
begin
writeln('Я Вас предупреждала уже семь раз');
writeln('Сейчас я говорю: "Это уравнение не имеет решений"');
Nitro PDF Trial
www.nitropdf.com

Страницы