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

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

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

1. Колебания потока направлены вертикально и совпадают с потоком естественной конвекции (рис. 77). В этом случае вто­ ричные вихревые течения и поток естественной конвекции совпа­ дают по направлению в верхней части цилиндра и направлены противоположно у нижней половины его поверхности.

Для малых амплитуд колебаний скорости и при малых значе­ ниях чисел Pr ^ 1 критериальные уравнения для средней тепло­ отдачи по поверхности цилиндра Nu и для теплоотдачи в окрест­ ности нижней лобовой точки Nu0 согласно работе [33] имеют вид

Nu2 =

0,4 Re»Pr -J- j/" R ^ P r2+

PrGr ,

(376)

где Re0 = A«2/(ov

(Ди — амплитуда

колебания

скорости);

Nuo =

) 2 = Re»Pr + У

Re^Pr2 + A

PrGr .

(377)

При малой интенсивности колебаний (малые значения

ReoPx)

значение числа Nu стремится к соответствующему значению для стационарной естественной конвекции, а при большой интенсив­ ности колебаний теплообмен в основном осуществляется посред­ ством колеблющегося потока. При этом число Nu стремится к соот­ ветствующему значению числа Nu для колеблющегося потока без учета влияния естественной конвекции.

2. Колебания потока направлены перпендикулярно плоскости симметрии направления потока естественной конвекции (рис. 77). В этом случае направление вторичных вихревых течений и на­ правление свободной конвекции совпадают у нижней половины поверхности цилиндра и противоположны у верхней его половины. Такая схема течений приводит к дополнительной турбулизации пограничного слоя в нижней части цилиндра в результате взаимо­ действия вторичных вихревых течений и потоков свободной кон­ векции, которые в данном случае в нижней части цилиндра на­ правлены навстречу друг другу. Согласно экспериментальным данным работы [51 ] интенсивность теплообмена оказывается значительно больше, чем в предыдущем случае, и для высоко-

Рис. 77. Схема течения около цилиндра в зависимости от иаправления колебаний:

/ — направление потока естествен­ ной конвекции; 2 —•направление вторичных вихревых течений; 3 — направление колебаний потока

171

Рис. 78. Изменение локального коэффициента теплоотдачи а (ккал/м.1°С) по поверхности цилиндра (см. рис. 77):

--------- расчеты

частотных колебаний малой ампли-

туды при Рг < 1 обобщается

зави-

30 симостью

 

Nu = 2,18(Re0Pr)°-25Gr°>278.

 

В этом случае с увеличением

интенсивности колебаний (Re„

уве­

личивается) и с увеличением числа

Gr интенсивность теплообмена

уве­

личивается.

 

Исследование теплообмена от горизонтального цилиндра при озвучивании его вдоль оси показало, что в стоячих акустических волнах коэффициент теплоотдачи увеличивается в два с лишним раза [51 ].

Распределение локальных коэффициентов теплоотдачи по по­ верхности цилиндра согласно опытам, описанным в работах [51, 33], представлено на рис. 78.

Теплообмен при колебаниях в замкнутом объеме

Экспериментальное исследование влияния колебаний в замкнутом объеме на естественную конвекцию проведено в ра­ ботах [27, 36]. Экспериментальная камера, образованная двумя вертикальными пластинами с различным отношением высоты Я к ширине зазора между пластинами В (Н/В = 9,4-5-42,7), подвер­ галась вибрации [36] в вертикальном направлении с частотами О—400 Гц и с ускорениями 0—110g. В результате визуального наблюдения пограничного слоя на горячей и холодной пластинах установлено, что в зависимости от частоты колебаний погранич­ ный слой на пластинах может быть как ламинарным, так и турбу­ лентным. В области частот, близких к первой резонансной гар­ монике, наблюдается турбулентный пограничный слой, при зна­ чительном отклонении от резонанса — ламинарный и смешанный (на определенном расстоянии ламинарный слой переходит в тур­ булентный). В работе получено существенное увеличение коэф­ фициента теплоотдачи при вибрациях в диапазоне резонансных частот колебаний. Причиной, вызывающей увеличение коэффи­ циентов теплоотдачи, вероятно, является развивающаяся тур­ булентность пограничного слоя по всей поверхности замкнутого объема, которая была тем значительней, чем ближе частота выну­ жденных колебаний совпадала с резонансом (собственной часто­ той колебаний столба жидкости в камере). Параметрами, оказы­ вающими влияние на теплоотдачу, являются частота колебаний

172

