Выбрать главу

⎦𝑟=𝑎

=0,

(26)

1

∂α

𝑟

∂θ

+

∂β

∂𝑟

-

β

𝑟

⎭𝑟=𝑎

=0,

∂α

∂𝑧

+

∂𝑤

∂𝑟

⎭𝑟=𝑎

=0.

(27)

Подставляя в эти условия значения 𝑝, α, β и 𝑤, задаваемые формулами (12) и (21), и исключая 𝐵/𝐴 и 𝐶/𝐴, получаем уравнения для определения 𝑏. Поскольку вычисления оказываются довольно длинными и приводят к очень громоздкому результату, мы не будем воспроизводить указанную процедуру точно, а ограничимся приближением, достаточным для наших целей.

В экспериментах численное значение величины | 𝑖𝑎𝑏 | оказывается малым, так как длина волны обычно велика по сравнению с диаметром струи: значение же величины | 𝑖𝑎𝑑 |, напротив, велико, так как мал коэффициент вязкости. (Во всех экспериментах | 𝑖𝑎𝑏 | < 0,24 и | 𝑖𝑎𝑑 | > 20.)

При всех значениях 𝑥 справедливо разложение

𝐽

𝑛

(𝑥)

=

𝑥𝑛

2𝑛⋅𝑛!

-

𝑥𝑛+2

2𝑛+2⋅1!(𝑛+1)!

+

𝑥𝑛+4

2𝑛+4⋅2!(𝑛+2)!

-…

(28)

Ряд (28) быстро сходится при малых 𝑥, но очень медленно — при больших 𝑥. Из разложения (28) следует

𝐽

'

𝑛

(𝑥)

=

𝑛

𝑥

𝐽

𝑛

(𝑥)

1-

𝑥2

2𝑛(𝑛+1)

-

𝑥4

23⋅𝑛(𝑛+1)2(𝑛+2)

-…

и далее с помощью (19)

𝐽

''

𝑛

(𝑥)

=

𝑛(𝑛-1)

𝑥²

𝐽

𝑛

(𝑥)

1-

𝑥²(2𝑛+1)

2(𝑛-1)𝑛(𝑛+1)

+

𝑥4

23(𝑛-1)𝑛(𝑛+1)2(𝑛+2)

.

Поэтому, вычисляя диссипативные члены в уравнении для определения 𝑏, мы будем полагать

𝐽

'

𝑛

(𝑖𝑎𝑏)

=-

𝑖𝑛

𝑎𝑏

𝐽

𝑛

(𝑖𝑎𝑏)

1+

𝑎²𝑏²

2𝑛(𝑛+1)

и

𝐽

''

𝑛

(𝑖𝑎𝑏)

=-

𝑛(𝑛-1)

𝑎²𝑏²

𝐽

𝑛

(𝑖𝑎𝑏)

1+

𝑎²𝑏²(2𝑛+1)

2(𝑛-1)𝑛(𝑛+1)

.

(29)

Для вычисления 𝐽𝑛(𝑥) при больших значениях 𝑥 мы воспользуемся асимптотическим выражением

𝐽

𝑛

(𝑥)

(2π𝑥)

[

𝑃

𝑛

(𝑥)

+

𝑖𝑄

𝑛

(𝑥)

]

exp 𝑖

𝑥-

2𝑛+1

4

π

[

𝑃

𝑛

(𝑥)

-

𝑖𝑄

𝑛

(𝑥)

]

exp -𝑖

𝑥-

2𝑛+1

4

π

(30)

где

𝑃

𝑛

(𝑥)

=

1-

(4𝑛²-1²)(4𝑛²-3²)

2!(8𝑥)²

+

(4𝑛²-1²)(4𝑛²-3²)(4𝑛²-5²)(4𝑛²-7²)

4!(8𝑥)4

-…

и

𝑄

𝑛

(𝑥)

=

4𝑛²-1²

1!8𝑥

-

(4𝑛²-1²)(4𝑛²-3²)(4𝑛²-5²)

3!(8𝑥)3

