Математическое моделирование процесса стратифицированного течения расплавов полимеров в осесимметричной постановке

ВЕСТНИК ПНИПУ

2014 Механика № 2

УДК 621.315

М.В. Бачурина, А.В. Казаков, Н.М. Труфанова

Пермский национальный исследовательский политехнический университет,

Пермь, Россия

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССА СТРАТИФИЦИРОВАННОГО ТЕЧЕНИЯ РАСПЛАВОВ ПОЛИМЕРОВ В ОСЕСИММЕТРИЧНОЙ ПОСТАНОВКЕ

Данная статья посвящена определению наиболее рациональной геометрии кабельной головки, позволяющей эффективно распределять потоки расплавов полимеров в каналах; оценке стабильности границ раздела потоков материалов с различными физико-реологическими свойствами для технологических диапазонов режимов переработки; расчету и построению полей скоростей, давлений, температур внутри каналов кабельной головки; экспериментальному определению зависимостей толщин накладываемых слоев изоляции и полупроводящих материалов от линейной скорости протяжки жилы и расходов для каждого канала. При моделировании процессов использовалась кабельная головка для совместного наложения трех слоев полимерного покрытия (слой полупроводящего экрана по жиле, изоляция, внешний полупроводящий экран), применяемая в производстве современных электрических кабелей среднего и высокого напряжения. Для анализа процессов тепломассопереноса в условиях стратифицированного течения в каналах кабельной головки реальные физические процессы были заменены математической моделью, представляющей собой систему нелинейных дифференциальных уравнений, отражающих основные законы сохранения. Система была дополнена граничными условиями и физико-реологическими свойствами перерабатываемых материалов. С целью упрощения модели предложен ряд допущений, позволивших, в частности, перейти к осесимметричной постановке задачи. Для решения сформулированной математической модели был применен численный метод, а именно метод конечных элементов, программно реализованный посредством комплекса АпБуБ. На основе полученных результатов была разработана более эффективная геометрия кабельной головки, устраняющая эффекты завихрения потоков полимера; получены распределения скоростей, давлений и температур в каналах кабельной головки; представлены возможные перегревы материала внутри каналов; оценено влияние некоторых технологических параметров процесса наложения многослойной полимерной изоляции на толщины изолирующих слоев.

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

M.V. Bachurina, A.V. Kazakov, N.M. Trufanova

Perm National Research Polytechnic University, Perm, Russian Federation

MATHEMATICAL MODELLING OF STRATIFIED FLOW OF POLYMER MELTS IN AN AXISYMMETRIC FORMULATION

This study is about definition of rational geometry of the cable die, which effectively distributes polymer melt flows in channels; estimation of stability boundary flows of materials with different physical and rheological properties for the technological range of modes of processing; calculation and construction of the velocity fields, pressures, temperatures inside the channels of cable die; experimental determination of the dependence of the thickness of imposed layers of insulation and semiconducting materials from the linear velocity of pulling wires and expenses for each channel. In the modeling of the processes was used the cable die for joint overlay three layers of the polymer coating (layer semicon-ductive conductor screen, insulation, outer semiconductive screen) used in the production of modern electric cables of medium and high voltage. To analyze the processes of heat and mass in terms of stratified flows in channels of cable die real physical processes have been replaced by a mathematical model, which is a system of nonlinear differential equations, reflecting the basic conservation laws. The system was supplemented of boundary conditions and physical and rheological properties of the materials processed. In order to simplify the model offers a number of assumptions that allowed us, in particular, go to the axisymmetric formulation of the problem. To solve the formulated mathematical model the numerical method was applied, namely the finite element method, implemented through the Ansys set of software. Basing on the received results was made more effective geometry of the cable die, eliminating the effects of twist the polymer streams is developed; velocity, pressure and temperature distribution in the channels of the cable die were obtained; possible overheating of the material inside the channels is presented; the effect of certain parameters of overlay multilayer polymeric insulation process on layer thickness was assessed.

Keywords: stratified flow, mathematical modeling, extrusion, rheology, polymers, abnormally viscous liquids, cable die, numerical solution, interface, axisymmetric, laminar.

