Условие задачи:
Газовая смесь, содержащая (% объёмных) хлора и (% объёмных) воздуха, обрабатывается в насадочном абсорбере водным раствором NaOH с целью удаления из смеси 99,5% хлора.
Объёмная скорость газового потока - . Процесс абсорбции проводится при температуре газового потока 17 и атмосферном давлении Р = 1 атм. Диаметр абсорбера составляет . Структуру газового потока в колонне можно интерпретировать моделью идеального вытеснения.
Требуется рассчитать высоту насадочного слоя в колонне, при которой достигается требуемая степень поглощения хлора, т.е. 99,5%, если массовая скорость потока жидкой фазы (водного раствора NaOH), орошающего слой насадки, достаточна для обеспечения рабочего режима процесса хемосорбции. Конечное уравнение представляет собой:
dy/dz=-0.5*S_a/G_vm*sqr(1-y)*ln(1/(1-y))*exp(0.8*ln((71*y/(1-y)+29)*G_vm/S_a));
Где S_a-площадь поперечного сечения абсорбера
G_vm- мольная скорость потока воздуха
Теперь сам готовый код:
Program A;
uses crt;
var
i:integer;
k1,k2,k3,k4,k5,h,z,y,S_a,G_vm:real;
Function f(y:real): real;
begin
G_vm:=0.00042*(100-50)*120;
S_a:=pi*sqr(0.45)/4;
f:=-0.5*S_a/G_vm*sqr(1-y)*ln(1/(1-y))*exp(0.8*ln((71*y/(1-y)+29)*G_vm/S_a));
end;
begin
clrscr;
h:=0.0001;
z:=0; y:=0.5; i:=0;
Repeat
z:=z+h;
i:=i+1;
k1:=h*f(y);
k2:=h*f(y+k1);
k3:=h*f(y+k1/2+k2/2);
k4:=h*f(y+3*k1/8+9*k3/8);
k5:=h*f(y+3*k1/2-9*k3/2+6*k4);
y:=y+1/6*(k1+4*k4+k5);
if ((i=1) or (10204-i<6)
or (i mod 500=0)) then
begin
writeln('z=',z:6:4,'y=',y:6:6);
readln;
end;
until(y<=0.005);
readln;
end.
