close

Вход

Забыли?

вход по аккаунту

- Вестник МГСУ

код для вставкиСкачать
Гидравлика. Инженерная гидрология. Гидротехническое строительство
ГИДРАВЛИКА. ИНЖЕНЕРНАЯ ГИДРОЛОГИЯ.
ГИДРОТЕХНИЧЕСКОЕ СТРОИТЕЛЬСТВО
УДК 532.542.1
А.Л. Зуйков
ФГБОУ ВПО «МГСУ»
АКСИАЛЬНЫЕ И РАДИАЛЬНЫЕ СКОРОСТИ
В ПОЛЗУЩЕМ ТЕЧЕНИИ В ТРУБЕ
Приведено аналитическое исследование полей аксиальных и радиальных скоростей в установившемся неравномерном ползущем течении на начальном участке цилиндрического канала. Показано, что поля скоростей такого течения являются
двумерными и определяются функцией тока. Распределения скоростей получены в
виде рядов Фурье — Бесселя.
Ключевые слова: ньютоновская жидкость, неравномерное ползущее течение, азимутальный вихрь, функция тока, аксиальные и радиальные скорости, разложение Фурье — Бесселя.
Настоящая статья является продолжением серии работ [1—3], где были
получены нормированные аналитические функции радиально-аксиальных
распределений таких характеристик установившегося неравномерного ползущего течения в цилиндрической трубе радиусом R, как азимутальный вихрь
∞
∂u&r ∂u& x
(1)
−
= 4r& + ∑ An J1 (λ n r&) exp(−λ n x& ), ∂x&
∂r&
n =1
и функции тока
2 
∞

& (r&, x& ) = r& 2 1 − r&  − r& ∑ An  J 0 (λ n ) J (λ r&) − r&J (λ r&) exp(−λ x& ), (2)
Ψ


1
n
0
n
n

2  2 n = 1 λ n  J1 (λ n )


связанные соотношением
& ∂2 Ψ
& ∂Ψ
&
∂2 Ψ
&θ, +
−
= − r&ω
(3)
r&∂r&
∂x& 2
∂r& 2
здесь r и х& — нормированные радиальная и продольная координаты (см. рис. 1
[3]), r = r R , х& = х R; u&r и u& x — нормированные радиальные и аксиальные
локальные скорости течения, u&r = u x V , u& x = u x V ; V — среднерасходная скорость течения, V = Q/pR2; Q — объемный расход; J 0 (λ n r&), J1 (λ n r&) — функции
Бесселя первого рода нулевого и первого порядков [4]; ln — константа разделения; An — постоянная, зависящая от граничных условий.
Из равенства (2) можно заключить, что J1 (λ n ) ≠ 0 , ибо стоит в знаменателе первого слагаемого в квадратных скобках, т.е. ln не является действительным нулем функции Бесселя первого рода первого порядка. Следовательно,
единственным приемлемым решением в данном случае, удовлетворяющим
граничным условиям (6) и (8) из [3], может быть только равенство нулю функции J 0 (λ n ) = 0, где ln — действительный нуль функции Бесселя первого рода
& θ (r&, x& ) =
ω
© Зуйков А.Л., 2014
127
5/2014
нулевого порядка. Отсюда можно вычислить неизвестную ранее постоянную
Ап, приведя (2) к виду
2 
2 ∞

& (r&, x& ) = r& 2 1 − r&  + r& ∑ An J (λ r&) exp(−λ x& ).
(4)
Ψ
0
n
n

2  2 n = 1 λ n

При этом на оси и стенках по всей длине трубы сохраняются корректные
 (0, x ) = 0 при r = 0 и Ψ
 (1, x ) = 0,5 при r = 1, а на бескоусловия (8) из [3]: Ψ
нечном удалении от входа в трубу (при x = ∞ ) имеем профиль Пуазейля, соответствующий условию (6) [3]

