Файл: Прикладная математика [сборник статей]..pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 30.10.2024

Просмотров: 57

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

- эч -

5 . Лебедев Н.Н. Специальные функции и их приложения. М., Физматгиз, 1963.

6 . Седов Л.И. Механика сплошной среды.- Г. I . И ., "Наука", 1970.

А.Я.ФЁДОРОВ, В.Я.ЛОТАРЕВ

РАССЕЯНИЕ И ПОГЛОЩЕНИЕ ЗВУКОВЫХ ВОЛН СФЕРОЙ В ВЯЗКОЙ СРЕДЕ

Когда звуковая волна встречает препятствие, некоторая ее часть отклоняется от своего первоначального направления.

Разность между действительной волной и волной, которая сущест­ вовала бы при отсутствии препятствия, называют рассеянной вол­

ной.

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

превращается в тепло.

 

 

 

1

 

 

Определим скорость, с которой часть энергии звуковой волны

поглощается вследствие вязкости, а часть

рассеивается

благода­

ря присутствию сферического препятствия.

 

 

 

 

Пусть сферические звуковые волны, исходящие от точечного

источника Е , распространяются в вязкой

среде и встречают

на

своем пути абсолютно твердую и неподвижную сферу радиусом

а ,

центр

которой находится

на расстоянии

6

от источниКа^вдка Е

При установившемся движении (с временным множителем

1)

зада

ча

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

вязкой среды приводится

к решению двух уравнений

Гельмгольца:

 

 

 

 

. дср4-Лг[ф=0 ;

д<Р+/^Ф=0,

 

 

(I)

удовлетворяющих граничным условиям на поверхности сферы:

^ = 0 ;

Г0= О ;

J V 0 ,

( 2)


где «

! - с ч

 

: V т

3 5 л

* =

“\ Г

 

 

 

 

 

К ,

r 8 , V v

 

- проекции скорости точки на оси сферической

 

Р , в , U)

 

системы

координат, связанной

со

сферой;

 

-

сферические координаты;

 

 

 

о

-

частота

колебаний.

 

 

 

Из работы [ij известно решение данной задачи:

 

 

 

Я=0

(k%r)?n (cosBj£>-6 ^ 4-

 

 

 

 

 

 

 

*ik i "Hi A n (In +-1)hn (ki&)hn(klr)Pn (co$0 )p'6ii'; (3 )

 

 

'

 

 

> r)

 

где

hft{z),jf](z)

- сферические

функции Ханкеля и Бесселя

 

 

J^rj (cos 9)

 

П -го порядка;

о

 

-

полиномы Лежандра;

 

 

Aft

и Bft -

постоянные

коэффициенты, определяемые

по

д

 

 

 

следующим формулам:_ .

 

_ - П (п+Р/п [к^а)Ьп {кгаЬк^<Х]п(кух) \kiahb(к^акЬ^а)]

Ал

kfik'n(k.ja) [к ^ К (кгakhn[кь^\-п (п+ 1)ЛЙ (к^)кп [кьа)

 

в

_ k ih n ik iS jjJ k ia l+ k iA n h n (k i& )h n (k ^ c t)

(4)

 

Л

к 2кп (кг6)[k0ia/}l,l(k1a)fhftikla}]

 

ГДе

Jn (z)=

 

M

Ш hn (z)-

 

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

Л {p?+Ps)(v? r ^ t )(/f,

где ? Г Р °7 Г

F On

•давление и радиальная скорость па-

э

= ,

дающей сферической в'олны; .

ps,=p Q^ £ £ ; y^r= -J& ’^

давление и радиальная скорость волн,

 

 

отраженных от сферического препят­

 

 

ствия.

 

Так как

механическая

энергия в замкнутом

пространстве

остается постоянной, то среднее значение этого

интеграла пред-


- 96 -

ставляет энергию, которая поглощается благодаря трению жид­ кости. Энергия, рассеянная препятствием, определяется ин­ тегралом

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

4 - И 1 Р , У£г * р ^ 9 е )с1$.

(5)

s

Интеграл (5) берется по поверхности сферы, радиус кото­

рой меньше б ,

но значительно больше величины

, ]*>'* -

толщина вязкого

пограничного слоя:

 

.

h j w *

(б)

В этом случае потенциалы скоростей падающей и отраженной

волн имеют вид

 

 

сp p = i k i ^ { l n + { ) h n i k i6 )jn (Ajе)Рп ( c o s 8 ) e~6 i l ;

 

ф < г=ik. £

