УДК 519.711.3
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ РЕАКЦИОННОЙ СИСТЕМЫ ПОЛИМЕРИЗАЦИИ ДИЦИКЛОПЕНТАДИЕНА
В.Г. Бондалетов, А.А. Ляпков, Е.И. Мельник, Г.В. Борисов
Томский политехнический университет E-mail: alexdes@tpu.ru
На основании математической модели полимеризации дициклопентадиена предложена математическая модель реактора-полимеризатора. Смоделирован реактор для производства полидициклопентадиена, работающий в адиабатическом и изотермическом тепловых режимах. Определена максимально допустимая концентрация мономера в растворе толуола. Показано, что продукты реакции полимеризации в разных реакторах имеют одинаковую молекулярную массу и различное молекулярно-массовое распределение.
Ключевые слова:
Полидициклопентадиен, адиабатический реактор, математическое моделирование, молекулярная масса, молекулярно-массовое распределение.
Решение проблемы утилизации побочных продуктов производств низших олефинов (этилена и пропилена), образующихся при пиролизе различного углеводородного сырья, является важной инженерно-экономической задачей при организации эффективного производства. Основным побочным продуктом, получающимся при пиролизе бензина, являются жидкие продукты пиролиза, содержащие в своем составе ценные продукты [1].
Один из важнейших компонентов жидких продуктов пиролиза - дициклопентадиен, содержание которого в различных фракциях может достигать 50 %. Полимер, получаемый на основе дициклопентадиена, является хорошим конструкционным пластиком, однако единственным способом получения изделий из этого вещества является реак-ционно-инжекционное формование [2].
Альтернативной технологией получения полидициклопентадиена является полимеризация его под действием каталитической системы А1Е12С1—Т1С14 в трубчатом реакторе вытеснения. При этом процесс полимеризации протекает по катионному механизму с раскрытием одной из двойных связей в кольцах дициклопентадиена и сопровождается образованием полимера линейного строения, способного к однократной переработке, так как сшивка и структурирование полимерных цепей в этом случае эффективно протекает при температурах выше 100...120 °С [3].
Реактор вытеснения рассматривают как модель идеального вытеснения. Условие идеальности такого аппарата состоит в том, что каждый элемент реакционной массы в данном поперечном сечении движется вдоль оси потока с одинаковой линейной скоростью (поршневой режим). Режим течения, близкий к идеальному вытеснению, имеет место в длинных трубах постоянного поперечного сечения при значениях критерия Рейнольдса, характерных для развитой турбулентности [4]. При этом отношение длины трубы к ее эквивалентному диаметру должно быть больше 20.
Характерным представителем такого типа реакторов является трубчатый реактор (рис. 1).
Рис. 1. Материальный баланс трубчатого реактора: г - координата длины; L - длина реактора
Полное математическое описание процесса представляется покомпонентным материальным балансом и тепловым балансом элементарной ячейки йУР для малого промежутка времени йт, так как параметры потока меняются по длине аппарата и во времени.
Конечный вид уравнения материального баланса реактора идеального вытеснения для вещества в размерности (моль/(л-с)) будет иметь вид:
СС СР Ст СУР
где г - скорость изменения концентрации вещества в результате химического превращения, моль/(л-с); ¥=уС=>БС - мольный поток вещества, моль/с; V -объемный расход потока, м3/с; С - концентрация вещества, моль/л; V - линейная скорость потока, м/с; 5 - площадь поперечного сечения реактора, м2. Решением этого уравнения является функция
с=/УТ сс
В случае стационарного режима работы — = 0
ёт
и уравнение материального баланса упрощается до
й¥
----= г. Решением его будет функция С=ДУр).
СУР
При УР=ео^ (мономолекулярная реакция) и Т =Т=Т :
О V ОІ.Т V •
dF
dVP
d(USC) _ UdC_ d (Sz) dz
dC dC
_ -U---------------+ r
dT dz
Решение этой системы уравнений - функция С=Дг,т).
и йС 0
В случае стационарного режима — = 0 и
йт
йт = — = (при постоянных и и Vпо длине реак-
ч йС г йС
тора) — = — и С=Дг) или — = г и С=Дт). В этих
й2 и йт
уравнениях I - длина реакционной зоны, м.
Результаты и их обсуждение
Исходя из кинетических данных, полученных для изучаемой системы (дициклопентади-ен-А1Е^С1-Т1С14-толуол) [5], составили математическую модель полимеризации, которая включает стадии равновесного образования комплексов каталитической системы, активной частицы, а также роста и обрыва цепи.
Стадия роста цепи заключается в последовательном присоединении молекулы мономера к активному центру, представляющему собой сольватно-разделенную ионную пару. Образование полимера происходит путем переноса цепи на мономер или катализатор:
DCPD + DEAH KOMP + TCT
± AC;
DCPD + AC —— GAC; GAC + DCPD —— POL + AC;
‘ POL + AC,
оксида пропилена, при раскрытии оксидного цикла которого протекают реакции дезактивации кислот Льюиса с образованием алкоксидов титана и алюминия, неактивных в катионной полимеризации.
Температура, достигаемая в адиабатическом реакторе вытеснения, определяется концентрацией дициклопентадиена в растворе толуола. Выбор правильной концентрации мономера, при которой максимальная температура в реакторе составит 60 °С, является необходимым условием для определения оптимальных значений расхода раствора мономера и диаметра реактора, поскольку при температурах свыше 60 °С образующийся полимер начинает сшиваться, что ухудшает растворимость продукта в ароматических и хлорированных растворителях. Кроме того, при повышенных температурах образуется большое количество коротких полимерных цепей, что приводит к общему снижению молекулярной массы образующихся продуктов, а также к расширению молекулярно-массового распределения.
В качестве растворителя использовали толуол. Концентрация мономера составляла 1530 моль/м3, а каталитического комплекса - 8,9 моль/м3. Расход мономера в реактор рассчитывали исходя из данных, полученных с производства «ЭП-300» (Производство мономеров ООО «Томскнефте-хим»), учитывая, что его содержание во фракции С9 составляет 16,6 %.
Математическое описание реактора состоит из уравнения материального баланса по мономеру
GAC + TCT
где DCPD - мономер, дициклопентадиен; TCT - тетрахлорид титана; DEAH - диэтилалюминийхло-рид; KOMP - сольватированная молекула катализатора, включая ценовые комплексы тетрахлорида титана; AC - первоначальный активный центр; GAC - растущие активные центры; POL - полимер, полидициклопенадиен.
Математическая модель была использована для проверки кинетической схемы реакции путем сравнения экспериментальных результатов с теоретическими, полученными расчетным путем. Решение обратной кинетической задачи позволило из данных кинетического эксперимента найти недостающие параметры кинетической модели и уточнить значения параметров, полученных экспериментальным путем.
Для моделирования нами был выбран реактор идеального вытеснения, который не отягощен продольным и обратным перемешиванием. Для такой модели несложно составить математическое описание, которое в общем случае должно включать в себя уравнения материального и теплового балансов, а также граничные условия, налагаемые на модель из практических соображений. Ограничение роста цепи проводили путем добавления дезактиватора компонентов каталитической системы -
UdCu
dT
_-k0e
C
'“'Л/
и уравнения теплового баланса
UdT_
dT
Ea
rt C
k0 e
Ea,
RT C
Поскольку при составлении математического описания реактора идеального вытеснения применяют систему координат, движущуюся с потоком (^=еоияі), можно записать
dCM
dT
_ - k0 e
AP
RT C
dT_
dT
AH„, -^
Lk0Pe
Cm +
k0 e RTCK
Для изотермического реактора Г=сопб1;.
Данную систему обыкновенных дифференциальных уравнений решали для различных начальных условий (расход мономера и диаметр реактора). Для трех значений V и четырех значений Л результаты моделирования представлены на рис. 2.
Режим течения, близкий к идеальному вытеснению, имеет место в длинных трубах постоянного поперечного сечения при значениях критерия Рейнольдса, характерных для развитой турбулент-
Е
A
а
/. м
б
I, М
I, М
Рис. 2. Результаты моделирования степени превращения в полимеризации дициклопентадиена при концентрации мономера 1530моль/м3 под действием каталитической системы AlEt2Cl~TiCl4 (1:1 моль, 2,0 мас. %) при расходе: 1) 4,1; 2) б,1; 3) 8,1 м3/ч для изотермического реактора диаметром: а) 0,05; б) 0,1; в) 0,15; г) 0,2 м
в
г
ности. При этом отношение длины трубы к ее эквивалентному диаметру должно быть больше 20. Из рис. 2 следует, что оптимальным расходом для реакторов различного диаметра является 4,1 м3/с, так как при этом значении наиболее быстро достигается максимальная мольная доля полимера, однако, учитывая критерий Рейнольдса, только расход 6,1 м3/с обеспечивает адекватное решение задачи для всех значений диаметра аппарата.
При одинаковых значениях расхода мономера и длины реактора, учитывая необходимый критерий, максимальная мольная доля полимера достигается в аппарате диаметром 0,2 м (рис. 3).
Важным условием для проведения полимеризации является расчет места ввода катализатора в реактор. В адиабатическом реакторе диаметром 0,1 м место ввода катализатора расположено на расстоянии 0,21 м от впрыска мономера и АШ2С1,
так как максимальная концентрация каталитического комплекса образуется через 1,5 с после начала процесса (рис. 4).
Проведен анализ молекулярной массы полученного полимера и молекулярно-массового распределения, которые находятся в прямой зависимости от скорости образования активных центров, возникающих под действием каталитического комплекса, и достижения системой стационарного состояния.
Анализ математической модели показал, что в изотермическом реакторе стационарное состояние достигается быстрее, чем в адиабатическом при равных концентрациях мономера и объемном расходе раствора (рис. 4), следовательно, молекулярно-массовое распределение полученного полимера в этом реакторе будет уже, в то время как средняя молекулярная масса полимера в обоих аппаратах будет одинаковой.
Рис. 3. Результаты моделирования полимеризации дициклопентадиена в изотермическом реакторе под действием каталитической системыы А1(С2Н5)2й~Т1С!4 при расходе мономера 6,1 м3/с и диаметре реактора: 1)0,2; 2) 0,15; 3) 0,1; 4) 0,05 м
Рис. 4. Кинетическая кривая образования первоначальных активных центров в адиабатическом (1) и изотермическом (2) реакторах диаметром 0,2 м, расходе раствора мономера 6,1 м3/ч при полимеризации ДЦПД концентрацией 1530 моль/м3 под действием каталитической системы АІЕі2СІ~ТІСІ4 (1:1 моль, 2,0 мас. %)
Выводы
Предложена математическая модель изотермического реактора вытеснения для проведения растворной полимеризации дициклопентадиена в толуоле под влиянием каталитической системы Т1С14-А1(С2Н5)2С1 (1:1 моль, 2,0 мас. %).
Смоделированный изотермический реактор полимеризации дициклопентадиена при достижении максимальной степени превращения имеет диаметр 0,2 м и длину 4,8 м. При этом оптимальными параметрами являются расход раствора мономера
6,1 м3/ч и концентрация дициклопентадиена 1530 моль/м3, при которых достигается максимальная степень превращения.
Полученные данные могут быть использованы при математическом моделировании и расчете кинетических параметров олигомеризации различных фракций жидких продуктов пиролиза, содержащих дициклопентадиен, а также для разработки технологических процессов получения полидициклопентадиена и нефтеполимерных смол.
СПИСОК ЛИТЕРАТУРЫ
1. Думский Ю.В., Но Б.И., Бутов Г.М. Химия и технология нефтеполимерных смол. - М.: Химия, 1999. - 312 с.
2. Vervacke D. An introduction to PDCPD. - Waarschoot: Product Rescue, 2008. - 129 p.
3. Ионова Е.И., Ляпков А.А., Бондалетов В.Г. Процессы окисления и структурирования олигомеров на основе дициклопентадиена // Известия Томского политехнического университета. - 2011. - Т. 318. - № 3. - C. 101-105.
4. Царева З.М., Товажнянский Л.Л., Орлова Е.И. Основы теории химических реакторов. - Харьков: Изд-во ХГПУ, 1997. -624 с.
5. Ионова Е.И. Олигомеризация непредельных компонентов жидких продуктов пиролиза под действием каталитической системы Е12А1С1—Т1С14: дис. ... канд. хим. наук. - Томск, 2012.- 126 с.
Поступила 19.09.2013 г.
UDC 519.711.3
MATHEMATICAL MODELING OF THE REACTION SYSTEM OF DICYCLOPENTADIENE POLYMERIZATION
V.G. Bondaletov, A.A. Lyapkov, E.I. Melnik, G.V. Borisov Tomsk Polytechnic University
Based on the mathematical model of dicyclopentadiene polymerization the authors have proposed the mathematical model of the polymerization reactor. The reactor for producing polydicyclopentadiene was designed. It is operated in adiabatic and isothermal temperature conditions. The maximum permissible monomer concentration in toluene solution was determined. It was shown that the polymerization products in different reactors have the same molecular mass and different molecular weight distribution.
Key words:
Polydicyclopentadiene, isotermic reactor, mathematic modeling, molecular mass, molecular weight distribution.
REFERENCES
1. Dumskiy Yu.V., No B.I., Butov G.M. Khimiya i tekhnologiya nef-tepolimernykh smol (Chemistry and technology of polymeric petroleum resin). Moscow, Khimiya, 1999. 312 p.
2. Vervacke D. An introduction to PDCPD. Waarschoot, Product Rescue, 2008. 129 p.
3. Ionova E.I., Lyapkov A.A., Bondaletov V.G. Bulletin of the Tomsk Polytechnic University, 2011. 318, 3, pp. 101-105.
4. Tsareva Z.M., Tovazhnyanskiy L.L., Orlova E.I. Osnovy teorii khimicheskikh reaktorov (Foundation of chemical reactor theory). Har’kov, HGPU, 1997. 624 p.
5. Ionova E.I. Oligomerizatsiya nepredelnykh komponentov zhidkikh produktov piroliza pod deystviem kataliticheskoy sistemy Et2AlCl-TiCl4 (Oligomerization of unsaturated components of liquid pyrolisis products under the effect of Et2AlCl-TiCl4 catalytic system). Dissertatsiya kandidata khimicheskikh nauk. Tomsk, 2012. 126 p.
УДК 667.6
ПОЛУЧЕНИЕ ПЛЕНКООБРАЗУЮЩИХ НА ОСНОВЕ ЦИКЛО-, ДИЦИКЛОПЕНТАДИЕНСОДЕРЖАЩИХ ФРАКЦИЙ ЖИДКИХ ПРОДУКТОВ ПИРОЛИЗА
А.А. Мананкова, Е.И. Задорожная, Н.В. Власова
Томский политехнический университет E-mail: nadezhda_hope89@mail.ru
Работа посвящена качественной переработке побочных жидких продуктов пиролиза путем их полимеризации с образованием нефтеполимерных смол. В качестве исходного сырья использовали цикло-, дициклопентадиенсодержащие фракции жидких продуктов пиролиза прямогонного бензина. Был проведен синтез нефтеполимерных смол с использованием каталитического комплекса на основе моноалкокситрихлорида титана Ti(OC3H7)Cl3 и диэтилалюминийхлорида Al(C2H5)2Cl. Моноалкокситрихло-рид титана получали взаимодействием TiCl4 и изопропилового спирта при мольном соотношении 1:1. Установлено, что использование данного каталитического комплекса позволяет получать нефтеполимерные смолы с улучшенными физико-химическими характеристиками. Исследованы свойства пленок на основе синтезированных нефтеполимерных смол, а также эксплуатационные характеристики пигментированных и масляно-смоляных покрытий. Показана возможность использования полученных пленкообразующих в лакокрасочной промышленности как самостоятельно, так и в различных композициях.
Ключевые слова:
Жидкие продукты пиролиза, нефтеполимерные смолы, пленкообразующие, переработка побочных продуктов, дициклопента-диеновая фракция, композиции, лакокрасочный материал.
Нефтеполимерные смолы (НПС) - низкомолекулярные полимерные продукты, которые образуются в результате полимеризации непредельных соединений при качественной переработке жидких продуктов пиролиза с температурой кипения 30...200 °С (фракция С5-С9). НПС применяются в лакокрасочной, мебельной, целлюлозно-бумажной, резинотехнической промышленности, в производстве печатных красок, искусственных кож, антикоррозийных и защитных покрытий, в строительстве и различных композициях.
Так, в работе [1] получена композиция гидро-фобизирующего состава, содержащая водную эмульсию НПС в органическом растворителе и анионактивное поверхностно-активное вещество, в качестве разбавителя воду, а также дополнительно водную дисперсию нефтяного парафина. Данная композиция обладает хорошими гидрофобными свойствами и используется для получения водоотталкивающих покрытий на керамических материалах, например кирпичах.