УДК 532.517.2
М. А. Кутузова, А. Г. Кутузов, Г. С. Кутузова,
Ф. Х. Тазюков
ТЕЧЕНИЕ ВЯЗКОУПРУГОЙ ЖИДКОСТИ FENE-P В ФОРМУЮЩЕЙ ГОЛОВКЕ ЭКСТРУДЕРА
Ключевые слова: вязкоупругая жидкость FENE-P, экструзионная головка, течение.
Рассматривается течение вязкоупругой жидкости FENE-P в канале, характерном для экструзионной головки. Исследуются условия течения полимеров вблизи твердой стенки и на линии трехфазного контакта.
Keywords: viscoelastic liquid FENE-P, extrusion head, flow.
Flow of a viscoelastic liquid FENE-P in the channel typical for extrusion head is considered. Conditions of polymers’ flow near the solid wall and on the line of the three-phase contact are researched.
Важным направлением изучения экструзии полимеров является исследование причин, приводящих к искажениям свободной поверхности экс-трудата. Образование регулярных искажений свободной поверхности полимерных жидкостей (эластической турбулентности) является следствием многих факторов [1-3]. Среди этих факторов особое место занимают условия течения экструдата вблизи твердых границ. При определении формы свободной поверхности также необходимо учитывать влияние линии контакта трех фаз, являющейся важнейшим фактором, способствующим пиковому росту напряжений.
Математическая постановка задачи. Рассматривается течение вязкоупругой жидкости FENE-P в плоском канале, характерном для формующей головки экструдера.
Основными уравнениями, описывающими медленное установившееся течение, являются следующие зависимости:
р (—) =-У Р +У-~ (1)
5 1
V -V = 0 (2)
где: р - плотность жидкости, V - вектор скорости, Р - давление, х - девиатор напряжения.
В соответствие с принципом расщепления напряжений:
~ (3)
Л 1
я
I
A
1 - (trA ) /(3 L2)
1 - (tr A ) /(3 L2)
v
+ я A =
1 - 1 /L2
1 - 1 / L2
(4)
(5)
х2 = 2^2 В . (6)
В выражениях (3)-(6): ^ 1 - вязкость полимерной составляющей жидкости, ^ - вязкость растворителя, X - характерное время релаксации, Ь2 -безразмерный параметр, характеризующий степень растяжения данной макромолекулы. Верхняя конвективная производная от тензора конфигурации определяется из выражения:
V 5 А ~ ~ ~ (7)
А = ------+ А - V V-А - А - (V? )т
51
Тензор скоростей деформации определяется соотношением:
D = —( Vv + Vv T ).
(8)
Запишем уравнения (1)-(6) в декартовой системе координат. Это приведет к следующей записи уравнений движения:
Эи Эи Эи ЭР Эт Хх Эт Ху э2и э2и (9)
Re(—+u— + v—)=—-+—— +--+р І- ,
дt дx Су дx дx Су дx2 ду2
ду dv dv дР дт xy дт Уу n.d2v d2v (10)
Re(—+u—+v— )=---+---+---+pl 1 v ’
' д t д x дУ д x д x дУ
2 +------т І
д x ду2
д u dv д x д у
■ 0 .
(11)
Уравнения эволюции тензора конфигурации в проекциях на оси координат запишутся в виде:
(12)
Х|ЗЧиЗАх+уЗАхЬя— а + 2-А-------+ ■
/^u ^+vCAxx W—A,x+ 2-A —4+^ +Л
[a dx ду ) dx^ ду1* 1 A.x+Ay 1-l
I---^-- 1 .2
dAyy dA dA . -.
я| ——yy+u—yy+v—yy |=2я— Ayy + 2я—A,, 1 ct dx 1 Уу ^ ■
dy
dAx,
dv
, dv ,
32
Ayy
(13)
dy
dx
.Ax+Ayy 1_____^
3L2 L2
я| —+u--------xy + v---xy |=2 я-Axx + 2я—A„
d x 1x1 Як УУ
д У
dv
Axy
(14)
dx
dy
■ _ Лсх+Лу Э1.2
Уравнения, связывающие компоненты полимерной части тензора напряжений с компонентами тензора конфигурации, запишутся в виде:
Л1
Л1
(15)
(1б)
(17)
В результате приведения уравнений (1)-(8) к безразмерному виду можно получить выражения чисел Рейнольдса, Вайссенберга, коэффициента ретардации и капиллярности:
т =
xx
я
т
УУ
я
Т =
xy
Ре =
Л :
риь X и
-—, We = —, р_
Л Ь Л1 + Л 2
Са =ни.
Граничные условия. На границе входа в расчетную область Б1 задаются профили скорости и
распределение напряжений: 1
и _ т(4 _ у 2 )>
4
V = 0,
т хх = 2We • (1 _р)(-|у)2, тху = We • (1 -Р)
Э и
Э у,
т уу _ 0-
На выходе из области течения Б2 задается установившийся однородный профиль скорости и
напряжений: у = о = о
’ 5 х ’
Э Т хх _ Э т ху _ Э т уу = 0, Р = 0
5 х 5 х 5 х
На границе Б3 задаются условия симметрии: у = о; х^ = 0 . На свободной поверхности
Б5 задаются кинематические и динамические условия:
V - п = 0 ; п - ~ - Т = а 0 / р - п , где п , Т - единичные векторы нормали и касательной на поверхности Б 5 ; а 0 - коэффициент поверхностного натяжения; р — радиус кривизны
свободной поверхности. На линии трехфазного контакта задается также контактный угол смачивания.
Решение уравнений. Аппроксимация уравнений (1)-(6) и вычисления проводились методом конечных элементов (МКЭ) второго порядка на нерегулярных сетках, сгущающихся к зоне истечения полимера из насадки. Так как искомые функции существенно меняются лишь в окрестности истечения полимера, то применение мелких в этой окрестности и разреженных вне ее МКЭ-сеток играет существенную роль в экономии вычислительных затрат. Для расчетов строилась последовательность сгущающихся сеток 9-узловых четырехугольных элементов (число узлов 2000, 8400). Для расчета напряжения использовались линейные четырехугольные элементы. Местоположение деформируемой свободной поверхности находилось из аппроксимации кинематического условия (9), затем сетка конечных элементов вблизи нее перестраивалась для получения решений уравнений (1)- (6), с помощью которых находилось поле скоростей, давлений и напряжений на новом временном слое. Стационарное положение формы выходящей струи находилось методом установления эволюционной задачи.
Результаты моделирования. Местоположение деформируемой свободной поверхности находилось из аппроксимации кинематического условия на сво-
бодной поверхности. На рис. 2 изображено распределение осевой составляющей скорости и по различным горизонтальным сечениям. Из этого рисунка видно, что при We _ 1.0, Ь2 _ 100.0 профиль скорости становится почти однородным на расстоянии ширины выходного канала.
Для анализа степени ориентации макромолекул полимера в потоке использован физический закон пропорциональности тензора напряжений тензору коэффициентов преломления или тензору оптической анизотропии.
Рис. 2 - Распределение скорости и для а) We = 1.0 и для Ь) We = 5.0 по продольным сечениям.
1 — у = 0.0, 2 — у = 0.5, 3 — у = 0.75, 4 —у = 0.95
В соответствие с этим законом, разность глав-ньк напряжений а — а2 =/N + 4хХу > где N1 =ххх —Хуу,
пропорциональна разности главных коэффициентов преломления: а 1 — а2 = С(п 1 — п2) . Таким образом, в области потока упруговязкой жидкости с предпочтительной ориентацией макромолекул разность скоростей прохождения обыкновенного и необыкновенного луча (эффект двойного лучепреломления) (п<1 — п2) количественно характеризует
степень ориентации макромолекул в потоке под действием развивающихся в области течения напряжений.
Следовательно, в соответствие с оптическим законом, разность главных напряжений с точностью до множителя С характеризует степень ориентации макромолекул в потоке.
На рис. 3 представлено распределение разности главных напряжений по этим же сечениям для Шв = 1.0 и Шв = 5.0 соответственно.
а) Ь)
Рис. 3 - Распределение а 1 _ а 2 для
We _ 1.0 (а) и We _ 5.0 (Ь) по продольным сечениям
1 — у _ 0.0, 2 — у _ 0.5, Э — у _ 0.75, 4 -у_ 0.95
а
0
а
Отметим также зону максимальных напряжений вблизи верхней выходной кромки канала. Увеличение времени релаксации полимера приводит к большему разбуханию экструдата и графики скорости и имеют некоторые особенности. Из рис.ЗЬ видно, что по мере приближения полимера к выходу из формующей головки происходит небольшое увеличение скорости на оси канала; это можно объяснить тем, что вблизи выхода образуется небольшая область замедленного течения жидкости, сужающая сечение выхода полимера из насадки. Наличие данной зоны приводит к некоторому замедлению скорости полимера вблизи стенки. Влияние «застойной» зоны сказывается и на распределении степени ориентации макромолекул полимера в потоке вблизи выходного сечения насадки. Из полученных результатов следует, что в застойной зоне, расположенной вблизи выходного сечения, степень ориентации макромолекул заметно уменьшается. Данное обстоятельство приводит к тому, что механические свойства экструдата резко меняются по сечению, что крайне нежелательно.
Влияние зоны замедленного течения видно и на распределении а1 — а2 по сечениям, представленного на рис.3, где отметим зону уменьшения а1 — а 2 перед выходом из головки вблизи стенки и более интенсивный пик при выходе. На рис.4 представлено распределение разности главных напряжений для различных значений числа Вайссенберга.
Видно, что возмущающее влияние выхода очень велико: оно приводит к значительному перераспределению концентрации разности главных напряжений вблизи кромок выходного сечения плоского канала формующего инструмента.
Для предупреждения образования застойной области и уменьшения пиков напряжений можно допустить некоторое постоянное проскальзывание экструдата по внутренней поверхности формующего инструмента. В работе показано, что скольжение экструдата вблизи выходного сечения уменьшает первую разность главных напряжений а1 — а2, что свидетельствует об уменьшении степени ориентации макромолекул, а также снижает пики напряже-
ний в окрестности линии трехфазного контакта.
bj
4x;j
Рис. 4 - Распределение разности главных напряжений а1 -а2 a) we = 1.0, L2 = 100.0 , b) We = 5.0, L2 = 100.0.
Однако при этом образуется дополнительная особая точка, расположенная вблизи твердой стенки и разделяющая области прилипания и скольжения.
В этой новой особой точке образуются пиковые напряжения, способные стать дополнительным фактором, приводящим к разрушению экструдата. Поэтому имеет смысл допускать проскальзывание по всей внутренней поверхности формующего инструмента. В этом случае на линии трехфазного контакта пики напряжений заметно снижаются, а дополнительная особая точка не образуется.
Литература
1. Migler K., Son Y., Qiao F., Flynn K. Extensional deformation, cohesive failure and boundary conditions during sharkskin melt fracture. // Journal of Rheol.- 2001, v.31, p.6134.
2. Miller E., Rothstein J. Control of the sharkskin instabilitity in the extrusion of polymer melts using induced temperature gradients. // Rheol. Acta.- 2004, v.44, p.160-173.
3. Снигерев Б.А. Течение упруговязкой жидкости со свободной поверхностью / Б.А. Снигерев, Ф.Х. Тазюков, А.Г. Кутузов, А. Аль-Раваш // Вестник Казанского технологического университета. - Казань: КГТУ, 2007.- № 1. - С.85-93.
© Г. С. Кутузова - канд. техн. наук, доц. каф. химии НХТИ КНИТУ; А. Г. Кутузов - д-р. техн. наук, проф. каф. автоматизации технологических процессов и производств НХТИ КНИТУ, [email protected]; М. А. Кутузова - канд. техн. наук, доцент той же кафедры; Ф. Х. Тазюков - д-р техн. наук, ст. науч. сотр. каф. теоретической механики и сопротивления материалов КНИТУ.