Примеры расчетов методом конечных разностей
В данном разделе приводятся примеры расчетов, выполненных с помощью МКР, что позволит пояснить области его возможного применения. Описание ограничивается расчетами течений с использованием явных разностных схем, так как именно они играют наиболее важную роль при моделировании процессов переработки пластмасс.
Сужающийся щелевой канал [16]
Изучение распределения скоростей и температур в щелевом канале с переменной высотой зазора показало, что профиль скорости в значительной мере зависит от геометрии канала и слабо — от температуры расплава (см. рис. 4.7).
Стенка Центр Стенка Приведенная высота канала |
Г--------------------- 111--------------------------- ч
60- |
1 |
1 мм
=1 „5" |
т |
"Г
„1" „2" „3" „4"
Форма канала (продольный разрез прямоугольной щели)
Рис. 4.7. Распределения скоростей и температур по поперечному сечению сужающегося щелевого канала при постоянной температуре стенок (200°С)
Профиль температуры достигает максимума вблизи стенки канала, так как высокая скорость сдвига приводит к повышенному тепловыделению за счет трения. С увеличением сужения канала этот максимум становится более выраженным вследствие роста диссипации энергии. Из-за низкой температуропроводности расплава выделяемое тепло может только медленно передаваться стенке канала с постоянной температурой (Tw= 200°С) или к более холодным слоям расплава в средней части канала.
На формующем участке канала (где его стенки параллельны) температура расплава в его средней части растет быстрее, при этом выделение тепла за счет трения (диссипация) уравновешивается теплоотводом через стенку канала и с утекающим расплавом. В результате этого максимум температуры смещается к средней части канала.
На рис. 4.8 показано влияние температуры стенок канала на профиль температуры при течении расплава в относительно коротком прямоугольном канале с сужающимся участком входа. Температура расплава на входе в канал (х = 0) предполагается постоянной (Тм = 170°С).
Из приведенного рисунка видно, что равномерного повышения или понижения температуры в канале за счет изменения температуры стенок не происходит. Поэтому необходимо добиваться максимально возможной термической однородности расплава, поступающего из экструдера в экструзионную головку (т. е. обеспечивать,
и
Tw= 160-С
160 171 |
Tw= 170-С |
т - 20 кг/ч; х = 0; Tw= 170'С |
Приведенная высота канала |
Н/2 |
Рис. 4.8. Профили температуры по поперечному сечению канала при различных температурах стенок |
—•— х=0 мм —о— ю » —+— 20 - —л— 30 - —о— 40 " —7— 50 •• |
по возможности, .минимальную температурную неоднородность по поперечному сечению расплава).
Угловая экструзионная головка для экструзии с раздувом [17]
В работе [17] анализируются распределение скоростей и температур в угловой экструзионной головке для экструзии с раздувом. На рис. 4.9 показан канал угловой экструзионной головки в разрезе и представлены расчетные профили скоростей, скоростей сдвига и температур. Температура массы (расплава) на входе (х = 0) постоянна и равна 200 °С; поверхность дорна предполагается адиабатической.
Расчеты позволяют выяснить, имеется ли температурная негомогенность в потоке расплава, какова ее величина, происходит ли выравнивание температуры расплава по сечению, и какое влияние оно оказывает на температурный профиль экструдата на выходе из головки. На рис. 4.10 показаны расчетные профили температуры экструдата на выходе из головки для двух рабочих режимов с экстремальными различиями в диссипации (предполагалось, что на входе имелось два различных температурных профиля с одинаковой средней температурой). Из этого примера видно, что даже предельные температуры на входе в экструзионную головку сами по себе не оказывают существенного влияния на распределение температур на выходе. Эта тенденция
Температура расплава ТМ. ”С Скорость сдвига у, с'1 Скорость v, мм/с |
Количество слоев тур в угловой экструзионной го- Внутренний Внешний ловке для экструзии с раздувом |
т — 40 кг/ч, s = 1,4 мм □ Температурный профиль 1 о Температурный профиль 2 А Тм = const = 200°С |
— прос) |
эиль 2 - |
|||
I /Темп |
It'S |
О w К |
1SX |
N |
^профиль 1/ |
о. 204 <и с 2 200 I |
Количество слоев |
220
■ Г С О д |
■ п =10 кг Темп< Темш ТМ‘ |
■ /Ч. s = 2 ?ратурнь? ратурнь const = |
- ,9 мм |й проф |й проф 200”С |
иль 1 иль 2 |
!=&s=i |
||||
210 |
200 |
Температурный — профиль 2 - |
consl
/Темпера L—турный профиль
О 2 4 6 8 10
Количество слоев
Рис. 4.10. Профили температуры на входе в экструзионную головку (нижний рисунок) и выходе из нее (верхний рисунок)
подтверждается и экспериментальными измерениями (см. [17]). Даже если в расчетах использовалась постоянная среднемассовая входная температура, обнаруживаются лишь незначительные различия температур. Необходимо принимать во внимание, что такие результаты всегда зависят от рабочих режимов.
Плоское течение при соэкструзии [11]
Численное моделирование можно применять не только к течению расплава одного материала, но и к многослойным течениям. В этом случае итерационными методами определяются не только такие параметры течения, как давления, скорости и температуры, но и расположение границ между слоями различных расплавов.
Рассмотрим асимметричное трехслойнос течение расплавов ПЭНП/ПЭВП/ПЭНП в простом прямоугольном канале размерами сечения Я х В = 12 х 200 мм. Температура массы на входе в канал (х = 0) для всех трех слоев составляет 235 °С, а температура стенки канала экструзионной головки на всей его длине — 260 °С. Зависимость вязкости от скорости сдвига описывается моделью Карро (см. раздел 2.1.1.2), а температурный сдвиг определяется в соответствии с уравнением Вильямса-Ланделла - Ферри (ВЛ Ф, см. раздел 2.1.1.3). Характеристики используемых для анализа материалов приведены в табл. 4.1. Расчет выполнялся в предположении постоянства теплопроводности и удельной теплоемкости.
Таблица 4.1. Характеристики материалов
|
На рис. 4.11 и 4.12 представлены расчетные профили скорости и температуры на входе в канал (х = 1 мм) и на различных расстояниях от него (х = 50 мм и х = 500 мм). Хотя повышенная температура стенки канала (рис. 4.11), зависящая от пути течения, не повлияла на температурные профили в середине канала, неоднородности в температурах слоев могут привести к возникновению нежелательного распределения температур в многослойном изделии, полученном методом соэкструзии (рис. 4.12).
Переходные слои Переходные слои
250 |
240 |
5 а» |
о
о
Рис. 4.11. Развитие профилей скорости и температуры при течении расплава с однородной температурой на входе и при постоянной температуре стенки
0 у/Н 1 0 у/Н 1 Относительная высота Относительная высота канала канала Н = 12 мм Материал т, кг/ч Тм. *с Tw. - С
|
Переходные слои на границе раздела различных расплавов |
Переходные слои на границе раздела различных расплавов |
Относительная высота канала |
и кг |
5 а> |
Скорость v, мм/с Скорость v, мм/с |
Рис. 4.12. |
Материал т, кг/ч |
ПЭНП ПЭВП ПЭНП |
25 100 55 |
Переходные слои на границе раздела различных расплавов |
о . CL гг- |
Относительная высота канала |
260 |
235 |
- 1 м 7 |
||
м — 50 |
ММ — |
|
0° |
К |
255 |
250 |
240 |
0 у/Н 1 Относительная высота канала Н = 12 мм |
Переходные слои |
О у/Н 1 Относительная высота канала |
Н= 12 мм |
235 235 235 |
260 |
260 |
'tv- |
Влияние различия температур расплава во внутреннем и внешних слоях на окончательное распределение температуры и скорости в экст - рудате |
ТМ.‘С |
Iw |
Материал т, кт/ч
260 |
Рис. 4.13. |
260 |
ПЭВП |
25 |
235 |
ПЭНП |
100 |
235 |
пэвп |
55 |
235 |
Влияние изменения материала на распределение скорости и температуры |
Влияние замены материала на профили скорости и температуры показано на рис. 4.13. Здесь течение, состоящее из слоев ПЭНП/ПЭВП/ПЭНГ1, было заменено течением, состоящим из слоев ПЭВП/ПЭНП/ПЭВГ1, при этом все остальные параметры (расход и температура) поддерживались на том же уровне, что и при рабочем режиме, рассматриваемом на рис. 4.11. По сравнению с описываемым случаем профили скорости и температуры резко изменились, что привело к изменению толщины слоев.
Компьютерная программа, предназначенная для расчета профилей температуры, скорости (и, следовательно, времени пребывания), профиля скорости сдвига, потерь давления, а также толщины слоев в процессе соэкструзии, позволяет моделировать процесс, изменяя параметры режима переработки и геометрические характеристики головки, определяющие толщину слоев. Таким образом, можно оценить влияние отдельных параметров на процесс в целом.
Следовательно конструкторы экструзионных головок и инженеры-технологи получают возможность оценить влияние на конечный результат вносимых в конструкцию экструзионной головки или в технологический режим изменений, во многих случаях обходясь без дорогостоящих натурных экспериментов. Кроме того, такое моделирование позволяет накапливать информацию, необходимую для конструирования экструзионных головок. Это справедливо и для тех случаев, когда принятые при расчетах допущения не позволяют получить достаточно точные конечные результаты.