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

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

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

При рассмотрении членов правой части этого уравнения, т. е. при оценке источников шума, введем дальнейшие упрощения отно­ сительно структуры потока в зоне смешения турбулентной струн применительно к области наиболее интенсивного излучения. Так, известно [65], что поток в этой области близок к параллельному со значительным градиентом средней скорости. Кроме того, как будет показано далее, расчет интенсивности шума струи сводится к сум­ мированию интенсивностей излучения' объемов, пульсации скорости в которых коррелированы, т. е. турбулентных вихрей. В пределах каждого турбулентного вихря можно принять, что градиент сред­ ней скорости постоянен, а изменение скорости вдоль оси струи пре­ небрежимо мало. По мере удаления от центральной области зоны смешения такая идеализация потока становится все менее реаль­ ной, но это не имеет особого значения при оценке суммарного шу- -ма струи. Поэтому для поля средних скоростей можно принять сле­ дующее упрощающее приближение:

U i = (J (х2)\ и 2= и 3= 0,

(2.32)

где направление Х\ совпадает с направлением истечения струи. Далее введем предположение о несжимаемости среды при оцен­

ке интенсивности излучения звука. Будем считать, что

плотность

среды в турбулентной струе постоянна Q = Q C и

мало отличается от

плотности окружающей среды Q0.

 

(2.31), описывающие ис­

Тогда члены правбй части уравнения

точники шума турбулентного потока, представляются

 

<*211;

:6c^ L ;

^ - = 2 Qc^

dUj

(2.33)

dxi dxj

dUi

dxidxj

dxi

dx i

dx j

 

где производная

вследствие

(2.32) равна

dUi

 

 

dxj

 

 

 

dx2

 

dniij

dxidxj содержит вторую пространственную про-

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

Второй член — L включает произведение градиента средней

ОХI

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

dPi

логии член — — выражает «сдвиговый» шум пульсации скорости

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

Действительно, поскольку принято,

что

dUi

dUi

ди\J

dxj

6x2

ТО дх i

представляет производную поперечной

пульсационнои скорости в

дч2

продольном направлении -----.Рассмотрим перемещение элемента dxi

жидкости вниз по потоку из точки А в точку Б под действием сред­ ней скорости Uj и поперечной пульсации скорости и2/ (рис. 2.1). Ввиду того, что существует значительный градиент средней скоро­ сти, величина скорости U\ в точке Б меньше величины скорости в- точке А. Тогда вследствие замедления, а затем ускорения, возни­ кает переменная сила, пропорциональная градиенту средней скоро­ сти и поперечной пульсации скорости.

Очевидно, что излучение «собственного» и «сдвигового» шум'а должно быть наиболее интенсивно вблизи центральной области зоны смешения, где пульсационные скорости и градиент средней скорости принимают наибольшие значения.

Вследствие линейности неоднородного волнового уравнения (2.31) его решение, определяющее «сдвиговый» шум объема тур­ булентного потока Q, имеет вид

 

 

р)

(у.

I* — yl )

 

Q - Q о =

 

«о

d V(y),

(2.34)

 

 

 

 

4лс% дхi

 

 

\Х — У\

 

 

а решение, определяющее «собственный» шум,

 

 

1

д2

(

' \Х ~ У '

 

 

г

 

 

----------

 

 

 

 

 

 

 

со

dV (у),

(2. 35)

6 б0 4 я cl

dxi dxj

)

 

-3-1

 

 

 

где у — координата элемента dV

жидкости; х — координата точ­

ки наблюдения в звуковом поле;

