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

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

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

Г л а в а 2

ШУМ СВОБОДНОГО ТУРБУЛЕНТНОГО ПОТОКА

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

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

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

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

тельно к турбулентной струе, а также приведены результаты экс­ периментальных исследований акустического поля. Авторы отчет­ ливо представляют необходимость дальнейших теоретических ис­ следований процесса образования шума в струйных течениях, и осо­ бенно, исследований формирования пространственного распреде­ ления шума. Действительно, постановка проблемы шума свободно­ го турбулентного потока в виде акустической аналогии Лайтхилла является точной и включает все эффекты взаимодействия звука с потоком. В принципе, неоднородное волновое уравнение (1.194) должно полностью описывать звуковое поле турбулентного потока,, в том числе и турбулентной струи. Трудность заключается в томг каким образом выделить необходимую информацию из этого уравнения. В то же время, постановка проблемы шума турбулент­ ного потока может быть осуществлена с помощью уравнения, в ко­ тором более отчетливо разделены процессы излучения и распрост­ ранения звука. Так, развитие общей теории с учетом влияния ха­ рактеристик среднего потока, как отмечалось в гл. 1, сначала было* осуществлено Филлипсом, а затем в более полном виде Лилли, ко­ торые выявили влияние градиентов плотности, температуры и сред­ ней скорости на излучение и распространение звука в турбулент­ ном потоке. В этой постановке проблемы аэродинамического шума считается, что источники шума не непосредственно взаимодейству­ ют с окружающей средой, а подвергаются воздействию окружаю­ щего потока. Однако полученное неоднородное конвективное волно­ вое уравнение не имеет пока общего решения, поскольку содержит комбинации сложных членов и включает трудности,' свойственные сжимаемым, вихревым, нелинейным потокам, преодоление которых идет в основном путем введения упрощенных моделей источников, шума и существенной идеализации структуры турбулентного по­ тока.

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

2.1.1. Неоднородное волновое уравнение. Акустическая аналогия Лайтхилла

Рассмотрим основные положения теории шума свободного тур­ булентного потока [80, 81, 83], в общем виде впервые сформулиро­ ванной Лайтхиллом, и те упрощающие предположения, которые были использованы при анализе основного уравнения с целью полу­ чения решения применительно к турбулентной струе.

Уравнение неразрывности в неподвижной системе координат при отсутствии источников массы имеет следующий вид:

dQ . й(ои|) = 0

(2 . 1)

dt т dxi

