Расчет акустического поля с учетом течений
В данной статье решена задача о распространении звука в неоднородной движущейся среде в виде рядов по степеням ε. В этом случае нулевое приближение соответствует задаче о распространении звука в неподвижной неоднородной среде, а поправка, обусловленная наличием течения, выражается через нулевое приб...
Gespeichert in:
| Veröffentlicht in: | Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану) |
|---|---|
| Datum: | 2008 |
| Hauptverfasser: | , , , |
| Format: | Artikel |
| Sprache: | Russisch |
| Veröffentlicht: |
Науково-технічний центр панорамних акустичних систем НАН України
2008
|
| Online Zugang: | https://nasplib.isofts.kiev.ua/handle/123456789/19217 |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Zitieren: | Расчет акустического поля с учетом течений / А.И. Гончар, Ю.А. Гончар, О.С. Голод, Н.С. Григорьева // Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану): Зб. наук. пр. — Запоріжжя: НТЦ ПАС НАН України, 2008. — № 5. — С. 34-72. — Бібліогр.: 10 назв. — рос. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraine| _version_ | 1859639474909085696 |
|---|---|
| author | Гончар, А.И. Гончар, Ю.А. Голод, О.С. Григорьева, Н.С. |
| author_facet | Гончар, А.И. Гончар, Ю.А. Голод, О.С. Григорьева, Н.С. |
| citation_txt | Расчет акустического поля с учетом течений / А.И. Гончар, Ю.А. Гончар, О.С. Голод, Н.С. Григорьева // Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану): Зб. наук. пр. — Запоріжжя: НТЦ ПАС НАН України, 2008. — № 5. — С. 34-72. — Бібліогр.: 10 назв. — рос. |
| collection | DSpace DC |
| container_title | Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану) |
| description | В данной статье решена задача о распространении звука в неоднородной движущейся среде в виде рядов по степеням ε. В этом случае нулевое приближение соответствует задаче о распространении звука в неподвижной неоднородной среде, а поправка, обусловленная наличием течения, выражается через нулевое приближение. Так как отношение характерного размера изменения свойств среды по вертикали к характерному размеру изменения свойств среды в горизонтальном направлении имеет порядок 10-3, то с помощью пространственно-временного лучевого метода решена задача по нахождению коэффициентов ряда по степеням малого параметра ε.
У даній статті вирішено задачу про поширення звуку в неоднорідному середовищі, що рухається, у вигляді рядів по ступенях ε. У цьому випадку нульове наближення відповідає завданню про поширення звуку в нерухливому неоднорідному середовищі, а виправлення, обумовлене наявністю течії, виражається через нульове наближення. Оскыльки відношення характерного розміру зміни властивостей середовища по вертикалі до характерного розміру зміни властивостей середовища в горизонтальному напрямку має порядок 10-3, то за допомогою просторово-тимчасового променевого методу вирішено задачу по знаходженню коефіцієнтів ряду по ступенях малого параметра ε.
This article solves the problem about the propagation of sound into nonhomogeneous moving medium in series in terms of power of ε. In this case zeros approximation corresponds to the problem of sound propagation into nonmoving nonhomogenenous medium and the allowance conditioned by availability of current is expressed by zeros approximation. Since the relation of characteristic dimension of change of medium properties in vertical position to characteristic dimension of change of medium properties in horizontal position has the order 10-3 the problem for estimation of ratios of series in terms of powers of small parameter ε was solved by the spatio-temporal ray-path method.
|
| first_indexed | 2025-12-07T13:19:52Z |
| format | Article |
| fulltext |
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
34
УДК 551.463.22
РАСЧЕТ АКУСТИЧЕСКОГО ПОЛЯ С УЧЕТОМ ТЕЧЕНИЙ
© А.И. Гончар, Ю.А. Гончар, О.С. Голод, Н.С. Григорьева, 2008
Государственный Северо-Западный технический университет, г. Санкт-Петербург
Научно-технический центр панорамных акустических систем НАН Украины, г. Запорожье
У даній статті вирішено задачу про поширення звуку в неоднорідному середовищі, що рухається, у
вигляді рядів по ступенях ε. У цьому випадку нульове наближення відповідає завданню про поширення звуку в
нерухливому неоднорідному середовищі, а виправлення, обумовлене наявністю течії, виражається через
нульове наближення. Оскыльки відношення характерного розміру зміни властивостей середовища по вертикалі
до характерного розміру зміни властивостей середовища в горизонтальному напрямку має порядок 10-3, то за
допомогою просторово-тимчасового променевого методу вирішено задачу по знаходженню коефіцієнтів ряду
по ступенях малого параметра ε.
В данной статье решена задача о распространении звука в неоднородной движущейся среде в виде
рядов по степеням ε. В этом случае нулевое приближение соответствует задаче о распространении звука в
неподвижной неоднородной среде, а поправка, обусловленная наличием течения, выражается через нулевое
приближение. Так как отношение характерного размера изменения свойств среды по вертикали к характерному
размеру изменения свойств среды в горизонтальном направлении имеет порядок 10-3, то с помощью
пространственно-временного лучевого метода решена задача по нахождению коэффициентов ряда по степеням
малого параметра ε.
This article solves the problem about the propagation of sound into nonhomogeneous moving medium in series
in terms of power of ε. In this case zeros approximation corresponds to the problem of sound propagation into
nonmoving nonhomogenenous medium and the allowance conditioned by availability of current is expressed by zeros
approximation. Since the relation of characteristic dimension of change of medium properties in vertical position to
characteristic dimension of change of medium properties in horizontal position has the order 10-3 the problem for
estimation of ratios of series in terms of powers of small parameter ε was solved by the spatio-temporal ray-path
method.
АКУСТИЧЕСКОЕ ПОЛЕ, СКОРОСТЬ ЗВУКА, РАСПРОСТРАНЕНИЕ ЗВУКА,
НЕОДНОРОДНАЯ ДВИЖУЩАЯСЯ СРЕДА, УРАВНЕНИЯ АКУСТИКИ
Акустические поля течений могут быть как в свободном пространстве, так и вблизи
разного рода препятствий, где образуются колебания частиц среды или условия, что
вызывают течения.
В условиях реального океана все основные характеристики среды (давление,
плотность, температура, скорость движения жидкости) – функции координат и времени
(диапазон изменения составляет, в частности: для температуры – от 4 до 40 0С, давления – от
0 до 1000 бар, плотности – от 1,00 до 1,04 г/см3, скорости движения жидкости – от 0 до
2,5 м/с, скорости распространения звука – от 1,45·103
до 1,60·103 м/с). При этом, так как
отношение ϑmax/cmin не превосходит 2·10-3, где ϑmax – максимальное значение скорости
течения, сmin – минимальное значение скорости звука, то при моделировании
распространения акустических волн в неоднородном океане с течениями всегда имеется
малый параметр ε.
Наша цель состоит в том, чтобы представить решение задачи о распространении звука
в неоднородной движущейся среде в виде рядов по степеням ε.
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
35
В этом случае нулевое приближение будет соответствовать задаче о распространении
звука в неподвижной неоднородной среде, а поправка, обусловленная наличием течения,
выразится через нулевое приближение.
Так как отношение характерного размера изменения свойств среды по вертикали к
характерному размеру изменения свойств среды в горизонтальном направлении имеет
порядок 10-2, то с помощью пространственно-временного лучевого метода [7, 8] решение
задач по нахождению коэффициентов ряда по степеням малого параметра ε удается
существенно упростить.
Как показывают экспериментальные исследования [1], наличие течений существенно
влияет на распространение звука в океане. Необходимость учета движения среды
подтверждается и теоретическими работами [2 - 6], в которых рассматривается
распространение звука в канале с горизонтальными граничными поверхностями. В [2] в
пределах канала скорость звука постоянна, а скорость течения ϑ
r
имеет постоянное
направление и величину: ;x yi jϑ ϑ ϑ= +
r r r
ϑx, ϑy=const. В [3] скорость звука также
предполагается постоянной, однако скорость течения линейно меняется с глубиной:
( ) ( )x yz i z jϑ ϑ ϑ= +
r r r
; ϑx(z), ϑy(z) – линейные функции z. В работах [4] и [5] скорость звука
линейно зависит от глубины, причем в [4] скорость течения имеет постоянное направление и
величину, а в [5]: ( ) 0
( ) ; ,z j z zϑ ϑ ϑ ϑ ν= = +
r r
где ϑ0 – скорость течения у источника, ν –
градиент скорости течения. Наконец, в [6] для лучевого представления поля получены
поправки к оптическим длинам лучей и амплитудам поля на них для произвольного профиля
скорости течений в предположении, что скорость звука в канале постоянна (все
окончательные формулы в [2 - 6] содержат только линейные члены относительно 0cϑ
r
, где
с0 – скорость звука у источника. Члены, порядка ( )2
0cϑ
r
отбрасываются, т.к. 0cϑ
r
имеет
порядок 10-3.
На рис. 1а представлен профиль изменения скорости звука с(z) по всей глубине океана
для района расположения приемного судна. Ось ПЗК находилась на глубине около 300 м.
Значение скорости звука у дна превышало его значение у поверхности на 40-50 м/с. Глубина
перемешивания приповерхностного слоя составила 35-45 м. Подробное гидрологическое
обследование трассы дальнего распространения было проведено несколько дней спустя
после окончания опыта. На рис. 1б для верхнего 1000-метрового слоя представлены профили
с(z), зарегистрированные в начале, в конце и на середине обследованной трассы. Всего вдоль
трассы гидрологическим зондом "Исток-3" было зарегистрировано 15 профилей изменения
температуры и солености (электропроводности) с глубиной. Максимальная глубина
зондирования составила 1000-1200 м. Шаг по расстоянию между соседними
гидрологическими станциями составил 25-30 км. По материалам этого обследования для
800-метрового подповерхностного слоя воды было построено поле значений скорости звука,
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
36
представленное на рис. 1в. В приповерхностном перемешанном слое воды скорость звука по
мере продвижения на юг вдоль обследуемой трассы изменялась от 1493 до 1510 м/с.
Наиболее резкие изменения скорости звука (температуры) с расстоянием в
приповерхностном слое наблюдались на участке трассы 130-160 км. Именно с этого участка
трассы и далее на юг происходило проникновение теплых вод на глубины до 300-400 м.
Трасса дальнего распространения в этом опыте пересекала северную границу зоны,
выходящую к поверхности на расстоянии -140-160 км от приемного судна. Средний наклон
линий постоянной скорости звука южнее этой границы составил 0.10-0.15°. По мере
продвижения на юг происходило их заглубление, составившее -100-150 м на каждые 50 км
(см. рис. 1в) [10].
а) б) в)
Рис. 1 – Тихий океан. Район субарктического фронта. Профили скорости звука с глубиной:
а) у приемного судна (по всей глубине океана), б) – в начале, в середине и в конце трассы
(для глубин до 1000 м), в) – поле скорости звука, построенное по материалам гидрологического
обследования трассы дальнего распространения (числа у изолиний – скорость звука в км/с) [10]
Уравнения акустики неоднородной движущейся среды. Переход к безразмерным
переменным
Общие уравнения движения сжимаемой жидкости выражают законы сохранения
вещества, импульса и энергии. Для того чтобы сформулировать эти законы, выберем
некоторую систему координат х, у, z, неподвижную относительно невозмущенной среды.
Пусть, далее, t – время, ( )1 2 3, ,ϑ ϑ ϑ ϑ=
r
- скорость жидкости в этой системе, ρ и µ -
соответственно плотность и вязкость жидкости, а р – скалярное давление в ней. В этих
обозначениях закон сохранения вещества, выражаемый уравнением непрерывности, и закон
сохранения импульса запишутся в форме [9]:
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
37
( ) 0div
t
ρ ρϑ∂ + =
∂
r
(1)
1
3
p div g
t
ϑρ µ ϑ µ ϑ ρ∂ = −∇ + ∆ + ∇ +
∂
r
r r r , (2)
где g – вектор ускорения силы тяжести, направленный всегда к центру Земли,
ρ g
r
– сила тяжести, действующая на единицу объёма жидкости,
d dtϑ
r
- полная производная скорости жидкости ϑ
r
по времени:
( ) 2
, ,
2
d rot
dt t t
ϑ ϑ ϑ ϑϑ ϑ ϑ ϑ
∂ ∂= + ∇ = +∇ +
∂ ∂
r r r
r r r r
.
Уравнение сохранения энергии, в свою очередь, можно записать в виде [9]:
dS QT
dt
λ
ρ ρ= ∆Τ +
%
, (3)
где Т - температура жидкости,
S - ее энтропия,
Q% - диссипативная функция
3
; 1 ik iki k
Q S ϑ
=
= ∑% , (4)
которая выражается через тензор деформаций ϑik :
1 2; 2 ; 2 ,
2 3
i k
ii iiik ik ik
ik
S S div
x x
ϑϑ
ϑ µϑ µϑ µ ϑ
∂∂
= + = = −
∂ ∂
r
(5)
λ=ρcvχ,
сv - теплоемкость жидкости при постоянном объёме,
χ - коэффициент теплопроводности жидкости;
( ),
dS S
S
dt t
ϑ∂= + ∇
∂
r
.
К выписанным трем законам гидродинамики (1) – (3) следует еще присоединить
уравнение состояния жидкости, связывающее давление р, плотность ρ и энтропию S:
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
38
P=F(ρ, S).
Пусть теперь в среде, состояние которой описывается величинами ϑ
r
, р, ρ, S,
распространяется звук. Исходное состояние среды будем считать устойчивым, а звук будем
рассматривать как малое колебание. Тогда при прохождении звука все величины ϑ
r
, р, ρ, S
получат малые приращения , , ,ξ π δ σ
r
, то есть ξ
r
- скорость звуковых колебаний, π - давление
звука, δ - изменение плотности среды, σ - изменение энтропии.
Чтобы получить уравнение акустики неоднородной движущейся среды, заменим в
(1) – (3) ϑ
r
на ϑ ξ+
r r
, р на р+π, ρ на ρ+δ, S на S+σ и, ограничиваясь линейным
приближением, отбросим члены второго и более высоких порядков относительно малых
величин , , , .ξ π δ σ
r
Влиянием вязкости и теплопроводности среды на распространение звука
пренебрежем. Это означает, что в линейных уравнениях для , , ,ξ π δ σ
r
отбрасываются члены
пропорциональные вязкости µ и теплопроводности λ (согласно (4), (5), диссипируемое в
жидкости тепло Q% также принадлежит к числу величин, пропорциональных µ), и мы
получаем [9]:
( ) 2
, , ,rot rot p
t
ξ π δξ ϑ ϑ ξ ϑ ξ
ρ ρ
∂ ∇
+ + + ∇ = − + ∇ ∂
r
r r r r r r
, (6)
( ) ( ), , 0div div div
t
δ ξ ρ ϑ δ ρ ξ δ ξ δ ϑ∂ + ∇ + ∇ + + + =
∂
r r r r r
, (7)
( ) ( ), , 0S
t
σ ϑ ξσ∂ + ∇ + ∇ =
∂
r r
, (8)
2c hπ δ σ= + , (9)
где с2 – адиабатическая скорость звука:
2 ,
S
p p
с h
S ρρ
∂ ∂ = = ∂ ∂
.
Пусть для рассматриваемой жидкости р0 и Т0 – характерные значения давления и
температуры, а ρ0 - значение плотности ρ при р=р0 и Т=Т0. Определим характерную
адиабатическую скорость звука формулой:
0
0
0
р
с γ
ρ
= ,
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
39
где γ=cp/cv – коэффициент Пуассона (cp, cv – теплоемкости среды при постоянном
давлении и объёме соответственно), и перейдём от размерных переменных х, у, z, t к новым
безразмерным переменным х*=х/l, у*=у/l, z*=z/r, t*=c0t/l. Здесь r – характерный размер
изменения свойств среды по вертикали, а l - характерный размер изменения свойств среды в
горизонтальном направлении.
Введем также новые безразмерные величины *, *, *,pϑ ρ
r
описывающие состояние
среды, *, *, *,ξ π δ
r
описывающие процесс прохождения звука через среду
0 0
0 0
0 0 0 0 0 0 0 0
* , max ; * ; * ; * , max ; * , * .
c cp
p
p p
ϑ ρ ξ π δϑ ϑ ϑ ρ ξ ξ ξ π δ
ϑ ρ ξ ξ ρ ξ
= = = = = = = =
r r
r r r r
Очевидно, что р*, ρ* - величины порядка единицы; * 1, * 1.ϑ ξ≤ ≤
r r
Предположим, что при распространении звуковых колебаний 0 0
0 0 0 0
и
с с
р
π δ
ξ ρ ξ
также
будут величинами порядка единицы (нетрудно убедиться, что в случае плоской волны это
условие действительно выполнено).
Если теперь в качестве уравнения состояния взять
0 0
v
S S
cp
p
γ
ρ
ρ
−
=
%
l , (10)
где S% - энтропия среды при р=р0, ρ=ρ0,
то, заменяя р на р+π, ρ на ρ+δ, S на S+σ и отбрасывая члены второго и более высоких
порядков малости относительно малых величин π, δ и σ, получим
0 0 vp c
π δ σγ
ρ
= + .
Значит, σ/сv, c0/ξ0 по порядку величин не превосходит единицы, поэтому определим
безразмерное изменение энтропии σ*, вызванное прохождением звука, равенством
0
0
*
v
c
с
σσ
ξ
= .
Для того, чтобы ввести безразмерную энтропию среды, обратимся к уравнению (3)
сохранения энергии. Его левая часть
dS
Т
dt
представляет собой прирост тепла единицы массы
вещества и, значит, этот прирост обусловлен исключительно теплопроводностью и работой
сил трения.
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
40
Если среда неподвижна, т.е. 0,ϑ =
r r
то диссипативная функция обращается в нуль, и
уравнение (3) сохранения энергии принимает вид:
dS S
T T
dt t
λ
ρ
∂= = ∆Τ
∂
.
В этом случае весь прирост тепла единицы массы вещества связан только с
теплопроводностью.
При 0ϑ =
r r
естественно предположить, что работа сил трения вызывает прирост тепла
единицы массы вещества по порядку величины равный Т (ϑ
r
, ∇S).
Если Т*=Т/Т0 – безразмерная температура, а
0
*
S S
S
S
−=
%
- безразмерное приращение
энтропии, то
( ) 0 0 0 1 * * *
, * *, * , где *
* * *
T S S S S
T S T S S i j q k
r q x y z
ϑϑ ϑ ∗ ∗
∂ ∂ ∂∇ = ∇ ∇ = + + ∂ ∂ ∂
rr r r r
,
Q=ℓ/r – отношение характерного размера изменения свойств среды в горизонтальном
направлении к характерному размеру изменения свойств среды по вертикали; в океане q
имеет порядок 102.
С другой стороны,
2
* *0
2
0
31
* , 1
Q
Sik ikr i k
µϑ ϑ
ρ ρ ρ
=
=
∑
%
.
Здесь
( )
** *
* * * * *
11 22 33
** * *
* * * *
12 21 13 31* * * *
* *
* * * * * *
23 32 * *
1 1
* , , ; ; ; ;
* * *
1 1 1
; ;
2 2
1 1 2 1
; 2 ; 2
2 3
yx z
x y z
yx x z
y z
ik ik ii ii
q x q y z
q y x z q x
S S div
z q y q
ϑϑ ϑϑ ϑ ϑ ϑ ϑ ϑ ϑ
ϑϑ ϑ ϑϑ ϑ ϑ ϑ
ϑ ϑϑ ϑ ϑ ϑ ∗
∂∂ ∂= = = =
∂ ∂ ∂
∂ ∂ ∂ ∂= = + = = + ∂ ∂ ∂ ∂
∂ ∂= = + = = − ∂ ∂
r r
*;ϑ
r
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
41
** *
*
* * *
.yx zdiv q
x y z
ϑϑ ϑϑ∗
∂∂ ∂= + +
∂ ∂ ∂
r
Поэтому введем S0 таким образом, чтобы
2
0 0 0 0
2
0
Т S
r r
ϑ µϑ
ρ
= ,
то есть положим 0
0
0 0
.S
T r
µϑ
ρ
=
После перехода в уравнениях (6) – (9) к безразмерным переменным и безразмерным
величинам *, *, *, *; *, *, *, *,p Sϑ ρ ξ π δ σ
r r
получаем
( )* *, * *, * *, ***
2
*
1 1 ** *.
* ***
t
rot rot
p
p
ξ ξ ϑ ϑ ξ ϑ ξε ε ε
δπ
ρ γ γ
∂ + ∇ ∂
∇ ∇
+ + =
= − +
r
r r r r r r
(11)
*
*
* **, * * 0*, *
*
div
t
divδ ρ ξξ ϑδ εδ∂ + ∇ + + ∂
=
rr r
, (12)
( ) ( )0*
*, * *, * 0
* v
S
S
t c
σ ϑ σ ξε ∗
∂ + ∇ + ∇ =
∂
r r
, (13)
π* = c*2δ* + h*ε*, (14)
где γ - коэффициент Пуассона;
ε = ϑ0/с0.
В океане ε≤2,5/1500<2·10-3, т.е. ε - безразмерный малый параметр задачи
* ** ** *
* * ;
* * * * * *
y yx xz zrot q i q j k
y z z x x y
ξ ξξ ξξ ξξ
∂ ∂ ∂ ∂∂ ∂= − + − + − ∂ ∂ ∂ ∂ ∂ ∂
rr r r
*2 * *
; *
* *
S
p p
c h
S ρρ
∂ ∂ = = ∂ ∂
.
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
42
Оценим параметр 0 0
0 0
.
v v
S
c T c r
µϑ
ρ
= Так как для воды µ=1·10-2 г/см·с, ρ0=1 г/см3,
Т0=3·102 град., сv=4·107 см2/с2
·град, то в случае, когда ϑ0=2·102 см/с, r=1·105
см, мы получим,
что 150 1,5 10
v
S
c
−≈ ⋅ и, значит, не умаляя общности, можно положить множитель S0/cv равным
нулю, а уравнение (13) заменить уравнением
( )*
*, * 0
*t
σ ϑ σε∂ + ∇ =
∂
r
. (15)
Если теперь в качестве уравнения состояния взять (10), то для коэффициента h*
получится следующее выражение
0*
*
* v
S S
Sp c
h
S c
ν
ρ
−
∂ = = ∂
%
l ,
так как
0
* .
S S
S
S
−=
%
Таким образом h* также можно положить равным нулю и заменить уравнение (14) на
π*=(с*)2δ*. (16)
При этом
∇*Ρ*=(с*)2∇*ρ*+h*∇*S*≈ (с*)2∇*ρ*. (17)
Система (11), (12), (15), (16) содержит малый параметр ε=ϑ0/с0, поэтому будем искать
*ξ
r
, π*, δ*, σ* в виде:
0 0 0 0
* ; * ; * ; *m m m m
m m m m
m m m m
ξ ξ ε π π ε δ δ ε σ σ ε
∞ ∞ ∞ ∞
= = = =
= = = =∑ ∑ ∑ ∑
r r
. (18)
Тогда нулевое приближение в (18) будет соответствовать решению задачи о
распространении звука в неоднородной неподвижной среде.
Представление решения задачи о распространении звука в неоднородной
движущейся среде в виде (18)
Подставим ряды (18) в систему (11), (12), (15), (16). Тогда в нулевом приближении по
ε получим:
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
43
0
* *
0
0
2
1
* *2
*
*
t
cξ
γρ
δ
π ρ
γρ
∂ = − ∇ ∇
∂
+
r
, (19)
( )0
0 * * 0, * * 0
*
div
t
δ ξ ρ ρ ξ∂ + ∇ + =
∂
r r
, (20)
0 0
t
σ∂ =
∂
, (21)
π0 = (с*)2δ0. (22)
Из равенства (21) следует, что σ0 – не зависит от t*, то есть
0 0 * 0
0
t
σ σ
=
= =
Уравнение (19), в свою очередь, с помощью формул (17) и (22) можно переписать в
форме:
0 0
*
1
* *t
ξ π
γ ρ
∂ = − ∇ ∂
r
.
Положив 0
0 *
π
ρ
Π = , мы получим для определения 0ξ
r
и П0 систему:
0
0
1
,
*t
ξ
γ ∗
∂ = − ∇ Π
∂
r
(23)
( )0
0 * * 02
*
, * * 0
( *) *
div
c t
ρ ξ ρ ρ ε∂Π + ∇ + =
∂
r r
. (24)
Исключим 0ξ
r
. Для этого продифференцируем уравнение (24) по t* и воспользуемся
формулой (25):
( )
2 *2 *2
* 0 * 0 **2
, *
*
c c
t
ρ
γ ρ γ
∂ Π = ∆ Π + ∇ Π ∇
∂
. (25)
Следует заметить, что (25) нетрудно преобразовать в уравнение типа волнового. С
этой целью достаточно вместо П0 ввести новую искомую функцию 0
~Π , положив
0 0*ρΠ = Π% .
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
44
Функция 0
~Π будет решением уравнения
2 *2
0
* 0 0*2
,
c
f
t γ
∂ Π = ∆ Π + Π
∂
%
%
где
( )
*2
* * **2
1
*, * 2 * *
4
c
f ρ ρ ρ ρ
γρ
= ∇ ∇ − ∆ . (26)
Для определения поправки, обусловленной наличием течения, подставим ряды (18) в
систему (11), (12), (15), (16) и приравняем коэффициенты при ε в левых и правых частях
полученных уравнений. Тогда для нахождения первого приближения 1, 1 1 1, ,ξ π δ σ
r
будем
иметь систему:
1
* *
1
1 0
2
1
* * *2
*
*
t
c
Lξ
γρ
δ
π ρ
γρ
∂ + ∇ − ∇
∂
=
r r
, (27)
( )1
1 * * 1 0, * *
*
div M
t
δ ξ ρ ρ ξ∂ + ∇ + =
∂
r r
, (28)
( )1
* 0*, 0
*t
σ ϑ σ∂ = − ∇ =
∂
r
, (29)
π1 = (с*)2δ1. (30)
Здесь введены обозначения:
( )* *
0 * 0 * 0, ,L rot ϑ ξ ϑ ξ = − − ∇
r r r rr
, (31)
( )* *
0 * 0 0 *,M divϑ δ δ ϑ= − ∇ =
r r
. (32)
Из соотношения (29) сразу же следует, что
1 1 * 0
0
t
σ σ
=
= = .
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
45
Сведем оставшиеся уравнения системы к неоднородному уравнению для давления π1.
Подставляя выражение для δ1 из (30) в уравнения (27), (28) и полагая
1
1 *
π
ρ
Π = ,
получаем
1
* 1 0
1
1
L
t
ξ
γ
∂ + ∇ Π =
∂
r
r
, (33)
( )1
1 * * 1 0*2
*
, *
*
div M
с t
ρ ξ ρ ρ ξ∂Π + ∇ + =
∂
r r
. (34)
Продифференцировав затем (16) по t* и воспользовавшись формулой (15), находим,
что П1 – решение уравнения
( )
2 *2 *2
1
* 1 * 1 * 0*2
, *
*
c c
F
t
ρ
γ γρ
∂ Π = ∆ Π + ∇ Π ∇ +
∂
, (35)
где
*2
0
0 * 0 0 ** ( , *)
* *
M c
F div L L
t
ρ ρ
ρ
∂= − − ∇ ∂
r r
.
Если же перейти к новой искомой функции 1 1 1* , тоρΠ = Π Π% % , функция будет
решением уравнения типа волнового (26):
2 *2
1
1 1 0*2
*x
c
f F
t
ρ
γ
∂ Π = ∆ Π + Π +
∂
% % .
Зная П1, с помощью (35) получаем :1ξ
r
* *
1 * 1 0
1
* *
t t
dt L dtξ
γ
= − ∇ Π +∫ ∫
r r
.
Вполне аналогично, для m–го приближения будем иметь (m=2,3…)
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
46
* *
* 1
* 0
1
*;
0.
t t
m m m
m m t
dt L dtξ
γ
σ σ
−
=
= − ∇ Π +
=
∫ ∫
r
Здесь Пm=πm/ρ* удовлетворяет уравнению
( )
2 *2 *2
* * * 1*2
, *
*m m m
c c
F
t
ρ
γ γρ −
∂ = ∆ Π + ∇ Π ∇ +
∂
.
Функции 1 1,m mF L− −
r
выражается через решения задачи на предыдущем шаге:
( )
( )
*
1 * 1 * * 1 * 1
*2
1
1 * 1 1 *
* *
1 * 1 1 *
, *, , ,
* , * ,
* *
, .
m m m m
m
m m m
m m m
L rot rot
M c
F div L L
t
M div
ξ ϑ ϑ ξ ϑ ξ
ρ ρ
ρ
ϑ δ δ ϑ
− − − −
−
− − −
− − −
= − − − ∇
∂= − − ∇ ∂
= − ∇ −
r r r r r rr
r r
r r
Все предыдущие рассмотрения никак не использовали тот факт, что исследуемая
задача имеет ещё один малый параметр 1/q=r/ℓ - отношение характерного размера изменения
свойств среды по вертикали к характерному размеру изменения свойств среды в
горизонтальном направлении. Наличие такого малого параметра позволяет существенно
упростить решение задач по нахождению П0 и последующих приближений Пm, m=1,2…
благодаря возможности использования аппарата пространственно-временного лучевого
метода [7, 8].
Применение пространственно-временного лучевого метода для нахождения
коэффициентов рядов (18)
Будем предполагать, что жидкость ограничена сверху свободной поверхностью z=0, а
снизу – твердым дном z=-Н(х,у). Пусть Р - постоянное давление на свободной поверхности.
Условие непрерывности давления на поверхности z=0 и кинематическое граничное условие
даются уравнениями
0 0
; 0zz z
P P ϑ
= =
= = .
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
47
На твердом дне z=-Н(х,у) должна обращаться в нуль нормальная компонента скорости
ϑ
r
:
0z x yx y
ϑ ϑ ϑ∂Η ∂Η+ + =
∂ ∂
.
Таким образом, коэффициенты рядов (18) (m=1,2…) должны удовлетворять
граничным условиям:
,* 0 * 0
0, 0m m zz z
π ξ
= =
= = , (36)
, , ,
1 * *
0 при
* *m z m x m yq x y
ξ ξ ξ ∂Η ∂Η+ + = ∂ ∂
( ) ( )1
* * , * * * , *z H x y H x y
r
= − = −l l l l . (37)
В частности, если m=0, то, продифференцировав равенство (37) по t* и
воспользовавшись уравнением (23), мы получим, что при * *( * , * )z H x y= − l l
0 0
2
1 * *
0
* * * * *z q x x y y
∂Π ∂Π ∂Η ∂Π ∂Η+ + = ∂ ∂ ∂ ∂ ∂
. (38)
Будем искать решение задачи (25), (36), (38), на основе пространственно-временного
лучевого метода, в виде ряда по степеням малого параметра iq-1:
( ) ( )( )∑
∞
=
=Π
0j
0
j
*t*,y*,xQ
i
q
0 *t*z*,y*,xAe , (39)
где Q(x*, y*, t*) и Aj
(0)(x*, y*, z*, t*); j=0,1,2 функции, подлежащие определению.
Подставляя разложение (39) в (25) и сокращая на экспоненциальный множитель,
получаем уравнение:
( )
( ) ( )2 0 02
02
*2
0 0 0
2
* * *
j j j
j j
j
j j j
A AQ i q Q i i
q A
t q i t t q t q
∞ ∞ ∞
= = =
∂ ∂ ∂ ∂− + + + ∂ ∂ ∂ ∂
∑ ∑ ∑
( ) ( )
22 *2 2
0 02
*2 2
0 0
2
*
j j
j j
j j
q Q i c Q i q Q
A q A
i t q x q i xγ
∞ ∞
= =
∂ ∂ ∂+ = − + × ∂ ∂ ∂
∑ ∑
(40)
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
48
( )
( )
( ) 20 022
0 2
*2 *2
0 0 0 *
j j j
j j
j
j j j
A Ai q Q i i Q
A q
x q i x q x q y
∞ ∞ ∞
= = =
∂ ∂ ∂ ∂× + + − × ∂ ∂ ∂ ∂
∑ ∑ ∑
( )
( )
( )
0 2
0 0
*2
0 0 0
2
* *
j j j
j
j j
j j j
Ai q Q i q Q i
A A
q i y y q i y q
∞ ∞ ∞
= = =
∂ ∂ ∂× + + + ∂ ∂ ∂
∑ ∑ ∑
( ) ( )0 02 2 *2
2
*2 *2
0 0
*
* * *
j j
j j
j j
A Ai i c q Q
q
y q z q x i x
ρ
γρ
∞ ∞
= =
∂ ∂ ∂ ∂ + + + × ∂ ∂ ∂ ∂
∑ ∑
( )
( )
( )
0
0 0
0 0 0
*
* * *
j j j
j
j j
j j j
Ai i q Q i
A A
q x q y i y q
ρ∞ ∞ ∞
= = =
∂ ∂ ∂× + + + ∂ ∂ ∂
∑ ∑ ∑
( ) ( )0 0
2
0 0
*
0
* * *
j j
j j
j j
A Ai i
q
y q z z q
ρ∞ ∞
= =
∂ ∂ ∂ + + = ∂ ∂ ∂
∑ ∑ .
Граничные условия (36), (37), в свою очередь, дают:
( )0
0
* 0
0
j
j
j
z
i
A
q
∞
=
=
=
∑ , (41)
( )
( )
( )
( )
( )
0 0
0
2
0 0 0
0
0
0 0
1
* * * *
*
0 * *.
* * *
j j j
j j
j
j j j
j j
j
j
j j
A Ai q Q i i
A
z q q x i x q x q
Aq Q i i
A при z H
y i y q y q
∞ ∞ ∞
= = =
∞ ∞
= =
∂ ∂ ∂Η ∂+ + + ∂ ∂ ∂ ∂
∂ ∂Η ∂ + = = − ∂ ∂ ∂
∑ ∑ ∑
∑ ∑
(42)
Приравнивая к нулю коэффициент при q2 в уравнении (40) и коэффициенты при q0 в
граничных условиях (41), (42), будем иметь
( )
( )
( )
0
20 02 (0)0
0 0 0 0*2
*
* * 0
* *
A
W A A A
z z c
γρρ ω ρ υ
∂∂ = + − = ∂ ∂
r
, (43)
( )
( )0
0 0
0
* 0
* *
0; 0
*z
z H
A
A
z=
=−
∂= =
∂
. (44)
Здесь введем обозначения:
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
49
,
* * *
Q Q Q
i j Q
t x y
ω υ ⊥
∂ ∂ ∂= − = + = ∇
∂ ∂ ∂
r rr
.
Задача (43), (44) – задача на нахождение собственных значений. Она определяет
некоторую связь между иω υr :
( )υω=ω
r
. (45)
Соотношение (45) называется дисперсионным уравнением. Это – аналог уравнения
Эйконала.
Учитывая, что , ,
*
Q
Q
t
ω υ ⊥
∂= − = ∇
∂
r
соотношение (45) можно записать в форме
уравнения Гамильтона-Якоби
*, *, , 0
* * *
Q Q Q
x y
t x y
ω ∂ ∂ ∂+ = ∂ ∂ ∂
, (46)
где роль функции Гамильтона играет функция ω.
Характеристиками уравнения (46) являются кривые х*=х*(t*), у*=у*(t*), вдоль
которых заданы
* * *, ,
* *t x y
Q Q Q
Q
t x y
υ υ∂ ∂ ∂= = =
∂ ∂ ∂
.
Эти кривые находятся из классической канонической системы, соответствующей
уравнению Гамильтона-Якоби (46):
*
* *
* * *
1, , , ,
*
x
x y
dt dx dy
S dS dS dS x
υω ω ω
υ υ
∂ ∂ ∂ ∂= = = =
∂ ∂ ∂ ∂
* *, 0
*
y t
d dQ
dS y dS
υ ω∂= =
∂
.
(47)
Вектор
* *
,гр
x y
ω ωϑ
υ υ
∂ ∂ = ∂ ∂
r
,
называется групповой скоростью. Первые три уравнения системы (47) показывают, что точка
х*=х*(t*), у*=у*(t*) движется с групповой скоростью.
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
50
Кривые t*=S, х*=х*(S), у*=у*(S) в пространстве х*, у*, t*, определяемые из системы
(47), называются пространственно-временными лучами.
Перейдем к нахождению амплитуды нулевого приближения А0
(0). Для этого
рассмотрим задачу для первого приближения А1
(0). Приравнивая нулю коэффициент при iq в
выражении (40) и коэффициенты при iq-1 в граничных условиях (41), (42), будем иметь (см.
также (43))
( ) ( )0 0
0 1 1 0W A W A = , (48)
где введено обозначение:
( )
( )
( )( ) ( ) ( ) ( ) ( )
0
0 0 0 0 00
1 0 1 0 1 0 0 1 0*2 *2
* *
2 2 * , *, *
* *
A
W A A A A QA
c t c t
γρ γρ ωω ρ υ ρ ν ρ∂ ∂ = + ∇ + ∇ + + ∆ ∂ ∂
r r
, (49)
2 2
*2 *2
,
Q Q
Q
x y⊥
∂ ∂∆ = +
∂ ∂
причем
( )
( )
( ) ( )
0
0 01
1 1 0
* 0
0; *, * *
*z
A
A A при z H
z
ν
=
∂= = ∇ Η = −
∂
r
. (50)
Задача (48), (50) – неоднородная задача Штурма-Лиувилля на собственном числе. Для
того, чтобы найти условие ее разрешимости, умножим обе части уравнения (48) на А0
(0),
проинтегрируем полученное равенство по z* от –Н* до нуля и воспользуемся уравнением
(43):
( )
( )
( ) ( )( )
0 00
0 0 00
0 0 1 0*2
* *
*
2 * 2 * , , *
*H H
A
A dz A A dz
c t
ρωγ ρ ν
− −
∂ + ∇ +
∂∫ ∫
r
( ) ( ) ( )
0 02 2
0 0 *
0 0*2
* *
*
*, * *
*H H
A dz A dz
c t
γρ ωρ ν ρ⊥ ⊥
− −
∂ + ∇ + ∆ + + ∂
∫ ∫
r
( ) ( ) 20
1 0
* - *
* *,
z H
Ñ Aρ ν
=
+ Η
r
.
(51)
Соотношению (51) нетрудно придать форму пространственно-временного закона
сохранения энергии
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
51
( ) ( )
20 0
0 0
0 0*2
* *
*
* * * 0
* H H
A dz div A dz
t c
ργω ν ρ⊥
− −
∂ + = ∂
∫ ∫
r
. (52)
Здесь
* *
bx by
div b
x y⊥
∂ ∂= +
∂ ∂
r
.
Пусть теперь А0
(0)(х*, у*, z*, t*) – решение уравнения (43), зафиксированное
соотношением:
( )
0 20
0*2
*
*
* 1
H
A dz
c
ρ
−
= ∫ % . (53)
Тогда
( ) ( ) ( ) ( )0 0 0
0 0 0*, *, *A a x y t A= % ,
где скаляр а0
(0)(х*, у*, t*) играет роль амплитудного множителя, и
( ){ } ( ) ( )
02 2
0 0 0
0 0 0
*
* * 0
* H
a div A dz a
t
γω ν ρ⊥
−
∂ + = ∂
∫
r % . (54)
Запишем равенство (54) в более удобной форме. Для этого обратимся к уравнению
(43) для ( )0
0А
~ . Умножим обе его части на ( )0
0А
~ и проинтегрируем по z* от –Н* до нуля. Тогда
мы получим, что
( ) ( )
( )
( )
0 0 0 02 220 0 02 0
0 0 0*2
* * *
*
* * * * * 0
* *H H H
A
A dz A dz A dz
c z z
ργω ν ρ ρ
− − −
∂∂ − + = ∂ ∂
∫ ∫ ∫
%r% % % ,
или
( )
( ) 20 0 022 02 0
0
* *
* * * * 0
*H H
A
A dz dz
z
γω ν ρ ρ
− −
∂ − − = ∂
∫ ∫
%r % .
Продифференцируем теперь последнее равенство по νr :
( ) ( )
( )0 0 02 20 0 0
0 0
* *
2 2 * * 2 * *
H H
A
A dz A dz
ωγω ν ρ ν ρ
ν ν− −
∂∂ − − − ∂ ∂∫ ∫
%r r% %
r r (55)
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
52
( ) ( )0 0 02
0 0
*
*
2 * * 0
*H
A A
dz
z z
ρ
ν−
∂ ∂− =
∂ ∂ ∂∫
% %
r .
С другой стороны, умножая обе части уравнения (43) для ( )0
0А
~ на ( )0
0А
~∂ и интегрируя
по z* от –Н* до нуля, будем иметь:
( ) ( )
( )
( )0 00 0 0
020 0 0
0*2
* *
*
* * *
* *
Н H
A A A
dz A dz
z z c
ρρ γω
ν ν− −
∂ ∂ ∂∂ + − ∂ ∂ ∂ ∂
∫ ∫
% % %
%
r r
( )
( )0 0
2 0 0
0
*
* * 0
H
A
A dzν ρ
ν−
∂− =
∂∫
%r %
r ,
или
( )
( ) ( ) ( )0 00 0 02
2 0 0 0 0
0
* *
* * * * 0
* *H H
A A A
A dz dz
z z d
ν ρ ρ
ν ν− −
∂ ∂ ∂− − =
∂ ∂ ∂∫ ∫
% % %r %
r r , (56)
так как
( )
( )
( )
0 00 2
0 00
0 0*2 *2
* *
* 1 *
* * 0
2H H
A
A dz A dz
с с
ρ ρ
ν ν− −
∂ ∂ = = ∂ ∂∫ ∫
%
% %
r r .
Таким образом (см. (55), (56)),
( )
0 20
0
*
* *
Н
A dz
ωγω ν ρ
ν −
∂ = ∂ ∫
r %
r ,
и значит
( )
0
20
0
*
1
, * *гр
x y Н
A dz
ω ω ν ωϑ ρ ν
ν ν ν ν γω −
∂ ∂ ∂ = = = ∂ ∂ ∂
∫
r
r r%
r r .
Используя полученное выражение для групповой скорости, можно соотношение (54)
записать в форме
( ){ } ( ){ }2 20 0
0 1 0 0
*
a div a гр
t
γω γω ϑ∂ + = ∂
r
. (57)
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
53
Проинтегрируем уравнение (57). Будем считать, что у нас имеется поле
пространственно-временных лучей, не имеющее особенностей в интересующей нас области.
Рассмотрим лучевую трубку – совокупность лучей, пересекающих плоскость t*=0 в точках
бесконечно малого прямоугольника х0
*
≤х
*
≤х0
*+dх0
*, у0
*
≤х
*
≤у0
*+dу0
*
Проинтегрируем соотношение (57) по объёму Ω лучевой трубки, вырезанному
плоскостями t*=0 и t*=const>0. Применяя формулу Остроградского, получим
( ) ( ) ( ){ }20
0 .cos , , 0грa n t n dγω ϑ∧ + Σ = ∫
rrr r
. (58)
Здесь ∂Ώ – граница объёма Ώ.
Интеграл по боковой поверхности лучевой трубки равен нулю, так как вектор
{ }.., ; , ,1гр грх уϕ ϑ ϑ=r
направлен вдоль луча, т.е. лежит в касательной плоскости к боковой
поверхности трубки, а подынтегральное выражение пропорционально ( ), nϕr r , где n
r - внешняя
нормаль. Таким образом, из соотношения (58) следует равенство интегралов по торцам
выделенной трубки и значит
( ) ( ) ( ) ( ) ( ) ( )2 20 0* * * *
0 * 0 0 0 0 0 0*, *, * *, *, * , ,0 , ,0tx y t a x y t d x y a x y dω ω Σ = Σ , (59)
или
( ) ( ) ( ) ( ) ( ) ( )2 20 0* * * * 0
0 0 0 0 0 0
*
*, *, * *, *, * , ,0 , ,0
t
d
x y t a x y t x y a x y
d
ω ω Σ = Σ
.
Пространственно-временной луч в рассматриваемом поле лучей определяется той
точкой х0
*, у0
*, в которой он пересекается с плоскостью t*=0, поэтому х*=х*(х0
*, у0
*, t*),
у
*=у*(х0
*, у0
*, t*)
Таким образом, равенство (59) можно переписать в виде:
( ) ( ) ( ) ( ) ( ) ( )2 20 0* * * *
0 0 0 0 0 0
1
*, *, * *, *, * , ,0 , ,0x y t a x y t x y a x y
J
ω ω = ,
* *
0 0
( *, *)
( , )
D x y
J
D x y
= .
Решение задачи (49), (50), в свою очередь, будем искать в форме:
( )0(0) (0) (0)
1 1 1 0А А а А= +% % .
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
54
Здесь ( )0
1А
~ - частное решение (49), удовлетворяющее условиям (50), а функция
а0
*(х*, у*, t*) определяется из условия разрешимости задачи для А2
(0) и т.д. Зная П0, находим
0ξ
r
:
( ) ( ) ( )
*
* * *, , 0
0 * 0
0
1 q
i
t
Q x y t ji
j q
j
dt eξ α
γ
∞
=
= − ∇ Π = ∑∫
r r
,
в частности,
( ) ( ) ( ) ( )
( )0
0 0 0 0 0
0, 0, 0 0, *
1
;x y z
Ai
i j A
z
α α ν α
γω γω
∂+ = =
∂
r r r
. (60)
Обратимся теперь к первому приближению по ε , то есть к поправке, обусловленной
наличием течения. Будем искать П1=π1/ρ* в виде:
( ) ( ) ( )
* * *, , 1 * * * *
1
0
, , ,
jq
Q x y t
i
j
j
i
q e A x y z t
q
α
∞
=
Π =
∑ , (61)
где Q(x*,y*,t*) – функция, входившая в формулу (39), а постоянная α подлежит
определению.
Подставляя разложение (61) в (52), получим, что при q→∞ (40)
( ) ( ) ( )* * *2 *2 *2 , ,* 21
* 1 * 1 **2 *
q
Q x y t
i
П c c
O q e
t
αρ
γ γρ
+∂ − ∆ Π − ∇ Π ∇ =
∂
.
Обратимся теперь к
( )
*2
*0
* * 0 0 ** *
,
M c
F div L L
t
ρ ρ
ρ
∂ = − − ∇ ∂
r r
,
где 0 0,L M
r
определяется равенствами (48), (49). Получим вначале выражение для 0 :L
r
Так как
( ) ( ) ( )( )
* * *, , 0* *
0
0
, ,
jq
Q x y t
i
j
j
i
e
q
ϑ ξ ϑ α
∞
=
=
∑
r r r r
,
то
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
55
( ) ( ) ( )
* * *,* * 0
* 0 *
0
,
jq
Q x y t
i
j
j
i
e
q
ϑ ξ ϑ α
∞
=
∇ = ∇ =
∑
r r r r
( )( )( ) ( )( )( )0 0*
1 *
0 0
, ,
q
i
q
Q j jQi ii
j jq q
j j
e qe k
z
ϑ α ϑ α
∞ ∞
= =
∂= ∇ + =∑ ∑ ∂
rr rr r
( )( )( ) ( )( )( )0 0* *
0 0
, ,
q q
i i
j jQ Qi i
j jq q
j j
q
Q e e
i
ϑ α ϑ α
∞ ∞
⊥ ⊥
= =
= ∇ + ∇ +∑ ∑
r rr r
( )( )0*
*
0
,
q
i
j
Q
j
j
i
qe k
z q
ϑ α
∞
=
∂+ ∂
∑
rr r
.
Далее
** **
* 0 0, 0,* * * *
, yx xz
z yrot q i
z x x y
ϑϑ ϑϑϑ ξ ξ ξ
∂ ∂ ∂∂
= − − − + ∂ ∂ ∂ ∂
r r r
* * ** * *
0, 0, 0,* * * * * *
y y yx z z
x z yq j q
x y y z y z
ϑ ϑ ϑϑ ϑ ϑξ ξ ξ
∂ ∂ ∂∂ ∂ ∂+ − − − + − − ∂ ∂ ∂ ∂ ∂ ∂
r
( ) ( )
** **
0
0, ,* * * *
0
q
i
jQ yx xz i
x j z q
j
q k e q
z x x y
ϑϑ ϑϑ ξ α
∞
=
∂ ∂ ∂∂ − − = − − ∂ ∂ ∂ ∂
∑
r
( ) ( ) ( ) ( )
* ** *
0 0
, ,* * * *
0 0
j jy yx xi i
j y j xq q
j j
i
x y x y
ϑ ϑϑ ϑα α
∞ ∞
= =
∂ ∂∂ ∂− − + − − ∂ ∂ ∂ ∂
∑ ∑
r
( ) ( ) ( ) ( )
* ** *
0 0
, ,* * * *
0 0
j jy yz zi i
j z j yq q
j j
q j q
y z y z
ϑ ϑϑ ϑα α
∞ ∞
= =
∂ ∂∂ ∂− − + − − ∂ ∂ ∂ ∂
∑ ∑
r
( ) ( )
* *
0
,* *
0
j
x z i
j x q
j
q k
z x
ϑ ϑ α
∞
=
∂ ∂ − − ∂ ∂
∑
r
.
Поэтому
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
56
( )( )( )20, 00, 0, 2 *
* 0 * * *
0
,
q
i
jQyx z i
j q
j
LL L
div L q e q Q
x y z
ϑ α
∞
⊥
=
∂∂ ∂
= + + = − − ∇ +∂ ∂ ∂
∑
rr r
( )( )( ) ( )( )( )0 0* *
* * * *
, ,
j j
i i
j jq q
j o j o
q Q q
i x x i y y
ϑ α ϑ α
∞ ∞
= =
∂ ∂ ∂ ∂Θ+ + + ∂ ∂ ∂ ∂
∑ ∑
r rr r
( )( ) ( ) ( )( ) ( )0 0* *
* * * *
, ,
j j
i i
j jq q
j o j o
q q Q
i x x i y y
ϑ α ϑ α
∞ ∞
= =
∂Θ ∂ ∂ ∂+ + +
∂ ∂ ∂ ∂∑ ∑
r rr r
( )( )( ) ( )( )( ) ( )( )( )
2 2 2
0 0 0* * 2 *
*2 *2 *2
, , ,
j j j
i i i
j j jq q q
j o j o j o
q
x y z
ϑ α ϑ α ϑ α
∞ ∞ ∞
= = =
∂ ∂ ∂+ + + +
∂ ∂ ∂∑ ∑ ∑
r r rr r r
( ) ( ) ( ) ( )
** **
0 0
, ,* * * * *
0 0
j jyx xz i i
j z j yq q
j j
q
q
i x z x x y
ϑϑ ϑϑ α α
∞ ∞
= =
∂ ∂ ∂∂∂Θ+ − − − + ∂ ∂ ∂ ∂ ∂
∑ ∑
( ) ( ) ( ) ( )
** **
0 0
, ,* * * * *
0 0
j jyx xz i i
j z j yq q
j j
q
x z x x y
ϑϑ ϑϑ α α
∞ ∞
= =
∂ ∂ ∂∂∂+ − − − + ∂ ∂ ∂ ∂ ∂
∑ ∑
( ) ( ) ( ) ( )
* ** *
0 0
, ,* * * * *
0 0
j jy yx zi i
j x j zq q
j j
q
q
i y x y y z
ϑ ϑϑ ϑα α
∞ ∞
= =
∂ ∂∂ ∂∂Θ+ − − − + ∂ ∂ ∂ ∂ ∂
∑ ∑
( ) ( ) ( ) ( )
* ** *
0 0
, ,* * * * *
0 0
j jy yx zi i
j x j zq q
j j
q
y x y y z
ϑ ϑϑ ϑα α
∞ ∞
= =
∂ ∂∂ ∂∂+ − − − + ∂ ∂ ∂ ∂ ∂
∑ ∑
( ) ( ) ( ) ( )
* ** *
0 0
, ,* * * * *
0 0
j jy xz zi i
j y j xq q
j j
q q q
z y z z x
ϑ ϑϑ ϑα α
∞ ∞
= =
∂ ∂∂ ∂∂ + − − − ∂ ∂ ∂ ∂ ∂
∑ ∑ ;
( )
* * * *
*
0 * 0, 0, 0,* * * *
,
q
i Q
x y zL L L qL e
x y z x
ρ ρ ρ ρρ ∂ ∂ ∂ ∂∇ = + + = − ×
∂ ∂ ∂ ∂
r
( )( )( ) ( )( )( )
* *
0 0* *
* * * *
0 0
, ,
j j
x zi i
j jq q
j j
q Q
q
i х x z x
ϑ ϑϑ α ϑ α
∞ ∞
= =
∂ ∂∂ ∂× + + − × ∂ ∂ ∂ ∂
∑ ∑
r r r
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
57
( ) ( ) ( ) ( )
* * *
0 0
, ,* * * *
0 0
q
i
j j Qy xi i
j z j yq q
j j
q Q
e
x y y i y
ϑ ϑ ρα α
∞ ∞
= =
∂ ∂ ∂ ∂× − − − × ∂ ∂ ∂ ∂
∑ ∑
( )( ) ( ) ( )( )( ) ( ) ( )
* *
0 0 0* *
,* * *
0 0 0
, ,
j j jy xi i i
j j j xq q q
j j jy x y
ϑ ϑϑ α ϑ α α
∞ ∞ ∞
= = =
∂ ∂∂× + + − − ∂ ∂ ∂
∑ ∑ ∑
r r
( ) ( ) ( )( )( )
** *
0 0*
,* * * *
0 0
q
i
j jQyz i i
j z jq q
j j
q q e q
y z z z
ϑϑ ρα ϑ α
∞ ∞
= =
∂ ∂ ∂ ∂− − − + ∂ ∂ ∂ ∂
∑ ∑
r r
( ) ( ) ( ) ( )
* ** *
0 0
, ,* * * *
0 0
j jy xz zi i
j y j xq q
j jz
q q
y z x
ϑ ϑϑ ϑα α
ϑ
∞ ∞
= =
∂ ∂∂ ∂ + − − − ∂ ∂ ∂ ∂
∑ ∑ .
Рассмотрим теперь ( )* *
0 * 0 0 *,M divϑ δ δ ϑ= − ∇ −
r r
. Так как
( ) ( ) ( )
* *, ** *, 0
0 0 0*2 *2 *2
0
1 q
Q x y t j
ii
j q
j
e A
c c c
ρ ρδ π
∞
=
= = Π = ∑ ,
то
** *
* * *0 0 0
0 0* * * * * *
yx z
x y zM q q
x y z x y z
ϑδ δ δ ϑ ϑϑ ϑ ϑ δ
∂∂ ∂ ∂ ∂ ∂= − + + + + + = ∂ ∂ ∂ ∂ ∂ ∂
( ) ( ) ( ) ( )
* *
0 0*
* *2 * *2
0 0
q
i
j jQ i i
x j jq q
j j
q Q
e A A
i x c x c
ρ ρϑ
∞ ∞
= =
∂ ∂= + + ∂ ∂
∑ ∑
( ) ( ) ( ) ( )
* *
0 0* *
* *2 * *2
0 0
j j
i i
y j j zq q
j j
q Q
A A q
i y c y c
ρ ρϑ ϑ
∞ ∞
= =
∂ ∂+ + + × ∂ ∂
∑ ∑
( ) ( ) ( ) ( )
** ** *
0 0
* *2 * * * *2
0 0
j jyx zi i
j jq q
j j
A q A
z c x y z c
ϑϑ ϑρ ρ∞ ∞
= =
∂ ∂ ∂∂ × + + + ∂ ∂ ∂ ∂
∑ ∑ ,
и значит
( ) ( ) ( ) ( )
* *
0 0*0
* * * *2 * *2
0 0
q
i
j jQ i i
x j jq q
j j
M q Q q Q
e A A
t i t i x c x c
ρ ρϑ
∞ ∞
= =
∂ ∂ ∂ ∂= − + + ∂ ∂ ∂ ∂
∑ ∑
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
58
( ) ( ) ( ) ( )
* *
0 0* *
* *2 * *2
0 0
j j
i i
y j j zq q
j j
q Q
A A q
i y c y c
ρ ρϑ ϑ
∞ ∞
= =
∂ ∂+ + + × ∂ ∂
∑ ∑
( ) ( ) ( ) ( )
** ** *
0 0
* *2 * * * *2
0 0
j jyx zi i
j jq q
j j
A q A
z c x y z c
ϑϑ ϑρ ρ∞ ∞
= =
∂ ∂ ∂∂ × + + + − ∂ ∂ ∂ ∂
∑ ∑
( ) ( ) ( ) ( )
* *
0 0*
* * *2 * *2
0 0
q
i
j jQ i i
x j jq q
j j
q Q
e A A
t i x c x c
ρ ρϑ
∞ ∞
= =
∂ ∂ ∂− + + ∂ ∂ ∂
∑ ∑
( ) ( ) ( ) ( )
* *
0 0* *
* *2 * *2
0 0
j j
i i
y j j zq q
j j
q Q
A A q
i y c y c
ρ ρϑ ϑ
∞ ∞
= =
∂ ∂+ + + × ∂ ∂
∑ ∑
( ) ( ) ( ) ( )
** ** *
0 0
* *2 * * * *2
0 0
j jyx zi i
j jq q
j j
A q A
z c x y z c
ϑϑ ϑρ ρ∞ ∞
= =
∂ ∂ ∂∂ × + + + ∂ ∂ ∂ ∂
∑ ∑ .
Таким образом,
( ) ( )
* **2 *
0 02 * *
0 0 0* *2 * * * * * *2
q
i Q z
x y
c Q Q Q Q
F e q A i A
c t x y t z c
ϑ ρρ ϑ ϑ
ρ
∂ ∂ ∂ ∂ ∂= + + + ∂ ∂ ∂ ∂ ∂
( )( ) ( )( ) ( ) ( )2 0 0 0 0* * *
0 0 0, 0,*2 * * * *
, , yx
z z
Q Q
Q i i
z x z y z
ϑϑρ ϑ α ϑ α α α⊥
∂∂∂ ∂ ∂+ − ∇ + − − − ∂ ∂ ∂ ∂ ∂
r rr r
( ) ( )( ) ( ) ( )
* ** * *
0 0 0 0*
0, 0 0, 0,* * * * * * *
,y y
x y xz z z z z z z
ϑ ϑρ ρ ρα ϑ α α α
∂ ∂∂ ∂ ∂ ∂ ∂ − + − − + ∂ ∂ ∂ ∂ ∂ ∂ ∂
r r
( ) ( ) ( )
* **2 *
0 02 *
0 0* *2 * *2
( ) ,
q q
i iQ Q zc
e Q q q e A i A
c t c
ϑ ρρ ω ϑ ν ω
ρ
∂+ = − − − − ∂
r r
( )( ) ( )( ) ( )
*
0 0 0* * * * *
0 0 0,* * * *
, , x
xz z z z
ϑν ρ ϑ α ρ ϑ α ρ α ∂∂ ∂ ∂ − + − + ∂ ∂ ∂ ∂
r rr r r
( ) ( )
* *
0 0* *
0, 0,* *
; ( )
q
i Qy
y zi e Q q
z z
ϑ ϑρ α α ρ ν
∂ ∂ + − + ∂ ∂
r
r
,
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
59
где введено обозначение:
( )* * *
* *
; x y
Q Q
x y
ϑ ν ϑ ϑ∂ ∂≡ +
∂ ∂
r r
.
Воспользуемся формулами (60). Тогда выражение для F0 можно будет записать в
форме:
( )
*2
02
0 0 0*
( )
q q
i i
Q Qc
F q e f A e Q q
γρ
= − + .
Здесь
( ) ( ) ( ) ( )
* **
0 0 0*
0 0 0 0*2 * *2
, zf A A i A
c z c
ϑ ρργω ϑ ν
∂ ≡ + + ∂
r r
( )
( ) ( )
( )0 0*
2 20* * * * *0 0
0* * *
1
, , z
A A
A i
z z z
ϑρ ν ρ ν ϑ ν ν ϑ ρ
ω
∂ ∂∂+ − + + − ∂ ∂ ∂
r
rr r r r
( )
( ) ( )0 0
* * * *0 0
* * * * *
, z
A A
i
z z z z z
ρ ϑ ν ρ ϑ
∂ ∂∂ ∂ ∂ − − ∂ ∂ ∂ ∂ ∂
r r
.
(62)
Обратимся теперь к граничным условиям (36), (37) при m=1. Продифференцировав
равенство (37) по t* и воспользовавшись уравнением (33), мы получим, что при z*=-Н*
* * * *
1 1 1
0, 0, 0,* 2 * * * * 2 * *
1
z x yL L L
z q x x y y q q x x
γ γ ∂Π ∂Π ∂Π∂Η ∂Η ∂Η ∂Η+ + = + + ∂ ∂ ∂ ∂ ∂ ∂ ∂
.
Следовательно, α≥0. Если бы α=0, то для нахождения А0
(1) мы получили бы
неоднородную задачу Штурма-Лиувилля на собственном числе (43), (62):
( ) ( )1 0
0 0 0 0W A f A = , (63)
( )
( )
( ) ( ) ( )
* * * ** *
1 0 0 02
* * *0 0 0 0
* * * * *2
1
, z z
z zz
A A A Ai
i
z z z z z
ϑ ν ϑ ϑ
ω ω
=−Η =−Η=−Η
∂ ∂ ∂ ∂∂= + = − ∂ ∂ ∂ ∂ ∂
r r
;
(64)
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
60
( )
*
1
0
0
0
z
A
=
= .
Коэффициент А0
(0) уже определен при рассмотрении приближения, соответствующего
отсутствию течения (ε=0), поэтому f0[А0
(0)] и
( )
* *
02
* 0
*2z
z
Ai
z
ϑ
ω
=−Η
∂−
∂
- известные функции.
Значит, неоднородная задача Штурма-Лиувилля (63), (64), вообще говоря, не имеет
решения. В связи с этим, будем искать П1 в виде
( ) ( ) ( ) ( )* * *, , 1 * * * *
1
0
, , ,
q
i Q x y t ji
j q
j
qe A x y z t
∞
=
Π = ∑ . (65)
Тогда для W0[A0
(1)] получится однородное уравнение W0[A0
(1)]=0 с однородными
краевыми условиями:
( )
( )
*
*
1
1 0
0 *0
0; 0
z
z
A
A
z=
=−Η
∂= =
∂
.
Таким образом, ( ) ( ) ( )1 1 0
0 0 0А а А= % , где ( )0
0А
% - решение задачи (43), (44), зафиксированное
соотношением (53).
Для определения а0
(1) надо рассмотреть задачу для А1
(1). Приравнивая нулю
коэффициент при iq2 в выражении, полученном в результате подстановки разложения (65) в
уравнение (35), будем иметь (см. (43), (49), (62)):
( ) ( ) ( )1 1 0
0 1 0 0 0 0
1
W A W A f A
i
= + , (66)
причем
( )
( ) ( )
( )
* *
01
1* * 01
1 0* * *
1
, z
z
AA
A
z z z
ν ϑ
ω
=−Η
∂∂ ∂= ∇ Η − + ∂ ∂ ∂
r
( )
( )
( ) ( )
( )
* ** *
0 02*
1* *0 0
1 0* *2
, , z
z
zz
A Ai
A
z z
ϑϑ ν ν
ω ω
=−Η=−Η
∂ ∂= ∇ Η − ∂ ∂
r r
,
( )
*
1
1
0z
A
=
=0.
(67)
Задача (66), (67) – неоднородная задача Штурма-Лиувилля на собственном числе. Для
того, чтобы найти условие ее разрешимости, умножим обе части уравнения (66) на А0
(0),
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
61
проинтегрируем полученное равенство по z* от -Н* до нуля и воспользуемся уравнением (43)
и граничными условиями (44), (67):
( ) ( ) ( ) ( ) ( ) ( ) ( )
( )
* * * *
0 0 02*
0 1 0 0 0 1* * * * 0
0 1 0 0 0 0 1 0 *2
, z
Н Н z
A
A W A dz i A f A dz A A
z
ϑρ ν
ω− − =−Η
∂ − = − ∇ Η − ∂
∫ ∫
r
. (68)
Примем теперь во внимание, что
( )
( )
( )
( )
1
1 00
0 00
0
a
A A
a
= .
Тогда
( ) ( ) ( ) ( ) ( )
* * *
0
0 1 0 1* * *
0 0 0 0,
Н z
A W A dz A Aρ ν⊥
− =−Η
+ ∇ Η = ∫
r
( )
( )
( )( ) ( ) ( )
*
0 1*
0 1 1* *0
0 0 0*2 *
2 2 , ,
A
A A A
c t
γρ ω ρ ν ρ ν⊥ ⊥
−Η
∂= − + ∇ + ∇ + ∂
∫
r r
( ) ( ) ( ) ( ) ( )
* *
*
1 1 0 1* * * *
0 1 0 0 0*2 *
,
z
A QA dz A A
c t
γρ ω ρ ρ ν⊥
=−Η
∂+ + ∆ + ∇ Η =∂
r
( )
( )
( )
( )
( )( ) ( )
*
01 0* *
0 0 0*0 0
0 0 0*2 * *2 *0
0
2 2 ,
a A
A A A
c t c ta
γρ ω γρ ωρ ν ⊥
−Η
∂ ∂= + ∇ + + ∂ ∂
∫
r
( )}
( )
( ) ( ) ( )
* *
1 20 0* * * *0
0 0 *0
0
, 2
z
a
QA dz A
ta
ρ ρ ν γω⊥ ⊥
=−Η
∂ + ∆ + ∇ Η + × ∂
r
( )
( )
( )
( )
( )
( )
* *
0 01 1* 2
0 0* * *0 0
0 0*20 0
0 0
2 ,
H H
a a
A dz A dz
ca a
ρ ν ρ⊥
− −
× + ∇
∫ ∫
r
,
и равенства (51), (53) дают
( )
( )
( )
( )
( )
*
01 120* *0 0
0* 0 0
0 0
1
, ,
H
a a
A dz
t a a
ρ ν ϕ
γω ⊥
−
∂ + ∆ = ∂
∫
r% (69)
где введено обозначение:
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
62
( ) ( ) ( ) ( )
( )
* * *
0 02* *
0 0 0* * * * 0
0 0 0 *2
1
, ,
2
z
H z H
A
x y t i A f A dz A
z
ϑ ρϕ
γω ω− =−
∂ = + ∂
∫
%
% % .
Как было показано,
( )
*
0
20* *
. 0
1
гр
H
A dzϑ ρ ν
γω −
= ∫
r r% ;
поэтому уравнение (69) можно переписать в форме:
( )
( )
( )
( ) ( )
1 1
* * *0 0
.* 0 0
0 0
, , ,гр
a a
x y t
t a a
ϑ ϕ⊥
∂ + ∇ = ∂
r
. (70)
На пространственно-временном луче t*=Q, x*=x*(t*), y*=y*(t*) справедливы равенства
* *
. ., , ,гр гр
dx dу
х у
dS dS
ϑ ϑ= = ,
и, значит, уравнение (70) вдоль луча принимает вид обыкновенного дифференциального
уравнения
( )
( ) ( )
1
* * *0
* 0
0
, ,
ad
x y t
dt a
ϕ
=
.
Решение задачи (66), (67), в свою очередь, будем искать в форме:
( ) ( ) ( ) ( )1 1 1 0
1 1 1 0A A a A= +% % ,
где ( )1
1A% - частное решение уравнения (66), удовлетворяющее условиям (67), а
функция а1
(1) определяется из условия разрешимости задачи для А2
(1) и т.д.
Зная П1, с помощью формулы
* *
*
* *
1 * 1 0
1 t t
dt L dtξ
γ ∇
= − ∇ Π +∫ ∫
r r
,
находим 1ξ
r
.
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
63
Следующие приближения по ε (m=2, 3…) рассматриваются аналогично. При этом
разложение
*
m
m
π
ρ
Π = в ряд по степеням i/q будем иметь вид:
( ) ( ) ( )( )* * *, , * * * *
0
, , , ;
q
i
jQ x y t mm i
m j q
j
q e A x y z t
∞
=
Π = ∑
Q(x*, y*, t*) – функция, входившая в формулу (39) и, значит, для π* получается
выражение:
( ) ( )* * *( , , )* * *
0 0 0
q
i
jmQ x y tm m i
m j q
m m m
e q Aπ ρ ε ρ ε
∞ ∞ ∞
= = =
= Π =∑ ∑ ∑ . (71)
Разложение (71) удобно тем, что его нулевое приближение по ε соответствует хорошо
изученной задаче о распространении звука в неподвижной среде. С другой стороны, новый
безразмерный параметр 0
0
0,2q q
c
ϑε = ≤ , поэтому при вычислении по приведенной схеме
звуковых полей в океане надо удерживать довольно много членов полученных рядов.
Выясним, не будет ли для расчетов предпочтительнее второй возможный подход к
решению рассматриваемой задачи, когда * * * *, , ,ξ π δ σ
r
пишутся в виде:
( ) ( ) ( )* * *
0 0 0
, ,
q q q
i i i
m m mG G Gi i i
m m mq q q
m m m
e B e D e Eξ π σ
∞ ∞ ∞
= = =
= = =∑ ∑ ∑
r r
, (72)
( )
( )
*
*
, , ,2*
, , ,m m x m y m zB B B B
c
πδ = =
r
,
здесь , , , 0,1...m m mB D E m =
r
- функции от х*, у*, z*, t*, подлежащие определению,
G=G(x*, y*, t*) также надо найти.
Последовательность рекурентных задач для определения коэффициентов рядов (72)
Подставим выражение (72) для * * * *, , ,ξ π δ σ
r
в уравнение (11), (12), (15). При этом в
уравнении (11):
( ) ( )
*
* *
0 0
;
q q
i i
m mG G mi i
m q q
m m
Bq G
e B e
t i t t
ξ ∞ ∞
= =
∂∂ ∂= +
∂ ∂ ∂∑ ∑
r r
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
64
{ ( ) ( ),* * *
* ,* *
0 0
,
q
i
m mG m x i i
z m zq q
m m
B q G
rot e q B
z i x
ξ ϑ ϑ
∞ ∞
= =
∂ ∂
= − − ∂ ∂
∑ ∑
r r
( ) ( ) ( ),, *
,* * *
0 0 0
q
i
m m mG m ym z i i i
y m yq q q
m m m
BB q G
e B
x i x x
ϑ
∞ ∞ ∞
= = =
∂∂ ∂− − + − ∂ ∂ ∂
∑ ∑ ∑
( ) ( ) {, *
,* * *
0 0
q
i
m m Gm xi i
m x q q
m m
Bq G q G
B xi x e
i y y i x
∞ ∞
= =
∂ ∂ ∂− − + × ∂ ∂ ∂
∑ ∑
r
( ) ( ) ( ) ( ), ,
, ,* * *
0 0 0 0
m m m mm y m xi i i i
m y m xq q q q
m m m m
B Bq G
B B
x i y y
∞ ∞ ∞ ∞
= = = =
∂ ∂ ∂× − − − −∂ ∂ ∂
∑ ∑ ∑ ∑
( ) ( ) ( ),,*
,* * *
0 0 0
q
i
m m mG m ym zi i i
z m z q q q
m m m
BBq G
e B q j
i y y z
ϑ
∞ ∞ ∞
= = =
∂ ∂ ∂ − + − + ∂ ∂ ∂
∑ ∑ ∑
r
( ) ( ) ( ),*
, ,* *
0 0 0
q
i
m m mG m yi i i
y m z m zq q q
m m m
Bq G
e B B q
i y z
ϑ
∞ ∞ ∞
= = =
∂ ∂+ + − − ∂ ∂
∑ ∑ ∑
[ ( ) ( ) ( ), ,*
,* * *
0 0 0
;
q
i
m m mG m x m zi i i
x m zq q q
m m m
B Bq G
e q B k
z i x x
ϑ
∞ ∞ ∞
= = =
∂ ∂ ∂− − − ∂ ∂ ∂
∑ ∑ ∑
r
( )
** **
*
* ,* * * *
0
,
q
i
mG yx xz i
m z q
m
rot e q B
z x x y
ϑϑ ϑϑϑ ξ
∞
=
∂ ∂ ∂∂ = − − − × ∂ ∂ ∂ ∂
∑
r r
( ) ( )
* ** *
, ,* * * *
0 0
q
i
m mG y yx zi i
m y m xq q
m m
B i e B q
x y y z
ϑ ϑϑ ϑ∞ ∞
= =
∂ ∂∂ ∂ × + − − − × ∂ ∂ ∂ ∂
∑ ∑
r
( ) ( )
* ** *
, ,* * * *
0 0
q
i
m mG y xz zi i
m z m yq q
m m
B j e q B q
y z z x
ϑ ϑϑ ϑ∞ ∞
= =
∂ ∂∂ ∂ × + − − − × ∂ ∂ ∂ ∂
∑ ∑
r
( ),
0
;
m
i
m x q
m
B k
∞
=
×
∑
r
( ) ( ) ( )
*
* * *
* , ,*
0 0
,
q
i
m mG x i i
m x x m xq q
m m
q G
e B B
x i x
ϑϑ ξ ϑ
∞ ∞
= =
∂ ∂∇ = + + + ∂ ∂
∑ ∑
r r
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
65
( ) ( ) ( )
*
,* *
, ,* * *
0 0 0
m m mym x i i i
x m y y m yq q q
m m m
B q G
B B
x x i x
ϑ
ϑ ϑ
∞ ∞ ∞
= = =
∂∂ ∂+ + + +
∂ ∂ ∂∑ ∑ ∑
( ) ( ) ( )
*
,* *
, ,* * *
0 0 0
m m mm y zi i i
y m z z m zq q q
m m m
B q G
B B
x x i x
ϑϑ ϑ
∞ ∞ ∞
= = =
∂ ∂ ∂+ + + +
∂ ∂ ∂∑ ∑ ∑
( ) ( )
*
,* *
,* * *
0 0
q
i
m mGm z xi i
z m x xq q
m m
B q G
i e B
x y i y
ϑϑ ϑ
∞ ∞
= =
∂ ∂ ∂+ + + × ∂ ∂ ∂
∑ ∑
r
( ) ( ) ( )
*
,* *
, ,* *
0 0 0
m m mym xi i i
m x x m y yq q q
m m m
B q
B B
y y i
ϑ
ϑ ϑ
∞ ∞ ∞
= = =
∂∂
× + + + ×
∂ ∂∑ ∑ ∑
( ) ( ) ( )
*
,*
, ,* * *
0 0 0
m m mm y zi i i
m y y m zq q q
m m m
BG
B B
y y y
ϑϑ
∞ ∞ ∞
= = =
∂ ∂∂× + + +
∂ ∂ ∂∑ ∑ ∑
( ) ( )
*
* *
, ,* *
0 0
q
i
m mG xi i
z m z m x xq q
m m
q G
B j qe B
i y z
ϑϑ ϑ
∞ ∞
= =
∂∂+ + + × ∂ ∂
∑ ∑
r
( ) ( ) ( )
* *
,, *
,* * * *
0 0 0
m m my m ym x zi i i
m y yq q q
m m m
BB
B
z z z z
ϑ ϑϑ
∞ ∞ ∞
= = =
∂ ∂∂ ∂× + + + ×
∂ ∂ ∂ ∂∑ ∑ ∑
( ) ( ),*
, *
0 0
;
m mm zi i
m z zq q
m m
B
B k
z
ϑ
∞ ∞
= =
∂
× + ∂
∑ ∑
r
( ) ( )*
** * * *
0 0
1 1 q
i
m mG mi i
m q q
m m
Dq G
e D i
i x x
π
ρ γ ρ γ
∞ ∞
= =
∂∂ − ∇ = − + + ∂ ∂
∑ ∑
r
( ) ( )* *
0 0
;
m m
mi i
m q q
m m
Dq G
D j q k
i y z
∞ ∞
= =
∂∂+ + ∂ ∂
∑ ∑
rr
( )
** * * *
* **
**2 *2 *2 * * *
0
1 1 1 q
i
mG i
m q
m
e i j q k D
x y z
ρδ ρ ρ ρρ π
γ ρ γ ρ γρ
∞
=
∇ ∂ ∂ ∂∇ = − = + + ∂ ∂ ∂
∑
rr r
,
а в уравнениях (12), (15):
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
66
*
*2 * *2 * *
0 0
1 1
;
q
i
m m
G m
m
m m
Dq G i i
e D
с t c i t q t q
π ∞ ∞
= =
∂∂ ∂ = + ∂ ∂ ∂
∑ ∑
( ) ( ) ( ) ( )
*
* *
* , , ,* * *
0 0
, ;
q
i
m m
m m mG i i i
m x m y m zq q q
m o m m
e B B q B
x y z
ρ ρ ρξ ρ
∞ ∞ ∞
= = =
∂ ∂ ∂∇ = + + ∂ ∂ ∂
∑ ∑ ∑
r
{
*
* * * * *
* * **2 * *2 *2 *
1 1 1
, , ,
q
i G
x
q G
e
c c c c i x
πϑ ϑ π ϑ π ϑ ∂ ∇ = ∇ + ∇ = × ∂
r r r
* *
* * *2 * *2
0 0
1 1q
i
m m
Gm
m x y
m m
Di i
D e
q x q x c y c
ϑ ϑ
∞ ∞
= =
∂ ∂ ∂ × + + + + ∂ ∂ ∂
∑ ∑
*
* *2
0
1
;
m
z m
m
i
q D
z c q
ϑ
∞
=
∂ + ∂
∑
( ) ( ),* * *
* ,* *
0 0
q
i
m mG m xi i
m x q q
m m
Bq G
div e B
i x x
ρ ξ ρ
∞ ∞
= =
∂ ∂= + + ∂ ∂
∑ ∑
r
( ) ( ) ( ), ,
,* * *
0 0 0
;
m m mm y m zi i i
m y q q q
m m m
B Bq G
B q
i y y z
∞ ∞ ∞
= = =
∂ ∂ ∂+ + + ∂ ∂ ∂
∑ ∑ ∑
( )
** **
*
**2 *2 * * *
0
1
;
q
i
mGyx z i
m q
m
div e D
с c x y z
ϑϑ ϑπ ϑ
∞
=
∂∂ ∂= + + ∂ ∂ ∂
∑
r
( ) ( )*
* *
0 0
;
q
i
m mG mi i
m q q
m m
Eq G
e E
i t t
σ
∞ ∞
= =
∂∂ ∂ = + ∂ ∂
∑ ∑
( ) ( ) ( )* * *
* * *
0 0
,
q
i
m mG mi i
x m q q
m m
Eq G
e E
i x x
ϑ σ ϑ
∞ ∞
= =
∂∂ ∇ = + + ∂ ∂
∑ ∑
r
( ) ( ) ( )* *
* * *
0 0 0
m m m
m mi i i
y m zq q q
m m m
E Eq G
E q
i y y z
ϑ ϑ
∞ ∞ ∞
= = =
∂ ∂∂+ + + ∂ ∂ ∂
∑ ∑ ∑ .
Приравнивая коэффициенты при q1 в левой и правой частях уравнений (11), (12), (15)
и сокращая на экспоненту, получаем:
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
67
0,* *
0 0, 0, 0,* * * * *
1 1 1 1x
z z y y x
BG G G G
B B B B i
i t z i x i x i y
ε ϑ ϑ
∂ ∂ ∂ ∂ ∂+ − − − + ∂ ∂ ∂ ∂ ∂
r r
0,* *
0, 0, 0,* * * *
1 1 1 y
x y x z z
BG G G
B B B j
i x i y i y z
ϑ ϑ
∂ ∂ ∂ ∂− − − + ∂ ∂ ∂ ∂
r
0, 0,* *
0, 0,* * * *
1 1y x
y z x z
B BG G
B B k
i y z z i x
ϑ ϑ
∂ ∂ ∂ ∂ + − − − + ∂ ∂ ∂ ∂
r
* ** *
0, 0, 0, 0,* * * *
y yx x
z z y xB i B j B B k
z z z z
ϑ ϑϑ ϑε
∂ ∂∂ ∂ + + − + + ∂ ∂ ∂ ∂
rr r
( ) ( )* * 0
0 0 0* * * *
1 1 1
, ,
D
B G B k D G k
i z z
ε ϑ ϑ
ρ γ γ ρ⊥ ⊥
∂ ∂ + ∇ + = ∇ − ∂ ∂
r rr rr r
,
(73)
( )
*
0,*
0 0, 0*2 * * * *2
1 1 1
, z
z
BG
D B B G
ic t z i z ic
ρ ρ ε⊥
∂ ∂ ∂+ + ∇ + + × ∂ ∂ ∂
r
( )
*
*
0 0* *2
, 0zD G D
z c
ϑϑ ε⊥
∂× ∇ + = ∂
r
,
(74)
( )* * 0
0 0* *
1 1
, 0z
EG
E E G
i t i z
ε ϑ ϑ⊥
∂∂ + ∇ + = ∂ ∂
r
, (75)
где ( )0 0, 0,* *
, x y
G G
B G B B
x y⊥
∂ ∂∇ = +
∂ ∂
r
Будем 0 0 0, и ,B D Е
r
в свою очередь, искать в форме рядов по степеням малого
параметра:
0 0 0 0 0 0
0 0 0
; ;j j j j j j
j j j
B B D D E Eε ε ε
∞ ∞ ∞
= = =
= = =∑ ∑ ∑
r
. (76)
Тогда в нулевом приближении по ε
( ) ( )0 0
0, 0* * *
1
x
G G
B D
t xρ γ
∂ ∂= −
∂ ∂
, (77)
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
68
( ) ( )0 0
0, 0* * *
1
y
G G
B D
t yρ γ
∂ ∂= −
∂ ∂
, (78)
( )
( )0
0 0
0,* * *z
DG i
B
t zγ ρ
∂ ∂= − ∂ ∂
, (79)
( ) ( ) ( )( )
( )0*
0 0 0 0,*
0 0, 0*2 * * *
1 1
, 0z
z
BG
D B B G
ic t z i z
ρ ρ ⊥
∂∂ ∂ + + ∇ + = ∂ ∂ ∂
r
, (80)
( )0
0*
1
0
G
E
i t
∂ =
∂
. (81)
Из равенства (81) сразу же следует, что Е0
(0)=0. Если же подставить выражение для
( ) ( ) ( )0 0 0
0, 0, 0,, ,х у zВ В В из (77) – (79) в (80), то для ( ) ( ) ( ) 10 0 *
0 0D D ρ
−
≡
)
мы получим уравнение:
( )
( ) ( )
0 *
20 0* *0
0 0* * *2
0
D
D i D
z z c
γρρ ρ
∂∂ + − = ∂ ∂
)
) r )
, (82)
где введены обозначения
* * *
,
G G G
x i i j
t x y
∂ ∂ ∂= − = +
∂ ∂ ∂
r r r
.
Будем по-прежнему предполагать, что жидкость ограничена сверху свободной
поверхностью z*=0, а снизу – твердым дном z*=-Н*. Тогда граничные условия запишутся в
виде (см. (72)):
( )( ) ( )( )* * * * * *
,
0 0
, ,0, 0, , ,0, 0,
m m
i i
m m zq q
m m
D x y t B x y t
∞ ∞
= =
= =∑ ∑
( ) ( ) ( )
* *
* *
, , ,* *
0 0
1
0
m m m
i i i
m z m x m yq q q
m o m m z H
B B B
q x y
∞ ∞ ∞
= = = =−
∂Η ∂Η+ + = ∂ ∂
∑ ∑ ∑ . (83)
Приравнивая нулю коэффициент при q0 в равенстве (83), получаем, что
* * * *0 0, 0,0 0
0; 0; 0z zz z z H
D B B
= = =−
= = = , (84)
и значит (см. (76), (79))
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
69
( )
( )
*
* *
0
0 0
0 *0
0; 0
z
z H
D
D
z=
=−
∂= =
∂
)
)
. (85)
Таким образом, (82), (85) – задача на нахождение собственных значений и
собственных функций, совпадающая с уже рассмотренной задачей (43), (44).
Перейдем к нахождению амплитуды d0
(0)(x*,y*,t*) нулевого приближения:
( ) ( ) ( ) ( ) ( ) ( )0 0 0* * * * * * * * * * *
0 0 0, , , , , , , , ;D x y z t d x y t D x y z t=
)
%
( ) −0
0D
~
решение уравнения (82), удовлетворяющее граничным условиям (85) и
зафиксированное, например, соотношением:
( )
*
0 * 2
0 *
0*2
1
H
D dz
c
ρ
−
= ∫ % .
Для этого рассмотрим задачу для первого приближения по ε. Так как Е0
(0)=0, то
уравнение (75) сразу же дает, что и Е0
(1)=0. Из уравнений (73), (74), в свою очередь, следует,
что
( ) ( )
( )1
1 1 0
0 0* * * *
1 1 DG i
B D G k
i t zρ γ γ ρ⊥
∂ ∂− ∇ + = ∂ ∂
rr
( )
( ) ( )
0
0 00,* *
0, 0,* * *
1 1x
z z y x
B G G
B B i
z i x i x
ϑ ϑ
∂ ∂ ∂ = − − − − ∂ ∂ ∂
r
( ) ( ) ( )
( )0
0 0 0 0,* *
0, 0, 0,* * * *
1 1 1 y
x y x z z
BG G G
B B B j
i x i y i y z
ϑ ϑ
∂ ∂ ∂ ∂
− − − − − ∂ ∂ ∂ ∂
r
( )
( ) ( )
( )
0 0
0 0, 00,* *
0, 0,* * * *
1 1y x
y z x z
B BG G
B B k
i y z z i x
ϑ ϑ
∂ ∂∂ ∂
− − − − −
∂ ∂ ∂ ∂
r
( ) ( ) ( ) ( )
* ** *
0 0 0 0
0, 0, 0, 0,* * * *
y yx x
z z y xB i B j B B k
z z z z
ϑ ϑϑ ϑ ∂ ∂∂ ∂− − + + − ∂ ∂ ∂ ∂
rr r
( )( ) ( )( )0 0* *
0 0*
1
, , ;B G B k
i z
ϑ ϑ⊥
∂− ∇ −
∂
rr rr r
(86)
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
70
( ) ( ) ( )( )
( )1*
1 1 1 0,*
0 0, 0*2 * * *
1 1
, z
z
BG
D B B G
ic t z i z
ρ ρ
∂∂ ∂ + + ∇ + = ∂ ∂ ∂
r
( ) ( ) ( )
*
0 0*
0 0*2 * *2
1
, zD G D
ic z c
ϑϑ ⊥
∂= − ∇ − ∂
r
.
(87)
Воспользовавшись равенствами (74) – (79):
( ) ( ) ( ) ( ) ( )
( )0
0 0 0 0 0 0
0, 0 0, 0 0,* * *
1 1
, ;x y z
DG G i
B D B D B
x x x y x zγ γ γ
∂∂ ∂= = =
∂ ∂ ∂
)
)
,
векторное уравнение (86) можно записать в виде трех скалярных:
( ) ( ) ( ) ( )
( )1
1 1 0 1 0
0, 0 0, 0* *
, ,x x
Di i G
B D R D
x xγ ρ
∂= − + = ∂
) )
(88)
( ) ( ) ( )1 1 0
0, 0 0,*
,y y
i i G
B D R
x yγ
∂= − + ∂
)
(89)
( )
( )
( )
1
1 00
0, 0,*z z
Di i
B R
x zγ
∂= − − + ∂
)
, (90)
где введены обозначения:
( )
( )
( ) ( )
( )0 0*
0 0 * *0 0
0, 0* * * * *
1
, ,x
x z
D DG G
R i i D
x z z x x z
ϑ ϑ τ ϑ
γ
∂ ∂ ∂∂ ∂ = − − + ∂ ∂ ∂ ∂ ∂
) )
r) r
( )
( )
( ) ( )
( )* 0 0
0 0 * *0 0
0, 0* * * * *
1
, ,y
y z
D DG G
R i i D
x z z y y z
ϑ
ϑ τ ϑ
γ
∂ ∂ ∂∂ ∂ = − − + ∂ ∂ ∂ ∂ ∂
) )
r) r
( )
( )
( )
( )0 0
0 * *0 0
0, * * *
1
,z z
D D
R i i
x x z z
ϑ τ ϑ
γ
∂ ∂∂ = − + ∂ ∂ ∂
) )
r r
.
Представим теперь выражения (88) – (90) для ( ) ( ) ( )1 1 1
0, 0, 0,, ,x y zB B B в (87). Тогда уравнение
для ( )1
0D
)
примет вид:
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
71
( )
( ) ( ) ( ) ( )
1 *
21 1 1 0* 2 *0
0 0 0 0* * *2
D
x D D R D
z z c
γρρ ρ τ
∂∂ + − = ∂ ∂
)
) ) )r
, (91)
( ) ( ) ( ) ( ) ( )
* **
1 0 0 0*
0 0 0 0*2 * *2
, zx
R D D ijx D
c z c
ϑ ργ ρ ϑ τ ∂ = + − ∂
r) ) )r
( )
( ) ( )0 0
* * * *0 0
* * * * *
1
, z
D Di
x z z x z z z
ϑ τ ρ ρ ϑ
∂ ∂∂ ∂ ∂− − − ∂ ∂ ∂ ∂ ∂
) )
r r
( )
( ) ( )
( )0 0* * * *
2 20 * *0 0
0* * *
2
, , z
D Di
D
x z z x x z
ρ ϑ ρ ρτ τ ϑ τ ϑ τ
∂ ∂∂− + + ∂ ∂ ∂
r) )
r)r r r r
,
причем (см. (84))
( ) ( )
* * *
1 1
0 0,
0
0, 0z
z z H
D B
= =−
= =
)
, (92)
и значит (см. (85), (90)),
( )
( )
( )
( )* *
* *
1 0 *
0 *0 0
0,* * *
1
, z
z
z H
z H
D D
R i
z x z z
ϑγ ϑ τ
=−
=−
∂ ∂ ∂= = − + ×∂ ∂ ∂
) )
r r
( ) ( ) ( )
* ** *
0 0 02 2
* *0 0 0
* *2 *2z z
z Hz H
D D Di
i
z z x z
ϑ ϑ
−=−
∂ ∂ ∂× + = −∂ ∂ ∂
) ) )
.
(93)
Задача (91) – (93) – неоднородная задача Штурма-Лиувилля на собственном числе.
Для того, чтобы найти условие ее разрешимости, умножим обе части уравнения (91) на ( )0
0D
)
,
проинтегрируем полученное равенство по z* от –Н* до нуля и воспользуемся соотношениями
(82), (85):
( ) ( ) ( ) ( )
( )
* * *
0 02
0 1 0 0* * * 0
0 0 0 0 *2z
H z H
Di
D R D dz D
z
ρ ϑ
γ− =−
∂ = ∂∫
)
) ) )
,
или
( ) ( ) ( ) ( ) ( )
( )
* * *
0 0220 0 1 0 0* * * 0
0 0 0 0 0 *2
0z
H z H
Di
d D R D dz D
z
ρ ϑ
γ− =−
∂ − = ∂
∫
)
) ) )
.
Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану), 2008 (№ 5)
72
Так как выражение, стоящее в фигурных скобках последнего равенства, вообще
говоря, в нуль не обращается, то задача (91) – (93) разрешена только в том случае, если
амплитуда d0
(0)(х*, у*, t*) ≡ 0, то есть D0
(0)(х*, у*, z*, t*)≡0. Аналогично, D0
(j)(х*, у*, z*, t*)≡0
при j=1,2,… Значит, при построении нулевого приближения по q функций * * *, ,ξ π σ
r
это
нулевое приближение нельзя представлять в виде рядов по степеням малого параметра ε.
Моделирование на ПЭВМ показывает, что, при наличии течений, подводный
волновод может образоваться в тех условиях, где стандартные методы (без учета течений)
его не обнаруживают.
Литература
1. Sanford T.B. Observations of strong current shears in the deep ocean and some implications on sound rays.
– J. Acoust. Soc. America, 1974, u.56, № 4, Р.1118-1121.
2. Stallworth L.A., Jacobson M.J. Acoustic propagation in an isospeed channel with uniform tidal current and
depth change. – J.Acoust. Soc. America,1970, v.48, № 1, Р.382-391.
3. Stallworth L.A., Jacobson M.J. Sound transmission in аn isospeed ocean channel with depth-dependent
current. J.Acoust. Soc. America, 1972, v.51, № 5, Р.1738-1750.
4. Stallworth L.A., Jacobson M.J. Acoustic propagation in аn uniformy moving ocean channel with depth-
dependent sound speed. – J.Acoust. Soc. America, 1972, v.52, № 1, Р.344-355.
5. Franchi E.R., Jacobson M.J. Ray propagation in a channel with depth-variable sound-speed and carrent. –
J.Acoust. Soc.America, 1972, v.52, № 1, Р.316-331.
6. Голод О.С., Григорьева Н.С. Оценка влияния течения на поле давления монохроматического
точечного источника в однородном океане. Акустический журнал, т. 28, вып. 6, 1982.
7. Барридж Р., Вейнберг Г. Горизонтальные лучи и вертикальные моды. – В кн.: Распространение волн
и подводная акустика, М., Мир, 1980, С.76-125.
8. Бабич В.М., Булдырев В.С. Молотков И.А. Метод возмущений в теории распространения волн
(пространственно-временной метод, линейные и нелинейные уравнения). – В кн.: Теория
распространения волн в неоднородных и нелинейных средах. М., 1979, С.28-143.
9. Блохинцев Д.И. Акустика неоднородной движущейся среды. М., Наука, изд. 2-е, 1981, 206 с.
10. Вадов Р.А. Дальнее распространение звука в районе субарктического фронта//Акустический журнал,
2008, том 54, №2, С.251-261
|
| id | nasplib_isofts_kiev_ua-123456789-19217 |
| institution | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| issn | 1815-8277 |
| language | Russian |
| last_indexed | 2025-12-07T13:19:52Z |
| publishDate | 2008 |
| publisher | Науково-технічний центр панорамних акустичних систем НАН України |
| record_format | dspace |
| spelling | Гончар, А.И. Гончар, Ю.А. Голод, О.С. Григорьева, Н.С. 2011-04-22T20:06:34Z 2011-04-22T20:06:34Z 2008 Расчет акустического поля с учетом течений / А.И. Гончар, Ю.А. Гончар, О.С. Голод, Н.С. Григорьева // Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану): Зб. наук. пр. — Запоріжжя: НТЦ ПАС НАН України, 2008. — № 5. — С. 34-72. — Бібліогр.: 10 назв. — рос. 1815-8277 https://nasplib.isofts.kiev.ua/handle/123456789/19217 551.463.22 В данной статье решена задача о распространении звука в неоднородной движущейся среде в виде рядов по степеням ε. В этом случае нулевое приближение соответствует задаче о распространении звука в неподвижной неоднородной среде, а поправка, обусловленная наличием течения, выражается через нулевое приближение. Так как отношение характерного размера изменения свойств среды по вертикали к характерному размеру изменения свойств среды в горизонтальном направлении имеет порядок 10-3, то с помощью пространственно-временного лучевого метода решена задача по нахождению коэффициентов ряда по степеням малого параметра ε. У даній статті вирішено задачу про поширення звуку в неоднорідному середовищі, що рухається, у вигляді рядів по ступенях ε. У цьому випадку нульове наближення відповідає завданню про поширення звуку в нерухливому неоднорідному середовищі, а виправлення, обумовлене наявністю течії, виражається через нульове наближення. Оскыльки відношення характерного розміру зміни властивостей середовища по вертикалі до характерного розміру зміни властивостей середовища в горизонтальному напрямку має порядок 10-3, то за допомогою просторово-тимчасового променевого методу вирішено задачу по знаходженню коефіцієнтів ряду по ступенях малого параметра ε. This article solves the problem about the propagation of sound into nonhomogeneous moving medium in series in terms of power of ε. In this case zeros approximation corresponds to the problem of sound propagation into nonmoving nonhomogenenous medium and the allowance conditioned by availability of current is expressed by zeros approximation. Since the relation of characteristic dimension of change of medium properties in vertical position to characteristic dimension of change of medium properties in horizontal position has the order 10-3 the problem for estimation of ratios of series in terms of powers of small parameter ε was solved by the spatio-temporal ray-path method. ru Науково-технічний центр панорамних акустичних систем НАН України Гідроакустичний журнал (Проблеми, методи та засоби досліджень Світового океану) Расчет акустического поля с учетом течений Розрахунок акустичного поля з врахуванням течій Article published earlier |
| spellingShingle | Расчет акустического поля с учетом течений Гончар, А.И. Гончар, Ю.А. Голод, О.С. Григорьева, Н.С. |
| title | Расчет акустического поля с учетом течений |
| title_alt | Розрахунок акустичного поля з врахуванням течій |
| title_full | Расчет акустического поля с учетом течений |
| title_fullStr | Расчет акустического поля с учетом течений |
| title_full_unstemmed | Расчет акустического поля с учетом течений |
| title_short | Расчет акустического поля с учетом течений |
| title_sort | расчет акустического поля с учетом течений |
| url | https://nasplib.isofts.kiev.ua/handle/123456789/19217 |
| work_keys_str_mv | AT gončarai rasčetakustičeskogopolâsučetomtečenii AT gončarûa rasčetakustičeskogopolâsučetomtečenii AT golodos rasčetakustičeskogopolâsučetomtečenii AT grigorʹevans rasčetakustičeskogopolâsučetomtečenii AT gončarai rozrahunokakustičnogopolâzvrahuvannâmtečíi AT gončarûa rozrahunokakustičnogopolâzvrahuvannâmtečíi AT golodos rozrahunokakustičnogopolâzvrahuvannâmtečíi AT grigorʹevans rozrahunokakustičnogopolâzvrahuvannâmtečíi |