Нелинейная стадия модуляционной неустойчивости тема автореферата и диссертации по физике, 01.04.02 ВАК РФ

На правах рукописи

Гелаш Андрей Александрович

Нелинейная стадия модуляционной неустойчивости

01.04.02 Теоретическая физика

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

4 ,Л 2014

Новосибирск 2014

005556393

005556393

Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Новосибирский национальный исследовательский государственный университет».

Научный руководитель: Захаров Владимир Евгеньевич

доктор физико-математических наук, академик РАН

Официальные оппоненты: Матвеев Владимир Борисович

доктор физико-математических наук Институт ИБМП ФГАОУ ВПО СПб ГУАП г. Санкт-Петербург, главный научный сотрудник

Габитов Ильдар Равильевич

кандидат физико-математических наук профессор университета Аризоны, г. Тусон, США (The University of Arizona, Department of Mathematics)

Ведущая организация: Федеральное государственное бюджетное

учреждение науки Математический институт им. В.А. Стеклова Российской академии наук

Защита диссертации состоится 26 декабря 2014 г. в 11 часов 30 минут на заседании диссертационного совета Д 002.207.01 при Федеральном государ- ствен-ном бюджетном учреждении науки Институт теоретической физики им. JI. Д. Ландау Российской академии наук по адресу: 142432, Москов- екая обл., г. Черноголовка, просп. Академика Семенова, д. 1-А, Институт теоретической физики им. JI. Д. Ландау РАН.

С диссертацией можно ознакомиться в библиотеке Института теоретической физики им. Л .Д. Ландау РАН и на сайте www.itp.ac.ru.

7 ¿Г //

Автореферат разослан _ __2014 г.

Учёный секретарь диссертационного совета,

доктор физико-математических наук Гриневич Пётр Георгиевич

Общая характеристика работы

Актуальность темы.

Фокусирующее нелинейное уравнение Шредингера (НУШ) является одной из важнейших моделей для изучения распространения квазимонохроматических волн в слабонелинейных средах. В частности НУШ описывает волны на глубокой воде [1], волны в оптическом волокне [2], Ленгмюров-ские волны в плазме [3] и Бозе-конденсат с притяжением [4].

НУШ имеет простое решение - монохроматическую волну с частотой, зависящей от амплитуды - конденсат. Конденсат неустойчив относительно малых возмущений. Это явление носит название модуляционная неустойчивость (МН). Линейная стадия МН представляет собой экспоненциальный рост начального возмущения. Но как выглядит поведение решения на больших временах - нелинейная стадия? Для пространственной размерности Б=2, 3 ответ известен - МН приводит к формированию сингулярных решений за конечное время - коллапсов. В размерности Э=1 коллапсы запрещены (см. |5]), однако теперь развитие МН может приводить к формированию волн экстремально большой амплитуды (т.н. "волн-убийц") [6, 7, 8]. Поэтому, изучение последствий развития модуляционной неустойчивости на больших временах является задачей, важной с практической точки зрения, ключевой для создания теории волн-убийц в океане и теории экстремальных волн в оптическом волокне.

НУШ является моделью первого приближения. Для поверхности жидкости эта модель описывает слабонелинейный волновой пакет с крутизной не более чем 0.15 [9]. В нелинейной оптике его приложения также ограничены [4]. В настоящее время разработано большое количество моделей, обобщающих НУШ. Для поверхностных волн это уравнения Дысте [10, 11], для волн в оптическом волокне это уравнения включающие третью производную по времени и более сложные формы нелинейности (см. [12, 13]). Также волны-убийцы в океане изучались с помощью численного моделирования в рамках точных уравнений Эйлера для потенциального течения со свободной границей [6, 14]. Однако, многие особенности решений НУШ сохраняются на качественном уровне и в более точных моделях, что делает исследование интегрируемого НУШ чрезвычайно актуальным. К тому же, в настоящее время уровень экспериментальных установок позволяет создавать условия, необходимые для выполнения НУШ. При этом анали-

тические решения НУШ на фоне конденсата воспроизводятся в лабораторных экспериментах по изучению волн в оптоволокне и водных бассейнах с превосходной точностью [16, 15, 17].

С 1972 г. известно, что НУШ является системой, которая может быть полностью проинтегрирована с помощью метода обратной задачи рассеяния [18]. С тех пор сотни статей и несколько монографий были посвящены этой теме (см.,например, [19, 3, 20, 21). ). Среди них можно найти и работы посвящённые применению метода обратной задачи к фокусирующему НУШ на фоне конденсата [22], [23], [26]. Тем не менее, важнейший вопрос о том какова нелинейная стадия модуляционной неустойчивости был не решен.

