Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,441

ФОРМИРОВАНИЕ НЕРАВНОВЕСНОЙ МОДЕЛИ ПРОЦЕССА РЕАКЦИОННОЙ РЕКТИФИКАЦИИ

Индыгашева Е.В. 1 Митянина О.Е. 1
1 ФГБОУ ВПО «Национальный исследовательский Томский политехнический университет»
Управление, моделирование и анализ совмещенных процессов осложнены наличием нескольких стационарных состояний. Существующие методы анализа множественности стационарных состояний основаны на математических моделях процесса.В статье сформирована фундаментальная неравновесная модель реакционно-ректификационного процесса, предназначенная для оптимизации, анализа, проектирования совмещенных процессов. Был проведен расчет изменения количества вещества за счет массопереноса,рассчитаны коэффициенты массопереноса, удельная поверхность контакта фаз и движущая сила процесса реакционной ректификации. Было проведено сравнение полученной неравновесной модели с данными расчета [8] и равновесной моделью. Отклонение результата расчета по неравновесной модели от равновесной модели незначительно, поэтому для проектирования и моделирования совмещенных процессов рекомендовано использовать равновесную модель.
реакционная ректификация
неравновесная модель
метил-трет-бутиловый эфир (МТБЭ)
массоперенос
1. Астарита Дж. Масссопередача с химической реакцией. – М.: Химия, 1974. – 224 с.
2. Динамика нестационарной ректификации, совмещенной с реакцией переэтерификации в жидкой фазе / В.В. Дильман, В.А. Лотхов, А.А. Глаголева, Н.Н. Кулов, С.Я. Квашнин // Теоретические основы химической технологии. – 2002. – Т. 36, № 3. – С. 248.
3. Анализ равновесной и неравновесной моделей реакционно-ректификационной колонны синтеза МТБЭ / М.А. Самборская, О.Е. Митянина, А.В. Кравцов, К.В. Дерина // Нефтегазопереработка: тезисы докл. междунар. научно-практическая конф.(Уфа, 24–27 Мая 2011 г.). – Уфа, 2011 – С. 200–201.
4. Серафимов Л.А. Технология основного органического синтеза. Совмещённые процессы. – М.: Химия, 1993. – 410 с.
5. Ульянов Б.А., Бадеников В.Я., Ликучев В.Г. Процессы а аппараты химической технологии в примерах и задачах. – Ангарск, 2006. – С. 275–285.
6. Doherty M.F., Fidkowski Z.T., Malone M.F., Taylor R. Perry’s chemical engineers handbook, 8th edition, section 13, Distillation. – USA, The McGraw-Hill Companies, 2008. – Р. 50–51.
7. Hauan S., Hertxberg T. Why Methyl tert-Butyl Ether Production by Reactive Distillation May Yield Multiple Solutions – Industrial & Engineering Chemistry Research, 1995, no. 34(3).
8. Huss R.S., Chen F., Malone, Doherty M.F. Reactive distillation for methyl acetate production // Computers and Chemical Engineering. – 2003. – № 27. – Р. 1855–1866.
9. Murphree E.V. Reactive Column Calculations // Ind. Eng. Chem. – 1925. – № 17. – Р. 747–750.
10. Wang, S.J., Wang D.S.H., Lee E.K. Effect of Interaction Multiplicity on Control System Design for a MTBE Reactive Distillation Column // J. Proc. Contr. – 2003. – № 13. – Р. 503–515.

В последние годы в химической промышленности все большее значение приобретает такой технологический вариант процесса, как совместное проведение химический реакции с разделением реакционной смеси в одном и том же аппарате. Рациональное использование явления переноса массы в момент протекания химической реакции обеспечивает дополнительные преимущества как в кинетическом, так и в термодинамическом аспектах [1].

На сегодняшний день реакционно-ректификационные процессы приобретают все большую актуальность в силу того, что они имеют ряд технологических и экономических преимуществ перед стандартными схемами разделения. Таким образом, проектирование, исследование и внедрение данных процессов в производство является актуальной задачей для российской экономики. С целью создания более экономичной технологии, а также технологии безотходных производств очень часто выгодно проводить несколько процессов в одном аппарате.

В настоящее время опубликовано большое количество научных работ, посвященных моделированию и анализу динамического поведения процесса получения МТБЭ в колонне реакционной ректификации, в которых отмечается наличие множественных стационарных состояний с сильно различающимися выходами целевых продуктов, а также неидеальность парожидкостного равновесия [4, 7]. Поэтому качественное проектирование, исследование и оптимизация реакционной ректификации является актуальной научной и технологической задачей, для решения которой можно использовать популярное специализированное программное обеспечение на основе математических моделей.

