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

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

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

Наименование граничных условий

Круглое отверстие диаметром 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

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