Об асимптотическом поведении пакетов линейных волн в многомерном случае. Эталонные решения
В. Г. Гневышев1^, С. И. Бадулин1,2,6
1 Институт океанологии имени П. П. Ширшова Российской академии наук.
Россия, 117997, Москва, Нахимовский пр., д. 36.
2 Новосибирский государственный университет.
Россия, 630090, Новосибирск, ул. Пирогова, д. 2.
E-mail: a avi9783608@gmail.com, 6 badulin@ioran.ru Статья поступила 12.09.2016, подписана в печать 14.10.2016.
Рассматривается классическая задача о распространении линейных волновых пакетов в диспергирующих средах. В терминах интегралов Фурье для волновых пакетов, имеющих в начальный момент времени гауссово распределение, построены асимптотики задачи Коши для случая двух пространственных измерений. Построенные асимптотики регулярны на каустике и описывают новые физические свойства пакетов — повороты в пространстве и формирование аномально слабо диспергирующего (квазидиспергирующего) волнового фронта.
Ключевые слова: линейные волны, волновая дисперсия, метод стационарной фазы, метод наискорейшего спуска, метод перевала, дисперсия волновых пакетов.
УДК: 532.591. PACS: 47.35.-i, 92.10.-H.
Введение
Задача о распространении линейных диспергирующих волн считается хорошо изученной. Мощный арсенал математических средств базируется на представлении общего решения в виде разложения по собственным функциям задачи, что позволяет решать разнообразные задачи распространения волн в неоднородной и нестационарной среде. Асимптотические методы стационарной фазы и перевала были предложены еще в XIX в. (П.-С. Лаплас, Дж. Г. Стокс, У. Кельвин) как подходы, существенным образом учитывающие волновые свойства: осциллирующий характер решений и/или локализацию этих решений в физическом пространстве или в пространстве волновых векторов. Именно эти методы дают возможность анализа основного физического эффекта, определяющего эволюцию пакетов линейных волн, — волновой дисперсии.
Случай одномерного распространения волн представляет особый интерес, поскольку позволяет сформулировать ряд важных результатов в виде явных аналитических выражений. Прежде всего, это результаты об асимптотическом поведении волновых пакетов при больших временах или на большом удалении от источника волн. Асимптотический анализ предсказывает затухание амплитуды как корня квадратного времени в системе координат, связанной с максимумом волнового возмущения, для волн произвольной физической природы за исключением особых случаев вырождения дисперсионного соотношения — обращения в нуль второй производной частоты по волновому числу. Естественное решение проблемы вырождения, включение в анализ высших дифференциалов дисперсионной зависимости приводит к решению в терминах функции Эйри и более
медленному дисперсионному затуханию амплитуды волны как A ~ t-1/3.
Обобщение асимптотических методов на многомерный случай довольно часто проводится «по аналогии», без должного учета многомерности задачи. Очень часто считается, что преобразованием системы координат можно свести задачу к одномерной и далее следовать рецепту одномерной задачи буквально. Вырождение квадратичной формы, связанной со вторым дифференциалом, трактуется в том же духе: собственный вектор, связанный с нулевым собственным значением, определяет направление, в котором задача может быть сведена к одномерной с привычными решениями в виде функции Эйри. Во многих волновых задачах подобная вольность не является существенной, а проведение математически строгого рассмотрения может рассматриваться как излишество, не имеющее отношения к решаемым физическим вопросам. Ниже мы покажем, что уже в двумерной задаче математически строгое рассмотрение обнаруживает важные физические эффекты и даже при вырождении квадратичной формы обращение к более высоким (кубическим) членам асимптотического разложения не является необходимым. При этом решением оказывается не степенная зависимость A ~ t-1/3, а A = const.
Геофизическая гидродинамика предлагает богатейший ряд задач распространения волн в среде, анизотропия которой связана с эффектами вращения и стратификации. Именно в этих задачах корректное применение асимптотических методов обнаруживает целый ряд замечательных эффектов, определяемых именно неодномерностью волновых пакетов и конечностью их объема. В настоящей работе рассмотрены основные эффекты, связанные с многомерностью, и примеры, в которых эти эффекты возникают.
Статья открывает серию исследований эффектов многомерности в волновых задачах.
1. Асимптотические методы в теории распространения точечных пакетов линейных волн
Рассмотрим классическую задачу о развитии возмущения в анизотропной, пространственно-однородной среде, допускающей существование волн малой амплитуды. Случаи пространственной неоднородности и нестационарности среды здесь не рассматриваются и будут представлены в последующих работах. Пусть источники волн и затухание отсутствуют, а начальное возмущение локализовано в некоторой области физического пространства (достаточно быстро стремится к нулю на бесконечности). Предположим, что в процессе решения задачу удалось свести к многомерному интегралу Фурье вида
П(х, ^ = (2п)-п/2
Fn(k)exp(ikx - iш(k)t) dk, (1)
где п — размерность пространства; Fn(k) — фурье-образ функции ^ в момент времени t = 0; k — многомерный волновой вектор; ш^) — дисперсионное соотношение, которое определяется типом рассматриваемых волн малой амплитуды и полагается известным. Действительная часть (1) Яе Д(х, 0 является решением физической задачи.
В общем случае интеграл (1), очевидно, аналитически вычислен быть не может. В то же время асимптотические методы могут быть эффективно использованы именно для интегралов вида (1). В современной литературе математический аппарат асимптотических методов изложен достаточно полно [1-4]. Многомерный метод стационарной фазы (см. классическую книгу Дж. Лайтхилла [4]) отражает особенности этих методов наиболее ярко. Большой параметр (время или расстояние) приводит к вычислению интеграла (1) от быстро осциллирующей функции. Старший член соответствующего асимптотического разложения связан с вкладами окрестностей точек экстремумов, стационарности фазы (отсюда название метода). Разложение в ряд Тейлора фазы в = kx - с точностью до квадратичных членов оказывается достаточным для вычисления старшего члена разложения. В том или ином виде соответствующая квадратичная форма используется в различных асимптотических подходах (например, классический метод стационарной фазы [2], метод канонического оператора Маслова [1]).
Присутствие формально большого параметра в фазе в (1) означает медленность изменения функции Fn(k), т.е. сильную локализацию волновой функции в физическом пространстве. В этом случае имеет смысл говорить о точечных волновых пакетах и точечных же асимптотиках решений (1). Разложение фазовой функции в = kx — в этом случае связывается с точкой-максимумом волнового возмущения в физическом пространстве. При вы-
рождении квадратичной формы в выделенной таким образом точке возникает сингулярность в асимптотическом выражении.
В настоящей работе мы исходим из иных физических предпосылок и предлагаем метод эталонных решений. Мы полагаем пакет локализованным в пространстве волновых векторов, а значит, рассматриваем волновые пакеты конечного размера, слабо модулированные в физическом пространстве. Такой подход также приводит к анализу поведения квадратичной формы для разложения фазовой функции в = kx — ш(^, но не относительно текущего максимума волновой амплитуды, а относительно несущей гармоники модулированного волнового пакета. Говоря об эталонных функциях и решениях мы лишь меняем способ описания кинематики пакета, переходя от волновой функции в координатном пространстве (аналогия с эйлеровым описанием в гидродинамике), к описанию волнового пакета-квазичастицы, маркированного волновым вектором несущей волновой гармоники (лагран-жев подход). Заметим, что очень близкий подход с использованием эталонных (модельных) гауссовых пакетов широко используется, например, в задачах радиофизики и волоконной оптики. Для обоих подходов исключительно важным является вопрос о пределах их применимости, который здесь не рассматривается. Физические аргументы, приведенные выше, представляются достаточными для широкого круга задач распространения волн.
Рассмотрим разложение Тейлора для частоты ш^) в окрестности точки ^ = (^, /0):
V) = 1о) + — + ш\^ — 1О) +
+ -ш''к ^ — ^)2 + ^шЦ 0 — 1о)2 + ш''1 ^ — АО)(1 — 1О) +... .
2
2
(2)
Квадратичной части разложения фазы в в (1) поставим в соответствие действительную симметричную двумерную матрицу и характеристическое уравнение [5]
A2 =
tЩ)Ш^{tk tkoloшll
Шошы Щш'ц
ёе^ — А^) = ёе1
— А 1,2 ^010ш'ы
tkoloш'kl Щш' [ — А 1,2
(3)
= 0.
Раскрывая определитель, получаем квадратное уравнение для собственных значений
А?,2 — А 1,2 Тг^ + ёе^н 0. (4)
Коэффициенты в (4) — определитель ёе^^ и след матрицы Тг^2) являются ортогональными инвариантами матрицы A2:
1^) = ЭДИ^ — (ш] )2],
Тг^2) = Щш1к + ^ {]. Отметим, что эйлерова кривизна — это след матрицы Тг^2) [6]. В силу симметричности
(5)
и вещественности матрицы A2 характеристические значения А1>2 также будут вещественны. Подчеркнем простой, но важный факт: для вычисления интеграла (1) с помощью метода стационарной фазы необязательно вычислять корни (4). Амплитуда волнового пакета определяется произведением А;А2, т.е., в соответствии с теоремой Виета, определителем A2. Асимптотика интеграла (1) /2(х, 0 « (1АхА2)-х/2 « (I A2)-1/2 имеет простой физический смысл: произведение А1А2 пропорционально физическому объему волнового пакета (А1 и А2 задают главные оси квадратичной формы A2). Мнимая единица задает сдвиг фазы ±п/4, а знак этого сдвига определяется условием затухания решения на бесконечности и может быть получен с помощью выхода в комплексную плоскость и различных способов замыкания контура интегрирования. В одном случае это приводит к интегралу Фурье, впервые вычисленному Эйлером [7]. Введение малого затухания и обход полюса в комплексной плоскости приводит к методу перевала [см., например, [4]] и, естественно, дает тот же результат.
Обобщение на случай произвольной размерности пространства п не представляет особых трудностей. Асимптотика интеграла (1) имеет вид
П(х, 0 « ^(Ы^ An]-1/2 ехр^х - ш^] х
х ехр
1п л 4 An
(6)
в пределе г ^ 1х; | ^ Волновой вектор k0 задается условием стационарности
X; дш
а условие
г дк\'
( д2ш \
\дЩ дЩ )
= 0
(7)
(8)
определяет достаточность разложения (2) с точностью до квадратичных по малым отклонениям волнового вектора членов. В (6) 8§п An — сигнатура матрицы, разность числа положительных и отрицательных собственных значений матрицы An.
Читается асимптотика (6) следующим образом. При больших г, х; и конечном отношении х;/г разрешим лучевые уравнения (7) относительно волновых чисел Щ = Щ(х;/г). Найденное значение волнового вектора ^(х/г) в (6) характеризует изменение масштаба (длины, периода) несущей волнового пакета при его распространении. Заметим, что нахождение ^(х/О, т.е. обращение лучевых уравнений может оказаться нетривиальным даже для волн в однородной и изотропной среде, не говоря о более общих ситуациях.
Асимптотику (6)-(8) можно назвать дисперсионной: при больших временах и отличном от нуля определителе (8) каждая пространственная координата дает затухание г-1/2, а амплитуда волны зату-
хает обратно пропорционально «объему» волнового пакета в п -мерном пространстве, т.е. как г-п/2. Говоря об «объеме» волновых пакетов, следует помнить об условности этого термина, а именно мы говорим об асимптотическом поведении точечного волнового пакета, распределение которого описывается в пространстве волновых чисел дельта-функцией. Отсюда в (6) появляется амплитуда фурье-гармоники Рп^0) в точке стационарности (7). Ниже мы используем разложение (2) для альтернативного метода вычисления асимптотики интеграла (1). Отличие сводится фактически к другому определению вектора ^ = (Щ0, /0).
Переход к задаче о волновых пакетах конечного объема кажется очевидным — следует ввести в рассмотрение пространственно-временную модуляцию волнового пакета и учесть возникающие эффекты. Один из способов сделать это — введение формально малого затухания для только что рассмотренных точечных волновых пакетов. Окончательное решение получается в пределе бесконечно малого затухания. Проблема с предельным переходом оказывается непреодолимой на каустиках, связанных с сингулярностью асимптотического выражения, когда определитель матрицы 8§п An обращается в нуль. В таких точках у известных асимптотических методов, в том числе и у метода канонического оператора Маслова, возникают проблемы, которые решаются построением нового асимптотического решения. Эти проблемы решаются по аналогии с одномерной задачей учетом кубического члена в разложении частоты. Возникающие таким образом решения в виде функций Эйри сшиваются с асимптотическими решениями вдали от каустики. Задача оказывается довольно трудоемкой.
Ниже мы покажем, как конечность физического объема пакета может быть учтена путем «размазывания» решения в пространстве волновых векторов. Ключевым отличием нашего подхода является то, что мы не требуем выполнения условия стационарности (7). Это позволяет обойтись без решения лучевых уравнений, необходимого для точечных пакетов. Кроме того, мы избегаем проблем на каустике при обращении в нуль определителя An (8). Это достигается ценой ограничения на вид начальных условий — мы ограничиваемся эталонными формами пакетов. На наш взгляд, такой физически наглядный прием не может существенно исказить эффекты трансформации волновых пакетов малого, но конечного объема, прежде всего эффект поворота пакетов.
2. Физические эффекты распространения волновых пакетов конечного размера
Пусть в начальный момент времени г = 0 волновой пакет имеет гауссово распределение
/2 (х, У, г = 0) = = 1
= ДхДу ехр
1к0х + й0у —
У
2(Дх)2 2(Ду)2
(9)
2
х
и фурье-образ
F2(h, l) = exp
(k - hp)2 (l - lo)2
' 2(Ah)2 2(Д1)2
(10)
где Дк = (Дх)-1, Д/ = (Ду)-1 характеризуют ширину пакета в координатном и фурье-пространствах. Соответствующие решения для распределений (9), (10) будем называть эталонными решениями. В пределе Дк ^ 0, Д/ ^ 0 распределение (10) стремится к 5-функции, однако интеграл (1) будет принципиально отличаться от представленных нами выше точечных асимптотик (6)-(7). При этом, что очень важно, у нас не возникает необходимости в предельном переходе!
Для подынтегрального выражения в (1) построим симметричную комплексную двумерную матрицу B2 квадратичного приближения
B
2=
i + tt (Д^ у
itДhДlш
hh
hl
itДhДlw'ft^ 1 + it^l)2y\
(11)
exp
-1 (Bx, x) - i(x, g)
dx =
= (2я-)л/2( det B) 1/2 exp
-1 (B-1i, i)
. (13)
Здесь B 1 — обратная матрица. Для решений, затухающих на бесконечности, имеем
VdetB-7 =| detB |-
1/2
exp
v
- 2 arg B
Для начальных условий (10) действительная часть матрицы B строго положительна (Яе B >0), что обеспечивает абсолютную сходимость интеграла. В нашем случае в (13) появляется безразмерный аргумент квадратичной формы ¿(¿х, ¿у)
gx = Дk(x - Cgr.xt), gy = Д1(у - Cgr.yt)
(14)
которая отличается от матрицы A2 в (3) единицей на диагонали. Еще одно важное отличие: в A2 мы использовали несущее волновое число только из соображений размерности, на последующие выкладки этот масштаб не оказывал никакого влияния. Введение новых пространственных масштабов ширины пакета Дк, Д/ снимает подобное вырождение, присутствующее в точечных решениях.
Поведение эталонных решений по-прежнему определяется инвариантами квадратичной формы (11), однако сами эти инварианты имеют более сложный вид
ёе! B2 = 1 + И[(Дк)2 ш'и + (Д/) Ч/ ] +
+ ;2(Дй)2(Д/)2[«/)2 - <«], (12) Тг^) = И[(Дк)2ш''к + (Д/)2^'/ ] + 2.
Эти зависимости не сводятся к произведениям комбинации вторых производных дисперсионной зависимости (ср. (4)), что подразумевает отсутствие масштаба пакета. Соответственно амплитуда А не обращается в нуль при вырождении квадратичной формы, определяемой единственно волновой дисперсией. Более общий вид аргумента экспоненты, содержащей линейные члены, принципиально не усложняет задачу. Преобразование Фурье экспоненты квадратичной формы (см., например, [2]) дает
определяющим отклонение от траектории пакета с волновым вектором несущей гармоники k0; cgr. x = Ц., cgr. y = У в (14) — компоненты групповой скорости.
Следует отметить важные физические отличия от метода стационарной фазы: переход к характеристикам типа (14) g(gx, gy) = 0 обеспечивается наличием большого параметра, «бесконечно большого» времени или расстояния (см. (6)). В нашем подходе эталонных решений малость аргумента g(gx, gy), определяющая отклонение от характеристики, связана с узостью спектра пакета (Дk ^ 0, Д1 ^ 0). Явные ограничения для координат и времени в данном подходе отсутствуют. В этом смысле наш подход можно считать более общим, чем привычные асимптотические методы.
Для определителя матрицы B в окрестности волнового вектора несущей k0 = (h0, l0)
D2 = |detB|2 = [1+t2fy-VxVy]2+t2(vx+Vy)2, (15) где для производных дисперсионного соотношения использованы обозначения
fix = Д)2и'ы, Vy = Д)2Цц, Vxy = ДШ^.
Волновая функция при этом имеет вид
f2(x, y, t) = ^2(x, y, t)cos [^(x, y, t)] (16) с огибающей (при t = 0 — гауссиана) 1
^2(x, y, t) = ~d X
xexp<^
{¿2 [ixo+t2«+v2y))+i2y(1+t2(v2xy+V))-
- 2gx iyt2Vxy (fx +
Vy) }
(17)
и фазовой функцией ф(х, у, 0 = к0х + /0у — ш(к0, /0)t — в2^) — ш2(х, у, £).
(18)
Поправка к несущей частоте ш(к0, /0) задается квадратичной формой
^2 (х, у, t) = 2П-2 ¿Хх{ ^Ну (Нху — Цх Ну) — tн■х) + + ¿у^Нх^ху — НхНу) — +
(1 + Нху — НхНу)) , (19) а сдвиг фазы также зависит от времени
(Нх + Ну)
6»2 = 2 arctg
(20)
1 + Нху — НхНу)_ Решение (14)-(20) справедливо при временах t ^ 0, в момент t = 0 совпадает с (8) и описывает дисперсионное расплывание, сопровождающееся поворотом пакета конечной ширины в физическом пространстве
на некоторый угол. Легко убедиться в том, что скорость вращения пакета конечна и асимптотически стремится к нулю при г ^
Числитель показателя экспоненты (17) можно переписать в следующей форме:
С ($х, Су ) = —г2[(^ху У Мх )2 + (£у Мху—£х Му)2] —^Хс—^У,
(21)
наглядно показывающей характер трансформации начальной (не зависящая от времени часть) формы пакета со временем. Зависящая от времени часть представляет собой эллипс с осями мху — мх = 0, СУМХУ — СхМу = 0. При больших временах она становится определяющей. Интересно, что первоначально изотропное распределение вовсе не означает отсутствия эффекта поворота пакета: при г ^ форма пакета определяется как его начальной шириной (ДЩ, Д1), так и свойствами волновой дисперсии (мх, му,Мху). Построенное решение (16)-(20) формально не требует решения лучевых уравнений, в отличие от рассмотренных выше точечных асимптотик. Волновой вектор несущей гармоники k = (Щ0,10) фиксирован, что кардинальным образом упрощает рассмотрение.
В важном частном случае ш'^ = 0 и изотропного дисперсионного соотношения направление групповой скорости совпадает с одной из осей эллипса. При этом поворота пакета не происходит. Матрица B и обратная ей B-1 диагональны, и двумерное решение (17) распадается на произведение двух одномерных диспергирующих решений
ф2(х, у, г) = Ф(х, г)Ф(у, г),
(22)
ф(х, г) =
1
Ф(у, г) =
(1 + г2мХ) 1
(1 + г2мУ)
1/4
1/4
ехр
ехр
х - с.
ёгх
о2
2(Дх)0(1 + м|г2)
(у — %г. у 0 2 2(Ду)0(1 + м2^2)
Поправка к фазе распадается на сумму арктангенсов [8]
Мх + Му 1 — г2Мх Му _
Мг) = — 2 аг^
= 2 [аге1е(г(ДЩ)2шЩЩ) + аге1е(г(Д/)2ш;'/)] . (23)
При г ^ то (23) дает предельные значения фазы в 2 (то) = ±п/2 или 0 в зависимости от знака производных ш//. В п-мерном случае разность числа положительных и отрицательных вторых производных отвечает сигнатуре вещественной матрицы. Таким образом, для в2 в (23) мы имеем корректный предельный переход к асимптотике точечных пакетов. Заметим, что многомерные решения типа (22,) (23) хорошо известны [9].
Качественные различия между асимптотиками точечными и для пакетов конечного размера обнаруживаются в случае выполнения условия
Для точечных асимптотик эта ситуация не может быть разрешена разложением фазы в окрестности точки стационарной фазы до квадратичных членов. Разложения высшего порядка приводят к решениям в терминах функций Эйри. При рассмотрении пакетов конечного размера никакой проблемы не возникает: в показателе экспоненты огибающей Ф2(х, у, г) старшие по времени члены сворачиваются в полный квадрат, а (23) описывает вырождение эллиптических линий постоянной фазы в прямую линию в пределе г^ то. В одном направлении размер пакета стремится к предельному значению, в перпендикулярном же бесконечно вытягивается. Решение в этом случае можно назвать квазидиспер-гирующим, поскольку дисперсионное расплывание в двумерном пространстве происходит так же, как и в одномерном случае. Из (17)-(20) имеем
& = хс\[\ш1\ — ус\[ш\ 5ЙП [шЩ/ (ш'Щщ + )] ,
хс — х xt, ус — у yt,
Ф2(х, у, г) =
В1/2
(25)
х ехр
г2 (ДЩ)2 (Д/)2 |мх + му 1С2 — + Ф
2В2
Ды = (шЩ/)2 — ш'Щши = 0 при шЩ/= 0.
(24)
Изолинии решения (25) — эллипсы, эксцентриситет и ориентация главных осей зависят от времени. Условие (24) обеспечивает вырождение этих эллипсов в семейство параллельных прямых при г ^ +то. Любое сечение пакета имеет гауссову форму, но ширина изменяется со временем по-разному: в направлении малой оси стремится к постоянному значению, по большой же оси происходит линейное по времени расплывание, что и обеспечивает затухание пакета, аналогичное одномерному случаю В-1/2 ~ г-1/2. Видно, что наше рассмотрение, учитывающее наипростейшим способом конечный размер волнового пакета, приводит к качественно новым эффектам (поворот пакета и изменение его формы) в сравнении с одномерным случаем. В то же время наш подход снимает затруднение, связанное с вырождением определителя (24), описывающего дисперсию точечных решений без учета высших членов асимптотических разложений.
Отметим, что для изотропного закона дисперсии ш(Щ, /) = П(д), д = (Щ2 + /2), условие каустики (шЩ)2 — ш'Щш'и = 0 сводится к + 2дП") = 0 и реализуется только для линейной зависимости П = (ад1/2 + Ь) (а, Ь — произвольные постоянные).
3. Пример. Пакет волн Россби конечного размера
В заключение представим пример, иллюстрирующий приведенный выше асимптотический анализ эволюции двумерных пакетов линейных волн конечного размера. Рассмотрим линейные волны Россби в океане с постоянной стратификацией, диспер-
х
сионное соотношение для которых имеет вид
" = (26)
Здесь ¡3 = 2Псоз#о/й — градиент параметра Ко-риолиса в северном направлении на широте во', й — радиус; П — частота вращения Земли; а2 = 1г2п2{}2/(Ат2Н2) — квадрат собственного значения краевой задачи по вертикальной координате; А^2 = —^ — частота Брента-Вяйсяля, полагаемая постоянной; Н — глубина океана; & — компонента волнового вектора в зональном, I — в меридиональном направлениях; п = 0, 1,... — порядковый номер моды. Для баротропной моды п = 0 масштаб а обращается в нуль. Условие (24) вырождения волновой дисперсии точечных пакетов в терминах безразмерных волновых чисел к = к/а, 1 = 1/а приводит к уравнению гиперболической лемнискаты (частный случай лемнискаты Бута), хорошо известной в теории плоских алгебраических кривых:
¿4 + 2*2/2 + /4-ЗЛ2 + /2 = 0.
(27)
Эта универсальная для рассматриваемого типа волн кривая показана на рис. 1, Ь вместе с дисперсионным соотношением (26) для бароклинных волн Россби (рис. 1 ,а). Видно, что интересующие нас эффекты реализуются, например, для достаточно высокочастотных волн (векторы кь к2) в окрестности максимума безразмерной частоты штах = 1/2.
1
м/н/г.^ 0.5
/ 0.45/
1 к 1
1 11.5 -1 I -0-5 0
-0.5
оц^оз/ -1
Н/2,«^-
____ -0 5
*
-2 \ -1.5 -1 -0-5 у
у/-0.5 -1-
Рис. 1. а — Изочастоты дисперсионного соотношения (26) в безразмерных координатах к = (&,/). Показаны изочастоты ш = 0.45, 0.4, 0.35, 0.3 и вол-
новые векторы к; = (— 1, /1), 1\ = у \/ 17/4 — 3/2;
кг = (— \/2, /9) для примеров рис. 2-4. б — Условие вырождения дисперсии (24) для дисперсионного уравнения бароклинных волн Россби (26). Показаны векторы для примеров рис. 2-4
Лемниската разбивает всю фазовое пространство волн Россби на две области. Внутри кривой (наиболее высокочастотные волны) определитель В2 положителен, т. е. дисперсия положительна. Внешняя область отвечает отрицательной дисперсии. Знак дисперсии принципиальным образом определяет динамику слабонелинейных диспергирующих волн [11]. В частности, формирование циклонов и антициклонов, которые часто интерпретируются как нелинейные волны Россби, естественно связать с соответствующими областями волновых чисел,
которые выделяются лемнискатой (циклоны — положительная дисперсия — внутри лемнискаты, антициклоны — отрицательная дисперсия — внешность лемнискаты). Корректное рассмотрение этого вопроса требует анализа трехмерных волн Россби, что будет проделано нами в последующих работах.
Рассмотрим эффект квазидисперсии пакета двумерных волн Россби с несущим безразмерным волновым вектором к; = ^-1, [(17/4)1/2 - 3/2] ^ . Анизотропия дисперсионного соотношения делает возможным проявление этого нетривиального эффекта даже для первоначально изотропной формы волнового пакета. На рис. 2 показаны изолинии огибающей пакета в разные моменты времени. При малых временах = 1, £ = 10) эволюция сводится к медленному почти изотропному расплыванию пакета. При больших временах = 100, 1= 1000) рас-плывание продолжается по закону близкому к линейному (внимание на масштаб осей). Расплывание сильно анизотропизуется: пакет вытягивается вдоль некоторого направления, которое асимптотически стремится к предельному на бесконечно больших временах. Очевидно, что качественно близкий эффект сильной анизотропизации первоначально симметричной формы волнового пакета будет наблюдаться и для близких волновых векторов. Иными словами, этот эффект является структурно устойчивым, грубым для конечной области пространства волновых чисел и рассматриваемого типа волн (волн Россби), а значит, физически значимым. В рамках предложенного нами подхода вопрос о структурной устойчивости решается элементарно. Заметим,
1 Ус а \ хс
1 1 1,5 0 1 0 5 1 1
4- б
3-
1- 1 1 хс
4! з 1г 1 0 1- 1 2 1 3 4
3
4-
400 Ус г
300
200
100
100 200 300 400
Рис. 2. Форма пакета волн Россби с безразмерным вол-
новым вектором ¿1 = ( — 1, и \/17/4 — 3/2 ) в разные
моменты времени. Изолинии отвечают уровням 0.75 и 0.5 от максимума амплитуды волны в рассматриваемый момент времени t= 1, 10, 100, 1000
Рис. 3. Форма пакета волн Россби для безразмерно-
го волнового
вектора к2 = \/Зз) /2
в разные моменты времени t = 0, 5, 25, 125. Изолинии отвечают уровням 0.75 и 0.5 от максимума амплитуды
Ь Ус
б Ус____
-1 / -ОЙ 0 / 0.5/ 1
V "°'5
----л.
20 Ус
г
10
Рис. 4. Форма пакета волн Россби для безразмерного волнового вектора кг = ^^—5 + \/Зз) /2
в разные моменты времени t = 0, 5, 25, 125. Начальный контур пакета ориентирован под прямым углом к таковому в случае рис. 3. Изолинии отвечают уровням 0.75 и 0.5 от максимума амплитуды
что при построении точечных асимптотик волновых пакетов вблизи каустик, т. е. при рассмотрении высших приближений метода стационарной фазы, вопрос о грубости соответствующих решений во многих случаях остается открытым.
Рис. 5. Сечения волнового пакета на уровне половины максимума, отвечающие случаю рис. 3
к2 = У (-5 + N/33) /2^ для времен t = 0 -
сплошная, t = Ъ — штрих, t = 2Ъ — штрихпунктир, t= 125 — пунктир
В завершение приведем еще две иллюстрации поведения пакета волн Россби при критическом значении безразмерного волнового вектора к[ = (-21/2, [(331/2 - 5)/2]). Аналогично рис. 2 на рис. 3, 4 показаны изолинии огибающей волнового пакета в моменты ^ = 0, 5, 25, 125 для начальных распределений в виде эллипса. На рис. 3 эллипс с отношением главных осей 2 вытянут вдоль оси х и в процессе эволюции меняет свою форму и ориентацию. Эволюция происходит на существенно меньших временах, чем в случае рис. 2, причем поворот пакета происходит гораздо быстрее, чем последующее квазидисперсионное расплывание. Эволюция пакета с теми же параметрами, но вытянутого вдоль оси ординат (рис. 4), происходит аналогичным образом.
Серию примеров завершим рис. 5, где показаны сечения волнового пакета на уровне половины максимума, отвечающие случаю рис. 3 для времен ^ = 0, 5, 25, 125.
Заключение
Мы рассмотрели приближенные решения для некоторых задачи распространения волн в многомерной среде с учетом конечности физического объема волновых пакетов. Нельзя рассматривать наш подход как критику хорошо разработанных асимптотических методов в теории распространения линейных волн. Известные асимптотические методы (стационарной фазы, перевала и т. п.) неявно содержат предположение о бесконечно малом объеме волнового пакета. Проблема сингулярности этих асимптотик в случае вырождения дисперсионных характеристик хорошо понята специалистами и совершенно определенно связывается именно с «то-чечностью» пакетов [3]. Мы вводим модуляцию
пакетов, т. е. рассматриваем пакеты конечного объема, с самого начала. В такой постановке серьезные трудности, возникающие в классических асимптотических методах, оказываются несущественными. Возможный спор о справедливости того или иного подхода сталкивается с проблемой физических реалий: насколько справедливо рассматривать пакеты как точечные в пространстве координат или же в пространстве импульсов. Для волн в океане, на наш взгляд, разрешение спора возможно на пути примирения обеих точек зрения: диапазон возможных условий настолько широк, что следует иметь в виду все физические эффекты, предсказываемые этими подходами.
Наибольший выигрыш при использовании нашего весьма грубого, рутинного подхода мы получаем в многомерном случае, где эволюция волнового пакета обретает дополнительную степень свободы — возможность изменения формы и ориентации огибающей волнового пакета. Мы считаем излишним обсуждать богатство физических эффектов, связанных с этим довольно простым фактом. Сильная анизотропия волновых движений в океане порождает множество теорий, призванных ее объяснить. Мы показываем, что одно из таких объяснений связано с собственной кинематикой пакетов линейных диспергирующих волн. Соответствующие примеры для волн разных типов мы планируем представить в следующих публикациях.
Работа выполнена при финансовой поддержке Российского научного фонда (проект № 14-2200174).
Список литературы
1. Маслов В.П. // Асимптотические методы и теория возмущений. М.: Наука, 1988.
2. Федорюк М.В. // Асимптотика: интегралы и ряды. М.: Наука, 1987.
3. Уизем Дж. // Линейные и нелинейные волны. М.: Мир, 1972. (Whitham G.B. // Linear and Nonlinear Waves. Wiley, 1974.)
4. Лайтхилл Дж. // Волны в жидкостях. М.: Мир, 1981. (Lighthill J. // Waves in Fluids. Cambridge University Press, 1978.)
5. Беклемишев Д.В. // Курс аналитической геометрии и линейной алгебры. М.: Физматлит, 2006; 2007.
6. Погорелое А.В. // Дифференциальная геометрия. М.: Наука, 1974. (Pogorelov A.V. // Differential Geometry. 1st edition. Gronningen, The Netherlands: P. Noord-hoff N.V., 1954.)
7. Лаврентьев М.А., Шабат Б.В. // Методы теории функции комплексного переменного. Изд. 5-е, испр. М.: Наука, 1987.
8. Абрамовиц М., Стиган И. // Справочник по специальным функциям. М.: Наука, 1979. (Abramowitz M., Stegun I. // Handbook of Mathematical Functions with Formulas, Graphs and Mathematical Tables. National Bureau of Standards. Applied Mathematics Series. 55. Tenth printing. December 1972.)
9. Пелиновский Е.Н., Слюняев А.В., Талипова Т.Г., Хариф К. // Изв. вузов. Радиофизика. 2003. XLVI, № 7. С. 499.
10. Педлоски Дж. Геофизическая гидродинамика. Т. 2. М.: Мир, 1984. (Pedlosky G. // Geophysical Fluid Dynamics. New York: Springer-Verlag, 1987.)
11. Кадомцев Б.Б. // Коллективные явления в плазме. М.: Наука, 1988.
On the asymptotics of multidimensional linear wave packets: reference solutions V.G. Gnevyshev1a, S.I. Badulin12b
1 Shirshov Institute of Oceanology, Russian Academy of Sciences. Moscow 117997, Russia.
2 Novosibirsk State University. Novosibirsk 630090, Russia. E-mail: a avi9783608@gmail.com, b badulin@ioran.ru.
The classic problem of linear wave-packet propagation in a dispersive medium is considered. Asymptotic equations of the Cauchy problem for two-dimensional Gaussian wave packets are constructed in terms of Fourier integrals. These asymptotic solutions are regular at the caustics and describe new physical features of wave-packet propagation: rotation in space and formation of a wave front with anomalously slow dispersion (quasi-dispersive).
Keywords: linear waves, wave dispersion, method of stationary phase, method of steepest descent, saddlepoint method, wave packet dispersion. PACS: 47.35.-i, 92.10.-H. Received 12 September 2016.
English version: Moscow University Physics Bulletin. 2017. 72, No. 4. Pp. 415-423.
Сведения об авторах
1. Гневышев Владимир Григорьевич — канд. физ.-мат. наук, ст. науч. сотрудник; тел.: (499) 124-75-65, e-mail: 978-36-08@mail.ru, avi9783608@gmail.com.
2. Бадулин Сергей Ильич — доктор физ.-мат. наук, вед. науч. сотрудник, зав. лабораторией; тел.: (499) 124-75-65, e-mail: sergei.badulin@gmail.com, badulin@ioran.ru.