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

Аэродинамические источники шума

..pdf
Скачиваний:
57
Добавлен:
15.11.2022
Размер:
11.44 Mб
Скачать

Величину т в дальнейшем будем называть фазой. Последнее выра­ жение можно получить непосредственно из уравнения (1.62), за­ писанного в естественных координатах

 

^

+ V { S ’ ф )1 7 = а

 

 

(1Л 35)

Первому

интегралу r\ = const для уравнения

(1.131)

соответствует

ф = const

в естественных

координатах для

уравнения

(1.135). Об­

щее решение уравнения

(1.135), т. е. решение, зависящее от произ­

вольной функции, в соответствии с (1.132)

и

(1.134),

можно пред­

ставить в виде

 

 

 

 

 

 

О= ^ а (0)(г, ф, г),

 

 

(1. 136)

где произвольная функция а(0) может быть определена, если возму­ щение энтропии задано как функция времени и координат на неко­ торой поверхности, пересекающей линии тока. Пусть это будет по­ верхность $= 5о(ф), тогда а(0)(т, ф, г) должна быть заданной функ­ цией.

Решение (1.136) имеет простой физический смысл: энтропийные возмущения переносятся вдоль линий тока с местной скоростью •основного течения без изменения величины. Другими словами, энтропия отдельных частиц жидкости в основном течении остается постоянной. В неоднородном течении скорость вдоль линии тока изменяется и, следовательно, изменяется расстояние между от­ дельными частицами жидкости, что приводит к изменению во вре­ мени пространственных производных энтропии в системе коорди­ нат, движущейся с местной скоростью отдельной частицы.

Первыми интегралами уравнений характеристик для системы (1.127) — (1.129) являются также T = const, *il) = const и z = const. В характеристических переменных система уравнений для возму­ щений завихренности (1.127) —(1.129) сводится к обыкновенным дифференциальным уравнениям, из которых возмущения определя­ ются интегрированием вдоль линии тока. Чтобы их проинтегриро­ вать, необходимо сначала вычислить производные о, входящие в правые части уравнений (1.127) —(1.129). Учитывая (1.136), а так­

же (1.134) для производной

ds

 

будем иметь

 

 

 

 

 

 

 

 

 

 

 

да

да^

дтз

 

dq(Q)

1

(1.

137)

ds = с р

дх

ds ■= — Сг дх

V(s, ф)

 

 

 

Несколько сложнее вычисляется производная до/дп:

 

 

 

да

^ Iда^

дф

, да^

дх \

(1

. 138)

дп

}\

дф

дп

дх IhT)

 

 

 

Для вычисления дх/дп вместо (1.134) удобнее использовать выра­ жение (1.133), помня, что нижний предел интегрирования в обоих

интегралах, вообще говоря, имеет зависимость $о(Ф) и Ы л)- В°с‘ Пользуемся также формулами (1.94) и (1.96)

дтз

1

 

г

 

?

(s'’0 ) d % ' = — -х

У (£о* ° )

dF0 j-

дп

 

 

 

 

dii

лт)(£. л)

<*]

J

^ 2(6'.

^(S o. ч) ^ 2(£о. т1>

 

1

 

to СП

 

 

 

,KW "

 

 

 

 

 

е

 

 

 

! _

« ^ - 5

> Л '.

 

 

+ 2 ------ —I

Лт)

1

л)

(1.139)

 

—г

<?1)

 

 

М ^ )

 

J

 

 

V

