CC BY
44
130 Вестник СамГУ — Естественнонаучная серия. 2009. № 6(72)
УДК 532.517.2:534.2
О ФОРМЕ ВЫХОДЯЩЕЙ СТРУИ ПРИ ВЫДАВЛИВАНИИ ВЯЗКОУПРУГОЙ ЖИДКОСТИ ИЗ ПЛОСКОГО ЩЕЛЕВОГО КАНАЛА
© 2009 Б.А. Снигерев]1 К.М. Алиев, Ф.Х. Тазюков2
Приведены результаты численного моделирования течения упруговязкой жидкости со свободной поверхностью на примере задачи об экструзии. В качестве реологической модели использована конститутивная модель упруговязкой жидкости ЕЕМЕ-Р, построенная на основе молекулярно-кинетической теории. Показано влияние на форму струи реологических свойств полимерной жидкости.
Ключевые слова: вязкоупругая жидкость, свободная поверхность, экструзия.
1. Предварительные сведения
Одним из направлений изучения течений со свободной поверхностью является исследование устойчивости поверхности экструдата к образованию волн. Образование регулярных искажений свободной поверхности полимерных жидкостей или эластической турбулентности является следствием многих факторов. Одним из важнейших факторов является проскальзывание экструдата на стенках канала при достижении критических значений сдвиговых напряжений в этой области. В таких задачах при определении формы свободной поверхности также необходимо учитывать влияние линии контакта трех фаз, являющейся фактором, способствующим развитию возмущений. Условия на линии контакта могут существенно влиять на движение жидкости. Несмотря на исключительную важность проблемы, количество публикаций в этой области недостаточно, а происходящие в окрестности межфазных границ процессы и их влияние на величину сдвиговых напряжений до сих пор недостаточно поняты. Экспериментальные
хСнигерев Борис Александрович (snigerev@mail.knc.ru), лаборатория моделирования технологических процессов института механики и машиностроения КазНЦ РАН, 420111, Россия, г. Казань, ул. Лобачевского, 2/31.
2Алиев Кулу Мехрали-оглы (alievmm@rambler.ru), Тазюков Фарук Хоснутдинович (tazyukov@mail.ru), кафедра теоретической механики и сопротивления материалов Казанского государственного технологического университета, 420025, Россия, г. Казань, ул. К. Маркса, 68.
исследования весьма затруднены из-за того, что все межфазные эффекты весьма чувствительны к примесям и физическому состоянию поверхности. Особое внимание ученых привлекает изучение причин появления так называемых ”спурт-эффектов” — нарушения стабильности истечения жидкости из капилляров и каналов [1—3]. В данной работе теоретически исследуется эластическое восстановление и разбухание струи экструдата. Возрастание диаметра или разбухание — общая особенность процесса экструзии при производстве синтетических волокон и пластмасс.
Целью настоящей работы является исследование течения, деформированного состояния вязкоупругой жидкости на выходе из экструдера в зависимости от расхода и релаксационных свойств полимера.
2. Постановка задачи
Исходной является система уравнений закона сохранения массы и импульса, для замыкания которой привлекается конститутивное уравнение FENE-P жидкости [3, 4] :
ду
в(+ V = -^Р + V- т, (2.1)
V- V = 0, (2.2)
где в — плотность жидкости, V — вектор скорости, Р — давление, т — де-виатор напряжения.
Реологическое уравнение используемой модели жидкости имеет вид
т = т 1 + т 2, т 2 = 2^2-0, (2.3)
т П1
т 1 _ т
1 - (гтЛ)/(гь2) 1 - 1/ь2
(2.4)
Л Л—Л _ I
1 - (ітЛ)/(3Ь2) + Х1і _ 1 - 1/Ь2' (2.5)
Здесь —:---верхняя конвективная производная от тензора конфигурации
Ы
Л; Пі — вязкость полимера; П2— вязкость растворителя; Л— время релаксации напряжений; ЬтЛ— след тензора Л; Ь— безразмерный параметр модели FENE-P жидкости, характеризующий степень растяжения данной макромолекулы. В данной модели Ь должно быть значительно больше 1 (Ь >> 1).
Верхняя конвективная производная от тензора конфигурации Л определяется как
—Л _ + V УЛ — Уь ■ Л — Л ■ (Уь)т, (2.6)
0І ОІ
а тензор скоростей деформации определяется соотношением О _ 1 (Уь + + Уьт).
В декартовых координатах система уравнений (2.1)—(2.5) имеет вид
Пі
Ах
+
1
Л \ 1 _ [Ахх+Ауу ] 1 — -1
\ ^ ь ь2
П1
А
УУ
+
Л \ 1 [Ахх +АУУ] 1 — -1 ’
V1 зь2 ь2 /
дАхх дАхх
+ и—-------+ V
ді
ді
+ и-
дх
дх
+ V
Тху Пі Л \1
дАхх \ ди
дУ ) -2 дХ
дАуу \ дv
дУ ) — 2Л~я~л дУ
А
хУ
зь2
ди
Ах
ду 1 [Ахх+Ауу] 1 _ і
у 1 - зь2 ь2
А
УУ
дх 1 _ [Ахх+АУУ] 1 — А
-1- ОГ2 ь2
+ ■
зь2
А
хУ
/ дАху + идАху + ^дАху\ ^ . ЛАууди
+ и^ + -ЛдХА + ЛА дУ- 1- [Ахз+Ауу]
Течение полимера в расчетной области (рис. 2.1) описывается уравнениями (2.1)—(2.5), на границе входа в расчетную область 51 задаются профили скорости и полные распределения напряжений V = «о (у),т = То (у), на выходе из области 52 задается установившийся однородный профиль скорости и напряжений =0, V = 0, =0, условия симметрии на границе 5з
и кинематическое и динамическое условие на свободной поверхности 5б :
V ■ п = 0, п ■ т ■ Ь = х/Н ■ п, (2.7)
где п,Ь — единичные векторы нормали и касательной на поверхности 5б, Н — радиус кривизны свободной поверхности, х — коэффициент поверхностного натяжения.
Аппроксимация уравнений (2.1)—(2.5) и вычисления проводятся методом конечных элементов (МКЭ) второго порядка на нерегулярных сетках, сгущающихся к зоне истечения полимера из насадки [5]. Поскольку искомые
1
1
хх
Л
1
Л
функции существенно меняются лишь в окрестности истечения полимера, применение мелких в этой окрестности и разреженных вне ее МКЭ-сеток играет существенную роль в экономии вычислительных затрат. Для расчетов строилась последовательность сгущающихся сеток 9-узловых четырехугольных элементов ( число узлов 2000, 8400). Для расчета напряжения использовались линейные четырехугольные элементы. Местоположение деформируемой свободной поверхности находится из аппроксимации кинематического условия (2.7), затем сетка конечных элементов вблизи нее перестраивается для получения решений уравнений (2.1)—(2.5), с помощью которых находится поле скоростей, давлений и напряжений на новом временном слое [6]. Стационарное положение формы выходящей струи находится методом установления эволюционной задачи с использованием традиционных для уравнений данного класса алгоритмов [7].
3. Результаты численного моделирования
Влияние релаксационных свойств экструдата, определяемых безразмерным критерием — числом Вайссенберга Шв (Ш в = ^и/Ь) на форму выходящей струи, показано на рис. 3.1, а. Здесь приведены кривые, соответствующие профилю свободной поверхности для чисел Шв = 0,1; 1, 0; 5, 0; 9, 0. Увеличение поперечного сечения струи по сравнению с сечением насадки, называемое эластичным восстановлением (ЭВ), составляет 1,6 для наибольшего числа Вайссенберга Шв = 9, 0. Влияние на эластичное восстановление параметра модели Ь показано на рис. 3.1, Ь. Численные эксперименты показывают, что с ростом времени релаксации (числа Ш в) высокоэластическое восстановление возрастает. Возрастание также наблюдается при увеличении параметра Ь — степени растяжения макромолекулы полимера, причем этот рост в этом случае носит асимптотический характер. Распределение давлений, нормальных и сдвиговых напряжений для Шв = 5,0 показано на рис. 3.2. На этих картинах видно, что как только полимерная жидкость проходит выходное сечение и попадает в зону свободного течения, давление и напряжения начинают немедленно релаксировать. Процесс релаксации происходит до тех пор, пока накопленная высокоэластичная деформация не уменьшится до значения, соответствующего эластической деформации полимера, находящегося в состоянии стационарного течения. Выделим две основные причины изменения сечения струи: перестройка профиля скоростей от параболического, соответствующего течению в канале, к прямоугольному, соответствующему движению струи как сплошного твердого стержня; релаксация высокоэластических деформаций растяжения, накопленных в материале струи за время ее прохождения по каналу.
Перестройка профиля скорости приводит к возникновению продольных деформаций растяжения или сжатия. Из требования выполнения закона сохранения количества движения следует, что струя, выходящая из канала,
должна сужаться. Этот эффект известен для низкомолекулярных жидкостей, не проявляющих высокоэластических свойств, и составляет 0,87. Для полимерных систем реализуется упругое разбухание, обусловленное высокоэластическими деформациями.
Ь 21—-------------------------------
0 -2 -4 -6
а
^ 3 * ‘
** ** — » * / ♦ / и . ^ 2 — 1 " * 8 !
1 0 -2 -4 -6 X
б
Рис. 3.1. Влияние релаксационных свойств экструдата на форму выходящей струи: а — зависимость ее формы от числа при Ь2 = 100, 1 — Ш в = 0,1;
2 — Ш в = 1, 0; 3 — Ш в = 3,0; 4 — Ш в = 9, 0; б — форма свободной поверхности экструдата при изменении параметра Ь при Шв = 1, 0;
1 — Ь =10,0; 2 — Ь = 100,0; 3 — Ь = 1000, 0
Рис. 3.2. Картина течения для Шв = 5, 0, Ь2 = 100 : а — изобары давления; б, в — линии уровня напряжений тхх,тху
Большой интерес представляет также распределение разности главных напряжений а 1 — а2 = у/N + 4т2у, где N1 = тхх — Туу. На рис. 3.3 представлено распределение разности главных напряжений для чисел Шв = = 1,0; 5,0. Видно, что возмущающее влияние выхода очень велико, оно приводит к значительному перераспределению скоростей и концентрации разности главных напряжений вблизи кромок выходного канала.
Рис. 3.3. Распределение разности главных напряжений о\ — :
а — Шв =1, 0, Ь2 = 100; б — Шв = 5,0, Ь2 = 100
Более детально данное изменение скорости представлено на рис. 3.4, где изображены распределения осевой составляющей скорости и по различным горизонтальным сечениям. Из них видно, что при Шв = 1,0, Ь2 = 100 профиль скорости становится почти однородным на расстоянии ширины выходного канала. Из рис. 3.5, где представлено распределение разности главных напряжений по этим же сечениям, видно, что напряжения полностью исчезают на удалении 1,5 ширины канала для данного числа Шв.
Отметим также зону максимальных напряжений вблизи верхней выходной кромки канала. Увеличение времени релаксации полимера приводит к большему разбуханию экструдата, и графики скорости и имеют некоторые особенности. Из рис. 3.4 видно, что по мере приближения полимера к выходу из головки происходит небольшое увеличение скорости в центральной части, это можно объяснить тем, что вблизи выхода образуется небольшая циркуляционная зона, сужающая зону выхода полимера из насадки. Наличие данной вихревой зоны приводит к небольшому замедлению скорости полимера вблизи стенки (рис. 3.4, сечение 4).
Рис. 3.4. Распределение продольной скорости для Ше = 5,0, 1? = 100 по продольным сечениям:
1 - у = 0, 0; 2 - у = 0, 5; 3 - у = 0, 75; 4 - у = 0, 95
N
1
20
10
-
4 ч, л ч"
\
—- 2 ! Л
\ 1
, , , ,
0
-4
X
Рис. 3.5. Распределение а\ — 72 для Ше = 5, 0, 12 = 100 по продольным сечениям:
1 - у = 0, 0; 2 - у = 0, 5; 3 - у = 0, 75; 4 - у = 0, 95
0
с.
Влияние вихревой зоны (рис. 3.5) видно и на распределении 71 — 72 по сечениям, где отметим зону уменьшения 71 — 72 перед выходом из головки вблизи стенки и более интенсивный пик при выходе. Протяженность зоны релаксации напряжений увеличивается до двух ширин канала. Изменение поверхности вытекающей из насадки струи связывают с началом неустойчивого течения. Более детальное изучение влияния релаксационных свойств полимера на форму выходящей струи и количественная информация о распределении напряжений при сравнении с экспериментальными данными помогут выявить критические значения скорости сдвига, при которых цилиндрическая форма струи экструдируемого полимера нарушает-
ся, и на ее поверхности появляются вмятины и утолщения [2]. Необходима более детальная информация об ориентации пристенных слоев расплава, вызывающих кристаллизацию (стеклование), и приводящая их к возникновению периодического проскальзывания. На основе этих знаний можно подобрать такие технологические режимы, которые помогут значительно уменьшить технологический дефект.
Выводы
Результаты численного моделирования истечения вязкоупругой жидкости из насадки на основе FENE-P модели жидкости показали, что для моделирования эффекта разбухания жидкости необходим выбор нелинейных вязкоупругих моделей, учитывающих первую и вторую разницу нормальных напряжений в сдвиговом течении. Существование нормальных напряжений также является одной из причин ’’разбухания” струи или Баррус-эффекта. Однако для более адекватного реологического описания данного явления необходимо учитывать также температурные факторы с помощью теории неравновесной термодинамики, поскольку известно, что эластическое восстановление и усадка резко возрастают при понижении температуры предшествующей деформации. Результаты данной работы получены с помощью более простой теории, однако они позволяют получить некоторые качественные оценки по форме свободной поверхности и распределению нормальных и сдвиговых напряжений в потоке экструдата. Полученная информация будет полезной при инженерном прогнозировании технологического процесса переработки полимеров. Явления эластического восстановления и разбухания струи экструдата очень сложны и многообразны, поэтому для адекватности математического описания необходимо также обращаться к эксперименту для более качественного прогнозирования данного явления.
Литература
[1] Чанг Дей Хан. Реология в процессах переработки полимеров. М.: Химия, 1979. 368 с.
[2] Торнер Р.В. Теоретические основы переработки полимеров. М.: Химия, 1977. 467 с.
[3] Bird R.B., Armstrong R.C., Hassager O. Dynamics of Polymeric Liquids. Vol. 1. Fluid Mechanics. 2nd ed. New York: John Wiley and Sons, 1987. 565 с.
[4] Bird R.B., Warner H.R., Evans D.C. Kinetic Theory and Rheology of Dumbell Suspensions with Brownian Motions // Adv. Polym. Sci. 1971. № 8. P. 1-90.
[5] High-order finite volume schemes for viscoelastic flows / B.A. Snigerev [et al.] // J. Comput. Phys. 2004. V. 199. P. 16-40.
[6] Течение упруговязкой жидкости со свободной поверхностью / Б.А. Снигерев [и др.] // Вестник Казанского техн. ун-та. 2007. № 1. С. 85-92.
[7] Коннор, Дж., Бреббиа К. Метод конечных элементов в механике жидкости. Л.: Судостроение, 1979. 264 с.
Поступила в редакцию 24/VT/2009; в окончательном варианте — 24/VT/2009.
ABOUT THE FORM OF OUTPUT JET AT SQUEEZING OUT THE VISCOELASTIC FLUID FROM PLANE GAP CHANNEL
© 2009 B.A. Snigerev3 K.M. Aliev, F.Kh. Tazyukov4
The results of the numerical modelling of flow of viscoelastic fluid with free surface are shown on the example of the task about extrusion.
As a rheological model, the constitutional model of a viscoelastic liquid FENE-P based on the molecular-kinetic theory is used. The influence of rheological property of polymeric fluid on the shape of a jet is shown.
Key words: viscoelastic fluid, free surface, extrusion.
Paper received 24/VT/2009. Paper accepted 24/V7/2009.
3Snigerev Boris Alexandrovich (snigerev@mail.knc.ru), Lab. of Modelling of Technological processes, Institute of Mechanics and Engineering of RAS, Kazan, 420111, Russia.
4Aliev Kulu Mekhrali-Ogli (alievmm@rambler.ru), Tazyukov Faruk Khosnutdinovich (Tazyukov@mail.ru), Dept. of Theoretical Mechanics and Strength of Materials, Kazan State Technological University, 420025, Russia.
