Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Тепловые и гидродинамические процессы в колеблющихся потоках

..pdf
Скачиваний:
9
Добавлен:
12.11.2023
Размер:
13.64 Mб
Скачать

Рис. 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

 

Г

 

 

 

(246)

Дтиго

Го

V

1+ *

{

утр!

гоcth

(1 -fi)

Для высокочастотных колебаний бк < г0 и cth

1(1 +

i) rJbK]—>

—» 1; следовательно,

 

 

 

 

 

 

 

 

 

 

 

А (Р“)о/

 

1

 

 

 

В этом случае безразмерный коэффициент потерь

 

 

т"

- г Ь - 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 и,

следовательно,

PQQ/V _

J___ГоРо

 

(258)

Voro

3

Vo

 

 

 

100

Соседние файлы в папке книги