На сегодняшний день продолжают оставаться востребованными исследования, посвященные изучению процессов многокомпонентных стратифицированных течений, сопровождающих и определяющих условия многих технологических процессов. Особое внимание следует уделить многослойным течениям жидкостей с аномально-вязкими свойствами, так как для подобных материалов режим течения будет определяться значительным числом влияющих факторов. Осложняет задачу взаимодействие нескольких потоков материалов с различными физико-реологическими свойствами. При этом имеются участки совместного течения: исследования показывают [1], что стабильность границы раздела оказывает существенное влияние на результат всего процесса и качество получаемого изделия. Этот процесс весьма затруднителен для непосредственного наблюдения, так как происходит внутри инструмента. Проследить распределение материалов можно,

например, посредством моментального охлаждения инструмента с последующим разбором и произведением продольного среза. Далее с помощью специальных приемов можно получить картину распределения отдельных потоков [2-6]. Однако такой путь весьма трудоемок и подходит к каналам простой формы. Поэтому для анализа процессов течения полимеров в каналах сложной формы было принято решение о численном моделировании [7].

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

Геометрия каналов кабельной головки представлена на рис. 1. Различными цветами обозначены расплавы полимеров с отличающимися физико-реологическими свойствами. Высота канала на выходе составляет 6,9 мм, длина горизонтального участка соприкосновения с подвижной стенкой (жилой) 27 мм, толщина канала 1 в зоне примыкания к подвижной жиле 0,7 мм.

Такую геометрию удобно рассматривать в цилиндрической системе координат (рис. 2).

Здесь различными цветами отмечены слои материалов с различными свойствами, при этом хорошо заметны участки совместного течения. Было получено расположение границ раздела потоков на выходе каналов, по которым можно судить о количестве материала, накладываемого на готовое кабельное изделие. Кроме того, в ходе решения

Рис. 1. Конфигурация кабельной головки

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

Рис. 2. Каналы истечения полимеров в цилиндрической системе координат

Математическое описание процессов тепломассообмена расплавов полимера при стратифицированном течении основано на законах сохранения [8].

Для упрощения задачи был сделан ряд допущений:

- процесс стационарный;

- среда несжимаемая, без упругих свойств [9];

- массовые силы значительно меньше поверхностных;

- течение осесимметричное;

- на границах канала реализуются условия прилипания и непроницаемости;

- теплофизические характеристики постоянны.

С учетом допущений система дифференциальных уравнений [8, 10] для каждого из слоев потока в скоростях имеет следующий вид:

1 д(гуг) + Эу^_ о г дг дС,

А (2ц э ди)+М. ди

дг I дг I г дг

д_ дг

ц э