[In+1)Д я {kjS)hn(к{г)Рп(co s 8 )e■6ii

(?)

 

A

n=0

 

 

При составлении суммы Ptffyg + PgVifr. необходимо учитывать только такие члены, которые содержат пространственные сфери­ ческие функции одинакового порядка. Для длинных волн

 

с высокой степенью точности можно ограничиться

двумя первыми

членами в разложении (?)..

 

Найдем значение интеграла (5) при монопольной составля­

ющей падающей

волны. Для этого в формулах (7) нужно поло­

 

жить п = 0 .

Тогйа для действительных частей потенциалов

скоро­

стей падающей

и отраженной волн будем иметь

 

=ki [/о(A^jsin 6t-ji0{ki6)ca$6t}j0 (k%r);

(8)

<ps =

k ^ 0j 0{ k s ) - 5 0Tia ( k ^ t i a e t -

[60j 0 { k ^ ^ n d k ^ ) J c o s e tj,

где

a 0=H0j 0{ki&)-Kon °

 

 

 

60 = H 0n 0{k^ + K o j'o ( k t 6 );

 

Aa=Hn+iKD=

j o ( k

j

O i .

 

 

K ik ^ a )


-97 -

п0{к%6) - сферическая функция Нейиана нулевого порядка. Значение интеграла (5) в этом случае

 

 

Ig =

 

 

^ С

 

 

 

.

(9)

Применяя асимптотические формулы для сферических функций

Бесселя

при к^С=и>«\

и учитывая

[3 ],

ЧТО

 

 

 

 

 

 

 

n„(zh

COS z

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

7^ = -|-fiJ>oC0-aA Jflr6 .

 

 

 

 

 

(Ю )

Выражая этот результат через поток энергии падающих сфе­

рических

волн

 

Ро к

 

,

получаем

величину

 

 

 

 

 

 

 

- f - 5 u ? 4 к , а ) \

 

 

'

^ и ) .

которая

имеет

место и при отсутствии вязкости

(если считать Aj

действительным

числом).

 

 

 

 

 

 

Найдем теперь

значение интеграла (5)

при дипольной состав­

ляющей падающей волны. Для этого в формулах (7) нужно при­

нять

Л

= I . Тогда

получим действительные

части потенциалов

скоростей падающей и отраженной волн:

 

 

 

 

(PprsS/frjCOS

 

 

 

 

 

 

 

J[6Jts i n(ktr)e f - +

 

,q>jt= 3 A t co s э | [a ,1/

i

(

 

(A

^

 

 

 

 

+аг1я 1( А 1г ) ] с о з б ^ ,

 

 

 

 

где

 

 

ai =HJl (ki 6)-Klni lkl6);

 

 

 

°

 

 

 

 

 

 

 

 

Ai=Nl +iK i.

 

 

 

 

 

(13)

 

 

 

 

 

 

 

 

 

Согласно

работе

[3]

 

 

 

 

 

 

Л ( /) =

3 m z ~zi

c —

;

 

h i W B jiW + ith b ) * ( - - j r

 

Ai jbai [kxa -4A-!1Jb+2A1jbtft]idTiA1tf+4A1JbflJeos£1a4

1 ‘\ К ^ р 1а г ^ 1 к ^ а г \-1{1\рга 1- 1 к \ р га ‘>- к \ а 1)


- 98

-

 

 

В этом случае значение, интеграла (5)

 

/ г=-6Яр0С(5'

1

 

(is)

1+ к\61.

Н%.

Рассиотрии значение этого выражения в зависимости.отJ5c. Если $ й велико, т . е . толщина вязкого пограничного слоя много меньше радиуса сферы, то влиянием вязкости можно пренебречь.

Действительно, в этом случае

1

Т ^ с В 'Ц к ^ а ) 1 1+ ~к\¥

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

к \ а Ч \ а г и к\6г

(17)

Этот результат в совокупности с выражением (II)

опреде­

ляет потерю энергии в падающих волнах за счет рассеивания на

сферическом препятствии.

 

 

 

^

 

 

 

Е|зли значение а. имеет порядок

]?>

или меньший,

чем {У

f то вследствие

малости к \&

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

4ра t+ра

 

 

 

и потеря энергии, выраженная

через

поток

энергии

в первоначаль

ных волнах, оказывается равной

 

 

 

 

 

 

ska

1+

1

1

1

Т\аг.

 

(18)

 

ра

 

|Ь сзг

 

'Н16Ч

 

 

 

В этом случае выражением (II)

в сравнении

с

полученным

выражением вполне можно пренебречь.

 

 

 

 

Таким образом, формула (18)

определяет ту

часть энергии

в падающих волнах, которая поглощается вследствие вязкости

среды.

Если в формулах/(1 7 ),

(18)

перейти к пределу при<5-—°°,

то полученные результаты

совпадут

с

результатами

работы [2 ],

где определялись потери энергии при дифракции плоских волн на сфере.

Л и т е р а- т у р а

I . Цой П.И. Дифракция сферических звуковых волн на сфере в вязкой среде. - "Механика жидкости и газа ", 1973? вып. 1„