CC BY
41
УДК 532.135
В. Н. Харитонов, В. М. Шаповалов
ВАЛКОВОЕ ТЕЧЕНИЕ ЖИДКОСТИ ОСТВАЛЬДА-ДЕ ВИЛЯ ПРИ ЗНАЧИТЕЛЬНОЙ ФРИКЦИИ
Предложен способ решения задач течения неньютоновских жидкостей. Способ может найти применение при анализе течений, в которых доминирует простое сдвиговое течение с известными параметрами.
Ключевые слова: давление, касательное напряжение, компоненты скорости, второй инвариант тензора скоростей деформации.
Для интенсификации процесса перемешивания в валковых машинах используют режим работы со значительной фрикцией. Реологические свойства значительной части полимеров и резин можно описать законом Освальда-де Виля. Однако использование степенного закона приводит к значительным техническим трудностям анализа модели [1]. В данной работе предпринята попытка упрощения задачи путём линеаризации второго инварианта тензора скоростей деформации.
Схема течения представлена на рисунке 1. Валки имеют одинаковые радиусы (И), но различную окружную скорость иг и и2 (и1>и2). Начало координат помещено в середине сечения минимального зазора. Начало зоны течения характеризует координата Хо, окончание - хг. Текущая полувысота зазора Ь(х). Минимальный зазор 2Но. Течение плоское, изотермическое. Среда несжимаемая. В силу условия х - X >> Ь давление
однородно по высоте зазора (др/ду = 0), следовательно, p=p(x).
Течение описывается системой уравнений [2]:
^ А' ^ ^
dx ду
(А: А) ■ 2
n—1
2 dv
(1)
ду
Харитонов Владимир Николаевич — кандидат технических наук, доцент (Волжский политехнический институт (филиал) Волгоградского государственного технического университета, Волжский); e-mail: vtm@volpi.ru.
Шаповалов Владимир Михайлович — доктор технических наук, профессор (Волжский политехнический институт (филиал) Волгоградского государственного технического университета, Волжский); e-mail: svm-5@meil.ru.
© Харитонов В. Н., Шаповалов В. М., 2016
87
Рис. 1. Схема течения
h
Q ={vxdy, (2)
-h
y = h, vx = u, (3)
y = -h, vx = u2, (4)
x = xo, p = 0, (5)
x = xi: p = 0, dp/dx = 0, (6)
где p — давление; x, y — координаты; т^ — касательное напряжение; vx, vy — компоненты скорости; ц , n — реологические константы, Q — расход жидкости, приходящийся на один метр ширины валков, (Q = const )•
Здесь (1) — уравнение движения и состояния; (2) — уравнение расхода; (3), (4) — условия прилипания; (5), (6) — условия для давления (в концепции Гаскелла).
Входящий в уравнение состояния (1) второй инвариант тензора скоростей деформации с учётом соотношений vz = 0, дvx /ду >>dvy/Эх,
h « х " х0, /ду >> dvx/дх и уравнения (2), имеет вид:
(А: А) = (Э02. (7)
2 t дУ J
88
С учётом соотношения (7) уравнение состояния (1) примет виц:
т = ur
xy г^С
dv
öy
öy
Представим осевую скорость в виде суммы:
Щ + u2 Щ - u2
- +
2h
y + V* {X y)
(8)
(9)
Первые два слагаемых правой части характеризуют простое сдвиговое течение, обусловленное фрикцией, а третье — неоднородностью давления. Задача сводится к определению функции у* (х,у). При этом
считаем правомерной оценку
dГU1+U2 | U1-U2y
dy
2h
>>
dvX {x,y)
dy
(10)
Подставив выражение (9) в (8), с учётом оценки (10) запишем выражение для касательного напряжения в форме
Txy =^0
n-V А * Л
u,-u, dv*
1 2 _L x
2h
2h dy
(11)
Согласно (11) эффекты аномалии вязкости обусловлены простой
сдвиговой компонентой течения. Эффективная вязкость
2h
од-
нородна по высоте зазора, но изменяется по его длине, поскольку Ь(х). Из граничных условий (3), (4) с учётом (9) получим условия для
* * /Л
функции Ух : у = +Ь, Ух = 0 .
Дважды проинтегрируем уравнение движения (1) с учётом выраже-
*
ния (11) и граничных условий для функции У . При этом получим:
1 Ui - U2
2Ц с 2h
^ {y2 - h2 )■ dx ;
(12)
Найдём осевой расход жидкости по формуле (2). Выполнив интегрирование с учётом (9), (12), получим:
Q = {ui -u2)h--
2h3 U1 - U 2
3Мч, 2h
d£. dx
(13)
С учётом граничного условия (6) определим расход:
n-1
V., =
x
2
U1-U2
n-1
u, -u
Vx =
89
Р = (и1 " и2 )Ь1' Ь1 = ь(х1 )•
Уравнение для давления (13) с учётом (14) запишем так:
1 Ь
dp dx
2
h1
h2
(14)
(15)
Примем параболическую аппроксимацию поверхностей валков
h = H0 (l + f2) и перейдём к безразмерным переменным:
Лее }_ {x,xo,x1}, р_
fo'f1}_ V2RH; p=
pH™ U1 - U 2
3^2RHo 2
i
, (ui * U2) . (16)
С учётом соотношений (16) уравнение (15) примет вид:
dP
1
___
df (l+f2)1+n (l+f2)2+n
(17)
Разделив переменные в (17) и проинтегрировав, с учётом граничного условия (5) в форме $ = $0, Р=0, получим выражение для безразмерного давления
P _
i-(1+ f2)
2n +1 2(1 + n) _
f df
i (1+f2 Г
(1+f2) 2(1 + n)
f
fo
.(1 + f2 )1+n (1 + f 2 )1+n
(18)
Условие (6) в безразмерной форме имеет вид $ = $1, Р=0. Используя это условие для (18), получим уравнение, связывающее параметры
$ о, $1
1 -(1 + f)
2n +1 2(1 + n )
df
(1 + f2)
(1 + f2 )1+n 2(1 + n)
f1
fo
(1 + f2)1+n -(1 + f2)1+n
_ o •
(19)
Параметры $0, $1 характеризуют безразмерные границы зоны течения. Согласно (19) их численные значения зависят от индекса течения (п), но не зависят от фрикции.
Из полученных расчётных выражений (18), (19) видно, что использованный способ выделения доминирующей компоненты скорости (простой сдвиг) к рассматриваемой задаче позволяет значительно упростить расчётные выражения; расчётные выражения по форме близки к расчётным выражениям, характерным для ньютоновской жидкости. Точность асимптотического приближения повышается с увеличением фрикции.
90
Список литературы
1. Розе Н. В. Течение полимерных материалов между двумя вращающимися цилиндрами // ДАН СССР. 1972. Т. 206. № 4. С. 834—837.
2. Шаповалов В. М. Валковые течения неньютоновских жидкостей. М.: ФИЗ-МАТЛИТ, 2011. 168 с.
•Jc -Jc -Jc
Kharitonov Vladimir N., Shapovalov Vladimir M.
THE SWATH OF FLUID FLOW OF OSTWALD-DE VILLE AT SIGNIFICANT FRICTIONS
(Volzhsky Polytechnical Institute (branch of) State Educational Institution of Higher Professional Education 'Volgograd State Technical University', Volzhsky, Volgograd region)
A method for solving the flow of non-Newtonian fluids has been offered. The method can be used in the analysis of trends, which is dominated by a simple shear flow with known parameters.
Keywords: pressure, shear stress, velocity components, the second invariant of the strain rate tensor.
References
1. Roze N. V. Techenie polimernykh materialov mezhdu dvumya vrashchayushchimisya tsilindrami [The flow of polymer materials between two rotating cylinders ], Doklady Akademii Nauk SSSR, 1972, vol. 206, no. 4, pp. 834 -837.
2. Shapovalov V. M Valkovye techeniya nen'yutonovskikh zhidkostey (Roller flow of non-Newtonian fluids), Moscow, FIZMATLIT Publ., 2011. 168 p.
* * *
91
