Расчет и оптимизация процессов и аппаратов химических и пищевых производств в среде MatLab. Дворецкий Д.С - 35 стр.

UptoLike

ВУЗ: 

ТГТУ | Тамбов

Составители: 

Рубрика: 

Для вычисления летучести как функции давления при заданной температуре используется зависи-
мость:
∫
−
=
P
dP
P
Z
P
f
0
1
ln
,
где Z – функция сжимаемости газа; P – давление; f – летучесть аммиака.
Вычислить летучесть и объем газа при заданных условиях.
9 Зависимость сжимаемости аммиака от давления
№
опыта
Давле-
ние
Функция
сжимае-
мости газа
Z
№ опы-
та
Давле-
ние
Функция
сжимае-
мости газа
Z
1 1 0,998 9 100 0,801
2 10 0,981 10 200 0,551
3 20 0,961 11 300 0,462
4 30 0,942 12 400 0,495
5 40 0,922 13 500 0,557
6 50 0,902 14 600 0,621
7 60 0,882 15 800 0,755
8 80 0,841
В функции ammiak1 записывается подынтегральная функция
function f = ammiak1(p, z)
% расчет зависимости летучести аммиака
% z – независимый параметр (летучесть)
f = (z-1)./p;
Процедура ammiak рассчитывает летучесть и объем газа, результаты расчета выводятся на графики
в одно окно.
% программа ammiak расчета летучести аммиака
% в зависимости от давления
% с использованием оператора for
R = 8.31514; % универсальная газовая постоянная
T = 473; % температура, К
% табличные данные
P = [0.001 1 10 20 30 40 50 60 80 100 200 300 ...
400 500 600 800];
Z = [0 0.998 0.981 0.961 0.942 0.922 0.902 0.882...
0.841 0.801 0.551 0.462 0.495 0.557 0.621 0.755];
% Расчет
for i=1:length(P)
I = quad('ammiak1', 1.0e-10, P(i), [], [], Z(i));
f(i) = P(i)*exp(I); % летучесть аммиака
v(i) = (Z(i)*R*T)/P(i); % объем газа
end
% графический вывод результатов
subplot(3,1,1)
plot ( P, Z, 'k-')

Страницы