книги / Тепловые и гидродинамические процессы в колеблющихся потоках
..pdf1. Колебания потока направлены вертикально и совпадают с потоком естественной конвекции (рис. 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