[х~У\ — расстояние от

эле­

мента объема жидкости до точки наблюдения, величины Р; и Пг^ определяются в момент времени t—г/с0, т. е. в тот момент, когда

звуковая волна должна излучиться, чтобы достичь точки наблюде­ ния в заданный момент времени (рис. 2.2).

Выполним дифференцирование правых частей выражений (2.34) и (2.35). Так, имеем

r = \ x - y \ = = V { x i - y if\

 

 

 

 

 

 

 

 

d r

_

 

x i ■ У1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дхi

 

Г

J

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X i У1 .

 

 

 

 

 

 

 

 

 

 

d ( t —

r/c„)

 

 

 

 

 

 

 

 

 

 

 

 

 

d x i

 

 

 

С0Г

 

 

 

 

 

 

 

д

ГPj(y<t —Г/Со)-!

 

1

Г дР,

d(t

г / CQ)

г

р

 

 

 

дх -, [

 

 

г

 

J

 

Г2

1

dt

 

dxi

 

 

 

дх,

\

 

 

д

/ рз

 

 

 

 

 

дР)

,

Р1

\

1

 

 

 

 

 

 

т (

1

 

(2. 36)

 

 

dxi

\

 

г

)

 

 

Г2

(

с0

dt

 

г

)

 

 

 

 

 

 

 

 

 

Аналогично получаем, используя

(2.36),

 

 

 

 

 

 

 

<Э2

г П

, ; ( у ,

t

Г/CQ)1

 

д

 

|-

Xj

уj

I |

дпи

 

П<7 \п

 

д2

М

 

 

 

 

г

д

Г

*}— У) /

1

дПц

 

д х [ д х j [

 

 

 

 

d x t L

г2

 

\ с 0

dt

 

г ) ]

 

 

■У]

 

1 d 2 n U

d ( t

г / CQ) |

1