Существует несколько основных типов моделей реакционной ректификации. Наиболее распространенной является равновесная модель, основным допущением которой является равновесие между жидкостью и паром. Подобные модели разрабатывались как за рубежом – в работах [8], так и в нашей стране [2]. Равновесные стационарная и динамическая модели, благодаря своей наглядности и простоте, успешно используются при проектировании, различных исследованиях и управлении реакционно-ректификационными аппаратами [10], авторы которых полагают, что для качественного исследования динамики реакционно-ректификационной колонны необходимо использовать модели, учитывающие неидеальность парожидкостного равновесия, то есть учитывающие массо- и теплоперенос в аппарате.

Задачами данной работы является формирование неравновесной математической модели реакционной ректификации, а также сравнение результатов расчета материального баланса колонны по равновесной и неравновесной модели.

На кафедре ХТТ и ХК ИПР ТПУ исследовалась равновесная модель реакционно-ректификационной колонны синтеза МТБЭ[3].

При расчете по неравновесной модели необходимо рассчитать изменение количества вещества на каждой из тарелок колонны за счет массопереноса [6].

Eqn61.wmf (1)

где Ni – изменение количества вещества за счет массопереноса; Ki – коэффициент массопереноса, м/с; А – удельная поверхность контакта фаз, отнесенная к единице массы, м2/м2; xi, yi – средние составы жидкости и паров, поднимающиеся с i-й тарелки и поступающие на нее; Eqn62.wmf, Eqn63.wmf – равновесный состав жидкости и пара i-ю тарелку.

Таким образом, для формирования неравновесной модели реакционной ректификации необходимо:

– определить коэффициенты диффузии;

– рассчитать коэффициент массопереноса;

– рассчитать площадь поверхности контакта фаз;

– рассчитать градиент концентрации.

Определение удельной поверхности контакта фаз на тарелках ректификационной колонны

Поверхность контакта фаз на тарелках ректификационной колонны трудно подается определению. Тем не менее разработаны методы и выполнены измерения поверхности. Ниже представлена формула, позволяющая рассчитать удельную поверхность, отнесенную к единице активной площади, для ситчатых тарелок [5].

Eqn64.wmf (2)

где А – удельная поверхность контакта фаз, отнесенная к единице площади тарелки, м2/м2; φг – газосодержание, м3/м3; ΔРст – сопротивление столба жидкости, Н/м2; ΔPст = h0∙ρж∙g; h0 – высота слоя жидкости, м; wk – скорость газа в сечении колонны, м/с.

Высота запаса жидкости (h0) для ситчатых тарелок определяется по уравнению:

Eqn65.wmf (3)

где Vж – объемный расход жидкости, м3/с; lw – длина сливной перегородки, м; σж, σв – поверхностное натяжение жидкости и воды; hw – высота сливной перегородки, м; hf – высота пены, м; fc – доля свободного сечения тарелки [5].

Eqn66.wmf (4)

Eqn67.wmf (5)

В колоннах реакционной ректификации часто используется ситчатый тип тарелок.

Данные о характеристике ситчатых тарелок взяты из [5].

Таблица 1

Исходные данные

Температура, °С

Молярная масса, г/моль

Вязкость, Па∙с

Поверхностное натяжение, Н/м

Диаметр колонны, м

dэ = d0, м

ρж, г/см3

ωk, м/с

46,85

32

0,4077

26,126

1,4

0,04

0,78

6695

Пользуясь формулами (2)–(5), рассчитали удельную поверхность контакта фаз. Полученные результаты представлены в табл. 2.

Таблица 2

Расчетные данные удельной поверхности

Высота пены, м

Запас жидкости на тарелках, м

Газосодержание, м3/м3

Сопротивление столба жидкости, Н/м2

Удельная поверхность контакта фаз, м2/м2

16,434

16,081

0,022

122,97

9,204

После проведенных расчетов, получили удельную поверхность контакта фаз, отнесенной к единице площади тарелки, равной 9,204 м2/м2.

Расчет коэффициента массопереноса

Производительность массопереноса на тарелке часто выражается путем безразмерной группы, называемой числом единиц переноса [6].

Эти безразмерные числа определяются с помощью уравнений:

Eqn68.wmf (6)

где hf – высота пены, м; Z – длина пути потока жидкости, м; Eqn69.wmf – объемная скорость потока жидкости, м3/с; us = V/Eqn70.wmf – поверхностная скорость пара, м/с; Аbub – кипящий слой на тарелке, м2; hL – высота жидкости на тарелке, м; a’ – межфазная поверхность на единицу поверхности пара, м2/м3; Eqn71.wmf – межфазная поверхность на единицу поверхности жидкости, м2/м3.

Эти межфазные поверхности связаны с межфазной поверхности на единицу объема пены следующим образом:

Eqn72.wmf (7)

где a = hL/hf – относительная плотность пены на тарелке. Также tL и tV время пребывания в паровой и жидкой фазах, определяются:

Eqn73.wmf (8)

Eqn74.wmf (9)

Для определения числа единиц переноса мы используем корреляционную модель AIChE. В этой модели выражение для определения числа единиц переноса выглядит следующим образом:

Eqn75.wmf (10)