и ускорение. Экспериментальные данные [36 ] аппроксимируются уравнением

Nu = 1,117

(378)

где Ra — число Релея (в качестве характерного размера выбрана ширина зазора между вертикальными пластинами В); s — функ­ ция параметров вибрации, определяемая из выражения

 

с -L

S =

, In (1 -f- gt) C2;

здесь константа Ct — коэффициент, зависящий от теплофизиче­ ских свойств жидкости и параметров вибрации; f — частота коле­ баний; /„ — резонансная частота колебаний, соответствующая первой (основной) резонансной частоте; Р* — коэффициент, ха­

рактеризующий затухание колебаний в жидкости (в случае воды Р* = 0,7327).

Коэффициенты С, и С2 определялись для воды из эксперимен­ тальной зависимости Nu = / (Ra) при отсутствии вибрации. Принимая

, = ____________ ///«____________

критериальное уравнение зависимости теплоотдачи при есте­ ственной конвекции в замкнутом объеме в присутствии верти­ кальных колебаний объема можно записать так:

Nu = 1,117 ( - ^ ) cizln(I + *') + S

(379)

где gt = AAo>2lg0 — безразмерное ускорение пластин

камеры;

АЛ — амплитуда колебаний; со — круговая частота; g0 — уско­ рение силы тяжести.

Для объема, образованного двумя вертикальными пластинами и заполняемого водой, критериальное уравнение для расчета

теплоотдачи имеет

вид

 

Nu =

1,117 (_R£.)0,03l06zln 0+«<)-И-3042^

(380)

Уравнение (380) можно иногда рекомендовать и для других жидкостей, так как параметры и Сх мало влияют на крите­ рий Nu. Результаты экспериментального исследования влияния колебания на теплоотдачу от горизонтального цилиндра в замкну­ том объеме (вода), когда вся система (и поверхность нагрева и жидкость) подвергалась вибрации, приведены в работе [27].

173

 

 

 

 

 

Рис. 79.

Зависимость

относи­

 

 

 

 

J

тельного

коэффициента

тепло­

 

 

 

 

отдачи К от

критериального

 

 

 

 

m

 

 

R e ^ p /6

#

 

 

7 ОМ

 

комплекса

 

 

 

 

 

M

 

 

 

Ra°’26Re?.’4

 

W

Ra^Ref

 

#

— ф =

0*15;

а

— Ф =

15*20;

 

 

 

Ао Ш у ‘ °

О

— ф =

20 * 25;

д

— ф =

25в? 30

 

 

 

 

 

 

 

 

 

 

 

 

 

U‘ 25

 

 

 

 

 

 

 

 

 

 

 

В

качестве

поверхности

 

 

 

 

 

нагрева

использовалась

о т о,oi

 

0,1

 

1,0 2,0

горизонтально

располо­

 

 

 

 

Ке&Рг9*

женная

платиновая

про­

 

 

 

 

'$1

волока диаметром 0,81 мм,

 

 

 

 

 

помещенная в бак с водой.

Экспериментальные данные по относительному коэффициенту

теплоотдачи

К обобщаются следующей зависимостью:

 

 

К

= 1 + 2 ,6 4

Re^Pr0-6

 

 

/0,46

 

 

(381)

 

Ra0,26Re®’4 - И

- Ю

' 8 Ra0,2IRe<U

 

 

Последние два члена в правой части^уравнения (381), характе­ ризующие эффекты относительного движения и ускорения, учи­ тываются с помощью критериев ReAu и Re,*.

Установлено, что влияние относительного движения гораздо больше, чем влияние ускорения.

Один из важных выводов, следующих из анализа уравне­ ния (381), заключается в том, что среднее значение чисел Нуссельта Nu при колебаниях уменьшается по мере увеличения числа Rea. Поскольку частота колебаний в числителе в большей степени, чем в знаменателе, Nu увеличивается с увеличением частоты колебаний. Фактором, способствующим уменьшению Nu, яв­

ляется

число Релея

(Ra = PrGr).

Максимальное

увеличение

числа

К, равное

 

Re^Pr*6

2,4, наблюдается при —^

—ду = 4 н- 5

 

y0,46‘

 

Ra

'

Re^,’

й при

0-5-30. Влияние

параметра

вибрационного

— g-jj—g-y =

Ra • Rew’

ускорения весьма небольшое и в условиях проведенных экспери­ ментов не превышало 12%. Основное влияние оказывает относи­ тельное движение.