2 (e', Т))

 

 

 

 

£о (Ч)

 

 

 

 

 

 

 

Перейдем теперь обратно к переменной s, используя при этом со­ отношение (1.113)

-2т (s, ф) ^

ds'db/ds'

(1. 140)

т ($', ф)К(5 '„ ф>

$0 (Ф)

где m{s, <1»)=е0(«. ф)^(5, Ф).

С учетом (1.140) и (1.110) уравнение (1.138) можно переписать в; виде

да

до(°>

m(s,

ф) - 1

dsp

дп

дф

 

 

 

 

+

2 m(s, ф) J

ds'db/ds'

(1. 141)

m (s', ф Щ * ' , ф )

 

 

 

 

•So (Ф)

Поскольку далее предстоит интегрировать вдоль линии тока, выраженин для производных с удобно записать, выделив явную зави­ симость от переменной 5 . Учитывая (1.136), (1.137) и (1.141), за­ пишем

 

да __ / да \ Vo

(1. 142)

 

~ д 7 ~ [ ~ д 7 ) о ~

да

 

 

 

(1

. 143)

 

 

дп

(1. 144)

-дг£ - = (\- 7дг- )/ 0• где индексом нуль отмечены значения соответствующих величин в

точке s= s 0 (i|)), а также введен

интегральный параметр

0 , равный

6(5,

Ф) =

2m 0V 0

'i

ds'db/ds'

(1. 145)

т (s', <P)V(s', ф )

 

 

So (Ф)

 

 

 

 

 

На характеристике, т. е. на линии тока a^const, 2 =const при

x=const, уравнения

для

возмущений завихренности

(1.127) —

(1.129) сводятся к обыкновенным дифференциальным уравнениям вдоль линии тока.

Уравнение (1.128) на линии тока с учетом (1.144) записывается в виде

(1. 146)

Уравнение (.1.127), учитывая (1.145) и (1.146), можно представить следующим образом:

 

 

д

 

 

1

 

 

 

 

 

 

ds

 

Q0V

m0V0

 

 

 

, Ь147)

 

 

 

н

 

 

 

 

Аналогично для уравнения

(1.129) с учетом

(1.142), (1.143) и

(1.145)

имеем

 

 

 

 

 

 

 

д

Г

V2

 

 

 

 

 

 

ds

[ Qo

2т0ср

 

 

 

 

( L m >

где

 

^ (s ,

 

Ф)=

j (s',

<]>) — т0]

1

(1.149)

 

 

Qo(s', 4)

 

 

 

 

 

So (Ф)

 

 

ds

 

 

 

 

 

 

 

 

 

Из уравнений

 

(1.146) —(1.148)

имеем недостающие первые инте­

гралы системы

уравнений

(1.127) —(1.129)

 

 

 

 

 

 

 

 

 

 

(1.150)

(1.151)

я*

6 о

const.

(1.152)

Используя эти выражения, запишем общее решение всей системы уравнений. Если на поверхности з = з 0 ('Ф), пересекающей линии то­ ка, заданы функции

2 , = 2 <0>(/, ф, 2 );

2„=2<0)(*, ф, z);

(1.153)

Qz = Q i 0)(t, ф, 2),

тогда общее решение для возмущений завихренности с учетом (1.150) — (1.152) может быть записано в виде

(1. 154)

(1.155)

| 1 ^

da№

1

dv^

dsp j? (1/)1

(1.

156)

Vo

dr d i / /

V0

d v

dty

J ’

 

 

где функции Q«» и tf<°> берутся в момент времени т, определяемый уравнением (1.134).

Функции Q(°) и а<°), вообще говоря, не могут быть заданы пол­ ностью произвольно, поскольку решение должно удовлетворять ус­ ловию соленоидальности (1.130). Подставляя решение (1.154) — (1.156) в уравнение (1.130), получим

< Н у й = - —— L

т о К дх

<?2 i0)

d s 0

дх е о ( ^ о ) dty

Q o (K 0)

Vi—V\

Vp

 

 

Ко

•Я7 +

 

VRp a ? 4

V

L

 

 

 

mVо

dQ<0)

 

v20

 

 

m^V

d z

V R 0 d z

 

 

-m C

d V

I

dSp

<?2a(°> Q

(1.157)

J Qo (K )

J d<\>

dzdx

 

 

Сформулируем теперь достаточные условия, которые надо нало­ жить на функции £2<°) и о(0), чтобы уравнениен (1.157) удовлетворя­ лось во всем потоке. На поверхности 5=5о(ф) это уравнение мож­ но переписать в виде

1 «<°> _

о£»

<

0)

т0

е0 (^о)

dsp

М г0( )

Ко

(jg<°>

dt

Rp

 

дф

 

 

 

d z

R Q

d z

(1. 158)

С учетом уравнения (1.158), в котором необходимо заменить t на т, уравнение (1.157) можно представить следующим образом:.

d i v a = ( i - i . ) X

! L ( ^ .

* 1W

dsp

Ко

do<°>

 

V

mpj

V \ Rp

dx

-9o(Vo) di/

Rp

dz

 

 

Qo (Ко)

v i —vl

v

d V )

dsp

Э2о<0) .

 

 

 

■m \

0.

(1.159)

 

 

 

Qo(V')

dty

dxdz

 

 

 

V0

 

 

 

 

При отличных от нуля Q л0) о<0) и их производных это уравнение бу­ дет удовлетворяться во всем потоке, если

Ro

0 .

(1. 160)

Эти условия заведомо удовлетворяются в удаленной от обтека­ емого тела области выше по потоку, где течение однородно. Таким образом, если начальные возмущения Я<°>, которые задаются в об­ ласти однородности потока на удаленной плоскости, перпендику­ лярной к линиям тока, удовлетворяют условию

1

dQ<0)

dQ{n0)

dQ<0)