где Q — плотность; t — время; щ — скорость течения жидкости в направлении х{.

Уравнение сохранения импульса' (количества движения) при отсутствии внешних сил представляется

 

 

 

dJj T +

£ i ( е ^ + р .7)=°-

м

где Pij=pbij-\-,t\

( — 0/у -|—

 

 

— тензор напряжений от сил

давления р и вязкости; г] — коэффициент сдвиговой вязкости;

 

да-

да/

да-

8 ,. =

(1

При

i =

j

Кронекера.

6 , ; = — --1------;

дХ)

{п

при

. ,

.— символ

"

dxj

dxt

"

(0

1 ф J

v

Дифференцируя уравнение неразрывности (2.1) по времени, а уравнение количества движения (2 .2 ) по пространственной коор­ динате Xi и вычитая один результат из другого, имеем

32Q д2 дР dxidxj

(Qu,Uj + P ,j).

с2

Вычитая

затем из обеих

частей

этого

уравнения величину

 

hr

получаем неоднородное волновое уравнение

 

0 dxidxj

 

 

 

 

 

 

 

 

 

 

d*Q_

2 - g g - _

d2TU

 

(2.3)

 

 

 

dt2

0

fix’]

dxidxj

 

 

 

 

где

со— скорость распространения звука в невозмущенной

среде,

тензор Tij^QUiUj-^riPij —рЪц)-{-(р — CoQ)bij

представляет

собой

разность напряжений в потоке и напряжений в покоящейся среде. Видно, что левая часть уравнения (2.3) описывает распростра­ нение звука в покоящейся среде, а правая часть характеризует ис­ точники звука, образующиеся вследствие действия поля напряже­

ний в жидкости.

Основная идея Лайтхилла базируется на формальном преобра­ зовании уравнений (2.1) и (2.2) в уравнение (2.3) и заключается в том, что точные уравнения движения жидкости могут быть записа-

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

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

Решение уравнения (2.3) рассматривается при условии прене­ брежения в Tij вторым и третьим членами по сравнению с пер­ вым членом. С физической точки зрения это означает пренебреже­ ние эффектом поглощения звука из-за вязкости и теплопррводности. Кроме того, плотность среды в потоке считается постоянной и равной плотности окружающей среды Q ~ QO. При таких упроще­ ниях, допустимых только в случае дозвукового изотермического по­ тока, может быть получено решение уравнения *(2.3) для плотнос­ ти или давления вне поля турбулентного потока в предположении, что тензор Tij^QoUiUj известен.

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

Введенная модель потока, как показывают практические при­ ложения, применима в ограниченном диапазоне углов наблюдений и частот, и не может быть использована при оценке высокочастот­ ного шума, распространяющегося в области вблизи оси струи. По мере увеличения скорости потока рефракция звуковых волн прояв­ ляется все в большей степени, и аналогия Лайтхнлла становится слишком грубой. При скоростях, соответствующих М>1, анализ неоднородного волнового уравнения еще более усложняется, по­ скольку уже нельзя использовать допущение о пренебрежении ча-

стью слагаемых в Тц. Кроме того, становится неоправданным предположение о постоянстве средней плотности среды в струе.

Для неизотермического потока, даже если принять, что величи­ на Рц может быть приближенно заменена /?6Zj, анализ Тц также не может быть проведен на основе введенных приближений, по­ скольку помимо учета изменения средней плотности среды при рас­ смотрении процессов излучения и распространения звука необхо­ димо учитывать изменение скорости звука в зоне смешения струи.' Так, даже в предположении, что давление и плотность являются изэнтропийными функциями, при использовании соотношения dp=^ = C2G?Q, где с — местная скорость звука, для исключения члена Pijc2q6ij необходимо в различных частях струи использовать раз­ личную величину с. Поэтому результаты оценки шума неизотерми­ ческих струй, а также струй с переменной плотностью, основанные на теории Лайтхилла, не всегда подтверждаются эксперименталь­ ными исследованиями.

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

2.1.2. Неоднородное конвективное волновое уравнение

Первая существенная попытка сформулировать основы теории пэродинамического шума с учетом взаимодействия движущегося потока и звука, генерируемого этим потоком, была предпринята Филлипсом в 1960 г. [93]. Идея его работы состоит в том, чтобы, исходя из уравнений движения и второго закона термодинамики, представить основное уравнение шума турбулентного потока в ви­ де неоднородного конвективного волнового уравнения, в котором правая часть, характеризующая источники шума, в отличие от представления Лайтхилла не содержит плотность среды в явном

виде.

Уравнение неразрывности (2.1) в неподвижной системе коорди­ нат имеет вид

J_

I dul - n

(2.4)

Q D t ' dx,

где D __ а ,

а

Dt

dt + U i dXi

С помощью выражения для энтропии единицы массы идеаль­ ного газа

S .= c v In-4 -

(2.5)

Q7

 

уравнение неразрывности может быть записано несколько иначе. Действительно, представим уравнение (2.4) следующим образом

D In Q | dti[ __Q

Dt дх[

Уравнение (2.5) может быть переписано в виде

1пд=— In р -----—S,

УсР

где y= c p/cv — отношение удельных теплоемкостей

единицы мас­

сы газа соответственно при постоянном

давлении

и постоянном

объеме. Из последнего уравнения имеем

 

 

PlnQ =

\

Р р ____1

DS

 

Dt

ур

Dt

ср

Dt

 

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

dui __

1

Dp ,

1

DS

 

dxi

уP Dt

Cp

Dt

 

Если использовать соотношение — — = — , где a = \n (p/p0);

р Dt Dt

Po — давление в окружающей среде, то уравнение (2.6) записы­ вается

ди^

1

Do j

1_

DS

дх i

у

Dt

ср

Dt

Уравнение количества движения (2.2) в неподвижной системе координат также имеет вид

Put

___ 1_

дРи

Dt

Q

дх i

Взяв от (2.8) производную по пространственной координате и используя соотношение

 

 

J _ P _ = _D _ д _ , д^ ' _ д _

(2 9)

 

 

dxi Dt

Dt дхi

dxi

dxj

 

можно получить

 

 

 

 

D

даi

dui

 

 

 

 

Dt

dxi

dxj dxi

d x i)~ d x i 1

Q

d x jl \ j

3

( 2. 10)

После применения к уравнению (2.7) оператора D/Dt и подста­ новки результата в (2.10) с учетом соотношения

1

др

с2

да

Q

d x t

у

дХ[

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

(2. 11)

Левая часть этого уравнения описывает распространение звука и соответствует волновому уравнению движущейся среды относи­ тельно переменной \n(plp0). Эффекты изменения местной скоро­ сти звука и взаимодействия звука с потоком, связанные с процес­ сом распространения звука, присутствуют в левой части уравне­ ния, что является естественным в противоположность уравнению Лайтхилла, где эти эффекты включены в правую часть.

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

( 2. 12)

D t 2 д л :Д d x i ) d x j d x i

Для адиабатного процесса

р/о* = Ро/Оо

уравнение

(2.12)

при-

нимает следующий также часто используемый вид:

 

 

D2 \nQ

д /

2 д In Q\ _

duj

d“j

,2

13)

D&

dxt \

дхi )

dxj

dxi

 

 

Если мгновенную скорость щ разложить на среднюю U{ и пульсационную и/ составляющие Ui = Ui + Ui\ то с учетом только чле­ нов, содержащих Пульсации скорости и представляющих нестацио­ нарные взаимодействия, правая часть уравнения (2.13) записыва­ ется

^ _ | _ 2

dli dCjrj

dxj dxL

дх i dxj

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

взаимодействие типа сдвиг — турбулентность. Физическая сущ­ ность членов правой части более подробно будет рассмотрена да­ лее.

В случае исследования распространения звука в потоке с уче­ том только распределений средних скоростей и температуры и пре­ небрежения рассеянием и конвекцией звука пульсационными сос­ тавляющими потока левая часть уравнения (2.13) существенно упрощается. Если также принять модель двухмерного турбулентно­ го потока, характеристики которого изменяются 'только в направ­ лении х2, нормальном направлению истечения хи U\ = U\(x2), U2 = = С/з= 0, с= с(х2), то уравнение (2.13) представляется

Пр" In Q

_д_

с 2 (х 2) д In Q

да]

du^ _|_2—

(2. 14)

* ^

Dt2

дх,-

dxt

dxj

dxi

dxi

dx2

где

Dt dt dxi

Для изотермического дозвукового потока уравнение (2.13) мо­ жет быть сведено к уравнению Лайтхилла. Так, члены левой час­ ти уравнения (2.13) записываются

In Q

Dt2

g i i l i + 2 и,

dnaQ

4- Uiuj

d2lnQ +

^

d№

 

 

dxidt

 

dxjdxj

Dt dxi

_d_ I

2 d In Q\

_ c2

& Q

(dQ_ \ 2

 

d x i \

d x i

j

Q

dx]

Q2 [ d x j ’

 

Du/

d In о

 

/o o \

где член —------ -с помощью уравнения количества движения (2.о)

Dt

dxt

&2

Id/Qо \2

представляется в виде

в2

[d x j

 

 

При преобразовании правой части уравнения (2.13) используем соотношение

dui du.j__ d2(UjUj) _

_d_ /ди Л _

du± dJJj

dxj dxL

dxidxj

1dxj \dxi)

dxLdxj

которое с помощью уравнения неразрывности (2.4) записывается

dui

ди] =

&(UjUj) _

/ £ Q \2

_d_

D Q\

dxj

dxL

dxtdxj

Q2 \ D t '

1dxt \ Q

Dt j

Тогда в результате использования полученных представлений членов левой и правой частей уравнения, (2.13) для случая М<^1 получается основное уравнение теории Лайтхилла.

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

кового потока постановку задачи в форме акустической аналогии Лайтхилла.

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

В соответствии с предположением Лилли [79] только член, со­ держащий .пульсационную скорость в квадрате и представляющий взаимодействие типа турбулентность — турбулентность, может рас­ сматриваться в качестве члена, характеризующего источники шума турбулентного потока. Второй же член, линейный относительно пульсационной скорости и характеризующий взаимодействие типа турбулентность — сдвиг, может быть представлен в виде функции переменной а и включен в число линейных членов левой части урав­ нения как член «сдвиговой рефракции».

Так, если мгновенную скорость, давление и температуру пред­ ставить в виде суммы средней и пульсационной составляющих, то

уравнение (2.12) может быть записано

 

 

Р2д'

- 2 y p

f ^ = L ^

0.

(2. 15)

DP

дх

[ oxj

 

 

где DjDt=d/dt + Uid/dXi-, с — средняя скорость звука; а' — флукту­ ирующая часть переменной а; при малой величине пульсационного давления р' справедливо а= In(1 -\-р'/р0) ~ р'/ро=°'; функция Lo(x, t) представляет собой совокупность членов, содержащих про­ изведения пульсационных составляющих. Отметим, что поскольку местная скорость звука пропорциональна корню квадратному из местной температуры с2~Т, то величина с2 характеризует среднюю температуру, а с'2 — пульсационную составляющую температуры.

Уравнение (2.15) является неопределенным, поскольку содержа­ щиеся в левой части величины а' и и' взаимосвязаны. Эта неопре­ деленность устраняется путем применения оператора D/Dt к обеим

частям конвективного волнового уравнения.

 

части

уравнения

В

результате дифференцирования

правой

(2.12)

использования

соотношения

(2.9)

и уравнения

количества

 

 

 

 

Da!

 

с2 да

 

 

движения, записанного в виде —----------- -— , получаем

 

 

 

 

D t

 

у

dxi

 

 

 

D_

/ daj_

 

D_

ди±

,

duj

D_

< 4 / _

 

 

D t

\ d x j

dxi /

dxi D t

dxj

 

dxj

D t

dx;

 

__

 

(&_

 

 

 

Г— (—

- \ 4 - ^

 

d x i l d x j

\ у d x J ' d x j

dxk\

d x j \ d x t \

Y

dxj) '

dxt dxk J

Тогда после переноса в левую часть члена правой части, содер­ жащего переменную а, уравнение (2.12) представляется в следую­

щем виде:

 

 

А (с2

+

— • (2-16)

D fi

D t

d x j

dxt d x j \ d x j

dxj dxk dx/

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

Dfi

± ( - 2 ^1)_1_2с2 ^ L - ^ =

2 Y ( —

дх^ дх^ дх{ I

Dt дх Д дх i)

дх2 дх±дх2

\dxj дх^дх-,

,(2. 17>

Правая часть уравнения (2.17) по меньшей мере квадратична пульсационной скорости и уже не является ее линейной функцией.

Если оператор D\Dt применить к уравнению (2.15), то в резуль­ тате имеем

D3*'

 

D

±

f a

дх-,)

i 2 ^ i

Dfi

Dt

дх,

\

dxj

dUj_

dU b

d u j

2

д а )

ъ _

/ д и л

dxj

дх,-

dxk

 

dxi

Dt

(2. 18)

 

\ d x j)

где функция Li(x, t)

содержит произведения и члены более высо­

кого порядка пульсационных величин и отражает взаимодействие пульсаций скорости, давления и температуры. Уравнение (2.18): представляет собой наиболее общую форму неоднородного кон­ вективного волнового уравнения и впервые было получено Лил­ ли [79].

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

о_ * (“>})

дх2 dxidxk

 

Dt

dxidxj

 

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

с принятой терминологией член, содержащий

только пульсационные скорости и отражающий

взаимодействие

турбулентность — турбулентность,

представляет

составляющую

«собственного» шума, а член, включающий взаимодействие турбу­ лентность — градиент средней скорости, представляет составляю­ щую «сдвигового» шума.

Отметим, что в более общем случае члены правых частей урав­ нений (2.16) и (2.18), содержащие пульсационную скорость в пер­