Научная статья на тему 'Решение прямой и обратной кинетических задач для процесса эмульсионной сополимеризации бутадиена со стиролом'

Решение прямой и обратной кинетических задач для процесса эмульсионной сополимеризации бутадиена со стиролом Текст научной статьи по специальности «Химические науки»

CC BY
245
68
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
сополимеризация / бутадиен / стирол / кинетика / обратная задача / молекулярно-массовое распределение / Copolymerization / Butadiene / styrene / kinetic / inverse problem / Molecular-weight distribution

Аннотация научной статьи по химическим наукам, автор научной работы — Мифтахов Э. Н., Насыров И. Ш., Мустафина С. А., Спивак С. И.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по химическим наукам , автор научной работы — Мифтахов Э. Н., Насыров И. Ш., Мустафина С. А., Спивак С. И.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

A mathematical model of emulsion copolymerization of butadiene with styrene is presented in the article. There is solved a direct kinetic problem for calculating characteristics of the molecular-weight distribution. There is also solved an inverse problem by fitting the conversion-dependence to search for the missing parameter. A comparative analysis of dependence of the experimental and calculated values of intrinsic viscosity, number-average and weight average molecular weight from time are studied.

Текст научной работы на тему «Решение прямой и обратной кинетических задач для процесса эмульсионной сополимеризации бутадиена со стиролом»

УДК 517.977.5

РЕШЕНИЕ ПРЯМОЙ И ОБРАТНОЙ КИНЕТИЧЕСКИХ ЗАДАЧ ДЛЯ ПРОЦЕССА ЭМУЛЬСИОННОЙ СОПОЛИМЕРИЗАЦИИ БУТАДИЕНА СО СТИРОЛОМ

© Э. Н. Мифтахов1*, И. Ш. Насыров2, С. А. Мустафина1, С. И. Спивак3

1 Стерлитамакская государственная педагогическая академия им. З. Биишевой Россия, Республика Башкортостан, 453I03 г. Стерлитамак, пр. Ленина, 49.

Тел./факс: +7 (3473) 43 I0 39.

2ОАО «Синтез-Каучук»

Россия, Республика Башкортостан, 453II0 г. Стерлитамак, ул. Техническая, I4.

Тел./факс: +7 (3473) 20 40 09.

3Башкирский государственный университет Россия, Республика Башкортостан, 450074 г. Уфа, ул. Заки Валиди, 32.

Тел./факс: +7 (347) 229 96 35.

E-mail: promif@mail.ru

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

Ключевые слова: сополимеризация, бутадиен, стирол, кинетика, обратная задача, молекулярномассовое распределение.

Производство полимерных материалов занимает одно из ведущих мест в современной химической промышленности. Многие полимерные продукты получают методом радикальной сополиме-ризации, когда исходными соединениями являются несколько (2 и более) мономеров. В частности, этот механизм лежит в основе технологии производства эмульсионного бутадиен-стирольного каучука.

Процесс сополимеризации бутадиена со стиролом в производстве синтетического каучука проводится непрерывным способом в батарее последовательно соединенных между собой полимеризаторов при температуре 5 °С [1]. Для регулирования молекулярной массы применяется трет-додецилмеркаптан. Процесс сополимеризации останавливают при конверсии мономеров 70%.

Кинетический метод моделирования полимери-зационных процессов заключается в составлении и численном решении кинетических уравнений для концентрации всех типов частиц, участвующих в процессе. При построении модели сополимеризации бутадиена со стиролом в эмульсии будем считать, что реакционная способность активного центра на конце растущей цепи определяется лишь природой концевого звена [2]. Рассмотрим по четыре элементарных реакции роста и обрыва цепи с участием двух мономеров M1 и М2 и двух типов растущих цепей, отличающихся природой концевого звена. Тогда кинетическую схему сополимеризации бутадиена со стиролом можно описать в виде следующих стадий:

1. Инициирование. На стадии инициирования образуются свободные радикалы. Их образование происходит в два этапа. На первом этапе распадаются молекулы инициатора (I), содержащиеся в исходной смеси:

тим, что в литературе часто встречается запись, где в результате распада инициатора происходит образование не одного, а двух свободных радикалов (2Д). Однако здесь мы несколько упрощаем вид стадии, а саму константу скорости распада инициатора рассчитываем из результатов его исчерпывания в ходе процесса.

На втором этапе инициирования в результате взаимодействия первичного радикала и мономеров M1 и M2 образуются мономерные радикалы, отличающиеся природой концевого звена: д + М1 —Ь® p^

R + M2

->P2

0,1

где к, - константа скорости реакции инициирования. 2. Рост цепи:

P1 + M1-

n,m

P1 + M2-

n,m

P2 + M1-

n,m

P2 + M2-

P1 ,

n+1,m

> P2.........

где кр - константа скорости реакции роста цепи, Р - активная («растущая») цепь сополимера длиной п + т, содержащая п звеньев мономера М1 и т звеньев мономера М2.

3. Передача цепи в результате взаимодействия с регулятором 5:

Р1 + 5 > О + д,

