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

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

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

яниях.г^Х, чтобы оно удовлетворяло волновому уравнению

.(1.168), должно быть записано следующим образом:

 

В' = 1

(Ы86)

г

 

Излучение, соответствующее первому члену в (1.186), называется

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

формулой

(1.185)

такой нестационарный завихренный поток излу­

чает звук

квадрупольно.

 

 

 

Рассмотрим далее подробно решение (1.185) в волновой зоне

г^>Х. Запишем его в соответствии с определением

(1.176) для ве­

личины Pij

 

 

J

 

4 (1Л87>

В' (*, t)

1

[XjXj

азУ

4пс\

г3

 

 

 

 

 

где квадратные скобки означают, что величина, Заключенная в них, берется в момент времени т =4—г/с0. Второй интеграл в (1.187) тождественно равен нулю. Действительно, вспоминая, что интегри­ рование в (1.187) сводится к области |y|^L<cA,, где движение можно рассматривать в приближении несжимаемости, оценим вто­ рой интеграл. Запишем уравнение движения (1.42), пренебрегая сжимаемостью и неоднородностями энтропии,

 

i £ - + v5 =

- [ 0 x«].

(1.188)

 

ot

 

 

где В — энтальпия торможения

в приближении

несжимаемости

B = P!Qо + — •

 

 

 

Умножая скалярно

обе части

уравнения (1.188) на v будем

иметь

 

 

 

-1-

—~ = —v\/B= —div {vB).

(1. 189)

Воспользовавшись (1.189), вычислим интеграл

 

 

 

|й 1 у [г»5]й?3г/.

(1.190)

Последний интеграл преобразуется в интеграл по бесконечно удаленной поверхности, где выражение в скобках стремится к ну­ лю. Если сходится интеграл

 

 

^ Ь> X й] d.3y<iоо,

то

можно показать, что скорость tf(Q = rot_fl) при \у\-+оо стремит­

ся

к нулю как 0 (|у |~ 3). В то же время В на бесконечности стре­

мится

к постоянной. Отсюда следует, что интеграл (1.190) равен

нулю.

Другими словами, доказано, что для вычисления энтальпии

торможения В'(х, t) на расстояниях г^>Х нужно пользоваться сле­ дующим выражением для тензора Pij\

(1. 191)

Для определения звукового давления р' воспользуемся форму­ лой (1.64), которую на расстояниях г^>Х, где отсутствует средний поток, можно записать в виде

др'

дВ'

(1.192)

dt ~ 6°

dt

 

где Qo — невозмущенная плотность

покоящейся среды. Из

(1.192)

следует, что

 

 

p' = Q0B‘.

(1. 193)

Сопоставляя уравнения (1.174), (1.176), (1.191) и (1.193), при­ ходим к следующему результату: с точки зрения определения зву­ кового давления на расстояниях г^$>Х задачи, связанные с решени­ ем неоднородного волнового уравнения:

(1. 194)

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

_

&Ttj

(1 . 195)

^SB

dxtdXj

(1. 196)

4^ = e 0 div[Q X v s],

где тензор 7\j определяется согласно Лайтхиллу [5Д

 

Tij = Q0vsivsj.

(1. 197)

Уравнение (1.194) с источником вида (1.195), (1.197) является

основным уравнением в теории Лайтхилла, а с

источником

(1.196) — в теории Пауэлла [9Ц. Тем самым доказано, что в волно­ вой зоне эти теории являются эквивалентными в смысле определе­ ния звукового поля, порождаемого завихренностью при малых чис­ лах М турбулентного потока. Обе эти теории связывают характе­ ристики звукового поля с аэродинамическими параметрами неста­ ционарного завихренного потока. Например, зная основные пара­ метры турбулентности изотермической затопленной струи, можнополностью определить создаваемое ей звуковое поле, решая клас­ сическое волновое уравнение Даламбера с правой частью, опреде­ ляемой по Лайтхиллу, либо Пауэллу. Необходимо отметить, что в соответствии с формулой (1.197) источник звукового поля выража­ ется через тензор турбулентных напряжений Рейнольдса. Сущест­ венно, что сам источник определяется уравнениями движения не­ сжимаемой жидкости, его принципиально можно определить неза­ висимо от порождаемого им звука.

Следует также подчеркнуть, что уравнения (1.194) —(1.197) позволяют определить только главный член разложения по числу-

М для звукового поля, когда М2 <с1. В принципе могут быть най^ дены и высшие члены такого разложения с помощью метода сра­ щиваемых асимптотических разложений [49]. Рассмотрение излуче­ ния звука турбулентным потоком при малых числах М, проведен­ ное здесь, представляет собой один из примеров рассуждений, про­ водимых более формально в методе сращиваемых асимптотичео ких разложений.

В случае генерации звука завихренным потоком при малых числах М мы убедились, что источник может быть представлен в-, двух различных формах согласно Лайтхиллу и Пауэллу, причем, было показано, что в области генерации, где отлична от нуля за­ вихренность, эти источники существенно различны. Однако созда­ ваемые ими звуковые поля в волновой зоне совпадают в главных, членах разложения по числу М. Различие источников в области ге­ нерации определяется величиной QQAU2/2, где под V понимается, соленоидальная составляющая скорости. Эта неоднозначность ис­ точника наводит на мысль, что не все в источнике существенно с точки зрения генерации звука. Источник однозначно определяет звуковое поле в волновой зоне, в то же время обратное утвержде­ ние — не верно. В этой же связи можно поставить вопрос, каков, должен быть источник, чтобы он вообще не излучал звук. Ответ на этот вопрос позволяет сформулировать понятие эквивалентных по звуковому полю источников. Очевидно, что эквивалентные источни­ ки различаются на источник, не излучающий звук.

Перейдем теперь к более точному рассмотрению этих вопросов. Для простоты рассмотрим случай, когда справедливо классическое волновое уравнение (1.194). Пусть источник q сосредоточен в трех­ мерной области Д характерный линейный размер которой равен L. Представим функцию q(x, t) в правой части волнового уравнения (1.194), а также возмущение давления р'(х, t) в виде спектрально­ го разложения Фурье

+ оо

4 - оо

 

 

 

(2я)4 I

rftoe_/a)/ [ d3kelkxq (й, ID);

+ 00

 

+ 00

(1 .198>

1

 

p ' (л:,/)= -^- \ du>e-‘alp(x,u>)

^

[ d3kelkxp (k , u>)

2it J

(2Jt)4

 

 

 

 

 

(1 .199>

и используем это представление в формуле (1.177), дающей реше­ ние волнового уравнения, тогда будем иметь

p'(x,<o)——

\ d 3y

е‘

C°R

(1 . 2 0 0 )

 

 

4 JT

J

/?

 

 

при r^>L с учетом (1.179)

формула

(1.200)

приобретает

вид

. ш

 

 

 

 

 

 

I — Г

—I —пу ~

 

 

С о

 

 

",=-ii r S * *

е

с°

 

4яг

СО

( 1. 201)

 

q(y, ш)=

— Л, to

тде tti = ----единичный вектор в направлении излучения. Таким

образом, на расстояниях r ^ L звуковое поле определяется только теми спектральными компонентами источника, для которых волно­

вой вектор k лежит на сфере \к\ = —

, называемой излучающей

со

-

сферой. Другие волновые числа никакого вклада в излучение не дают, хотя они определяют пространственное распределение источ­ ника q. Отсюда, в частности, следует, что из измерений звукового поля на расстояниях r^>L нельзя сделать практически никаких суждений о пространственном распределении источника. Только

двухмерное многообразие |Л| = — в трехмерном пространстве водно­ го

вых чисел существенно с точки зрения излучения звука. Два источ­ ника, совпадающих на излучаюшей сфере, но совершенно различ­ ных при других волновых .числах, будут давать одинаковое звуко­ вое поле при r^>L. Излучения вообще не будет, если спектральная

плотность источника на излучающей сфере равна нулю q

/г, со^=

= 0 .Примером такого источника может служить источник вида

 

оо

 

q (x ,t)=

| \ - У 2)т<?гп(хЛ

(1 -2 0 2 )

 

т=*1

 

где все функции qm(х, t) быстро стремятся к нулю вне области D

(заведомо быстрее, чем убывает звуковое поле).

Спектральная

плотность, соответствующая источнику (1.2 0 2 ), имеет вид

оо

 

q(k, <о) = ^ ] (A2- ' J ) ’* £»(*.“).

(!• 203)

m=*l ^

•а соответствующее ему решение волнового уравнения (1.194) мож­ но записать

оо

(1.204)

т. е. оно весьма быстро стремится к нулю, не давая никакого вкла­ да в поток звуковой энергии через сферу достаточно большого ра­ диуса.

Предположим теперь, что размеры источника, т. е. области D, в которой отлична от нуля функция q, много меньше характерной длины волны X. Производной по времени в формуле (1.202) тогда можно пренебречь по сравнению с оператором Лапласа, так как <c0ldJdt=O(X- 1), a V = 0 (L _1) и формула (1.202) приобретает вид

<7(-М) = 2

t).

(1.205)

m= 1

 

 

Легко видеть, что источники в теориях Лайтхилла и Пауэлла различаются как раз на величину, подобную (1.205), т. е. на неиз­

лучающий источник. В этом случае ^ = —- - у - < 7т = 0 (т > \).

Функция q выражается через гидродинамические поля, связанные с несжимаемой жидкостью, и потому быстро убывает при rS>L.

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

Формула (1.201) позволяет также сделать некоторые суждения об интенсивности аэродинамических источников звука таких, как турбулентность при малых числах М. Действительно, если ограни­ читься линейным приближением для источника подобно тому, как: это делалось в разд. 1.4, то спектр волновых чисел источника, бу­ дет весьма похожим на спектр турбулентности. Максимум спектра’ турбулентности, очевидно, относится к волновым числам порядка k = (i)/Vc, где Vc = 0(V) — скорость конвекций вихрей и энтропийных неоднородностей; V — характерная скорость потока. Другими сло­ вами, при малых числах М излучающая сфера k=<d/c0 лежит да­

леко от максимума функции q(k , со), т. е. только малая часть энер­ гии турбулентного потока может превращаться в звук.

1.7. ИСТОЧНИКИ ЗВУКА В УРАВНЕНИИ БЛОХИНЦЕВА — ХОУ

Вернемся к рассмотрению линейных уравнений для вихревых,, энтропийных и акустических возмущений, распространяющихся в: безвихревом основном потоке с постоянной энтропией. Эти уравне­ ния были получены в разд. 1 .4 , а в разд. 1.5 найдены сами возму­ щения, относящиеся к вихревой и энтропийной компонентам дви­ жения, и исследовано их распространение в плоском двухмерном’ основном потоке. Общее решение, приведенное в разд. 1.5, позволя­ ет рассмотреть акустическую часть общей задачи о малых трех­ мерных возмущениях. В отличие от вихревых и энтропийных воз­ мущений, где решение удается получить без дальнейшей конкре­ тизации основного потока и вида самих возмущений, в случае акус­ тических Возмущений в общем виде можно только математически сформулировать задачу, т. е. определить все величины, входящие в уравнение Блохинцева — Хоу (1.63). Левая часть этого уравне­ ния, т. е. Волновой оператор, полностью определена, поскольку ос­ новной поток Но предположению известен. Остается определить правую часть, т. е. источник звука, связанный с вихревыми и энтро­ пийными возмущениями, чему посвящен данный раздел.

Правую часть уравнения (1.63) можно представить в виде

 

<7= div А -----!^-(Лу) V2,

 

 

4

2с2 v v'

 

где

A = [W x K

] - r 0va.

(1.206)

Следует заметить, что отсутствие начальных возмущений за­ вихренности и энтропии ft(°)= 0 и а<0) = 0 приводит вообще к отсут­ ствию этих возмущений во всем потоке, т. е. общее движение с уче­ том возмущений будет безвихревым и изэнтропийным й = 0 и а = 0 . Источник акустического излучения q в этом случае тождественно равен нулю, а уравнение (1.63) переходит в однородное уравнение Блохинцева. Этот случай нетривиален только с точки зрения рас­ пространения звука в основном потоке, возникновение которого не связано с самим потоком. Эта задача имеет самостоятельный ин­ терес.

В связи с тем, что общее решение задачи о возмущениях завих­ ренности и энтропии может быть представлено в виде суммы ре­ шений, из которых одно соответствует нулевым начальным возму­ щениям энтропии, а другое — нулевым начальным возмущениям

.завихренности, источник q также можно представить в виде суммы

q = qQJtq°-

(1.207)

Источник q2 соответствует нулевым начальным возмущениям энтропии, a q<j — нулевым начальным возмущениям завихренности. Обратим внимание на одно формальное отличие источников q& и Источник q<z вообще не зависит от возмущений энтропии, по­ скольку сг=0 , если сг<0) = 0 . Однако источник q зависит от возму­ щений завихренности даже при Q(0) = 0 , поскольку начальные воз­ мущения энтропии порождают возмущения завихренности, что сле­

дует из решения в разд. 1.5.

В силу линейности уравнения (1.63), а также граничного усло­ вия (1.65) справедлив принцип суперпозиции звуковых полей, по­ рождаемых источниками q& и qaj а стало быть, можно рассматри­ вать уравнение (1.63) с источниками q$> и qa отдельно друг от дру­ га. В связи с тем, что различие между обоими случаями с точки зрения общей структуры теории не носит принципиального харак­ тера, в дальнейшем рассмотрим только источник 9 2 , соответствую­ щий излучению звука возмущениями завихренности при их взаимо­ действии с неоднородным основным потоком.

С учетом изложенного источник q2 определяется выражением

(1.206),

где

 

 

 

 

Л = [Й'хК],

(1.208)

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

Q'

в соответствии

с формулами

(1.154),

(1.155) и (1.162) в отсутствие начальных возмущений энт­

ропии можно представить в виде

 

 

 

Q's=—

[2i0) +

BQ(n0)];

(1. 209)

 

(Г. 210)

а ; = Со(Ю Q(0)

( 1. 2 1 1 )

ео(к0)

v

причем начальные возмущения Q(°) являются произвольными функ­

циями т, ф, г [см. (1.134)], удовлетворяющими условию

(1.161).

На основании приведенных в разд. 1.5 формул векторного анализа в естественных координатах источник q2 можно представить сле­ дующей формулой:

[dQ.' dQ' \ т/з ,

Обратим внимание на то обстоятельство, что источник q& вооб"- ще не зависит от возмущений £2/, что впрочем следует уже из фор­ мул (1.208) и (1.206). Это означает, что винтовые вихри, т. е. вих­ ри, у которых вектор завихренности £2' коллинеарен с вектором ско­ рости V основного потока, не дают никакого вклада в излучение звука и это тем более интересно, поскольку имеется тенденция к; образованию винтовых вихрей в неоднородном потоке. В соответст­ вии с формулой (1.209) возмущения £2/ отличны от нуля даже в> том случае, если начальные возмущения £2S<°) отсутствуют, посколь­ ку они порождаются в результате взаимодействия поперечной сос­ тавляющей £2Пс неоднородностями основного потока.

С учетом (1.209) — (1.2 1 2 ) источник q можно выразить непос­ редственно через начальные возмущения завихренности Я<°>. После­ довольно простых преобразований имеем

 

 

Qo(VO^

V2

 

 

 

 

 

. QO(KO)

с2

 

 

 

+

2

Ко

«Эф

Ко к

0)

(1.213)

 

т0

дг ’

 

где m — Q0(V)V\- т0 и Vo — параметры

набегающего однородного

потока; R — радиус кривизны линии тока; до'(У)

— означает про­

изводную плотности go основного потока

по скорости V, а осталь­

ные величины определяются согласно формулам

(1.26),

(1.157),

(1.145).

 

 

 

 

 

 

Таким образом, получено самое общее выражение для источни­

ка звука в

уравнении

Блохинцева—Хоу

(1.63) в

отсутствие на­

чальных возмущений энтропии.

Когда скорость V основного потока много меньше скорости зву­

ка (V 7c«l), первым

членом в выражении (1.213) можно прене­

бречь по сравнению

с остальными членами, поскольку QO( V)

~QOO= const. Источник в этом случае можно записать в виде

аа<0)

QoVo

Ко

< 0)

(1.214)

дг

дг

 

В однородном потоке (V=V0) источник q2 приобретает вид

Л<2(°)

<зд(°)

(1.215)

<72 = 60^ 0 —^-----— = K0 rotQ'

 

дф

dz

-

Согласно результатам, приведенным в разд. 1.4, этот источник не приводит к излучению звука.

Наиболее сильное отличие источников уя и qQ наблюдается в области существенной неоднордности основного потока и оно сос тавляет ту эффективную часть источника q$f которая связана с из­ лучением звука. Отсюда следует, что источник звука эффективно локализован в области существенной неоднородности основного потока, т. е. вблизи обтекаемого тела на расстояниях порядка его характерного размера. В этой области радиус кривизны линии то­ ка достигает своих минимальных значений, а интегральный пара­ метр 0 достигает максимума.

Согласно физической постановке задачи начальные возмущения задаются выше по потоку от обтекаемого тела в области однород­ ности основного течения, где возмущения завихренности удовлет­

воряют более простому, чем

(1.60), уравнению

 

^

+ (1/oV)«o=0 .

(1.216)

dt

 

 

Начальные возмущения завихренности можно представить в ви­ де интеграла Фурье по пространственно-временным гармоникам вида ехр (—кot+ikx), т. е. в виде разложения по плоским вол­ нам. В связи с этим рассмотрим эволюцию возмущения заданного выше по потоку в виде отдельной гармоники и исследуем вклад, который вносит возмущение в источник звука q2 .

Пусть в набегающем на обтекаемое тело однородном потоке на­ чальные возмущения завихренности заданы в виде плоской гармо­ нической волны

tii= Q ,ne - l"t+lkx

(1.217)

где k x = k lx+k2y + k3z.

на бесконеч­

Будем считать, что вектор скорости У0 однородного

ности потока совпадает с осью х. Условие соленоидальности возму­

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

d iv Q o = 0 требует, чтобы

постоянный вектор

Я*п был поперечным,

т. е.

 

 

 

 

Qink = 0 .

(1.218)

Подстановка выражения (1.217)

в уравнение

(1.216) дает

закон

дисперсии малых возмущений завихренности

 

 

 

^ = V 0k{.

(1.219)

Рассмотрим, как изменяются

начальные

возмущения

вида

(1.217) в области неоднородности основного потока. В рамках об­ щего решения, представленного формулами (1.209) — (1.211), не­ обходимо выбрать произвольную векторную функцию Я<°) так, что­ бы это решение в области однородности выше по потоку переходи­

ло в решение вйДа (1.217). Нетрудно показать, что в этом случае­ произвольная векторная функция й (0) от переменных т, ф, 2 долж­ на быть выбранЯ следующим образом:

 

k y 0x-\-^ / n o - f k3z,

(1 . 2 2 0 )'

где

(1. 2 2 1 )

a есть вектор» составляющие которого постоянны относительно локального репера естественных координат es, en, ez:

+

(1. 2 2 2 )

и совпадают с составляющими вектора Й<я в формуле (1.217), при­ чем удовлетворяется условие поперечности (1.218). Заметим, что вектор й гп поворачивается вместе с локальным базисом так, что его составляющие остаются постоянными. Выбор векторной функ­

ции й<°> вида (1.2 2 0 ) согласуется

с условием

соленоидальности.

(1.161). Из формул (1.209) — (1.211), (1.220),

и (1.218) следует

й 'у ф =

0,

(1.223)

что есть более общее условие поперечности вихревых волн в неод­

нородном потоке. Можно непосредственно видеть, что

в области

однородности потока выше обтекаемого тела фаза волны

в ( 1.2 2 0 )

совпадает с фазой в формуле (1.217). Таким образом,

решение

(1.209) — (1.211), (1.220), (1.221) соответствует начальным воз­ мущениям вида (1.217).

Рассмотрим теперь, какой источник звука qs будет соответство­ вать вихревым возмущениям, заданным в однородном потоке в виде

плоских гармонических волн. Для этого необходимо

подставить

формулы (1.2 2 0 ) и (1.2 2 1 ) в общее выражение для

источника

(1.213) через начальные возмущения завихренности. В результате будем иметь следующее выражение для источника звука qs:

 

' 6р (Р ) V

V2

 

. Qo (^о)

С2

- м

mWo2 {k2- b k x)Q(г0)

 

то

 

( 0)

(1.224)

lk3v (jQ'n

где Q(°) определяется формулой (1.2 2 0 ).

Приведенные выражения для источника q$> завершают полное определение всех величин, входящих в неоднородное уравнение Блохинцева—Хоу в случае отсутствия начачьных возмущений энт­ ропии. Эти выражения могут быть далее конкретизованы при рас­ смотрении различного рода частьых случаев течений.

В общем случае решение уравнения Блохинцева—Хоу постро­ ить нельзя. В конце разд. 1.4 было уже указано на эти затрудне­ ния, причина которых заключается в более сложной структуре вол­ нового оператора Блохинцева по сравнению с классическим опе­ ратором Даламбера. Мощным средством решения классического волнового уравнения с источником в правой части является исполь­ зование функции Грина, которую удается построить в ряде практи­

чески важных случаев. Построение, функции Грина для уравнения Блохинцева—Хоу без каких-либо дополнительных предположений не представляется возможным, и только использование всякого ро­ да упрощающих обстоятельств может привести к успеху. Так, на­ пример, при использовании теоремы взаимности Хоу [74] построил низкочастотную асимптотику функции Грина для источника, нахо­ дящегося вблизи обтекаемого тела, когда точка наблюдения уда­ лена. В высокочастотном пределе использование метода геометри­ ческой акустики [5] и его различных модификаций также позволяет получить решение уравнения Блохинцева—Хоу. Таким образом, в отношении получения решений этого уравнения и извлечения фи­ зических следствий предстоит еще довольно значительная работа.