CC BY
105
УДК 539.219.3
РАСЧЕТ ЭФФЕКТИВНОГО КОЭФФИЦИЕНТА ДИФФУЗИИ ПО СТРУКТУРНЫМ ХАРАКТЕРИСТИКАМ ПОЛИМЕРОВ
С.П. Рудобашта1, А.М. Климов2, Ю.А. Тепляков3, В.М. Нечаев4
Кафедры: «Теплотехника и энергообеспечение предприятий»,
ФГБОУ ВПО «Московский государственный агроинженерный университет имени В.П. Горячкина» (1); «Технологии пищевых продуктов» (2), «Прикладная геометрия и компьютерная графика» (3), «Технологические процессы и аппараты» (4), ФГБОУ ВПО «ТГТУ»;
Представлена членом редколлегии профессором С.В. Мищенко
Ключевые слова и фразы: нитроцеллюлоза; пористость; степень кристалличности; структура; сушка; экстрагирование; эффективный коэффициент диффузии.
Аннотация: Исследованы структурные характеристики нитроцеллюлозы и внутренний массоперенос при сушке от этанола. Предложена зависимость для расчета эффективного коэффициента диффузии с учетом структурных характеристик нитроцеллюлозы.
Одним из проблемных вопросов расчетной практики массообменных процессов является знание физических закономерностей массопереноса и информация о кинетических коэффициентах. Разработаны различные математические модели [1-6], описывающие влияние структуры и массосодержания (концентрации) распределяемого вещества на эффективный коэффициент диффузии в твердой фазе. Тем не менее имеются еще значительные трудности в расчете последнего, что вызывает необходимость его экспериментального определения.
Подавляющее большинство методов расчета эффективного коэффициента диффузии разработано и используется для капиллярно-пористых материалов, если известны значения коэффициента свободной диффузии в газовой среде или растворе и ряд характеристик пористой структуры материала: общая пористость ем; размер пор гк; интегральная и дифференциальная функции распределения пор по радиусам; коэффициент извилистости пор Ти; коэффициент формы (гофрировки) пор ц г. Наши исследования [7, 8] и литературные данные других авторов [4, 5] показали, что такие факторы, как Ти и цг, могут быть рассчитаны из дифференциальной функции распределения пор по размерам и по основному параметру материала - общей пористости ем.
Массоперенос распределяемого вещества в непористых (в кинетическом и сорбционном отношениях) полимерных материалах отличается двумя особенностями: 1) материал, как правило, состоит из аморфной и кристаллической фаз, резко различающихся по своим диффузионным свойствам. Диффузия молекул осуществляется, в основном, через аморфную часть полимера, тогда как диффу-
зионная проводимость кристаллитов пренебрежимо мала; 2) миграция распределяемого компонента осуществляется путем диффузии его молекул не по жестко фиксированным порам, а через «окна» в матрице полимера, которые хаотически появляются и исчезают вследствие теплового движения сегментов макромолекул, в соответствии с законами вероятности. Температурная зависимость коэффициентов диффузии в непористых полимерных материалах хорошо аппроксимируется уравнением Аррениуса
D э u = D даu exP(—ED, u/R*T), (1)
где Dэ u - эффективный коэффициент диффузии в материале, м2/с; D,, даu - эффективный коэффициент диффузии в материале формально при T ^ да, м2/с; Ed u -
энергия активации диффузии в полимере, Дж/кмоль; R* = 8,314 Дж/(моль • К) -универсальная газовая постоянная.
Проведенные нами исследования [9] показывают, что основным параметром структуры непористых полимеров, оказывающим влияние на диффузию распределяемого вещества в нем, является объемная доля аморфной фазы еам которую можно оценить как еам = 1 — х, где х - степень кристалличности полимера. Эту величину можно рассматривать в качестве некоторого, хотя и неполного, аналога общей пористости в капиллярно-пористых материалах ем. В соответствии с этим эффективный коэффициент диффузии низкомолекулярного вещества в частично закристаллизованном полимере при его массосодержании u можно рассчитать по соотношению
Dэu ~ Dэu амеам ~ Dэu ам (1 — х), (2)
где Dэ u ам - эффективный коэффициент диффузии в аморфном полимере с нуле-
вой степенью кристалличности, м2/с.
Для учета влияния степени кристалличности х и массосодержания полимера на параметры Ed, u и Dx u в уравнении (1) предложены эмпирические зависимости [10], полученные при обобщении опытных данных по коэффициентам диффузии воды в полимерах:
Ed,u = 7,8•Ю-4 [1 — 0,475 (1 — х)](1 — 0,16 u/umaxс); (3)
Dдаu = exP[—11 (1 — Х) — 0,92 u/umaxс], (4)
где umax с - максимальное сорбционное массосодержание полимера, кг/кг твердой фазы.
Диффузия молекул какого-либо вещества в полимере осуществляется путем их перескока из одного положения локального равновесия в другое через появляющиеся «окна» в матрице полимера. Для осуществления этого акта вблизи молекулы должно появиться «окно» площадью f > 0,785 d диф.в, где dfl^.H - диаметр молекулы диффундирующего вещества. В соответствии с этими представлениями эффективный коэффициент диффузии, при прочих равных условиях, должен быть обратно пропорционален диаметру молекул диффундирующего вещест-
2 /2
ва и, следовательно, можно записать: D.3uв = аDэh2o u, где а = dH Q/ dfl^.H ; Dэu h2o - эффективный коэффициент диффузии молекул воды, м2/с; dH2o -диаметр молекул воды, м.
Коэффициенты D^, u в рассчитывали по уравнению (5), в котором величины D^, да u и Ed u определяли по соотношениям (3) и (4). Поскольку они получены для диффузии молекул воды, то для учета размеров других диффундирующих молекул введен корректирующий множитель а = dН о jdj^ в, который при сушке
от воды смачивающихся полимеров, находящихся в высокоэластическом состоянии, равен единице. Было найдено, что при экстрагировании низкомолекулярного соединения из поликапроамида а = 0,54.
С учетом корректирующего множителя а перепишем уравнение (1) в виде
Du = ^эда u exP (^—RTf) (5)
В качестве примера на рис. 1 приведено сравнение расчетных и экспериментальных концентрационно-температурных зависимостей эффективных коэффициентов при сушке от воды полиамида 66 (х = 0,7), поликапроамида (х = 0,35) и экстрагировании из поликапроамида.
В данной работе приведены результаты изучения закономерностей внутреннего массопереноса в коллоидных капиллярно-пористых полимерных материалах при сушке их от жидкостей, являющихся активными растворителями по отношению к твердой фазе. Типичным примером может служить сушка нитроцеллюлозы от этанола. На первом этапе с целью идентификации структуры материала и его состояния изучили структурно-сорбционные и термомеханические свойства. Исследование структуры сорбционным методом (по поглощению паров воды, не абсорбируемых матрицей полимера) и методом ртутной порометрии выявило наличие в материале пор. Пористость нитроцеллюлозы по сорбционным данным составляет 0,031 м3/м3, а по ртутно-порометрическим, в процессе которых происходит проникновение ртути в замкнутые поры - 0,036 м3/м3, то есть объем замкнутых пор может быть оценен в 0,005 м3/м3, что составляет 14 % от общей пористости. Определяющий радиус транспортных пор по данным сорбционных измерений равен ~ 30 -10 10 м, а их радиус в устье (у поверхности тела) согласно ртутно-порометрическим измерениям равен ~ 55 • 10 10 м. Емкость адсорбционного монослоя по воде —2
um = 0,75-10 кг/кг сорбента. Удельная поверхность открытых пор, измеренная по низкотемпературной ад-
D3 u -1010, м2/с
u • 102, кг/кг
Рис. 1. Зависимости Бэи = /(и) при сушке от воды (1, 5) и при экстрагировании (2):
1 - полиамид 66 (Т = 373 К); 2 - поликапроамид (Т = 412 К); 3 - поликапроамид (Т = 423 К);
----------экспериментальные данные;
----------расчет по формулам (3) - (5)
сорбции азота с расчетом по методу БЭТ, равна 26,6 м2/г. Из чего можно сделать заключение, что нитроцеллюлоза может быть отнесена к классу коллоидных капиллярно-пористых тел по классификации [6] с микрокапиллярами (определяю-
_7
щий размер пор гопр р 10 м).
Известно, что характер диффузии низкомолекулярных веществ в полимере во многом зависит от того, в каком физическом состоянии последний находится -стеклообразном, высокоэластическом или вязкотекучем [11]. Ввиду отсутствия этих данных для нитроцеллюлозы были проведены соответствующие измерения термомеханическим методом [12].
На рисунке 2 приведены зависимости температур стеклования и деструкции нитроцеллюлозы от массовой доли этанола в материале. С увеличением концентрации спирта обе температуры снижаются, что объясняется ослаблением межмо-лекулярных связей полимера. Кривая 1 разграничивает области застеклованного (Г < Гст) и высокоэластического (Г > Гст) состояний полимера. Кривая 2 определяет верхнюю допустимую температуру нагрева материала в процессе сушки, при ее превышении в структуре полимера возникают пузыри пара. Исследования показали, что температурным условиям сушки нитроцеллюлозы в производстве соответствует высокоэластическое и стеклообразное состояние этого полимера.
Проведены исследования структуры этого материала микроскопическим методом. Обнаружено двухфазное строение полимера (кристаллическая и аморфная фазы) и подтверждено наличие замкнутых микропор в листовой нитроцеллюлозе. Расчеты кристаллической фазы полимера планиметрическим методом, по микрофотографии структуры нитроцеллюлозы, увеличенной в 600 раз, показали, что степень кристалличности изменяется в пределах х= 0,4 к 0,78, при изменении массосодержания от 0,11 до 0,01 (кг / кг твердой фазы).
Для определения эффективных коэффициентов диффузии этанола в нитроцеллюлозе снимали поля массосодержаний при сушке на пластинах толщиной 3 и 4 мм при температурах 303 и 313 К. Эксперимент заключался в сушке пластин одинаковой толщины в воздушном термостате при заданной температуре. В нужный момент времени одну из пластин извлекали, в центре ее на штампе вырубали три таблетки, которые разрезали на восемь слоев параллельно плоскости листа и каждый слой анализировали на содержание этанола. Сушка листовой нитроцеллюлозы происходит крайне медленно ввиду большой величины внутридиффузи-
онного сопротивления и низких температур процесса.
На рисунке 3 приведены поля массосодержаний для одного из опытов. Как видно, они имеют обычный для полей массосодержаний при сушке вид, необычным является лишь медленное - в течение всего времени сушки - понижение поверхностного массосодержания в образце, тогда как в соответствии с условием Ыт ^ Ж
оно должно сразу после начала процесса принять равновесное значение и сохранять его в ходе всего процесса сушки, проводимой при постоянных параметрах внешней среды.
Медленное установление десорб-ционного равновесия на поверхности
и-102, кг/кг
Рис. 2. Зависимость температуры стеклования (1) и деструкции (2) от массосодержания этанола в нитроцеллюлозе
раздела фаз следует объяснить напряженным состоянием полимера: под действием набухшего центрального ядра его поверхностные слои испытывают растягивающие усилия, что препятствует десорбции растворенных в матрице полимера молекул этанола, возникают диффузионные затруднения внутри полимера. Из полей массосодержаний этанола при сушке нитроцеллюлозы рассчитали эффективные коэффициенты диффузии, на рис. 4 приведены опытные значения Дэ и = / (и). Было найдено, что в области температур 303 и 313 К и массосодержа-
ний 0,01...0,11 кг/кг твердой фазы температурная зависимость коэффициента диффузии этанола в нитроцеллюлозе подчиняется закону Аррениуса (1).
Коллоидные капиллярно-пористые полимерные материалы по своим структурным свойствам занимают промежуточное положение между непористыми и капиллярно-пористыми материалами, что предопределяет промежуточность их диффузионных свойств. В общем случае эффективный коэффициент диффузии для конкретной системы является функцией концентрации распределяемого вещества и, температуры Т, структуры материала и может быть представлен в виде
Оэи = I[и, Т, ем, х(и)]. (6)
Представляет интерес проверить применимость зависимостей (3) - (5) для расчета Дэ и при диффузии этанола в нитроцеллюлозе, находящейся в высокоэластическом и стеклообразном состояниях, по двум определяемым параметрам: сте-
и-102, кг/кг
х-103, м
Рис. 3. Поля массосодержания при сушке листовой нитроцеллюлозы от этанола:
(2 5 = 4-10-3 м ) Т = 313 К т = 10-4 с: 1 - 4,32; 2 - 8,64; 3 - 17,28; 4 - 25,92; 5 - 34,56; 6 - 43,20;
7 - 60,48; 8 -86,40; 9 - 161,28
Бэ и -1011, м2/с
и-102, кг/кг
Рис. 4. Зависимости Бэ и = / (и) при Т = 303 К (1) и Т = 313 К (2):
сплошные линии - расчет по формулам (3), (4), (7)
пени кристалличности полимера х, максимальному сорбционному массосодержа-нию итахс с корректирующим множителем а в уравнении (5).
Из сопоставления эффективных коэффициентов диффузии, рассчитанных по уравнениям (3) - (5) и экспериментально полученных из полей массосодержаний (см. рис. 4) было найдено значение а = 7,0 и аппроксимирующая формула для диффузии этанола в нитроцеллюлозе, таким образом, может быть записана в виде
Средняя относительная погрешность вычислений Дэи по уравнениям (3),
(4), (7) составляет ± 20 %. В связи с тем что общая пористость нитроцеллюлозы незначительна, то ее не выделяли в уравнении (6) в качестве отдельного параметра.
Корректирующий множитель а в уравнении (7), в данном случае учитывает не только отличие диаметра диффундирующих молекул этанола от диаметра молекул воды, но и наличие пор в материале, облегчающих диффузию.
Сравнение экспериментальных и рассчитанных эффективных коэффициентов диффузии (см. рис. 4) позволяет распространить расчетные зависимости (3), (4), (7) на застеклованный полимер (см. рис. 2) (и < 0,075 кг/кг при Т = 303 К; и < 0,04 кг/кг при Т = 313 К): как следует из рис. 4, уменьшение массосодержа-ния ниже значений ист при указанных температурах не приводит к изменению зависимостей Дэи = I(и). Это объясняется сравнительно небольшой разницей температур (Тст - Т) < 12 К.
Уравнения (3), (4), (7) имеют вполне определенный физический смысл, носят полуэмпирический характер и требуют минимальной информации о структуре полимера и его сорбционных свойствах. Вычисленные по этим уравнениям концентрационно-температурные зависимости коэффициентов диффузии удовлетворительно аппроксимируют экспериментальные данные.
1. Проанализировано влияние структуры твердой фазы на эффективный коэффициент диффузии, показана возможность его расчета по основным параметрам: общей пористости ем для капиллярно-пористых материалов и степени кристалличности х для непористых полимерных материалов соответственно.
2. Определены структурные характеристики нитроцеллюлозы. Установлена величина общей пористости, определяющая радиус пор, удельная поверхность открытых пор. Выявлено, что структурные характеристики нитроцеллюлозы изменяются в зависимости от массосодержания этанола. Полученные данные по структуре этого полимера позволяют идентифицировать полимер как коллоидный капиллярно-пористый материал и определять температурный режим сушки.
3. Экспериментально исследованы поля массосодержаний этанола в нитроцеллюлозе при сушке, из которых определены концентрационные зависимости эффективного коэффициента диффузии при различных температурах.
4. Предложена зависимость, позволяющая прогнозировать значение эффективного коэффициента диффузии при сушке нитроцеллюлозы от активного растворителя - этанола по двум параметрам структуры полимера: степени кристалличности и максимальному сорбционному массосодержанию полимера.
(7)
Заключение
1. Хейфец, Л.И. Многофазные процессы в пористых средах / Л.И. Хейфец, А.В. Неймарк. - М. : Химия, 1982. - 319 с.
2. Аксельруд, Г.А. Массообмен в системе твердое тело - жидкость / Г.А. Ак-сельруд. - Львов : Изд. Львов. ун-та, 1970. - 188 с.
3. Рейтлингер, С.А. Проницаемость полимерных материалов / С.А. Рейт-лингер. - М. : Химия, 1974. - 272 с.
4. Саттерфилд, Ч.Н. Массопередача в гетерогенном катализе / Ч.Н. Саттерфилд ; пер. с англ. А.Р. Брун-Цехового. - М. : Химия, 1976. - 240 с.
5. Кришер, О. Научные основы техники сушки / О. Кришер ; под ред. А.С. Гинзбурга. - М. : Изд-во иностр. лит., 1961. - 539 с.
6. Рудобашта, С.П. Массопренос в системах с твердой фазой / С.П. Рудо-башта. - М. : Химия, 1980. - 248 с.
7. Rudobasta, S.P. Computation of the Pore Sinuosity and Shape Factors in Diffusion Through Homogeneously Porous Materials / S.P. Rudobashta, Yu.A. Teplyakov // Fluid Mechanics Soviet Research. - 1987. - Vol. 16, No. 5. -P. 101-108.
8. Внутренний массоперенос веществ при экстрагировании из однороднопористых материалов / Ю.А. Тепляков [и др.] // Вестн. Тамб. гос. техн. ун-та. -2010. - Т. 16, № 2. - С. 319-325.
9. Внутренний массоперенос веществ при экстрагировании и сушке / А.М. Климов [и др.] // Изв. вузов. Химия и хим. технология. - 2011. - Т. 54, вып. 5. - С. 117-120.
10. Рудобашта, С.П. Диффузия в химико-технологических процессах // С.П. Рудобашта, Э.М. Карташов. - 2-е изд., перераб. и доп. - М. : Колос С, 2010. -478 с.
11. Роджерс, К. Проблемы физики и химии твердого состояния органических соединений : пер. с англ. / К. Роджерс. - М. : Мир, 1968. - 475 с.
12. Степанов, В.П. Деформация и разрушение полимеров / В.П. Степанов // Механика полимеров. - 1975. - № 1. - С. 95-100.
Calculation of Effective Diffusion Coefficient on Structural Properties of Polymers
S.P. Rudobashta1, A.M Klimov2, Yu.A. Teplyakov3, V.M. Nechaev4
Departments: “Heat and Power Supply of Enterprises”,
Moscow State Agroengineering University named after V.P. Goryachkin (1);
“Technology of Food Products” (2), “Applied Geometry and Computer Graphics” (3), “Technological Processes and Apparatuses” (4), TSTU; [email protected]
Key words and phrases: degree of crystallinity; drying; effective diffusion coefficient; extraction; nitrocellulose; porosity; structure.
Abstract: The paper studies structural characteristics of nitrocellulose and internal mass transfer in drying from ethanol. The dependence to calculate the effective diffusion coefficient given the structural characteristics of nitrocellulose has been proposed.
Berechnung des wirksamen Diffusionskoeffizienten nach den strukturellen Charakteristiken der Polymere
Zusammenfassung: Es sind die strukturellen Charakteristiken des Cellulosenitrats und die innere Massenübertragung beim Trocknen vom Äthylalkohol untersucht. Es ist die Abhängigkeit für die Berechnung des wirksamen Diffusionskoeffizienten unter Berücksichtigung der strukturellen Charakteristiken des Cellulosenitrats vorgeschlagen.
Calcul du coefficient efficace de la diffusion d’après les caractéristiques de structure des polymères
Résumé: Sont étudiées les caractéristiques de structure de la nitrocellulose et le transfert de masse lors du séchage à partir de l’alcool éthylique. Est proposée une dépendance pour le calcul du coefficient efficace de la diffusion compe tenu des caractéristiques de structure de la nitrocellulose.
Авторы: Рудобашта Станислав Павлович - доктор технических наук, профессор, заведующий кафедрой «Теплотехника и энергообеспечение предприятий», ФГБОУ ВПО «Московский государственный агроинженерный университет им. В.П. Горячкина», г. Москва; Климов Анатолий Михайлович - кандидат технических наук, профессор кафедры «Технологии пищевых продуктов»; Тепляков Юрий Александрович - кандидат технических наук, доцент кафедры «Прикладная геометрия и компьютерная графика»; Нечаев Василий Михайлович -кандидат технических наук, доцент кафедры «Технологические процессы и аппараты», ФГБОУ ВПО «ТГТУ».
Рецензент: Килимник Александр Борисович - доктор химических наук, профессор, заведующий кафедрой «Химия», ФГБОУ ВПО «ТГТУ».