∂Ψ
∂   r 2  
u x = =  r 2 1 −   =2 (1 − r 2 ) . (5)
2 
r∂r r∂r  
Рассмотрим трансформацию профилей аксиальных и радиальных скоростей неравномерного ползущего течения по длине цилиндрической трубы. Для
этого, прежде всего, определим граничные условия задачи.
Кроме профиля Пуазейля (5) при x = ∞ и условия сохранения объемного
расхода по длине трубы
1
∫ u& x 2r&dr& = 1 для
x ≥ 0, (6)
0
принятыми в работе [3], дополнительными граничными условиями положим:
равенство нулю радиальных скоростей при x = ∞ соответствующем равномерному движению жидкости по Пуазейлю;
равномерное распределение аксиальных скоростей по сечению на входе в
трубу при x = 0 , при котором u& x = 1;
ввиду осевой симметрии течения равенство нулю на оси трубы (при r = 0 )
радиальных скоростей и частной производной от аксиальной скорости по радиусу, отвечающее условию мягкого экстремума функции u& x ;
равенство нулю скоростей движения жидкости на стенках трубы (при r = 1 )
ввиду вязкости жидкости и непроницаемости стенок.
Дополнительные граничные условия описываются системой
u&r = 0
при x& = ∞ для 0 ≤ r& ≤ 1; 
при x& = 0 для 0 ≤ r& ≤ 1; 
u& x = 1

(7)
∂u& x

&
&
&
ur = 0,
= 0 при r = 0 для x ≥ 0; 
∂r&

u&r = 0, u& x = 0
при r& = 1 для x& > 0. 
Переписав функцию (4) на входе в трубу (при x = 0 ) в виде
2 
2 ∞

& (r&,0) = r& 2 1 − r&  + r& ∑ An J (λ r&),
Ψ
0
n

2  2 n = 1 λ n

где, следуя второму граничному условию системы (7)
&
∂Ψ
u& x = 1 =
,
r&∂r&

Ψ
r
2
 = rdr
 = r . 

после интегрирования ∫ d Ψ
,
получим
Ψ
∫0
2
0
128
(8)
(9)
ISSN 1997-0935. Vestnik MGSU. 2014. № 5
Гидравлика. Инженерная гидрология. Гидротехническое строительство
Тогда выражение (8) после очевидных сокращений перепишется как
r&3 − r& =
∞
A
∑ λn
n =1
r&J 0 (λ n r&). (10)
n
Воспользуемся условиями ортогональности функций Бесселя [4], согласно
которым, если li и lj — два нуля функции J m (λ n ) = 0 — действительные, то

если i ≠ j; 

0
(11)

1
1
2

& & = [ J m′ (λ i ) ] , если i = j .
∫0 J m (λi r&)J m (λ j r&)rdr

2

Тогда, умножая правую и левую части (10) последовательно на J 0 (l1r)dr,
J 0 (λ 2 r&)dr& , …, J 0 (λ п r&)dr& и далее до J 0 (l ∞ r)dr, и интегрируя по интервалу от
0 до 1, согласно условиям ортогональности получим систему равенств
1
∫J
m
& & = 0,
(λ i r&)J m (λ j r&)rdr


0
0

1
1
1

A
3
2











r
J
r
dr
rJ
r
dr
J
r
J
r
rdr
l
−
l
=
l
l
(
)
(
)
(
)
(
)
;
0
2
0
2
∫0 0 2
∫0 0 2
l 2 ∫0




(12)

1
1
1
An

3
 ;
∫0 r J 0 (l n r)dr − ∫0 rJ 0 (l n r)dr= l n ∫0 J 0 (l n r) J 0 (l n r)rdr





1
1
1
A∞
3
 .
∫0 r J 0 (l∞ r)dr − ∫0 rJ 0 (l∞ r)dr= l ∞ ∫0 J 0 (l∞ r) J 0 (l∞ r)rdr


