CC BY
18
УДК 539.374+539.377
DOI: 10.14529/mmph170305
О ЗАДАЧЕ ИДЕНТИФИКАЦИИ НАЧАЛЬНОГО ГИДРОДИНАМИЧЕСКОГО УЧАСТКА ПРИ ЛАМИНАРНОМ ТЕЧЕНИИ НЬЮТОНОВСКОЙ ЖИДКОСТИ В ГОРИЗОНТАЛЬНОМ КОЛЬЦЕВОМ КАНАЛЕ
А.В. Ряжских
Воронежский государственный технический университет, г. Воронеж, Российская Федерация E-mail: ryazhskihav@bk.ru
В рамках физической линеаризации уравнений Навье-Стокса об однонаправленном аксиальном ламинарном течении ньютоновской жидкости найдено аналитическое решение задачи идентификации поля скоростей во входной области кольцевого канала, что позволило получить в явном виде соотношение для расчета длины начального гидродинамического участка.
Ключевые слова: кольцевой канал; ньютоновская жидкость; длина начального гидродинамического участка.
Кольцевая геометрия, наряду с цилиндрическими трубами и плоскими каналами, наиболее часто используется в различных технических системах [1]. В [2] показано, что проведение процесса теплопередачи в режиме начального гидродинамического участка является одним из способов интенсификации теплообмена. Поэтому выявление закономерностей формирования структуры поля скоростей на начальном участке и определение его длины в направлении течения жидкости актуально.
Впервые такая задача для кольцевого канала была рассмотрена в [3] на основе погранслойно-го приближения с применением линеаризации Лангхаара [4]. Детальный обзор последующих теоретических и экспериментальных исследований этой проблемы [5] показал, что различие в определении характеристик начального гидродинамического участка существенно. Уточнение длины начального гидродинамического участка на основе численного интегрирования уравнений Навье-Стокса, например [6], вызывает вопросы корректности используемых вычислительных процедур, связанных с существованием и единственностью решения самих уравнений. Предпринимаются также дальнейшие попытки отыскания аналитических решений, например, использующих идеализацию о «ползущем» течении [7].
В связи с этим в данной работе с позиций физической линеаризации уравнения пограничного слоя на основе гипотезы однонаправленности течения формулируется и аналитически анализируется задача отыскания длины начального гидродинамического участка в горизонтальном кольцевом канале.
Постановка задачи. Математическая формулировка задачи в цилиндрической системе координат имеет вид [8]:
dvr dvr dvrJ dp
dvr dvr
> 1 - + vz r
r dr z dz
dvz
- + v, -- + vz
dr
2 d i r v 1 - +1 (v +dv
r dr V dr J r2 dz V dr dz
(1)
I dvz +v dvz +v dvz j _ dp + I 1 d
dt r dr z dz J dz I r dr
dvz + dvr
dr dz
+ 2^1= (2)
1 d,, dvz
rvr )+-oz=0; (3)
r or dz
vr (r, z,0) = vz (r, z,0) = 0, vr (r,0, t ) = 0, vz (r,0, t ) = v0 = const; (4)
dvr (r,l,t) dv_ (r,l,t) , ч , \ r, i \ i \ г, /^ч
rV ' = z V ' 7 = 0, vr (ri, z, t ) = vz (ri, z, t ) = 0, vr (r2, z, t ) = vz (r2, z, t ) = 0, (5) dz dz
где vr, vz - радиальная и аксиальная компоненты вектора скорости жидкости; p - избыточное давление; p,h - плотность и коэффициент динамической вязкости жидкости; t - время; r1, r2 -радиусы внутренней и внешней соосных цилиндрических труб, составляющих кольцевой канал; l - длина кольцевого канала (l » r2).
r
Как правило, для замыкания системы (1)-(5) дополнительно рассматривается уравнение относительно давления [9]
1d fr dp 1+ о
r dr V dr J dz2
i d f rv du
. r dr V r dr
1 d
+--
r dr
d( rvr)
dz
0 f1 d2 f du
+2'i Vr
dr
++dv _ 1 d fv 1+1Л
dz3 r dr V r J r drdz
+_d (v du
dz V r dr
dvz dvr
r I —^ + ■ r
d ( dvz
+ — I vz—z
dz V dz
+
dr dz
с краевыми условиями
p (r,0, t ) = p0 = const, dp (rt) = 0,
dz
dp (r1, z, t) dp (r2, z, t)
= 0.
(6)
(7)
(8)
Эг Эг
Пусть течение по всей длине канала является однонаправленным, т. е. ьг ° 0, тогда из (3) следует, что ь2 =ь2 (г, t), причем (1) вырождается в Эр/Эг = 0, а (2) трансформируется в
dv_ dp h d f dv_ p—- = —- + ——I r—z dt dz r dr V dr
с граничными условиями
v (r,0 ) = vo, v (ri, t ) = vz (r2, t ) = 0 и выполнением балансового соотношения
2
2 2 r _ r
'2
J rvz (r, t) dr = v .
(9) (10) (11)
При этом из (6) следует постоянство аксиального градиента давления, т. е. Эр/dz = const. В безразмерном виде система (9)-(11) запишется
dV=с+±1A f r dV
dZ Re R dR { dR V (R,0) = 1; V (R1, Z ) = V (R2, Z ) = 0;
R2 _ R1 R
2
J RV (R, Z)dR = 1,
(12)
(13)
(14)
где Т = = = Z; R = ; V = ь2 /ь0 ; Р = р/(рЬ); R1,2 = г12/й ; Re = ь0йр/Т] - число Рейнольдса; й = 2(г2 — г1) - гидравлический диаметр; С = — dP|dZ .
Решение. Применим к (12), (13) одностороннее интегральное преобразование Лапласа по переменной Z [10]:
d 2V 1 dV,
C
-- +--: — sReVт + I — +1 I Re = 0;
dR2 R dR : К 5 )
V:. (Rl, 5 ) = VL (R2, 5 ) = 0,
где s,VL - изображение Zи Vсоответственно. Решение (15), (16) есть
VL (R,5) = (С + s)/s2 + р(R,5)/у(5),
где р(R,5) = (С + 5){[К0 (R1^¡sR£) — К0 (R24sR£)] 10 (R^/R) +
+[ 10 (R24R)—10 (R14R)] к0 (R^/R)};
у(5) = 52у(5), у(5) = 10 (R1,¡sR£)К0 (R24sR£) —10 (R24sR£)К0 (R1^fsR£) 10, К0 - функции Бесселя и Макдональда мнимого аргумента [11]. Оригинал первого слагаемого правой части (17) равен [12]
(15)
(16)
(17)
2
Г1 [(С + 5 )/5 2 ] = С2 +1, (18)
где Г-1 - обратный оператор преобразования Лапласа.
Для нахождения оригинала второго слагаемого правой части (17) применим вторую теорему Ващенко-Захарченко [13]. Корни знаменателя у/( 5) = 0 таковы: 1) 5 = 0 - кратности два; 2) бесконечное счетное множество однократных корней 5п уравнения /(5п) = 0 . Для корня 5 = 0
р(Я,0)
.. . = lim
у/( 0)
Z exp (sZ) ((R,s) + exp (sz ) - exp (sZ ()Y( S)
Y(s) Y(s) у (s)
Так как I'0 (z) = I, (z) и K'0( z) = -K, (z), то (p's (R, s) = [ K0 (yfsRe) - K0 (R2JsRe )] I0 (RyfsRe) + [ I0 (JsRe) -10 (JsRe )] K0 ((/¡Re) +(C + s)4Re {[R2K1 (R24sRe) -R1K1 (Rл/sRe )]I0 (RyfsRe) + +[K0 (yfsRe) -K0 (R2yfsRe )]RI1 (RyfsRe) + [R211 ((yfsRe) - RJI, (( yfsRe )]K0 (RyfsRe) + +[I0 (yfsRe)-10 (R2*JsR£)]RK1 (RyfsRe )}/(2>/J), Y (s ) = JRe [ RlIl (VsRe )K0 (-JsRe)- R210 (JsRe )) (JsRe )--R2I, (r2 VsRe )K0 ((VsRe) + RjI0 (yfsRe )) ((VsRe )]/(2>/s),
+
тогда
= -CZ + J CRe [-(R2 -R2 ) + (R2-R2 )ln (R/R )ln (R/ R)] .
Для корней sn имеем
sn ) ( Z)
^-^TV exP (snZ),
(19)
(20)
п=1 ¥ (5п)
откуда следует, что физическому смыслу задачи удовлетворяют чисто мнимые корни, поэтому
45пКе = Мп1 или 5п = —ЦКе .
п п
Так как 10 (1г) = J0 (г); К0 (1г) = 1—^ (г); 11 (1г) = и1 (г); К1 (1г)=—Ы1 (г), где J0 (г),
J1 (г) - функции Бесселя первого рода; Ы0 (г), Ы1 (г) - функции Неймана [14], то
(
Г цП > ( R,-ц
v Re / v
2 Л
C-
¥
ц
Re
( M ^
iП{[N (RiMn)-N0 (R2M )]Л (RMn)+[Л (RMn)-Л (RiMn )]N0 (RMn )};(21)
ц\ П
R i 7 [-RJl (RiMn )N (R2Mn ) - R2 Л (RiMn )Ni (R2Mn ) Re 4 L
v Re /
+R2 J (RM )N (RiMn)+Ri Л (RM )Ni (RiMn )]. Суммируя (18)-(20), с учетом (21) и (22), получим
+
(22)
V (R, Z )=1 CRe 4
-(r 2 - R2 иМ-^ inR
in R2 Ri
Ri
-2^(CRe-цЦ){[N0 (RiMn) -N0 (R2Mn)]J (RMn)
+
n=1
+
Mn
[ J (R2M) - J0 (RiMn N0 (RMn)} exp ^ - R^Z
/ {mM [-Ri Ji (RiMn )N0 (R2Mn) - R2 J0 (RiMn )Nj (R2M)
+R2 Jl (R2Mn )N0 (RiMn ) + Ri J0 (R2Mn )Ni (RiMn , где Mn - корни уравнения
+
(23)
Ряжских А. В.
R1 =1R /(1 — R), R2 = I /(1 — R), R = г1/ г2 ;
R = 1 г2 ; константа С определена из (14) при Z ® <
С = 32 (1 — к )21[1 + к + (1 — к 2)/ 1п я ] ,
(24)
совпадающая с результатом в [15].
Анализ. Характерная структура гидродинамического поля безразмерных скоростей на входном участке (рис. 1) показывает, что при малых значениях R наблюдается ассиметрия профиля скорости со сдвинутым максимальным значением в сторону поверхности соосного цилиндра
0.1 I I I I I ' 0.6 0.7 0.8 0 9
Рис. 2. Вид функции f (Е )
Рис. 1. Поля скоростей, примыкающих ко входному сечению кольцевого канала при различных числах Re и отношениях радиусов внутренней и внешней соосных цилиндрических труб Ё: а) Не = 100,Ё = 0,99 ; б) Не = 100,Ё = 0,01 ; в) В.е = 10, Л = 0,5 ; В.е = 1000, Л = 0,5 ; г) В.е = 1000, Л = 0,5
меньшего радиуса, а увеличение числа Re практически линейно увеличивает длину гидродинамического начального участка.
Предполагая, что дрейфом радиальной координаты максимальной скорости на начальном гидродинамическом участке можно пренебречь и считать ее постоянной и равной
1
Рис. 3. Изменение максимальной безразмерной скорости вдоль кольцевого канала при Ре = 1: □ - расчет; о - эксперимент
R* =-
R +1
'1пR , найденной из условия 1 — V(R*,Z)/V(R*,^) = е,
242\\R—1
где относительная точность е принимается, как правило, равной 0,02 [16], и ограничиваясь регулярным режимом (в решении (23) учитывается только первое слагаемое ряда), найдено выражение длины начального гидродинамического участка в кольцевом канале:
ZL = f (R) Re
(25)
где
f(й) = 1п{<(<)/е(R\«)]}; р(<) = 2(с& —т){[N (Ж) — N0 (R2ml)]Л () + [Л (R2Ml) — Л (АМ)]N0 (R^м )}/ {м3 [—RЛ (RlMl) N0 ()—R2 Л (Ам) N1 (R2Ml)+R2 ¿1 (R2Ml) N0 (ж)+^ ¿0 (R2Ml) N1 ()]}
Как следует из (25), длина начального участка пропорциональна числу Re (константа C обратно пропорциональна Re ), что согласуется с оценками [17]. Отметим, что при R ® 0 и R ® ¥ получаем частные случаи для круглой трубы и плоского канала соответственно (рис. 2). Количественное сравнение полученных результатов проведено с классическими экспериментальными данными [3] (рис. 3). Различие рассчитанных и экспериментальных значений скорости в области, примыкающей ко входному сечению, может быть объяснено тем, что в модели не учитывается расход кинетической энергии потока жидкости на нивелирование неоднородности давления по сечению канала. Тем не менее, это не оказывает существенного влияния на длину гидродинамического участка.
Заключение. Синтезированная модель гидродинамики во входной области кольцевого канала характеризуется качественной и количественной адекватностью и может быть применена для расчета длины начального гидродинамического участка при ламинарном течении вязких несжимаемых жидкостей.
Литература
1. Massoud, M. Engineering termofluids: thermodynamics, fluid mechanics and heat transfer / M. Massoud. - Springer-Verlag Berlin Heidelberg, 2005. - 1120 p.
2. Jiji, Latif M. Heat convection / Latif M. Jiji. - Springer-Verlag Berlin Heidelberg, 2009. - 543 p.
3. Heaton, H.S. Heat transfer in annular passages-simultaneous development of velocity and temperature fields in laminar flow / H.S. Heaton, W.C. Reynolds, W.M. Kays // Int. J. Heat and Mass Transfer. - 1964. - Vol. 7, № 7. - P. 763-781.
4. Кэйс, В.М. Конвективный тепло- и массообмен / В.М. Кэйс. - М.: Энергия, 1972. - 448 с.
5. The development lengths of laminar pipe and channel flows / F. Durst, S. Ray, B. Unsal, OA. Bayoumi // Journal of Fluids Engineering. - 2005. - Vol. 127, no. 6. - P. 1154-1160.
6. Shaker, A.A. A numerical study of low Reynolds number incompressible flow of entrance and disturbed regions of concentric circular pipes / A.A. Shaker // J. of Engineering and Development. -2012. - Vol. 16, no. 2. - P. 16-33.
7. Ryazhskikh, A.V. Hydrodynamic entrance length for high-viscosity Newtonian fluid flow in an annular channel / A.V. Ryazhskikh, S.V. Ryabov // Journal of Engineering Physics and Thermophysics.
- 2013. - Vol. 86. - Issue 2. - P. 396-401.
8. Poole, R.J. Development length requirements for fully developed laminar flow in concentric annuli / R.J. Poole // Journal of Fluids Engineering. - 2010. - Vol. 132, № 6. - P. 501-504.
9. Берд, Р. Явления переноса / Р. Берд, В. Стьюарт, Е. Лайтфут. - М.: Химия, 1974. - 688 с.
10. Дёч, Г. Руководство к практическому применению преобразования Лапласа и z-преобразования / Г. Дёч. - М.: Физматгиз, 1971. - 288 с.
11. Янке, Е. Специальные функции / Е. Янке, Ф. Эйде, Ф. Лещ. - М.: Наука, 1968. - 344 с.
12. Диткин, В.А. Справочник по операционному исчислению / В.А. Диткин, А.П. Прудников.
- М.: Высшая школа, 1965. - 465 с.
13. Беляев, Н.М. Методы теории теплопроводности. В 2 ч. Ч.1. / Н.М. Беляев, А. А. Рядно -М.: Высшая школа, 1982. - 327 с.
14. Korenev, B.G. Bessel functions and their applications / B.G. Korenev. - NY: Taylor & Francis, 2002. - 288 p.
15. Идельчик, И.Е. Справочник по гидравлическим сопротивлениям / И.Е. Идельчик. - М.: Машиностроение, 1992. - 672 с.
16. Слезкин, Н.А. Динамика вязкой несжимаемой жидкости / Н.А. Слезкин. - М.: ГИТТЛ, 1955. - 519 с.
17. McComas, S.T. Hydrodynamic entrance length for ducts of arbitrary cross section / S.T. McComas // J. Basic Eng. - 1967. - Vol. 89, no. 4. - P. 847-850.
Поступила в редакцию 26 февраля 2017 г.
Ряжских А. В.
Bulletin of the South Ural State University Series "Mathematics. Mechanics. Physics" _2017, vol. 9, no. 3, pp. 34-40
DOI: 10.14529/mmph170305
ON THE IDENTIFICATION OF ENTRANCE HYDRODYNAMIC REGION IN CASE OF LAMINAR FLOW OF NEWTONIAN FLUID IN HORIZONTAL ANNULAR CHANNEL
A.V. Ryazhskikh
Voronezh State Technical University, Voronezh, Russian Federation E-mail: ryazhskihav@bk.ru
In the frameworks of physical linearization of the Navier-Stokes equations in a cylindrical coordinates system on the one-way axial force-feed laminar flow of Newtonian fluid, a mathematical model of the flow development in the entrance region of a horizontal annular channel is formulated. The unknown constant gradient of pressure along the channel is connected with the equation of continuity written in an integral form of stability of liquid flow in any cross section of a channel. Use of the one-way integral Laplace transformation along the longitudinal coordinate allowed to obtain an analytical expression of the local hydrodynamic field at the entrance region and determine pressure losses coincided with the classic data. Analysis of the characteristic structure of the hydrodynamic field of dimensionless velocities at the entrance region showed that for small values of the ratio of the radii of the inner and outer coaxial cylindrical tubes constituting the annular channel, asymmetry of the longitudinal velocity profile is observed with a shift of the maximum value towards the surface of a coaxial cylinder of smaller radius, and an increase in the Reynolds number practically linearly increases the length of the hydrodynamic entrance region. Assumption about the absence of drift of the radial coordinate of the maximum velocity in the entrance hydrodynamic region, limited to the so-called "regular" regime, made it possible to identify the length of the entrance hydrodynamic region in the annular channel by the completed expression in an explicit form that correlates with the classical estimates obtained as a result of computational experiments. It is noted that when the ratio of the radii of the inner and outer coaxial cylinders approaches zero or infinity (corresponding to particular cases of a circular tube and a flat channel), the known results for the lengths of the entrance hydrodynamic regions are obtained. Difference between velocity values calculated by the proposed model and experimental values in the region adjacent to the entry section is explained by the fact that kinetic energy of the liquid flow is not accounted for by leveling the pressure inhomogeneity along the channel cross-section. Nonetheless, it is shown that it does not have a significant influence on the length of hydrodynamic entrance region.
Keywords: annular channel; Newtonian fluid; length of the entrance hydrodynamic region.
References
1. Massoud M. Engineering termofluids: thermodynamics, fluid mechanics and heat transfer. Springer-Verlag Berlin Heidelberg, 2005, 1120 p. DOI: 10.1007/b138870
2. Jiji Latif M. Heat convection. Springer-Verlag Berlin Heidelberg, 2009, 543 p. DOI: 10.1007/978-3 -642-02971 -4
3. Heaton H.S., Reynolds W.C., Kays W.M. Heat transfer in annular passages-simultaneous development of velocity and temperature fields in laminar flow. Int. J. Heat and Mass Transfer, 1964, Vol. 7, no. 7, pp. 763-781. DOI: 10.1016/0017-9310(64)90006-7
4. ^ys WM. Convective heat and mass transfer. McGraw-Hill Science, 1993, 480 p.
5. Durst F., Ray S., Ünsal B., Bayoumi O.A. The development lengths of laminar pipe and channel flows. J. Fluids Eng., 2005, Vol. 127, no. 6, pp. 1154-1160. DOI: 10.1115/1.2063088
6. Shaker A.A. A numerical study of low Reynolds number incompressible flow of entrance and disturbed regions of concentric circular pipes. J. of Engineering and Development, 2012, Vol. 16, no. 2, pp. 16-33. (http://www.iasj.net/iasj?func=fulltext&aId=67504)
7. Ryazhskikh A.V., Ryabov S.V. Hydrodynamic entrance length for high-viscosity Newtonian fluid flow in an annular channel. Journal of Engineering Physics and Thermophysics, 2013, Vol. 86, Issue 2, pp. 396-401. DOI: 10.1007/s10891-013-0847-6
8. Poole R.J. Development length requirements for fully developed laminar flow in concentric annuli. J. Fluids Eng., 2010, Vol. 132, no. 6, pp. 501-504. DOI: 10.1115/1.4001694
9. Bird R.B., Stewart W.E., Lightfoot Е.N. Transport phenomena. New York, John Wiley, 2007, 928 p.
10. Doetsch, G. Anleitung zum praktischen gebrauch der Laplace-transformation und der Z-transformation. München-Wien: R. Oldenbourg Verlag, 1967, 256 p.
11. Janke Е., Emde F., Lösch F. Tafeln höherer Funktionen. Stuttgart, B.G.Teubner Verlagsgesellschaft, 1960, 318 p.
12. Ditkin V.A., Prudnikov A.P. Spravochnik po operatsionnomu ischisleniyu (Reference book on operational calculus). Moscow, Vysshaya shkola, 1965, 465 p. (in Russ.).
13. Belyaev N.M., Ryadno A.A. Metody teorii teploprovodnosti. V 2 ch. Ch.1. (Heat conductivity theory techniques. In two parts. Part 2). Moscow, Vysshaya shkola, 1982, 327 p. (in Russ.).
14. Korenev B.G. Bessel functions and their applications. NY: Taylor & Francis, 2002, 288 p.
15. Idel'chik I.E. Spravochnik po gidravlicheskim soprotivleniyam (Reference book on hydraulic resistance). Moscow, Mashinostroenie Publ., 1992, 672 p. (in Russ.).
16. Slezkin N.A. Dinamika vyazkoy neszhimaemoy zhidkosti (Dynamics of viscous incompressible liquid). Moscow, GITTL publ., 1955, 519 p. (in Russ.).
17. McComas S.T. Hydrodynamic entrance length for ducts of arbitrary cross section. J. Basic Eng., 1967, Vol. 89, no. 4, pp.847-850. DOI: 10.1115/1.3609713
Received February 26, 2017
