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

Основы_акустики_Гринченко_Вовк

.pdf
Скачиваний:
105
Добавлен:
10.04.2023
Размер:
16.81 Mб
Скачать

2h(1)(ka ) = h(1)

(ka )h(1)

(ka ).

(8.81)

n

n 1

n +1

 

 

Тогда для мягкой сферы коэффициенты ряда (8.76) имеют такие при- ближенные значения:

 

 

 

 

 

C0

 

=

j0(ka)

 

1

 

 

 

 

 

 

 

 

 

 

 

 

h0(1)(ka)

1 ka

 

 

 

 

 

 

 

 

 

 

С

 

=

 

jn (ka)

 

(ka)2n +1

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

hn(1)(ka)

 

2n +1

 

 

 

 

 

= ka,

(8.82)

C0 , n > 0.

Отсюда, для мягкой сферы при (ka) << 1 в ряде (8.76) можно оставить только первый член; тогда

σ = 4πa2

(8.83)

s

 

и, как и в случае малого мягкого цилиндра, имеем рассеяние моно- польного типа. Но для сферы, как тела конечных размеров, величина σs не стремится к бесконечности, как это характерно для цилиндра бесконечной длины.

Рассмотрим теперь жесткую сферу при (ka) << 1. Соответствую- щие приближенные значения коэффициентов Cn имеют вид

C1 = j1(ka)

h1(1)(ka)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j1(ka)

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C0

 

 

 

=

 

 

j0(ka)

 

=

 

 

 

(ka) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

h(1)(ka)

h(1)(ka)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

j (ka)

j

2

(ka)

 

 

 

 

 

 

 

ka

(ka)2

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3ka

 

3 5

 

 

 

 

=

 

 

ka 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ka)

,

(8.84)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

h

(1)(ka)

h(1)(ka)

 

1

 

 

 

 

3

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ka

1

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

ka (ka)2

 

 

 

 

(ka)3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С

 

 

=

 

 

 

 

 

=

 

 

 

jn 1(ka) jn +1(ka)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jn (ka)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

h(1)(ka)

 

 

h(1)

(ka) h

 

 

(1)(ka)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

n 1

 

n +1

 

 

 

 

 

 

 

 

jn 1(ka)

 

 

 

 

 

 

(ka)2n +1

C

 

,C ,

n >1.

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h(1)

(ka)

 

 

 

 

(2n +1)

!(2n 1)!

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

n +1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Как и в случае малого жесткого цилиндра, в ряде (8.76) нужно удер- жать два первых члена:

σs

4πa2

 

+

3 (ka)6

=

7

πa2(ka)4.

(8.85)

(ka)2

1

4

 

9

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

521

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

Как видим, полное сечение рассеяния малой мягкой сферы в че- тыре раза больше площади сечения сферы, тогда как сечение рассея- ния малой жесткой сферы во много раз меньше ее геометрического сечения. Согласно (8.85) полное сечение рассеяния жесткой сферы пропорционально частоте в четвертой степени, что обусловливает усиление рассеяния в области высоких частот.

При рассеянии коротких волн (ka >> 1) как для мягкой, так и для жесткой сфер, полное сечение рассеяния

σ = 2πa2.

(8.86)

s

 

Расчеты величины σs /πa2 по точной формуле (8.77) приведены на рис. 8.7.

Отметим, что полное сечение рассеяния для мягкой сферы на низ- ких частотах в два раза больше, чем на высоких частотах (формулы (8.83) и (8.86)). Это можно объяснить тем, что на низких частотах па- дающая звуковая волна возбуждает колебания, которые равномерно распределены на поверхности сферы таким образом, что на преобра- зование мощности плоской волны в мощность рассеянной волны ис- пользуется площадь, равная площади поверхности сферы. На высо- ких частотах рассеянная мощность возникает вследствие неравно- мерного возбуждения сферической поверхности плоской волной. По- этому рассеянная волна разделяется на две отраженную и теневую волны; для каждой из них сечение рассеяния равняется поперечному сечению сферы.

8.8. Рассеяние плоской волны на звукопроницаемой сфере