Сравнение экспериментальных данных в работе [27] с рас­ четами по уравнению (381) для относительного коэффициента теплоотдачи приведено на рис. 79.

Г Л А В А V

ТЕПЛОВЫЕ И ГИДРОДИНАМИЧЕСКИЕ ПРОЦЕССЫ В ТУРБУЛЕНТНЫХ КОЛЕБЛЮЩИХСЯ ПОТОКАХ

1. ВОЗДЕЙСТВИЕ ВЫНУЖДЕННЫХ КОЛЕБАНИЙ НА УСТОЙЧИВОСТЬ ЛАМИНАРНОГО ТЕЧЕНИЯ

Колебания давления и скорости, наложенные на осредненное течение, могут служить причиной возмущающего воздей­ ствия на устойчивость и структуру течения. В случае ламинар­ ного течения эти воздействия могут ускорить или замедлить пере­ ход к турбулентному режиму течения, а в случае турбулентного потока — существенно изменить гидродинамические характери­ стики потока, в частности спектр турбулентности, и даже при­ вести к образованию вихревых течений.

Вопрос о влиянии возмущений на устойчивость течения в экспе­ риментальном и теоретическом аспектах рассматривался в ра­ боте [41].

Основная информация об устойчивости 'пульсирующих те­ чений получена в виде экспериментальных зависимостей и ча­ стных теоретических решений.

Наиболее разработанной является теория устойчивости погра­ ничного слоя. Поле скоростей пульсирующего течения в первом приближении обладает свойствами, характерными для колеблю­ щегося пограничного слоя, который возникает как при есте­ ственных, так и при вынужденных возмущениях в процессе перехода.

Исследование устойчивости возмущенного движения можно проводить двумя методами: энергетическим и методом малых колебаний.

Первым методом определяется изменение во времени энергии возмущающего воздействия. Вопрос об устойчивости в этом слу­ чае решается в зависимости от того, увеличивается или умень­ шается во времени энергия возмущающего воздействия. Во вто­ ром методе исследуется развитие во времени возмущенного поля скоростей. Если со временем амплитуда возмущения скорости нарастает, — течение окажется неустойчивым, если затухает, — течение устойчиво. Более широкое распространение получил второй метод.

Наиболее разработанной считается линейная теория неустойг чивости, которая предполагает, что наложенные на стационарное ламинарное течение возмущения параметров потока малы по сравнению с величинами осредненного потока.

175

Для двумерного возмущенного поля скоростей вводится функ­ ция тока

 