( д П и

d ( t — r / c Q)

Г —

 

 

г»

 

 

с0

dfi

 

 

d x L

г2

(:\ dt

 

 

dxi

 

 

 

 

 

 

 

 

 

 

 

 

П

 

дг

 

 

 

1

дПц

П/;

\

д (г2)

 

 

 

 

 

 

и

дх[ ) _

г2 — ( —

 

 

 

 

 

 

 

 

 

£2

/

П/у \

 

 

 

Ч со

 

 

 

 

 

 

 

 

 

( x i — yi)(xj — yj)

(

1

d2Xlij

 

 

3

dll/y

 

ЗП/у

. (2.37)

dxidxj

\

г )

 

 

 

гЗ

 

 

 

 

 

 

c0r

ctf

 

Г2

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

Для дальнего звукового поля решение уравнения (2.31), опреде­ ляющее «сдвиговый» шум, представляется

\ - Q o = — — \

( x i У д

д Р j ( У » t — r / c о) dV (у),

(2. 38)

4 п с 6 J

г2

dt

 

О2

 

 

 

я решение, определяющее «собственный» шум,

_J__ Г(x i— yi)(xj — yj)

д * П и ( у, t — г / с 0)

dV (у). (2.39)

Q-Qo

 

dt2

4ш?о \

г3

 

 

 

Проведенное дифференцирование правых частей выражений '(2.34) и (2.35) применимо, вообще говоря, только к турбулентным лотокам, движущимся с малыми дозвуковыми скоростями. В этом

случае источники шума можно считать практически покоящимися и поэтому подынтегральные выражения можно рассматривать в не­ подвижной системе координат. -

2.2.2. Учет конвекции источников шума

Поскольку в реальных турбулентных потоках наблюдается кон­ векция источников шума, то при вычислении подынтегральных вы­ ражений, содержащихся в (2.34) и (2.35), необходимо перейти от неподвижной системы координат к системе координат, связанной с движущимися источниками [80, 105]. Так, далее будет показано, что в турбулентной струе перемещение источников шума происхо­ дит со скоростями, распределение которых подобно распределению* средних скоростей.

Введем систему координат т|, перемещающуюся со скоростью

конвекции источников Uu. Учтем, что за время х = \ х —у\/с0, в те­

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

от источника,,

находящегося в точке у, до точки наблюдения х,

источник переме­

щается на расстояние xUu = b\K\x—y\> где MK= t/K/c0. Следователь­ но, неподвижная и движущаяся системы координат связаны соот­ ношением

Т1= у +

Мк|лг — у\.

 

 

(.2*40)

Тогда с учетом того, что

 

 

 

 

 

dV {i\) = dV (у) [l -

]

,

 

выражения (2.34) и (2.35), определяющие

«сдвиговый» и «собст­

венный» шумы, представляются

 

 

 

 

 

д

С

t'jW’t —r/Co)

 

(2.41)

б - 6 о =

J

Г — Мк ( X I (/;)

 

4лс%

11

 

2

П,7(у, t

 

 

 

1

f

г/ср)

 

(2.42)

е— е0= + cl дхiдх]

О

r - t \ * ( X i - y i )

 

 

 

Проведем дифференцирование

(2.41) и

(2.42)

по

аналогии с

тем, как провели дифференцирование выражений (2.34) и (2.35). Для этого необходимо определить производную dxi в системе координат т|

dr

_ x i У1 Л __ дуЛ

dxi

г \

dx-J

Из соотношения

(2.40)

при постоянной величине г\ находим

 

Х 1'

У1f 1

ду/ \

р.

ду1’

Мк (Xj у{)

дх^

к

г

\

dxi)

dxi

Мк (Х[ у{) г

Подставляя найденное - ^ L получаем dxi

dr _

Xj — у i

dxi

r — MK(xi i/,)

В результате дифференцирования подынтегральных выраже­ ний, содержащихся в (2.41) и (2.42), получаются выражения, по­ добные (2.36) и (2.37) соответственно. Для области дальнего зву­ кового поля решения неоднородного волнового уравнения, пред­ ставляющие «сдвиговый» и «собственный» шум, имеют соответст­ венно следующий вид:

1

Г

(*1 —yj)

dPjiyt, t —/-/с0)

dV(r[)\

(2.43)

4 я с30

M K ( X i — (/,)]2

дх

 

 

1

С

{ Xj — У i ) ( X j — Уj )

(?2 I I ,; (1), t — r/c0)

(2.44)

4ясд £

tr Мк(*,-—</<)Р

дт2

dV{x\).

 

 

Введем угол 0 между направлением конвекции источника шума со скоростью t/K=CoMK и направлением на точку наблюдения. По­ скольку индекс i в знаменателе выражения (2.44) в этом случае соответствует направлению движения источника, то можно запи­ сать

0 = arOCOS ( X i У г ) / Г .

Для дальнего звукового поля r= |JC—у | — |jf| =х, х,-—#,~х* и выражения (2.43) и (2.44) принимают следующий вид:

__1________ xt______г dPj(r!,t —r/c0) -dV{xi)- (2.45)

 

 

cl

 

r2 0 — М к cos 0)2

J

dv

 

 

U

 

 

Q

 

 

 

1

 

 

X j X j

f d2Jl ij(l ),

t — r/c 0)

 

 

4яс„

 

r3 (1 — MK cos 0)3 J

dV{x\). (2.46>

 

 

 

 

 

 

u

 

 

 

2

 

 

Тогда без учета членов взаимной корреляции получаем интен­

сивность / =

(QQ0)2 co/Qoдля «сдвигового» и «собственного» шумов-

соответственно

 

 

 

 

 

/ =

1

 

X i X l

Ц J

(T,, 0) ^ ( Л + l . T)rfK(4)j d V №

 

 

 

 

 

1бя2 Q0 CQ

x * ( l —

 

M K COS 0 )4 ^

 

 

/ =

I

XiXjXnXi

г

I (

д~ 11ij

(л ,0 )^ (л - Ы ,

,

* 6 ( 1 — Мк c o s 0)6

,

,

dx2

 

1 6 я 2 б 0 с “

J

J

’ dx * V

 

 

 

V L 00

 

 

X d V ( t|),

(2.47)

x ) d V ( x [ ) I x

" ]1

(2.48)

где черта сверху означает усреднение во времени; %— пространст­ венное разделение.

Рис. 2.3. Схема излучения зву­

Рис.

2.4. Схема

излуче­

ка конвектируемым протяжен­

ния

звука конвектируе­

ным источником

мым

точечным

источни­

ком

Интегрирование выражений (2.47) и (2.48) должно быть осуще­ ствлено, вообще говоря, по бесконечному пространству. Однако вследствие ограниченности области смешения турбулентного пото­ ка с окружающей средой при повторном интегрировании практи­ чески можно ограничиться эффективным объемом V этой области.

Из выражений (2.47) и (2.48) видно, что в направлении, об­ разующем острый угол с направлением движения потока, интен­ сивность звука больше, чем в направлении, образующем тупой угол с направлением движения потока. Поскольку при этом увели­ чение интенсивности звука больше, чем ослабление, то происходит увеличение суммарной акустической мощности, определяемой ин­ тегрированием по сферической поверхности фактора направленно­ сти (1—MHcos 0)"п, где /г= 4,6.

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

t'cQ= t c 0-\-t'N[Kc0cos б,

или t'/t = (\ — Мк cos б)-1.

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

Оценку

влияния

движения

Неподвижный

источника

на эффективность

 

излучения

можно

провести

 

также исходя из графического

 

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

(рис.

2.4).

 

Пусть точечный источник кон-

 

вектируется

со

 

скоростью

 

UK= c0Мк.

Для

того

чтобы

 

звук, излученный из точки у,

 

достиг точки х в тот

 

же

мо­

 

мент времени, что и излучен­

 

ный из точки у+Ау,

 

звук

из

 

точки у должен быть излучен

 

раньше. Разница

во

 

времени

 

излучения равна At=Ay cos0/co

Рис. 2.5. Влияние конвекции на эффек­

или

 

 

 

 

 

 

 

 

 

 

 

тивный объем источника и время излу-

ду с0

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

= и неподвижным вихрями удобно представить в координа­ тах t п у (рис. 2.5). За промежуток времени At вихрь перемеща­ ется на расстояние А// = с0МкД^ т. е. связь между координатами t и у определяется

dt =

1

ду

с0Мк

Координаты излучающего объема, из которых звук приходит и точку наблюдения одновременно, изображаются на схеме прямой Ы, тангенс угла наклона которой к оси определяется соотношением dt/dy = cos 0/со. Оценим изменение эффективного объема у1{ и вре­ мя излучения tu конвектируемого .вихря по сравнению со временем

излучения 10 и объемом уо неподвижного вихря.

иметь уку0 =

Используя приведенные соотношения, будем

^=tKCoMKi г/KCOS 0/соТогда

чисто геометрически

получаем, что

=

(1 — Мк cos 0)-1,

(2.49)

кУо

т.е. эффективный объем и время излучения вихря за счет конвек­ ции изменяется на множитель (1—MKcos0)_1.

Поскольку интервал времени в движущейся системе координат

изменяется пропорционально (1—M K CO S '0 ) - 1 , т о производные по времени в неподвижной и движущейся системах координат связа­

ны следующим образом: -^ - = (1 — MKcos 0)-1-^ -

dt

дх

Обратимся к выражениям (2.38), (2.39) и (2.45), (2.46), пред­ ставляющим решение волнового уравнения для неподвижных и конвектируемых источников шума соответственно. Видно, что изме­ нение плотности, определяющее «сдвиговый» шум или дипольные источники шума, пропорционально элементу объема турбулентного потока и первой производной по времени от Pj, а изменение плот­ ности, определяющее «собственный» шум или .квадрупольные ис­ точники шума,— элементу объема и второй производной по време­ ни от Пij. Элемент объема в движущейся системе' координат равен элементу объема в неподвижной системе, умноженному на (1—Мк cos О)-1.

Таким образом, эффект конвекции в случае дипольных источни­ ков шума, изменяя эффективный объем и интервал времени излу­ чения на множитель (1—MKcos О)-1, вызывает изменение плотно­ сти пропорционально (1—MKcos0)-2, а интенсивности излучения {1—MHcos0)-4. В случае квадрупольных источников шума эффект конвекции, изменяя эффективный объем и дважды интервал вре­ мени излучения на множитель (1—MHCOS B)- 1, приводит к измене­ нию плотности пропорционально (1—MKcosO)“3, а интенсивности излучения (1—MHcos0)-6.

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

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

Схематически влияние огрниченности области излучения на звуковое поле источников представим также в координатах t и у {рис. 2.6). Общее количество источников шума выражается дли­ ной отрезка ys. Количество источников, звуковые сигналы от кото­ рых приходят в точку наблюдения одновременно, представляется длиной отрезка у. Так же, как и ранее, из чисто геометрического рассмотрения видно, что отношение числа источников, дающих вклад в звуковое поле в определенный момент времени, к общему числу источников равно у1уъ = 1—MKcos 0.

Качественное объяснение изменения фактора направленности продемонстрируем на примере упрощенной модели турбулентной струи (рис. 2.7). Так, источники шума существуют только внутри некоторого объема и перемещаются относительно окружающей сре­ ды со скоростью UKy равной половине скорости распространения звука с0. Точка наблюдения в звуковом поле находится на оси

Рис.

2.6. Влияние ограничен-

Рис. 2.7. Влияние конвекции на величи-

ности

области излучения

на

ну эффективного объема турбулентного,

количество эффективных

ис-

потока

 

точников

 

 

струи. Положение источников, ‘ вуковые сигналы от которых дости­ гают точку наблюдения одновременно, изображаются линией N В момент t0 звук, излученный от среза сопла (сечение 1—/), достигает точку наблюдения в какой-либо определенный момент времени. Пусть за время t\10 фронт звуковой волны N перемес­ тился на расстояние, равное половине длины области излучения. Источники, которые в момент времени tQ>находились на срезе соп­ ла, за это время переместились вниз по потоку в сечение 2—2 на расстояние, равное половине расстояния от сечения 1— 1 до сече­ ния 33. Источники, находящиеся в момент t\ на линии N, в мо­ мент 4 располагались в сечении 22. Заштрихованная область по^ тока между сечениями 1— 1 и 22, т. е. половина объема, прой­ денного фронтом волны, определяет интенсивность бегущей звуко­ вой волны. Область потока между сечениями 2— 2 и 33 еще не

успела излучить звук, распространяющийся с фронтом N.

В момент t2 фронт звуковой волны достиг границы области из­ лучения турбулентного потока. Видно, что только половина всего объема области между сечениями 1— 1 и 33 внесла вклад в зву­ ковое поле. Другая половина объема (между сечениями 33 и 4—4) не успела излучить звук, приходящий в точку наблюдения вместе с фронтом N, так как источники шума уже вышли из обла­ сти излучения. Источники шума, пришедшие в рассматриваемую область излучения после момента t0, также не успевают излучить звук, который приходит в точку наблюдения в выбранный момент времени, т. е. вместе с фронтом звуковой волны N. Поэтому число источников шума, дающих вклад в звуковое поле, уменьшается на 0,5, что также следует из выражения (1—MKcos0) при Мк = 0,5, cos 0= 1.

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

муле для интенсивности шума дипольных источников будет иметь вид (/—MKcos0)-3 вместо (/—MKcos0)~4, а для интенсивности шума квадрупольных источников (1—MKcos0)-5 вместо (/—MKcos 0)~6. Еще раз подчеркнем, что происходит это не вслед­ ствие изменения эффективности излучения источников шума, а вследствие изменения числа источников, дающих вклад в звуковое воле.

Тогда выражения (2.47) и (2.48), определяющие интенсивность «сдвигового» и «собственного» шумов турбулентной струи соответ­ ственно, принимают вид

XiXt

 

 

L(,1+ 1’ x)dvw>1 ^ (4 );

\[ \

^ (ть 0 )^

1бл2р0Со х4(1 — MKCOS0)3 ^

| J

дх

 

 

 

 

(2. 50)

X i X j X k X t

 

дзпы

/ = •

(П.

x)d V (4) d V (IJ).

16л2д0СдЛ:б(1 — MKcos0)5 \ [

 

 

кL ° °

(2.51)

Следует отметить, что приведенные ранее рассуждения относи­ тельно изменения эффективности акустического излучения турбу­ лентного потока подразумевают идеальные нераспадающиеся во времени источники шума. При сверхзвуковых скоростях истечения в случае М„ cos 0 = 1 решения волнового уравнения (2.45) и (2.46) теряют смысл, так как формально интенсивность излучения в на­ правлении под углом 0 к направлению движения должна возрас­ тать до бесконечности. Это объясняется увеличением эффективного объема и времени излучения источников, поскольку, излучая звуко­ вую волну, источники продолжают перемещаться вместе с волной. Поэтому в какой-либо определенный момент времени в точку на­ блюдения должен приходить звук, излучаемый в течение неогра­ ниченного промежутка времени.

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

2.2.3. Сдвиговый и собственный шум струи

Обратимся к выражениям для интенсивности «сдвигового» и «собственного» шумов и покажем, что в случае стационарного слу­ чайного процесса, произведение первых производных по времени в (2.50) можно заменить второй производной по времени, а произ­ ведение вторых производных в (2.51) заменить четвертой произ­ водной.