Пусть в среду с параметрами ρ и с внесено сферическое препятствие радиусом а. Сфера наполнена средой с параметрами ρ1 и с1. Принимая во внимание опыт построения решения задачи рассея- ния плоской волны на идеальной сфере, можно решить задачу рас- сеяния плоской волны на таком сферическом препятствии.

Поле плоской волны определяется по формуле (8.65), а рассеянное поле можно записать в виде (8.71). Выделение множителя (2n + 1)in ничего не изменяет при решении задачи, ведь все будет откорректи- ровано при определении коэффициентов Cn, но при этом мы имеем возможность, использовать формулу (8.76) для сечения рассеяния. Итак, звуковое поле в среде, окружающей сферу имеет вид

 

 

 

n

 

 

(1)

 

 

 

p = p

 

+ p =

(2n +1)i

j

(kr ) C h

 

(8.87)

0

 

 

(kr ) P (θ).

 

s

n =0

n

n n

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

522

В отличие от идеальной сферы в данном случае звук проникает вглубь сферического препятствия. Поле внутри сферы запишем так:

p =

(2n +1)inD j

(k r )P (θ).

(8.88)

1

n =0

n n

1

n

 

 

 

 

 

 

Снова отметим, что выделение множителя (2n + 1)in осуществляется только ради унификации решения. Отметим, что k = ω/c, k1 = ω/c1. По- ле p1 строится только с учетом функции jn(k1r), сферические функции Неймана yn(k1r) во внимание не принимаются, поскольку при r 0 они стремятся к (вспомните, подобная ситуация сложилась при изу- чении колебаний круговой мембраны).

На поверхности сферы должны выполняться условия сопряжения

полей снаружи и внутри сферы:

 

 

 

 

 

 

 

 

 

 

p0 + ps = p1, r = a,

 

 

 

 

1

p

p

 

 

1

p

 

 

 

 

 

 

0 +

s

 

=

 

1

,

r = a.

(8.89)

 

 

 

 

iωρ

r

r

 

 

iωρ1 r

 

 

 

Подставив (8.87) и (8.88) в систему (8.89), получим следующие равен- ства:

j (ka) C h

(1)(ka) = D j (k a),

 

 

 

 

n = 0,1,2,...,

(8.90)

n

 

 

n

n

 

 

 

 

n n

1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

(1)

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ka) =

 

 

 

 

 

 

 

(k1a).

(8.91)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jn (ka) Cnhn

 

 

 

ρ1c1

Dn jn

 

ρc

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поделив (8.90) на (8.91), запишем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

(ka) C h(1)(ka)

 

 

 

 

ρ c

 

j

(k a)

 

 

 

 

 

 

 

n

 

 

n

 

n

 

 

 

=

 

1 1

 

n

1

 

 

.

 

 

 

 

 

 

 

 

 

 

(1)

(ka)

 

 

 

 

 

ρc

 

 

 

 

(k1a)

 

 

 

 

 

 

 

jn

(ka) Cnhn

 

 

 

 

 

 

 

 

jn

 

 

 

Отсюда искомые коэффициенты Cn

равны

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ c

 

j (k a) j

 

(ka)

 

 

 

 

 

 

 

 

 

 

 

 

11 1

 

n

 

1

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρc

j

(k a)

 

jn (ka)

 

 

 

 

j (ka)

 

 

 

Cn =

 

 

 

n

 

1

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

.

(8.92)

 

 

ρ c

 

j

(k a)

h

(1)(ka) h(1)(ka)

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1 1

 

n

 

 

1

 

 

 

n

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

ρc

j (k a) h

(1)(ka)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

1

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда коэффициенты Dn определяются из уравнения (8.90) или (8.91). Итак, решение найдено. Он дает возможность вычислить звуковое поле снаружи и внутри сферы.

Как пример, на рис. 8.10 показано распределение нормированной (к амплитуде падающей волны) амплитуды давления снаружи и внут-

523

ри сферы. Окружающей средой является вода (ρ = 103 кг/м3, с = 1,5

103 м/с), параметры сферы: ρ1 = 1,1 103 кг/м3, с1 = 0,95 103 м/с,

