- •А.В. Кравцов, н.В. Ушева, е.А. Кузьменко,
- •1. Моделирование кинетики химических реакций
- •1.1. Лабораторная работа №1 Моделирование кинетики гомогенных химических реакций
- •Этапы развития химической кинетики
- •Основные понятия химической кинетики
- •Общие представления одношаговых методов решения обыкновенных дифференциальных уравнений
- •Варианты заданий
- •1.2 Лабораторная работа №2 Моделирование кинетики гетерогенных химических реакций
- •Кинетика гетерогенных химических реакций
- •Варианты заданий
- •2. Моделирование структуры потоков в аппаратах
- •2.1. Лабораторная работа №3 Исследование гидродинамики насадочного абсорбера
- •Типовые математические модели структуры потоков в аппаратах
- •М f tатематическое описание гидродинамики насадочного абсорбера
- •Варианты заданий
- •Содержание отчета
- •2.2. Лабораторная работа №4 Моделирование процесса смешения технологических потоков нефтехимического производства
- •Математическое описание процесса смешения
- •Варианты заданий
- •3. Моделирование тепловых процессов химической технологии
- •3.1. Лабораторная работа №5 Моделирование теплообменных аппаратов в стационарном режиме
- •Моделирование и интенсификация работы теплообменной аппаратуры
- •Моделирование теплообменных процессов
- •Варианты заданий
- •Пример результатов расчетов
- •4. Моделирование массообменных процессов химической технологии
- •4.1. Лабораторная работа №6 Исследование процесса разделения многокомпонентной смеси в газовом сепараторе
- •Разделение газожидкостных потоков в химико-технологических процессах
- •Расчет однократного испарения многокомпонентной углеводородной смеси
- •Варианты заданий
- •4.2 Лабораторная работа №7 Математическое моделирование процесса ректификации
- •Описание объекта моделирования
- •Основные уравнения модели
- •Варианты заданий и исходные данные для расчета процесса ректификации
- •5. Моделирование химических реакторов
- •5.1. Лабораторная работа №8 Моделирование гомогенных химических реакторов
- •Классификация реакторов
- •Математическая модель реактора идеального перемешивания
- •Математическая модель реактора идеального вытеснения
- •Исследование химического процесса, протекающего в гомогенном реакторе идеального смешения
- •Исследование химического процесса, протекающего в реакторе идеального вытеснения в стационарном режиме
- •Литература
- •Приложения Приложение а Программы расчета кинетики гомогенных химических реакций
- •Расчет кинетики химических реакций методом Эйлера
- •Расчет температурной зависимости скоростей химических реакций с использованием метода Эйлера
- •Расчет кинетики химических реакций методом Рунге-Кутта
- •Приложение б Программа расчета кинетики гетерогенных химических реакций
- •Приложение в Программа расчета гидродинамики насадочного абсорбера
- •Приложение г Программа расчета смесителя
- •Приложение д Программа расчёта теплообменника
- •Приложение е Программы расчета гомогенных химических реакторов
- •Программа расчёта реактора идеального вытеснения
- •Файл с исходными данными
- •Программа расчёта реактора идеального смешения
- •Файл с исходными данными
- •Приложение ж
- •1. Моделирование кинетики химических реакций 3
- •Математическое моделирование химико-технологических процессов
Приложение г Программа расчета смесителя
Обозначения, принятые в программе:
N – число компонентов в потоке;
C1[1..N] – содержание компонентов в 1-м входном потоке, масс. доли;
C2[1..N] – содержание компонентов во 2-м входном потоке, масс. доли;
С3[1..N] – содержание компонентов в выходном потоке, масс. доли;
G1,G2 – расходы входных потоков, кг/ч;
G – расход выходного потока, кг/ч;
CP1,CP2 – удельные теплоемкости входных потоков;
CP – удельная теплоемкость выходного потока;
a[1..N], b[1..N], c[1..N], d1..N] – коэффициенты для расчета значений теплоемкости компонентов потока в зависимости от температуры;
Т1,Т2 – температуры входных потоков;
Т – значение температуры выходного потока;
Тx – значение температуры выходного потока на предыдущей итерации;
Е – заданная точность вычисления температуры выходного потока.
Program smesit;
var
g,cp,cp1,cp2,e,g1,g2,t1,t2,t,tx:real;
i,k:integer;
c1,c2,c3,c4,a,b,c,d:array[1..4] of real;
f1,f2:text;
begin
assign (f1,'smesit.dta');
assign(f2,'smesit.rez');
reset(f1);rewrite(f2);
for i:=1 to 4 do read(f1,a[i]);
for i:=1 to 4 do read(f1,b[i]);
for i:=1 to 4 do read(f1,c[i]);
for i:=1 to 4 do read(f1,c1[i]);
for i:=1 to 4 do read(f1,c2[i]);
read(f1,e,g1,g2,t1,t2);
g:=g1+g2;
for i:=1 to 4 do c3[i]:=(g1*c1[i]+g2*c2[i])/g;
cp1:=0; cp2:=0;
for i:=1 to 4 do begin
cp1:=cp1+(a[i]+b[i]*t1+c[i]*t1*t1+d[i]*t1*t1*t1)*c1[i];
cp2:=cp2+(a[i]+b[i]*t2+c[i]*t2*t2+d[i]*t2*t2*t2)*c2[i];
end;
t:=(t1+t2)/2;
repeat
tx:=t;cp:=0;
for k:=1 to 4 do
cp:=cp+(a[k]+b[k]*t+c[k]*t*t+ d[k]*t*t*t)*c3[k];
t:=(g1*cp1*t1+g2*cp2*t2)/g/cp;
until abs(tx-t)<e;
writeln(f2,t:8:3,cp:8:3,cp1:8:3,cp2:8:3,g:8:3);
writeln(f2,c3[1]:8:3,c3[2]:8:3,c3[3]:8:3,c3[4]:8:3);
close(f2);
end.
Приложение д Программа расчёта теплообменника
Обозначения, принятые в программе:
d– диаметр внутренней трубы теплообменника;
pi – ;
v1, v2 – объёмные скорости теплоносителя и хладоагента;
ro1, ro2 – плотности теплоносителя и хладоагента;
Cp1, Cp2 – теплоёмкости теплоносителя и хладоагента;
h– шаг по длине теплообменника для вывода результатов;
hl – шаг по длине теплообменника для расчёта температур и вывода соответствующих точек на графике;
K– коэффициент теплопередачи;
L – текущая длина теплообменника;
lt – координата текущей длины теплообменника на графике;
t1, t2 – температуры теплоносителя и хладоагента;
tt1, tt2 – текущие координаты температур теплоносителя и хладоагента;
n– число вычислений между выводом результатов в раздел;
gd, gm – параметры графического режима;
f1 – файловая переменная, соответствующая разделу с результатом.
Program Teplo;
Const d=0.03; v1=2.3e-4; v2=5.1e-4; pi=3.14; ro1=900.0; ro2=1000.0;
cp1=3.35e3; cp2=4.19e3; h=0.2; hl=0.01; k=4900.0;
Var n:integer;
t1,t2,lk,l,b1,b2:real;
f1:text;
Begin
assign(f1,'rteplo.pas');
rewrite(f1);
t1:=200.0; t2:=35.0; lk:=3.0;
b1:=k*pi*d/(ro1*cp1*v1);
b2:=k*pi*d/(ro2*cp2*v2);
l:=0.0;
writeln(f1,'Расчёт процесса теплообмена ');
writeln(f1,' Длина, Т1, Т2,');
writeln(f1,' м град.С град.С');
writeln(f1,l:5:1,t1:10:2,t2:10:2);
n:=0;
while l<=lk do
begin l:=l+hl;
n:=n+1;
t1:=t1-b1*(t1-t2)*hl;
t2:=t2+b2*(t1-t2)*hl;
if n=round(h/hl) then
begin
writeln(f1,l:5:1,t1:10:2,t2:10:2);
n:=0;
end; end;
Close(f1);
end.