
книги / Тепловые и гидродинамические процессы в колеблющихся потоках
..pdfНаименование граничных условий
Круглое отверстие диаметром d в экране, стоящем поперек трубы диаметром D (без учета сопротивления излучения и трения)
Ряд круглых отверстий площадью <л =
= -5^- в бесконечной перегородке (без
учета сопротивления трения) при di <
< 0,2 "j/ " (Fi — площадь на одно от
верстие)
Ряд (я) круглых отверстий диаметром
di > 0,7 (Fi — площадь одного от
верстия), которые в общем занимают кру-
говую зону площадью F = |
nd2 |
—г— |
|
! |
4 |
Схематическое |
Электрический |
|
изображение |
||
эквивалент |
||
граничных |
||
импеданса |
||
условий |
п
|
* |
3 ' |
[ |
|
Ф6 |
[ |
|
|
[ |
|
|
[ |
о [о ]о |
ООО
О о ) т
потверстии п ветвей
4:i< j 1
Продолжение табл. 1
Акустический импеданс
-F А р
Z s s А и
' ( т г ) “ ( 1 _ м , т г + а д ‘ 5 +
d6 |
\ —1 |
|
+ °. °7 |
Ч------) |
(формула Фока); |
при |
d - * D Е |
->• 0 |
Z & /coMj/г,
где
а
Й
[Наименование граничных условий
Импеданс на конце трубы F при при соединении в центре импеданса za , сосре-
nd2
доточенного на площади а — ——
Имцеданс входного отверстия (площадью F) в емкости объемом V, в которую вклю чен импеданс za (на площади о) при А >
> V v
Импеданс трубы (сечением F), развет вленной на ряд трубок площадью ai, а2, а3, . . . с импедансами zi, z2> z3
|
|
|
|
|
Продолжение табл. 1 |
|
Схематическое |
Электрический |
Акустический импеданс |
||||
изображение |
||||||
граничных |
эквивалент |
|
7 _ |
F A Р |
|
|
условий |
импеданса |
|
|
|
А и |
|
|
|
|
|
|
||
|
|
1 ** о5- (Z<T+ 1(0"2tT) _ |
1? (Z<T + 1<оМ): |
|||
|
|
при |
d < |
0,2-2 "j/" |
||
|
|
при |
d > |
0,7-2 |
|
|
|
|
|
М^О |
|
||
|
|
|
1 _ |
|
1 |
1 |
|
|
z ~ |
|
z^ |
|
|
h |
|
|
_ |
|
pa2F2 . |
|
|
|
Z v _ |
icoF |
’ |
||
|
|
|
F*_ |
_°i r |
_ai + _ai |
+ |
Z |
zx + |
z2 ^ z3 |
^ |
Решая систему уравнений (172) при граничных условиях (175), получим систему уравнений для определения констант Cll2:
Ci [оо (0) |
г|5, (0)— 2 (0)] Fl (0) + |
С2 [сю(0) ^ |
(0) - 2 (0)] F*2(0) = Ф0; |
С, [оо (L) t i (Ц — 2 (L)] Fi* (L) + |
С2 [fl0 (С) afc (L) — 2 (L)] F2 (L) = 0; |
||
здесь FIj2 |
— прямая и обратная волны, |
|
f b " T ^ w exp№ ' 2W)-
Решая систему относительно констант Ci и С2, получим
Ci |
- |
Ф 0 [оо (L) тр2 (L) — 2 (L)] F2 (L) ^ |
; |
С2 |
= |
Ф 0 [оо (L) o|)i (L) — 2 (L)] Fi (L) |
, |
где А — определитель системы, который можно выразить уравне* нием
А = [а0(0) t i (0) - 2 (0)] [а0 (L) |
(L) - |
|
|
— 2 (L)] Fi (0) F l (0) — |
[оо (0) ip2(0) — |
2 (0)1 х |
|
X [а0(L) % (L) — 2 (L)] Fl (L) Fl (0). |
(176) |
||
При этом следует иметь в виду, что |
|
|
|
К 2(0) = |
Vа0(0) ' |
|
|
Подставляя полученные значения констант Ci и С2, получим следующие выражения для определения комплексных амплитуд
колебания давления Ф (х) и |
массовой скорости F (х) |
по длине |
канала: |
|
|
Рф} =~i-U«o (L) |
(L) - 2 (L)] Fl (L) Fl (х) - |
|
— [а0(L) % (L)— 2 (L)] Fl (L) Fl (*)); |
|
|
-Ф 7 = -^{№o(L)b (L) - 2 (L)} Fl (L) х |
(177) |
|
X atf (х) Fl (х) — [оо (L) я|?1(L) — 2 (L)} х |
|
|
X Fl(L)aob(x) Fl(x)}. |
|
В качестве примера рассмотрим ряд частных случаев.
73
Предположим, что на входе заданы амплитуда колебания дав
ления Ф (0) = |
Ф0 и импеданс Z (0) = |
|
0; тогда |
|
|||
|
[«О (L) % (L) - Z (L)] F'2 (L) F\ (х) - |
[а0 (L) % (L) - |
|
||||
FJM = |
|
_________ -Z (L )\F \(L )F *2 {X) ______________. |
|
||||
Ф° |
Та0 |
(0) Ь (0)] [а0 (L) >|>2 (0) - |
Z (L)] |
(0) F*2 (L) - |
’ |
||
|
- |
[«о (0) % (°)] [«о (L) |
(L) - |
Z (L)] F\ (L) F\ (0) |
|
||
|
[«0 (D % (L) ~ z (Щ F\ (L) o fa |
(x) Fl (X) - |
|
||||
Ф(х) _ |
- [ a 0 ( L ) % ( L ) - Z ( L ) ] F \( L ) a 0% (x)Fl(x) |
|
|||||
фо |
[a0 (0) ^ (0)] [a0 (L) i|>2 (L) - |
Z (L)] F\ (0) F*2 (L) - |
' |
||||
|
~ |
K (0 ) % (0)] [a0 (L) ^ |
(L) - |
Z (L)] F\ (L) F*2 (0) |
|
В некоторых случаях при расчете конкретных задач приведенное выражение удобно представить в виде, когда в качестве мас
штаба выбрана амплитуда колебания давления на выходе из ка нала:
[a0( L ) % ( L ) - Z i i . )
F1 (L)
f ( ,, |
_ |
|
|
|
|
|
|
|
|
Ф{1) |
К |
(L) i|>2 (L) - |
Z (L)] a0 (L) ^ |
(L) - |
* |
||||
|
— К |
(L) 4>i (L) - |
Z (L)] a0 (L) i|>2 (L) |
||||||
|
|
|
|
|
|
|
|
(178a) |
|
|
|
|
|
|
|
|
f j |
lx) |
|
|
[a„ (L) % (L) - Z (L)] a0 (x) t|>x (x) |
- |
|||||||
|
|
|
|
|
|
|
Fi (L) |
||
|
- |
[<*o (L) Ь |
(L) - |
Z (L)] a0 (x) % (x) |
FUx) |
||||
Ф (x) |
F2 (L) |
||||||||
_ |
[a0 (L) |
(Z.) - |
Z (L)] a0 (L) |
|
|||||
Ф (L) |
|
(L) - |
|||||||
|
|
- |
(ao (L) 4>i (L) - |
Z (L)) a,, (L) % (L) |
При сравнительно небольших градиентах скорости звука 4р
идля больших значений Rea согласно выражениям (120) и
(121)имеем
k = т |
Г м„ ± Уд |
dx |
tl, 2 (*) |
м0 ± Уд . |
||
J |
! - М0 |
М*) ’ |
1 MQ ’ |
|||
|
||||||
|
|
|
21Л? |
|
1/4 |
|
t i W — |
= l-Mg(L) |
q -X,i{x)q^ ( L ) « l . |
74
Тогда выражения (17&а) для акустически закрытого канала [Z (0 = 0 ] примут вид
Ф (х) |
_ |
~ i f а (а) |
1 — Мц (L) |
[М0 (х) — V д (х) |
||
0(L) |
~ |
У а0 (L) |
2 V qjL ) |
|
{ 1 — М$(х) |
|
X exp [ki (X) — fcj (L)] |
М0 (*) + К у (х) |
|||||
|
l- M g ( x ) |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
(179) |
|
|
X exp [k2{x) — k2(L)\, |
||||
F(x) |
|
i — |
i - м а ю |
|||
Ф(Ц |
|
V a0 (x)a„(L) |
|
= |
7— {exp [kt (X) |
|
|
2 1Л? (L) |
|||||
|
|
— kt (L)] — exp [k2(x)- |
■кЛЩ- |
Действительная часть полученных выражений характеризует амплитуду колебаний, а мнимая — соответствующую фазу коле баний. В качестве примера рассмотрим случай малых значений коэффициента ослабления р. Можно принять
|
|
|
|
|
|
|
|
(180) |
|
|
|
|
Фь2 = ± |
т ± М0 |
„— 1/4 ~ 1. |
|
|
Подставляя приближенные значения kU2 и г|;1>а в |
выраже |
|||||||
ния |
(178а), |
при |
|
Z (0) = 0 получим |
|
|
||
|
|
|
|
|
|
. ехр [*! (х) — kx (L)] — |
|
|
|
F (л) |
_ |
1 ~ М0 (L)______— В exp [k2 (х) — k2(/.)] |
|
||||
|
|
V а0 (х) а0 (L) |
1 — М0 (L) + |
В (1 + М 0 (L)) |
(181) |
|||
|
|
_____ |
, |
|
(1 — М0 (х>) exp [kt (х) — kx (L)] + |
|||
Ф(х) |
I / |
а0 (х) |
1 — м о (О |
+ |
В (1 + Мр (х)) exp [k2(х) — kx (L)] |
|
||
Ф(Ю |
У |
а0(L) |
|
1 — Мл (х) |
|
1 - М 0 ( Ю + б ( 1 + M 0 (L)) |
|
|
здесь параметр |
|
|
|
|
|
|
||
|
|
|
|
, ______ ао (L) |
|
|||
|
|
|
|
R _ |
|
Z(L) ( 1 - М 0 (Ю) |
(182) |
|
|
|
|
|
1 | |
М Ц ______ |
|||
|
|
|
|
|
+Z (L) (1 — М„ (L))
Ограничимся рассмотрением случая, когда на выходе из канала имеется только активная составляющая импеданса Z (L), т. е^ Z (L) является величиной действительной.
75
Отделив действительную и мнимую части в выражении (181), полупим распределение амплитуд колебания давления и массовой скорости:
|
1 S X1 |
|
l-Mg(L) |
X |
|
||
|
|
Ф{1) |
|
V«о (*) ао(£•) |
|
||
|
|
jexp |
|
2со | |
р dx J |
+ |
|
|
|
|
|
|
|
10,5 |
|
+ В2ехр [ М |
Р Л ] - |
26 cos [ 2ш \1 Г ~ Мо О)) а0 (*) |
|
||||
X |
|
1 - M 0(L) + B (1 + M 0(D) |
|
||||
|
|
|
|||||
Ф(х) |
|
/ |
Щ, (х) |
1'- M 20(L) |
(183) |
||
Ф (L ) |
| |
f |
аa0 (L) |
\1.— |
д о X |
||
|
|
X |(1 |
— |
(* ))2 х |
|
|
|
X exp J^— 2(0 J |
p d * J + 5 * (1 + М0 (л:))2ехр ^2со J p < te j + |
|
|||||
|
|
|
|
|
х |
|0,5 |
|
X |
|
|
|
[ |
|
|
|
1 |
+ Mo |
(L) -|- В (1 + М0 (L)) |
|
||||
|
|
Для акустически закрытого канала Z (L) —» оо, поэтому В — 1 и выражение (183) преобразуется к виду, характерному для аку стически закрытого канала.
В заключение этого раздела рассмотрим колебания идеальной жидкости с постоянными физическими свойствами при отсутствии
осредненного движения |
= 0, ф112 ~ ± 1; q = klt 2 = + |
x'j . |
|||
Согласно формуле (176) распределение амплитуд колебания F \х) |
|||||
и Ф (х) будет иметь вид |
|
|
|
|
|
la0 + Z(L)] ехр |
( х— |
|
— Z ( L ) ] X |
|
|
|
Xexp[ ~ ^ (JC~ L)] |
|
|
||
[а0 — Z (0)] [а0 + |
Z (Z.)] exp |
i |
L J — |
|
|
— [flo + 2 (°)1 k |
- Z (L)] exp |
[^— |
L j |
|
76
|
|
|
[а0 + Z (L)] exp |
|
|
L)j — |
|
|||||
|
|
|
|
— [a0 — Z (L)] exp [ - ^ - ( x - Z-)] |
|
|
||||||
|
|
|
[a0 — Z (0)] [a0 + |
Z (£)] exp [ t |
L J — |
|
||||||
|
|
|
— [a„ + Z (0)] [a0— Z (L)Jexp |
|
L j |
|
||||||
Для акустически закрытого канала на выходе |
Z (L) —>оо и |
|||||||||||
открытого на |
входе Z (0) = 0 получим |
|
|
|
|
|
||||||
F(x) |
___i_ |
Sin[ j - (* ~ L)] . |
Ф(х) |
C°S [ ^ |
(X- |
L)] |
(184) |
|||||
Ф (0) |
a0 |
CS(0 iL) |
Ф(0) |
|
|
c o s ( | - L ) |
||||||
|
|
|
||||||||||
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
||||
Соответственно |
для |
акустически |
закрытого |
канала с обоих |
||||||||
концов |
|
|
|
|
|
|
|
|
|
|
|
|
|
Ф(х) |
|
|
Ф(х) |
C0S[~^ |
{X~ |
L)] |
|
|
|||
|
Ф(0) Z(0) = |
|
F\0) = ш * |
. |
/ |
CO |
\ |
|
|
|||
|
|
|
|
|
|
|
sin |
\ |
— L |
) |
|
|
|
|
|
|
|
|
|
|
«о |
/ |
|
|
|
|
|
F(x) |
2(0) = |
F(x) |
sin [ ^ |
(>~ ° ] |
|
(184a) |
||||
|
|
Ф(0) |
F(0) _ |
|
|
|
|
|
|
sln ( l7 0
Как следует из приведенных соотношений, при определенных значениях частоты cos
при Z(0) = 0 |
и Z(L) = oo |
c o s ( ^ - Z ,) = 0 ; |
|||
|
а0 |
л (2/г— 1) . |
|
|
|
s — |
L |
2 |
’ |
|
(185) |
при Z (0) —>оо |
и Z (L) —>оо |
|
sin |
||
|
Lj = 0; |
||||
f>s = |
х |
|
|
|
где /г = (1, 2, 3 . . .).
Следовательно, приведенное решение существует не при всех значениях частот. Если со = cos, то знаменатель в выражениях (184) обращается в нуль и решение задачи без начальных условий не существует, т. е. не существует установившегося режима [ампли туды колебаний Ар или А (р«) по времени неограниченно возра стают, так как не учитывается трение ]. Это соответствует резонанс ному режиму колебаний, при котором частота возмущений коле баний совпадает с частотой собственных колебаний газового столба в канале, величина п = 1, 2, 3 . . . определяет номер резонансной гармоники. Как следует из приведенных соотношений, собствен-
77
Щи) |
Рис. 9. Распределение амплитуды ко |
||||||
|
лебания давления и массовой ско |
||||||
|
рости по длине идеальной стоячей |
||||||
|
волны |
|
|
|
|
|
|
|
ная частота зависит от гра |
||||||
|
ничных |
условий. |
В |
качестве |
|||
|
примера |
рассмотрим |
случай |
||||
|
резонансных |
колебаний |
столба |
||||
|
жидкости в акустически закры |
||||||
|
том канале. |
В |
случае, |
если |
|||
|
частота |
колебаний |
со5 |
соответ |
|||
|
ствует первой резонансной гар |
||||||
|
монике |
(tl |
1) |
о)'s |
= |
- |
, на |
|
длине |
канала |
. укладывается |
||||
|
половина длины |
волны AI2=L |
(рис. 9), т. е. образуется стоячая волна. На входе и выходе в этом случае располагаются максимумы (пучности) колебаний давле ния Ар и минимумы (узлы) колебания А (ры), а в центральной части канала — максимумы (пучности) колебания Ар. Причем в стоячей волне идеальной жидкости колебания давления и мас совой скорости А (ры) сдвинуты по фазе на 90° как по времени, так и по длине канала. В момент времени, когда амплитуда колебания давления Ар в пучности давления максимальная, колебания мас совой скорости А (ры) по длине канала равнынулю, и, наоборот, когда колебания давления по длине канала равны нулю, ампли туда колебания массовой скорости в пучности скорости максималь ная. При этом каждый участок канала длиной Л/4, заключенный между узлом давления и ближайшим к нему узлом скорости, не обменивается энергией с соседними участками. Таким образом, в каждом таком участке происходит дважды за период колебаний превращение кинетической энергии колебаний, сосредоточенной преимущественно около пучности скорости (узла давления),
впотенциальную (энергию деформации), сосредоточенную преи мущественно около пучности давления (узла скорости).
При резонансных колебаниях идеальной жидкости распреде ление амплитуд давления и скорости по длине волны носит синусо идальный характер, а амплитуда их может расти беспредельно. В неизотермическом потоке вязкой жидкости установится конечное значение амплитуд колебания скорости и давления и исказится синусоидальный закон их распределения вследствие диссипации энергии по длине волны. При этом пучности скорости сместятся
всторону входного сечения канала. Поскольку возмущения дав ления малой амплитуды распространяются со скоростью звука, которая переменна в неизотермическом потоке, то это также при ведет к дополнительному смещению прочности скорости в сторону меньших значений температур газа.
Г Л А В А III
ТЕПЛОВЫЕ И ГИДРОДИНАМИЧЕСКИЕ ПРОЦЕССЫ В ЛАМИНАРНЫХ КОЛЕБЛЮЩИХСЯ ПОТОКАХ
1. МЕТОДЫ РАСЧЕТА ГИДРОДИНАМИКИ ДВУМЕРНЫХ КОЛЕБЛЮЩИХСЯ ПОТОКОВ
Существующие методы решения нестационарных урав нений гидродинамики в колеблющихся потоках основаны на упро щающих систему уравнений допущениях, достоверность которых требует экспериментальной проверки. Одним из наиболее распро страненных приближенных методов анализа нестационарной гидродинамики является метод последовательных приближений. Рассмотрим этот метод на примере плоского нестационарного по граничного слоя.
Нестационарное течение в плоском пограничном слое описы вается уравнениями неразрывности и движения
ар |
, д(ри) |
d(pv) |
_ п . |
а/ |
дх |
"Т" ду |
~ U’ |
где и, v — проекции скорости на оси лс и у, х, у — продольная
ипоперечная координаты.
Вобщем случае граничные условия для нестационарного те чения можно записать в следующем виде:
|
и = |
Щр(лс, |
t)\ |
v = vw (х, |
t) |
при |
у = 0; I |
|
|
и = |
и» (х , |
ty, |
v = Ооо |
(х , |
t) |
при |
(187) |
|
у — о о , / |
|||||||
здесь |
(лс, t) |
и vw (х, t) — есть |
продольная и поперечная ско |
рости движения поверхности тела, если оно движется; ит (х, t) —
скорость невязкого внешнего |
течения, |
связанная с давлением |
||||
в пограничном слое соотношением |
|
|
||||
1 |
др |
_ |
dam |
дит |
(188) |
|
р |
дх |
~ |
dt " г Ысо |
дх |
||
|
Соотношение (188) получается из уравнения (186), если в по следнем пренебречь членом, зависящим от трения, и поперечной
С К О Р О С Т Ь Ю Ооо.
Граничные условия в уравнениях (187) зависят от выбора си стемы координат. В дальнейшем будем рассматривать систему координат, жестко связанную с телом; следовательно, в такой си стеме координат % = 0 и % = 0. Однако в случае колеблющейся поверхности тела и стационарного внешнего течения предпочти-
79
тельнее использовать систему координат, в которой внешнее те чение является стационарным. Для несжимаемой среды различные системы координат равноценны.
При использовании метода последовательных приближений решение нестационарных уравнений пограничного слоя можно
записать в виде следующего ряда: |
|
|
|
|
||
и (х, у, t) = |
и0(х, у) + |
Ды(,) (х, |
у, 0 + |
' |
||
+ АиГ)(х, у, () + |
Аыс”)(х, |
у, |
t) Н------; |
|
||
V (х, у, t) = |
v0 (х, |
у) + |
До(,) (х, у, 0 + |
’ |
||
4-До(' ) (х, |
у, 0 + |
A v ^ ix , |
у, |
0 4 ------; |
|
здесь индекс ноль соответствует стационарному потоку жидкости, а штрихи означают порядок малости колеблющихся составляю щих скорости.
Рассмотрим случай несжимаемой жидкости (т. е. р = const)
спостоянными физическими свойствами. Подставляя выражения
(189)в уравнения (186) и пренебрегая в первом приближении не-
u |
и |
ди |
и |
v |
да |
получим: |
линеиными членами |
|
|
уравнения стационарного пограничного слоя в нулевом при
ближении |
|
dun . |
дил |
|
|
|
|
du„ |
|
|||
|
|
|
|
|
|
|
||||||
|
|
|
*/> |
|
|
CO |
|
|||||
|
|
“•Tto |
|
|
“ v |
|
|
|
(190) |
|||
|
|
|
du0 |
i |
|
|
= 0 - |
|
|
|||
|
|
|
dy |
|
|
|
|
|||||
|
|
|
dx |
+ |
|
U’ |
|
|
|
|||
уравнение пограничного слоя в первом приближении |
||||||||||||
|
|
|
дА «<'> |
|
. а2Д«<'> |
_ |
дАих ' |
|
|
|||
|
|
|
dt |
|
V |
ду* |
~ |
|
dt |
’ |
(191) |
|
|
|
|
дА «<'> |
, |
дА г><'> |
л |
|
|||||
|
|
|
|
дх |
+ |
|
ду |
~ |
U’ |
|
|
|
уравнения пограничного-слоя во втором приближении |
||||||||||||
ад «<"> |
v |
ааДи<"> |
. |
ЭАиа — Аи' |
ад ^ ^ |
— Av' |
дАиУ' |
|||||
dt |
ду2 |
~ A"°° |
|
дх |
|
|
|
дх |
|
ду |
||
|
|
|
дА и{,,) |
. д Av("} |
__ |
|
|
(192) |
||||
|
|
|
~ д х |
|
|
1 |
ду |
~ |
|
|
|
Граничные условия для уравнений в системе координат, жестко связанной с телом, примут вид:
при у — 0 и0 = Аи(П = Дыг) =••.== 0; ■
v0 = Av<') = Avr >= . . . = |
0; |
при у-* оо и0= «о.» (х); |
(193) |
Дц' = Ды0~ (х, 0» Ди' = Дм" = •••= 0.
80