а/λ = 0,9 (λ длина волны в среде). Графики иллюстрируют распре- деление амплитуды давления вдоль лучей θ = 180° (z/a < 0) и θ = 0 (z/a > 0). На графике можно видеть достаточно резкий максимум давления, которое отвечает белому пятну на карте распределения ам- плитуды давления. Можно сказать, при выбранных параметрах мы имеем жидкую сферическую линзу.

Рис. 8.10. Диаграммы распределения (справа) и графики (слева) для норми- рованной амплитуды давления снаружи сферы и внутри сферы:

ρ = 103 кг/м3, с = 1,5 103 м/с, ρ1 = 1,1 103 кг/м3, с1 = 0,95 103 м/с, а/λ = 0,9

Рассмотрим интересный случай сферы малого волнового размера, т.е. ka << 1 и k1a << 1. Для этого согласно формулам (8.78)—(8.81) найдем асимптотику коэффициентов Cn при ka << 1 и k1a << 1. Пусть n = 0:

j0(k1a) = − j0(k1a) ≈ − 3 , j0(k1a) j1(k1a) k1a

h0(1)(ka) = −h1(1)(ka) = − 1 , h0(1)(ka) h0(1)(ka) ka

тогда

 

(ka)

 

 

 

 

 

j1(ka)

 

 

 

 

 

j0

= −

 

= −ka

,

 

j0

(ka)

j(ka)

 

 

 

 

 

 

3

 

 

j0(ka)

 

 

1

 

 

= ika,

 

h(1)(ka)

 

i /(ka )

0

 

 

 

 

 

 

 

 

 

 

524

 

 

ρ c

3

ka

 

 

ρ c2

 

χ

 

 

 

 

 

 

 

1 1

 

 

 

 

 

 

 

1

1 1

 

 

 

 

 

1 ρc2

1

1

 

 

 

ρc

k a

3

 

 

χ

 

3

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

C0 =

 

 

 

 

ika = i

 

 

ka = i

 

 

 

(ka)

, (8.93)

 

ρ1c1

3

1

 

ρ1c1 ka 1

2

 

χ1

1

 

 

 

1 3

ρc k1a

 

 

 

(ka)

3

 

 

 

ρc

k1a

 

ka

 

(ka)2

χ

 

 

где χ = ρ c2

упругость вещества сферы, χ = ρc2 упругость среды

1

1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вокруг сферы. Пусть n = 1:

 

 

 

j1(k1a)

 

=

 

 

 

 

j1(k1a)

 

 

 

 

 

 

=

 

 

 

 

 

1

 

 

 

 

 

(k a),

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j1(k1a)

 

 

 

j1(k1a)

j2(k1a)

 

 

1

 

 

 

(k a)2 3

1

 

 

 

 

 

 

 

 

 

 

 

 

 

k a

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(k1a)

 

3

5(k1a)

 

 

 

 

(ka)

 

 

 

1

 

 

h(1)(ka)

 

 

 

 

1

 

 

 

h(1)(ka)

 

 

1

 

 

2

2

 

 

j0

 

 

,

 

1

 

=

 

 

 

 

2

 

 

 

=

 

3 (ka)

= −

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j0(ka)

 

ka

h(1)(ka)

ka

h

(1)(k a)

ka

ka

 

 

 

 

 

 

 

 

 

 

 

 

(ka)3 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j (ka)

 

 

 

 

ka (ka)2

 

 

(ka)3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

=

 

3 (i)

 

= i

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h1(1)(ka)

 

 

 

3

 

 

 

 

 

 

 

 

 

тогда

 

 

 

 

 

 

 

 

 

 

1 ρ1c1

k1a

 

 

 

 

 

 

 

 

 

 

ρ1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ka)3

 

 

 

1

 

 

(ka)3

 

 

 

 

 

 

 

 

 

C

 

=

 

 

 

 

pc

 

 

ka

 

i

= i

 

 

 

ρ

 

 

.

 

 

(8.94)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1 +

ρ1c1 k a

2

 

 

 

3

 

 

 

 

1 + 2 ρ1

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pc

1

 

 

ka

 

 

 

 

 

 

 

 

 

ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть n > 1, проводя аналогичные выкладки, нетрудно убедиться, что

 

j

(k a)

 

j (ka)

 

j

(k a)

h(1)(ka)

 

 

 

k a

 

 

n

1

 

 

n

 