НУШ обладает большим количеством аналитических решений. Поэтому естественно надеяться, что нелинейное развитие модуляционной неустойчивости описывается некоторыми из них. Нас интересует неустойчивость, возникающая из локализованных возмущений конденсата. Исторически, первое подобное решение было найдено Перегрином в 1983 [27] (т.н. солитон Перегрина). Это полностью гомоклинический сценарий развития модуляционной неустойчивости. Развитие неустойчивости начинается с конденсата при 4 -»■ -оо, и возвращается к нему при 4 -* +оо. При этом максимальное значение амплитуды решения в три раза превышает амплитуду невозмущённого конденсата. В последнее время это решение привлекает к себе много внимания [28], [29]. Было экспериментально подтверждено его воспроизведение в гидродинамике и оптике [16, 17]. В 1985 был найден двухсолитонный аналог этого решения [30], а в настоящее время В. Б. Матвеевым, П. Дюбардом и их соавторами активно изучаются мульти-Перегриновские решения (см. [31], [32]).

Однако, результаты численного моделирования [25], [11] демонстрируют формирование осциллирующих структур - бризеров. Поэтому, мы полагаем, что гомоклинический сценарий развития модуляционной неустойчивости является не самым правдоподобным. К тому же, решение Перегрина является весьма специальным сценарием развития МН и существует необходимость в построении теории, позволяющей описывать эволюцию произвольного начального возмущения.

Цель диссертационной работы. Целью работы является описание нелинейной стадии модуляционной неустойчивости с помощью точных решений НУШ на фоне конденсата. А также, построение общего /У-

солитонного решения НУШ на произвольном фоне с помощью метода одевания, исследование солитонных решений НУШ на фоне конденсата и обобщение результатов на векторный случай - систему Манакова.

Научная новизна, заключается в следующих основных результатах, выносимых на защиту:

1. С помощью метода одевания построено общее N - солитонное решение НУШ на фоне конденсата. Применение униформизации позволило получить явные выражения для солитонных решений в достаточно простом и ясном виде, удобном для анализа. Солитонные решения НУШ на фоне конденсата описаны с единой точки зрения. Введено понятие регулярного решения, т.е. решения которое не меняет фазу конденсата и тем самым может описывать модуляционную неустойчивость, развивающуюся из локализованных возмущений.

2. Обнаружен широкий класс 2ДГ - солитонных (суперрегулярных) решений, которые описывают нелинейную стадию модуляционной неустойчивости, развивающуюся из малого локализованного возмущения. Эти решения формируют бесконечномерное функциональное пространство. Мы полагаем, что любое неустойчивое возмущение конденсата может быть в пределе аппроксимировано нелинейной комбинацией суперрегулярных решений.

3. Исследованы вырожденные солитонные решения НУШ на фоне конденсата, в частности, описаны вырожденные суперрегулярные решения.

4. Построено ЛГ-солитонное решение НУШ на фоне солитона Перегрина. Этот результат позволяет избежать вычисление предела всякий раз при исследовании новой нелинейной комбинации солитонов и тем самым демонстрирует преимущества метода одевания. В дальнейшем это позволит исследовать взаимодействие солитона Перегрина и суперрегулярных решений.

5. Показана возможность наблюдения суперрегулярных решений в экспериментах с гравитационными волнами на поверхности жидкости. Указаны параметры при которых целесообразно проводить данные эксперименты. Представлены результаты, демонстрирующие хорошее согласие теории и эксперимента.

6. Построено обобщение метода одевания для НУШ на фоне конденсата на векторный случай - систему Манакова. Описано общее одно-солитонное решение. Установлено, что выбор пары полюсов вблизи разреза с противоположными значениями параметра а не приводит к появлению суперрегулярных решений, за исключением тривиального обобщения решений скалярного случая.

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

Апробация диссертации. Результаты диссертации докладывались: на XX и XXI научных сессиях Совета по нелинейной динамике РАН (Москва, 2011 и 2012), на всероссийской конференции молодых ученых в рамках XII научной школы «Нелинейные волны - 2012» (Нижний Новгород, 2012), на международной конференции «SOLITONS, COLLAPSES AND TURBULENCE» (Новосибирск, 2012), на международной конференции «American mathematical society meeting» (Тусон, США, 2012), на международной конференции «European Geosciences Union General Assembly» (Вена, Австрия, 2013), на международной конференции «Nonlinear Waves Theory and Applications» (Пекин, Китай, 2013), на международной конференции «SOLITONS, COLLAPSES AND TURBULENCE» (Черноголовка, 2014), на международной конференции «Наука будущего» (Санкт-Петербург, 2014), а также на теоретическом семинаре в ИТФ им. JI. Д. Ландау РАН и теоретическом семинаре в ИЯФ СО РАН.