Отсюда для произвольного п-го частного решения в результате интегрирования получим
An
2
1
1
2
(13)
J 2 (l n ) −
J 3 (l n ) −
J1=
(l n )
[ J 0′ (l n )] , 2
ln
ln
ln
2l n
где J2(ln), J3(ln) — функции Бесселя первого рода второго и третьего порядков.
Далее в результате использования рекуррентных соотношений [4], связывающих цилиндрические, в т.ч. функции Бесселя, из равенства (13) находим
8
An = − 2
.
(14)
λ п J1 (λ n )
Обратимся теперь к уравнениям (1) и (3), приведя их согласно (14) к виду
∞
& ∂2 Ψ
& ∂Ψ
&
J (λ r&)
∂2 Ψ
& 2 + 8r& ∑ 1 n exp(−λ n x& ), 4
r
−
=
−
+
(15)
2
2
2
r&∂r&
∂r&
∂x&
n = 1 λ п J1 (λ n )
и замечая далее, что
& ∂Ψ
&
& 
∂u&
∂  ∂Ψ
∂2 Ψ
 = r& x ,
−
= r& 
2
r&∂r&
∂r&
∂r&  r&∂r& 
∂r&
1
1
3
∫ r J 0 (l1r)dr − ∫ rJ 0 (l1r)dr=
1
A1
 ;
J 0 (l1r) J 0 (l1r)rdr
l1 ∫0
Hydraulics. Engineering hydrology. Hydraulic engineering
129
5/2014
а по (4) и (14)
&
∂2 Ψ
∞
J (λ n r&)
exp(−λ n x& ) ,
(
λ
)
∂х
n
1
n
n =1
с учетом записанного находим
2
∑ λ 0J
= −4r& 2
∞ 
∂u& x
2 J (λ r&) r&J (λ r&) 
= −4r& + 4 ∑  2 1 n + 0 n  exp(−λ n x& ).
∂r&
 λ п J1 (λ n ) λ п J1 (λ n ) 
n = 1
(16)
Разделяя переменные и интегрируя в пределах по радиусу от r до 1 и по
аксиальной скорости от u& x до 0 на стенке трубы, получим
∞
4  r&J (λ r&) 2 J (λ r&) 
(17)
u& x (r&, x& ) = 2(1 − r& 2 ) − ∑ 2 1 − 1 n + 0 n  exp(−λ n x& ).
J1 (λ n ) λ n J1 (λ n ) 
n =1 λп 
Заметим, что в соответствии с распределениями (16) и (17) аксиальные
скорости имеют мягкий экстремум на оси трубы при r = 0, ибо в этом случае производная ∂u x ∂r =0, и одновременно u& x = 0 при r = 1. Таким образом,
граничные условия два и три системы (7) выполнены.
Проинтегрируем теперь распределение (17) для проверки условия (6)
1
∞


 &