и

 

n

1

 

 

n

 

 

 

имеют величину порядка

1

, то-

 

j (k a)

jn (ka)

j

(k a)

h(1)(ka)

 

 

 

 

 

 

 

 

ka

 

 

n

1

 

 

 

 

 

n

1

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jn (ka)

 

 

2n +1

 

 

 

гда как отношение

 

 

 

(ka)

 

. Итак, именно это отношение

h(1)(ka)

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

будет определять порядок величины коэффициентов Cn. Действи- тельно, если коэффициенты C0 и C1 имеют порядок (ka)3, то C2 уже имеет порядок (ka)5.

Таким образом, в случае произвольного препятствия малого вол- нового размера в ряде (8.71) для рассеянной волны ps достаточно ос- тавить лишь первые два члена, т.е.

ps = −(C0h0(1)(kr ) + 3iC1h1(1)(kr )cos θ).

(8.95)

В дальнем поле, где kr >> 1, согласно асимптотическим формулам (8.74) выражение (8.95) будет иметь вид

525

p

=

 

i

[C

0

+ 3C

cos θ]exp(ikr ).

(8.96)

 

 

s

 

kr

1

 

 

 

 

 

 

 

 

Таким образом, произвольная сфера малого волнового размера созда- ет одновременно и монопольное, и дипольное рассеяние. Другими со- ставляющими рассеянного поля можно пренебрегать. Вообще приве- денный результат действителен для малого препятствия произвольной формы.

Подставляя в (8.96) асимптотические формулы для коэффициен- тов C0 и C1 малой сферы, получаем следующую формулу для рассеян- ного поля:

 

 

exp(ikr )

1 − χ

 

/χ

 

 

1 − ρ /ρ

 

3

 

 

 

 

 

1

 

 

 

 

1

 

 

p

= −

 

 

 

 

 

+

 

 

 

cos θ

(ka) .

(8.97)

kr

(ka)2

3χ /χ

 

 

+ 2ρ1 /ρ

s

 

 

 

1

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

Приняв во внимание, что обычно (ka)2 << 3χ1 /χ,, (8.97) можно пере- писать в виде

 

exp(ikr ) 1 − χ1 /χ

 

1 − ρ1 /ρ

 

3

 

ps =

 

 

3χ /χ

 

cos θ

(ka) .

(8.98)

kr

1 + 2ρ /ρ

 

 

 

1

 

1

 

 

 

Неравенство (ka)2 << 3χ1 /χ не выполняется только для сферы, кото- рая подобна по своим свойствам мягкой сфере, у такой сферы значе- ние упругости χ1 очень мало. Формула (8.98) впервые была получена Релеем .

Используя формулу (8.76) в случае рассеянного поля в виде (8.97) и (8.98), получаем следующие значения для полного сечения рассеяния:

 

 

 

 

 

 

1 − χ1 /χ

 

 

 

2

 

1 1 − ρ1 /ρ

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

4

 

 

σ = 4πa

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

(ka)

,

(8.99)

 

 

 

2

 

 

 

3 1+

2ρ /ρ

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ka) 3χ /χ

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

1

 

 

1 − ρ1 /ρ

 

2

 

 

 

 

 

 

σs = 4πa

2

 

1 − χ1 /χ

 

+

 

 

 

 

 

4

 

(8.100)

 

 

 

3χ /χ

 

 

3

 

1+ 2ρ

 

/ρ

 

 

(ka) .

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

Итак, основными свойствами рассеяния звука малой сферой (ka<< 1, k1a << 1; и вообще препятствиями малых волновых разме- ров) являются:

• наличие монопольной и дипольной составляющих в рассеянном поле; если препятствие отличается от среды только сжимаемостью

Рэлей (Rayleigh) Стретт Джон Уильям (1842—1919)английский физик, лауреат Нобелевской премии (1904).

526

(ρ1 = ρ), то имеем рассеяние монопольного типа, а если только плотно- стью (χ1 = χ), то дипольного типа;

• сечение рассеяния пропорционально частоте в четвертой степе- ни, что обуславливает увеличение рассеивающей способности пре- пятствия с ростом частоты. Такую особенность рассеяния звука на малых препятствиях называют рэлеевским рассеянием (в честь Рэлея, который открыл этот закон в 1898 г.).

