CC BY
281
УДК 519.7:541.12
Ю. А. Комиссаров, Дам Куанг Шанг
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ И АЛГОРИТМ РАСЧЕТА ПРОЦЕССА РЕКТИФИКАЦИИ МНОГОКОМПОНЕНТНЫХ СИСТЕМ
Ключевые слова: математическая модель, алгоритм, ректификация, многокомпонентные системы, точка ки-
На основе метода независимого определения концентраций использованием трехдиагональной матрицы [1, 2, 3], авторы создали эффективный алгоритм для моделирования и расчета процесса ректификации многокомпонентных смесей. Алгоритм позволяет рассчитать процесс ректификации с питаниями в виде кипящей жидкости, пара или обеих этих состояний. В этой работе коэффициент равновесия системы пар-жидкость и энтальпия фаз рассчитаны с помощью уравнения состояния Соава-Редлиха-Квонга (СРК).
Key words: mathematical model, algorithm, distillation, multi-component systems, bubble point, the equation of state.
On the basis of the independent determination of the concentration method by using tridiagonal matrix [1, 2, 3], the authors have created a efficient algorithm for modeling and calculation of
the distillation of multi-component systems. Algorithm allows to calculate process of distillation with the feeds in the form of boiling liquid, vapor or both of these states. In this paper, the equilibrium ratio of the vapor-liquid phase and the enthalpy ofphases are calculated by using the equation of state Soave-Redlich-Kwong (SRK).
1. Математическое описание процесса ректификации
Рассмотрим процесс ректификации, схема потоков которого приведена на рис. 1. Ректификационная колонна содержит N тарелок, смесь состоит из C компонентов. Тогда математическая модель в соответствии со схемой рис.1 для этой системы будет иметь вид:
- уравнение покомпонентного материального баланса
пии H , количества F , составы z жидкой и паровой фаз питания; 1 < i < С; 1 < j < N;
пения, уравнение состояния.
(1)
уравнение фазового равновесия
(2)
- уравнения стехиометрических соотношении
для жидкой фазы:
С
(3)
для паровой фазы:
С
соответственно энталь-
(4)
(5)
L = 0; Vn+i = 0; Fh = о.
Рис. 1 - Структурная схема парожидкостных и тепловых потоков в ректификационной колонне
В приведенных уравнениях
кц = ки (Т ■Р . X У), (6)
Н = Н (Т. Р. X) , (7)
Н = НV (Т. РУ) • (8)
где X = (( Чі......ХСі ); У = (. У2і.Усі )■
Тогда задача сводится к решению N(3+20) независимых уравнений, перечисленных
ниже:
Уравнения
Количество
Уравнение покомпонентного материального баланса КхС
Уравнение фазового баланса КхС
Уравнение стехиометрического соотношения 2хК
Уравнение теплового баланса N
Сумма N(3+2^
2. Выбор независимых переменных
Из раздела 1 имеется система N(3+20) независимых уравнений, поэтому необходимо выбрать N(3+20) независимых переменных. В табл. 1 показан список выбранных независимых переменных. Остальные переменные являются исходными данными.
Таблица 1 - Список исходных и независимых переменных
Переменный Число Исход. Независ.
Потоки жидкости и пара ( Ц; V .) 2хК ЦУ, 2хN - 2
Нагрузка теплоты ( 0]) N N - 2
Питание на каждую тарелку ( р. ) N N
Давление питания ( Р'р ) N N
Температура питания ( ТІ ) N N
Состав питаний ( 2(. .) ^С ^С
Боковые отборы ( и і ) 2хN 2хN
Температура на тарелке (Т ■) N N
Давление на каждой тарелке ( Р.) Состав жидкой и паровой фаз на каждой тарел- N N 2х^С
ке ( хи'; У і і ) 2х^С
Сумма ^(3С+10) ^(С+7) ^(2С+3)
3. Метод независимого определения концентраций с использованием трехдиагональной матрицы для расчета процесса ректификации
Идеей этого метода является разделение системы уравнений на три группы (уравнения покомпонентного материального баланса, уравнения фазового равновесия и уравнения теплового баланса), которые решаются не одновременно, а последовательно по группам итерационным методом [1, 2, 3]. Однако его недостатками являются использование простых термодинамических моделей для расчета энтальпии, и не возможность расчета при условии, когда питание содержит жидкую и паровую фазы. В работе [4] предложена модификация этого метода, позволяющая рассчитать процесс ректификации ширококипящих смесей. В данной статье авторы добавили процедуру расчета состава фаз и энтальпии питаний при заданных давлении и температуре [1, 3, 5], а также использовали уравнение состояния СРК [2, 6] для определения коэффициента равновесия пар-жидкость и энтальпии фаз. Это позволяет рассчитать процесс ректификации с питаниями в виде кипящей жидкости, пара или обеих этих состояний, а также обеспечивает высокую точность. На первой итерации коэффициент парожидкостного равновесия компонентов определяется с помощью модели Чао-Сидера.
3.1 Расчет состава фаз и энтальпии в питании
Для каждого питания мы имеем количество ^ и состав 2, 1 при давлении Рур и температуре Т ^ . Тогда необходимо определить количество Ру, Р 1, концентрации 2^^, 2^ и энтальпии Н^, Нр1 паровой и жидкой фаз питания соответственно. Для этого используется метод расчета парожидкостного равновесия [1, 3, 5].
3.2 Начальный профиль температуры и потоков по паровой фазе
Задаем начальный профиль температур линейной интерполяцией между предполагаемой температурой верхней и нижней частями колонны:
ДТ _ Тниз Тверх .
_ N -1 .
т _ Тв
верх
TJ _ у + ЛТ, (2 <у< N).
(9)
Начальный профиль потока по паровой фазе вычисляется из предположения постоянного молярного расхода:
V2 = V + + Uл); 1 = V, (2 <у< N -1). (10)
3.3 Расчет профиля концентрации жидкой фазы
Уравнение общего материального баланса сечения колонны относительно верха имеет
вид:
У-1
(11)
V. -L. 1 + FV +У (FV + FL - 1/У - U V- V. = 0
1 1-11 ¿—1 \ т т т ^т ) 1 •
т=1
Отсюда:
Ц-, = V. + FJ + £(С + Г - 1^т -ит)-V.
т=1
Заменяя у, у и Ll в уравнении (1) соответственно К! ¡Х ■,. (уравнение 2), а также используя выражение (11), получим:
А1Х! , 1-1 + В1Х! 1 + С1Х! , 1+1 = В1 ,
где
1 -1
А. = V- + Fv +У (FV + FL -М - и )-У1
1 1 1 /.и \ т т т т) 1
(12)
т_1
В _-
і+І +£ (С + ^ - и„ )-и, + иJ +(^ + ^ )к„
т_1
У j І j+1 і, j+1
С _ V+к,у+1 (1 < У < N -1); й: _ и -
х,0 _0; х N+1 _0; ^ _0; +1 _0.
(2 < і < N);
(1 < І < N); (1 < І < N);
Таким образом, мы получили систему линейных уравнений в виде трехдиагональной матрицы:
" в С1 0 0 " хі,1" " О "
Аг Вг Сг 0 X ,2 °2
0 Аз Вз Сз 0 х/,з Оз
X _
А В N-1 '-N-1 СМ-1 х/м-1 DN-1
0 Аы BN_ _ XN _ _
При решении этой системы методом Томаса [1, 3], основанным на приведении к диагональному виду путем элементарных преобразований по рекуррентным формулам, получим Х, 1 (1 < 1 < N ). Такой процесс расчета повторяется для каждого компонента.
121
3.4 Профиль концентрации жидкой фазы по высоте колонны
Рассчитанный состав жидкой фазы необходимо пронормировать по формуле:
х
(норм) _ і, і Хі, і С
(1 < і < N ; 1 < і < С).
к _1
3.5 Расчет нового профиля температуры и состава паровой фазы
Для расчета новой температуры и состава паровой фазы можно использовать алгоритм расчета точки кипения [1, 4] на каждой тарелке при известных составе х, і и давлении Рі
жидкой фазы. Уравнение расчета точки кипения дляу-ой тарелки выражается:
¿Пі-1 =ІЛЛі -1 _ 0, (14)
і_1 і_1
где К і ■ - коэффициент равновесия пар-жидкость по уравнению СРК.
3.6 Расчет нового профиля потоков паровой фазы
Значение парового потока У2 (рис. 1) получено по уравнению общего материального баланса для первой тарелки (конденсатор):
у _ у + +и1.
Заменяя Ц-1 и Ц в уравнении теплового баланса (5) выражением (9), получим:
у¡-1 -аі-У-1
V _
Р'і-1
(У = 3
(15)
где а, _ Н|-1 - ну ; ^ _ Н^ - Н\;
У/ _
і-1
ІН - - и„) +ру - у
т_1
(Н - н|_1)+р;, н - н* р
+?Ч « - "? — («Г - нЧ )—
3.7 Расчет профиля потоков жидкой фазы
Из уравнения (11) получим:
Ч=V.+1+^+£ «+« - и т - м ) - V (1 < у < N)
(16)
где К+1 _ 0; ^ _ 0.
3.8 Проверка сходимости расчета
Процедура заканчивается при удовлетворении условия:
-і2
N
І
і _1
Т (к) - т (к-1) і і Т (к)
N
+ І
і _1
У (к) - У (к-1) і і
V(к)
< а.
где е- заданная ошибка. В работе [3] Ванга и Хенка предложен следующий критерий:
N Г_І[Т(к) -т(к-1)] < 0 ,0^, і _1
(17)
(18)
где Т(к) , Т(к 1) - соответственно значения температуры нау-ой тарелке к-ой и (к-1)-ой итераций. 3.9 Расчет отвода и подвода теплоты 01 и Оы
После того, как критерий сходимости выполняется, отвод и подвод теплоты рассчитываются по формулам:
о, = уд*-[уД +(у + и,) Д- ],
N N-1
о^ = I (Г + р Д - и НУ - w ду) -1 - УД* - ¿*Д .
(19)
(20)
7=1 7=1
На рис. 2 приведена блок-схема алгоритма метода независимого определения концентраций с использованием трехдиагональной матрицы для расчета процесса ректификации многокомпонентных систем.
Рис. 2 - Блок-схема алгоритма расчета по методу независимого определения концентраций с использованием трехдиагональной матрицы
4. Поверочный пример алгоритма расчета по методу независимого определения концентраций с использованием трехдиагональной матрицы
В соответствии с алгоритмом была составлена программа расчета на языке Visual Ba-sic.Net. На рис. 3 приведен пример поверочного расчета процесса ректификации смеси, состоящей из 5 компонентов: этана (1), пропана (2), н-бутана (3), н-пентана (4), н-гексана (5). Данные физических свойств этих компонентов взяты из работы [6]. В работе получены результаты: профиль Lj, Vj, Tj, Xjj, y\ j, а также Q-i, Qn (представлены в табл. 2). Данные табл. 2, находящиеся в скобках, представляют собой исходные данные. На рис. 4 представлен график профиля концентрации паровой фазы по высоте колонны.
Рис. 3 - Поверочный пример алгоритма расчета по методу независимого определения концентраций с использованием трехдиагональной матрицы
Таблица 2 - Результаты расчета ( жи, у, ■ в мольных долях)
Номер тарелок L (моль/ч) V (моль/ч) Т (К) О (КДж/ч) и (моль/ч) № (моль/ч)
1 (68039,9) (6803,9) 315,3 957733,5 (2268) 0
3 77799,5 76125,8 327,5 (212000) (1360,8) 0
13 118356,8 84562,2 392,4 0 0 (16782,9)
16 18143,6 96648,9 412,4 -1663362,3 0 0
Номер тарелок Уіі У2І Узі Х4і У4і *5] У5і
1 0,0744 0,1719 0,8805 0,8116 0,0449 0,0165 0,0002 0,0000 0,0000 0,0000
3 0,0182 0,0474 0,7801 0,8595 0,1982 0,0923 0,0035 0,0008 0,0000 0,0000
13 0,0000 0,0001 0,0130 0,0291 0,5656 0,6901 0,3896 0,2694 0,0318 0,0113
16 0,0000 0,0000 0,0013 0,0033 0,2595 0,3926 0,6248 0,5460 0,1144 0,0580
Профиль концентрации паровой фазы
1.0
о.э
0.3
0.7
/■“Ч Ч 0 Р
Ч
и
и.
ч 0.5
и
.2.
0.4
0.3
0.2
0.1
0.0
А \
\ V
\ г \
\ / \ к г
\ / \
У / \ \ У Д.
/ / \ V >
\ У / V ■ ш
К > С -■ ■ ■ -Л—
- У(ч -У(2) -У[3) ■ У(4) -У(5)
1 2 3 4 5 6 7
9 10 11 12 13 14 15 16
Номер тарелок
Рис. 4 - График профиля концентрации паровой фазы по высоте колонны
Вывод
Авторы предложили модификацию алгоритма расчета по методу независимого определения концентраций использованием трехдиагональной матрицы и программное обеспечение для моделирования и расчета многокомпонентной ректификации. Использование процедуры определения состава фаз, уравнения состояния СРК и модели Чао-Сидера для расчета коэффициента равновесия системы пар-жидкость, а также энтальпии фаз, позволяет применить этот алгоритм в программном обеспечении расчета колонны ректификации с различным питанием сырья и отбором жидкой и паровой фаз в разных точках колонны (рис.1).
Литература
1. Кафаров, В.В. Разделение многокомпонентных систем в химической технологии. Методы расчета: учебное пособие / В.В. Кафаров, Ю.Н. Ковалев, В.Н. Ветохин - М.: МХТИ. 1987. - 84 с.
2. Комиссаров, Ю.А. Научные основы процессов ректификации / Ю.А. Комиссаров, Л.С. Гордеев, Д.П. Вент. - М.: Химия. 2004. Т. 2. - 416 с.
3. Ernest, J. Equilibrium-Stage Separation Operations in Chemical Engineering / J. Ernest, Henley, J. D. Seader - New York: John Wiley & Sons Inc. 1981. - 742 pp.
4. Rosendo, M.L. Simulation of Multicomponent Multistage Vapor Liquid Separations. An Improved Algorithm Using the Wang Henke Tridiagonal Matrix Method / M. L. Rosendo - Ind. Eng. Chem. Res. 2003, 42, pp 175-182.
5. Holland, C.D. Fundamentals of multicomponent distillation / C.D. Holland. New York: McGraw-Hill Book Company. 1994. - 626 pp.
6. Robert, C. R. The Properties of Gases and Liquids / C. R. Robert, M. P. John, E. P. Bruce - New York: “McGraw- Hill, Inc”. 1987. - 741 pp.
© Ю. А. Комиссаров - д-р техн. наук, проф., зав. каф. электротехники и электроники РХТУ им. Д.И. Менделеева, [email protected]; Дам Куанг Шанг - асп. той же кафедры.
