3115
.pdfстает вдоль положительного направления оси *, причем форма профиля пленки б (.х) — неизвестная функция. Тепло передается от нагретой пластины, которая имеет постоянную температуру Т0, к поверхно сти раздела фаз с температурой Т = Тт (предполагается, что поли мер является частично-кристаллическим веществом). Как было показано в разд. 9.2, аморфные полимеры в области Тё не имеют четко выраженного перехода от хрупкого состояния к вязкому. Поэтому выбор значения Тт не является очевидным. Можно произ вольно выбрать температуру Тт (выше Тб), при которой полимер наверняка становится вязкой жидкостьюД В противоположность этому, как недавно предположили Сандстром и Ло [27 ], температура стеклования совместно с уравнением ВЛФ [см. уравнение (6.4-10)] может быть использована для определения ТП1 расчетным путем.
ЩтгКу будем искать аналитическое решение, определяющее ско рость плавления и распределение температур в образующемся слое расплава. Ясно, что эти переменные являются функциями физиче ских свойств твердого полимера, температуры и скорости нагретой
пластины и ширины бруса Первый, кто рассматривал эту задачу в приложении к изучению
механизма плавления полимерной пробки в червячном экструдере, был Тадмор с сотр. [28—30]. Позднее Вермеулен [31] и Сандстром [32, 37 ] исследовали этот вопрос экспериментально и теоретически,
Маунт |
[33] экспериментально |
определил скорость плавления, |
а Пирсон [34] дал теоретический |
анализ этой проблемы. Восполь |
|
зуемся |
анализом Пирсона. |
|
Чтобы детально разобраться в механизме плавления при описан ных условиях, рассмотрим свойства твердого полимерного стержня. Для абсолютно твердого, несжимаемого тела, надвигаемого на нагре тую пластину без вращения, скорость плавления на поверхности раздела фаз не должна зависеть от координаты х , потому что скорость твердой фазы в любом сечении х одинакова. Следовательно, б (х), Р (.х) и поля скоростей и температур в пленке расплава должны принимать значения, которые будут удовлетворять как этому тре бованию, так и уравнениям движения и энергии при соответствую щих граничных условиях. Однако в тонких пленках сильновязких полимеров при больших скоростях сдвига более приемлемым является предположение о постоянстве давления в пленке. Это в свою оче редь дает основание предполагать, что при установившемся режиме скорость плавления в общем случае зависит от х, хотя эта зависи мость может быть очень слабо выражена.
Переменная скорость плавления означает, что твердая фаза под вергается или деформации, или вращению, или тому и другому вместе. Твердые полимеры, в частности в виде пробки спрессованных гранул или порошков (как это обычно наблюдается в процессах пере работки), можно считать деформируемыми. Расплав, образующийся в очаге плавления, проникает внутрь пустот между твердыми ча стицами пробки, позволяя им скользить и перестраиваться в области, прилегающей к поверхности раздела фаз. Физическая сущность деформации твердой пробки состоит в следующем. Медленно дефор-
282
мируемый твердый стержень придвигается к движущейся нагретой поверхности. Граница стержня имеет незначительную скорость в отрицательном направлении оси у , которая может несколько изме няться при изменении х. Однако можно считать, что при опреде ленном значении х стержень обладает достаточной прочностью, чтобы выдержать сдвиговые напряжения, возникающие в пленке расплава и тем самым препятствовать возникновению компоненты скорости в направлении оси х.
Исходя из этого, Изведем допущения, позволяющие упростить исходное дифференциальное уравнение: 1) теплофизические свойства постоянны; 2) расплав — несжимаемая жидкость; 3) на стенках нет проскальзывания; 4) справедлив степенной закон течения ненью
тоновской |
жидкости с вязкостью, зависящей |
от |
температуры: |
|
|
т = т0е “ (Т Т'^ |
|
|
(9.8-1) |
5) процесс |
плавления — установившийся; |
6) гравитационные силы |
||
пренебрежимо малы; 7) течение ламинарное; |
8) |
толщина пленки |
||
во много раз меньше ее ширины (6/№ |
1). Поэтому с учетом ма |
лости числа Рейнольдса в пленке расплава справедливо применение для этой задачи смазочной аппроксимации, позволяющей пренеб речь эффектом выхода (при х = W) и не задавать точных условий
входа (при х = 0). |
|
|
сводятся в этом случае |
||||
Уравнения неразрывности и движения |
|||||||
к виду: |
|
|
|
|
|
|
|
dv^ |
^dvy_ |
(9.8-2) |
дР |
_ |
dTjp |
(9.8-3) |
|
дх ^ |
ду |
дх |
~ |
ду |
|||
|
|
Поскольку предполагается, что течение в пленке является чисто сдвиговым, то уравнение (9.8-3) дополнительно упрощается:
атУ* = о |
(9.8-4) |
ду |
|
Выразив напряжение сдвига в уравнении (9.8*4) через локаль ный градиент скорости, получим:
(9.8-5)
Интегрирование уравнения (9.8-5) по у дает:
дУх |
с а (Т - Т т)/п |
(9.8-6) |
|
ду |
1 |
||
|
Таким образом, если а — 0 (т. е. вязкость не зависит от темпе ратуры), то профиль скорости будет линейным как для ньютонов ской, так и для степенной жидкости. Если, однако, а ф 0, то ло кальный профиль скорости становится функцией температуры. Так как температура меняется резко вдоль у, то можно ожидать, что профиль скоростей в направлении оси у будет существенно нелинейным. Кроме того, вследствие конвекции температура Т 9дабо зависит от х, что приводит к появлению слабой зависимости
профиля скоростей от х. Следовательно, уравнения движения и энергии должны решаться одновременно.
Уравнение энергии в данном случае имеет вид:
где Pm., Ст и km — теплофизические свойства также, что теплопередача в направлении оси оси у, и что наибольший вклад в вязкостную напряжений, которые содержат компоненты
(9.8-7)
расплава полимера. Предполагается х много меньше, чем в направлении диссипацию вносят те члены тензора
тух.
Определим граничные условия для пленки расплава:
Т (0) = 7\>; vx (0) = V„; |
vu (0) = 0 |
(9.8-8а) |
Т (6) = Tmvx (6) = |
0 |
(9.8-86) |
Скорость vy (б) в любой точке х определяется скоростью плав ления на поверхности раздела фаз и может быть получена из урав нения теплового баланса:
km {-'%~)у=б ~ Pm [~ Vy (6)1 Х + ks ^ f ) y =б |
(98'9) |
где А,— удельная теплота плавления; ks и Ts — коэффициент теплопроводности и температура твердой фазы.
Левая часть уравнения представляет собой поток тепла, посту пающий за счет теплопроводности от горячей поверхности пленки
кочагу ппавления (к единичной поверхности раздела). Первый член
вправой части — произведение интенсивности плавления на еди ничной поверхности раздела на теплоту плавления; второй член поток тепла, отводимого от поверхности раздела. Так как происходит плавление, то дТ/ду < 0. Поэтому левая часть уравнения положи тельна и выражает тепло, которое, как это видно из правой части уравнения, используется для двух целей: плавления расплава при
Т= Тт и нагрева полимера до температуры плавления.
Второй чаен правой части уравнения может быть получен из ре шения, описывающего температурный профиль в слое твердого полимера. Рассмотрим элементарный участок, выделенный на оси х, в пленке расплава и слое твердой фазы (см. рис. 9.13]. Предполо жим, что твердая фаза занимает область у > 6 (где о — толщина пленки в данном месте) и движется с постоянной скоростью vxy
кповерхности раздела фаз. Задача, таким образом, сводится к ре шению стационарной одномерной задачи теплопроводности конвек цией. В твердой фазе устанавливается экспоненциальное распреде ление температуры подобно тому, как это имело место в задаче, описанной в разд. 9.5. Уравнение энергии в данном случае сводится
квиду:
п Г>.. 01 S |
и |
° 1S |
(9.8-10) |
PsCsVsy - Щ - |
■ - ks |
дуЪ |
где Ps, Cs и ks — теплофизические свойства твердой фазы полимера,
254
Уравнение (9.8-10) легко интегрируется при граничных усло виях Ts (б) = Тт и Та (оо) = Ts0 и дает следующий профиль тем пературы:
|
Т |
Tso + |
(Т1П- |
Tso) exp [ Vsy {y~ |
h) ] |
(9.8-11) |
||
Так как |
vsy < 0, |
то |
уравнение |
(9.8-11) удовлетворяет обоим |
||||
граничным |
условиям. |
Поток |
тепла |
от поверхности раздела фаз |
||||
с учетом того, что |
v8yp8 = vy (б) pm, определяется |
как |
||||||
|
^s ( |
|
ду )^=б = |
^ m |
Vy ^ |
PmCs |
(9.8-12) |
|
Таким образом, |
уравнение (9.8-9) может быть |
записано так: |
||||||
где |
|
|
|
(~ду~)//=б = |
PmVy ^ |
|
(9.8-13) |
|
|
|
** = * + |
Cs (Tm — Ts0) |
|
(9.8-14) |
|||
|
|
|
|
А,* — суммарное количество тепловой энергии, необходимое для нагрева твердой фазы от начальной температуры Ts0 до темпера туры Тт и плавления при этой температуре. Сандстром и Юнг [32] решили эту систему уравнений численным методом, заменив урав нения в частных производных уравнениями в обыкновенных произ водных на основе методов теории подобия. Пирсон [34] использо вал аналогичный подход и получил ряд аналитических решений для более простых счучаев. Он использовал безразмерные перемен ные, которые полезны, как это будет далее показано, при физиче ской интерпретации результатов:
|
о |
Т '— Тт |
|
(9.8-15) |
|
|
|
Т0 - Т т |
|
||
|
|
|
|
||
|
l = x/W-, г| = у/6 |
|
(9.8-16) |
||
“* *= |
и» = |
l/0 (60V f |
(9.8-17) |
||
(значение б0 будет объяснено ниже). |
|
|
|
||
Перепишем граничные условия: |
|
|
|
||
0(0) — 1; и*(0)- |
1; «»(0)-.0 |
(9.8-18з) |
|||
|
в (1) = |
0; и*(1) = |
0 |
(9.8-186) |
|
Условие плавления |
на поверхности |
раздела фаз |
[уравнение |
||
(9.8-13)] сводится к виду: |
|
|
|
|
|
km{T t-T m)W |
|
|
|
(9.8-19) |
|
ь'РпУЛ |
|
|
|
||
|
|
|
|
||
Это соотношение |
позволяет |
сделать обоснованный вы |
бор 60. Поскольку заданные условия определяют конкретный
физический процесс |
плавления, |
то |
безразмерная |
группа |
||
km (Т0 — Тт) Щ^*р„(У0бо) должна |
быть порядка единицы. |
Следо |
||||
вательно, можно выбрать |
р0 следующим |
образом: |
|
|||
с |
/ |
km (T0- T m)W |
у /2 |
(9.8-20) |
||
б° ~ \ |
Щ М , |
I |
||||
|
нении энергии можно пренебречь. Безразмерный параметр b — мера температурной зависимости вязкости в рассматриваемом диапазоне температуры (величина, пропорциональная энергии активации вяз кого течения).
Полное решение уравнений, описанных выше, как указывалось ранее, сложно. В дополнение к численному решению [32] Пирсон [34] предложил аналитический подход. Исследование природы плавления при удалении расплава за счет сдвига может быть осуще ствлено путем рассмотрения нескольких специальных случаев, для которых можно получить аналитические решения. Эти упрощенные случаи дают очень полезную информацию для моделирования, для процессов переработки.
Ньютоновская жидкость, вязкость которой не зависит от температуры (конвекция пренебрежимо мала)
Для ньютоновской жидкости при условиях, близких к изотер мическим (п = 1, | b | 1), и при пренебрежимо малой конвекции (М > 1) уравнение движения имеет вид:
д2их |
= 0 |
(9.8-27) |
дг\2 |
|
|
Для граничных условий, определяемых уравнением |
(9.8-18), |
|
его решение записывается следующим образом: |
|
|
Их-= |
1 — л |
(9.8-28) |
Уравнение энергии, которое в этом случае может быть решено |
||
независимо, имеет вид: |
|
|
д2е |
дих \2 |
|
ди2 |
дц ) = 0 |
(9.8-29) |
Подставляя уравнение (9.8-28) в (9.8-29) и интегрируя получен
ное уравнение, имеем для |
профиля температуры: |
|
||||
|
0 = 1 |
— Л |
+ - ^ г 1( 1 - л ) |
(9.8-30) |
||
Средняя |
температура 0 определяется из уравнения (9.8-30) как |
|||||
|
|
1 |
|
|
|
|
|
|
[ нх-0 dr\ |
|
|
||
|
* = h ------- = 4 |
+ -БГ |
(9'8-3,) |
|||
|
|
| |
UXrfl] |
|
|
|
|
|
о |
|
|
|
|
Теперь можно решить уравнение (9.8-19) относительно |
иу (1), |
|||||
подставляя |
(дО/дт)),,^ = —(1 |
+ |
Вг/2) из (9.8-30): |
|
||
|
«ио |
|
- |
4 4 |
1 + - ^ - ) |
(9-8'32) |