CC BY
332
УДК 004.65:541.1
Ю. А. Комиссаров, Дам Куанг Шанг
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРИ РАСЧЕТЕ ПАРОЖИДКОСТНОГО РАВНОВЕСИЯ МНОГОКОМПОНЕНТНЫХ СИСТЕМ
Одной из важнейших задач, которые приходится решать при математическом моделировании любых процессов разделения смесей, в том числе и процессов ректификации, является исследование и математическое описание условий равновесия в парожидкостных смесях. Изучение результатов исследований в этой области позволило нам создать программное обеспечение для расчета и анализа парожидкостного равновесия многокомпонентных систем при изменении входных параметров.
В многокомпонентных системах парожидкостное равновесие для каждого /-го компонента выражается уравнением
У/ = К/Х/, (1)
где К/ - коэффициент фазового равновесия / -го компонента.
В общем случае К/ является функцией, зависящей от давления, температуры и фазового состава системы, т. е.
К/ = I/ робщ, Хь ..., Хс, Уъ ..., Ус),
где Т - температура; Робщ - общее давление системы; х1, ..., Хс, у1, ...,Ус - мольные доли компонентов в жидкой и паровой фазах соответственно; С - число компонентов в системе.
Ниже перечислены наиболее распространенные модели для расчета коэффициента фазового равновесия.
Расчет равновесия пар-жидкость для идеальной смеси (законы Рауля и Дальтона) [1-3]:
К/ = у- • (2)
^общ
где р - давление насыщенного пара / -го компонента при температуре Т .
Модель Чао - Сидера используется для смеси неполярных компонентов (в первую очередь углеводородов) [1, 3, 4]:
К/ =и0ф, (3)
Ф/
где и0 - коэффициент летучести (фугитивности) чистого / -го компонента в жидкой фазе; у, -коэффициент активности / -го компонента в жидкой фазе, определяемый по уравнению (6)
Скэтчарда - Гильдебранда; - коэффициент фугитивности / -го компонента в паровой фазе,
определяемый по уравнению состояния Редлиха - Квонга.
Расчет К1 по уравнению состояния для обеих фаз [1, 3]:
Ф ^
К/-фг • (4)
где Ф^, Фг - коэффициенты фугитивности жидкой и паровой фаз соответственно, рассчитанные по одному из уравнений состояния: Соава - Редлиха - Квонга или Пенга - Робинсона [5].
Эти модели ориентированы на углеводородные смеси и легкие газы (водород, углекислый газ, азот и др.). Они не рекомендованы для неидеальных смесей при низком давлении. Кроме этого, уравнение состояния должно хорошо прогнозировать коэффициент фугитивности как жидкости, так и пара.
Расчет К1 с использованием коэффициента активности [2, 3, 6]:
К/ -У/у-, (5)
общ
где у- - коэффициент активности, определяется по одному из уравнений: Маргулеса, Ван-Лаара, Вильсона и НРТЛ [4, 7].
Эта модель применяется для неидеальных систем при невысоком давлении.
Расчет коэффициента парожидкостного равновесия обычно сводится к вычислению коэффициента активности жидкой фазы. Из существующего множества моделей для расчета коэффициента активности наиболее часто применяются следующие.
Идеальная модель: используется для идеальной системы. При этом коэффициент активности всегда равен единице.
Регулярная модель [1, 2, 7]: коэффициент у^ определяется на основе теории регулярных растворов по уравнению Скэтчарда - Гильдебранда:
1п у - . (6)
' ЯТ
При определении у/ необходимо знать значения двух констант для каждого компонента:
мольного объема жидкости у/1 и параметра растворимости 5^. Величина 5 является усредненным значением параметра растворимости смеси:
_ 2 Х/у^5/
5 —
Величину 5; можно определить по формуле Г ильдебранда:
/ ч1/2
ЛТТ
5/ -
Ш,
V у
где АиI - энергия, необходимая для того, чтобы изотермически перевести мольный объем у//
/-й жидкости из состояния насыщения в состояние идеального газа; у// - мольный объем /-й жидкости, см3/моль.
При значениях температуры значительно ниже критических
А и »АНу/ - ЯТ,
где АНу/ - мольная теплота парообразования чистой /-й жидкости при температуре Т .
Модификацию уравнения Скэтчарда - Гильдебранда, предложенную П. Дж. Флори и М. Хаггинсом, называют расширенным уравнением Скэтчарда - Гильдебранда или уравнением Скэтчарда - Гильдебранда - Флори - Хаггинса:
1пУ/ - (Я-5)2 + 1п0/ +1 -0/, (7)
где 0/ ----/--
Х/у/
2:
Модель Маргулеса [2, 7]:
1п У/ - (а/} + 2(а}/ - А )х/ К2. (8)
Модель Ван-Лаара [2, 7]:
1п у, =-
(9)
1 +
Аух1_ АлХ} ;
В моделях Маргулеса и Ван-Лаара параметры Ау, А]1 рассчитываются по данным для бинарных пар у . Модели Маргулеса и Ван-Лаара непригодны для представления многокомпонентных смесей при наличии только параметров бинарного взаимодействия [7]. Они используются только для бинарных смесей (/ = 1, у = 2 и / = 2, у = 1).
Модель Вильсона [2, 3, 7]:
1п Уг =- 1п
Г N
Е х, Л і ^^'=1
+1-
хк Л кг
N
'Е>
к=^ х, Л
г к
у=1
(10)
где
Л- = -
VI
V1
ехр
Г 1 у -1,у Ї
ят
VI
VI
-ехр
£ у
ят
Л
В модели Вильсона для описания многокомпонентных смесей требуется задать только параметры ^ , определяемые по экспериментальным данным для бинарных пар, входящих
в данную многокомпонентную смесь, и физико-химические свойства чистых компонентов.
Модель двух несмешивающихся жидкостей (НРТЛ) [3, 7]:
N
Ех Ах
1п У, =-
N
- +
N
х,а,
N
Е а
у,ху
Е акухк
N
Е Хк]ак]хк
к_____________
' N
Е акухк
(11)
где
ят
а
у
схр (-ау,ху,);
*й = 0; Ои = 1,
где а у, - определяемые константы, зависящие от природы компонентов. В [7] показано, что значения коэффициента активности незначительно изменяют а^ , которая находится в пределах [-1; 0,3].
Опыт показывает, что равновесие систем пар-жидкость в многокомпонентных смесях может быть рассчитано с достаточной для инженерной практики точностью по уравнениям Вильсона или НРТЛ, но для этого необходимо иметь навык определения параметров бинарного взаимодействия [3]. Модель НРТЛ дает самые лучшие результаты для водных систем [7].
Алгоритм парожидкостного равновесия и решаемые при этом задачи
Задача определения температуры кипения. Часто встречающейся задачей является определение температуры кипения и состава паровой фазы для заданных составов жидкости и общего давления. Из условия парожидкостного равновесия для каждого компонента сумма значений у, должна быть равна единице. Исходя из этого, получаем следующую систему нелинейных уравнений:
2
X
к
к
Еі = уі - Кіхі = 0 і = 1, С;
С
Б¥ = ^ Уі -1 = 0.
і=1
(12)
Задача определения температуры точки росы. В этой задаче заданы состав паровой фазы и общее давление, требуется определить температуру точки росы. Аналогично первой задаче имеется:
Е, = Уі
■ Кіхі = 0
і = 1, ..., С;
БХ = ^ х, -1 = 0.
(13)
і =1
Задача расчета равновесия системы пар-жидкость в общем виде. В общем виде задача расчета парожидкостного равновесия для процесса многокомпонентной ректификации формулируется следующим образом. Задан состав исходной смеси. Требуется определить общие количества пара, жидкости и составы фаз (в молях), находящихся в равновесии при заданных температуре и давлении. При этом имеются следующие уравнения:
— покомпонентный материальный баланс:
V
(14.1)
где zi — мольная доля /-го компонента в исходной смеси;-----отношение пара; V — количество
F
паровой фазы; F — количество исходной смеси;
— уравнения фазового равновесия:
Е = Кх - Уі = 0;
стехиометрическое соотношение:
(14.2)
і=1
Следовательно, решается система уравнений:
7) у'+ - 7)х'- 2' =0'
Кх - уг = 0,
С
2 (Уі - хі )= 0.
і=1
(14.3)
(15)
Так, имеется 2С + 1 уравнений и 2С + 3 (Т , Р, V/F , х1, ..., хс, у1, ..., ус ) переменных. Для решения вышеприведенной системы нелинейных уравнений надо знать две независимые переменные. Обычно выбирают следующие параметры:
— давление Р и температуру Т (РТ);
— давление Р и отношение пара V|F или жидкости V|F (PVF).
Метод решения
Для решения вышеуказанных систем нелинейных уравнений (12), (13) и (15) используется метод Ньютона - Рафсона [2]. В качестве независимых переменных выбирают следующие параметры:
— для расчета температуры кипения: состав паровой фазы и температура (или давление):
X = (хь х2, ..., хс, Т) или X = (хь ^ ..., хс, Робщ ),
вектор функций невязок:
F = (El,E2,...,Ec, SY);
— для расчета температуры точки росы: состав жидкой фазы и температура (или давление):
X = (Уъ У2,..., yc, T) или X = (Уl, У2,..., yc, робщ), вектор функций невязок: F = (El, E2,..., Ec, SX);
— для расчета равновесия системы пар-жидкость в общем виде: состав паровой и жидкой фаз и температуру (или отношение пара):
X = (xl ^...,xc,Уl,У2,...,yc, T) или X = ^x1, x2,..., xc,Уl,У2,...,yc,, вектор функций невязок:
F = (M1,M2, ...,Mc,E1,E2, ...,Ec, SXY).
Программа расчета равновесия системы пар-жидкость
На основе полученных математических моделей мы создали компьютерное обеспечение расчета равновесия CMPR на языке Visual Basic. Эта программа позволяет решить задачи парожидкостного равновесия для многокомпонентных смесей с различными термодинамическими свойствами и входными параметрами. Она состоит из несколько модулей, представленных на рисунке.
Алгоритм расчета программы CMPR
— модуль «Вход-Выход»: позволяет вводить входные параметры и выводить результаты;
— модуль «Операция»: позволяет выбрать тип задачи равновесия системы пар-жидкость и термодинамическую модель;
— модуль «Термодинамическая модель»: рассчитывает коэффициенты фазового равновесия по выбранной термодинамической модели;
— модуль «Математический»: решает систему нелинейных уравнений;
— модуль «База данных»: содержит важную информацию о более чем 600 индивидуальных компонентах пара и жидкости и параметры бинарного взаимодействия более чем 300 пар веществ в уравнениях состояния (Соава - Редлиха - Квонга и Пенга - Робинсона) или уравнениях вычисления коэффициента активности (Маргулеса, Ван-Лаара, Вильсона и НРТЛ).
Для работы с программой CMPR необходимы следующие требования к персональному компьютеру: оперативная память - от 512 МБ, свободное пространство жесткого диска — 30 МБ, операционная система — Microsoft Windows XP и выше.
В табл. 1-5 приведены экспериментальные и расчетные данные равновесия пар-жидкость, полученные с помощью программы CMPR для нескольких систем при разных условиях. Экспериментальные данные парожидкостного равновесия взяты из [8]. В этих таблицах составы жидкой и паровой фаз выражены в мольных долях. Сравнение этих данных показало, что с помощью программы CMPR хорошо прогнозируются составы паровой и жидкой фаз, температура кипения и давление.
Таблица 1
Значения парожидкостного равновесия бинарной системы ацетон - бензол при температуре 318 К (по модели Вильсона)
Номер смеси Состав жидкой фазы x1, % Состав паровой фазы у*экс, % Робщ , бар у*рас, % Робщ, бар Ошибка
Ay —г—, % уГ AP , % Р экс общ
1 4,70 14,44 0,330 15,63 0,335 8,24 1,52
2 40,11 61,39 0,500 61,84 0,510 0,73 2,00
3 95,29 96,99 0,663 97,00 0,668 0,01 0,75
Таблица 2
Значения парожидкостного равновесия бинарной системы пропилен - изобутан при давлении 2,04 МПа (по уравнению Соава - Редлиха - Квонга)
Но- мер смеси Состав жидкой фазы x1, % Состав паровой фазы у*экс, % T экс , K у*рас, % T рас K Ошибка
Ay —г—, % У1экс AT , % T ЭКС
1 20,0 32,1 364,1 30,46 361,9 -5,11 -0,60
2 51,5 67,7 346,4 66,74 344,3 1,42 -0,61
3 84,8 92,7 333,0 91,83 328,9 -0,94 -1,23
Таблица З
Значения парожидкостного равновесия многокомпонентной системы: метиловый спирт (1) - этиловый спирт (2) - вода (3) при давлении 0,1 МПа
(по модели Вильсона)
Номер смеси Номер компо- нента x , моль, % у*экс моль, % Tэкс, K у*р“, % T рас , K Ошибка
Ay . , % у *экс AT -=—, % T ЭКС
1 1 15,66 23,61 348,8 23,74 348,4 0,55 -0,11
2 83,06 75,37 74,92 -0,60
3 1,28 1,02 1,34 31,37
2 1 27,53 37,74 346,92 39,19 346,64 3,84 -0,08
2 71,43 61,36 59,85 -2,46
3 1,04 0,9 0,96 6,67
3 1 49,08 63,60 343,9 62,35 343,6 -1,97 -0,09
2 48,95 34,87 36,23 3,90
3 1,97 1,53 1,42 -7,19
4 1 78,29 85,43 340,1 86,38 339,8 1,11 -0,09
2 19,96 13,69 12,68 -7,38
3 1,75 0,88 0,94 6,82
5 1 69,61 79,71 341,1 80,13 341,0 0,53 -0,03
2 28,17 18,99 18,68 -1,63
3 2,22 1,30 1,18 -9,23
Таблица 4
Значения парожидкостного равновесия многокомпонентной системы: этан (1) - бутан (2) - пентан (3) при температуре 338,6 К (по уравнению Соава - Редлиха - Квонга)
Номер смеси Номер компо- нента x , % y*3“, % Робщ, бар y *рас, % Рорбщ , бар Ошибка
Ay —г—, % y*3KC AP , % Р экс общ
1 1 85,8 0,91,9 59,75 91,42 59,47 -0,52 -0,47
2 3,0 1,7 2,11 24,12
3 11,2 6,4 6,47 1,09
2 1 70,3 91,7 48,38 91,75 47,24 0,05 -2,36
2 3,9 1,8 1,69 -6,11
3 25,8 6,5 6,56 0,92
3 1 74,0 87,1 52,3 86,97 51,37 -0,15 -1,78
2 20,8 11,0 11,13 1,18
3 5,2 1,9 1,90 0,00
4 1 51,0 82,8 35,01 81,94 35,24 -1,04 0,66
2 37,8 14,9 15,65 5,03
3 11,2 2,3 2,41 4,78
5 1 71,3 88,6 50,18 87,99 49,04 -0,69 -2,27
2 17,5 7,7 8,51 10,52
3 11,2 3,7 3,49 -5,68
Таблица З
Значения парожидкостного равновесия многокомпонентной системы: метан (1) - этан (2) - пропан (3) - бутан (4) - пентан (5) при температуре 310,8 К (по уравнению Свона - Редлиха - Квонга)
Номер смеси Номер компо- нента x , % у*экс, % Рбщ , бар y*-“, % Ррщ ,бар Ошибка
Ay y *экс AP , % Р экс общ
1 1 41,99 84,42 97,3 85,17 94,6 0,89 -2,77
2 7,52 6,56 6,05 -7,77
3 6,81 2,85 2,81 -1,40
4 10,55 2,07 2,24 8,21
5 33,13 4,10 3,73 -9,02
2 1 38,38 85,61 88,2 85,29 86,38 -0,37 -2,06
2 7,56 6,15 6,2 0,81
3 7,05 2,77 2,81 1,44
4 11,29 2,15 2,2 2,33
5 35,72 3,32 3,5 5,42
Выводы
В результате исследований нами были созданы математические модели и программное обеспечение парожидкостного равновесия многокомпонентных систем, что позволяет анализировать и прогнозировать парожидкостное равновесие смесей при изменении входных параметров. База данных, соединенная с CMPR, содержит 27 термодинамических свойств более 600 индивидуальных компонентов пара и жидкости и параметры бинарного взаимодействия более 300 пар веществ в уравнениях состояния и уравнениях вычисления коэффициента активности. Программа позволяет использовать ее для моделирования и расчета химико-технологических процессов (ректификация, абсорбция и др.) с целью анализа и синтеза сложных систем разделения.
СПИСОК ЛИТЕРАТУРЫ
1. Холланд Ч. Д. Многокомпонентная ректификация. - М.: Химия, 1969. - 351 с.
2. Комиссаров Ю. А., Гордеев Л. С., Вент Д. П. Научные основы процессов ректификации. - М.: Химия, 2004. - Т. 1. - 270 с.
3. Комиссаров Ю. А., Гордеев Л. С., Вент Д. П. Научные основы процессов ректификации. - М.: Химия, 2004. - Т. 2. - 416 с.
4. Уэйлес С. Фазовые равновесия в химической технологии. - М.: МИР, 1989. - Т. 2. - 360 с.
5. Reid R. C., Prausnitz J. M., Poling B. E. The Properties of Gases and Liquids. - N. Y.: McGraw-Hill, Inc., 1987. - 741 p.
6. Кафаров В. В., Ветохин В. Н., Бояринов А. И. Программирование и вычислительные методы в химии и химической технологии. - М.: Наука, 1972. - 486 с.
7. Уэйлес С. Фазовые равновесия в химической технологии. - М.: МИР, 1989. - Т. 1. - 304 с.
8. Коган В. Б., Фридман. В. М., Кафаров В. В. Равновесие между жидкостью и паром. - М.; Л.: Наука, 1966. - Т. 2. - 795 с.
Статья поступила в редакцию 24.12.2010
MATHEMATICAL MODELING AT THE CALCULATION OF THE VAPOR-LIQUID EQUILIBRIUM OF MULTICOMPONENT SYSTEMS
Yu. A. Komissarov, Dam Quang Sang
A mathematical description of equilibrium conditions in vapor-liquid mixtures with the use of the equations of Wilson, NRTL, Soave - Redlich - Kwong, Peng - Robinson and others is given in the paper. Software is developed in order to calculate and analyze the vapor-liquid equilibrium of multicomponent systems when input parameters (temperature, pressure, composition of the liquid and vapor phases) are changed. A database is created that contains important thermodynamic properties of more than 600 components for gas and liquid, parameters of the binary interaction (more than 300 pairs) in the above mentioned equations. This database was used in calculating the vapor-liquid equilibrium of multicomponent systems.
Key word: mathematical modeling, vapor-liquid equilibrium, multicomponent mixture, Newton - Raphson method, software, thermodynamic properties, database.