Приведенные соотношения позволяют записать полученные преж- де формулы для полного сечения рассеяния жесткой и мягкой сфер. Действительно, в случае малой жесткой неподвижной сферы имеем χ1 / χ → ∞ и ρ1 /ρ → ∞; поэтому из формулы (8.100) получаем значение σs, которое согласуется с (8.85):

σs = 4πa

2

 

1

+

1

 

1

4

=

7

πa

2

4

;

(8.101)

 

 

9

3

4

 

(ka)

9

 

(ka)

 

 

 

 

 

 

 

 

 

 

 

 

 

в случае малой мягкой сферы имеем χ1 / χ 0 и ρ1 /ρ 0, поэтому из формулы (8.99) получаем значение σs, которое совпадает с (8.83):

 

 

 

1

+ 1

 

 

 

σ

= 4πa2

 

(ka)4

4πa2.

(8.102)

4

s

 

 

3

 

 

 

 

 

 

(ka)

 

 

 

Подставляя в формулу (8.96) асимптотические значения коэффициентов С0 и С1 при ka << 1, нетрудно убедиться, что малая мягкая сфера имеет ненаправленное рассеяние (монопольный тип), а малая жесткая неподвижная сфера создает направленное рассеянное

поле, которое определяется соотношением:

 

 

p

 

 

 

1 3 cos θ

 

.

(8.103)

 

 

 

 

 

 

 

s

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

В завершение приведем следующие соображения относительно частотной зависимости сечения рассеяния в формулах (8.99) и (8.100), из которых следует, что с ростом частоты имеет место суще- ственное увеличение рассеивающей способности малой сферы. Релей в 1898 г. в своей фундаментальной работе [50, т. 2, с.153] писал: “Ес- ли первичный звук является сложной музыкальной нотой, то разные ее составляющие тона рассеиваются в неодинаковой пропорции. Ок- тава, например, в шестнадцать раз сильнее по сравнению с основ- ным тоном во вторичном звуке, чем это было в первичном. Нетрудно, таким образом, понять, каким образом эхо, отраженное от такого препятствия, как группа деревьев, может стать повышенным на ок- таву”. (Под октавой понимается тон с частотой в два раза большей, чем частота основного тона.)

527

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

Отношение длины волны красного к длине волны синего света при- близительно составляет 1,44. Возводя это число в четвертую степень, получаем приблизительно 4,3. Таким образом, согласно закону Рэлея интенсивность синего света, рассеянного в атмосфере, в четыре раза превышает интенсивность красного света. И потому в слое воздуха высотой в десятки километров в значительной мере доминируют си- ний и голубой цвета.

Мы видим небо синим, а не фиолетовым, хотя закон Релея прогно- зирует доминирование фиолетового цвета. Это расхождение обуслов- лено двумя причинами. Первая, в спектре солнечного света фиолето- вых лучей значительно меньше, чем синих. Вторая, чувствительность глаза человека зависит от длины волны света; на фиолетовые лучи глаз реагирует значительно слабее, чем на сине-зеленые.

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

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

Предлагаем читателям провести такой опыт [24, с. 339]. Налейте в стеклянную посуду воду и направьте на сосуд сбоку луч света от кар- манного фонаря. Посмотрите на свет, который рассеян под углом

90 , и на фонарь сквозь сосуд с водой. Прибавьте к воде несколько капель молока и перемешайте. Продолжайте смотреть, постепенно добавляя молоко. Обратите внимание на голубой оттенок рассеянного света и на желтоватый или красный оттенок у света, который прошел сквозь сосуд. Объясните это явление. Если добавить еще молоко, то го-

528

лубой оттенок рассеянного света исчезнет. Рассеянный свет будет поч- ти белым, а воображаемый закатбудет становиться все более крас- ным. Объясните это. В конце концов, вода в сосуде станет настолько мутной, что мы не сможем видеть сквозь ее толщу лампочку фонаря; рассеянный свет станет белым, и мы не будем видеть луч света. “Воз- духпревратится в белую тучу”.

