CC BY
32
УДК 532.135
Д.В. Ананьев, Е.К. Вачагина, Г.Р. Халитова
ИССЛЕДОВАНИЕ ГИДРОДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК ПОТОКА ПРИ ЛАМИНАРНОМ ТЕЧЕНИИ РЕОЛОГИЧЕСКИ СЛОЖНОЙ ЖИДКОСТИ
В КОНВЕРГЕНТНОМ КАНАЛЕ
D.V. Ananyev, E.K. Vachagina, G.R. Halitova
INVESTIGATION OF HYDRODYNAMIC CHARACTERISTICS OF LAMINAR FLOW FLUID WITH TRANSVERSE VISCOSITY IN THE CONVERGENT CHANNEL
Исследуется гидродинамика при ламинарном течении водного раствора гидроксиэтилцеллю-лозы, которая представляет собой нелинейно-вязкую жидкость с поперечной вязкостью, в конвергентном канале. В качестве реологической модели используется нелинейная модель, определяющая жидкость Рейнера — Ривлина.
ГИДРОДИНАМИКА. ГИДРОКСИЭТИЛЦЕЛЛЮЛОЗА. НЕЛИНЕЙНО-ВЯЗКАЯ ЖИДКОСТЬ. ПОПЕРЕЧНАЯ ВЯЗКОСТЬ. КОНВЕРГЕНТНЫЙ КАНАЛ. ЖИДКОСТЬ РЕЙНЕРА-РИВЛИНА.
We study the hydrodynamics in the laminar flow of an aqueous solution of hydraxyethyl cellulose, is a non-linear viscous fluid with a viscosity shear, in a convergent channel. A nonlinear rheological model defining the Reiner — Rivlin fluid is used.
HYDRODYNAMICS. HYDROXYETHYL CELLULOSE. NONLINEAR VISCOUS FLUID. TRANSVERSEVISCOSITY. CONVERGENT CHANNEL. REINER — RIVLIN FLUID.
Широкое применение реологически сложных сред в химической, нефтехимической и пищевой промышленности обусловливает возрастающий интерес к изучению гидродинамики потока в сложных каналах. К числу жидкостей Рейнера — Ривлина относятся коллоидные растворы, высокополимерные соединения, некоторые геоматериалы.
Как известно, наличие конвергентных и дивергентных частей канала позволяет значительно интенсифицировать процессы теплообмена. Наряду с другими причинами интенсификации, большое значение имеет изменение конфигурации профиля составляющих вектора скорости в соответствии с наличием нелинейно-вязких эффектов поведения жидкостей Рейнера—Ривлина. Именно наличие поперечной вязкости определяет особенности течения в каналах со сложной геометрией. И это диктует необходимость исследования течений такого рода жидкостей.
Известно множество работ, посвященных течению сред в каналах со сложной геометрией.
В работе [1] численно было исследовано полностью развитое течение внутри двумерного периодически сужающегося-расширяющегося волнистого канала с пульсирующим потоком ньютоновской жидкости. Решение задачи теплообмена при течении вязкой несжимаемой жидкости между двумя непараллельными пластинами рассмотрено в [2]. Профиль скорости определяется с помощью эллиптических функций Якоби, а расчет температурных полей осуществляется численно. Аналитическое решение задачи Джеффри — Хамеля для случая скольжения на границах канала представлено в [3]. Проведено сравнение со случаем прилипания на границах.
Поведение жидкости с поперечной вязкостью при течении в каналах, имеющих сложную геометрию, имеет особенности, отличающие его от поведения ньютоновских сред. Однако в литературе практически нет работ, описывающих эти особенности. Поэтому в данной статье численно решается задача течения жидкости с поперечной вязкостью в конвергентном канале.
Постановка задачи
При построении математической модели были приняты следующие допущения:
течение среды — установившееся, ламинарное, со сформировавшимся профилем вектора скорости на входе в канал;
реологические свойства среды характеризуются наличием нелинейно-вязких эффектов;
плотность среды в ходе процесса меняется незначительно;
силы тяжести пренебрежимо малы. Геометрическая область течения представляет собой конвергентный канал, образованный наклонными пластинами, расстояние между которыми много меньше ширины канала. Тогда можно предположить, что распределение основных гидродинамических характеристик в любом продольном сечении не меняется. Все вышесказанное позволяет свести задачу к двумерной постановке. Введем декартову систему координат, ось 0Х которой направлена по оси симметрии канала в направлении течения жидкости. Продольное сечение такой геометрической области, представлено на рис. 1.
Для описания реологических свойств жидкости будем использовать нелинейную модель Рейнера — Ривлина [4, 5]:
Т = _р1 + ф (/2,Т )В + Ф2 (12,Т )В2;
В = 2Б; Б = 2(gradv + <га(у )Т).
12 = 4^ (Б2 ) = ^ (В2),
(1) (2)
(3)
где р — давление; I — единичный тензор; Б — тензор скоростей деформаций; V — вектор скорости; ф1 (12), ф2 (12) — материальные функ-
ции; 12 — второй инвариант тензора скоростей деформаций; индекс «Т» — символ транспонирования; <4г» — операции взятия следа тензора.
Математическая модель рассматриваемой задачи получена на основе уравнений механики сплошных сред и имеет вид
Эух Эух л
V -— + V -—
Ух дх Уу ду
Эр + д_Г
дх дх
2Ф1 <12 £
+
+-
ду
Ф1 (12 )
'дух+дУу
Эу дх
V
Эу
+
(у<2 )); (4)
УУ
Эу л
V -^ + V -у
Ух дх Уу
ду дх
Ф1 (12 )
дУ У
'дух+д!у ^
ду дх
+
УУ
+
ду
/ ч Эуу 2ф1 (12
эу
\\
+
(у (2 ))
У У
(5)
(6)
Эух + дуу = о дх ду
где р — плотность жидкости; (у(ф2В2)) и (у (ф2В2)) — декартовы компоненты вектора &у (ф2В2); ух и уу — компоненты вектора скорости.
Эллиптический характер используемых уравнений математической модели диктует необходимость постановки граничных условий на входе и выходе из канала. Рассмотрим вместо действительной геометрической области определения переменных задачи (ЛВСБ на рис. 1)
К Уi 1 Л
Н
Ь
х1 — ■
Б
Рис. 1. Геометрическая область исследуемого конвергентного канала:
Н — наибольшая высота канала; Ь — длина боковой поверхности канала; а — угол наклона пластин, образующих канал; х1 — длина конвергентной части канала
Ь
х
видоизмененную теоретическую область, полученную из данной увеличением длины канала (плоские границы) до входа (AKи DL на рис. 1) и после выхода (BE и CF на рис. 1). Причем BE и CF имеют длину, при которой происходит стабилизация скоростных полей. На входе в канал (LK) задается распределение скорости, являющееся решением стационарной задачи течения в канале, образованном двумя плоскостями, расстояние между которыми равно Н (H — наибольшая высота канала). На стенках принимаются гидродинамические граничные условия прилипания — v = 0 . На продолженных границах принимаются гидродинамические граничные условия прилипания — v = 0 . В качестве граничных условий на выходе (EF) принимаются граничные условия стабилизации искомых полей по длине: dvx/dx = 0, vy = 0 .
Метод решения
Для решения поставленной задачи использовался метод конечных элементов, реализованный в коммерческом пакете Femlab. Разработанная модель была описана с помощью опций Femlab и введения дополнительных членов, связанных с наличием поперечной вязкости, что привело к созданию новой пользовательской модели в рамках программы Femlab. Для решения задач гидродинамики при течении жидкости с поперечной вязкостью использовались нелинейные элементы. Метод решения применяемой программы основан на использовании интегральных соотношений метода Бубнова — Га-лёркина.
Тестирование метода решения проводилось путем сравнения результатов решения задачи течения гидроксиэтилцеллюлозы в цилиндрическом канале квадратного поперечного сечения, полученного в [6] В.Г. Литвиновым и в данной работе (рис. 2). Разница между распределением продольной составляющей вектора скорости vz на половине диагонали квадрата, которое получено В.Г. Литвиновым (кривая 1), и распределением, которое получено методом, использованным в данной работе (кривая 2), составляет не более 4 %.
Также было проведено сравнение результатов решения задачи течения ньютоновской жидкости в конвергентном каналеиз работы [7] с решением, которое получено с использованием
vz, м/с
0
0,004
0,008
r, м
Рис. 2. Распределения у2: 1 — получено В.Г. Литвиновым, 2 — получено авторами данной работы
программного пакета Реш1аЪ (рис. 3, б). Рассматривалось течение воды в конвергентном канале (а = 0,01). Схема течения и расположение осей цилиндрической системы координат задачи, рассмотренной в [7], изображены на рис. 3. Кривая 1 на рис. 3, б — из работы [7], кривая 2—решение, полученное с помощью программного пакета Реш1аЪ.
Результаты
Для численных расчетов был выбран водный раствор гидросиэтилцеллюлозы, эффективная и поперечная вязкость которой определяются следующим образом [6]:
Ф1 (12) = 2 при 0<7^2<4; 4
ф1 (12 ) = ущ при 4 <VV2 < 103; Ф2(12) = -0,32 при 0<VV2 < 1;
при 1 <VV2
(7)
Ф2 (2 ) =
-0,32
(8)
(12/2 )
3/4
Конвергентный канал можно однозначно определить, задав угол а наклона плоскостей, образующих канал, длину плоскостей Ь и расстояние Н между плоскостями на входе в канал (рис. 1). Как следует из уравнений вязкости (7), (8), гидродинамическая картина течения жидкости с поперечной вязкостью в конвергентном канале зависит от изменения вектора скорости.
Рис. 3. Схема течения, принятая в [7], (а) и распределения уг/угтах (б): уг — радиальная составляющая вектора скорости; уг тах — ее максимальное значение
Распределение составляющих вектора скорости представлено на рис. 4. Расчет был проведен для канала со следующими геометрическими характеристиками: а =0,2; H|Ь =0,8. Режим течения определялся расходом, равным 0,029 м2/с , что соответствовало значению числа Яе=57 на входе в канал, т. е. в самой широкой части канала.
Из рис. 4 видно, что наибольшие абсолютные величины максимумов поперечной и продольной составляющих вектора скорости — в узкой части канала.
Изменение второго инварианта тензора скоростей деформации, от которого явно зависят эффективная и поперечная вязкости, представ-
лено на рис. 5, a. Распределения !2 показывают, что изменение !2 в узкой части канала (кривая 1) характеризуется наибольшим значением !2 на границе, а дальше монотонно убывает до нуля на оси симметрии канала.
Распределение функции Ф2 в поперечном сечении канала представлено на рис. 5, б. Из рисунка видно, что абсолютное значение поперечной вязкости изменяется от своего наибольшего значения в центре канала до наименьшего на границе, при этом в узкой части наименьшее значение Ф2 меньше, чем в более широкой части, что обусловлено более высоким значением !2 (кривая 1).
а)
ух 0,20
б)
-0,004
-0,008
-0,012
0,04 0,08 0,12 0,16 0,20 у
0,04 0,08 0,12 0,16 0,20 у
0
V
у
0
Рис. 4. Распределение составляющих V (а) и V (б) вектора скорости в поперечном сечении канала:
х у
а) 1 — х = х{; 2 — х = 2хх/3; 3 — х = хх/3; 4 — х = 0; б) 1 — х = 5хх/6; 2 — х = 2хх/3; 3 — х = хх/2; 4 — х = хх/3; 5 — х = хх/6
а)
к
50 40 30 20 10 0
б)
0 0,04 0,08 0,12 0,16 0,20 у
0 0,04 0,08 0,12 0,16 0,20 у
Рис. 5. Распределение второго инварианта /2 тензора скоростей деформации (а) и поперечной вязкости ф2 (б) в поперечном сечении канала: 1 — х = х1; 2 — x = 2x1/3 ; 3 — x = x1/3; 4 — x = 0
а)
-0,004
-0,008
-0,012
0,04 0,08 0,12
б)
0 0,04 0,08 0,12
-0,1-
у
у
0
у
у
Рис. 6. Распределение поперечной составляющей вектора скорости (а) и поперечной вязкости (б) при х = 2х1у/3 в зависимости от числа Яе: 1 — Яе =1 5; 2 — Яе = 38; 3 — Яе = 57
На рис. 6, а представлено распределение поперечной составляющей скорости в поперечном сечении канала при х = 2 хх/3 в зависимости от числа Яе. Как видно, с увеличением числа Яе происходит рост поперечной составляющей скорости. Расчеты проведены для канала со следующими геометрическими характеристиками: а =0,2; Ь=0,5; Н/Ь =0,8.
Режим течения, на который существенное влияние оказывает расход, определяет в конечном итоге изменение поперечной вязкости, что представлено на рис. 6, б. При увеличении числа Рейнольдса происходит уменьшение абсо-
лютного значения поперечной вязкости, что также обусловлено увеличением величины 12 вблизи границы канала с увеличением Яе и псевдопластичными свойствами рассматриваемой жидкости.
Представленные результаты численных расчетов позволяют описать особенности течения, возникающие из-за наличия поперечной вязкости и различных размеров конвергентного канала.
Характер течения жидкости с поперечной вязкостью существенно зависит от геометриче-
ских параметров конвергентного канала, а также гидродинамического режима, определяемого числом Яе. При увеличении числа Рейнольдса происходит рост поперечной составляющей вектора скорости. С увеличением угла наклона растет абсолютное значение составляющих вектора
скорости. Абсолютное значение поперечной вязкости уменьшается с уменьшением отношения Н/Ь.
Работа выполнена при финансовой поддержке РФФИ и Академии наук РФ (грант № 12-08-97034-р_поволжье_а).
СПИСОК ЛИТЕРАТУРЫ
1. Gong-Nan Xie. Heat transfer enhancement inside converging-diverging wavy channel by pulsating flow [Текст] / Xie Gong-Nan, Wang Qiu-Wang, Zeng Min [et al.] // Journal of Chemical Engineering of Chinese Universities.— 2006.— Vol. 20, № 1.— Р. 31-35
2. Millsapst, K. Thermal Distributions in Jeffery-Hamel Flows Between Nonparallel Plane Walls [Текст] / K. Millsapst, K. Pohlhausent // Journal of the Aeronautical Sciences (Institute of the Aeronautical Sciences.— Vol. 20, Is. 3.
3. Dorrepaal, J.M. Slip flow in converging and diverging channels [Текст] / J.M. Dorrepaal // Journal of Engineering Mathematics.— 1993.— Vol. 27.— Р. 343-356.
4. Брутян, М.А. Гидродинамика неньютоновских
жидкостей [Текст] / М.А. Брутян, П.Л. Крапивский // Итоги науки и техники. Сер.: Комплексные и специальные разделы механики.— М.: ВИНИТИ, 1991.— Т. 4.— С. 3-98.
5. Astarita, G. Principles of non-newtonianfluid mechanics / G. Astarita, G. Marrucci.— McGRAW-HILL, 1974.— 309 p.
6. Литвинов, В.Г. Движение нелинейно-вязкой жидкости [Текст] / В.Г. Литвинов.— М.: Наука, 1982.— 374 с
7. Vachagina, E.K. Fourier method for heat transport equation in the convergent channel [Текст] / E.K. Vachagina, D.V. Ananyev // International Journal of Heat and Mass Transfer.— 2013.— Vol. 57.— Р 148-154.
REFERENCES
1. Xie Gong-Nan. Heat transfer enhancement inside converging-diverging wavy channel by pulsating flow [Текст] / Xie Gong-Nan, Wang Qiu-Wang, Zeng Min [et al.] // Journal of Chemical Engineering of Chinese Universities.— 2006.— Vol. 20, № 1.— Р 31-35.
2. Millsapst, K. Thermal Distributions in Jeffery-Hamel Flows Between Nonparallel Plane Walls [Текст] / K. Millsapst, K. Pohlhausent // Journal of the Aeronautical Sciences (Institute of the Aeronautical Sciences).— Vol. 20. Is. 3.
3. Dorrepaal, J.M. Slip flow in converging and diverging channels [Text] / J.M. Dorrepaal // Journal of Engineering Mathematics.— 1993.— Vol. 27.— Р 343-356.
4. Brutian, M.A. Gidrodinamika nen'iutonovskikh zhidkostei [Tekst] / M.A. Brutian, PL. Krapivskii // Itogi nauki i tekhniki. Ser.: Kompleksnye i spetsial'nye razdely mekhaniki.— M.: VINITI, 1991.— T. 4.— S. 3-98. (rus.)
5. Astarita, Dzh. Osnovy gidromekhaniki nen'iutonovskikh zhidkostei [Text] / Dzh. Astarita, Dzh. Maruch-chi.— M.: Mir, 1978.— 309 s.
6. Litvinov, V.G. Dvizhenie nelineino-viazkoi zhidkosti [Tekst] / V.G. Litvinov.— M.: Nauka, 1982.— 374 s. (rus.)
7. Vachagina, E.K. Fourier method for heat transport equation in the convergent channel [Tekst] / E.K. Vachagina, D.V. Ananyev // International Journal of Heat and Mass Transfer.— 2013.— Vol. 57.— P. 148-154.
СВЕДЕНИЯ ОБ АВТОРАХ /AUTHORS
АНАНЬЕВ Дмитрий Владиславович — кандидат технических наук научный сотрудник лаборатории ТФИ Исследовательского центра проблем энергетики Казанского научного центра Российской академии наук; 420111, ул. Лобачевского, 2/31, Казань, Россия; e-mail: dima211083@yandex.ru ANANYEV Dmitry V. — Research Center of Power Engineering Problems of RAS Institution of the Kazan Scientific Centre of RAS; 420111, P.B. 190, 2/31 Lobachevsky str., Kazan, 420111, Russia; e-mail: dima211083@ yandex.ru
ВАЧАГИНА Екатерина Константиновна — доктор технических наук заведующая лабораторией ТФИ Исследовательского центра проблем энергетики Казанского научного центра Российской академии наук; 420111, ул. Лобачевского, 2/31, Казань, Россия; e-mail: evachagina@mail.ru
Энергетика. Электротехника -►
VACHAGINA Ekaterina K. — Research Center of Power Engineering Problems of RAS Institution of the Kazan Scientific Centre of RAS; 420111, P.B. 190, 2/31 Lobachevsky str., Kazan, 420111, Russia; e-mail: evachagina@mail.ru
ХАЛИТОВА Галия Рахимовна — старший научный сотрудник ЭТПИЭ Исследовательского центра проблем энергетики Казанского научного центра Российской академии наук; 420111, ул. Лобачевского, 2/31, Казань, Россия; e-mail: halitova-galia@mail.ru
HALITOVA Galia R. — Research Center of Power Engineering Problems of RAS Institution of the Kazan Scientific Centre of RAS; 420111, P.B. 190, 2/31 Lobachevsky str., Kazan, 420111, Russia; e-mail: halitova-galia@mail.ru
© Санкт-Петербургский государственный политехнический университет, 2013
