Аэродинамические источники шума
..pdfВеличину т в дальнейшем будем называть фазой. Последнее выра жение можно получить непосредственно из уравнения (1.62), за писанного в естественных координатах
|
^ |
+ V { S ’ ф )1 7 = а |
|
|
(1Л 35) |
Первому |
интегралу r\ = const для уравнения |
(1.131) |
соответствует |
||
ф = const |
в естественных |
координатах для |
уравнения |
(1.135). Об |
|
щее решение уравнения |
(1.135), т. е. решение, зависящее от произ |
||||
вольной функции, в соответствии с (1.132) |
и |
(1.134), |
можно пред |
||
ставить в виде |
|
|
|
|
|
|
|
О= ^ а (0)(г, ф, г), |
|
|
(1. 136) |
где произвольная функция а(0) может быть определена, если возму щение энтропии задано как функция времени и координат на неко торой поверхности, пересекающей линии тока. Пусть это будет по верхность $= 5о(ф), тогда а(0)(т, ф, г) должна быть заданной функ цией.
Решение (1.136) имеет простой физический смысл: энтропийные возмущения переносятся вдоль линий тока с местной скоростью •основного течения без изменения величины. Другими словами, энтропия отдельных частиц жидкости в основном течении остается постоянной. В неоднородном течении скорость вдоль линии тока изменяется и, следовательно, изменяется расстояние между от дельными частицами жидкости, что приводит к изменению во вре мени пространственных производных энтропии в системе коорди нат, движущейся с местной скоростью отдельной частицы.
Первыми интегралами уравнений характеристик для системы (1.127) — (1.129) являются также T = const, *il) = const и z = const. В характеристических переменных система уравнений для возму щений завихренности (1.127) —(1.129) сводится к обыкновенным дифференциальным уравнениям, из которых возмущения определя ются интегрированием вдоль линии тока. Чтобы их проинтегриро вать, необходимо сначала вычислить производные о, входящие в правые части уравнений (1.127) —(1.129). Учитывая (1.136), а так
же (1.134) для производной |
ds |
|
будем иметь |
|
|
|
|||
|
|
|
|
|
|
|
|
||
да |
да^ |
дтз |
|
dq(Q) |
1 |
(1. |
137) |
||
ds = с р |
дх |
ds ■= — Сг дх |
V(s, ф) |
||||||
|
|
|
|||||||
Несколько сложнее вычисляется производная до/дп: |
|
|
|
||||||
да |
^ Iда^ |
дф |
, да^ |
дх \ |
(1 |
. 138) |
|||
дп |
}\ |
дф |
дп |
дх IhT) |
|||||
|
|
|
Для вычисления дх/дп вместо (1.134) удобнее использовать выра жение (1.133), помня, что нижний предел интегрирования в обоих
интегралах, вообще говоря, имеет зависимость $о(Ф) и Ы л)- В°с‘ Пользуемся также формулами (1.94) и (1.96)
дтз |
1 |
|
г |
|
? |
(s'’0 ) d % ' = — -х |
У (£о* ° ) |
dF0 j- |
||
дп |
|
|
|
|
dii |
|||||
лт)(£. л) |
<*] |
J |
^ 2(6'. |
^(S o. ч) ^ 2(£о. т1> |
||||||
|
1 |
|
to СП |
|
|
|
,KW " |
|
|
|
|
|
|
е |
|
|
|
! _ |
« ^ - 5 |
> Л '. |
|
|
+ 2 ------ —I |
Лт) |
1 |
л) |
(1.139) |
|||||
|
—г |
<?1) |
|
|||||||
|
М ^ ) |
|
J |
|
|
V |
2 (e', Т)) |
|
|
|
|
|
£о (Ч) |
|
|
|
|
|
|
|
Перейдем теперь обратно к переменной s, используя при этом со отношение (1.113)
-2т (s, ф) ^ |
ds'db/ds' |
(1. 140) |
т ($', ф)К(5 '„ ф> |
$0 (Ф)
где m{s, <1»)=е0(«. ф)^(5, Ф).
С учетом (1.140) и (1.110) уравнение (1.138) можно переписать в; виде
да |
до(°> |
m(s, |
ф) - 1 |
dsp |
|
дп |
дф |
|
|||
|
|
|
|||
+ |
2 m(s, ф) J |
ds'db/ds' |
(1. 141) |
||
m (s', ф Щ * ' , ф ) |
|||||
|
|
|
|
•So (Ф)
Поскольку далее предстоит интегрировать вдоль линии тока, выраженин для производных с удобно записать, выделив явную зави симость от переменной 5 . Учитывая (1.136), (1.137) и (1.141), за пишем
|
да __ / да \ Vo |
(1. 142) |
|
|
~ д 7 ~ [ ~ д 7 ) о ~ ’ |
||
да |
|
|
|
|
(1 |
. 143) |
|
|
|
дп
(1. 144)
-дг£ - = (\- 7дг- )/ 0• где индексом нуль отмечены значения соответствующих величин в
точке s= s 0 (i|)), а также введен |
интегральный параметр |
0 , равный |
||||
6(5, |
Ф) = |
2m 0V 0 |
'i |
ds'db/ds' |
(1. 145) |
|
т (s', <P)V(s', ф ) |
||||||
|
|
So (Ф) |
|
|||
|
|
|
|
|||
На характеристике, т. е. на линии тока a^const, 2 =const при |
||||||
x=const, уравнения |
для |
возмущений завихренности |
(1.127) — |
(1.129) сводятся к обыкновенным дифференциальным уравнениям вдоль линии тока.
Уравнение (1.128) на линии тока с учетом (1.144) записывается в виде
(1. 146)
Уравнение (.1.127), учитывая (1.145) и (1.146), можно представить следующим образом:
|
|
д |
|
|
1 |
|
|
|
|
|
|
ds |
|
Q0V |
m0V0 |
|
|
|
, Ь147) |
|
|
|
н |
|
|
|
|
||
Аналогично для уравнения |
(1.129) с учетом |
(1.142), (1.143) и |
|||||||
(1.145) |
имеем |
|
|
|
|
|
|
|
|
д |
Г |
V2 |
|
|
|
|
|
|
|
ds |
[ Qo |
2т0ср |
|
|
|
|
( L m > |
||
где |
|
^ (s , |
|
Ф)= |
j [т (s', |
<]>) — т0] |
1 |
(1.149) |
|
|
|
Qo(s', 4) |
|||||||
|
|
|
|
|
So (Ф) |
|
|
ds |
|
|
|
|
|
|
|
|
|
|
|
Из уравнений |
|
(1.146) —(1.148) |
имеем недостающие первые инте |
||||||
гралы системы |
уравнений |
(1.127) —(1.129) |
|
||||||
|
|
|
|
|
|
|
|
|
(1.150) |
(1.151)
я*
6 о
const.
(1.152)
Используя эти выражения, запишем общее решение всей системы уравнений. Если на поверхности з = з 0 ('Ф), пересекающей линии то ка, заданы функции
2 , = 2 <0>(/, ф, 2 );
2„=2<0)(*, ф, z); |
(1.153) |
Qz = Q i 0)(t, ф, 2),
тогда общее решение для возмущений завихренности с учетом (1.150) — (1.152) может быть записано в виде
(1. 154)
(1.155)
| 1 ^ |
da№ |
1 |
dv^ |
dsp j? (1/)1 |
(1. |
156) |
||
Vo |
dr d i / / |
V0 |
d v |
dty |
J ’ |
|||
|
|
где функции Q«» и tf<°> берутся в момент времени т, определяемый уравнением (1.134).
Функции Q(°) и а<°), вообще говоря, не могут быть заданы пол ностью произвольно, поскольку решение должно удовлетворять ус ловию соленоидальности (1.130). Подставляя решение (1.154) — (1.156) в уравнение (1.130), получим
< Н у й = - —— L
т о К дх
<?2 i0) |
d s 0 |
дх е о ( ^ о ) dty |
|
Q o (K 0) |
Vi—V\ |
Vp |
|
|
Ко |
•Я7 + |
|
|
VRp a ? 4 |
V |
L |
|
|||
|
|
|||||
mVо |
dQ<0) |
|
v20 |
|
|
|
m^V |
d z |
V R 0 d z |
|
|
||
-m C |
d V |
I |
dSp |
<?2a(°> Q |
(1.157) |
|
J Qo (K ) |
J d<\> |
dzdx |
|
|
Сформулируем теперь достаточные условия, которые надо нало жить на функции £2<°) и о(0), чтобы уравнениен (1.157) удовлетворя лось во всем потоке. На поверхности 5=5о(ф) это уравнение мож но переписать в виде
1 «<°> _ |
о£» |
< |
0) |
т0 |
е0 (^о) |
dsp |
М г0( ) |
Ко |
(jg<°> |
|
dt |
Rp |
|||||||||
|
дф |
|
|
|
d z |
R Q |
d z |
(1. 158)
С учетом уравнения (1.158), в котором необходимо заменить t на т, уравнение (1.157) можно представить следующим образом:.
d i v a = ( i - i . ) X |
! L ( ^ . |
* 1W |
dsp |
Ко |
do<°> |
|
||
V |
mpj |
V \ Rp |
dx |
-9o(Vo) di/ |
Rp |
dz |
|
|
|
Qo (Ко) |
v i —vl |
v |
d V ) |
dsp |
Э2о<0) . |
|
|
|
■m \ |
0. |
(1.159) |
|||||
|
|
|
Qo(V') |
dty |
dxdz |
|||
|
|
|
V0 |
|
|
|
|
При отличных от нуля Q л0) о<0) и их производных это уравнение бу дет удовлетворяться во всем потоке, если
Ro |
0 . |
(1. 160) |
Эти условия заведомо удовлетворяются в удаленной от обтека емого тела области выше по потоку, где течение однородно. Таким образом, если начальные возмущения Я<°>, которые задаются в об ласти однородности потока на удаленной плоскости, перпендику лярной к линиям тока, удовлетворяют условию
1 |
dQ<0) |
dQ{n0) |
dQ<0) |
(1. 161) |
- ------оН----------------------— |
||||
V0 |
dt |
дф |
dz |
|
то поле возмущений завихренности будет соленоидальным во всем? потоке. Начальные возмущения энтропии на удаленной плоскости могут быть заданы произвольно.
С учетом только что сказанного последние члены в уравнении (1.156) должны, очевидно, быть опущены. После этого оно приоб ретет вид
_ QoflO |
) ( 0) |
V2 |
да"» |
~Qo{V) |
V 2 -V \ |
|
2Vl |
(1. 162) |
|||
Qo (Vo) |
|
dv |
|
2 |
Кроме того, параметры, отмеченные индексом нуль, исключая р0* приобретают смысл параметров в набегающем потоке.
Так же как и возмущения энтропии, малые возмущения завих ренности Q переносятся вдоль линии тока с местной скоростью ос новного потока, однако в отличие от возмущений энтропии их вели чина не остается постоянной для отдельных частиц жидкости ос новного потока. Изменение вдоль линии тока компонент завихрен ности определяется изменением параметров основного потока: ло кальным параметром — скоростью и интегральным параметром 0 . Величина завихренности £2П в соответствии с уравнением (1.155) определяется только локальным параметром V, поэтому завихрен ность жидкой частицы будет иметь равные значения в точках линии тока с равными скоростями в моменты времени, соответству ющие прохождению жидкой частицы этих точек. Этого нельзя ска зать о компонентах Qs и S2Z, поскольку они определяются интеграль ным параметром 0 . Завихренность в направлении линии тока Qs линейно зависит от завихренности Qn с коэффициентом пропорцио нальности, равным интегральному параметру 0 .
Важной особенностью решения (1.154) — (1.156) является так же то, что в неоднородном потоке возмущения завихренности зави сят от возмущений энтропии. Если даже начальные возмущения за вихренности отсутствуют, т. е. й<°)=0 , то это не означает, что за вихренность будет отсутствовать во всем потоке, напротив, она будет индуцироваться начальными возмущениями энтропии. Все эти особенности решения для возмущений завихренности не проти воречат известной теореме Кельвина о сохранении циркуляции по замкнутому жидкому контуру.
Обратим внимание еще на одно обстоятельство, связанное с ли нейностью системы уравнений (1.60), (1.61), (1.62) для возму щений завихренности и энтропии. Произвольное решение этой сис темы уравнений может быть представлено в виде линейной супер позиции двух решений, одно из которых определяется начальными возмущениями завихренности Q<°> при s ^ = 0 , а другое — возмуще ниями энтропии при й<0)^=0 .
1.6. ГЕНЕРАЦИЯ ЗВУКА ТУРБУЛЕНТНЫМ ПОТОКОМ ПРИ МАЛЫХ ЧИСЛАХ МАХА.
ТЕОРИЯ ЛАЙТХИЛЛА И ПАУЭЛЛА
Применим полученные в разд. 1.3 уравнения к случаю аэроди намической генерации звука турбулентным потоком при малых числах М2«с1. Поток предполагается неограниченным в том смыс ле, что в нем отсутствуют твердые или податливые границы, либо
эффекты, связанные с их присутствием, пренебрежимо малы. Такой поток называется свободным. Предположим также, что завихрен ность и неоднородности энтропии отличны от нуля в ограниченной области, характерный линейный размер которой определяет макро скопический масштаб турбулентности L. Примером такого потока, очень важным в практическом отношении, является струя воздуха, истекающая в спокойную атмосферу (затопленная струя) или в не однородный спутный поток. В этом случае характерным масшта бом L может служить диаметр сопла.
Турбулентный поток при малых числах М является акустически компактным источником, т. е. масштаб L много меньше характер ной длины волны акустического излучения. Действительно, из со ображений подобия следует, что число Струхаля, образованное по
частоте f, соответствующей |
основному энергонесущему масштабу |
|
L и характерной |
скорости |
потока У, Sh = fL/V является, вообще |
говоря, функцией числа Рейнольдса Re=VL/v и числа М =V/c |
||
|
Sh = /(R e,M ). |
|
В нашем случае, |
поскольку |
рассматривается турбулентный поток, |
числа Re весьма велики, а число М мало, так что можно допустить, что число S h = / ( o o , 0) является постоянной величиной порядка единицы для всех чисел Re, больших некоторого критического. В экспериментах со струями это число равно 0,3, причем под ха
рактерной скоростью потока понимается |
скорость на срезе сопла, |
|||
а частота соответствует максимуму спектра излучаемой |
звуковой |
|||
мощности. Таким образом, частота |
имеет порядок f~V/L, а соот- |
|||
ветствующая длина волны Х |
с |
L* |
много больше |
характер |
|
^ > 1 |
ного размера источника.
Поместим начало координат где-нибудь в области турбулентно го потока, в случае струи, скажем, — на срезе сопла. При выполне нии неравенств M2 <Cl, L<^X в области г<СА, поток можно рассмат ривать в приближении несжимаемой жидкости, где г — расстояние от начала координат. Волновое уравнение Блохинцева — Хоу (1.48) в соответствии со сделанными ранее предположениями мо жно значительно упростить, пренебрегая малыми членами.
Сделаем необходимые порядковые оценки для величин, входя щих в уравнение (1.48). Рассмотрим сначала область г< 0 , где за ведомо отличны от нуля завихренность и неоднородности энтропии. Для производных по координатам и времени в этой области имеем djaxi = 0 (L-1), d/dt = 0 (f)= 0 (VL~X). Следовательно, для матери альной производной имеем оценку D/Dt = 0 (V *L~X) , отсюда также получаем Dv/Dt=0(V2L~l) .
Если поле скорости v представить в виде суммы |
соленоидаль- |
ной и безвихревой составляющей |
|
<o=vs-\-vf1 |
(1. 163) |
гд е d iv t f * = 0 ; r o t f l ^ O , |
|
46
то в области |
безвихревой составляющей скорости можно пре |
||
небречь. Действительно, из ургЦшения непрерывности |
|
||
|
J - |
-2 £ \{ -d iv « '= 0 |
(1.164) |
|
QC2 |
Dt |
|
с учетом сделанных ранее оценок, а также принимая во внимание, что P ~ QV2, легко получить следующую оценку vr 0 (V М2). Таким образом, если пренебречь всеми членами, пропорциональными М2 в обеих частях уравнения (1.48), то на расстояниях г<^Х для эн тальпии торможения получим уравнение Пуассона
—v2/? = dlv{[fl X Vs] — TyS}, |
(1. 165) |
при этом также предполагается справедливым уравнение адиабатичности DSJDt = 0.
Рассмотрим теперь область r^>Ly в которой по предположению
отсутствуют завихренность и неоднородности энтропии. |
Следова |
||
тельно, правая часть уравнения |
(1.48) в этой области обращается |
||
в нуль. Для производных по координатам на расстояниях |
r » L |
||
справедливо d/dxi = О (Аг1), а |
для производной по времени оста |
||
ется справедливой прежняя оценка. При записи левой части |
(1.48) |
||
можно различать два случая: |
1) на расстояниях r^>L |
среда по |
коится, как это имеет место в случае затопленной струи; 2 ) имеет ся внешний, вообще говоря, неоднородный потенциальный поток (струя в спутном потоке). Вторым членом в левой части уравнения (1.48) в большинстве случаев можно пренебречь. В первом случае от него остается только часть, пропорциональная пульсационной скорости, которая относится к членам, описывающим нелинейное распространение звука. Для потоков с малыми числами М нели нейные эффекты, связанные с распространением звука, который порождается самим потоком, пренебрежимо малы. Во втором слу чае, когда имеется внешний поток, этот член имеет порядок
1 |
Dv |
п |
/ |
/Ц2М2 |
\ |
. 166) |
-------- v = |
o |
I |
АХ |
(1 |
||
с2 |
Dt |
|
) |
|
где Л — масштаб неоднородности внешнего потока, поле скорости в котором будем обозначать через Vo; т — параметр, равный отно шению скорости внешнего потока V0 к характерной скорости тур булентного потока V. В случае струи в спутном потоке этот пара метр называют параметром спутности. В обоих случаях, о которых
говорилось ранее, материальная производная |
имеет одинаковый |
порядок |
|
— — = 0(Х-!). |
(1.167) |
сDt
Однако в первом случае эта производная |
записывается в виде |
D/Dt = d/dt, а во втором случае следует |
учитывать скорость У0 |
внешнего потока D/Dt=d/dt-\-Vо/ ^ “ 7 • |
|
Таким образом, в области r 3 >L с учетом сделанных оценок эн тальпия торможения В удовлетворяет в первом случае уравнению Даламбера
(iS^vJ)B=0' <1л«»
где с0 — скорость звука в покоящейся среде, а во втором случае — конвективному волновому уравнению
< 1 Л 6 9 )
Легко видеть, что уравнение для величины В, справедливое равно мерно во всем пространстве, можно записать в первом случае в виде
^ - V 2) 5 = div{[Qx ^ ] - 7 ’vS}, |
(1.170) |
а во втором случае
+ |
B=d,v |
(1.171) |
В самом деле, в области r^>L правая часть обоих уравнений (1.170) и (1.171) обращается в нуль, т. е. справедливы уравнения (1.168) и (1.169), а в области гС ft, с уравнениями (1.170) и (1.171) можно повторить проведенные ранее рассуждения, с помощью ко торых можно убедиться в справедливости уравнения (1.165). Заме
тим также, что в промежуточной области L<r<^K=— , энтальпия
М
торможения удовлетворяет уравнению Лапласа
ДД= 0. |
(1. 172) |
Кстати, существование самой этой промежуточной области свя зано с малостью числа М. На расстояниях, больших по сравнению с L, но малых по сравнению с ft, (L<^r<^%) с учетом убывания В с расстоянием (точнее той ее части, которая зависит от времени), можно записать общий вид решения уравнения (1.172)
в=в0+ |
г |
OXi \ Г I OXi OXj \ г ] |
(1. 173) |
|
|
|
A ( a f f l ) + |
- £ - ( S i £ W ) + ..., |
|
где величины р, Diy |
Qijy ...от координат |
не зависят, В0— постоян |
ное значение энтальпии торможения вне турбулентного потока. Ограничимся далее в этом разделе рассмотрением случая, ког
да, в области r^$>L среду можно рассматривать как покоящуюся, которая, как показано ранее, описывается уравнением (1.170). Будем также для простоты предполагать, что неоднородностями энтропии можно пренебречь. В частности, это предположение явля ется вполне реальным в случае изотермической затопленной струи.
Иными словами, далее будем рассматривать излучение звука за вихренностью, которое описывается уравнением
( 4 г ~ - V2) fi= d iv [Q X v*] = q. |
(1. 174) |
Это уравнение может быть преобразовано к виду, который впервые был получен в работе Лайтхилла [80] на основе соображе ний, приведенных в разд. 1.1. Правую часть уравнения (1.174) с учетом формулы (1.141) для соленоидальной составляющей скоро сти можно переписать следующим образом:
<7 = d iv [Q х г>]= — ^ ^ |
----—v2^2> |
(1. 175) |
dxidxj |
2 |
|
при этом индекс 5, обозначающий соленоидальную составляющую скорости, далее будем только подразумевать. Учитывая выражение (1.175) для источника, представим его следующим образом:
(1.176)
dxidxj ’
где Pij = vivj —
Решение уравнения (1.174) для безграничного пространства, удов летворяющее условию излучения, может быть записано в виде
В (х, t ) - B 0 = |
B'(x, t) = - ± - \d * y |
' ’ D с |
|
|
4я |
J |
R |
_ - L f |
dyidyj |
R |
(1. 177) |
4 я J |
|
где B0— постоянное значение энтальпии торможения в отсутствие завихренности; 7? = |де—у |. Интегрирование по У (1.177) эффектив но сводится к области, в которой отлична от нуля завихренность, т. е. \м\ < L. Важное свойство интеграла (1.177) состоит в том, что производные по координатам можно вынести за знак интеграла
’ d‘yd' =
На расстояниях г=\х\, удовлетворяющих условию r^>L, выраже ние для R можно приближенно представить в виде
R = \ x - y \ ~ r — ^iji—\-... = r-\-0(L), |
(1.179) |
т. е. интеграл (1.178) можно переписать
•у ( Л « - г г ^ { т ! М * ' - £ Н (г>>ц- (,Л80>
Характерное время Т, за которое заметно изменяется величина Pij(y> *), равно T = 0(L/V), поэтому выражение (1.180) можно за писать также в виде
(г»д (,Л81>
пренебрегая различием в запаздывании во времени от различных точек области \у\< L, в которой отлична от нуля функция Рц. Оно не превышает величины порядка L/cQи поэтому при М<с1 мало па
сравнению с Т. При |
формулу (1.181) можно представить |
||||
в форме (1.173), поскольку |
при — <^ — = — = Т |
вообще мож- |
|||
но пренебречь запаздыванием |
с0 |
С0 |
V |
|
|
|
|
|
|
||
5'(д:^) = - 4 |
- ( |
— |
) |
(£«г«Х ), |
(1-182) |
OXiOXj |
\ г |
] |
|
|
|
где |
|
|
|
|
|
Qi}{t) = - ^ |
\ p u(y ’ t)d3y- |
(Ы 83) |
Из сравнения с формулой (1.173) можно заключить, что в случае генерации звука завихренностью функции ц(/) и Di(\t) тождествен но равны нулю. С помощью рассуждений, аналогичных сделанным, в разд. И и 73 из [21], можно показать, что функция ц(^) отличня от нуля, если в нестационарном завихренном потоке имеются по датливые тела, объем которых может изменяться, a Di(t) отлична от нуля, если в потоке имеются твердые тела, не изменяющие сво его объема. В свободном турбулентном потоке разложение (1.173)> называемое мультипольным разложением, начинается с Qzj.
На расстояниях |
пренебрегая членами порядка |
в формуле (1.181) |
г~х можно вынести за знак дифференцирования* |
а производные по Х{ вычислить согласно следующей формуле:
д _ |
\ |
дг |
д |
__ |
Xj |
J_ |
d |
dxi |
CQ |
dxt |
dt |
— |
r |
Co |
dt |
В результате будем иметь
B '( x ,t) = |
1 - |
* |
Q>A t |
- - ) ■ |
|
_ 1 |
|
Г |
dxidxj |
|
Co / |
X i X j |
<Э2 Г |
Ри (у,{ — r |
OC |
||
4лсд |
гЗ |
dfi j |
\ |
c0 |
|
(1. 184)
»Х). (1. 185)
Решение общего вида, соответствующее формуле (1.173) на рассто-
50