Аэродинамические источники шума
..pdfРис. 2.14. Спектры продольной пульсационной скорости (ширина полосы фильт ра 1/10 октавы)
Перейдем к переменной г\ = (у—0,5D)/x, связанной с переменной Tji= (у—у1)/6 соотношением
rli= (Tl+ tg ai)/(tga i + tga2) или Т11= 3,62г) + 0,362.
После подстановки в формулу (2.76) получаем представление обоб щенной функции в следующем виде:
/(г 1)=1-1,35т1-81г12+320т13. |
(2.77) |
Результаты расчета распределения пульсаций скорости в на чальном участке струи по формуле (2.77) находятся в соответствии с данными экспериментальных исследований (см. рис. 2.11). Мак симальная величина пульсаций скорости, определяемая из условия Г)!= 1/3 или riSO, наблюдается вблизи линии, проходящей через кромку сопла параллельно оси y = 0,5D. Некоторое превышение экспериментальных данных над результатами расчета наблюдает ся лишь вблизи внутренней границы струи, где определенное вли яние на величину пульсаций скорости оказывает начальная турбу лентность. Следует отметить, что при большой начальной турбу лентности >3% происходит увеличение и максимальной величины пульсаций скорости в зоне смешения струи.
В переходном и основном участках в области между внешней границей струи и линией y=0,5D продольные пульсации скорости
также можно представить с'помощью функ ции f(r\):
У и ' ~ |
= |
\' |
( г 7 |
г ) |
/(Т 1 ) |
( 2 - 78): |
и с |
|
и ,С |
/ шах |
|
|
для переходного участка (хн
У 7 Я |
( V I |
l / f a . /(л ) |
(2.79). |
£7С |
\ t / c |
'm ax ' X |
|
Рис. 2.15. Распреде ление пульсационных скоростей в попереч ном сечении струи:
/—продольные пульсации скорости; 2—поперечные пульсации скорости
для основного участка (л;п^ * )- В области между осью струи и линией
y=0,5D с удалением от плоскости среза сопла профили пульсационных скоростей: выравниваются так, что практически в ос новном участке величины пульсаций скоро сти в этой области можно считать постоян ными и равными значению на оси.
Спектры пульсаций скорости представ* лены на рис. 2.14, г д е ] / ^ — среднеквадра
тичная пульсационная скорость в полосе частот шириной 1/10 ок тавы. По мере удаления точки измерения от среза сопла или от оси струи в спектре турбулентности возрастает доля низкочастот ных составляющих. Увеличение скорости истечения струи приво дит к возрастанию в спектре турбулентности высокочастотных со ставляющих.
Приведенные данные указывают на то, что характерная часто та пульсаций скорости определяется зависимостью f ~ U / L , где в качестве характерной скорости следует взять среднюю местную* скорость потока, а в качестве характерного пространственного раз мера L — некоторый линейный размер источника, величина кото рого постоянна в каждом поперечном сечении струи и растет при' увеличении расстояния от среза сопла; в качестве такого харак терного размера обычно используется масштаб турбулентности. Например, частота максимума в спектрах турбулентности, изме ренных в третьоктавных полосах частот на линии г/= 0,5Д при уве личении расстояния от плоскости среза сопла подчиняется следую щему соотношению /тах=1,35£/с/х [6]:
Отмеченные ранее закономерности справедливы для всех сос тавляющих пульсационной скорости. Различие заключается лишь- в абсолютной величине пульсаций скорости. Так, практически мож но считать, что величины поперечных v' и тангенциальных со' пуль саций скорости составляют приблизительно 0,7 величины продоль ных пульсаций скорости и\ которые являются максимальными-
(рис. 2.15), например, )/"v'*/Uc^ 0,7 j /"и'7 и с.
Последнее утверждение равносильно тому, что, как это следует из соотношения (2.74), при определенной величине градиента сред ней скорости, или в каждой точке потока, величина пути, смещения*.
например, поперечных пульсаций скорости |
составляет |
-—'0,7 вё- |
личины пути смешения продольных пульсаций скорости |
/„, т. е. |
|
t v— 0,7/м. |
|
|
2.3.2. Пространственно-временные характеристики потока
Пространственно-временная структура потока определяется на основе исследований корреляционных зависимостей между пуль сациями скорости в различных точках потока. Так, статистические свойства поля пульсационных скоростей описываются с помощью пространственно-временных корреляционных функций
x)= u'i (0, 0)и}(£, х).
Для удобства представления результатов измерений величины корреляции целесообразно нормализовать посредством среднеквад ратичных величин пульсаций скорости с тем, чтобы получаемый та ким образом коэффициент корреляции не превышал по абсолютной величине единицы,
п \ ( 0, 0)иу(е, |
т) |
£/;(£» х) = |
(2. 80) |
V и\ (°) V |
а) (5) |
С помощью коэффициента корреляции £(£), определенного при нулевом времени задержки т = 0, 1исследуются пространственные характеристики поля пульсационных скоростей. При таких измере ниях обычно один из датчиков остается неподвижным, а другой дат чик в случае определения продольного коэффициента корреляции перемещается в направлении, параллельном оси струи, в случае определения поперечного коэффициента корреляции датчик пере мещается нормально оси.
В качестве характерного результата экспериментальных иссле дований на рис. 2.16 приведены продольные коэффициенты корре ляции продольных kux и поперечных kvx пульсаций скорости. Из сравнения видно, что наибольшими являются коэффициенты корре ляции продольных пульсаций скорости. Максимальная величина корреляции между различными компонентами пульсаций скорос ти достигается при минимальном сближении точек измерения, а с увеличением расстояния между датчиками величина корреляции убывает. Причем, чем дальше от среза сопла расположен непод вижный датчик, тем менее интенсивно падает величина коэффи циента корреляции с увеличением расстояния между точками из мерения (рис. 2.17).
Пространственные характеристики поля пульсационных скорос тей определяются в виде интегральных масштабов турбулентности, характеризующих объемы турбулентных вихрей или коррелирован
ные объемы,
Ео
L = \ k { \ ) d ^
Рис. 2.16. Зависимость про дольных коэффициентов корре ляции от расстояния между точками измерения
Рис. 2.17. Зависимость про дольного коэффициента корре ляции от расстояния между точками измерения:
y=0,5D; |
l—xlD= 6; |
2—xJD= 3; |
3— |
x/D= 2; |
4'—x,fD = 1; |
5-x/D = 0,5; |
D= |
|
= 40 M M ; UC—100 м/с |
|
где £o— расстояние между точками измерения, при котором коэф фициент корреляции принимает первое нулевое значение.
На примере результатов исследований пространственных харак теристик продольных пульсаций скорости видно, что как продоль ные, так и поперечные масштабы турбулентности незначительноизменяются в поперечном сечении (рис. 2.18). С увеличением рас стояния от среза сопла х, как это следует, например, из рис. 2.17г пространственные масштабы турбулентности возрастают. Для м ас штабов продольной турбулентности эти зависимости приближенно* могут быть представлены в следующем виде [67, 77]\
LaX={ 0,12... |
0,14) JT, |
8 |
£„„=(0,03... |
0,036)*. |
' ] |
Отметим, что понятие поперечного масштаба турбулентности LUiy>.
определяющего эффективный поперечный размер перемешиваю щихся объемов турбулентного потока, близко по смыслу понятию* пути смешения /ц., определяющего среднее расстояние в попереч
ном потоку направлении, на котором пульсационная скорость со храняет свое значение. Из сравнения на примере продольной тур
булентности результатов измерений пути смешения |
1и= 0,027л; и |
поперечного масштаба турбулентности Luy= (0,03 |
0,036)x вид |
но, что практически наблюдается хорошее соответствие величин и Luy.
Из результатов измерений пространственных коэффициентов: корреляции, например, рис. 2.16, следует, что величина соотноше ния пространственных масштабов поперечной и продольной турбу
лентности составляет приблизительно LV/LM~ 0,6 |
0,65. Это соот |
ношение также соответствует приведенному ранее |
соотношению |
для путей смешения поперечных и продольных пульсаций скорости /г^0,7/и.
Рис. 2.13. Изменение интетральных масштабов продоль ной турбулентности в попереч ном сечении струи
Коэффициент корреляции k(x) при нулевом .пространственном разделении £= 0 или коэффициент автокорреляции определяет вре менные характеристики поля турбулентного потока. Характер изме нения временных, так же как и пространственных, коэффициентов •корреляции с изменением аргумента (времени т или пространст венного разделения g) одинаков; с увеличением времени задержки •коэффициент автокорреляции уменьшается. Так же как и в случае пространственных характеристик, величина временного коэффици ента корреляции поперечных пульсаций скорости меньше времен ного коэффициента корреляции продольных пульсаций скорости. Аналитическая зависимость временных коэффициентов корреляции от аэродинамических и геометрических параметров струи будет да на позже после приведения результатов исследования закономер ностей движения вихревой структуры.
В турбулентном потоке с течением времени происходят как рас пад турбулентных вихрей, так и их перемещение в пространстве. Поэтому наряду с информацией о структуре турбулентного потока
в неподвижной системе координат необходима также информация
отурбулентности в системе координат, движущейся «вместе с тур булентными вихрями».
Если зафиксировать расстояние между датчиками Ах, то тур булентный вихрь пройдет это расстояние с некоторой скоростью конвекции UK за определенный промежуток времени т* = Дх/£/ь. Момент прохождения вихря через выбранную точку потока опреде ляется моментом наблюдения в точке максимальной величины кор реляции при изменении времени задержки подачи сигналов с дат чиков.
Практически определение скорости конвекции вихрей осущест вляется на основе анализа результатов измерений пространствен но-временных коэффициентов корреляции пульсаций скорости {рис. 2.19). При каждом определенном расстоянии между датчика ми Ах зависимость коэффициента корреляции от времени задерж ки х имеет вид кривой с максимумом, величина которого умень шается с увеличением Ах вследствие разрушения турбулентных вихрей. Скорость конвекции определяется исходя из установленных
значений расстояния Ах и времени задержки ттах, соответствую щего максимальной величине коэффициента корреляции,
iJ к = Ах/Хтах.
Рис. 2.19. Зависимость коэффициента пространственно-временной корреляции пульсаций скорости от времени задержки (£>='60 мм; U c = &2 м/с; #/0,5Z) = 1)*
Приведенные результаты экспериментальных исследований так же можно представить в виде кривых постоянных значений коэф фициентов корреляции на диаграмме взаимной зависимости рас стояния между датчиками и времени задержки сигналов [67] (рис.. 2.20). Тогда скорость конвекции можно определить как тангенс уг ла наклона прямой, соответствующей наиболее медленному изме нению коэффициента корреляции и характеризующей движущуюся: «вместе с вихрями» систему координат.
Исходя из данных таких исследований установлено, что изме нение скорости конвекции в зоне смешения струи как для продоль ных, так и для поперечных пульсаций скорости, подобно измене нию средней скорости потока (рис. 2.21); по мере приближения к внешней границе величина £/к уменьшается. Причем в области су ществования максимальных пульсаций скорости, т. е. при г) = 0, ве личина скорости конвекции составляет UKc^0,6t/c и мало отлича ется от величины местной средней скорости. В общем виде рас-
Рис. 2.20. Кривые постоянных значе |
Рис. 2.21. |
Распределение |
скорости |
||
ний коэффициент^ |
корреляции |
пуль |
конвекции |
в поперечном |
сечении |
саций скорости |
(D = 25 мм; |
0 с ~ |
|
струи |
|
«150 м/с; */£>= 1,5; y = 0 ,5 D )
ттределение скорости конвек_ции в начальном участке струи записывается.
и к/ и с= ср(л);
% |
II |
О.;
------------1
У =0,5П
-ср(.п)= 0,6 — 2,5^1 — Ют)2 + 50г13. |
Г- • |
• - 0=25м м ; Uс=75м/с |
|||||
|
|
(2,82) |
|
||||
|
|
|
о - |
46 |
125 |
||
Анализ |
изменения корреля |
|
*-П = ¥ 6м м ; UC = 6 2 M/ C |
||||
ционных функций в движущей |
7 200 ¥00 |
800 |
1600 |
3200 ft Гц |
|||
ся «вместе с вихрями» системе |
Рис. 2.22. Зависимость скорости конвек |
||||||
:координат |
показал, |
что |
ана |
||||
литическое |
представление |
вре |
ции |
от частотн |
|
||
менного масштаба |
турбулент |
приближенно имеет вид |
|
||||
ности по данным одних а в т о р о в |
|
||||||
|
|
|
t = Lxl0Л67£/К, |
|
|
(2.83) |
;а согласно сведениям ряда других авторов p=3jc/4t/c.
Оба приведенных соотношения означают, что характерная час тота пульсаций скорости /П=1/'Р обратно пропорциональна рас стоянию от среза сопла и прямо пропорциональна скорости потока. Такая же закономерность отмечается при рассмотрении спектров •турбулентности.
Представленные результаты исследований характеризуют уро вень современного понимания мелкомасштабной турбулентности е зоне смешения струн. Однако они дают все-таки несколько упро щенное представление о структуре турбулентного потока, посколь ку турбулентное движение не сможет считаться полностью подоб ным для различных частотных компонент.
Так, в результате узкополосного частотного анализа простран ственно-временных корреляционных функций установлено [15, 67], что высокочастотные составляющие пульсаций скорости конвектируются с большей скоростью по сравнению с низкочастотными со ставляющими (рис. 2.22). Кроме того, с увеличением частоты на блюдается уменьшение временных и пространственных масштабов турбулентности. Последнее обстоятельство указывает на тот факт, •что разрушение высокочастотных компонент происходит быстрее, чем низкочастотных компонент.
Изложенный экспериментальный материал дает возможность
.принципиально представить общий характер изменения микрост руктуры турбулентного потока в зоне смешения струи. На оснозе •принятого подхода к процессу излучения звука, как результату взаимодействия нестационарных турбулентных вихрей, и использо вания приведенных экспериментальных данных о турбулентности возможно провести оценку распределения источников шума в зоне смешения, спектральных характеристик шума и суммарного шума струи.
2.3.3.Крупномасштабная квазиупорядоченная структура
Взаключение необходимо отметить, что наряду с рассмотрен ной мелкомасштабной турбулентностью в зоне смешения струи су ществует крупномасштабная квазиупорядоченная структура, кото рая характеризуется наличием когерентных связей низкочастотных компонент в достаточно удаленных друг от друга точках потока. Крупномасштабная упорядоченная структура образуется в резуль тате смешения потоков со сдвигом средней скорости и проявляется
ввиде волн неустойчивости, постепенно развивающихся и скручи вающихся в вихри. Существование такой структуры доказано да же при низких уровнях начальной турбулентности в широком диа пазоне чисел Рейнольдса. Физически механизм развития крупно масштабных вихрей в сдвиговом слое смешения можно представить, следующим образом [88].
Вблизи среза сопла возмущения .потока первоначально вызыва ют возбуждения сдвигового слоя, из которых развиваются волнынеустойчивости (рис. 2.23). По мере перемещения вдоль сдвигового слоя волны неустойчивости усиливаются и скручиваются в вихри,, которые, захватывая окружающий поток, образуют крупномасш табное перемешивание. Такой естественный процесс в целом мож но рассматривать как нарастание крупномасштабной волны неус тойчивости из мелкомасштабного возмущения в тонком начальном: сдвиговом слое. Взаимное внедрение потоков происходит с по мощью вихрей, увеличивающих общую поверхность взаимодейст вия и способствующих тем самым более эффективному развитию мелкомасштабной диффузии. При удалении от среза сопла крупно масштабные вихри, образовавшиеся вследствие возбуждения сдви гового слоя, постепенно разрушаются.
Попытки предсказать шум турбулентной струи исходя из рас смотрения крупномасштабной структуры, основывались как на ана лизе неустойчивости границы струи, так и на представлении струн в виде ряда взаимодействующих вихрей. В настоящее время наи более распространенное мнение заключается в следующем: в до звуковых струях такая структура в основном управляет процессом генерации шума турбулентными пульсациями скорости, а сама не является непосредственным источником интенсивного акустического излучения [88]. Однако в ряде работ, например [106], предполагает ся, что неустойчивая крупномасштабная структура может генери ровать определенную часть шума струи, в частности, в области низ ких частот и в направлениях, образующих небольшой острый угол с направлением истечения. Окончательная точка зрения на роль крупномасштабной структуры в общем шуме струи еще не выра ботана.
В этой связи следует выделить также работы по определению реакции сдвигового слоя и звукового поля струи на искусственное акустическое возбуждение, осуществляемое часто на практике по средством подачи звукового сигнала из камеры избыточного дав ления перед соплом. Действительно, поскольку процесс смешения
Границы зоны смешен
Сиппо \ струи
I \-Cnapu8nHuir
Вихрей LОВоиэоСание Вихря
'-Вовлечение окружающего Воздуха.
-Осцилляция сдВигабого слоя
Рис. |
2.23. Схема |
разви |
Рис. 2.24. Изменение |
широкополосного |
||
тия |
сдвигового |
слоя |
шума |
струи при акустическом |
возбуж |
|
|
|
|
|
дении: |
|
|
|
|
|
______ |
невозбужденная |
с т р у я ;----------струя |
|
|
|
|
при |
низкочастотном воздействии |
(Shw0,3); |
|
|
|
|
------------- струя при высокочастотном |
воздей |
||
|
|
|
|
ствии (Sh«3) |
|
потоко-в со сдвигом можно представить как постепенное развитие крупномасштабной волны неустойчивости, то спусковым механиз мом такого развития может являться звуковая волна возбуждения. В ряде случаев при экспериментальных исследованиях в результа те низкочастотного воздействия наблюдалось значительное увели чение шума струи (рис. 2.24).
Вследствие того, что существование акустического излучения от волны неустойчивости не доказано, то ее можно рассматривать как усилитель акустической мощности струи при определенной час тоте возбуждения, а струю в делом как усилитель внутреннего акустического возбуждения. Причем ответная реакция волны не устойчивости при акустическом возбуждении небольшой интенсив ности является линейной и описывается с помощью линейной тео рии неустойчивости сдвигового слоя. При более интенсивном воз буждении реакция волны неустойчивости становится нелинейной, поскольку отбираемая от среднего потока энергия частично пре образуется в энергию турбулентного движения. Описание нелиней ного поведения сдвигового слоя заключается в следующем.
По мере возрастания уровня возбуждения увеличивается амп литуда движения сдвигового слоя вследствие увеличения интенсив ности Полны неустойчивости. Это вызывает возрастание величины локального мгновенного сдвига и, следовательно, увеличение ин тенсивности мелкомасштабной турбулентности. При этом происхо дит передача энергии от волны неустойчивости к широкополосной
турбулентности. |
Кроме того, |
при |
возбуждении |
возрас |
||
тает |
скорость |
расширения сдвигового |
слоя |
и, |
следова |
|
тельно, |
эффективная площадь |
поперечного |
сечения зоны |
смещения, что приводит к увеличению общей энергии турбулентно сти. J3 результате действия этих механизмов энергия волны неус тойчивости образующаяся благодаря энергии среднего потока, передается энергии турбулентного движения, определяющей звуко
вое поле струи. Следовательно, факт увеличения шума струи следу ет рассматривать как результат увеличения интенсивности излуче ния звука мелкомасштабной турбулентностью, а не как следствие появления нового источника шума.
Таким образом, воздействуя на волну неустойчивости с целью ее усиления, можно добиться более интенсивного развития погра ничного слоя и получить значительное увеличение широкополосно го шума струи. В то же время, воздействуя на волну неустойчиво сти с целью разрушения с момента ее образования, т. е. вблизи среза сопла, можно добиться существенного снижения шума струи [10].
Так, при -высокочастотном акустическом облучении создаются условия, препятствующие образованию крупномасштабной струк туры, и наблюдается менее интенсивное развитие пограничного слоя, выражающееся в уменьшении его ширины, снижении интен сивности турбулентных пульсаций скорости и менее резком паде нии осевой скорости при удалении от плоскости среза сопла. Эф фект уменьшения интенсивности смешения истекающей струи с окружающей средой при высокочастотном облучении наблюдается совместно со снижением широкополосного шума струи (рис. 2.24)
Различным характером возбуждения сдвигового слоя можно от части объяснить появление «избыточного» шума струи и расхож дение результатов измерений шума струи на различных установ ках. Причем возбуждение волны неустойчивости может происхо дить как вследствие акустических, так и турбулентных и темпера турных пульсаций. Практически внутренние возмущения различ ного рода присутствуют в любых установках, только при различ ных обстоятельствах превалируют те или иные возмущения.
Полная совокупность источников возбуждения содержится в турбореактивных двигателях. Так, камера сгорания и турбина создают температурные, акустические и турбулентные пульсации. Кроме того, во время течения по каналу двигателя, при отрыве по тока, при обтекании различных стоек и центральных тел образу ются турбулентные и акустические пульсации. Каждый из этих ис точников внутренних возмущений может привести к пульсациям тяги и массового расхода на срезе сопла, вызвать нестабильность сдвигового слоя и изменить звуковое поле струи. Это означает, что дополнительный или «избыточный» шум двигателя по сравнению с шумом «чистой» выхлопной струи обусловливается не только внут ренним шумом двигателя, но также и увеличением интенсивности излучения самой струи вследствие эффекта возбуждения. Отсюда, кстати, следует возможный способ снижения шума выхлопа турбо реактивных двигателей путем уменьшения интенсивности внутрен них возмущений [13].
Таким образом, неустойчивая крупномасштабная структура определяет процесс смешения струи с окружающей средой и обра зование мелкомасштабной турбулентности, являющейся генерато ром широкополосного шума. С помощью крупномасштабной струк туры как бы раскрывается внутренний механизм развития турбу