(1. 161)

- ------оН----------------------—

V0

dt

дф

dz

 

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

С учетом только что сказанного последние члены в уравнении (1.156) должны, очевидно, быть опущены. После этого оно приоб­ ретет вид

_ QoflO

) ( 0)

V2

да"»

~Qo{V)

V 2 -V \

 

2Vl

(1. 162)

Qo (Vo)

 

dv

 

2

Кроме того, параметры, отмеченные индексом нуль, исключая р0* приобретают смысл параметров в набегающем потоке.

Так же как и возмущения энтропии, малые возмущения завих­ ренности Q переносятся вдоль линии тока с местной скоростью ос­ новного потока, однако в отличие от возмущений энтропии их вели­ чина не остается постоянной для отдельных частиц жидкости ос­ новного потока. Изменение вдоль линии тока компонент завихрен­ ности определяется изменением параметров основного потока: ло­ кальным параметром — скоростью и интегральным параметром 0 . Величина завихренности £2П в соответствии с уравнением (1.155) определяется только локальным параметром V, поэтому завихрен­ ность жидкой частицы будет иметь равные значения в точках линии тока с равными скоростями в моменты времени, соответству­ ющие прохождению жидкой частицы этих точек. Этого нельзя ска­ зать о компонентах Qs и S2Z, поскольку они определяются интеграль­ ным параметром 0 . Завихренность в направлении линии тока Qs линейно зависит от завихренности Qn с коэффициентом пропорцио­ нальности, равным интегральному параметру 0 .

Важной особенностью решения (1.154) — (1.156) является так­ же то, что в неоднородном потоке возмущения завихренности зави­ сят от возмущений энтропии. Если даже начальные возмущения за­ вихренности отсутствуют, т. е. й<°)=0 , то это не означает, что за­ вихренность будет отсутствовать во всем потоке, напротив, она будет индуцироваться начальными возмущениями энтропии. Все эти особенности решения для возмущений завихренности не проти­ воречат известной теореме Кельвина о сохранении циркуляции по замкнутому жидкому контуру.

Обратим внимание еще на одно обстоятельство, связанное с ли­ нейностью системы уравнений (1.60), (1.61), (1.62) для возму­ щений завихренности и энтропии. Произвольное решение этой сис­ темы уравнений может быть представлено в виде линейной супер­ позиции двух решений, одно из которых определяется начальными возмущениями завихренности Q<°> при s ^ = 0 , а другое — возмуще­ ниями энтропии при й<0)^=0 .

1.6. ГЕНЕРАЦИЯ ЗВУКА ТУРБУЛЕНТНЫМ ПОТОКОМ ПРИ МАЛЫХ ЧИСЛАХ МАХА.

ТЕОРИЯ ЛАЙТХИЛЛА И ПАУЭЛЛА

Применим полученные в разд. 1.3 уравнения к случаю аэроди­ намической генерации звука турбулентным потоком при малых числах М2«с1. Поток предполагается неограниченным в том смыс­ ле, что в нем отсутствуют твердые или податливые границы, либо

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

Турбулентный поток при малых числах М является акустически компактным источником, т. е. масштаб L много меньше характер­ ной длины волны акустического излучения. Действительно, из со­ ображений подобия следует, что число Струхаля, образованное по

частоте f, соответствующей

основному энергонесущему масштабу

L и характерной

скорости

потока У, Sh = fL/V является, вообще

говоря, функцией числа Рейнольдса Re=VL/v и числа М =V/c

 

Sh = /(R e,M ).

В нашем случае,

поскольку

рассматривается турбулентный поток,

числа Re весьма велики, а число М мало, так что можно допустить, что число S h = / ( o o , 0) является постоянной величиной порядка единицы для всех чисел Re, больших некоторого критического. В экспериментах со струями это число равно 0,3, причем под ха­

рактерной скоростью потока понимается

скорость на срезе сопла,

а частота соответствует максимуму спектра излучаемой

звуковой

мощности. Таким образом, частота

имеет порядок f~V/L, а соот-

ветствующая длина волны Х

с

L*

много больше

характер­

 

^ > 1

ного размера источника.

Поместим начало координат где-нибудь в области турбулентно­ го потока, в случае струи, скажем, — на срезе сопла. При выполне­ нии неравенств M2 <Cl, L<^X в области г<СА, поток можно рассмат­ ривать в приближении несжимаемой жидкости, где г — расстояние от начала координат. Волновое уравнение Блохинцева — Хоу (1.48) в соответствии со сделанными ранее предположениями мо­ жно значительно упростить, пренебрегая малыми членами.

