Численное моделирование стационарного течения жидкости Балкли - Гершеля в канале с внезапным расширением
Исследуется установившееся течение неньютоновской жидкости в плоском канале с внезапным расширением. Реологическое поведение среды описывается моделью Балкли - Гершеля. Для нахождения стационарных полей скорости и давления используется численный алгоритм, в основе которого лежат метод установления и процедура SIMPLE. С целью ускорения расчета применяется технология параллельных вычислений MPI. Приводятся распределения основных характеристик потока с выделением зоны двумерного течения в области скачка сечения. Показано влияние основных параметров задачи на картину распределения квазитвердых ядер в потоке жидкости.
Numerical simulation of the steady-state Herschel - Bulkley fluid flow in a channel with sudden expansion.pdf Течения реологически сложной жидкости в каналах с различными геометрическими особенностями реализуются в технологических процессах во многих отраслях промышленности. В частности, течения в каналах с внезапным расширением или сужением широко распространены в технологии формования полимерных изделий на этапах транспортировки. Сложное реологическое поведение многих полимерных композиций сильно осложняет математическое моделирование их течений. За последнее десятилетие было выполнено множество теоретических [1-8] и экспериментальных [9-13] исследований течений неньютоновских жидкостей в каналах различной геометрии. Авторами [1-3] численно решались задачи о течении вязкопластичных жидкостей в различных системах. Исследование течений ньютоновских, степенных, вязкопластичных и вязкоупругих жидкостей в каналах с внезапным сужением/расширением проводилось в работах [4-8]. Обзор экспериментальных данных по исследованиям жидкостей с пределом текучести приведен в [9]. В данном случае исследуется ламинарное стационарное течение жидкости Балкли - Гершеля в канале с внезапным расширением. Сформулированная задача решается численно с помощью алгоритма SIMPLE, при этом для обеспечения устойчивого расчета в областях малых скоростей деформаций используется регуляризация реологической модели. Постановка задачи Рассматривается установившееся ламинарное течение неньютоновской жидкости в плоском канале с внезапным расширением. Жидкость считается несжимаемой. Течение описывается системой, состоящей из уравнений движения и неразрывности, записанных в безразмерной форме „ . dUU dVU Re I -- + dp J d3U 5 2U ^ = -^- + B dx ax2 dy dx dy SB dU dB (dU 5V, + 2--+-|-+-|; (1) dx dx dy ^ dy dx „ . dUV dVV Re I-+ = -d-P+B dy dx dy ( d 2V d 2V ^ dx2 dy2 „ dB 5V dB ( 5V 5U \ + 2--+-|-+-|; (2) dy dy dx \dx dy dU 5V „ -+-= 0. (3) 5x dy Реологические свойства жидкости подчиняются закону Балкли - Гершеля, согласно которому выражение для эффективной вязкости B определяется формулой [14] (dU 5V + |-+ - ^ 5y dx Se+Ak A dU dx dV_ ~dy A = B = + 2 (4) Здесь: U и V - поперечная и продольная проекции вектора скорости u на оси декартовой системы координат (x, y), p - давление, k - степень нелинейности, A - U 02-klk безразмерная интенсивность тензора скоростей деформаций, Re = р д - безразмерный параметр вязкоплат lk Se = V обобщенное число Рейнольдса, ^U 0 стичности, р - плотность жидкости, то - предел текучести, д ческого закона. В качестве масштабов длины, скорости и давления используются величины: l - полуширина входного сечения, U0 - среднерасходная скоконстанта реологиy I Г3 рость во входном сечении, комплекс д ветственно. Область течения схематично показана на рис. 1, отношение поперечных размеров входного и выходного каналов равняется 2. Жидкость подается через входное сечение Г2 с постоянным расходом, при этом профиль скорости совпадает с профилем, характерным для установившегося течения рассматриваемой жидкости в плоском бесконечном канале. В выходном сечении Г3 используются «мягкие» граничные условия для продольной скорости и равенство нулю поперечной скорости. На твердых границах Г4 выполняются условия прилипания, на плоскости симметрии -условия симметрии. Входная и выходная границы находятся на достаточном удалении от уступа во избежание влияния последнего на характер течения в окрестностях этих границ. соот- Г t Г t Согласно такой постановке задачи, граничные условия записываются в виде dV Г1: - = 0,U = 0; dx Г2: V = V1(x), U = 0; dV Г3: - = 0,U = 0; dy Г4: V = 0,U = 0. Профиль скорости во входном сечении V1(x) находится из решения одномерной задачи об установившемся течении жидкости Балкли - Гершеля в плоском бесконечном канале с заданным постоянным расходом и имеет вид ' 1 k k+1 (-C - Se) k при x < x1; C k +1 1 k (9) при x > x1, C k +1 V1 (x ) = k+1 k+1 (-Cx- Se)~ -(-C - Se)~F" где x1 = -Se / C - уравнение плоскости, разделяющей области квазитвердого и сдвигового течений; C < 0 находится из решения алгебраического уравнения k+1 k (-C - Se )Т" [k Se-(k + 1)C]-C2 (k + 1)(2k +1) = 0. Распределения скорости во входном сечении для значений Se = 0 и Se = 1 в зависимости от степени нелинейности представлены на рис. 2. (5) (6) (7) (8) V1(x) 1.4 -1.2 -1 0.8 - 0.6 0.4 0.2 0 0 7 0.2 0.4 0.6 0.8 V1(x)-| 1.4 1.2 1 0.8 0.6 0.4 0.2 0 0 0.2 0.4 0.6 0.8 x x Рис. 2. Профили скорости Vi(x) при Se = 0 (а) и Se = 1 (б). 1 - k = 0.5, 2 - k = 1, 3 - k = 1.5 Таким образом, решение поставленной задачи сводится к отысканию стационарных полей скорости и давления, удовлетворяющих системе уравнений (1) - (4) и граничным условиям (5) - (9), а также линий тока, иллюстрирующих картину течения. Метод решения Задача решается численно с помощью конечно-разностного метода. Для нахождения стационарных полей скорости и давления в расчетных узлах разнесенной сетки используется метод установления [15], который подразумевает добавление производных искомых функций по времени в уравнения (1) и (2) соответственно, с последующей записью системы (1) - (4) в разностном виде и её решением с применением процедуры SIMPLE [16]. С целью ускорения расчета применяется технология параллельных вычислений MPI. Рассматриваемая реологическая модель обладает особенностью «бесконечной» вязкости при стремлении к нулю второго инварианта тензора скоростей деформаций. С целью обеспечения сходимости методики расчета применяются различные способы регуляризации реологической модели [17]. Суть регуляризации состоит в том, что в реологический закон вводится малый параметр, который ограничивает рост вязкости в областях малых скоростей сдвига, слабо изменяя ее величину в остальной части потока, например: 1. B = (Se + (A + е )k)j (A + e) - простая регуляризация; (V A2 + e2 )k y(V A2 + e2) - модель Беркавиер - Энглмэна (Bercovier - Engleman model). Рис. 3 иллюстрирует влияние малого параметра на профиль вязкости, рассчитанный для случая установившегося течения жидкости в плоском бесконечном канале. Сплошная линия соответствует аналитическому решению, пунктирная -простая модель, штрих-пунктир - модель Беркавиер - Энглмэна. Предлагаемые модификации, допуская предельный переход при е ^ 0 к модели Балкли - Гершеля, 2. B = | Se + B-1614 12 10 8 6 4 2 0 ~~г- 0.6 0.2 0.4 0.8 0 Рис. 3. Распределение вязкости при Se = 0, k = 0.5 и е = 10 обеспечивают возможность сквозного расчета течений с наличием квазитвердых ядер или застойных зон. Выбирая величину е заведомо большей ошибок аппроксимации, но достаточно малой для того, чтобы не исказить характер течения, можно сгладить профили эффективной вязкости в областях малых скоростей деформаций и, в то же время, получить решения, близкие к решениям с использованием исходной модели. Параметрические исследования показали, что лучшую сходимость и скорость расчета обеспечивает модель Беркавиер - Энглмэна с использованием значения параметра е ~ 10-2. В качестве условия выделения зон квазитвердого движения вязкопластичных жидкостей используется неравенство BA < Se, которое является безразмерным аналогом условия т < т0. Решение рассматриваемой задачи подразумевает формирование профиля скорости в сечении Г3, соответствующего решению одномерной задачи об установившемся течении жидкости в плоском бесконечном канале. В таблице представлены значения максимальной относительной ошибки рассчитанного профиля скорости в выходном сечении по сравнению с аналитическим решением для различных реологических параметров в зависимости от шага сетки h. Из таблицы видно, что наблюдается аппроксимационная сходимость численной методики. Все дальнейшие расчеты проводились на сетке с шагом 1/40. Максимальная ошибка по скорости в выходном сечении, % h 1/10 1/20 1/40 Ае (Se = 0, k = 0) ньютоновская жидкость 0.18 0.12 0.11 Ае (Se = 0, k = 0.8) псевдопластичная жидкость 0.23 0.13 0.11 Ае (Se = 1, k = 1) вязкопластичная жидкость 39.86 7.44 2.85 Результаты На рис. 4 показана типичная картина течения ньютоновской жидкости. Безразмерная вязкость в этом случае во всей области равна единице. По характеру распределения представленных характеристик поток можно условно разделить на три области: область одномерного движения вблизи входной границы, соответствующую течению в плоском бесконечном канале; область двумерного движения с образованием циркуляционной зоны в окрестности уступа; область одномерного движения в окрестности выходной границы. Длина области течений, представленных на рисунках, практически полностью охватывает двумерную зону, возмущенную за счет скачка сечения. Размеры этой зоны зависят от параметров задачи. Расчеты показали, что при увеличении числа Рейнольдса продольный размер циркуляционной зоны L увеличивается (рис. 5), что объясняется усилением инерционных эффектов. При этом, как видно из графика на рис. 6, данная зависимость удовлетворительно аппроксимируется линейной функцией в рассматриваемом диапазоне изменения Re, как в случае ньютоновского поведения жидкости, так и в случае течения реологически сложной жидкости (псевдопластичной k < 1, дила-тантной k > 1). 14 12 0.5 4 0.3 0.1 2 0.1 0 0.068 0.044 14- 14 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 1.5 12- 12 12 0.028 10 0.012 0.004 0 -0.004 4 -0.02 -0.044 -0.06 1.4 1.2 10 1 8 0.8 0.7 6 10- 10 8 8 8 6 6 6 4 4 2 2 2 0 0 0 0 0 1 2 0 1 2 0 1 2 Рис. 4. Распределения характеристик течения при Re = 20, k = 1 и Se = 0: а - линии тока, б - поле продольной скорости, в - поле поперечной скорости, г - поле давления Рис. 5. Структура течения в окрестности циркуляционных зон при Se = 0 и k = 1: а - Re = 1, б - Re = 10, в - Re = 20, г - Re = 40 10 987654321 0 L, усл.ед. 54- 2- 1 - 3 10 20 30 ~I Re Рис. 6. Зависимость длины циркуляционной зоны от числа Рейнольдса при Se = 1: 1 - k = 1, 2 - k = 0.6, 3 - k = 1.4 0 На рис. 7 представлена картина течения для псевдопластичной жидкости, качественно совпадающая с таковой для ньютоновской жидкости. Видно, что в центре циркуляционной зоны и на плоскости симметрии (в областях одномерного течения) локализуются зоны повышенной вязкости. Кинематический характер течения качественно совпадает с ньютоновским случаем. Усиление псевдопластичных свойств приводит к росту интенсивности поперечного движения в области двумерного течения, наряду с ее увеличением. 0.100 0.090 0.080 0.070 0.060 0.050 0.040 0.030 0.020 0.010 0.004 0.002 0.000 -0.005 -0.010 0 I 16 14 12 0 1 2 16- 16 14 12 10 8 6 4 16 1.3 1.2 1.1 1 0.9 0.8 0.7 0.65 0.62 0.55 0.4 0.2 0 -0.05 0 14121086 4 2 0 14 12 10 8 6 4 2 10 8 1 2 10.2 2-Щ - И 0 1 ■ 0 0 1 2 1 2 0 1 2 Рис. 7. Картина течения при Re = 20, k = 0.4 и Se = 0: а - линии тока, б - поле продольной скорости, в - поле поперечной скорости, г - поле вязкости, д - поле давления Влияние степени нелинейности жидкости на структуру течения исследовалось в диапазоне значений 0.4 < k < 2 и показано на рис. 8. На рис. 9 представлена аппроксимирующая кривая длины циркуляционной зоны L от степени нелинейности. Рис. 8. Структура течения в области двумерного течения: а - k = 0.4, б - k = 0.7, в - k = 1, г - k = 1.5, д - k = 2 10987654321 0 012012012012012 L, усл.ед. _ 321 - Рис. 9. Зависимость длины циркуляционной зоны от степени нелинейности при Re = 20 и Se = 0 0 0.4 0.8 1.2 1.6 k На рис. 10 показаны распределения продольной скорости и вязкости в сечении канала, проходящем через центр циркуляционной зоны для различных степеней нелинейности. Для дилатантной жидкости с ростом k кинематика течения в области двумерного течения меняется слабо. В случае течения псевдопластичной среды с уменьшением k в области расширения наблюдается значительный рост значений вязкости, наряду с уменьшением величины скорости. При k = 0.4 жидкость практически покоится. Зависимость длины циркуляционный зоны от k имеет нелинейный характер. Рис. 10. Распределение скорости (а) и вязкости (б) в сечении, проходящем через центр циркуляционной зоны: 1 - k = 0.4, 2 - k = 0.7, 3 - k = 1.0, 4 - k = 1.5 Далее проводились исследования течений жидкости Шведова - Бингама, для которой характерно образование квазитвердых ядер в областях малых скоростей деформаций. На рис. 11 показаны квазитвердые ядра и линии тока для различных параметров вязкопластичности Se. Местоположение квазитвердых ядер соответствует разделению потока на зоны одномерного и двумерного течений. С увеличением Se наблюдаются формирование застойной зоны в окрестности скачка сечения и увеличение размеров квазитвердых ядер в областях одномерного течения, при этом продольные размеры зоны двумерного движения уменьшаются. При этом, начиная со значения Se = 1, дальнейшее его увеличение практически не влияет на размеры застойной зоны. Регуляризация реологической модели позволяет проводить расчеты без явного выделения границы ядер, что в данном случае приводит к формированию ненулевого поля скоростей в застойной зоне. Однако максимальное значение модуля скорости не превышает 0.002. 15 ^^^ | | j ^ 10- -- -- -- -- -- -9- -- -------------- ± Г:: HFRHF 3- а - - б - - в - - г - -II д - -II е - 2- -- -- -- -- -- 1 --Ill .......... H I .......... 1 1 1 1 1 1 J 1 J 1 F 1 1 1 1 а б в г е 1 1 1 1 д 012 012 012012 012 012 Рис. 11. Распределение линий тока и квазитвердых ядер при Re = 20 и к = 1: а - Se = 0.1, б - Se = 0.3, в - Se = 0.5, г - Se = 1.0, д - Se = 1.5, е - Se = 2.0 0 На рис. 12 показаны размеры и местоположение квазитвердых ядер для различных чисел Рейнольдса. Наблюдается увеличение зоны двумерного течения и продольных размеров застойной зоны при увеличении Re. Влияние параметра нелинейности k в случае течения жидкости с пределом текучести показано на рис. 13. 0 12 0 12 0 12 Рис. 12. Зависимость формы квазитвердых ядер Рис. 13. Зависимость формы квази-от Re при k = 1 и Se = 1: а - Re = 1, б - Re = 10, твердых ядер от k при Re = 20 и в - Re = 20, г - Re = 40 Se = 1: а - k = 0.8, б - k = 1, в - k = 1.2 201918171615141312 11 10 9 8 7 6 5 4 3 2 1 0 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 i i 0 1 2 12 0 12 Заключение Сформулирована математическая постановка задачи о стационарном течении неньютоновской жидкости в канале с внезапным расширением и предложена численная методика ее решения. Показано наличие в потоке трех зон: одномерного движения в окрестности входа; одномерного движения в окрестности выхода; двумерного движения в области скачка сечения, которое характеризуется наличием циркуляционной зоны. Продемонстрирована линейная зависимость продольных размеров последней в зависимости от числа Рейнольдса в диапазоне 1 < Re < 40. В случае течения степенной жидкости зависимость L(k) имеет нелинейный характер при изменении k в пределах от 0.4 до 2. Для течений жидкостей с пределом текучести показано наличие квазитвердых ядер вблизи плоскости симметрии в зонах одномерного движения и застойных участков в области скачка сечения. Исследовано влияние параметра вязкопластичности в диапазоне 0 < Se < 5 на картину течения.
Ключевые слова
течение жидкости,
канал с внезапным расширением,
численное моделирование,
модель Балкли - Гершеля,
квазитвердые ядра,
fluid flow,
channel with sudden expansion,
numerical simulation,
Herschel-Bulkley model,
dead zoneАвторы
Борзенко Евгений Иванович | Томский государственный университет | кандидат физико-математических наук, доцент кафедры прикладной газовой динамики и горения | borzenko@ftf.tsu.ru |
Хегай Ефим Игоревич | Томский государственный университет | студент физико-технического факультета | efim_h@ftf.tsu.ru |
Всего: 2
Ссылки
Frigaard, C. Nouar. On the usage of viscosity regularisation methods for visco-plastic fluid flow computation // J. Non-Newtonian Fluid Mechanics. 2005. No. 127. P. 1-26.
Патанкар С. Численные методы решения задач теплообмена и динамики жидкости. М.: Энергоатомиздат, 1984. 152 с.
Годунов С.К., Рябенький В.С. Разностные схемы. М.: Наука, 1977. 440 с.
Смольский Б.М., Шульман З.П., Гориславец В.М. Реодинамика и теплообмен нелиней-но-вязкопластичных материалов. Минск: Наука и техника, 1970.
Maillard M., Boujlel J., Coussot P. Flow characteristics around a plate withdrawn from a bath of yield stress fluid // J. Non-Newtonian Fluid Mechanics. 2015. No. 220. P. 33-43.
Ovarlez G., Cohen-Addad S., Krishan K., Goyon J., Coussot P. On the existence of a simple yield stress fluid behavior // J. Non-Newtonian Fluid Mechanics. 2013. No. 193. P. 68-79.
Stephane Mossaz, Pascal Jay, Albert Magnin. Experimental study of stationary inertial flows of a yield-stress fluid around a cylinder // J. Non-Newtonian Fluid Mechanics. 2012. No. 189-190. P. 40-52.
Perez-CamachoM., Lopez-Aguilar J.E., Calderas F., Manero O., WebsterM.F. Pressure-drop and kinematics of viscoelastic flow through an axisymmetric contraction - expansion geometry with various contraction-ratios // J. Non-Newtonian Fluid Mechanics. 2015. No. 222. P. 260-271.
Coussot P. Yield stress fluid flows: A review of experimental data // J. Non-Newtonian Fluid Mechanics. 2014. No. 211. P. 31-49.
Lober Hermany, Daniel Dall'Onder dos Santos, Sergio Frey, Monica F. Naccache, Paulo R. de Souza Mendes. Flow of yield-stress liquids through an axisymmetric abrupt expansion-contraction // J. Non-Newtonian Fluid Mechanics. 2013. No. 201. P. 1-9.
Dhinakaran S., Oliveira M.S.N., Pinho F.T., Alves M.A. Steady flow of power-law fluids in a 1:3 planar sudden expansion // J. Non-Newtonian Fluid Mechanics. 2013. No. 198. P. 48-58.
Primo"z Ternik. Planar sudden symmetric expansion flows and bifurcation phenomena of purely viscous shear-thinning fluids // J. Non-Newtonian Fluid Mechanics. 2009. No. 157. P. 15-25.
Fernanda B. Link, Sergio Frey, Roney L. Thompson, Monica F. Naccache, Paulo R. de Souza Mendes. Plane flow of thixotropic elasto-viscoplastic materials through a 1:4 sudden expansion // J. Non-Newtonian Fluid Mechanics. 2015. No. 220. P. 162-174.
Ferras L.L., Afonso A.M., Alves M.A., Nobrega J.M., Carneiro O.S., Pinho F.T. Slip flows of Newtonian and viscoelastic fluids in a 4:1 contraction // J. Non-Newtonian Fluid Mechanics. 2014. No. 214. P. 28-37.
Yiolanda Damianou, Georgios C. Georgiou, Irene Moulitsas. Combined effects of compressibility and slip in flows of a Herschel - Bulkley fluid // J. Non-Newtonian Fluid Mechanics. 2013. No. 193. P. 89-102.
Patel S.A., Chhabra R.P. Steady flow of Bingham plastic fluids past an elliptical cylinder // Journal of Non-Newtonian Fluid Mechanics. 2013. No. 202. P. 32-53.
Ying-Hsin Wua, Ko-Fei Liu. Start-up flow of a Bingham fluid between two coaxial cylinders under a constant wall shear stress // J. Non-Newtonian Fluid Mechanics. 2015. No. 223. P. 116-121.