Публикации. Результаты диссертации опубликованы в трёх статьях в научных журналах, входящих в перечень ВАК. Список работ приведен в конце автореферата.

Структура диссертации.

Диссертация состоит из введения, четырёх глав, заключения и списка литературы.

Краткое содержание работы Введение

Во введении обоснована актуальность рассматриваемой темы, научная новизна исследований, а также сформулированы цели и приведены основные результаты работы.

Глава 1. Метод одевания для нелинейного уравнения Шредингера

В первой главе строится схема метода одевания для НУШ, которое мы записываем в следующей форме:

т - ^и - (М2 -\а\2)Ч> = 0- (1)

с неубывающими граничными условиями \Л\2 при х -*■ ±оо. Где А -

произвольная константа, которую без ограничения общности можно считать действительной. Уравнение (1) является условием совместности следующей переопределенной системы на матричную функцию Ф [18]

Лф _

(2)

ох

= + (3)

си

Где

_ 1 / 1<р12 - А2 ч>х и = 1А + и, ЛУ=- 1 , , ,2 „о

2\ VI ~Ы 2 + А2

('•".)• -и :)■

(4)

Здесь А - спектральный параметр. Идея метода одевания [34] состоит в том, что мы выбираем некоторое решение ¡ро НУШ (1), для которого известно соответствующее матричное решение Фо системы уравнений (2), (3). Затем выбираются N комплексных параметров А„ (п = 1,..., Лг; ДеА„ > 0), которые есть координаты полюсов одевающей функции (см. [3] или основной текст диссертации), а также N произвольных комплексных

констант Сп. Обозначим Гп = Фо(х,-Л*) и определим N векторов с помощью следующего соотношения:

(5)

В результате, новое Л^-солитонное решение НУШ записывается в виде отношения детерминантов:

д,

(р = Ро + 2-

М

(6)

Где М есть определитель матрицы

Мпт =

А„ + Л*

а определитель М„р записывается так:

М„я =

О

1*1, а

91,/3

мг

(7)

Затем рассматривается случай, когда (ро = А - конденсат. Тогда в плоскости спектрального параметра А существует разрез от —А до А. Для упрощения вычислений и дальнейшего анализа в работе предложена, т.н. униформи-зация плоскости спектрального параметра, состоящая в применении преобразования Жуковского:

(8)

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

= ЛГ1схр(и*п), Сп = схр(10п + цп).

(9)

Где - координата полюса п-го солитонного решения, а С„ - произвольные комплексные константы. Для удобства записи окончательных формул солитонных решений мы также будем пользоваться дополнительной параметризацией Лп = ехр(гп).

Для случая конденсата явно вычисляются векторные функции q„ и соответствующие солитонные решения, выводятся общие соотношения, в частности групповые скорости.

В заключение вводится понятие регулярного решения, играющего ключевую роль в дальнейшем. Регулярным называется решение фазы которого при х —► ±оо одинаковы. Т.к. в работе мы полагаем, что модуляционная неустойчивость развивается из локализованных возмущений, на роль кандидатов для её описания мы рассматриваем только регулярные решения. Критерий того, что ЛГ-солитонное решение является регулярным, особенно просто выглядит в униформизованных переменных. Для этого необходимо, чтобы сумма угловых параметров солитонов была равна нулю или ±7г/2:

О! +- + £*„= 0, ±|. (10)

Глава 2. Солитонные решения нелинейного уравнения Шредингера на фоне конденсата

В этой главе с единой точки зрения изучаются солитонные решения НУТТТ па фоне конденсата. Глава начинается с рассмотрения односолитон-ного решения, которое мы записываем в следующем явном виде:

(И)

(12)

(13)

В зависимости от выбора положения полюса различают несколько принципиально различных типов решений. Когда полюс расположен на действительной оси вне разреза Ля > А, решение представляет собой солитон

Где

<p = -/l(cosh 2: cosh 2« + cos a cos2i>)~1 x (cosh 2 cos 2a cosh 2u + cosh 2z cos a cos 2v +i (cosh .z sin 2a sinh 2u + sinh2zcosasin2u)).

и = sex - jt + ц/2, a; = A sinh 2 cos a,

к = A cosh z sin a, Асимптотики (11) следующие

tp -Aexp(±2ia),

M2 = A2,

v = kx-Ljt- 6/2 7

A2

— cosh 2г sin 2a, 2

A2

cu = —sinh 2z cos 2a. 2

x -» ±oo, x -*■ ±oo.

Кузнецова [22|. Это локализованное решение, которое стоит на месте и осциллирует. Период осцилляций растет при Ад -* А и становиться бесконечным при Ад = А. В этом пределе солитон Кузнецова переходит в решение Перегрина [27].

Если полюс расположен на разрезе 0 < Ад < А, то решение периодично в пространстве и локализовано во времени - это т.н. бризер Ахмедиева [30]. Это решение почти гомоклинично. Оно начинается с конденсата и возвращается также к конденсату, но уже с другой фазой.

В общем случае полюс может быть расположен в любой точке правой полуплоскости. При этом решение движется с определённой групповой скоростью и осциллирует. При х —>- ±оо оно переходит в конденсат обладающий одинаковыми амплитудами, но различными фазами. В явном виде такое решение было впервые представлено в работе А.Р. Итса, A.B. Рыбина и М.А. Салля [35].

Все эти решения подробно обсуждаются в диссертации, приводятся соответствующие явные формулы. Особую важность для дальнейшего рассмотрения играет случай, когда полюс расположен близко к разрезу. Это решение представляет собой бризер большого размера L к ^.c'osa; Движущийся с большой групповой скоростью Vgroup и - и малой фазовой скоростью Vvh « Л sin о2" • называем такое решение "квази-Ахмедиевский"бризер, т.к. в пределе Д -> 1 он переходит в бризер Ахмедиева. Пример данного решения представлен на Рис. 1. Далее подробно

jUa-

Рис. 1: Квадрат модуля "квази-Ахмедиевского"бризера </> при ( = 0 с различными значениями параметра а. Слева: Я = 1.02, а = тг/3,[1 = 0, 9 = 0, справа: Д = 1.02, а = тг/11, ц = 0, в = 0.

рассматривается двухсолитонное решение. Особое внимание уделено регу-

лярным двухсолитонным решениям, которые не меняют фазу конденсата. Изучаются солитонные атомы, всевозможные нелинейные комбинации двух различных типов солитонов. Например, такие комбинации как соли-тон Кузнецова плюс бризер Ахмедиева.

Отдельная часть посвящена построению А^-солитонное решения НУШ на фоне солитона Перегрина. Этот результат позволяет избежать вычисления предела всякий раз при исследовании новой нелинейной комбинации солитонов. В дальнейшем это позволит исследовать взаимодействие соли-тона Перегрина и суперрегулярных решений.

Глава оканчивается обсуждением того, что в пределе R -*■ 1 пары квази-Ахмедиевских бризеров с противоположными значениями параметров а полностью аннигилируют, т.е. соответствующее решение есть невозмущённый конденсат.

Глава 3. Суперрегулярные решения

В этой главе содержится основной результат работы. Факт полной аннигиляции двух бризеров Ахмедиева с противоположными значениями параметра а, делает возможным построение специального и очень важного класса 2Лг-солитонных решений, описывающих эволюцию сколь угодно малых возмущений конденсата. Эти решения представляют собой подкласс регулярных решений. Мы называем их суперрегулярными солитонными решениями.

Вначале подробно описывается случай N = 1, т.е. двухсолитонное суперрегулярное решение. Для того, чтобы его построить необходимо выбрать Ri и Дг близкими к единице:

Ri = 1 + е, Д2 - 1 + ае, (14)

е > 0, а > 0, |е| « 1. Такие два полюса мы называем суперрегулярной парой. Тогда, в определённый момент времени (для определённости, t = 0), решение представляет собой слабовозмущённый конденсат:

1p = A + 6ip. (15)

Где 5ip - малое локализованное возмущение конденсата, амплитуда которого пропорциональна е. В работе подробно изучен возможный вид функции

¿¡р(х). В простейшем случае, когда полюса выбираются симметричным образом (а = 1) и константы = = 0, в\ + 02 = тг, функция ¡¡р в первом порядке такова:

5(р & 41еА

созасоз(2кх - (в1 - в2)/2) соэЬ( 2е А сое ах)

(16)

Откуда видно, что амплитуда возмущения прямо пропорциональна, а характерная ширина обратно пропорциональна величине параметра е. Вна-

Рис. 2: Суперрегулярные двухсолитонные решения с параметрами: Д1 = Дг = 1-2, сц = 7г/3, <*2 = —тг/З, /Их = дг = 0, 6\ = 62 = 7г/2 (верхние рисунки) и Д1 = 1.1, Лг = 1.05, «1 = 7г/3, «2 = —7г/3, = /л2 = 0, 61 = #2 = т/2 (нижние рисунки). Квадрат модуля решения (/5 в моменты времени < = 0 (рисунки слева) и í = 10 (рисунки справа).

чале такое возмущение растет экспоненциально, а затем формируется пара квази-Ахмедиевских бризеров, распространяющихся в противоположных направлениях. В симметричном случае а = 1 эти бризеры одинаковы и решение симметрично. В общем случае а Ф 1 решение несимметрично и описывается с помощью пяти параметров а, е, а, 81,62, а также может быть сдвинуто по х и t.

1.15

. 1.05

1.10

— 1.00

ч 0.9

0.8

0.95

-10

0

10

20

-60 - 40 -20 0 20 40 60

Рис. 3: Возмущения конденсата, соответствующие суперрегулярным решениям в начальный момент времени 1 = 0. Рисунок слева соответствует двухсолитонному решению представленному на решению Рис. 2, а рисунок справа - четырехсолитонному решению из Рис. 4.

Типичные примеры развития двухсолитонных суперрегулярных решений в симметричном и несимметричном случаях представлены на Рис. 2, а пример начального возмущения на Рис. 3 слева.

Чрезвычайно важным является тот факт, что данные рассуждения можно обобщить на любое количество пар полюсов. Мы доказали, что в первом порядке начальное возмущение соответствующее N суперрегулярным парам полюсов есть линейная сумма функций 5<рт, соответствующих отдельным двухсолитонным решениям со своими параметрами:

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

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

Гравитационные волны на глубокой воде можно описывать с помощью

N

4>п = А+ 5(рп.

(17)

771=1

-300 -200 -100

100 200 300 -300 -200 -100 0 100 200 300

Ц

-300 -200 -100

100 200 300 -300 -200 -100

too 200 300

Рис. 4: Развитие четырехсолитонного суперрегулярного решения Квадрат модуля ip в разные моменты времени. Параметры: Ri = 1.05, R3 = 1.05, аг = тг/5, а3 = -тг/5, щ = Цз = 0, дг = в3 = тг/2; Д2 = 1-075, Д4 = 1.075, ct2 = 7г/5, а4 = —тг/5, ß2 = ßi = 0, 02 = 64 = 7г/2.

НУШ впервые полученного В.Е. Захаровым [1] в 1968 г:

(да а \ ш0 д2а шо^п, ,о „ , ,

Где шо = \/дко - закон дисперсии для гравитационных волн в линейной теории. сдг = 2%\к=к0 ~ f^ " групповая скорость. Величина, которую мы измеряем в эксперименте - это высота поверхности жидкости r/(x,t), которая даётся выражением (см. [1]):

r/(x,t) = Äe{a(a;,i)exp[i(fco:E-Ci)o£)]} (19)

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

момент времени, т.е. тем больше решение годится на роль в описании развития неустойчивости из малого возмущения. В качестве компромисса нами были выбраны параметры а = 0.4 и Л = 1.15. Соответствующее решение представлено на Рис. 5. А. Чабчоуб в сотрудничестве с нами начал работу

У

1 I

Кг

Рис. 5: Двухсолитонное суперрегуляное решение с параметрами, подобранными для эксперимента: Й1 = 1.15, Дг = 1.15, а\ = 0.4, аг = -0.4, = /¿2 = О, в1 = 02 = тг/2;

по воспроизведению суперрегулярных решений в бассейне. Эксперимент показал отличное согласие с нашими теоретическими предсказаниями. На Рис. 6 представлен один из экспериментальных графиков на котором приведено сравнение теории и эксперимента для конечной стадии роста суперрегулярных возмущений.

О 10 20 30 40 50 60

Кв)

Рис. 6: Сравнение теории (красным) и эксперимента (синим) для конечной стадии роста суперрегулярных возмущений, (датчик, расположен на расстоянии 31.5 метра от генератора волн). Высота поверхности жидкости.

Глава 4. Векторное нелинейное уравнение Шредингера на фоне конденсата

В этой главе проведено обобщение результатов работы на векторное нелинейное уравнение Шредингера на фоне конденсата (ВНУШ) - систему Манакова [36]:

iVit " ~ (|</=>i[2 " Ы2 + ¡V>212 - М2)«/^ = О,

i<¿>21 ~ \чъ*х ~ (\Ч>1? ~ ki|2 + lv>2¡2 ~ Ы2)¥>2 = 0. (20)

Рассматривается задача с неубывающими граничными условиями \ч>\,2\ -»■ lai^l при -»• оо. Где ai и а2 - произвольные константы, которые мы полагаем действительными. Схема одевания аналогична той, что мы использовали для скалярного случая. Однако теперь решения ещё более многообразны. Униформизация вновь существенно упрощает вычисления и в результате мы находим общее односолитонное решение ВНУШ на фоне конденсата в следующем виде:

-a2N\ + a\N2

ifi = ai - e cos a cosh z-----,

exp(2«o + z) + 4 cosh 2и cosh z + 4 cos 2v cos a

o u aiNi+a2N2

V2 = o,2 - 8 cos a cosh z-----. (21)

exp(2uo + z) + 4 cosh 2u cosh z + 4cos2i>cosa

Где

JVi = exp(uo + z/2 + i^o + ia/2)[cos(u - a/2) cosh(u - z/2) -

isin(?; - a/2) sinh(u - z/2)], N2 = cos a cosh 2u + i sin a sinh 2и + cos 2v cosh z + i sin 2v sinh z, (22)

u0 = as0x - 70i + fJ.o/2, v0 = k0x - ui0t - d0/2,

a2

aeo = a cosh z cos a, 70 = sinh 2z sin 2a, a2

kon = a sinh z sin a, u>on = — cosh 2z cos 2a,

u = eex -7t + n/2, v = kx-ut- 9/2, a2

as = a sinh z cos a, 7 =--cosh 2z sin 2a,

2

a2

k = a cosh z sin a, ui = — sinh 2z cos 2a. (23)

В добавок к основным параметрам а и Л, решение(21) имеет дополнительные четыре параметра /¿о, во, в. Они соответствуют сдвигам по времени и пространству, а также относительным сдвигам компонент решения.

Мы проводим подробное исследование решения (21). В общем случае это солитон, движущийся из правого полупространства х > 0 и в точке ж = 0,4 = 0 разделяется на два солитона движущихся в противоположных направлениях, как показано на Рис. 7. В частных случаях мы исследуем векторные аналоги солитона Кузнецова и бризера Ахмедиева, которые в векторном случае обладают рядом интересных особенностей. Общее одно-

Рис. 7: Векторный аналог общего односолитонного решения. Квадрат модуля решения в зависимости от х и Параметры: ах = 0.5, аг = 1, II = 1, а = тг/4, цо = 0, /х = 0, в0 = 0, в = 0.

солитонное решение ВНУШ на фоне конденсата недавно было получено в работе [37]. Однако мы считаем, что использование униформизации и запись решения в виде (21) является наиболее удобной, что позволило нам качественно понять поведение решения.

Далее записывается общее двухсолитонное решение. Полное изучение его поведения выходит за рамки диссертации. Мы останавливаемся лишь на вопросе о существовании векторных аналогов суперрегулярных решений.

Как было установлено ранее, математическое объяснение возникновения двухсолитонных решений, которые являются малыми возмущениями в определённый момент времени, состоит в полной аннигиляции квази-Ахмедиевских бризеров. При этом числитель двухсолитонного решения тождественно обращается в ноль при £1 = е1а, £2 = е-"*, тогда как знаменатель остаётся конечным.

=0 5

кя1 :

Мы показали, что в векторном случае числитель тождественно равен нулю только если две из трёх независимых комплексных констант решения выбраны равными нулю. Однако это соответствует тривиальному обобщению решений скалярного НУШ на векторный случай, которое может быть получено с помощью следующего известного преобразования:

ipi(x,t) = aiip(a2t,ax), ip2(x,t) = a2'fi(a2t,ax). (24)

(где ip - решение скалярного НУШ). Таким образом, в векторном случае малое возмущение при выборе полюсов £1 = е1а, £2 = е~1а появляется только в указанном случае и не приводит к качественно новым эффектам. Вопрос о нелинейной стадии модуляционной неустойчивости в рамках ВГТТУТТТ оказался существенно более запутанным и безусловно интересен для дальнейших исследований. На данный момент не ясно, существует ли иной со литонный механизм развития малых локализованных возмущений. Интересно рассмотреть обобщение на В НУШ с тремя компонентами, которое также находит физические применения [38, 39].

Заключение

В Заключении сформулированы основные результаты диссертации.

Работы автора по теме диссертации

[1] Zakharov, V. Е. Nonlinear stage of Modulation Instability / V. E. Zakharov and A. A. Gelash // Physical Review Letters -2013 - Vol. 111, № 5. P. 054101-1 - 054101-5.

[2] Zakharov, V. Freak waves as a resuit of modulation instability / Vladimir Zakharov and Andrey Gelash // Procedia IUTAM -2013.- Vol. 9C, № 5. P. 165-175.

[3] Gelash, A. A. Superregular solitonic solutions: a novel scénario for the nonlinear stage of modulation instability / A. A. Gelash and V. E. Zakharov // Nonlinearity -2014,- Vol. 27, № 4. P. R1-R39.

Цитируемая литература

[1] Zakharov, V. E. Stability of periodic waves of finite amplitude on the surface of a deep fluid / V. E. Zakharov //J. Appl. Mech. and Tech. Phys. -1968.- Vol. 9, № 2. - P. 190-194.

[2] Akhmediev, N. Dissipative Solitons / N. Akhmediev and A. Ankiewicz // Springer-Verlag, Berlin, Heidelberg, 2005.

[3] Sulem, C. The nonlinear schrodinger equation. Self focusing and wave collapse / C. Sulem and P-L. Sulem // New York, Berlin, Heidelberg: Springer-Verlag, 1999.

[4] Kivshar, Y. S. Optical Solitons: From Fibers to Photonic Crystals / Y. S. Kivshar and G. Agrawal // Academic Press, 2003.

[5] Falkovich, G. Fluid Mechanics. A Short Course for Physicists / Gregory Falkovich // Cambridge University Press, 2011.

[6] Zakharov, V. E. Freak waves as nonlieral stage of Stokes wave modulation instability / V. E. Zakharov, A. I. Dyachenko and A. O. Prokofiev // Eur. J. Mech. B/Fluids. -2006,- Vol. 25, № 5. - P. 677-692.

[7] Zakharov, V. E. New method for numerical simulation of a nonstationary potential flow of incompressible fluid with a free surface / V. E. Zakharov, A. I. Dyachenko and O. A. Vasilyev // Eur. J. Mech. B/Fluids. -2006.- Vol. 23, № 1. P. 283-291.

[8] Pelinovsky, E. Extreme Ocean Waves / E. Pelinovsky and C. Harif // Springer, 2008.

[9] Dyachenko, A. I. On the formation of freak waves on the surface of deep wather / A. I. Dyachenko and V. E. Zakharov // JETP Lett. -2008.- Vol. 88, № 5. P. 307-311.

[10] Dysthe, K. B. Note on breather type solutions of the NLS as a model for freak-waves / K. B. Dysthe and K. Trulsen // Physica Scripta -1999,-Vol. T82. P. 48-52.

[11] Zakharov, V. E. About shape of giant breather / V. E. Zakharov and A. I. Dyachenko // Eur. J. Mech. B/Fluids. -2010- Vol. 29 P. 127-131.

[12] Zakharov, V. E. Optical solitons and quasisolitons / V. E. Zakharov and E. A. Kuznetsov // JETP. -1998.- Vol. 86, № 5. P. 1035-1046.

[13] Balakin, A. A. Structural features of the self-action dynamics of ultrashort electromagnetic pulses / A. A. Balakin, A. G. Litvak, V. A. Mironov and S. A. Skobelev // JETP. -2007.- Vol. 10, № 3. P. 363-378.

[14] Chalikov, D. Modeling of extreme waves based on equations of potential flow with a free surface / D. Chalikov and D. Sheinin //J. Comp. Phys. -2005,- Vol. 210. P. 247-273.

[15] Kibler, B. Observation of Kuznetsov-Ma soliton dynamics in optical fibre / B. Kibler, J. Fatome, C. Einot, G. Millot, G. Genty, B. Wetzel, N. Akhmediev, F. Dias and J. M. Dudley // Nature Sci. Reports -2012.-Vol. 2. P. 463.

[16] Kibler, B. The Peregrine soliton in nonlinear fibre optics / B. Kibler, J. Fatome, C. Finot, G. Millot, F. Dias, G. Genty, N. Akhmediev, and J. M. Dudley // Nature Phys. Lett. -2010,- Vol. 6. P. 790-795.

[17] Chabchoub, A. Rogue wave observations in wather wave tank / A. Chabchoub, N. P. Hoffman and N. Akhmediev // Phys. Rev. Lett. -2011.-Vol. 106, 204502.

[18] Zakharov, V. E. Exact Theory of Two-dimensional Self-focusing and One-dimensional Self-modulation of Waves in Nonlinear Media / V. E. Zakharov and A. B. Shabat // Sov. Phys.-JETP. -1972,- Vol. 34, № 1. P. 62-69.

[19] Faddeev, L. D. Hamiltonian Methods in the Theory of Solitons / L. D. Faddeev and L. A. Takhtajan // Berlin: Springer-Verlag, Berlin, Heidelberg, 2007.

[20| Kharif, C. Rogue Waves in the Ocean / C. Kharif, E. Pelinovsky and A. Slunyaev // Springer, 2009.

[21] Zakharov, V. E. Theory of solitons. The inverse scattering method / V. E. Zakharov, S. V. Manakov, S. P. Novikov and L. P. Pitaevskii // New York, London: Plenum Press, 1984.

[22] Kuznetsov, E. A. Solitons in a parametrically unstable plasma / E. A. Kuznetsov // Sov. Phys.-Dokl. (Engl. Transl.) -1977.- Vol. 22. P. 507-508.

[23] Ma, Y-C. The Perturbed Plane-Wave Solutions of the Cubic Schrodinger Equations / Y-C. Ma // Stud. Appl. Math. -1979,-Vol. 60. P. 43-58.

[24] Chiao, R. Y. Self-Trapping of Optical Beams / R. Y. Chiao, E. Garmire and C. H. Townes // Phys. Rev. Lett. -1964,- Vol. 13. P. 479-482.

[25] Agafontsev, D. Rogue waves statistics in the framework of one-dimensional Generalized Nonlinear Schrodinger Equation / D. Agafontsev and V. Zakharov // arXiv:1202.5763. -2012.

[26] Kawata, T. Inverse Scattering Method for the Nonlinear Evolution Equations under Nonvanishing Conditions / T. Kawata and H. Inoue // Journ. Phys. Soc. Japan. -1978.- Vol. 44. P. 1722-1729.

[27] Peregrine, D. H. Water waves, nonlinear Schrodinger equations and their solutions / D. H. Peregrine // J. Aust. Math. Soc. Ser. B -1983-Vol. 25. P. 16-43.

[28] Shrira, V. I. What makes the Peregrine soliton so special as a prototype of freak waves? / V. I. Shrira and V. V. Geogjaev // J. Eng. Math. -2010.-Vol. 67. P. 11-22.

[29] Akhmediev, N. How to excite a rogue wave / N. Akhmediev, J. M. Soto-Crespo and A. Ankiewicz // Phys. Lett. A -2009.- Vol. 80, 043818.

[30] Akhmediev, N. N. Generation of periodic trains of picosecond pulses in an optical fiber: exact solutions / N. N. Akhmediev,V. M. Eleonskii and N. E. Kulagin // Sov. Phys.-JETP. -1985.- Vol. 89. P. 894-899.

[31] Dubard, P. Multi-rogue waves solutions of the focusing NLS equation and the KP-I equation / P. Dubard, P. Gaillard, C. Klein and V. .B. Matveev // Eur. Phys. J. Special Topics -2010.- Vol. 185. P. 247-258.

[32] Dubard, P. On multi-rogue waves solutions of the NLS equation and positon solutions of the KdV equation / P. Dubard and V. B. Matveev // Nat. Hazards. Earth. Syst. Sci. -2011.- Vol. 11. P. 667-672.

[33] Dubard, P. Multi-rogue waves solutions: from the NLS to the KP-I equation / P. Dubard and V. B. Matveev // Nonlinearity -2013.- Vol. 26. P. R93-R125.

[34] Zakharov, V. E. Integration of nonlinear equations of mathematical physics by the method of inverse scattering II //V. E. Zakharov and A. B. Shabat, Functional Analysis and Its Applications -1980.- Vol. 13, № 3. P. 166-174.

[35] Its, A. R.. On exact integration of nonlinear schrodinger equation / A. R. Its, A. V. Rybin and M. A. Sail // Theor. Math. Phys. -1988.- Vol. 74, № 1. P. 20-32.

[36] Manakov, S. V. On the theory of two-dimensional stationary self-focusing of electromagnetic waves / S. V. Manakov // Sov. Phys.-JETP. -1974,- Vol. 38, № 2. P. 248-253.

[37] Zhao, Li-Chen Localized nonlinear waves in a two-mode nonlinear fiber / Li-Chen Zhao and Jie Liu // J. Opt. Soc. Am. B -2012,- Vol. 29, № 11. P. 3119.

[38] Yomba, E. Modulational instability and exact solutions for a three-component system of vector nonlinear Schrodinger equations / Emmanuel Yomba and George R. Sell // J. Math. Phys. -2009.- Vol. 50, 053518.

[39| Zhao, Li-Chen Rogue-wave solutions of a three-component coupled nonlinear Schrodinger equation / Li-Chen Zhao and Jie Liu // Phys. Rev. E -2013.- Vol. 87. P. 013201.