
книги / Тепловые и гидродинамические процессы в колеблющихся потоках
..pdfДля болыиихзначений частоты колебаний скорости (бK/r0 |
1) |
согласно выражениям (257) и (258) получим
(259)
хт го 3 % о °f дх / 3 го ’
Как следует из последнего выражения, касательное напряже ние на стенке канала складывается из соответствующего квазистационарного значения и пульсационного:
|
адut |
|
|
|
|
Of |
|
|
|
'W о |
\ AV дх / |
2v |
|
(260) |
1 = Р«о/ |
u0f |
V (О |
|
|
|
|
|||
Здесь в качестве характерного размера принят эквивалентный |
||||
диаметр плоского канала d3= 4г0. Величина £„ = т^0/ |
Р“о |
96 |
||
|
|
/ |
8 |
Re. |
является квазистационарным значением коэффициента сопроти
вления |
трению |
для |
плоского |
канала, |
последний член |
характе |
|
ризует |
влияние |
колебательного движения |
жидкости |
в канале |
|||
на осредненную по времени силу трения на стенке канала. |
|||||||
Для стоячей |
волны Auof |
А 0 sin (kx) и, |
следовательно, |
||||
|
|
/ |
dAuo t \ |
1 |
|
|
|
|
|
\ А и 0/ ~ ^ г / |
= ~2 -kAlsin 2kx. |
|
|||
Поскольку k = 2л/А = а>/а0, |
|
|
|
||||
то |
|
|
|
|
|
|
|
|
6 = |
|
- f e + 4 v "5 / T |
^ |
sln(2fa:)' |
|
|
|
P«0f = |
|
Таким образом, коэффициент сопротивления трению при воз буждении высокочастотных колебаний жидкости стоячей волной обратно пропорционален корню квадратному из акустического
а 2
числа Рейнольдса Rea = — .
3. НЕКОТОРЫЕ ОСОБЕННОСТИ СТРУКТУРЫ КОЛЕБЛЮЩИХСЯ ПОТОКОВ
Одной из основных особенностей колеблющихся пото ков является образование стационарных вихревых течений как в свободном колеблющемся потоке среды, так и вблизи препят ствий, помещенных в колеблющуюся среду, или вблизи колеб лющихся тел.
Как уже упоминалось выше, в 'теории колеблющихся потоков существенную роль играет толщина колеблющегося пограничного
Щ
слоя бк = ^ ^ У ’5 . Отношение толщины колеблющегося слоя к длине волны колебаний Л
А |
- \ 2 ( ±иЛ ( - V Л 0 ,5 |
Л |
" [ \(о Л А ДиЛ ) |
Поскольку о)А ~ 2nW и для малых колебаний скорость рас пространения волны W равна скорости звука (W ~ а), то послед нее соотношение можно представить в виде
•J [
Ля MA/ReA]°,S
где Мд = —----- колебательное число Маха; ReA = —------- ко
лебательное число Рейнольдса.
При ReA > Мд длина волны колебаний больше, чем толщина колеблющегося пограничного слоя, и в этом случае течение потока 'можно разделить на течение в пограничном слое, скорость которого зависит от вязкости, и течение вне пограничного слоя, скорость которого не зависит от вязкости среды. Когда ReA < Мд, четко выраженного пограничного слоя нет, и во всем пространстве ско рость течения зависит от вязкости.
Возникновение вихревых течений в колеблющихся потоках формально учтено нелинейными конвективными членами в урав нениях Навье-Стокса, значение которых может быть вычислено посредством определения функции F (х, у) в уравнении (197). Как следует из выражения (198), возникновение вихревых тече ний в значительной степени зависит от градиента скорости внеш него потока. Градиент скорости внешнего потока может быть обусловлен стоячей волной, например резонансными колебаниями или обтеканием криволинейных поверхностей шара, цилиндра и т. д. Влияние градиента скорости на структуру колеблющегося пограничного слоя определим методом последовательных прибли жений. В этом случае для анализа удобно внести функции тока
для пульсационных составляющих: |
|
|
|
||
Аы(,) II |
ДоП -- |
дг|?<Г) |
’ |
||
д Г |
|||||
|
|
|
|||
Ди<"> “ |
Д „ (" ) |
_ |
дх |
|
|
ду ’ |
-- |
|
|||
|
LAU' ' |
|
|
Решая систему уравнений (191) и (192) по методу последова тельных приближений, получим
ф' = бкАмооо (*) o' (л) ехр (Ш)\
(262)
ФП = бк |
[аП (л) ехр ( 2 Ш) + а р (л )], |
102
где t] = t//8K, а |
значения функций |
o(,) (л), |
(л), |
o p |
(rj) соот |
|||
ветственно |
|
|
|
|
|
|
|
|
an (tj) = - |
iL p il + |
n + |
-Ц-L exp [ - (1 + |
0 4l; |
|
|||
|
|
(л) = |
|
exp [ — 1/2(1 +0Л 1 + |
|
|
||
|
+ |
“ Г Ч |
е х р ( - |
( 1 + |
О ЧЛ — |
’• |
|
(263) |
of* (л) |
= |
— j |
ехр (— |
2л)— ехр (— |
л) cos л — |
|
||
— ехр (— л) sin л ---- у- ехр (— л) sin л — |
4 Ч + |
-j- |
• |
Как следует из приведенных соотношений (262) и (263), во вто ром приближении в колеблющемся пограничном слое появляются
колебания |
удвоенной частоты, |
и |
соответственно |
осредненные |
|||||||
составляющие скорости |
|
|
|
|
|
|
|
|
|
||
|
(Ды<">) |
< ? г |Л |
1 |
/ . |
|
д Д « о а > |
\ |
даР |
|
(264) |
|
|
ду |
(о |
\ |
ЫОо° |
|
дх |
/ |
дх\ |
’ |
||
|
|
|
|
||||||||
|
|
сК|>П |
|
|
6К |
1 |
д2(Аы0со) |
|
(265) |
||
(ДцП) = — ~ д Г |
|
|
со |
2 |
д *2 |
|
|
||||
где |
|
|
|
|
|
|
|
|
|
|
|
&>Р |
ехР (— 2л) + 4* |
— Ч) ехР (— Л)cos Л + |
|
||||||||
дг) |
4 |
|
|||||||||
|
+ |
4 ( 4 + |
Ч)ехр(— л)51пл — 4 - |
|
|
|
Осредненное по времени поле скоростей существенно зависит от градиента колебательной составляющей скорости внешнего
дАUQQQ
потока; при — = 0 (Ди<">) = 0. Поперечная скорость
(ДцО) достигает максимума в точках, где градиент |
!(д“о») |
|
дх2 |
||
|
максимален.
Зная осредненную по времени скорость в пограничном слое, можно определить и осредненное по времени пульсационное ка сательное напряжение на поверхности:
|
|
а(АцП)1 |
|
|
|
|
ду J < /= o |
|
|
_М_ |
Диош |
dA“jco |
(266) |
|
собк |
дх |
|||
|
|
103
где |
|
|
|
|
|
|
|
|
dVP |
|
i |
|
2т1) + |
|
|
||
~ д ^ Г = — 2 exP( |
|
|
||||||
+ (ri— 2) exp (— ri) sin T] + |
exp (— |
cos ц. |
(267) |
|||||
Отсюда находим |
|
|
|
|
|
|
|
|
|
( |
|
\ |
_ |
_1_ |
|
|
|
И |
\ *1* |
/ п=о~ |
2 |
|
|
|
||
|
|
|
|
|
|
|
|
|
<АтР> = |
^ i |
4 |
- \ Д«о» ^ |
г у |
= |
|
||
|
1 |
/ |
|
дДи |
\ |
|
|
|
= |
fiK- f |
P \ |
|
AU“ T |
|
/ - |
|
<268) |
Таким образом, во втором приближении касательное напря жение силы трения на поверхности пропорционально толщине колеблющегося слоя и градиенту колебательной составляющей скорости внешнего потока.
Рассмотрим плоскую пластину, обтекаемую высокочастотной стоячей волной. При этом будем считать, что стационарное тече ние отсутствует. В частном случае для гармонической стоячей волны
Апоа> = — A cos kx,
где k = -д----- волновое число.
Вэтом случае компоненты осредненных по времени-скоростей
впограничном слое
(Аы° ) = |
i - 4 - sin (2^ ) ^ Г |
(269) |
|
(ДуП) = |
М а A2cos (26л:)стр. |
||
|
|||
|
ю |
|
|
При 1] = у!6К<С 1, с точностью до вели |
|||
чин первого |
порядка малости |
<то ^ 0; |
|
— -g- г], и тогда получим |
|
||
(Ди<'>) = |
г]sin (2kx); (on ) = |
0. (270) |
Согласно выражению (270) скорость те чения (Ды(//)) вблизи стенки зависит от ки нематической вязкости («*v-0*5), причем ско-
Рис. 21. Распределение функции Сто (у/вк) по толщине колеблющегося пограничного слоя
104
рость равняется нулю в узлах и пучностях стоячей волны. Вид
функции а р ^у - ^ , определяющей скорость течения в колеблю
щемся пограничном слое стоячей волны, показан на рис. 21. На правление скорости течения вдоль стенки зависит от знака произ
водной ^ (264). Эта производная (рис. 21) при у <=* 6К меняет
знак на обратный, что приводит к образованию в пограничном слое стоячей волны вихрей, размеры которых имеют порядок —(Л/4) 28к.
Результаты, полученные для плоского пограничного слоя, при годны также и тогда, когда поверхность изогнута, если местный радиус кривизны поверхности много больше длины вязкой волны. Рассмотрим течение вокруг круглого цилиндра. При ReA > Мд, т. е. при толщине колеблющегося пограничного слоя, много мень шей по сравнению с длиной волны колебаний внешнего потока (8К С Л и kr0 < 1), для определения течения вблизи цилиндра можно воспользоваться решением (264) для плоского погранич ного слоя [33].
Для этой цели координату х направим вдоль поверхности ци линдра; тогда скорость потенциального внешнего потока
ДИо, (х) = 2Дыосо sin — = 2Д«осо sin <р,
Г0
где ф — угловая координата; г0 — радиус цилиндра.
Из выражений (269) следует, что компоненты скорости в по граничном слое во втором приближении соответственно
(Д «<">} |
ДЫО оо 1 |
доГ . |
|
|
a |
krQ Sin 2ф д (Уо/&к) ’ |
|
|
|
kb* |
(271) |
(Дц<">) |
|
cos 2ср а<' |
|
|
|
а |
|
здесь а — скорость звука; |
k = 2я/Л — волновое число. |
Характер течения вокруг цилиндра показан на рис. 22 для слу чая R 0/bK = 7. Как и в случае обтекания плоской пластины стоя чей волной, вихрь в пограничном слое вращается в направлении, противоположном вихрю вне пограничного слоя. Центр враще
ния внешнего вихря находится на |
бесконечности, |
а ско |
|
рость в вихре |
|
|
|
(ДиО <*, СС» = - А Аио„ % = = - |
-I- ^ |
s,„ ( £ ) . |
(272) |
Аналогичная картина возникает и в случае обтекания сферы колеблющимся потоком; в этом случае стационарное вторичное течение набегает на экватор сферы (рис. 23), а продольная компо
нента скорости вторичных течений вне пограничного слоя |
|
||
Attn |
Sin |
Ov |
(273) |
д иП(х, ОО) = - 1,4 |
, |
где R о — радиус сферы.
105
90° 600
|
|
|
|
|
|
Рис. 22. Картина линий тока |
||||
|
|
|
|
|
|
стационарного вторичного те |
||||
|
|
|
|
|
|
чения |
в |
окрестности ци |
||
|
|
|
|
|
|
линдра: |
|
|
|
|
|
|
|
|
|
|
а — колебание |
малой |
ампли |
||
|
|
|
|
|
|
туды; б |
— |
колебание большой |
||
|
|
|
|
|
|
амплитуды; |
в — |
малые |
значе- |
|
ati,n О |
2 |
4 |
6 |
8 Ю |
н |
ния числа |
R e. |
|
|
|
|
|
А |
|
|
||||||
|
|
|
В) |
|
|
|
|
|
|
|
Подобные вихревые вторичные течения возникают не только у неподвижных препятствий, помещенных в колеблющуюся среду, но также и около колеблющихся тел: пластин, цилиндров, шаров, газовых пузырей и т. д.
Таким образом, малые колебания среды или тела приводят к возникновению не только колеблющегося пограничного слоя, но и добавочного вторичного течения на большом расстоянии от тела, и притом такого, которое не зависит от вязкости среды.
Следует иметь в виду, что второе приближение пригодно для малых колебательных чисел Маха (Мд) и при больших значениях параметра kr0ReA, т. е. для малых амплитуд колебания скорости. При больших амплитудах колебания скорости и при условии, что течение остается ламинарным, харак тер обтекания тел изменяется. На рис. 22, б приведена картина обте кания круглого цилиндра с точно стью до величины четвертого по рядка малости для случая г0/бк = 7,
Мд/6г0 = 1 0 . |
В этом |
случае про- |
Рис. 23. Картина линий тока |
внешних тече- |
|
ний в окрестности |
сферы |
|
106
исходит деформация формы как внутреннего, так и внешнего вихрей [15].
При ReA ;£= Мд толщина колеблющегося пограничного слоя сравнима, или больше длины волны колебаний. В этом случае течение характеризуется тем, что параметр k\ R 0\ ReA мал и модель течения, основывающаяся на уравнениях пограничного слоя, ста новится недостаточно точной.
По мере уменьшения ReA размеры пограничного вихря ста новятся все большими, причем эти размеры могут значительно превышать толщину колеблющегося пограничного слоя 8К, и на конец, при очень малых ReA пограничное течение занимает все пространство. Такая картина течения представлена на рис. 22, в (соответствует случаю R 0/8K— 1,4 и &RoReA<10). Аналогичная вихревая картина движения жидкости может наблюдаться в ка нале при возбуждении колебаний стоячей волной.
Рассмотрим течение в плоском канале между двумя параллель ными плоскостями шириной 2г0, вызываемое стоячей волной (фронт волны перпендикулярен к плоскости). Распределение ско ростей в первом приближении определяется согласно выраже ниям (250) и (251) для гармонической стоячей волны уравнением
Аи0(0, х) = — Ам0 (0) cos kx\
распределение скоростей во втором приближении носит вихревой характер. В ядре потока вдали от стенки поле скоростей во втором приближении [15]
<д„°> = - 1 1 ^ 1 1 cos 2 b [ ( l — f - ) - ( l - i ) ! ] . j
Рис. 24. Картина течения в плоском канале, обусловленная стоячей волной:
а — схема вихревых течений в погра» ичном слое; б — схема вторичных вихревых тече ний в ядре потока
107
у = 0,423г0;
"Г " II оо| >
у - 1,577/v,
Рис. 25. Картина вторичных вихревых течений в цилиндрическом канале при поперечных резонансных акустических колебаниях
Характер течения согласно выра жению (274) показан на рис. 24. Течение вне пограничного слоя состоит из чередующихся вихрей с периодом Л/2 вдоль оси канала ох. Вращение вихрей происходит вокруг центров с координатами
|
|
|
|
Ю |
|
|
у = |
|
у = |
II |
00ЯI |
> |
|
0,423го; |
0,423г0; |
|
(275) |
|||
/ |
З А \ |
|
|
|
|
|
( ' — г ) 1' - - ; ( * = < V > A)‘ |
|
|||||
У = |
1,577г0; |
у = |
1,577г0. |
|
|
В отличие от течения в колеблющемся пограничном слое ско рость течения вне пограничного слоя не зависит от вязкости. Однако образование вихревого движения вне пограничного слоя обусловлено вязкостью среды. Вихри, образовавшиеся в колеб лющемся пограничном слое (рис. 24, б), возникают вследствие вязкости среды, а вихри вне пограничного слоя (рис. 24, а) воз никают в результате взаимодействия с вихрями в пограничном слое. Вращение вихрей в пограничном слое происходит в на правлении, противоположном вращению вихрей вне пограничного слоя. Аналогичная картина возникает в цилиндрическом канале. При возбуждении в канале поперечных резонансных акустических колебаний, направленных вдоль радиуса канала, возникают вто ричные вихревые течения, как и в случае продольных колебаний. Вращение вихрей осуществляется в плоскости поперечного сече ния канала (рис. 25). Методика расчета таких течений приведена
вработе [39].
4.ТЕПЛООБМЕН В ЛАМИНАРНЫХ КОЛЕБЛЮЩИХСЯ ПОТОКАХ
Анализ температурного поля в колеблющихся потоках показывает, что основное изменение колебаний температуры, как и в случае рассмотрения поля скоростей, наблюдается вблизи поверхности тела. Вблизи поверхности тела образуется тепловой колеблющийся пограничный слой. Из сравнения уравнения энер гии и уравнения движения толщина колеблющегося теплового 1Щг граничного слоя
108
\ Метод последовательных приближений можно использовать только в том случае, когда толщина теплового стационарного пограничного слоя много больше, чем толщина колеблющегося теплового слоя, т. е.
6к. т |
l / J Z - c * V |
°-08 |
1- |
||
6к |
У |
Шбт0 |
У |
Для анализа процесса теплообмена в колеблющихся погранич
ных слоях удобно |
уравнение энергии записать |
относительно |
функции тока ф и |
безразмерной температуры |
Т — т |
0 = =------. |
||
|
|
‘w 1 <о |
С учетом диссипации кинетической энергии для жидкости с по стоянными физическими свойствами уравнение энергии для по граничного слоя и граничные условия при постоянной температуре поверхности и внешнего потока, имеют вид
ае . |
аф |
ае |
аф |
ее |
е*е |
v |
/ д * у у |
dt + |
ду |
дх |
дх |
ду |
~ а д у * ^ |
Cp(Tw - T ^ ) |
\д у * ) » |
|
|
|
|
при |
у — 0 0 |
= 1 ; |
(276) |
|
|
|
|
|
|||
|
|
|
|
при у —►оо 0 |
= 0. |
|
Это уравнение относительно пульсационных и осредненных по времени величин можно записать, в виде:
для пульсационного движения
аде |
, |
афоаде _аф0аде _ |
е*де |
,______ |
v |
|
/а 2дф\2 |
|||||
dt |
|
ду дх |
дх |
ду |
~ й |
ду* |
cp (Tw |
~ Tо°) |
' дУ2 |
' |
||
|
|
д^ д % , д _ ^ д ь , _ д м ? ш |
адфаде |
, |
|
|||||||
|
|
ду дх |
дх ду |
ду дх |
' |
дх |
ду |
' |
|
|||
|
|
|
I_______у |
|
9 5*Фо а*дф. |
|
|
(277) |
||||
|
|
|
^ |
op (Tv - T |
a ) |
ду* |
ду*~' |
|
|
|||
|
|
|
|
|
|
|||||||
для |
осредненного |
движения |
|
|
|
( д2Ф0\ 2_ |
|
|||||
|
афо аво _ |
афо а§о __ |
д*% |
, |
|
у |
|
|
||||
|
ду |
ду |
дх ду |
|
а ду* |
^ |
ср (Tw — Т„) |
\ ду* ) |
|
|||
/адф аде\ |
/адф аде\ |
|
|
v |
|
/ / |
аадф \2\ |
|||||
\ д |
у |
д х / |
\ |
дх |
ду / |
' |
Ср (Tw - T |
/ |
\ V |
ду* ) |
/ + |
|
|
|
|
|
|
|
|
|
|
|
|
|
(278) |
Решение системы уравнений для сравнительно малых амплитуд колебания скорости внешнего потока будем искать в виде суммы ряда по отдельным гармоникам:
Д0 = 2 8<'>Л«> (л, 0; |
ч = у у |
^ \ |
е = t Sh. |
(279) |
1=1 |
Г |
VX |
|
|
109
Рассмотрим влияние колебаний скорости внешнего Потока с постоянной амплитудой колебаний на тепловой пограничный слой в предположении, что диссипацией кинетической энергии можно пренебречь. Это допущение может быть оправдано для сравнительно небольших амплитуд колебания скорости. Прене брегая в первом приближении влиянием нелинейных членов как в пульсационном, так и в осредненном по времени уравнениях энергии-и-лспользуя выражение (277), получим уравнения тепло вого пограничного слоя для степенного закона изменения скорости
«о = Ахт относительно безразмерных |
параметров: |
|
||||
для |
осредненного |
движения |
|
|
||
|
|
J_ |
т + 1 f д9о — п |
(280) |
||
|
|
Pr drf |
‘ |
2 '°дг\ |
|
|
для |
пульсационного движения |
|
|
|||
|
' 1 d2h , |
m +1 |
, dh |
— |
|
|
|
Pr 5r|2 |
2 |
dr] |
' |
' Зт) * dr| |
|
|
m +- 1 |
o0o |
|
d0o dh |
(281) |
|
|
—T~ g |
5r) |
— (1— m)C dn\ dt, |
Решение уравнений теплового пограничного слоя (279) можно, как и в случае анализа уравнения движения, представить в виде степенного ряда (215) и (216) для больших и малых значений чисел Струхаля. В частности, для малых значений чисел Sh
fibi, Sh) = £ с*л*(л); fe=0
тогда значение числа Нуссельта будет определяться степенным рядом
-77= - = — [9' (0) + е ехр (Ш) h'k(0)]. V Re
Для плоской пластины при т = 0 и Рг = 0,72 для малых зна
чений чисел Струхаля |
|
|
= 0,291 + е ехр (Ш) {0,148 + |
0,125 Sh2Н------ |
|
----- h (/0,02 Sh Н------ |
}. |
(282) |
Из этого выражения следует, что при малых значениях Частоты
иамплитуды колебания скорости внешнего потока, т. е. при е —* 0
иSh —*0, коэффициент теплоотдачи стремится к квазистационарному значению. При увеличении частоты наблюдается сдвиг по фазе между колебаниями скорости внешнего потока и коэф фициентом теплоотдачи. В общем случае коэффициент теплоотдачи является величиной комплексной, мнимая часть которого харак теризует фазу колебания.
ПО