Сделаем необходимые порядковые оценки для величин, входя­ щих в уравнение (1.48). Рассмотрим сначала область г< 0 , где за­ ведомо отличны от нуля завихренность и неоднородности энтропии. Для производных по координатам и времени в этой области имеем djaxi = 0 (L-1), d/dt = 0 (f)= 0 (VL~X). Следовательно, для матери­ альной производной имеем оценку D/Dt = 0 (V *L~X) , отсюда также получаем Dv/Dt=0(V2L~l) .

Если поле скорости v представить в виде суммы

соленоидаль-

ной и безвихревой составляющей

 

<o=vs-\-vf1

(1. 163)

гд е d iv t f * = 0 ; r o t f l ^ O ,

 

46

то в области

безвихревой составляющей скорости можно пре­

небречь. Действительно, из ургЦшения непрерывности

 

 

J -

-2 £ \{ -d iv « '= 0

(1.164)

 

QC2

Dt

 

с учетом сделанных ранее оценок, а также принимая во внимание, что P ~ QV2, легко получить следующую оценку vr 0 (V М2). Таким образом, если пренебречь всеми членами, пропорциональными М2 в обеих частях уравнения (1.48), то на расстояниях г<^Х для эн­ тальпии торможения получим уравнение Пуассона

—v2/? = dlv{[fl X Vs] — TyS},

(1. 165)

при этом также предполагается справедливым уравнение адиабатичности DSJDt = 0.

Рассмотрим теперь область r^>Ly в которой по предположению

отсутствуют завихренность и неоднородности энтропии.

Следова­

тельно, правая часть уравнения

(1.48) в этой области обращается

в нуль. Для производных по координатам на расстояниях

r » L

справедливо d/dxi = О (Аг1), а

для производной по времени оста­

ется справедливой прежняя оценка. При записи левой части

(1.48)

можно различать два случая:

1) на расстояниях r^>L

среда по­

коится, как это имеет место в случае затопленной струи; 2 ) имеет­ ся внешний, вообще говоря, неоднородный потенциальный поток (струя в спутном потоке). Вторым членом в левой части уравнения (1.48) в большинстве случаев можно пренебречь. В первом случае от него остается только часть, пропорциональная пульсационной скорости, которая относится к членам, описывающим нелинейное распространение звука. Для потоков с малыми числами М нели­ нейные эффекты, связанные с распространением звука, который порождается самим потоком, пренебрежимо малы. Во втором слу­ чае, когда имеется внешний поток, этот член имеет порядок

1

Dv

п

/

/Ц2М2

\

. 166)

-------- v =

o

I

АХ

(1

с2

Dt

 

)

 

где Л — масштаб неоднородности внешнего потока, поле скорости в котором будем обозначать через Vo; т — параметр, равный отно­ шению скорости внешнего потока V0 к характерной скорости тур­ булентного потока V. В случае струи в спутном потоке этот пара­ метр называют параметром спутности. В обоих случаях, о которых

говорилось ранее, материальная производная

имеет одинаковый

порядок

 

— — = 0(Х-!).

(1.167)

сDt

Однако в первом случае эта производная

записывается в виде

D/Dt = d/dt, а во втором случае следует

учитывать скорость У0

внешнего потока D/Dt=d/dt-\-Vо/ ^ “ 7

 

Таким образом, в области r 3 >L с учетом сделанных оценок эн­ тальпия торможения В удовлетворяет в первом случае уравнению Даламбера

(iS^vJ)B=0' <1л«»

где с0 — скорость звука в покоящейся среде, а во втором случае — конвективному волновому уравнению

< 1 Л 6 9 )

Легко видеть, что уравнение для величины В, справедливое равно­ мерно во всем пространстве, можно записать в первом случае в виде

^ - V 2) 5 = div{[Qx ^ ] - 7 ’vS},

(1.170)

а во втором случае

+

B=d,v

(1.171)

В самом деле, в области r^>L правая часть обоих уравнений (1.170) и (1.171) обращается в нуль, т. е. справедливы уравнения (1.168) и (1.169), а в области гС ft, с уравнениями (1.170) и (1.171) можно повторить проведенные ранее рассуждения, с помощью ко­ торых можно убедиться в справедливости уравнения (1.165). Заме­

тим также, что в промежуточной области L<r<^K=— , энтальпия

М

торможения удовлетворяет уравнению Лапласа

ДД= 0.

(1. 172)

Кстати, существование самой этой промежуточной области свя­ зано с малостью числа М. На расстояниях, больших по сравнению с L, но малых по сравнению с ft, (L<^r<^%) с учетом убывания В с расстоянием (точнее той ее части, которая зависит от времени), можно записать общий вид решения уравнения (1.172)