Рассеяние света и звука, в сущности, происходит по одинаковым законам. Хотя то, что длина световой волны в видимом диапазоне час- тот имеет порядок 5 10–5 см, обуславливает значительную разницу в рассеянии света и звука в воде. В отличие от воздуха, где флуктуации концентраций молекул не только возможны, а и постоянно возникают благодаря тепловому движению молекул, молекулы воды касаютсяи потому размещены в среде довольно равномерно. Рассеяние волн про- исходит на маленьких твердых частицах, которые находятся в естест- венной водной среде. Обычно эти частицы значительно превышают по размерам длину световой волны, поэтому почти для всех твердых час- тиц в море сечение рассеяния для света определяется их геометриче- ским сечением (формула (8.86)). Те самые частицы оказываются зна- чительно меньшими, чем длина звуковой волны в море (обычно 1см и больше), и потому звук на них испытает рэлеевское рассеяние (форму- ла (8.99)). Этим объясняется тот факт, что для света морская вода мут- ная, а для звука прозрачная.

8.9. Рассеяние звука газовым пузырьком в жидкости

Рассмотрим рассеяние на препятствиях, которые, несмот- ря на малые размеры по сравнению с длиной звуковой волны (ka << 1), при соответствующих частотах звуковой волны имеют вы- сокую эффективность рассеяния. Еще раз подчеркнем, что речь идет именно о малом препятствии. Природа этого аномального рассеяния связана с резонансными явлениями. Интересным и важным для практики примером такого препятствия является газовый пузырек в жидкости. Пузырек - препятствие, которое имеет значительно мень- шие значения плотности и упругости, чем жидкость, которая ее ок- ружает.

В начале приведем следующие физические соображения. Анали- зируя колебания пузырька, следует отметить, что в колебаниях участ- вует не только газо-паровая смесь, заполняющая объем пузырька, но и часть жидкости, которая окружает пузырек и движется с ним практически синфазно. Эта жидкость образует присоединенную мас- су. Газо-паровая смесь и поверхностное натяжение образуют упру- гость системы, а потери энергии на излучение и на теплообмен с ок-

529

ружающей средой эквивалентны присутствию трения в системе. По- этому, с точки зрения теории механических колебаний, пузырек эк- вивалентен массе на пружине при наличии трения, т. е. закономер- ности колебаний для этих двух случаев одинаковы. Известно, что в такой колебательной системе имеет место резонанс. Амплитуда коле- баний пузырька будет максимальной на той частоте падающей гар- монической волны, которая совпадает с частотой резонанса. На этой частоте компенсирующее движение будет наибольшим и также наи- большим будет рассеяние.

Теперь перейдем к количественным соотношениям. Пусть на газо- вый пузырек радиусом а, который находится в жидкости с плотно- стью ρ и со скоростью звука с, падает плоская гармоническая волна с частотой ω (рис. 8.8). Частота звука такая, что выполняется условие ka << 1. Будем считать, пренебрегая при этом вязкостью и теплопро- водностью, что газ внутри пузырька описывается линеаризованным адиабатическим уравнением состояния (4.16):

pп = cп2ρп,

(8.104)

где pп звуковое давление; ρп переменная плотность; cп

ско-

рость звука для газа, который заполняет пузырек.

Выразим давление pп через колебательную скорость υп нормаль- ных смещений поверхности пузырька. Поскольку

ρп

= −

dV

, то

ρп = −ρп0

dV

,

(8.105)

ρп0

V

V

 

 

 

 

 

где ρп0 средняя плотность газа; V объем пузырька; dV изме- нение объема, которое соответствует ρп . Перепишем (8.105), учиты-

вая формулы для объема сферы V = 4πa3 /3 и, соответственно, dV = 4πa2da:

ρп = −3ρп0 da .

(8.106)

a

 

Подставив (8.106) в (8.104) и продифференцировав это уравнение по времени, с учетом того, что dpп /dt = −iωpп, а da /dt = υп является ко-

лебательной скоростью нормальных смещений поверхности пузырька, получим искомое соотношение:

 

 

 

3ρ

c2

 

3χ

п

 

 

p

п

= −i

 

п0

п

υ = −i

 

υ ,

(8.107)

ωa

 

ωa

 

 

 

п

п

 

где χп = ρп0cп2 упругость газа в пузырьке.

530