Программа расчета кинетики гетерогенных химических реакций

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

Ca0,Cb0,Cc0,Cd0,Cz0,Czh0,Czc0 – начальные значения концентраций;

C – массив текущих значений концентраций;

f – массив правых частей уравнений математической модели;

k1,k2,k3,k4,k5 – константы скорости химических реакций;

w [i],i=1,5 – скорости стадий гетерогенной химической реакции;

tk – конечное время (время протекания реакций), с;

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

h – шаг по времени, с;

xa – степень превращения водорода;

i,j – переменные цикла;

n – счётчик шагов;

r – имя файловой переменной.

Program get;

Type

Mas1 = array [1..7] of real;

Mas2 = array [1..5] of real;

Var

f,c: mas1; w: mas2;

n, i, j: integer;

r: text;

V, t, xa: real;

Const

K1=0.47; k2=0.33; k3=0.21; k4=0.035; k5=0.11; h=0.5;

tk=60; Ca0=0.61; Cb0=0.47; Cco=0.0; Cd0=0.0; Cz0=1.0;

Czh0=0.0; Czc0=0.0;

Procedure PR (k1,k4,k2, k3, k5: real; w: mas2; var f: mas1;C: mas1);

Begin

W[1]:= k1*C[1]*C[5];

W[2]:= k2*C[2]*C[6];

W[3]:= k3*C[7];

W[4]:= k4*C[6];

W[5]:= k5*C[3]*C[4]*C[5];

f[1]:= -w[1]+w[4];

f[2]:= -w[2];

f[3]:= w[3]-w[5];

f[4]:= w[3]-w[5];

f[5]:= -w[1]+w[4]+w[3]-w[5];

f[6]:= w[1]-w[4]-w[2];

f[7]:= w[2]-w[3]+w[5];

BEGIN

assign (r, ‘rget.pas’);

rewrite(r);

writeln (r);

writeln (r, ‘ Таблица 1 – Концентрации реагирующих веществ’);

writeln (r, ‘ Время C_H2 C_C7H8 C_CH4 C_C6H6 C_Z C_ZH2 C_{ZC7H8 * H2}’);

C[1]:= Ca0; C[2]:=Cb0; C[3]:=Cc0; C[4]:=Cd0; C[5]:=Cz0; C[6]:=Czh0;

C[7]:=Czc0; t:=0.0; n:=0;

Writeln (r, t:5:2, C[1]:6:3, C[2]:6:3, C[3]:6:3, C[4]:6:3, C[5]:6:3, C[6]:6:3,

C[7]:6:3);

Repeat

t:=t+h;

n:=n+1;

pr (k1, k4, k2, k3, k5, w, f, C);

for i:=1 to 7 do

C[ i ]:=C[ i ]+h * f [ i ];

If n=10 then begin

Writeln (r, t:5:2, C[1]:6:3, C[2]:6:3, C[3]:6:3, C[4]:6:3, C[5]:6:3, C[6]:6:3,

C[7]:6:3);

n:=0; end;

until t > tk;

xa:= (Ca –C[ 1 ]) / Ca0;

writeln (r, ‘Степень превращения водорода Х_H2:= ‘, xa:5:3);

close (r);

End.

Приложение В