&
& 
2
& &.
&
&
&
&
=
−
r
−
u
2
rdr
2
1
) ∑ λ42 1 − rJJ1 ((λλn r)) + 2λJJ0 (λ(λn r))  exp(−λ n x&)  2rdr
∫0 x
∫0  (
n =1 п 
1
n
n 1
n 


Можно видеть, что условие сохранения объемного расхода по длине трубы не
выполняется, ибо имеет место равенство
1
∞
1
∫ u& x 2r&dr& = − ∑
4
2
n =1 λп
0
exp(−λ n x& ). (18)
Этот дефицит расхода невелик, но точное решение требует его устранения,
для чего используется дифференциальное уравнение распределения аксиальных скоростей
∂ 2 u& х
∂ 2 u& х
∞
J (λ r&)
∂u& x
= −8 + 8 ∑ 0 n exp(−λ n x& ), r&∂r&
λ J (λ n )
n =1 п 1
(19)
∂x& 2
∂r& 2
которое легко получить из уравнения азимутального вихря (1) и уравнения неразрывности
+
+
∂ (r&u&r ) ∂u& x
(20)
+
= 0. r&∂r&
∂x&
Причем для получения требуемого решения следует использовать однородную часть уравнения (19) в виде
∂ 2 u& х
∂u& x
= 0,
r&∂r&
∂x&
∂r&
поскольку неоднородное уравнение дает результат, соответствующий (17).
Принимая как и ранее дополнительные (компенсационные) осевые скорости произведением двух функций u& х (r&, x& ) = φ(r&)ϕ( x& ) , одна из которых φ(r)
зависит только от текущего радиуса, а вторая ϕ(x ) — только от текущей аксиальной координаты, находим
2
130
+
∂ 2 u& х
2
+
ISSN 1997-0935. Vestnik MGSU. 2014. № 5
Гидравлика. Инженерная гидрология. Гидротехническое строительство
2
 ∂ 2 φ ∂φ 
 2


 = 0, или, деля на φϕ, имеем 1 ∂ ϕ = − 1  ∂ φ + ∂φ  .
+
ϕ
+
 ∂r& 2 r&∂r& 
ϕ ∂x& 2
φ  ∂r& 2 r&∂r& 
∂x& 2


Последнее равенство, в котором правая часть является произведением
функций от x , а левая — произведением функций только от r, возможно только в случае если одновременно ни правая, ни левая части не зависят ни от x ,
ни от r , т.е. обе части этого равенства представляют собой одну и ту же постоянную h. Следовательно,
φ
∂2 ϕ
1 ∂2 ϕ
1  ∂ 2 φ ∂φ 
=
−
+
= η,
ϕ ∂x& 2
φ  ∂r& 2 r&∂r& 
что позволяет записать систему двух дифференциальных уравнений

∂ 2ϕ
− hϕ = 0;

2
∂x

∂ 2 φ ∂φ
+
+ hφ = 0.

∂r 2 r∂r
(21)
Можно видеть, что при η = λ2n > 0 решением первого уравнения системы
φn Bn J 0 ( l n r ) .
(21) будет являться ϕn = С exp ( −λ n x& ) , а решением второго =
Таким образом, получаем дополнительные (компенсирующие дефицит
расхода) аксиальные скорости в виде ряда Фурье — Бесселя
u& x (r&, x& ) =
∞
∑ Bn J 0 (λ n r&) exp(−λ n x& ), (22)
n =1
где константа С вошла в константу Вп.
Можно видеть, что эти дополнительные скорости не изменят нулевых
значений аксиальных скоростей у стенок трубы, ибо J 0 (λ n ) = 0 , однако расход от этих скоростей должен быть обратным равенству (18), т.е. должен
составлять
1
∞
1
∫ u& x 2r&dr = 2 ∑ Bn exp(−λ n x& ) ∫ r&J 0 (λ n r&)dr& =
∞
∑
4
exp(−λ n x& ).
2
λ
n =1
n =1 п
0
0
Результаты этих вычислений позволяют получить значение константы Вп
Вn =
2
.
λ п J1 (λ n )
(23)
Подставляя (23) в (22), находим компенсационную скорость
∞
J 0 (λ n r&)
exp(−λ n x& ), λ J (λ n )
n =1 п 1
u& x (r&, x& ) = 2 ∑
(24)
вводя которую в выражение (17), находим окончательное распределение аксиальных скоростей в неравномерном ползущем течении в начале трубы
u& x (r&, x& ) = 2(1 − r& 2 ) −
 4  r&J1 (λ n r&)   4
 2 J 0 (λ n r&) 


−
1
1
exp(−λ n x& ). (25)
+
−
  2
2 
 λ n J1 (λ n ) 
λ
J
(
)

1
n
  λп
n = 1 λ п 


∞
∑
Hydraulics. Engineering hydrology. Hydraulic engineering
131
5/2014
Найдем окончательное распределение функции тока, взяв интеграл
ψ
&
r
∫ Ψ& = ∫ u& x r&dr ,
0
0
в результате получаем
2 
∞

& (r&, x& ) = r& 2 1 − r&  − 2r& ∑ 1 r& − J1 (λ n r&) + 2r&J 0 (λ n r&)  exp(−λ x& ). (26)
Ψ

n
2 

J1 (λ n ) λ n J1 (λ n ) 
2 
n =1 λп 

Радиально-аксиальное распределение радиальных скоростей найдем дифференцированием функции тока (26) по аксиальной переменной

J1 (λ n r&) 2r&J 0 (λ n r&) 
(27)
+
r& −
 exp(−λ n x& ). J1 (λ n ) λ n J1 (λ n ) 

Теперь полученные радиально-аксиальные распределения характеристик
неравномерного ползущего течения в круглой трубе (25)—(27) отвечают всем
граничным условиям (5)—(7).
Отметим также, что, как показывают выполненные расчеты, гидравлические потери при равномерном ползущем течении ньютоновской жидкости соответствуют формуле Пуазейля — Гагена
u&r (r&, x& ) = −
∞
&
1
Ψ
= −2 ∑
&rdx&
λ
n =1 n
64 (28)
,
Re
где l — коэффициент гидравлического сопротивления по длине трубы; Re —
число Рейнольдса, Re = Vd/n; d — диаметр трубы, d = 2R; n — кинематическая
вязкость.
Нормированные графики аксиальных и радиальных скоростей исследованного течения приведены на рис. 1 (здесь профиль Пуазейля показан штриховой кривой). Расчеты выполнены для створов, расположенных от входа в
трубу на расстояниях х = 0,1R, 0,2R, 0,4R и 0,8R.
λ=
Рис. 1. Радиально-аксиальное распределение аксиальных u x и радиальных ur
скоростей в неравномерном ползущем течении на начальном участке трубы
Можно видеть, что согласно правому рисунку в ползущем течении аксиальные скорости u& x приобретают профиль Пуазейля в самом начале трубы на
расстоянии x близком к 0,8 ее радиуса R. При этом радиальные скорости ur ,
показанные на левой части рисунка, всегда отрицательны, т.е. направлены в
сторону оси трубы.
132
ISSN 1997-0935. Vestnik MGSU. 2014. № 5
Гидравлика. Инженерная гидрология. Гидротехническое строительство
На рис. 2 показаны изолинии нормированной функции тока.
 в неравномерном ползуРис. 2. Нормированное поле изолиний функции тока Ψ
щем течении на начальном участке трубы
В целом отметим, что структура неравномерного ползущего течения весьма близко напоминает структуру неравномерного ламинарного течения на
начальном участке трубы. Разница лишь в том, что если длина начального
участка, на котором имеет место неравномерное ламинарное течение, ограничивается длиной от начала трубы, равной x = 0,175ReR [1], то для ползущего течения вне зависимости от числа Рейнольдса она составляет порядка
x = 1,0R.
Библиографический список
1. Орехов Г.В., Зуйков А.Л., Волшаник В.В. Контрвихревое ползущее течение //
Вестник МГСУ. 2013. № 4. С. 172—180.
2. Моделирование и расчет контрвихревых течений / В.К. Ахметов, В.В. Волшаник,
А.Л. Зуйков, Г.В. Орехов ; под ред. А.Л. Зуйкова. М. : МГСУ, 2012. 252 с.
3. Зуйков А.Л. Азимутальный вихрь и функция тока в ползущем течении в трубе //
Вестник МГСУ. 2014. № 4. С. 150—159.
4. Korn G.A., Korn T.M. Mathematical Handbook for Scientists and Engineers:
Definitions, Theorems, and Formulas for Reference and Review. New York : General
Publishing Company, 2000. 1151 p.
Поступила в редакцию в январе 2014 г.
О б а в т о р е : Зуйков Андрей Львович — доктор технических наук, профессор кафедры гидравлики и водных ресурсов, Московский государственный строительный
университет (ФГБОУ ВПО «МГСУ»), 129337, г. Москва, Ярославское шоссе, д. 26,
8(495) 287-49-14 вн. 14-18, [email protected]
Д л я ц и т и р о в а н и я : Зуйков А.Л. Аксиальные и радиальные скорости в ползущем
течении в трубе // Вестник МГСУ. 2014. № 5. С. 127—134.
Hydraulics. Engineering hydrology. Hydraulic engineering
133
5/2014
A.L. Zuykov
AXIAL AND RADIAL VELOCITIES IN THE CREEPING FLOW IN A PIPE
The article is devoted to analytical study of transformation fields of axial and radial
velocities in uneven steady creeping flow of a Newtonian fluid in the initial portion of the
cylindrical channel. It is shown that the velocity field of the flow is two-dimensional and
determined by the stream function. The article is a continuation of a series of papers,
where normalized analytic functions of radial axial distributions in uneven steady creeping flow in a cylindrical tube with azimuthal vorticity and stream function were obtained.
There is Poiseuille profile for the axial velocity in the uniform motion of a fluid at an infinite
distance from the entrance of the pipe (at x = ∞), here taken equal to zero radial velocity. There is uniform distribution of the axial velocity in the cross section at the tube inlet
at x = 0, at which the axial velocity is constant along the current radius. Due to the axial
symmetry of the flow on the axis of the pipe (at r = 0), the radial velocities and the partial
derivative of the axial velocity along the radius, corresponding to the condition of the soft
function extremum, are equal to zero. The authors stated vanishing of the velocity of the
fluid on the walls of the pipe (at r = R, where R — radius of the tube) due to its viscous
sticking and tightness of the walls. The condition of conservation of volume flow along
the tube was also accepted.
All the solutions are obtained in the form of the Fourier — Bessel. It is shown that
the hydraulic losses at uniform creeping flow of a Newtonian fluid correspond to Poiseuille — Hagen formula.
Key words: Newtonian fluid, uneven creeping flow, azimuthal vorticity, stream function, axial and radial velocities, decomposition of Fourier — Bessel.
References
1. Orekhov G.V., Zuykov A.L., Volshanik V.V. Kontrvikhrevoe polzushchee techenie
[Creeping Counter Vortex Flow]. Vestnik MGSU [Proceedings of Moscow State University of
Civil Engineering]. 2013, no. 4, pp. 172—180.
2. Akhmetov V.K., Volshanik V.V., Zuykov A.L., Orekhov G.V. Modelirovanie i raschet
kontrvikhrevykh techeniy [Modeling and Calculation of Counter Vortex Flows]. Мoscow,
MGSU Publ., 2012, 252 p.
3. Zuykov A.L. Azimutal'nyy vikhr' i funktsiya toka v polzushchem techenii v trube [Azimuthal Vorticity and Stream Function in the Creeping Flow in a Pipe]. Vestnik MGSU [Proceedings of Moscow State University of Civil Engineering]. 2014, no. 4, рp. 150—159.
4. Korn G.A., Korn T.M. Mathematical Handbook for Scientists and Engineers: Definitions, Theorems, and Formulas for Reference and Review. New York, General Publishing
Company, 2000, 1151 p.
A b o u t t h e a u t h o r : Zuykov Andrey L'vovich — Doctor of Technical Sciences, Professor, Department of Hydraulics and Water Resources, Moscow State University of Civil
Engineering (MGSU), 26 Yaroslavskoye shosse, Moscow, 129337, Russian Federation;
+7 (495) 287-49-14, ext. 14-18; [email protected]
F o r c i t a t i o n : Zuykov A.L. Aksial'nye i radial'nye skorosti v polzushchem techenii v trube
[Axial and Radial Velocities in the Creeping Flow in a Pipe]. Vestnik MGSU [Proceedings of
Moscow State University of Civil Engineering]. 2014, no. 5, pp. 127—134.
134
ISSN 1997-0935. Vestnik MGSU. 2014. № 5
1/--страниц
Пожаловаться на содержимое документа