Программа расчета гидродинамики насадочного абсорбера

Обозначения, принятые в программе:

Ha – высота абсорбера, м;

S – площадь поперечного сечения абсорбера, м²;

V – объёмный расход исходной газовой смеси, м³/ч;

C0 – концентрация извлекаемого компонента в исходной газовой смеси, моль/л;

N – число ячеек в модели;

tau – время пребывания, с;

С – массив концентраций внутри ячеек;

Сe – массив экспериментальных значений концентраций на кривой отклика, моль/л;

Сr – массив расчетных значений концентраций на кривой отклика, моль/л;

tk – время переходного процесса, с;

t – текущее время, c;

h – шаг интегрирования;

F – массив сумм квадратов отклонений расчетных значений концентраций на кривой отклика от экспериментальных;

Fmin – минимальное значение суммы квадратов отклонений расчетных значений концентраций на кривой отклика от экспериментальных;

N0 – оптимальное число ячеек.

Program absor;

Type mas=array[1..10] of real;

const Ha=11.5;S=1.8;V=10000.0;C0=0.18;h=1;

Ce:mas=(0.002,0.05,0.1,0.12,0.14,0.15,0.155,0.159,0.161,0.166);

var Cr,F:mas;

C:array[0..10]of real;

tau,Fmin:real;

i,N,N0,t:integer;

f1:text;

begin

assign(f1,'rabsor.pas');

rewrite(f1);

tau:=Ha*S/V*3600;

Writeln(f1,'Расчет гидродинамики абсорбера');

for N:=1 to 10 do

begin

C[0]:=C0;

for i:=1 to N do

C[i]:=0;

F[N]:=0;

for t:=1 to 10 do

begin

for i:=1 to n do

C[i]:=C[i]+h*N/tau*(C[i-1]-C[i]);

Cr[t]:=C[N];

F[N]:=F[N]+sqr(Ce[t]-Cr[t]);

end;

end;

Fmin:=F[1];

N0:=1;

For i:=2 to 10 do

if F[i]<Fmin then

begin

Fmin:=F[i];

N0:=i;

end;

C[0]:=C0;

for i:=1 to N0 do

C[i]:=0;

writeln(f1,'Оптимальное число ячеек N=',N0, '.. Fmin=', Fmin:11);

Writeln(f1,'Время,с Сэкс Срасч');

for t:=1 to 10 do

begin

for i:=1 to n0 do

C[i]:=C[i]+h*N0/tau*(C[i-1]-C[i]);

Cr[t]:=C[N0];

writeln(f1,t:4,Ce[t]:9:3,Cr[t]:9:3);

end;

close(f1);

end.