close

Вход

Забыли?

вход по аккаунту

?

Моделирование процесса сополимеризации бутадиена и изопрена на катализаторе Циглера-Натта TiCl 4Al(i – c 4H 9) 3.

код для вставкиСкачать
Вестник СГТУ. 2012. № 1 (64). Выпуск 2
УДК 541.64.11+519.91
Э.Р. Гиззатова, С.И. Спивак
МОДЕЛИРОВАНИЕ ПРОЦЕССА СОПОЛИМЕРИЗАЦИИ БУТАДИЕНА
И ИЗОПРЕНА НА КАТАЛИЗАТОРЕ ЦИГЛЕРА-НАТТА TICL4- AL(I – C4H9)3
В работе рассмотрена кинетическая схема процесса сополимеризации
бутадиена и изопрена на титансодержащем катализаторе. Построена
математическая модель процесса, поставлена прямая и обратная
кинетические задачи.
Сополимеризация,
кинетическая задача
катализаторы
Циглера-Натта,
диены,
прямая
E.R. Gizzatova, S.I. Spivak
MODELLING OF SOPOLYMERIZATION PROCESS OF BUTADIENE
AND IZOPREN ON ZIGLER-NATTA CATALYSTS TICL4- AL(I − C4H9)3
This paper is shown kinetics scheme of sopolymerization process of
butadiene and izopren on titan-based Zigler-Natta catalysts. The mathematical
model also direct and inverse kinetics problems are given.
Sopolymerization, Zigler-Natta catalysts, direct kinetics problem
Изучение процессов сополимеризации явилось логическим продолжением исследований
процессов гомополимеризации. Механизм процессов подобный, однако, совместное участие мономеров может привести к образованию сополимерного продукта с нехарактерными для участвующих мономеров свойствами. Процессы гомополимеризации довольно хорошо изучены и
накоплен большой литературный материал. В нашей стране существует несколько научных
школ, которые занимаются различными проблемами процессов полимеризации. К рассматриваемым ими вопросам можно отнести исследование кинетической неоднородности катализатора,
дезактивацию активных центров в процессе, появление периода индукции для активации центров и т.д. Решение этих актуальных задач требует создания эффективной методики исследования и предложения оптимального алгоритма. Успешное решение любой научной задачи невозможно без соединения теории и практики, поэтому сочетание лабораторного эксперимента и его
дальнейшая математическая обработка позволяют достичь поставленных целей.
В работе рассматривается процесс сополимеризации бутадиена и изопрена на катализаторе TiCl4-Al(i – C4H9)3, предполагающий три элементарные стадии. Стадия инициирования
активных центров, которая проводится до начала процесса, считается прошедшей. В этом
случае за активные центры, которые будут осуществлять рост цепи, принимаются активные
цепи с одним звеном мономера – бутадиена или изопрена. Каждая элементарная стадия характеризуется своей константой скорости. Поскольку процесс считается изотермическим, то
константы скоростей стадий не меняют своего значения в ходе процесса.
Рост сополимерной цепи осуществляется согласно
k
k
Ri + B → RiB+1 Ri + B → RiB+1 .
Перенос цепи на мономер
B
p
66
B
p
Новые материалы и технологии
BB
kM
BI
kM
RiB + B → Qi + R1B RiB + I → Qi + R1I ,
IB
(1)
II
kM
kM
RiI + B →
Qi + R1B RiI + I →
Qi + R1I .
Перенос цепи на алюминийорганическое соединение
B
I
kA
kA
RiB + A →
Qi + R1B RiI + A →
Qi + R1I .
Здесь RiB и RiI – активные цепи длиной i с концевыми звеньями бутадиена и изопрена соответственно; Qi – неактивная сополимерная цепь длины i с произвольным концевым звеном; B и I –
мономеры бутадиена и изопрена соответственно; А – алюминийорганическое соединение (АОС).
Построение математической модели можно вести двумя способами: статистическим и
кинетическим. Применение первого подхода обосновано в тех случаях, когда необходимо
рассчитывать полимеры с очень большой молекулярной массой на основании заданных вероятностей наступления каждой стадии. Во втором случае математическая модель строится на
основании кинетической модели.
Согласно основным законам химической кинетики, скорость каждой элементарной
стадии определяется по следующей формуле:
w = k[ A1 ][ A2 ]...[ An ] ,
(2)
где w – скорость реакции; [ A1 ],[ A2 ],...,[ An ] – концентрации молекул A1 , A2 ,..., An в моль/л;
k – константа скорости реакции, называемая константой скорости реакции.
Тогда концентрации реагентов процесса C (t ) будут вычисляться по формуле (3):
m
dC dt = ∑ j =1 b ⋅ w j ,
(3)
в которой коэффициент b определяется как 1, если реагент образуется, (– 1), если реагент
расходуется и 0, если реагент в реакции не участвует.
В каждой элементарной стадии схемы (1) участвуют активные цепи с различной длиной цепи. Теоретически процесс образования предполагается бесконечным. Однако быстрое
расходование мономера (глобул полимера) или частичная гибель активных центров, или другие условия среды ведут к прекращению процесса, то есть, всегда будет существовать максимально допустимая длина цепи imax . Следовательно, количество участвующих реагентов будет ограничено.
Учитывая формулы (2), (3) и сделанные выше выводы, можно выписать систему
обыкновенных дифференциальных уравнений, в которой каждое уравнение показывает скорость изменения отдельно взятого реагента во времени:
∞
∞
dB dt = − k pB BC aB − k pB BC aI − k MBB B ∑2 ( RiB ) − k MIB B ∑ ( RiI ),
2
∞
∞
dI / dt = − k IC − k IC − k I ∑2 ( R ) − k I ∑2 ( R ), dA / dt = −k AB AC aB − k AI AC aI ,
I
p
I
a
I
p
B
a
II
M
I
i
BI
M
B
i
∞
∞
dR1B / dt = − k pB BR1B − k pB BR1I − k MBI IR1B + k MIB BC aI + k MBB B ∑2 ( RiB ) + k AB A∑2 ( RiB ),
dRiB / dt = − k pB BRiB + k pB BRiB−1 − k pI IRiB + k pB BRiI − k MBI IRiB − k MBB BRiB − k AB ARiB , i = 2,3...,
∞
∞
(4)
dR / dt = −k IR − k IR − k BR + k IC + k I ∑2 (R ) + k A∑2 (R ),
I
1
I
p
B
1
I
p
I
1
IB
M
I
1
BI
M
B
a
II
M
I
i
I
A
I
i
dRiI / dt = −k pB BRiI + k pI IRiB−1 − k pI IRiI + k pI IRiB − kMIB BRiI − kMII IRiI − k AI ARiI , i = 2,3...,
dQ1 / dt = kMBB BR1B + k MBI IR1B + kMIB BR1I + k MII IR1I + k AB AR1B + k AI AR1I ,
dQi / dt = kMBB BRiB + k MBI IRiB − kMIB BRiI + kMII IRiI + k AB ARiB + k AI ARiI , i = 2,3... .
Начальные данные для нее имеют вид:
67
Вестник СГТУ. 2012. № 1 (64). Выпуск 2
B(0) = B (0) , I (0) = I (0) , A(0) = A (0) , R1B (0) = C aB (0), R1I (0) = C aI (0), Q1 (0) = 0,
(5)
RiB (0) = RiI (0) = Qi (0) = 0, i ≥ 2.
Здесь ненулевыми концентрациями являются концентрации бутадиена B(0), изопрена
I(0), АОС А(0) и концентрации активных центров с одним звеном бутадиена R1B (0) и одним
звеном изопрена R1I (0) , поскольку именно эти реагенты экспериментаторы считают необходимыми для процесса. Концентрации остальных веществ в начальный момент времени считаются равными нулю.
Тогда можно поставить прямую и обратную кинетические задачи. Прямая задача будет
заключаться в решении системы дифференциальных уравнений (4) с начальными данными
(5) и известными значениями констант скоростей элементарных стадий. Обратная задача сводится к поиску констант скоростей элементарных стадий на основании экспериментальных
данных по средним молекулярным характеристикам сополимера.
В работе [2] приведена приведена методика решения прямой кинетической задачи для
процессов гомополимеризации, при этом получены количественные оценки молекулярномассового распределения (ММР). Однако, из-за больших значений констант скоростей было
невозможно рассчитать и построить ММР.
Наличие мощных компьютеров, способных хранить и обрабатывать большие объемы
данных позволяют справиться с этой проблемой и увидеть молекулярно-массовое распределение сополимера, для этого нами создан алгоритм и написано программное обеспечение. Результаты модельного вычислительного эксперимента приведены на рис. 1, 2.
На рис. 1. показано изменение среднечисленной молекулярной массы (левый график)
и среднемассовой молекулярной массы (правый график). Такое поведение кривых характерно
как для гомополимера, так и для сополимера.
Рис. 1. Зависимость средних молекулярных масс от времени сополимеризации
68
Новые материалы и технологии
Следующий рис. 2 иллюстрирует кривые молекулярно-массового
распределения для больших времен
сополимеризации.
В рамках работы можно сделать следующие выводы:
– построена
математическая
модель процесса сополимеризации
бутадиена и изопрена на одноцентровом катализаторе TiCl4- Al(i – C4H9)3;
– предложена методика решения прямой кинетической задачи
нахождения значений средних молекулярных масс и построения молекулярно-массового распределения;
– полученные графические зависимости для средних молекулярных
Рис. 2. Кривые молекулярно-массового
масс и молекулярно-массового расраспределения для 1, 5 и 10 минут
пределения согласуются имеют вид,
сополимеризации
аналогичный с экспериментально получаемым.
Работа выполнена при поддержке ФЦП «Научные и научно-педагогические кадры инновации России» на 2009-2013 гг. ГК № 02.740.11.0648 от 29.03.2010.
ЛИТЕРАТУРА
1. Кафаров В.В., Ветохин В.Н., Тихомиров С.Г. Моделирование кинетики процесса
полимеризации полиизопренового каучука // Доклады АН. 1989. Т.305. № 6. С.1425-1429.
2. Усманов Т.С., Максютова Э.Р., Спивак С.И. Математическое моделирование процесса полимеризации бутадиена на лантанидсодержащих катализаторах // Доклады АН. 2002.
Т.387, № 6. С.793-796.
Гиззатова Эльвира Раисовна –
кандидат физико-математических наук, доцент кафедры «Физическая химия и химическая
экология» Башкирского государственного университета
Спивак Семен Израилевич –
доктор физико-математических наук, заведующий кафедрой «Математическое моделирование» Башкирского государственного университета
Статья поступила в редакцию 6.03.12, принята к опубликованию 12.03.12
69
Документ
Категория
Без категории
Просмотров
5
Размер файла
167 Кб
Теги
моделирование, бутадиен, сополимеризация, циглера, катализаторы, процесс, изопрена, 4al, натти, ticl
1/--страниц
Пожаловаться на содержимое документа