- •Бийский технологический институт (филиал) г.В. Багров, г.И. Севодина, о.Ю. Кравцова общая химическая технология
- •Содержание
- •Введение
- •Цель работы
- •Краткие теоретические сведения
- •Математическое описание процесса
- •Алгоритм решения
- •Задание
- •Контрольные вопросы
- •Лабораторная работа № 2. Исследование процесса восстановления двуокиси углерода углем
- •2.1 Цель работы
- •2.2 Краткие теоретические сведения
- •2.3 Описание кинетики процесса
- •2.4 Задание
- •2.5 Контрольные вопросы
- •Лабораторная работа № 3. Сравнение эффективности работы изотермических реакторов
- •3.1 Цель работы
- •3.2 Краткие теоретические сведения
- •3.3 Математическое описание изотермических реакторов
- •3.4 Задание
- •3.5 Контрольные вопросы
- •Лабораторная работа № 4. Исследование влияния условий проведения газофазной реакции на объем реактора
- •4.1 Цель работы
- •4.2 Краткие теоретические сведения
- •4.3 Математическое описание кинетики процесса
- •4.4 Задание
- •4.5 Контрольные вопросы
- •Лабораторная работа № 5. Расчет неизотермических реакторов
- •5.1 Цель работы
- •5.2 Краткие теоретические сведения
- •5.3 Задание
- •5.4 Контрольные вопросы
- •Лабораторная работа № 6. Реакторы с неидеальной структурой потока
- •6.1 Цель работы
- •6.2 Краткие теоретические сведения
- •6.3 Задание
- •6.4 Контрольные вопросы
- •Лабораторная работа № 7. Обогащение минерального сырья. Флотация
- •7.1 Цель работы
- •7.2 Краткие теоретические сведения
- •7.3 Описание лабораторной установки
- •7.4 Методика проведения работы
- •7.5 Обработка данных
- •7.6 Контрольные вопросы
- •Лабораторная работа № 8. Жесткость воды и еЕ определение. Методы снижения жесткости воды
- •8.1 Цель работы
- •8.2 Краткие теоретические сведения
- •8.3 Описание лабораторной установки
- •8.4 Определение жесткости воды
- •8.5 Обработка данных
- •8.6 Контрольные вопросы
- •Литература
- •А2 Программа № 2
- •А3 Программа № 3
- •А4 Программа № 4
- •А5 Программа № 5
- •А6 Программа № 6
- •А7 Программа № 7
- •А8 Программа № 8
- •А9 Программа № 9
А5 Программа № 5
program lab_3_l;
var a,b,c,d,a1,b1,kl,k2,t,ca0,cb0,xa,xariv,xaris:real;
var tau, dtau, x, y, y1,ca,cb,cr,cs,cq,cd,friv,fris:real;
var i,j,n:integer;
label 1,11,10,100;
function f(x:real):real;
begin
f:=A*exp(3*ln(x))+b*sqr(x)+c*x+d;
end;
begin
tau:=100; t:=tau;
kl:=1.692E-3; k2:=1.708E-3;
ca0:=l; cb0:=l;
writeln ('______________________riv_______________________');
writeln;
dtau:=0.01*tau;
ca:=ca0; cb:=cb0; cr:=0; cq:=0;
writeln('i ca cb cr cq’);
for i:=l to l00 do
begin
ca:=ca-dtau*(kl *ca+k2*ca*cb);
cb:=cb-dtau*k2*ca*cb;
cr:=cr+dtau*k1*ca;
cq:=cq+dtau*k2*ca*cb;
if (i=10)or(i=20)or(i=30)or(i=40)or(i=50) then writeln (i,ca,cb,cr,cq);
if (i=60)or(i=70)or(i=80)or(i=90)or(i= 100) then
writeln(i,ca,cb,cr,cq);
end;
xariv:=(ca0-ca)/ca0;
writeln ('tau=i*dtau;',' dtau=',dtau);
friv:=cr/ca0;
writeln (' xariv=',xariv,' friv=',friv);
writeln ('____________________ris_______________________');
writeln;
a:=sqr(tau)*kl*k2;
b:=t*k2+tau*kl+tau*k2*cb0;
c:=l-tau*k2*ca0;
d:=-ca0;
al:=0.0l; bl:=cb0;
1:x:=(al+bl)/2;
y:=f(x);
if abs(al-bl)<0.001 then goto 10;
yl:=f(al);
if y*yl<0 then bl:=x else al:=x; goto 1;
10:ca:=x; cb:=cb0/(l+tau*k2*ca);
cr:=tau*kl*ca; cs:=cr;
cq:=Tau*k2*ca*cb; cd:=cq;
xaris:=(ca0-ca)/ca0;
fris:=cr/ca0;
writeln('ca=',ca,'cb=',cb,'cr=’,cr,'cq=',cq);
writeln (' xaris=',xaris,' fris=',fris);
xa:=0.1; while xa<xariv do begin tau:=tau+l;
a:=sqr(tau)*kl*k2;
b:=t*k2+tau*k1+tau*k2*cb0;
c:=l-tau*k2*ca0;
d:=-ca0;
al:=0.01; b:=eb0;
ll:x:=(al+bl)/2;
y:=f(x);
if abs(al-bl)<0.001 then goto 100;
yl:=f(al);
if y*yl<0 then bl:=x else al:=x; goto 11;
100:ca;=x; cb:=cb0/(l+tau*k2*ca);
cr:=tau*kl*ca; cs:=cr;
cq:=Tau*k2*ca*cb; cd:=cq;
xa:=(ca0-ca)/ca0;
end;
fris:=cr/ca0;
writeln (‘________________xa riv = xa ris__________________');
writeln;
writeln ('ха=’,хаriv,' tau riv=',t,' friv=',friv);
writeln (‘ tau ris=',tau,' fris=',fris);
end.
А6 Программа № 6
program lab_4;
var bn0,t,P,k,cb0,ca0,e,xp,x,xx,h,s,Vr,r:real;
var i,j:integer;
function f(x:real):real;
begin
f:=exp(3*LN(l+e*x))*exp(3*LN(T))/k/(sqr(l-x)*(Cb0-
0.5*Ca0*x)*exp(3 *LN(P))); end;
begin
bn0:=0.757;
T:=273; P:=5; k:=1.4E+4;
Cb0:=0.08; Ca0:=0.09; e:=-0.45; xp:=0.8; x:=0;
while (x<xp) do begin
x-x+0.1;
{T:=298+(12.15*x-35.25)/(7.005-0.143*x);
k:=1.4E+4*exp(1.7*LN(293/T));}
h:=0.01*x;{integral}
s:=f(0)+f(x);
xx:=0;
for j:=l to 99 do
begin
xx:=xx+h;
s:=s+2*f(xx);
end;
s:=h/2*s; VR:=bn0*s;
writeln ('x=',x,'Vr=',VR);
end;
end.
А7 Программа № 7
program lab_5;
var y,y1,a1,b1,cp,ca,cb,ca0,cb0,cr0,kt,tau,dtau,ro:real;
var t,t0,tt,tt0,ttk,xa,k,k0,ft,e,r,H,V,x,vt,cpt,cpa,cpr:real;
var i,j,m:integer;
label 1,10;
function f(vt:real):real;
begin
t:=t0; tt:=ttk; ca:=ca0; cb:=cb0; dtau:=0.01*tau;
for i:=l to 100 do begin
k:=exp(ln(k0)-E/(8.31*t));
{ writeln(‘k=’,k,’ca=’,ca,’t’t,'tt=',tt);}
r:=k*ca;
ca:=ca-dtau*r;
cb:=cb+dtau*r;
t:=t+dtau*(H*r+FT*kt/V*(tt-t))/(cp*ro);
tt:=tt+dtau*(ft*kt/(tau*vt*cpt)*(tt-t));
end;
f:=tt-284;
end;
begin
t0:=301; ca0:=1.8; cb0:=0; cr0:=10.8;
cpa:=62; cpr:=53; ro:=1000;
k0:=1.84E+5; e:=4.49E+4; h:=-448E+4; tt0:=284;
ttk:=293; kt:=0.320;
tau:=360; V:=0.283; FT:=14.16; cpt:=4.190;
cp:=ca0*cpa+cr0*cpr;
al:=0,l; bl:=4.5;
l:x:=(al+bl)/2;
y:=f(x);
{ writeln (‘x=’,x,'y=’,y);}
if abs(y)<0.1 then goto 10;
yl:=f(al);
if y*yl<0 then bl:=x else al:=x; goto 1;
10:
writeln (‘t0=',t0);
vt:=x;
t:=t0; tt:=ttk; ca:=ca0; cb:=cb0; dtau:=0.01*tau;
for i:=l to 100 do begin
k:=exp(ln(k0)-E/(8.31*t));
r:=k*ca;
ca:=ca-dtau*r;
cb:=cb+dtau*r;
t:=t+dtau*(H*r+FT*WV*(tt-t))/(cp*ro);
tt:=tt+dtau*(ft*kt/(tau*vt*cpt)*(tt-t));
if (i=10)or(i=30)or(i=50)or(i=100) then writeln(i,ca,t,tt);
end;
xa:=(ca0-ca)/ca0;
writeln(‘vt=',vt,'xa=',xa);
end.