CC BY
23
2016 Математика и механика № 6(44)
УДК 532.542
DOI 10.17223/19988621/44/5
Ю.Б. Банзула, Е.И. Борзенко, С.В. Карязов, Г.Р. Шрагер
КИНЕМАТИКА ТЕЧЕНИЯ ВЯЗКОЙ ЖИДКОСТИ В КАНАЛЕ С ЗАТВОРОМ1
Представлены результаты математического моделирования течения ньютоновской жидкости в канале с конструктивным элементом, моделирующим наличие затвора. Математическая постановка задачи формулируется в переменных вихрь - функция тока. Численное решение рассматриваемой задачи осуществляется конечно-разностным методом на основе схемы переменных направлений. Представлены результаты параметрических исследований кинематики течений как для полностью открытого, так и частично закрытого затвора.
Ключевые слова: течение, вязкая жидкость, затвор, численное моделирование, кинематика течения.
Во многих технологических процессах при переработке материалов в жидком состоянии реализуется транспорт жидкости по трубопроводу, который является частью технологической оснастки. Транспортный трубопровод, как правило, содержит конструктивные элементы специального назначения (соединения типа сужение/расширение, диафрагма, клапаны, краны и т.п.). При расчете расходно-напорных характеристик трубопровода необходимо учитывать потери энергии в потоке за счет конструктивных элементов [1]. Начиная с середины прошлого столетия в этом направлении ведутся интенсивные теоретические и экспериментальные исследования как для различных режимов течения, так и для каналов разной геометрии, в том числе с использованием сложных реологических моделей поведения жидкой среды. Обзор подобных исследований представлен в работах [1-4]. Утверждается, что исследования в данной области являются актуальными и в настоящее время. Это обусловлено не только многочисленными практическими приложениями, но и тем, что имеющиеся экспериментальные данные во многих случаях не согласуются не только между собой, но и с результатами аналитических и численных исследований [5]. Широко распространенным конструктивным элементом, используемым для открытия/закрытия трубопровода, является шаровой затвор. Конструкция шарового затвора позволяет осуществлять полное и частичное перекрытие сечения трубопровода. Течение жидкости в окрестности шарового затвора имеет сложный пространственный характер и его качественное и количественное исследование возможно лишь с использованием современных вычислительных технологий, реализуемых на высокопроизводительных ЭВМ. Имеются попытки численного моделирования течений жидкости в трубопроводах с шаровыми кранами и в системах с шаровыми клапанами с помощью пакетов прикладных программ, например «FLUENT» [6, 7].
Целью данной работы является численное моделирование течения жидкости в плоскости симметрии канала с затвором в положении как полного, так и частичного перекрытия сечения.
1 Работа выполнена при финансовой поддержке РФФИ (проект № 15-08-03935-а).
Постановка задачи
Рассматривается ламинарное, стационарное течение ньютоновской жидкости в плоском бесконечном канале шириной В. Область решения в декартовой системе координат схематично изображена на рис. 1. На некотором удалении от входной границе Г1 и от выходной границы Г5 в канале расположен затвор диаметра В1, который может менять ширину проходного сечения в зависимости от положения запорного элемента Г4. Значения Ь1 и Ь2 подбираются такими, чтобы исключить влияние границ Г1 и Г5 на характер течения в окрестности затвора. Геометрия стенки Г2 описывается непрерывной функциейДх), стенки Г3 - (-Дх)).
Г5
Рис. 1. Область течения Fig. 1. Flow region
Г
X
Задача формулируется в плоской постановке с использованием декартовой системы координат, начало которой располагается во входном сечении. Течение описывается уравнениями движения и неразрывности, которые в безразмерных переменных принимают вид
( ди ди Л Яе др
Яе I и — + V— 1 =---— + Ди , (1)
\ дх ду) 2 дх
„ ( дv дv Л Яе др
Яе I и — + V— 1 =---- + Дv, (2)
^ дх ду) 2 ду
^ + ^ = о. (3)
дх ду
На твердых стенках (Г2, Г3) и запорном элементе (Г4) выполняются условия прилипания; на входной границе (Г1) жидкость поступает с постоянным расходом, распределение скорости при этом соответствует установившемуся течению жидкости в бесконечном канале с заданным постоянным расходом; на выходной границе (Г5) используются мягкие граничные условия для продольной скорости и равенство нулю поперечной. Таким образом, граничные условия записываются в виде
(Г1): v=0, и=1.5(1 - у2); (4)
(Г2, Гз, Г4): v=0, и=0; (5)
(Г5): V = 0, ^ = 0. (6)
дх
Здесь: и, V - продольная и поперечная компоненты вектора скорости соответственно; р - давление; Яе = рПЬ/ц - число Рейнольдса (£=0/2); р - плотность; ц -
вязкость. В качестве масштабов обезразмеривания выбраны следующие величины: длина - полуширина входного канала 0/2; скорость - среднерасходная скорость во входном сечении П; давление - величина рП2/2 .
Используя безразмерные функцию тока и вихрь, определяемые выражениями
ду ду ¿V ди
дх' ду' дх ду'
перепишем постановку (1) - (6) в переменных функция тока - вихрь
дю дю 1
и--+ V— = —Дю ,
дх ду Яе
V > = -ю :
(Г1): у = 1 +1.5
,.зЛ
ю = 3 у;
(7)
(8) (9)
(10)
(Г2, Г4): у = 2,= 0;
дщ
(Гз): у = 0, ^ = 0;
дщ
(Г5): ^ = 0, ^ = 0.
дх дх
(11) (12) (13)
В уравнениях (11), (12) щ соответствует направлению внешней к границе нормали. Функция Дх) в безразмерном виде записывается следующим образом:
I (х) =
Я2-\х - Ь -1
0.5
х е[/ Ь, Ь1/ Ь + 21 ],
1, х Ь, Ь1/ Ь + 21 ],
I =
Уд2 - О2
О
= 4яГ-1 .
Таким
где Я = д - безразмерный радиус затвора,
образом, решение задачи определяется критериями подобия Яе, геометрическими параметрами Я и - степень закрытия запорного элемента (рис.1).
Метод расчета
Для решения сформулированной задачи в переменных функция тока - вихрь используется метод установления [8]. Физическую область течения с криволинейной границей Дх) преобразуем в прямоугольную введением новых координат
| = х, п = у /1 (х).
В преобразованной области решения строится квадратная разностная сетка с шагом И0.И ={^ = ¡И, п = /И, I = 0,...,} = 0,...,И2}. Уравнения (8), (9) с учетом преобразования координат записываются в разностном виде с применением
схемы продольно-поперечной прогонки [9]. Конвективные слагаемые в уравнении переноса вихря аппроксимируются с помощью схемы против потока [10].
Значение вихря на стенке определяется выражением
д 2ш
® =--2Т ,
дп2
для разностного представления которого используется формула Тома [22].
Профиль функции Л(х) в малой окрестности точек Ь11Ь, Ь1 /(Ь + 21) сглаживается полиномом третьей степени, обеспечивающим непрерывность самой функции, а также ее первой и второй производных.
Для проверки аппроксимационной сходимости методики расчета была проведена серия вычислительных экспериментов на последовательности сеток. Сходимость демонстрируется по распределению скорости в поперечных сечениях х = Ь1 /Ь и х = Ь1/(Ь +1) (рис. 2).
Рис. 2. Распределение продольной скорости в сечениях x = Lj/L (а) и x = L1/(L +1) (b):
Re = 1, n = 1, R = 2 (кр. 1 - 1/10, кр.2 - 1/20, кр. 3 - 1/40, кр. 4 - 1/80, кр. 5 - 1/160) Fig. 2. Distribution of the axial velocity in cross-sections (а) x = L1 /L and (b) x = L1 /(L +1) for Re = 1, n = 1, R = 2 ((1) 1/10, (2) 1/20, (3) 1/40, (4) 1/80, (5) 1/160)
Анализ результатов показывает аппроксимационную сходимость метода. Максимальные расхождения между результатами, полученными на сетках с шагом 1/80 и 1/160, в значениях скорости составляют 0.02 %. Все дальнейшие вычисления проведены на сетках с шагом 1/80.
Результаты расчетов
Расчеты показали, что в случае установившегося течения жидкости в рассматриваемом канале с полностью открытым затвором распределения характеристик симметрично относительно плоскости >>=0, при этом можно выделить три зоны в области течения. В окрестности входной и выходной границ наблюдаются зоны одномерного течения, характерного для установившегося движения в бесконечном канале с профилем скорости (4), в окрестности затвора движение имеет двумерный характер. Поля компонент скорости в случае малых значений чисел Рей-нольдса демонстрируются на рис. 3, а распределения линий тока в зависимости от Яе представлены на рис. 4.
Рис. 4. Распределения изолиний функции тока (a,b,c,d,ef,g - Re = 0.1, 1, 5, 10, 20, 50, 80) Fig. 4. Streamline distributions for Re = (a) 0.1, (b) 1, (c) 5, (d) 10, (e) 20, (f 50, and (g) 80
Видно, что распределения линий тока и продольной скорости при малых Яе симметричны также и относительно плоскости х = хс, где хс - координата центра затвора. Увеличение числа Яе приводит к возникновению в области затвора зон циркуляционного движения. Последние формируются в левой части затвора (рис. 4, d) и по мере увеличения Яе занимают всю полость (рис. 4, g). Распределения характеристик течения в области двумерного течения при Яе = 50 представлены на рис. 5. Видно, что внутри затвора ширина эффективного сечения, через которое протекает жидкость, уменьшается, что приводит к росту продольной скорости в окрестности плоскости симметрии. Параметрические расчеты показывают, что при Яе < 1 распределения кинематических характеристик практически не меняются.
-0.13 -0.09 -0.05 -0.01 0.03 0.07 0.11
-0.2 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8
Рис. 5. Картина течения Re = 50 (а, b - распределения v, u)
Fig. 5. Flow pattern for Re = 50 (distributions of the velocity components (a) v and (b) u)
Влияние радиуса затвора на картину течения демонстрирует рис. 6, где представлены распределения линий тока при различных значениях R для двух чисел Рейнольдса. В случае Re=1 с ростом размеров затвора характер течения не меняется, линии тока плавно огибают его контур.
Частичное закрытие затвора приводит к формированию в зоне двумерного течения несимметричной картины. Распределения составляющих вектора скорости для различных значений Re в случае, когда запорный элемент перекрывает половину сечения канала, представлены на рис. 7 и 8.
Влияние степени закрытия затвора на картину течения демонстрирует рис. 9. При малых значениях Re линии тока плавно огибают контуры канала, а по мере закрытия затвора за ним формируется циркуляционная зона (рис. 9, a,b,c). Перед элементом возникает зона циркуляционного движения, размеры которой практически не меняются при изменении числа Рейнольдса. При Re=10 практически вся жидкость за затвором вовлечена в циркуляционное движение, а размеры циркуляционной зоны в нижней части потока уменьшаются по мере закрытия затвора (рис. 9, d,ef).
1.
0.502
d
Рис. 6. Распределения изолиний функции тока при различных R (a,b,c - Re = 1, d,ef - Re = 10; a,d - R = 1.6, b,e - R = 1.8, cf - R = 2.0) Fig. 6. Streamline distributions for various R (Re = (a,b,c) 1, and (d,ef 10; R = (a,d) 1.6, (b,e) 1.8, and (cf 2.0)
12
1
1
-1.1 -0.9 -0.7 -0.5 -0.3 -0.1 0.2 0.4 0.6
-0.1 0.3 0.7 1.1 1.5 1.9 2.3 2.7
-i-1-1— i-1-1— -2"l-1-1-1 ' I1 -1-r-
7 8 9 10 11 12 6 7 8 9 10 11 12
Рис. 7. Картина течения Re = 1, s = 0.5 (a, b - распределения v, u) Fig. 7. Flow pattern for Re = 1, s = 0.5 (distributions of (а) v and (b) u)
-1.4-1.2 -1 -0.8-0.6-0.4-0.2 0.1 0.3 0.5 0.7 0.9
-0.6 -0.2 0.2 0.6 1 1.4 1.8 2.2 2.6 3 3.4
a ---\
tllip
'6 7
1 01 -1
I I I-1-1 -2--г-T-J I -1-1-
9 10 11 12 13 26 7 8 9 10 11 12 13
Рис. 8. Картина течения Re = 10, s = 0.5 (a, b - распределения v, u) Fig. 8. Flow pattern Re = 10, s = 0.5 (distributions of the velocity components (a) v and (b) u)
2
2 1
0
-1
-2
2
1
0
-1
-2
2
1
0 -1 -2
6 7 8 9 10 11 12 6 7 8 9 10 11 12
Рис. 9. Распределения изолиний функции тока при различных s (a,b,c - Re = 1, d,ef- Re = 10; a,d - s = 0.25, b,e - s = 0.5, cf- s = 0.75) Fig. 9. Streamline distributions for various s (Re = (a,b,c) 1, and (d,ef 10; s = (a,d) 0.25, (b,e) 0.5, and (cf) 0.75)
Заключение
В результате проделанной работы выполнено математическое моделирование течения ньютоновской жидкости по каналу с конструктивным элементом типа затвора. В результате параметрических исследований получены картины течения в зависимости от числа Рейнольдса в диапазоне изменения от 0.1 до 80, радиуса затвора - от 1 до 2, степени закрытия запорного элемента - от 0 до 0.5. Показано разделение потока на три характерные зоны: зону двумерного течения в окрестности затвора и зон одномерного течения вне его. Продемонстрировано влияние значений определяющих параметров на формирование и размеры выделяемых зон в потоке.
ЛИТЕРАТУРА
1. Fester V., Slatter P., Alderman N. Resistance coefficients for non-Newtonian flows in pipe fittings // Rheology. InTech. 2012. P. 151-186.
2. Sisavath S., JingX., Pain C.C., Zimmerman R.W. Creeping flow through axisymmetric sudden contraction or expansion // Journal of Fluids Engineering, Transactions of the ASME. 2002. V. 124. Iss. 1. P. 273-278. DOI: 10.1115/1.1430669.
3. Pienaar V.G. Non-Newtonian fitting losses. Unpublished M Tech. thesis. Cape Technikon. Cape Town, 1998.
4. Fester V.G., Kazadi D.M, Mbiya B.M., Slatter P.T. Loss coefficients for flow of Newtonian and non-Newtonian fluids through diaphragm valves // Chemical Engineering Research and Design. 2007. V. 85. Iss. 9A. P. 1314-1324. DOI: 10.1205/cherd06055.
5. Pienaar V.G. Viscous flow through sudden contractions. Dissertation submitted in fulfilment of the degree Doctor technologiae in the Faculty of Engineering. Cape Technikon, 2004.
6. Andhale V.A., Deshmukh D.S. Investigation of ball valve design for performance enhancement // Pratibha: International Journal of Science, Spirituality, Business and Technology. 2016. V. 4. Iss. 2. P. 105-112.
7. Zhang S.C., Zhang Y.L., Fang Z.M. Numerical simulation and analysis of ball valve three-dimensional flow based on CFD // IOP Conf. Series: Earth and EnvironmentalScience. 2012. V. 15. Iss. 5. P. 1-6. DOI: 10.1088/1755-1315/15/5/052024.
8. Годунов С.К., Рябенький В.С. Разностные схемы. М.: Наука, 1977. 440 с.
9. Яненко Н.Н. Метод дробных шагов решения многомерных задач математической физики. Новосибирск: Наука, 1967. 197 с.
10. Roache P.V. Fundamental of computational fluid dynamics. New Mexico: Hermosa Publishers, 1998. 648p.
Статья получена 17.11.2016 г.
Banzula Yu.B., Borzenko E.I., Karyazov, S.V., Shrager G.R. THE KINEMATICS OF A VISCOUS FLUID FLOW IN A CHANNEL WITH A VALVE. Tomsk State University Journal of Mathematics and Mechanics. 6(44). pp. 54-63
DOI 10.17223/19988621/44/5
In this paper, results of mathematical modeling of a steady flow of a Newtonian fluid in a flat channel are presented. At some distance from the inlet and outlet of the channel, a valve having a specified diameter and capable to change the width of the flow cross section depending on the position of the locking element is situated. On the walls, the adhesion conditions are performed; at the input, the profile corresponding to the steady flow of a predetermined constant flow is assigned, and the Neumann boundary conditions are used at the output. The mathematical statement of a problem is formulated using a dimensionless stream function and vorticity variables. The numerical solution of the problem is implemented with the application of an original finite-difference method based on the scheme of alternating directions. The physical flow area with a curved boundary is converted into a rectangular domain by introducing new coordinates, and the stationary solution of the problem is obtained using the relaxation method. The calculations demonstrate that in the vicinity of the inlet and outlet boundaries, the one-dimensional flow areas corresponding to a steady fluid flow in an infinite channel occur, and in the vicinity of the shutter, a two-dimensional flow with recirculations is observed. The effect of the Reynolds number on the distribution of the kinematic characteristics is investigated. The results of the parametric studies of the flow pattern in relation to the opening position and diameter of the locking element are presented.
Keywords: flow, viscous fluid, valve, numerical simulation, kinematics of the flow.
BANZULA Yuriy Borisovich (Doctor of Technical Sciences, Federal State Enterprise the Federal Center for Dual-Use Technologies Soyuz, Dzerzhinskiy, Moscow region, Russian Federation)
BORZENKO Evgeniy Ivanovich (Candidate of Physics and Mathematics, Tomsk State University, Tomsk, Russian Federation) E-mail: [email protected]
KARYAZOVSvyatoslav Vladimirovich (Candidate of Technical Sciences, Federal State Enterprise the Federal Center for Dual-Use Technologies Soyuz, Dzerzhinskiy, Moscow region, Russian Federation)
SHRAGER Gennady Rafailovich (Doctor of Physics and Mathematics, Professor, Tomsk State University, Tomsk, Russian Federation) E-mail: [email protected]
REFERENCES
1. Fester V., Slatter P., Alderman N. (2012) Resistance coefficients for non-Newtonian flows in pipe fittings. Rheology. InTech. pp. 151-186.
2. Sisavath S., Jing X., Pain C.C., Zimmerman R.W. (2002) Creeping flow through axisymmet-ric sudden contraction or expansion. Journal of Fluids Engineering, Transactions of the ASME. 124(1). pp. 273-278. DOI 10.1115/1.1430669.
3. Pienaar V.G. (1998) Non-Newtonian fitting losses. Unpublished M Tech. thesis. Cape Technikon. Cape Town.
4. Fester V.G., Kazadi D.M, Mbiya B.M., Slatter P.T. (2007) Loss coefficients for flow of Newtonian and non-Newtonian fluids through diaphragm valves. Chemical Engineering Research and Design. 85(9A). pp. 1314-1324. DOI 10.1205/cherd06055.
5. Pienaar V.G. (2004) Viscous flow through sudden contractions. Dissertation. Cape Technikon.
6. Andhale V.A., Deshmukh D.S. (2016) Investigation of ball valve design for performance enhancement. Pratibha: International Journal of Science, Spirituality, Business and Technology. 4(2). pp. 105-112.
7. Zhang S.C., Zhang Y.L., Fang Z.M. (2012) Numerical simulation and analysis of ball valve three-dimensional flow based on CFD. IOP Conf. Series: Earth and Environmental Science. 15(5). pp. 1-6. DOI 10.1088/1755-1315/15/5/052024.
8. Godunov S.K., Ryabenkii V.S. (1987) Difference schemes. North-Holland: Elsevier Science Ltd.
9. Yanenko N.N. (1971) The method of fractional steps: the solution of problems of mathematical physics in several variables. Springer-Verlag.
10. Roache P.V. (1998) Fundamental of computational fluid dynamics. New Mexico: Hermosa Publishers.