п,т *~-п,т '

P2 + S-

-»Q + R,

*~-n,m '

-»R,

где

к, -

константа скорости распада инициатора,

R - инициирующий (первичный) радикал. Заме-

где к - константа скорости реакции передачи цепи, О - неактивная («мертвая») цепь сополимера длиной п + т, содержащая п звеньев мономера М1 и т звеньев мономера М2.

к

к

к

к

к

* автор, ответственный за переписку

ISSN 1998-4812

Вестник Башкирского университета. 2011. Т. 16. №2

337

4. Обрыв цепи в результате рекомбинации:

P1 + P1 — > Q

n,m r ,q

P1 + P2 ■

r ,q

P1 + P1

P1 +P1 ■

r ,q

Ж

и диспропорционирования:

P1 + P1 > Q + Q ,

n,m r,q Z-'njm z-'r,q5

P1 + P2 kd12

P2 + P1

r ,q

2 kd.

> Q„,m + Qr ,q ,

> Q„,m + Qr,q ,

Pn2,m + Pr,q n,m r,q

->Q + Q ,

n,m r,q

где к , к - константы скоростей реакций рекомбинации и диспропорционирования.

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

Упростим систему за счет использования принципа квазистационарности [3], заключающегося в постоянстве концентраций активных цепей. Тогда скорость образования радикалов (активных цепей) на стадии инициирования равна скорости их гибели на стадии обрыва:

СРХ сР2

п,т _ о п,т __ о

с ’ а

Относительно Р1 и Р 2 получим:

п,т п,т •'

к 11МР1 + к 12М 2Р' - к 11М1Р11 - к 21М1Р21 + к 15Р1 +

р 1Г п,т р 1^ п,т р 11 п—1,т р21 п—1,т ге§1 п,т

.. .../P1 + k 12M 2P' - к „М 'P1, - к М 'P2, ...

plr n,m p1„ n,m pll n—\,m p 21 n—1,m re

+ к,,C .P1 + к C 2P' + kd11C .P1 + k.12C 2P' _0,

r11 M1 n,m r12 M2 n,m d11 M1 n,m d12 M2 n,m

(I)

+ к „С 2Р2 + к „С 1Р2 + к.22С 2Р2 + к.21С 1Р2 = 0,

г 22 м2 г 21 м1 п,т й 22 м п,т й 21 м1 п,т ~

где С и С - концентрации активных цепей,

м1 м2

отличающиеся природой концевого звена:

С 1 = У Р1 ,

м1 п,т ’

п,т=0

(2)

Принципиальное отличие полимеризационных процессов от всех других химических процессов состоит в том, что вместо продукта с фиксированной молекулярной массой получают некоторое молекулярно-массовое распределение полимера (ММР). Статическая теория полимеризации исходит из возможности анализа ММР по особого рода средним молекулярным массам м , м , называемых

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

среднемассовой и среднечисленной. Для анализа вводятся понятия моментов активных цепей:

V? = Е Е (1™а + т^Ь ) кР1п,т ,

п=0 т=0

У

ЕЕ (nwa + mwb ) kK„,

ka

n _0 m _0

и неактивных цепей:

Ук _ ЕЕ (﹫ + ^ ) к °п,т ,

п_0 т_0

где к _ 0,1,К,¥.

Перепишем систему (1) в виде:

(кр11М 1 + кр12М 2 + кг«15 + (кг11 + кС11 )СМ1 + (кг12 + кС12 )^М2 )Рп1.т _

к „М'Р\ + к ,М'Р\ ,

р 11 п-1,т р21 п- 1,т“

