Аэродинамические источники шума
..pdfГ л а в а 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, где с — местная скорость звука, для исключения члена Pij—c2q6ij необходимо в различных частях струи использовать раз личную величину с. Поэтому результаты оценки шума неизотерми ческих струй, а также струй с переменной плотностью, основанные на теории Лайтхилла, не всегда подтверждаются эксперименталь ными исследованиями.
Наблюдающееся противоречие теоретических предсказаний и данных экспериментальных исследований шума струи можно объ яснить упрощением механизма образования шума. В реальных ус ловиях источники шума турбулентного потока не погружены в по коящуюся среду, а подвержены эффектам взаимодействия с дви жущейся средой так, что мощность их излучения определяется также условиями локального окружения. Поиску формулировки проблемы шума турбулентного потока, обеспечивающей четкое разделение процессов генерирования шума и его распространения в окружающей среде и пригодной для изотермических и неизотер мических потоков в широком диапазоне скоростей, было уделено большое внимание в течение прошедшего двадцатилетия.
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), содержащие пульсационную скорость в пер