CC BY
72
УДК 664 ББК 36.81 М-52
Меретуков Заур Айдамирович, кандидат технических наук, докторант кафедры технологии, машин и оборудования пищевых производств Майкопского государственного технологического университета, e-mail: zamer@radnet. ru;
Кошевой Евгений Пантелеевич, доктор технических наук, профессор, заведующий кафедрой машины и аппараты пищевых производств Кубанского государственного технологического университета, e-mail: Koshevoi@,kubstu. ru.
ИССЛЕДОВАНИЕ ТЕЧЕНИЯ НЕНЬЮТОНОВСКОЙ ЖИДКОСТИ МЕЖДУ ДВИЖУЩЕЙСЯ И НЕПОДВИЖНОЙ ПЛАСТИНОЙ
(рецензирована)
Для развития теории процесса экструзионной агломерацииЦелью данной работы является решение задачи ламинарного течения фосфатидного концентрата в щелевом канале между двумя плоскими параллельными стенками, одна из которых неподвижна, а другая движется с заданной постоянной скоростью.
Ключевые слова: фосфатидный концентрат, ламинарное течение, щелевой канал, течение Куэтта, неньютоновская жидкость.
Meretukov Zaur Aydamirovich, Candidate Of Technical Sciences, assistant professor of Technology, Machines and Food Industry Equipment Department, Maikop State Technological University, e-mail: [email protected];
Koshevoi Eugene Panteleevich, Doctor Of Technical Sciences, Professor, Honoured Scientist of the Russian Federation, head of the Department of machines and equipment for food industry, Kuban State Technological University, e-mail: Koshevoi@,kubstu. ru.
STUDY OF NON-NEWTONIAN FLUID RUNNING BETWEEN THE MOVING AND
FIXED PLATE
The aim of this work is to solve the problem of laminar flow of phosphatide concentrate into a slot channel between two flat parallel walls, one of which is fixed, while the other moves with a given constant speed.
Key words: phosphatidic concentrate, laminar flow, slotted channel, Couette flow, non-Newtonian fluid.
Прежде чем описать течения фосфатидного концентрата, обладающего свойствами неньютоновской жидкости, в развернутом канале вала одношнекового экструдера проанализируем ряд вспомогательных задач.
Исследуем частный случай ламинарного течения несжимаемой ньютоновской жидкости в щели между двумя плоскими параллельными стенками, одна из которых неподвижна, а вторая движется с заданной постоянной скоростью Wпрeд - так называемое течение Куэтта [1]. Для одномерного потока уравнение неразрывности имеет вид:
0 Wj,lxt=u „гЫ = |> D
Уравнение Навье - Стокса в этом случае примет вид:
д д2
— р\х\= и —Wjx) (2)
«с сх
Учитывая постоянство скоростного напора, для установившегося режима, имеем возможность приравнять левую часть уравнения (2) постоянной величине ДР - градиенту давления по оси z, Па/м.
Используя прямое преобразование Лапласа, преобразуем уравнение (3) к виду:
В этом случае получаем обыкновенное дифференциальное уравнение второго порядка:
(3)
(4)
Разрешая полученное уравнение (4) относительно изображения искомой функции Цв) находим полученную зависимость, считая в простой переменной:
(5)
(6)
Используя обратное преобразование Лапласа, преобразуем уравнение (5) к виду:
1 і 1 4
м^і.х.і = — — ■ АР-х2 +х - р. ■ С2 + р. ■
м и .
Постоянные интегрирования Сі и С2 находим из условия прилипания вязкой жидкости к поверхности твердого тела на подвижной (х = ії) и неподвижной (х = 0) пластинах:
(7)
Решая систему уравнений (7) установили, что С1 = 0, а С 2 определяется формулой:
.2
_ 1 А
2 1 к-р.
(8)
Уравнение скорости для ламинарного течения несжимаемой ньютоновской жидкости в щели между двумя плоскими параллельными стенками, одна из которых неподвижна, а вторая движется с заданной постоянной скоростью \Упред имеет вид:
х др.х-А-ДР-ха + /л-2-И^
нд.х! =------------------------------— (9)
2 ц-Н
Зная распределение скорости по сечению щели можно рассчитать удельный расход
жидкости:
*-{ т АРхк-АРх2 + и-21¥^Л
(10)
Окончательно удельный расход жидкости в плоской щели шириной Ь можно представить формулой:
12 О1)
Представленные уравнения (10) и (11) превращаются в формулы Куэтта для безградиентного течения жидкости (АР = 0). Таким образом, установлено, что в зависимости от величины гидродинамического напора (АР) и скорости движения пластины возможны различные режимы ламинарного течения жидкости в щелевом канале (ри. 1). На рисунке 1 видно, что при изменении перепада давления (от отрицательного до положительного значения) возможно существование обратных течений вблизи неподвижной стенки.
Рис. 1. Моделирование изменения профиля относительных скоростей по сечению щели для различных градиентов давления при постоянной скорости Wпред
В реальном щелевом канале профиль скоростей при допущении потенциальности течения представляет собой двумерный профиль, который при условии симметрии поля скоростей можно представить как суперпозицию одномерных полей, на основе полученного уравнения (10). В этом случае решение задачи можно представить как произведение частных решений по оси х и по оси у [2].
Учитывая, что по оси у значения скоростей на противоположных стенках равны 0, а в центре такого канала скорость максимальна имеем возможность получить распределение скоростей на основе суперпозиции задачи Куэтта и напорного течения в плоском канале.
Рассмотрим решение задачи для случая напорного движения жидкости между неподвижными пластинами.
Исследуем частный случай ламинарного течения несжимаемой ньютоновской жидкости в щели между двумя неподвижными плоскими параллельными стенками. Для одномерного потока уравнение неразрывности имеет вид аналогичный уравнению (1), а уравнение Навье - Стокса в этом случае аналогично уравнению (2).
Постоянные интегрирования С1 и Сг находим из условия прилипания вязкой жидкости к поверхности твердого тела на неподвижных стенках (у = +Ь) и условии симметрии (у = 0) в центре потока:
1 М
У^у'Ь * — V т-ДР™ Ь2+Ь-
1-Ь ' = 1 Гг ■ АР ■ 1- + -О 1
,2
^ ■С2 + {А ■ Сг
= 0
Ь'№'С2+рС1
= о
(12)
■ | 0 1 = 0
Решая систему уравнений (12) установили, что Сг= 0, а С1 определяется формулой:
(13)
2 ^
Уравнение скорости для ламинарного течения несжимаемой жидкости в щели между двумя плоскими параллельными стенками имеет вид:
VI/
.(>■)=
АР.
!/
ь2)
(
2/1 ' (И)
В центре скорость потока максимальна и определяется градиентом давления ДРнап:
(15)
Учитывая, что скорость в центре потока задачи Куэтта и напорного течения по условию симметрии совпадают, имеем возможность, определить скорость по уравнению (15), приравняв её к скорости по оси х , взятой из задачи Куэтта (9):
(16)
Выразим напорную составляющую давления ДРнап как функцию решения задачи
Куэтта:
(17)
Получаем распределение скоростей напорного движения жидкости между неподвижными стенками:
(
(18)
Используя распределение скоростей, полученное по уравнению (18) проводили расчет модели скоростного напора в канале прямоугольного сечения с подвижной стенкой. Профили скоростей для противодавления (ДР = -100), безнапорного движения жидкости (ДР = 0) и напорного движения (ДР = 100) представлены (рис. 2) для сравнения в виде совмещенной поверхностной диаграммы, а для практически важного случая с противодавлением (ДР = -100) также построена диаграмма (рис. 3).
Практическое значение полученные зависимости имеют для моделирования процесса экструзионной обработки пищевой массы с реологическими свойствами и неньютоновского типа как первое приближение, однако более сложный характер реологических свойств реальной жидкости (фосфатидного концентрата) потребует применения численных методов.
Рас. 2. Поля скоростей для противодавления (АР = -100) безнапорного движения жидкости (АР = 0) а напорного движения (АР = 100)
Рис. 3. Поле скоростей для противодавления (АР = -100 Па/м)
Литература:
1. Романков П.Г., Курочкина М.И. Гидромеханические процессы химической технологии. - 3-е изд., перераб. - Л.: Химия, 1982. - 288 с.
2. Олдройд Дж.Г. Неньютоновское течение жидкостей и твердых тел //Реология, теория и приложения / Под ред. Эйриха Ф.М. М.: ИЛ. - 1962. - 763 с.
