𝑑ψ
𝑑𝑟
+ψ²+ψ
⎧
⎪
⎩
1
𝑟
+
2
∂𝑝
0
𝑝
0
∂𝑟
-
2
∂𝑊
𝑊
∂𝑟
⎫
⎪
⎭
-
2
∂𝑝
0
∂𝑊
𝑝
0
𝑊
∂𝑟
∂𝑟
=0.
(9)
Положим теперь, что 𝑊=𝑐+σ, где 𝑐 — средняя скорость струи, а σ — величина, малая по сравнению с этой скоростью. В этом случае ψ мало, и, пренебрегая членами порядка (σ/𝑐)², мы находим из (9)
𝑑ψ
𝑑𝑟
+ψ
⎧
⎪
⎩
1
𝑟
+
2
∂𝑝
0
𝑝
0
∂𝑟
⎫
⎪
⎭
-
2
∂𝑝
0
𝑑σ
𝑐𝑝
0
∂𝑟
𝑑𝑟
=0.
(10)
Численное значение |𝑖𝑘𝑟|, соответствующее условиям эксперимента, будет очень малым (длина волны велика по сравнению с диаметром струи). Поэтому, чтобы не усложнять формулы, мы будем при вычислении ψ использовать только первый член разложения 𝐽𝑛(𝑖𝑘𝑟). При этом (1/𝑝0)∂𝑝0/∂𝑟 и решение уравнения (9) принимает вид
ψ=
2𝑛
𝑐
𝑟
-(2𝑛+1)
⎧
⎪
⎩
𝑟
∫
0
𝑑σ
𝑑𝑟
𝑟
2𝑛
𝑑𝑟
+𝐶
⎫
⎪
⎭
.
Учитывая конечность ψ при 𝑟=0, находим, что 𝐶=0. Интегрируя по частям, получаем
ψ=
2𝑛
𝑐𝑟
σ-
4𝑛2
𝑐𝑟2𝑛+1
𝑟
∫
0
σ𝑟
2𝑛-1
𝑑𝑟
.
(11)
Предположим, что уравнение поверхности имеет вид
𝑟-𝑎=ζ=
𝐵
𝑒
𝑖𝑛θ+𝑖𝑘𝑧
.
Общие граничные условия дают при этом
𝐷
𝐷𝑡
(𝑟-𝑎-ζ)
=
⎧
⎪
⎩
α
∂
∂𝑟
+β
∂
∂θ
+𝑤
∂
∂𝑧
⎫
⎪
⎭
(𝑟-𝑎-ζ)
=0,
откуда мы получаем, пренебрегая величинами того же порядка, что и в случае уравнений (3)
α-𝑊
∂ζ
∂𝑧
=0,
ζ=
𝑖
𝑊𝑘
α.
Подобным же образом, обозначая через 𝑅1 и 𝑅2 главные радиусы кривизны, мы имеем далее
1
𝑅1
+
1
𝑅2
=
1
𝑎
-
ζ
𝑎²
-
1
𝑎²
∂²ζ
∂θ²
-
∂²ζ
∂𝑧²
=
1
𝑎
-α
𝑖(𝑛²-1+𝑘²𝑎²)
𝑎²𝑊𝑘
.
Обозначая коэффициент поверхностного натяжения через 𝑇, мы находим следующее динамическое условие на поверхности:
𝑇
⎧
⎪
⎩
1
𝑅1
+
1
𝑅2
⎫
⎪
⎭
-𝑝
=const.
Отсюда находим (в тех же приближениях, что и раньше), используя формулу (4):
⎡
⎢
⎣
𝑇(𝑛²-1-𝑘²𝑎²)
ρ𝑎²𝑘²𝑊²
∂𝑝
∂𝑟
-𝑝
⎤
⎥
⎦𝑟=𝑎
=0.
(12)
Из (12) с помощью (7) и (8) получаем
𝑘²
=
𝑇
𝑖𝑎𝑘
𝐽
'
𝑛
(𝑖𝑎𝑘)
ρ𝑎³ 𝐽
𝑛 (𝑖𝑎𝑘)
(𝑛²-1-𝑘²𝑎²)
⎡
⎢
⎣
1
𝑊²
⎧
⎪
⎩
1+ψ
𝑝0
∂𝑝0/∂𝑟
⎫
⎪
⎭
⎤
⎥
⎦𝑟=𝑎
.
(13)
Наконец, из формулы (13) с помощью (11) и в том же приближении, которое было использовано при вычислении ψ, имеем
𝑘²
=
𝑇
𝑖𝑎𝑘
𝐽
'
𝑛
(𝑖𝑎𝑘)
ρ𝑐²𝑎³ 𝐽
𝑛 (𝑖𝑎𝑘)
(𝑛²-1-𝑘²𝑎²)
⎡
⎢
⎣
1-
4𝑛
𝑐𝑎2𝑛
𝑎
∫
0
σ𝑟
2𝑛-1
𝑑𝑟
⎤
⎥
⎦
.
(14)
Это соотношение, за исключением последнего члена, совпадает с решением, полученным Рэлеем. Мы видим, таким образом, что влияние различия в величинах скоростей концентрических слоёв струи состоит в том, что в формуле для длины волны λ=2π/𝑘 выражение для средней скорости струи 𝑐 заменяется «эффективной средней скоростью»
𝑐'=𝑐+
2𝑛
𝑎2𝑛
𝑎
∫
0
σ𝑟
2𝑛-1
𝑑𝑟
(15)
Из выражения (15) мы видим, что чем больше и, тем ближе будет эффективная средняя скорость совпадать со скоростью движения поверхности. Это объясняется тем, что чем больше значение 𝑛 (т. е. число волн, укладывающееся на длине окружности сечения струи), тем быстрее будет уменьшаться скорость колебательного движения частиц жидкости при переходе от поверхности к оси струи.
Теперь можно показать, что 𝑐' будет меньше 𝑐, если скорость струи имеет наибольшее значение в центре и постепенно уменьшается к поверхности (это и будет иметь место в условиях эксперимента). Поскольку 𝑐 — средняя скорость струи, то
𝑎
∫
0
σ𝑟
𝑑𝑟
=0,
и в рассматриваемом случае
𝑎
∫
0
σ𝑟
𝑑𝑟
>0,
если 𝑎 > 𝑟 > 0. Отсюда получаем для 𝑛 ≥ 2 (в экспериментах 𝑛 = 2)
𝑐'-𝑐
=
2𝑛
𝑎2𝑛
𝑎
∫
0
σ𝑟
2𝑛-1
𝑑𝑟
=-
2𝑛(2𝑛-2)
𝑎2𝑛
𝑎
