
книги / Тепловые и гидродинамические процессы в колеблющихся потоках
..pdfРис. 14. Распределение амплитуды |
(а) и фазы (б) колебаний скорости |
||
в пограничном слое |
|
|
|
Для этого случая уравнение (214) будет иметь вид |
|
||
<Pg |
- * ь й - 2 |
- С , |
(219) |
дт)8+ / о | ^ — (2/о+£) |
|||
где ДиГ> = Axeg' (rj) exp (Ш); vt = — У Avtg (т]) exp (Ш). |
урав |
||
Для малых частот колебания, т. е. при Sh < |
1, решение |
нения (219) можно искать в виде ряда (215). Для больших значе ний частоты колебания решение задачи можно искать в виде ряда
(216). Используя метод ВКБ |
[57] для больших значений чисел Sh, |
||||
получим асимптотическое решение |
|
||||
|
|
8 (Л) = (л — y = j + у = exp г (£), |
|
||
где |
z ( Q = — V l |
+YjTjP jМлj• |
|
||
|
Зная значение g" (0) при t] = 0, определим касательное на |
||||
пряжение на |
поверхности |
|
|
||
|
~ |
+ |
0,225t Sh + |
0,020 Sh2 (для малых Sh); |
(220) |
1 |
Atцт”! _ |
|
|
|
|
|
0 ,8 1 1 3 /rS h + |
(для больших Sh). |
(220а) |
91
А |
|
|
|
<Р' |
|
|
|
|
|
|
|
|
|
|
Рис. |
15. Зависимость амплитуды |
|||
|
|
|
|
|
А = 1/е | Axw/rWo | фазы ф колебания |
||||
|
|
|
|
|
касательного |
напряжения в окрест |
|||
|
|
|
|
|
ности |
критической |
точки от^ ча |
||
|
|
|
|
|
стоты |
(Sh): |
|
|
|
|
|
|
|
|
1 — уравнение (220); |
2 — |
уравнение |
||
|
|
|
|
|
(220а) |
|
|
|
|
здесь Vo = |
'<К\ |
_ |
стационарное |
значение касательного |
|||||
fi(-a ; |
|
||||||||
напряжения; |
* |
|
/аДи<'> \ |
— пульсационное |
|
значение |
|||
Атш = р |
^ |
|
|
||||||
касательного |
напряжения; |
Shъп = |
0) |
/ <4 со |
|
л |
|||
—-j- — w / |
-----число |
Стру- |
|||||||
халя, определенное по градиенту средней скорости. |
|
пульса- |
|||||||
Зависимости амплитуды А |
и фазы ср |
относительного |
ционного трения на поверхности от частоты показаны на рис. 15. С увеличением частоты (числа Sh) амплитуда и фаза колебания касательного напряжения на поверхности увеличиваются, и при больших значениях Sh фаза колебаний достигает 45°.
Система уравнений во втором приближении соответственно за пишется:
дДи(2* . |
дДм<2) |
+ |
, |
А т |
ди0 . |
аД«42> |
. |
||
dt |
+ |
ал: |
|
Аы<2> |
дх |
+ |
У0- s f - + |
|
|
|
|
|
|
~ |
ду |
|
|||
|
дШ1' |
|
|
|
|
|
|
|
|
+ Ди<1> |
дх |
Дц£2> ^ |
+ |
Д ц < 1 > ^ = |
|||||
|
_ |
1 ГадиЗ«,] . |
, а2Ди<2> . |
(221) |
|||||
|
|
||||||||
|
“ |
2 [ |
дх |
J " r V |
а</2 |
’ |
|
||
|
|
ади«> |
|
, адг»«2) |
_ л |
|
|||
|
|
дх |
|
^ |
ду |
~ и - |
|
||
Граничные условия: Дм(2) = До(2) = |
0 при у = 0; Ды<2) = О |
при у — ОО.
Как следует из приведенной системы уравнений, нелинейные члены во втором приближении приводят к появлению гармоники
удвоенной |
частоты. |
|
|
Ах можно |
||
Решение системы уравнений (221) для случая, ыоо, = |
||||||
искать в виде |
|
|
|
|
|
|
|
Дutf) = |
Axe2{Go (г)) -f- G' (ri) ехр (2mt)}; |
1 |
|
||
|
Ду(2) = |
— YAve2{G0 (rj) + |
G(rj) exp (2mi)}. ) |
(222) |
||
Функция |
G0 (л) характеризует |
влияние |
пульсационного |
|||
движения |
на |
осредненное. Подставляя выражения |
(222) в |
92
уравнения (221), |
получим уравнение |
относительно |
функ |
ции G0 (п): |
< г ,+ /< Г о - 2 /;о ;+ /ч |
= |
|
|
|
||
= I |
№ ) 2 + ( g t f - g r g r - g i g l - 1}, |
(223) |
где gr и gt — соответственно мнимая и действительная части функ ции g (rj).
Граничные условия: Go = Go = 0 при t] = |
0; Go = |
0 при |
Г] —» оо. |
|
|
Для малых значений частоты колебания |
решение |
уравне |
ния (223) можно искать в виде ряда |
|
|
G0= t v l)Gi{44).
Для больших значений частоты колебаний в правую часть урав-' нения (223) необходимо подставлять функцию g (rj).
На, рис. 16 приведены расчетные значения функции Go, харак теризующие деформацию профиля осредненного движения под действием колебаний в пограничном слое. Как следует из при веденного рисунка, вблизи поверхности под действием колебаний скорость (осредненная по времени) потока жидкости увеличи вается, что приводит к увеличению градиента скорости, а следо вательно, и силы трения на стенке. С увеличением частоты коле баний (или критерия Sh) максимум осредненной по времени ско рости уменьшается и смещается в сторону поверхности стенки.
Выражение для осредненного (по времени) касательного на пряжения можно представить в виде £57]
= 1 + e 2K(Sh), |
(224) |
Vo
2 я
где (% s) = 2^- J p ( - ^ ) y=0d (©0 — осредненное по времени каса
тельное напряжение на поверхности; |
|
|
|||
F (Sh) = |
-jg — 0,0033Sh |
(для небольших Sh); |
|
(225) |
|
0.2868 . |
0,4932 . |
0,8313 , |
, |
(225а) |
|
|
Y s b + |
Sh- + |
(Для больших Sh). |
Рис. 16. Распределение функции G'0 в колеблющемся пограничном слое для различных значений Sh
93
F(Sh)
Рис. 17. Зависимость F (Sh) от частоты |
Рис. |
18. |
Зависимость |
F (Sh) |
от час |
||
в окрестности |
критической точки: |
тоты |
для |
плоской пластины: |
|||
/ — уравнение |
(225;) 2 — уравнение (225а); |
/ — уравнение |
(226); |
2 — |
уравнение |
||
3 — два члена |
уравнения (225а) |
(226а) |
|
|
|
|
|
Численные значения функции F (Sh) показывают, что с уве личением частоты (критерия Sh) влияние пульсаций скорости на осредненное по времени касательное напряжение уменьшается (рис. 17), а с увеличением относительной амплитуды колебания это влияние увеличивается, т. е. сила трения на поверхности уве личивается.
При обтекании колеблющимся потоком плоской пластины с по стоянной осредненной скоростью внешнего потока = const осредненное касательное напряжение на поверхности аппрокси мируется выражением (224), где
( |
— 0,621 (Sh)2 (для небольших Sh); |
(226) |
F (S h)= j |
|
|
( 0,1875 (Sh)-2 (для больших Sh). |
(226а) |
Значение функции F (Sh) представлено на рис. 18. С увеличе нием частоты колебания влияние колебаний потока на осредненную силу трения на поверхности уменьшается.
Сравнивая выражения (225) и (226), можно сделать вывод, что при одинаковых параметрах е и Sh влияние колебаний потока на трение на поверхности увеличивается с увеличением градиента скорости осредненного внешнего потока. Это связано с измене нием толщины стационарного пограничного слоя под действием градиента скорости.
В самом деле,
Т - = |
2’4 V Ь |
= 2)4 V т к (ДЛЯ “°“ = Ах); |
k = |
5 V k b ~ |
= 5 V ^ k (для M<>“ = const)- |
Влияние колебаний скорости внешнего потока с конечной амплитудой на осредненное движение в пограничном слое можно определить по уравнению (197),-в котором функция F (*, у), харак-
94
теризующая влияние пульсационного движения на осредненное, определяется согласно уравнению (203). /
Для случая линейного изменения скорости внешнего потока (ИОсо = Ах) функция
F (х, у) = - i. |
F. ( • £ ) - 4 - |
( £ ) . |
Тогда уравнения осредненного движения в пограничном слое можно записать так:
“•& + * % - М ' + -T»*. ( i ) ] + '-3?':
М -°-
Для решения этих уравнений удобно ввести переменные Блазиуса
«о = “о «.Go (л) = AxG'o (л); v0 = — VAvG (л).
Тогда уравнение движения относительно функции G0 (л) можно записать в виде
о; + G f i l - (G'Y + 1 + -i-eaF0 ( / Х ^ л ) = 0. |
(227) |
Граничные условия: Go = Go = 0 при л = 0; Go = 1 при л = °°- Безразмерную координату у/6к можно выразить через пере
менную Блазиуса л =
t - y / V
Таким образом, согласно уравнению (227) осредненное поле скоростей является функцией относительной амплитуды колеба ния скорости е и критерия Sh.
Осредненный по времени профиль скорости G' (л) =
и 0 со
в зависимости от е и Sh согласно численным расчетам приведен на рис. 19.
На рис. 20 приведена зависимость функции Go (0) для л = 0, характеризующей касательное напряжение на стенке канала. С увеличением относительной амплитуды колебания скорости внеш него потока и уменьшением числа Sh трение на поверхности уве личивается. При этом данные, полученные по методу, изложенному в работе [67], с погрешностью 0,5—3% совпадают с приведенным выше расчетом методом последовательных приближений.
При анализе колебания жидкости в канале можно использовать метод последовательных приближений, как и в случае расчета пограничного слоя. Рассмотрим двумерную модель течения жид.
95
Рис. 19. Распределение функций G0 (*)) и F0 (*]) в пограничном слое для различных значений относительной амплитуды и числа критерия Sh:
~ G'0 ( n ) ; -------------------------- F 0 (ц)
кости в канале, для которой уравнения движения и неразрыв ности (9) и (10) имеют вид
ди . |
ди |
, |
ди |
|
др . |
1 |
д (у«т). |
РаГ + Ры лГ + Р1'Ж |
= |
- Ж + |
Уп |
ду ’ |
|||
др |
, |
д_(ри) |
, |
_1_ |
д (упру) |
_ п |
(228) |
dt |
~Г |
дх |
Ч" |
уп |
ду |
~ |
’ |
где х — продольная координата; у — поперечная координата, на правленная вдоль радиуса канала (п = 0 — соответствует пло
скому каналу, п — 1 — цилиндрическому); т = Л ^ — напря
жение трения.
Будем считать, что давление и плотность постоянны по попереч ному сечению канала. Тогда, дифференцируя уравнения движе ния (228) по у и интегрируя уравнение неразрывности по попереч ному сечению канала с использованием граничных условий
у = 0 при у = 0; |
и — V = 0 при у = г0, J
3 |
6 |
9 |
12 15 Sh |
0 |
1 |
2 3 |
У е |
|
Рис. 20. |
Зависимость |
функции |
Go (0) |
от е |
при Sh = |
10 (о) |
||
и Sh = |
20 (б) |
и |
от Sh |
при е = |
2,0 (в) |
|
|
|
96
получим |
дх _ |
д |
Г |
1 д (упт) ~| ■ дР |
|
||
1 |
|
||||||
v |
dt |
ду |
L |
Уп |
ду |
J 'f 'P 'd y ’ ’ |
|
|
|
д£1 |
, |
д(ра); |
0; |
(230) |
|
|
|
dt |
~r |
~ d x ~ = |
|
||
|
дри . |
|
|
|
|
(231) |
|
|
+ |
|
|
|
|
||
|
|
|
|
|
|
||
здесь F — функция, характеризующая нелинейные члены в урав |
|||||||
нении движения |
|
|
|
|
|
|
|
|
F — |
|
( “ £ |
+ |
’> %) ■ |
<232> |
Индексом f обозначены средние по сечению значения плот ности р? и массовой скорости (gu)f.
1
Р/ = Х |
|
|
|
F |
|
Го |
(233) |
(р«), = T J |
(vu) dF |
tl -f- 1J |
(ры) уп dy. |
>o+1 о |
|
||
о |
|
|
Приведенную выше систему уравнений для пульсационного движения и осредненного по времени можно записать в виде:
для пульсационного движения
дДт |
|
д |
Г |
I |
д (i/лДт) 1 |
, |
дАР. |
|
(234) |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
дД (ри) |
, |
1 |
д (у*Д (рр)) |
_ |
дД (p«)f . |
|
(235) |
||||||
дх |
|
уп |
|
ду |
|
|
|
|
дх |
’ |
|
||
|
|
|
|
|
|
|
|
||||||
АГ |
/ А дДи |
, |
* <ЭД« |
+ |
|
ад“ |
|- , |
дД« \. |
(236) |
||||
AF = - |
( Аи |
|
|
|
|
и0Ж + |
о„ |
) , |
|||||
для осредненного |
движения |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
<2 3 7 > |
|
д^ |
1 |
= |
О |
[(ри), = |
const]; |
|
|
(238) |
||||
<’ — |
|
[ < * * т г > |
+ |
< 4' |
т г > ] ! |
(239) |
|||||||
|
|
ди0 |
, |
J _ |
|
|
|
_ |
|
|
|
|
|
|
|
д* |
' |
уп |
дУ |
|
|
|
|
|
|||
—Ж *£+*‘2 |
|
-+* ^ > ]- |
(240> |
7 Б. М. Галицейский |
^ |
При рассмотрении колеблющегося потока в канале при отсут ствии стационарного расхода необходимо принять (pu)of = 0.
Если плотность жидкости вдоль оси канала не изменяется, то уравнение неразрывности относительно пульсационных скоростей с точностью до линейных членов можно записать так:
д&и |
1 dynAv _ |
1 |
дА (рц)£ |
,2 4 П |
|
дх ‘ уп |
ду |
р0 |
дх |
' |
' |
Тогда уравнение неразрывности для осредненного движения можно записать в виде
ди0 . J _ dynv0 дх ' уп ду
Для синусоидальных колебаний градиента давления или мас совой скорости в канале пульсационное трение на стенке канала в первом приближении
Дтг = Д%0 ехр (Ш). |
|
(243) |
Пренебрегая нелинейными членами |
= 0 |
решение |
уравнения (234) относительно амплитуды колебания напряжения трения Дт0 с учетом граничного условия в центре канала (при
у — 0 т = р. |^- = 0^ для случая плоского канала можно записать
в виде
А ,. |
*Ь [ 1' + 0 | г ] |
(244) |
|
АХГ. |
* [ „ + „ £ ] • |
||
|
Распределение амплитуды колебания массовой скорости и сред нее ее значение по сечению канала определим посредством инте грирования выражения (244) с учетом граничных условий (229):
|
Д (ри)о |
_ |
f |
Д*о |
*У_ _ |
. J |
вк |
|
|
|
|
Ат yg |
|
J |
Дтц7о |
V |
v |
I —|—т |
|
|
|
|
|
|
г. |
|
|
|
|
|
|
|
|
,<!h[ (1 +i)^]— cth [(1 -ft)-* -] |
; |
(245) |
|||||||
|
sh [ ( 1 + 0 ^ |
] |
|
|
|
|
|
|||
A (P“)of |
1 |
1 |
|
6к |
Г |
|
|
|
(246) |
|
Дтиго |
Го |
V |
1+ * |
{ |
утр! |
гоcth |
(1 -fi) |
|||
Для высокочастотных колебаний бк < г0 и cth |
1(1 + |
i) rJbK]—> |
||||||||
—» 1; следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
А (Р“)о/ |
|
1 |
6К |
|
|
|
В этом случае безразмерный коэффициент потерь |
|
|
|||
т" |
- г Ь - V Т |
- - Т |
( - |
‘ f ) . |
(248) |
|
towо го |
|
|
|
|
где ReM= |
<ша"осil |
число |
Рейнольдса |
/ |
4 р |
—------колебательное |
( “э = |
~jf — |
= 4г0 — эквивалентный диаметр).
Для цилиндрического канала при высокочастотных колебаниях
коэффициент будет определяться также по формуле (248),
если в качестве характерного размера использовать эквивалент
ный диаметр d3 == ~ — d0.
Пульсационное поле поперечной скорости в первом прибли жении определим из уравнения неразрывности (235) с учетом гра ничных условий (229):
(249)
Для плоского канала в случае высокочастотных колебаний, когда 8К г0, получим
д ( й « “ >— « р [— 0 + 0 '* £ '■ ] ;
* £ £ 7 — г Ы ' - “ PO + O V ] - 1 <250>
Используя поле пульсационных скоростей в первом прибли жении, определим среднюю по времени нелинейную функцию
. . |
1 |
/ |
д (Р«)о/ |
(Р“ ) о / \ |
. |
(251) |
F (у!го)— |
о2 |
\\ |
п?Ро |
яг д* / |
ро (У/го)- |
Для этой цели выделим действительную часть в системе (250):
A (pU){ — cos a>t— exp [— — ^ j cos |
; (252) |
д (рц) _ V i |
|
£ C O S (^G)t— |
|
ад (p«)0f dx
(253)
Осредняя по времени выражение (251) и учитывая выраже
ния (252) и (253), получим |
' |
|
F0 (У/г0) = 1 + ехр [— 2 ^ = ^ ] — ехр [— г-°-“^ - ] |
х |
|
х [ 2Н |
т г ) + 81п( т г ) ] - |
(254) |
99
Распределение осредненного по времени касательного напря жения и скорости по радиусу канала определим посредством инте грирования уравнения движения (237). Полагая, что касательное напряжение на стенке канала известно, т. е.
|
при |
у = г0 |
т = % 0, |
|
ы = 0; |
(255) |
|
|
у — 0 |
|
ди0 |
||
|
при |
т = 0, |
|
ду»- = 0, |
|
|
в результате |
интегрирования |
уравнения (237) |
получим |
|||
То |
|
У |
|
|
1 |
|
гоРо |
1 j F (Y) Yn dY - |
|
Y j F (Y) Yn dYJ ; |
|||
Vo |
|
|
|
|
|
|
UjT = ~ |
f ( lT - ) dY = 4 V “ |
*") + |
|
|||
roro |
J \ w / |
L |
|
|
|
|
|
4- M i |(! __ Y 2) j F (Y) YndY — |
(256) |
||||
|
|
Интегрируя выражение (256) по поперечной площади канала, получим связь между касательным напряжением на поверхности канала и средней скоростью в канале:
PoVv |
= — (я -f- 1) |
|
YndY = |
|
V o ro |
|
|
|
|
|
|
U J O- \ Y "+UY = |
1 |
|
|
|
ti -j- 3 |
||
roPo |
|
F(Y)YndY — Y ^ F (Y )Y n+2dY] . (257) |
||
Vo |
|
|
|
|
Для плоского |
канала n = 0 и, |
следовательно, |
||
PQ“Q/V _ |
J___ГоРо |
|
(258) |
|
Voro |
3 |
Vo |
|
|
|
|
100