+…

Если в формуле (30), справедливой при положительной вещественной части 𝑥, взять несколько членов, то она будет давать очень хорошее приближение для 𝐽𝑛(𝑥) при больших 𝑥. При использовании формулы (30) 𝑥 можно записать в виде 𝑎-𝑖𝑏, где 𝑎 и 𝑏 — большие положительные числа. При этом член с 𝑒𝑖𝑥 будет преобладающим; пренебрегая членом с 𝑒-𝑖𝑥, имеем

𝐽

'

𝑛

(𝑥)

=

𝑖𝐽

𝑛

(𝑥)

1+

𝑖

2𝑥

-

4𝑛²-1

8𝑥²

+

𝑖(4𝑛²-1)

8𝑥³

-…

и, согласно (19),

𝐽

''

𝑛

(𝑥)

=

-𝐽

𝑛

(𝑥)

1+

𝑖

𝑥

-

2𝑛²-1

2𝑥²

-

𝑖(4𝑛²-1)

8𝑥³

-…

.

Исходя из сказанного, мы в дальнейших расчётах положим

𝐽

'

𝑛

(𝑖𝑎𝑑)

=

𝑖𝐽

𝑛

(𝑖𝑎𝑑)

1+

1

2𝑎𝑑

+

4𝑛²-1

8𝑎²𝑑

и

𝐽

''

𝑛

(𝑖𝑎𝑑)

=

𝐽

𝑛

(𝑖𝑎𝑑)

1+

1

𝑎𝑑

+

2𝑛²+1

2𝑎²𝑑²

(31)

Теперь из соотношений (27) с помощью (29) и (31) мы получаем

𝐴

1

𝑐ρ

𝐽

𝑛

(𝑖𝑎𝑏)

2𝑛(𝑛-1)

𝑎²𝑏

1+

𝑎²𝑏²

2(𝑛-1)(𝑛+1)

+

+

𝐵𝐽

𝑛

(𝑖𝑎𝑑)

2𝑛𝑏

𝑎𝑑

1+

3

2𝑎𝑑

+

4𝑛²-1

8𝑎²𝑑²

-

-

𝐶

𝑛

𝐽

𝑛

(𝑖𝑎𝑑)

𝑖𝑑²

𝑛

1+

2

𝑎𝑑

+

2𝑛²+1

𝑎²𝑑²

=0

(32)

и

𝐴

1

𝑐ρ

𝐽

𝑛

(𝑖𝑎𝑏)

2𝑛

𝑎

1+

𝑎²𝑏²

2𝑛(𝑛+1)

+

+

𝐵𝐽

𝑛

(𝑖𝑎𝑑)

𝑑

1+

1

2𝑎𝑑

+

4𝑛²-1

8𝑎²𝑑²

-

𝐶𝐽

𝑛

(𝑖𝑎𝑑)

𝑖𝑏

𝑎

=0.

(33)

Из соотношений (32) и (33) находим

𝐵𝐽

𝑛

(𝑖𝑎𝑑)

𝐴

1

𝑐ρ

𝐽

𝑛

(𝑖𝑎𝑏)

2𝑛

𝑎𝑑

1+

𝑎²𝑏²

2𝑛(𝑛+1)

1-

1

2𝑎𝑑

-

12𝑛²-8𝑛-3

8𝑎²𝑑²

и

𝐶𝐽

𝑛

(𝑖𝑎𝑑)

𝐴

1

𝑐ρ

𝐽

𝑛

(𝑖𝑎𝑏)

𝑖2𝑛²(𝑛-1)

𝑎²𝑑²𝑏

1+

𝑎²𝑏²

2(𝑛²-1)

1-

2

𝑎𝑑

-

2𝑛²-3

𝑎²𝑑²

.

(34)

Формула (26) с учётом равенств (12), (21), (29), (31), (34) и (13) теперь даёт

𝑏²-𝑖𝑏

μ

ρ

4𝑛(𝑛-1)

𝑎²𝑐

1+

𝑎²𝑏²