(кр2^ 2 + кр21М 1 + кг,я 25 + (кг 21 + кС21 См1 + (к 22 + кС 22 См> У-.т _

(3)

k „M2P' , + k „M2P2 ,,

p12 n,m-1 p 22 n,m-1“

п, т > 1.

Получившуюся систему уравнений (3) решим с использованием производящих функций вида:

G( s, и) _ ЕЕ s

n+munwa +mwb

n,m ?

n_0 m_0

(4)

F (s, w) _ ЕЕ s

\ ’ \ 1 ^n+^un^a +mWbP 2

n,m '

n_0 m_0

n+m nwa +mwl

s u

Домножив левую и правую части (3) на

Ь и просуммировав получившиеся выражения по п, т, получим линейную систему уравнений относительно 0(5, и) и Р(5, и) :

(кр„м1 + кр12м2 + к^ + (кт + кЛ1 )См1 +(кГ12 + км )См2 )о =

(кр1^м-0 + крпм 2Р + Р,0 )т'-, (5)

(кр22м 2 + крпм 1 + кге^2^ + (кг21 + кй21 )См1 + (кг22 + кй22 )См2 )р =

(к^м ’о+кр22м2 р+а,0 )ш'ь ,

Выражая 0(5, и) и Р(^, и) из системы (5), получим формулы для расчета моментов молекулярно-массового распределения (ММР) активных цепей сополимера [4]:

й0(5,и)

Ук _-

У _

duk dkF (s, u)

du

s _ u _ 1,

s _ u _ 1,

(6)

где yM1 ,yM2 - моменты

к -го порядка активных цепей сополимера, отличающихся природой концевого звена; Wa, ^ - молекулярная масса мономеров М1 и М2 соответственно.

Для расчета средних молекулярных масс сополимера необходимо знание моментов до 2-го порядка включительно, тогда система дифференциальных уравнений относительно моментов ММР сополимера примет вид:

dI _ dt

dM1 dt

dM2 dt

— _-k ,SC , -k 2SC 2,

dt reg. M1 2 M2 ’

-ki I,

_-k„M 1C , - L,M‘C

-k p12M2C 1 - k p22M2C 2

p12 M1 p22 m2

(7)

У _ Aa(kr, к,, k„g ,yf,yf,[S 1[I ], [m 1 ^[M2 ^Cm,,Cm2)

і а п с7

Моменты

(а _ 0,1,2) рассчитываются по формулам (6)

же активных цепей yM , ум

аа

k

k

к

к 22M2P" + к „М1P* - к „їм2P‘ , - к „їм2P* , + к 2SP" +

p 22 n,m p21 n,m p12 n,m—1 p 22 n,m—1 reg2 n,m

n,m_0

Начальные условия для системы (7) имеет вид: I(0) = /0, м40) = м01, м2 (0) = м2, 5(0) = 50, (5)

у (0) = ум1 (0) = ум2 (0) = 0, Р _ 0,1,2.

Найденные значения молекулярно-массового распределения сополимера подставим в формулы для нахождения среднечисленной () и среднемассовой (м®) молекулярных масс:

M

(9)

MQ =Ж.

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

Ы_ 5.4 х 10—4 м%066. (10)

Кинетические параметры сополимеризации бутадиена со стиролом в эмульсии (при 5 °С) были взяты из работы [6], однако не все найденные значения соответствуют условиям и механизму процесса эмульсионной сополимеризации. Потому для уточнения некоторых кинетических параметров (констант скоростей реакций рекомбинации и дис-пропорционирования) проводится решение обратной задачи путем аппроксимации конверсионной зависимости.

Стратегия поиска кинетических параметров сополимеризации заключается в поиске минимума функции соответствия (11) экспериментально определяемых параметров и вычисленных на основе математической модели, описывающей процесс сополимеризации (7)-(8):

Н(г) = £ (к™ - К?асч. )2, г = к, кг), (11)

1=1

где кэксп' - значения конверсии, определенные в ходе эксперимента, красч' - расчетные значения конверсии, определяемые по формуле:

рм 1 + рм 2 К расч-_1 — ,_____,

1 ~ м1 + м2 ’

где рм 1 и рм2 - значения концентрация мономеров м1 и м2:

рм1 _ м1 • ехр((— крп • Си 1 — крг22 • См2 / Г ) • X, ),

рм2 _ м2 • ехр((— крп • См 1 /Г — кр22 • См2 )• X,)

Для уточнения кинетических параметров путем минимизации функционала (11) использовался метод Хука-Дживса.

Применяя неявный метод Адамса-Мултона 4го порядка для численного решения жесткой системы дифференциальных уравнений (7)-(8) и подставляя значения моментов в (9)-(10), определим

зависимости м®, м® и [^] от времени сополимеризации (рис.). Из приведенных зависимостей видно, что значения усредненных молекулярных масс, вычисленных по формулам (9), а также значение

характеристической вязкости, вычисленной по формуле (10), удовлетворительно описывают данные результатов эксперимента, проведенного в центральной заводской лаборатории ОАО «Синтез-Каучук» (г. Стерлитамак), что подтверждает адекватность используемой нами модели.

8 t, ч

8 t, ч

t, ч

Рис. Зависимость экспериментальных (точки) и расчетных (сплошная линия) значений характеристической

вязкости (а), среднечисленных (б) и среднемассовых (в)

молекулярных масс от времени.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

ЛИТЕРАТУРА

1. Кирпичников П. А., Береснев В. В., Попова Л. М. Альбом технологических схем основных производств промышленности синтетического каучука. Л.: Химия, 1986. 224 с.

2. Семчиков Ю. Д. Высокомолекулярные соединения. М.: Академия, 2003. 368 с.

3. Подвальный С. Л. Моделирование промышленных процессов полимеризации. М.: Химия, 1979. 256 с.

4. Ray W. H., Douglas T. I., Godsalve E. W. Molecular weight distribution in copolymer systems. II. Free radical copolymerization // Macromolecules. 1971. V. 4. P. 166-174.

5. Booth C., Beason L. R., Bailey J. T. Incremental Modification of styrene-butadiene rubber with fer/-dodecyl mercaptan // J. Appl. Polymer Sci. 1961. V. 13. P. 116-123.

6. Brandrup J., Immergut E. H., Grulke E. A. Polymer Handbook. Wiley, 1998. 2317 p.

Поступила в редакцию 19.01.2011 г. После доработки - 10.05.2011 г.

i Надоели баннеры? Вы всегда можете отключить рекламу.