Ф(*. у,

t) = ф (у) е{ <«*-•«,

 

(382)

где ф =

Фг + 1ф,- — комплексная амплитуда

функции

тока

а =

= 2п!Х

(А, — длина

волны

возмущения);

р = Рг +

i‘P, — ком­

плексный коэффициент (Рг — круговая

частота отдельного

коле­

бания,

р,- — коэффициент

нарастания

колебания).

< 0

коле­

В случае рг > 0

течение неустойчиво, в случае рг

бание затухает и ламинарное течение устойчиво. Амплитуда возмущающего воздействия ф зависит от у.

Составляющие скорости возмущающего воздействия имеют вид

и =

дф

 

 

 

~W

 

(383)

,

дф

.

,

v = —

=

iaw,„\e‘

 

»аф(й)е<

Уравнения движения относительно возмущений в линейном приближении можно записать в следующем виде:

Подставив уравнение (383) в формулу (384) и приведя пара­ метры к безразмерному виду, получим уравнение Ора-Зоммер- фельда:

с) (ф" — а*ф) — й'ф =

(ф"" — 2а*ф" + а 4ф), (385)

где с = р/а = с, + ict (сг — фазовая распространения волны возмущения; щая судить о поведении колебания:

скорость, т. е. скорость ct — величина, позволяю­ ct > 0 — колебание нара­

стает, Ci < 0 — колебание затухает); и — безразмерная скорость; Ree — число Рейнольдса, определенное по толщине погранич­ ного слоя (Ree = иб/v).

Индексы (штрихи) у ф и и означают порядок производной. Уравнение Ора-Зоммерфельда решается при известных гранич­ ных условиях. Члены левой части уравнения (385) получены из инерционных членов уравнения движения, а члены правой части— из членов, учитывающих трение. Если вязкие силы малы, т. е. значения чисел Ree велики, то уравнение (385) можно упростить, отбросив в нем все члены правой части. Правомерность такого упрощения в первом приближении обоснована экспериментально, так как значения критических чисел Рейнольдса, при которых достигается предел устойчивости ct = 0, достаточно велики. Для

176

плоского ламинарного течения амплитуда малых возмущений будет удовлетворять дифференциальному уравнению возмущен­ ного движения без трения

 

(и — с) (<р"— а*<р)— й"<р = 0.

(386)

Полученное дифференциальное уравнение второго порядка

может

удовлетворять

двум граничным условиям:

ф = 0 при

у = 0

и ф = 0 при

у = оо.

 

Исследование устойчивости сводится к определению собствен­ ных значений дифференциального уравнения (386) при извест­ ных*граничных условиях. Найдя собственное значение функции ф и собственное значение коэффициента с, можно отыскать условия, соответствующие нейтральному '(безразличному) воздействию, т. е. условия, при которых возмущения не затухают и не усили­ ваются (С[ = 0). Эти условия будут соответствовать границе между устойчивой и неустойчивой областями течения.

Первым и весьма важным результатом, полученным на осно­ вании решения уравнения (386), является теорема о влиянии на устойчивость точки перегиба на профиле скоростей. Согласно этой теореме профили скорости, имеющие точку перегиба, яв­ ляются неустойчивыми. Расчеты на основе уравнения (386) для

профилей (см.

работу

[41 ])

показали, что выпуклые

профили

(рис. 80, а, б)

обеспечивают

устойчивость движения,

а

профили

скорости

(рис. 80, в, г)

приводят к

неустойчивости.

устойчивы,

При

R e —>оо профили скорости

типа а, б,

д, е

профили скорости типа в, г,

ж — неустойчивы.

Профиль типа д

получается при падении давления, профиль типа е — при постоян­ ном давлении, типа ж — при повышении давления в направлении течения. Расчеты в работе [41 ] дали основание считать, что про­ фили скорости с точками перегиба являются неустойчивыми

(рис. 80, ж, з). Позднее

это было доказано

и с учетом

вязкости.

Уравнение движения несжимаемой жидкости, выраженное

через напряжение Рейнольдса, имеет вид

 

 

—__ 1

дР

I

1

д („ dui

(387)

dt ~T~Ui дх,

р

дх{

^

р

ад**

дх

 

Напряжение Рейнольдса (щи, ) как дополнительное напряжение к силам давления и вязкого напряжения оказывает дополнитель­ ное влияние на осредненное течение. Если напряжение передает энергию от основного течения к возмущению, то это может вы­ звать неустойчивость. В работе [41 ] показано, что наличие этого напряжения благоприятствует переходу энергии осредненного движения в энергию возмущенного течения. Обмен энергией между основным течением и наложенными возмущениями яв­ ляется одним из физических механизмов, который используется как в теории турбулентности, так и в теории устойчивости лами­ нарных течений.

12 б . М. Гал.чцейский

177

-/(I

Рис. 80. Профили скоростей в ламинарном^ пограничном слое:

и т

скорость внешнего течения; б —

толщина пограничного слоя; б* — толщина вытес-

нения;

W — точка перегиба профиля

скорости

Характеристики устойчивости пограничного слоя на пластине приведены на рис. 81. По оси ординат отложена безразмерная длина волны возмущения аб* (б* — толщина вытеснения погра­ ничного слоя), а по оси абсцисс — число Рейнольдса, опреде­ ленное по толщине вытеснения; Точки, лежащие в областяхвнутри нейтральных кривых, определяют состояние движения, соответствующее неустойчивым колебаниям, точки вне нейтраль­ ных кривых — состояние, соответствующее устойчивым колеба­ ниям, а точки, лежащие на самих нейтральных кривых, — со­ стояние, соответствующее нейтральным колебаниям. При зна­ чительном увеличении чисел Рейнольдса Ree обе ветви нейтраль­ ных кривых приближаются к оси абсцисс. Наименьшее число Рейнольдса, при котором нейтральное возмущение возможно, (Ree)KP = 420.

Таким образом, при возмущениях, когда длина волны возму­ щения находится за некоторым нижним пределом (аб*)пред, всегда существует конечный интервал чисел Ree, в котором тече­ ние неустойчиво. В частности, движение полностью устойчиво при малых числах Рейнольдса.

178

Рис.

81.

Нейтральные

 

кривые,

■а$*

 

 

определяющие

функциональную

 

 

 

 

зависимость между длиной волны

 

 

 

 

возмущений аб* и числом Рей­

 

 

 

 

нольдса Re, для пограничного слоя

 

 

 

 

на плоской пластине

 

 

 

 

 

 

 

Положение

нейтральной

 

 

 

 

точки и . точки перехода ла­

 

 

 

 

минарного

течения

в турбу­

 

 

 

 

лентное определяется

интен­

 

 

 

 

сивностью

нарастания

неус­

 

 

 

 

тойчивых

возмущений

и сте­

 

 

 

 

пенью турбулентности внеш­

 

 

 

 

него

течения.

 

 

 

 

 

 

 

 

С учетом вязкости область

 

 

 

 

неустойчивости

на

рис.

81

 

 

 

 

сокращается и располагается

 

 

Автор работы*

[101

внутри

зоны

«невязкой»

неустойчивости.

определил

собственную -функцию

ф(4,> для

нейтральных

коле­

баний

(см. рис. 81,

точка /). Определенные по значению

этой

функции

линии

тока возмущенного

движения выглядят

так,

как

показано

 

на

рис.

82.

Многие

экспериментальные

работы при вынужденных возмущениях пограничного слоя пока­ зали хорошую сходимость результатов, характерных естествен­ ным возмущениям, с данными рис. 81. Однако при этом'выясни­ лось, что нарастание неустойчивых волн приводит к явно выра­ женной трехмерной структуре течения.

Влияние акустических колебаний на развитие и переход ламинарного пограничного слоя в турбулентный исследовалось экспериментально. В условиях акустических колебаний как при естественных, так и при вынужденных возмущениях в процессе

X

 

 

 

 

 

критический

 

 

 

 

 

C/WU\

 

 

----------W.

 

х т

]*■«•

«

f

'■

_

} ’

■ №

Рис. 82. Картина линий тока и распределение скоростей для нейтрального коле­ бания в пограничном слое на продольно-обтекаемой плоской пластине:

и у

— основное

течение; и

« ' (<*.

У» t ) * возмущенное

поле скоростей;

Re ®

= и

6*/v =893 — число Рейнольдса; X =

4,0; б* — длина волны

возмущения; А и

0

= 0,36

 

СО'

распространения

волны.возмущении,

 

и т — скорость

 

 

 

12*

179

Рис. 83. Зависимость влия­ ния частоты акустических колебаний на устойчи­ вость пограничного слоя:

в ) Д с , = F (f); б) i = F (« oe * / v )

перехода ламинарного пограничного слоя в турбулентный воз­ никает волновой пограничный слой с периодическими колебаниями.

При озвучивании пограничного слоя на

частотах,

лежащих

в области неустойчивости (ct > 0 ), т. е. на

частотах,

близких

к частоте естественных вихрей, которые и являются причиной возникновения волнового пограничного слоя, наблюдается син­

хронизация вихрей

с вынужденной частотой f,

что

приводит

к усилению возмущений (средняя часть кривой

Дct =

F (f) на

рис. 83). В области

низких частот (зона 1 и Дс, =

F (/)) акусти­

ческие колебания малой амплитуды не оказывают заметного влия­ ния на развитие пограничного слоя, а в области частот выше области нестабильности (зона 2) наблюдается подавление вихреобразования (зона 3), что замедляет переход ламинарного со­ стояния пограничного слоя. Увеличение уровня звукового сиг­ нала, соответствующего области неустойчивости, приводит к уве­ личению синхронизирующего воздействия внешних колебаний. При этом синхронизация наблюдается при большей разности частот колебаний естественных вихрей и вынужденной частоты. Исследование влияния энергии и спектра акустических возму­ щений на переход ламинарного пограничного слоя в турбулент­ ный показало, что при воздействии распространяющихся вдоль потока звуковых колебаний достаточно большой интенсивности резко возрастает турбулентность продольной составляющей ско­ рости, в то время как интенсивность турбулентности поперечных составляющих пульсаций скорости остается неизменной. Резуль­ таты влияния акустических возмущений на переход ламинарного пограничного слоя в турбулентный показаны на рис. 84, где пред­ ставлены спектры продольных составляющих скорости в фикси­ рованном положении пространства пограничного слоя. Воздей­ ствие звука ускоряет наступление перехода в пограничном слое при определенных условиях. При уровне звуковых возмущений (УЗД), равном 110—115 дБ, звуковое возмущение не влияет на переход при УЗД = 120-5-125 дБ; такое влияние обнаружено лишь на частоте 800 Гц, а при УЗД = 130-5-135 дБ, влияние зву­ ковых возмущений на переход имело место на частотах 200, 400 и 800 Гц. Наконец, при УЗД = 140 дБ переход ускоряется уже

180

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