где hw – высота перегородки, м; ScV – число Шмидта, равное ScV = μV/(ρVDV); Fs – фактор, кг/мс; Fs = usV(ρV)0,5; QL – объемная скорость потока жидкости, м3/с, Eqn69.wmf [6].

Ниже представлена таблица с данными для расчета коэффициентов массопереноса.

Таблица 3

Данные для расчета коэффициента массопереноса

Величина

Обозначение

Значение

Длина пути потока жидкости

Z, м

0,4835

Длина перегородки

W, м

1,2

Высота жидкости на тарелке

hL, м

0,079

Высота пены

hf, м

0,0715

Высота сливной перегородки

hw, м

0,05

Площадь поперечного сечения колонны

Ас, м2

1,538

Данные для расчетов по тарелкам взяты из справочника программы Хайсис. Пользуясь формулами (6)–(11), рассчитали поэтапно коэффициенты массопереноса, средние значения коэффициентов массопереноса по тарелкам: kV = 0,491, kL = 147 м/с.

Определение движущей силы процесса

Движущую силу процесса рассчитываем, учитывая эффективность контактного устройства, используя закон Мерфри [9]:

Eqn76.wmf (11)

где уk, уk-1 – средние составы паров, поднимающиеся с k-й тарелки и поступающие на нее; Eqn77.wmf – состав пара, равновесный с составом жидкости, покидающей k-ю тарелку; хk, хk-1 – средние составы жидкости, покидающие k-ю тарелку; Eqn78.wmf* – состав жидкости, равновесный с составом пара, поднимающегося на k-ю тарелку.

После преобразований получили уравнение для расчета рабочей концентрации компонентов в паровой и жидкой фазах:

Eqn79.wmf (12)

Эффективность контактного устройства приняли среднее для ситчатых тарелок: Е = 0,8;

При расчете движущей силы процесса использовались равновесные концентрации на тарелках колонны синтеза МТБЭ, представленные в [7].

Расчет массопереноса

Сформированная неравновесная модель представлена ниже:

Eqn80.wmf (13)

где i – индекс компонента; j – индекс номера тарелки; x, y – мольная доля в жидкой и паровой фазах соответственно; V, L, F – расход пара, расход жидкости, расход орошающей жидкости соответственно, моль/с; U – удержание жидкости на тарелке, м/с;vi,ρ– стехиометрический коэффициент j-го вещества в реакции с номером ρ; ri,r – скорость химической реакции, моль/(м3·с).

Результаты расчета по неравновесной модели в сравнении с данными, представленными в [7] и с результатами расчета по равновесной модели для жидкой и паровой фазы, приведены на рис. 1 и 2 соответственно.

pic_60.wmf pic_61.wmf
pic_204.wmf pic_205.wmf

Рис. 1. Сравнение неравновесной модели с экспериментальными данными в жидкой фазе

Как видно из рис. 1, среднее отклонение расчетов в жидкой фазе по неравновесной модели от данных расчета [7] не превышает 3 %.

Как видно из рис. 2, отклонение расчета неравновесной модели от данных расчета [7] в паровой фазе превышает отклонение от данных расчета [7] в жидкой фазе, и можно сделать вывод, что неравновесная модель адекватно описывает процесс.

Сравнение результатов расчета неравновесной модели с результатами расчета по равновесной модели для жидкой и паровой фазы приведены на рис. 3 и 4 соответственно.

pic_62.wmf pic_63.wmf
pic_64.wmf pic_65.wmf

Рис. 2. Сравнение неравновесной модели с экспериментальными данными в паровой фазе

pic_66.wmf pic_67.wmf
pic_68.wmf pic_69.wmf

Рис. 3. Сравнение неравновесной и равновесной модели в жидкой фазе

Как видно из рис. 3 и 4, неравновесная модель практически не отличается от значений равновесной модели. Отклонения от значений равновесной модели составили не более чем 4 %, поэтому целесообразно использовать для расчетов, моделирования, проектирования, анализа множественности стационарных состояний и оптимизации равновесную модель.

pic_70.wmf pic_71.wmf
pic_72.wmf pic_73.wmf

Рис. 4. Сравнение неравновесной и равновесной модели в паровой фазе

Рецензенты:

Ивашкина Е.Н., д.х.н., доцент кафедры химической технологии топлива и химической кибернетики, ФГБОУ ВПО «Национальный исследовательский Томский политехнический университет», г. Томск;

Иванчина Э.Д., д.т.н., профессор кафедры ХТТиХК, ФГБОУ ВПО «Национальный исследовательский Томский политехнический университет», г. Томск.

Работа поступила в редакцию 01.07.2013.


Библиографическая ссылка

Индыгашева Е.В., Митянина О.Е. ФОРМИРОВАНИЕ НЕРАВНОВЕСНОЙ МОДЕЛИ ПРОЦЕССА РЕАКЦИОННОЙ РЕКТИФИКАЦИИ // Фундаментальные исследования. – 2013. – № 8-3. – С. 615-620;
URL: http://fundamental-research.ru/ru/article/view?id=31968 (дата обращения: 31.10.2020).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074