(ду^ ЗуЛ

дг +дС

+

ц э

2 г д -—ц э уг

г дС,

ц э

+-

(ду^ ду.) дг + дС

(

дУг 1

2ц г

дг др

р'С'

(

дт дт) 1 д

уг эТ+у<

дС.

дГ

, гХ— II „

г дг ^ дг I дС, ^

дС I дС

д (, , дТЛ

+ Ф,

где индекс г определяет номер слоя в потоке и соответствующим ему материал; г, £ - радиальная и продольная цилиндрические координаты; уг, у^ - компоненты вектора скорости; Ф - диссипативный источник тепла; Р - давление; Т - температура; р - плотность; С - теплоемкость; X - теплопроводность; - эффективная вязкость, являющаяся функцией скорости сдвига и температуры [11]:

цэ _цое-р(Т-Т0)2) 2

где ц0 - начальная вязкость при Т0; в - температурный коэффициент вязкости; п - коэффициент аномалии вязкости; 12 - второй инвариант тензора скоростей деформации.

Система дифференциальных уравнений (1)-(5) замыкалась следующими граничными условиями:

- на неподвижных стенках компоненты скорости равны нулю;

- на границе контакта с подвижной жилой продольная компонента скорости равна скорости движения жилы (0,183м/с);

- на входе в каналы задавались давления;

- на выходе - ГУ второго рода по скорости и температуре;

- температура неподвижных стенок каналов 443 К;

- температура расплавов полимеров на выходе из экструдеров 423 К;

- температура подвижной жилы 383 К.

На границах раздела потоков задавались следующие граничные условия [12]:

п—1

по напряжениям

(Хгг С08(Г,И) + Х^ С08(С,П)) = (Хгг С08(Г,и) + Х^ С08(С,П)) |2 (Т^ С0Б(г,п) + Т?? С0Б(С,и)) ^ = (Т^ С0Б(г, и) + Т?? С0Б(С,и)) |2 (тгг С0Б(Г,П) + Х^ С°8(С,п)) |2 = (тгг С0Б(г,п) + Х^ С08(С,п)) |з (Х^ С0Б(г,п) + Т?? С0В(С,И)) |2 = (Т^ С08(г,и) + С0Б(С,й)) |з

по скоростям

по потокам тепла

Ч,; ^ = Ч? |2,

по температуре на границах раздела потоков

1 11 1 12'

' 12 ' 1з>

где индексы 1, 2, з - номера потоков стратифицированного течения (см. рис. 2); тгг хг^ - нормальные и касательные напряжения, определяемые как

£ т,

^ э

(Эу^ д^О

ЭС + дг

где Чг, Ч? - потоки тепла в направлении соответствующих осей.

Сформулированная таким образом система дифференциальных уравнений (1)-(4), дополненная реологическими уравнениями состояния (5), граничными условиями (6)-(10) и свойствами перерабатываемых материалов (таблица) решалась численно, с применением метода конечных элементов в программном комплексе ЛиБуБ.

Для оценки стабильности решения было проведено два расчета: один в режиме заданного расхода, другой - в режиме заданного перепада давления в каналах.

Для сравнения модели с заданными расходами (назовем ее «модель Q») с результатами, полученными для модели с заданными перепадами давления (назовем ее «модель Р»), были рассчитаны поля тем-

ператур, скоростей, распределения потоков, а также построены эпюры скорости по некоторым поперечным сечениям каналов (рис. 3).

Физико-реологические свойства исследуемых материалов

Материал Начальная вязкость Мо при Т = 433К, Пас Температурный коэффициент вязкости в, 1/К Показатель аномалии вязкости п Плотность р, кг/м3 Коэффициент теплопроводности X, Вт/

Изоляция (слой 2) 14946 0,0168 0,542 779 0,182 2500

Экран (слои 1, 3) 38523 0,0027 0,251 1080 0,182 2500

Для модели Р на входах в каналы задавались давления в узлах в виде прямоугольных эпюр.

Расход в модели Q задавался через узловые скорости для каждого из трех входов, вектор скорости ориентирован по нормали к линии, ограничивающей вход каждого канала. Эпюры скорости имели прямоугольный вид.

Рис. 3. Сравнение эпюр скоростей для различных сечений и моделей:--«модель

Q»;----«модель Р»; сечение А и сечение Б обозначены на рис. 2

Анализ эпюр скорости показал максимальное отклонение результата модели Р от модели Q в сечении Б (т.е. ближе к выходу из расчетной области) 8,4 % (при этом среднее отклонение скорости в этом сечении не превышает 3,9 %). Для сечения А максимальное отклонение скорости составило 2,3 %, среднее - 0,9 %.

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

Эпюра скорости расплава полимера в узком канале, образованном головкой и подвижной жилой, в поперечном сечении имеет явно выпуклый характер (скорость по центру этого канала более чем в 2,5 раза превышает скорость движения жилы), что объясняется малым поперечным сечением кольцевого канала и высоким перепадом давления. По мере слияния потока 1 и совместного потока 2 и 3 на расстоянии между сечением А и сечением Б происходит смещение границы раздела этих потоков в направлении от жилы, тем самым поперечное сечение кольцевого канала увеличивается, скорости потока 1 выравниваются и приближаются к скорости движения жилы (см. рис. 2).

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

Анализ полученных кривых показал, что для стабилизации решения достаточно 1500-2000 счетных итераций. Дальнейшее увеличение количества итераций на характер распределения потоков материала в установившемся стратифицированном течении практического влияния не оказывало. Следовательно, в расчетах при той же геометрии модели и той же сетке разбиения на конечные элементы можно ограничиться установленным количеством итераций.

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

мм 1,0 0,9

0,8 0,7 0,6 0,5 0,4

1 к

с

—X- -

й2, мм 1,6 1,4 1,2 1,0

0,8 0,6 0,4

0 200 400 600 800 1000 1200 1400 1600 1800 2000 I

Рис. 4. Изменение толщин слоев материала при увеличении количества счетных итераций (/): —д— внешний экран по изоляции (А1); —к— внутренний экран

по жиле

И1, мм

1,5 1,4 1,3 1,2 1,1 1

0,9

- \

V —

-4+ ч N._4

-А

Н2, мм 0,83 0,805 0,78 0,755 0,73 0,705 0,68

0

100000

200000

300000

400000

500000

N

Рис. 5. Изменение толщин слоев материала при увеличении количества элементов

Из рис. 5 видно, что при увеличении количества конечных элементов от двухсот тысяч и выше характер распределения потоков не менялся, на что указывала стабилизация границ раздела потоков разнородных материалов на выходе расчетной области.

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

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

Подобное явление достаточно хорошо известно и описано, например, в [14], где автор отмечает влияние такого вихревого течения на качество получаемого изделия.

Рис. 6. Завихрение потока в канале исходной кабельной головки

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

Рис. 7. Отсутствие завихрений материала в канале предложенной

геометрии

Так, после предварительного анализа существующей кабельной головки были определены недостатки ее геометрии и предложена новая конструкция канала истечения полимера, формирующая более рациональные потоки материала.

На рис. 8 приведен фрагмент каналов кабельной головки, на котором характеристически представлено отличие исходной кабельной головки от предложенного варианта.

Штриховой линией показано предложенное утонение канала, препятствующее возникновению завихрения.

С целью определения степени влияния температурной зависимости вязкости на возникновение завихрения материала были проведены исследования, подробно описанные в [15]. В целом они показали, что эффект завихрения обусловлен именно наличием сложной модели вязкости, учитывающей, помимо всего прочего, и температуру материала. В модели, где вязкость зависела исключительно от скорости сдвига, подобного завихрения не образовывалось.

Рис. 8. Различия в геометрии каналов исходной и предложенной кабельной головки

Предложенная геометрия кабельной головки была численно исследована на возможность переработки материалов, физико-реологические свойства которых отличаются от представленных в данной работе (плотности в диапазоне 400-1400 кг/м3, п = 0,1...0,7, цо = = 2000.80000). Результаты позволяют сделать вывод о возможности такой переработки без изменения предложенной конструкции. Однако для определения возможности переработки материалов, свойства которых существенно отличаются от приведенных, или описываются иными законами, или моделями (например, ПВХ), требуются дополнительные исследования.

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

На рис. 9 представлено поле давления внутри кабельного канала. Видно, что наибольшее давление создается в канале, по которому истекал полимер, образующий внутренний полупроводящий слой по жиле. Это связано с наименьшим поперечным сечением данного канала и с достаточно высокой вязкостью полимера.

Contours of Static Pressure

Рис. 9. Поле давлений в каналах кабельной головки

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

Contours of Velocity Magnitude

Рис. 10. Поле скоростей в каналах кабельной головки

Анализ температурного поля, приведенного на рис. 11, показал достаточно равномерный характер прогрева материала в канале. При этом на выходе полимер практически сохранял температуру расплава, получаемого от экструдера (рис. 12), - 423 К, а средняя температура на выходе 426 К. На неподвижных стенках температура соответствовала температуре подогрева кабельной головки, что улучшало прохождение материала через сечение канала.

Contours of Static Temperature

Рис. 11. Поле температур в каналах кабельной головки

T, K 440 430 420 410 400 390 380

Y < 2 5

\ >— \ 1 \ 1 N. I

1 1

1 1 1

1 1 1

1 1 1

1 1 1 1

0

1

2

3

6

7

А103, м

Рис. 12. Распределение температуры в сечении Б (И - высота канала): 1 - граница полупроводящего экрана по изоляции; 2 - граница полупроводящего экрана по жиле; 3 - токопроводящая жила

Рис. 11-12 показывают невозможность оперативного контроля вероятных диссипативных перегревов полимера в кабельной головке при помощи измерительных инструментов на выходе из головки, так как температуры выравнивались и уловить перегрев не представлялось возможным. Таким образом, предлагаемый способ контроля внутренних перегревов на основе численного моделирования представляется одной из возможных альтернатив контроля качества изоляции на предмет наличия следов термической деструкции материала.

Рис. 13. Распределение потоков материала внутри каналов кабельной головки

Рис. 13 демонстрирует траектории движения расплава полимеров по коническо-цилиндрическим каналам кабельной головки и отражает взаимодействие потоков при их слиянии. На рис. 14-15 отдельно показаны границы потоков расплавов полимеров. Можно отметить стабильность границ на всем участке взаимодействия потоков, что обеспечивает хорошее качество получаемой многослойной изоляции.

Рис. 14. Границы раздела потоков

Границы потоков, изображенные на рис. 15, позволили определить толщины накладываемых на жилу полимерных полупроводящих и изолирующих слоев.

Рис. 15. Границы раздела потоков на выходе из кабельной головки

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

Графики, приведенные на рис. 16, показывают изменение толщины накладываемых слоев материала при изменении массового расхода в первом (01), втором (02) и третьем ^3) каналах. При увеличении расхода Q1 в 10 раз толщина экрана по жиле возрастала в 6,8 раз, толщина изоляции уменьшалась на 20 %, а толщина экрана по изоляции уменьшалась на 30 %. При увеличении расхода Q2 в 8 раз толщина экрана по жиле уменьшается в 2 раза, толщина изоляции возрастает в 3,7 раза, а толщина экрана по изоляции падает в 5,7 раза. При увеличении расхода Q3 в 10 раз толщина экрана по жиле падает на 20 %, толщина изоляционного слоя падает на 40 %, а толщина экрана по изоляции возрастает в 5 раз. Анализ графиков позволил сделать вывод о минимальном влиянии расхода первого канала на два других.

На рис. 17 показана зависимость толщин слоев от линейной скорости. При увеличении Ул в технологическом диапазоне в 8 раз толщина экрана по жиле уменьшалась в 3,5 раза, толщина изоляции и толщина экрана по изоляции возрастали соответственно на 12 и 28 %.

а

б

в

Рис. 16. Изменение толщин слоев материала при изменении расхода: а - в канале 1; б - в канале 2; в - в канале 3; —д— экран по изоляции; ---♦-- изоляция; —х— экран

по жиле

Толщины Толщина

Рис. 17. Изменение толщин слоев материала при изменении линейной скорости ¥л: —&— экран по изоляции; ------- изоляция; —*— экран по жиле

Толщины слоев контролировались на выходе кабельной головки. Таким образом, суммарная толщина всех слоев фиксирована и ограничена высотой канала на выходе. При увеличении линейной скорости эпюра скорости материала на выходе теряла выпуклый характер: материал увлекался подвижной стенкой (жилой), и характер эпюры смещался ближе к течению Куэтта. Скорости второго и третьего потоков уменьшались, следовательно, для сохранения требуемого расхода было необходимо увеличение ширины каждого из потоков.

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

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

Результаты данного исследования с успехом могут быть применены при изучении процессов неизотермических многослойных течений в кабельных головках; при проектировании систем автоматизированного управления экструзионными линиями [16, 17].

Работа выполнена при финансовой поддержке РФФИ (грант №13-08-96034).

Библиографический список

1. The Finite Element Simulation of Polymer Coextrusion Based on the Slip Boundary / M. Zhang, C. Huang, S. Sun, Y. Jia // Polymer-Plastics Technology and Engineering. - 2009. - No. 48. - P. 754-759.

2. Visualisation and Analysis of LDPE Melt Flows in Coextrusion Geometry / M.T. Martyn, T. Gough, R. Spares, P.D. Coates, M. Zatloukal // SPE ANTEC. - 2002. - No. 60. - P. 937-941.

3. Experimental Observations of LDPE Melt Flow in Coextrusion Geometries / M.T. Martyn, T. Gough, R. Spares, P.D. Coates, M. Zatloukal // SPE ANTEC. - 2004. - No. 62. - P. 205-209.

4. Theoretical and experimental instabilities in coextrusion analysis of interfacial flows / M. Zatloukal, W. Kopytko, A. Lengalova, J. Vlcek // J. Appl. Polym. Sci. - 2005. - No. 98 (1). - P. 153-162.

5. Martyn M.T., Coates P.D., Zatlouka M. Visualisation and Analysis of Polyethylene Coextrusion Melt Flow // AIP Conference Proceedings, 7/24/2009. - 2009. - Vol. 1152. - Iss. 1. - P. 96-109.

6. Imaging and analysis of wave type interfacial instability in the coextrusion of low-density polyethylene melts / M.T. Martyn, R. Spares, P.D. Coates, M. Zatloukal // Journal of Non-Newtonian Fluid Mechanics. -2009. - No. 156. - P. 150-164.

7. Mitsoulis E., Heng F.L. Numerical simulation of coextrusion from a circular die // J. Appl. Polym. Sci. - 1987 - No. 34 (4). - P. 1713-1725.

8. Раувендаль К. Экструзия полимеров - СПб.: Профессия, 2008. - 786 с.

9. Khan A.A., Han C D. Trans.Soc.Rheol. - 1976. - Vol. 20. - No. 4. -P. 595-621.

10. Математическое моделирование соэкструзии длинномерных кольцевых изделий из резиновых смесей / П.П. Юрыгин [и др.] // Научно-технический вестник Поволжья. - 2013. - № 3. - С. 35-37.

11. Yankov V.I., Trufanova N.M., Shcherbinin A.G. Nonisothermal flow of polymer solutions and melts in channels of constant cross section // Theoretical Foundations of Chemical Engineering. - 2004. - Vol. 38. -No. 2. - P. 179-188.

12. Zhao R., Macosko W. Slip at molten polymer-polymer interfaces // MRS Symp. - 2000. - P. 629.

13. Казаков А.В., Савченко В.Г., Труфанова Н.М. Расчет влияния геометрии каналов технологического инструмента кабельной головки на возникновение вихревых потоков при наложении изоляции // Кабели и провода. - 2010. - № 2 (321). - С. 11-13.

14. Tanner R.I. Some Experiences Using Finite Element Methods in Polymer Processing and Rheology // Proceedings of the Seventh International Congress on Rheology. - Gothenburg, Sweden, 1975. - P. 140.

15. Казаков А.В., Савченко В.Г., Труфанова Н.М. Моделирование процессов тепломассопереноса полимера в головке экструдера с учетом и без учета зависимости вязкости от температуры // Интеллектуальные системы в производстве. - 2010. - № 1. - С. 130-134.

16. Казаков А.В., Труфанова Н.М. Численные исследования режимов стратифицированного течения и методика управления процессом экструзионного наложения многослойной изоляции // Известия Том. политехн. ун-та. - 2012. - Т. 320, № 4. - С. 167-171.

17. Kazakov A.V., Trufanova N.M. A system for adaptive monitoring of the process of polymer insulation production // Russian Electrical Engineering. - 2012. - Vol. 83. - Iss. 11. - P. 640-643.

References

1. Zhang M., Huang C., Sun S., Jia Y. The Finite Element Simulation of Polymer Coextrusion Based on the Slip Boundary. Polymer-Plastics Technology and Engineering, 2009, no. 48, pp. 754-759.

2. Martyn M.T., Gough T., Spares R., Coates P.D., Zatloukal M. Visualisation and Analysis of LDPE Melt Flows in Coextrusion Geometry. SPE ANTEC, 2002, no. 60, pp. 937-941.

3. Martyn M.T., Gough T., Spares R., Coates P.D., Zatloukal M. Experimental Observations of LDPE Melt Flow in Coextrusion Geometries. SPE ANTEC, 2004, vol. 62, pp. 205-209.

4. Zatloukal M., Kopytko W., Lengalova A., Vlcek J. Theoretical and experimental instabilities in coextrusion analysis of interfacial flows. J. Appl. Polym. Sci, 2005, no. 98 (1), pp. 153-162.

М.В. Eanypuna, A.B. Ka3aKoe, H.M. Tpy$anoea

5. Martyn M.T., Coates P.D., Zatlouka M. Visualisation and Analysis of Polyethylene Coextrusion Melt Flow. AIP Conference Proceedings, 7/24/2009, 2009, vol. 1152, iss. 1, pp. 96-109.

6. Martyn M.T., Spares R., Coates P.D., Zatloukal M. Imaging and analysis of wave type interfacial instability in the coextrusion of low-density polyethylene melts. Journal of Non-Newtonian Fluid Mechanics, 2009, no. 156, pp. 150-164.

7. Mitsoulis E., Heng F.L. Numerical simulation of coextrusion from a circular die. J. Appl. Polym. Sci, 1987, no. 34 (4), pp. 1713-1725.

8. Rauvendal' K. Ekstruziia polimerov [Polymer extrusion]. St. Petersburg: Professiia, 2008. 786 p.

9. Khan A.A., Han C.D. Trans. Soc. Rheol, 1976, vol. 20, no. 4, pp. 595-621.

10 Yurygin P.P. [et al.] Matematicheskoe modelirovanie soekstruzii dlinnomernykh kol'tsevykh izdelii iz rezinovykh smesei [Mathematical modeling of long annular coextrusion products from rubber compounds]. Nauchno-tekhnicheskii vestnikPovolzh'ia, 2013, no. 3, pp. 35-37.

11. Yankov V.I., Trufanova N.M., Shcherbinin A.G. Nonisothermal flow of polymer solutions and melts in channels of constant cross section. Theoretical Foundations of Chemical Engineering, 2004, vol. 38, no. 2, pp. 179-188.

12. Zhao R., Macosko W. Slip at molten polymer-polymer interfaces. MRS Symp, 2000, p. 629.

13. Kazakov A.V., Savchenko V.G., Trufanova N.M. Raschet vliianiia geometrii kanalov tekhnologicheskogo instrumenta kabel'noi golovki na vozniknovenie vikhrevykh potokov pri nalozhenii izoliatsii [The calculation of the geometry of the channel process tool cable head on the occurrence of eddy currents when applied insulation]. Kabeli iprovoda, 2010, no. 2 (321), pp.11-13.

14. Tanner R.I. Some Experiences Using Finite Element Methods in Polymer Processing and Rheology. Proceedings of the Seventh International Congress on Rheology. Gothenburg, Sweden, 1975. P. 140.

15. Kazakov A.V., Savchenko V.G., Trufanova N.M. Modelirovanie protcessov teplomassoperenosa polimera v golovke ekstrudera s uchetom i bez ucheta zavisimosti viazkosti ot temperatury [Simulation of heat and mass transfer processes in polymer extrusion head with and without tem-

perature dependence of viscosity]. Intellektual'nye sistemy v proizvodstve, 2010, no. 1, pp. 130-134.

16. Kazakov A.V., Trufanova N.M. Chislennye issledovaniya rezhi-mov stratifitsirovannogo techeniya i metodika upravleniya protsessom ek-struzionnogo nalozheniya mnogosloinoi izolyatsii [Numerical studies of stratified flow regimes and governance process extrusion multilayer insulation overlay]. Izvestiia Tomskogo politekhnicheskogo universiteta, 2012, vol. 320, no. 4, pp. 167-171.

17. Kazakov A.V., Trufanova N.M. A system for adaptive monitoring of the process of polymer insulation production. Russian Electrical Engineering, 2012, vol. 83, iss. 11, pp. 640-643.

Об авторах

Бачурина Мария Владимировна (Пермь, Россия) - аспирант кафедры конструирования и технологии в электротехнике Пермского национального исследовательского политехнического университета (614990, г. Пермь, Комсомольский пр., 29, e-mail: ktei@pstu.ru).

Казаков Алексей Владимирович (Пермь, Россия) - кандидат технических наук, доцент кафедры конструирования и технологии в электротехнике Пермского национального исследовательского политехнического университета (614990, г. Пермь, Комсомольский пр., 29, e-mail: ktei@pstu.ru).

Труфанова Наталия Михайловна (Пермь, Россия) - доктор технических наук, профессор, заведующая кафедрой конструирования и технологии в электротехнике Пермского национального исследовательского политехнического университета (614990, г. Пермь, Комсомольский пр., 29, e-mail: ktei@pstu.ru).

About the authors

Maria V. Bachurina (Perm, Russian Federation) - graduate student of Design and Technology in Electrical Engineering, Perm National Research Polytechnic University (29, Komsomolsky av., 614990, Perm, Russian Federation, e-mail: ktei@pstu.ru).

Alexey V. Kazakov (Perm, Russian Federation) - Ph.D. in Technical Sciences, Assistant Professor of Design and Technology in Electrical Engineering, Perm National Research Polytechnic University (29, Komsomolsky av., 614990, Perm, Russian Federation, e-mail: ktei@pstu.ru).

Natalia M. Trufanova (Perm, Russian Federation) - Doctor of Technical Sciences, Professor, Head of Design and Technology in Electrical Engineering Department, Perm National Research Polytechnic University (29, Komsomolsky av., 614990, Perm, Russian Federation, e-mail: ktei@pstu.ru).

Получено 13.03.2014

Просьба ссылаться на эту статью в русскоязычных источниках следующим образом:

Бачурина М.В., Казаков А.В., Труфанова Н.М. Математическое моделирование процесса стратифицированного течения расплавов полимеров в осесимметрич-ной постановке // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2014. - № 2. - С. 102-124.

Please cite this article in English as:

Bachurina M.V., Kazakov A.V., Trufanova N.M. Mathematical modelling of stratified flow of polymer melts in an axisymmetric formulation. PNRPU Mechanics Bulletin. 2014. No. 2. P. 102-124.