в=в0+

г

OXi \ Г I OXi OXj \ г ]

(1. 173)

 

 

A ( a f f l ) +

- £ - ( S i £ W ) + ...,

 

где величины р, Diy

Qijy ...от координат

не зависят, В0— постоян­

ное значение энтальпии торможения вне турбулентного потока. Ограничимся далее в этом разделе рассмотрением случая, ког­

да, в области r^$>L среду можно рассматривать как покоящуюся, которая, как показано ранее, описывается уравнением (1.170). Будем также для простоты предполагать, что неоднородностями энтропии можно пренебречь. В частности, это предположение явля­ ется вполне реальным в случае изотермической затопленной струи.

Иными словами, далее будем рассматривать излучение звука за вихренностью, которое описывается уравнением

( 4 г ~ - V2) fi= d iv [Q X v*] = q.

(1. 174)

Это уравнение может быть преобразовано к виду, который впервые был получен в работе Лайтхилла [80] на основе соображе­ ний, приведенных в разд. 1.1. Правую часть уравнения (1.174) с учетом формулы (1.141) для соленоидальной составляющей скоро­ сти можно переписать следующим образом:

<7 = d iv [Q х г>]= — ^ ^

----—v2^2>

(1. 175)

dxidxj

2

 

при этом индекс 5, обозначающий соленоидальную составляющую скорости, далее будем только подразумевать. Учитывая выражение (1.175) для источника, представим его следующим образом:

(1.176)

dxidxj

где Pij = vivj —

Решение уравнения (1.174) для безграничного пространства, удов­ летворяющее условию излучения, может быть записано в виде

В (х, t ) - B 0 =

B'(x, t) = - ± - \d * y

' ’ D с

 

J

R

_ - L f

dyidyj

R

(1. 177)

4 я J

 

где B0— постоянное значение энтальпии торможения в отсутствие завихренности; 7? = |де—у |. Интегрирование по У (1.177) эффектив­ но сводится к области, в которой отлична от нуля завихренность, т. е. \м\ < L. Важное свойство интеграла (1.177) состоит в том, что производные по координатам можно вынести за знак интеграла

’ d‘yd' =

На расстояниях г=\х\, удовлетворяющих условию r^>L, выраже­ ние для R можно приближенно представить в виде

R = \ x - y \ ~ r — ^iji—\-... = r-\-0(L),

(1.179)

т. е. интеграл (1.178) можно переписать

•у ( Л « - г г ^ { т ! М * ' - £ Н (г>>ц- (,Л80>

Характерное время Т, за которое заметно изменяется величина Pij(y> *), равно T = 0(L/V), поэтому выражение (1.180) можно за­ писать также в виде

(г»д (,Л81>

пренебрегая различием в запаздывании во времени от различных точек области \у\< L, в которой отлична от нуля функция Рц. Оно не превышает величины порядка L/cQи поэтому при М<с1 мало па

сравнению с Т. При

формулу (1.181) можно представить

в форме (1.173), поскольку

при — <^ — = — = Т

вообще мож-

но пренебречь запаздыванием

с0

С0

V

 

 

 

 

 

5'(д:^) = - 4

- (

)

(£«г«Х ),

(1-182)

OXiOXj

\ г

]

 

 

где

 

 

 

 

 

Qi}{t) = - ^

\ p u(y ’ t)d3y-

(Ы 83)

Из сравнения с формулой (1.173) можно заключить, что в случае генерации звука завихренностью функции ц(/) и Di(\t) тождествен­ но равны нулю. С помощью рассуждений, аналогичных сделанным, в разд. И и 73 из [21], можно показать, что функция ц(^) отличня от нуля, если в нестационарном завихренном потоке имеются по­ датливые тела, объем которых может изменяться, a Di(t) отлична от нуля, если в потоке имеются твердые тела, не изменяющие сво­ его объема. В свободном турбулентном потоке разложение (1.173)> называемое мультипольным разложением, начинается с Qzj.

На расстояниях

пренебрегая членами порядка

в формуле (1.181)

г~х можно вынести за знак дифференцирования*

а производные по Х{ вычислить согласно следующей формуле:

д _

\

дг

д

__

Xj

J_

d

dxi

CQ

dxt

dt

r

Co

dt

В результате будем иметь

B '( x ,t) =

1 -

*

Q>A t

- - ) ■

_ 1

 

Г

dxidxj

 

Co /

X i X j

<Э2 Г

Ри (у,{ — r

OC

4лсд

гЗ

dfi j

\

c0

 

(1. 184)

»Х). (1. 185)

Решение общего вида, соответствующее формуле (1.173) на рассто-

50