Динамика солитонов Давыдова-Скотта при пространственно-скоростном смещении компонент

ИНФОРМАТИКА И УПРАВЛЕНИЕ В ТЕХНИЧЕСКИХ И СОЦИАЛЬНЫХ СИСТЕМАХ

УДК 537.86

Л.Г. Бляхман, В.П. Морозов, В.В. Тютин

ДИНАМИКА СОЛИТОНОВ ДАВЫДОВА-СКОТТА ПРИ ПРОСТРАНСТВЕННО-СКОРОСТНОМ СМЕЩЕНИИ КОМПОНЕНТ

Национальный исследовательский университет «Высшая школа экономики»

Цель работы: Исследована динамика двухкомпонентных (векторных) солитонов Давыдова-Скотта (ДС) при пространственно-скоростном рассогласовании высокочастотных (ВЧ) и низкочастотных (НЧ) компонент. Рассмотрение проведено в рамках Захаровского типа системы двух связанных уравнений для ВЧ и НЧ поля. В этой системе ВЧ поле описывается линейным уравнением Шредингера с переменным во времени и пространстве потенциалом, вызванным НЧ компонентой. НЧ компонента в этой системе описывается уравнением Бус-синеска/Кортевега-де-Вриза со слагаемым квадратичного (стрикционного) воздействия ВЧ поля. Результат: Аналитически с применением метода интегральных моментов волнового поля найдена частота пространственно-скоростных колебаний ВЧ и НЧ компонент солитона ДС. Аналитические результаты подтверждаются численным счетом.

Научный подход: Исследование проведено как численно, так и аналитически.

Новизна: Результаты исследования новые и могут иметь приложение для исследования и описания взаимодействия нелинейных НЧ и ВЧ импульсов.

Ключевые слова: линейное уравнение Шредингера, уравнение Буссинеска, уравнение Кортевега-де-Вриза, солитон, дисперсия, нелинейность, численное моделирование.

Интерес к солитонам (устойчивым нелинейным локализованным волнам) обусловлен их способностью распространяться на значительные расстояния с сохранением своей формы. Солитонные решения возникают во многих моделях при описании распространения волн в нелинейных средах с дисперсией: при распространении поверхностных волн на глубокой и мелкой воде, внутренних волн в стратифицированной жидкости, электромагнитных, ленгмюровских и ионно-звуковых волн в плазме, оптических импульсов и пучков в волоконных линиях связи [1-8].

Динамика интенсивных высокочастотных (ВЧ) волн, например поверхностных волн на глубокой воде, электромагнитных и ленгмюровских волн в плазме, описывается вторым приближением теории дисперсии нелинейных волн, базовым уравнением которого является нелинейное уравнение Шредингера [9, 10]. ВЧ солитон огибающей в этом уравнении существует в результате баланса эффектов линейной дисперсии второго порядка и кубичной нелинейности.

Динамика интенсивных низкочастотных (НЧ) волн, например внутренних волн в стратифицированной жидкости, поверхностных волн на мелкой воде и ионно-звуковых волн в плазме, описывается двунаправленным уравнением Буссинеска (Бк) или (как частный случай) однонаправленным уравнением Кортевега-де-Вриза (КдВ). Солитонное решение в этих уравнениях возникает в результате баланса собственной нелинейности и дисперсии высокого порядка НЧ волн.

© Бляхман Л.Г., Морозов В.П., Тютин В.В., 2016.

В ряде случаев распространение ВЧ волн происходит в присутствии интенсивных НЧ волн. В этом случае взаимодействие ВЧ волн с интенсивными НЧ волнами исследуется в рамках Захаровского типа системы, состоящей из линейного уравнения Шредингера и уравнения Бк/КдВ. ВЧ поле в этой системе описывается линейным уравнением Шредингера с переменным во времени и пространстве потенциалом, обусловленным НЧ волной [11-13]. Для электромагнитных и ленгмюровских волн в плазме роль потенциала играет изменение концентрации плазмы, вызванное присутствием ионно-звуковой волны; для поверхностных волн роль потенциала играет течение, создаваемое внутренней волной на поверхности стратифицированном жидкости. Интенсивная НЧ компонента в этой системе описывается двунаправленным уравнением Бк [14, 15] или однонаправленным уравнением КдВ [16], учитывающим нелинейность и дисперсию НЧ волн и квадратичное (керровское) слагаемое воздействия ВЧ поля на НЧ волну. Данная система исследована довольно подробно в следующих случаях: 1) при соотношении коэффициентов нелинейности и дисперсии НЧ волн (Р и у, соответственно), равном Р / у = 1, найдено двухкомпонентное солитонное решение Давыдова - Скотта (ДС) [17], проведено исследование взаимодействия подобных солитонов; 2) при Р / уф 1 найдены асимптотические двухкомпонентные солитонные решения в приближении достаточно малой амплитуды ВЧ поля по сравнению с НЧ полем [13]; 3) в пренебрежении стрик-ционного воздействия ВЧ на НЧ волны, когда НЧ поле, описываемое уравнением Бк/КдВ, представляет собой солитон «разряжения», найдены локализованные стационарные состояния линейного ВЧ поля [16].

В данной работе в рамках системы Захаровского типа, состоящей из линейного уравнения Шредингера и уравнения Бк/КдВ, исследуются колебания ВЧ и НЧ компонент соли-тонов ДС при их пространственно-скоростном рассогласовании. Показан периодический характер изменения параметров ВЧ и НЧ компонент. Найдена частота подобных колебаний.

Модельная система уравнений Шредингера и Буссинеска

Рассмотрим динамику медленной огибающей и (х, X) комплексного ВЧ поля и(х, X) ехр(/к0х - гщХ) и НЧ поля п(х, X), представляющего для ВЧ компоненты эффективный потенциал (коэффициент преломления). Если ВЧ и НЧ волны представляют поверхностные и внутренние волны в океане, соответствующая система Захаровского типа представляет собой систему линейного уравнения Шредингера для поверхностных волн и уравнение Бк для внутренних волн, связанных стрикционным квадратичным слагаемым [18, 19]:

„/аи „ диЛ д2и 0 тт п

Ч"а +2пи=0-

^-г2 л+браЙ л =8аМ

а2 ^ ах2 р ах2 у дх4 ах2 '

где Унр = (дюн^ / дк)к - групповая скорость ВЧ волн; У1Ж = / к№ - фазовая скорость НЧ

волн; Р - коэффициент нелинейности НЧ волн; у - коэффициент дисперсии высокого порядка НЧ волн; е - коэффициент стрикционного воздействия ВЧ поля на НЧ волны. В очевидном случае при отсутствии ВЧ компоненты (и = 0) система (1) сводится к уравнению Бк, имеющему однопараметрическое (А - свободный параметр) решение в виде солитона

у 2 Г х - УХ Л 2 Тг2 у п =--— Бееп2 1 -т^2 1

ПД2 д ,, У = У2Р +72.

рА2 ^ А ) А2

При ип Ф 0 солитонные решения системы (1) анализировались в следующих случаях:

1) при у /Р =1 найдено точное двухкомпонентное (векторное) решение солитона ДС с

двумя свободными параметрами А, А [17]:

и = АБееЬ[л/Я(х - Уг)]ехр -^(я - ^ - V2)г + /(гнР - V

п = -ЯБееЬ2 [/Я

2) при у / Р ^ 1 в приближении достаточно слабого ВЧ поля |и| << п найдены асимптотические решения многогорбых солитонов [13];

3) при е = 0, когда НЧ компонента описывается уравнением Бк, а ВЧ поле описывается линейным уравнением Шредингера с переменным во времени и пространстве потенциалом, обусловленным присутствием НЧ поля. При НЧ волне в виде волны «разряжения» (п < 0) найдены стационарные локализованные состояния линейного ВЧ поля как решение стационарного линейного уравнения Шредингера с неоднородным потенциалом [16].

Колебания солитонов в системе уравнений Шредингера и Буссинеска

Система (1) при нулевых условиях на бесконечности (п,1) ^ ^ 0 имеет следующие соотношения для интегральных моментов волнового поля:

пт 1 +да ^2 +да

%  = 0, £ |пЛ =

7 +да +да

-да -да

где и = |и| ехр (/ф), К = дф/ дх - волновое число. В адиабатическом приближении решение

системы (4)-(6) представим в виде солитоноподобных (БесЬ-Нке) импульсов с пространственно-скоростным смещением ВЧ компоненты относительно НЧ компоненты при их фиксированных амплитудах и протяженностях:

и = А8ееь[л/Цх - У - х(г))]ехр - 1- (я - УНР - V2) + / (унр - V)х + /к(г)х

п = -ЯБееЬ2 [л/Я

+да

где х (г) = и1 j х||| -х - центр «масс» ВЧ компоненты; к (г) = К + V,

-да

V2 = V22p + 4Я-еА2 / Я. В начальный момент времени соотношение (7) отвечает солитону (3) со смещением ВЧ компоненты. Соотношение (7) содержит два свободных параметра: центр «масс» х(г) и волновое число к (г). Подставляя (7) в (5) и (6), обозначая ^ = х - VI и считая

малой вариацию центра «масс» пакета И << 1/>/Я, получим

= 0.

—2 15

Из (8) следует, что ВЧ компонента осциллирует около невозмущенной НЧ компоненты с частотой

ю = 2Ял/2/15 .

-да

Модельная система уравнений Шредингера и Кортевега-де Вриза

Для ряда сред величины дисперсии и нелинейности НЧ волн и коэффициента стрик-ционного воздействия ВЧ волн на НЧ компоненту в уравнении Бк достаточно малы. В соответствии с предположением о распространении однонаправленной волны уравнение Бк может быть сведено к уравнению ^В типа. В этом случае наиболее эффективное взаимодействие интенсивных ВЧ и НЧ волн возникает при групповом синхронизме ВЧ и НЧ волн, когда групповая скорость ВЧ волн УНР = (дюНР / дксовпадает с фазовой скоростью НЧ волн

У1Р = /к1Р , распространяющейся в том же направлении: УНР = У1Р = Уш. Подобная ситуация возникает, в частности, при взаимодействии поверхностных и внутренних волн в стратифицированном океане [16]. Оценим длины поверхностных и внутренних волн (Л^. и А1Р соответственно) при групповом синхронизме. При величине частоты Брента-Вяйсяля юву ~ 0.01 Гц и длине внутренней волны порядка Л№ ~ 100 м имеем для скорости линейных внутренних волн У^ ~ 15 см/с [18]. Характерное значение частоты поверхностных волн на глубокой воде при глубине в несколько сотен метров [19] составляет порядка ~1 Гц, что на два порядка больше частоты внутренних волн и позволяет провести разделение описания взаимодействия поверхностных и внутренних волн на ВЧ и НЧ волны соответственно. В классическом дисперсионном соотношении поверхностных волн на глубокой воде юнР = (g - ускорение свободного падения) групповой синхронизм при выбранных параметрах осуществляется при длине поверхностных волн порядка Лнр ~ 2 м, что на два порядка меньше длины волны внутренних волн. В этом случае система (1), при выполнении условия группового синхронизма в системе отсчета, движущейся с групповой скоростью У$ш, ^ = х — Х = X, в пренебрежении второй производной по времени, что возможно при достаточно малых величинах нелинейности и дисперсии внутренних волн, примет вид [ 16]

_. ди д2и _ тт Л .„ „ .

2/—---— + 2пи = 0,

&

дп д(п2) д3п д(и|2) дХ д^ д^3

2Узм — — 6Р -П2 + У — = .

При у / Р = 1 система (10) имеет двухкомпонентное (векторное) солитонное решение ДС с двумя свободными параметрами Ш, А [17]:

и = А8ееЬ[л/Ц^ — У )]ехр

где У = (4уШ — 8А2/ Ш)/(2У81Н). Скорость движения солитона (11) в лабораторной системе отсчета при достаточно малых величинах дисперсии и нелинейности НЧ волн и коэффициента стрикционной силы воздействия ВЧ волн на НЧ компоненту по сравнению с квадратом скорости НЧ волн близка к скорости солитона (3): У « У1Р + 2уШ/У1Р — 8А2 /(2У1Рш).

Колебания солитонов в системе уравнений Шредингера и Кортевега-де Вриза

Система (10) при нулевых граничных условиях на бесконечности (п, и^ 0 имеет следующие соотношения для динамики интегральных моментов волнового поля:

% " & /М 2&^ =0  & I <^=0 •

Л Ж * &

00

—о

—о

7 +да

-\{sK\U\2 —Fsinn2—5 = 0,

—да

- j ^ = —J n |•

-ГУ1 -(-У1 ^

— dt

-да

+да

Решение системы (12)-(15) представим в виде солитона ДС (11) с пространственно-скоростным смещением ВЧ компоненты относительно НЧ компоненты и изменением амплитуд и протяженностей ВЧ и НЧ компонент:

U = a(t )Asech

"VAfe—Vt—5(i )) "

A

exp [— ïkt — iV5 + ik (t )5], n = —b(t )Asech2

V^

A

где 5 = N 1 J5U —5 - центр «масс» ВЧ компоненты; k(t) = К + V, a2(t)A(t) = const,

—да

b(t)A(t) = const ; a(o) = b(0) = A(0) = 1. Соотношение (16) удовлетворяет интегральным соотношениям (12)-(15) и содержит три свободных параметра: координату центра «масс» 5 if) • протяженность A(t ) и волновое число k (t ). Подставляя (16) в (13), получим первый интеграл

a2 = b = - = 1 + 3VnA2

A 2A2

где k0 = k(0). Далее, подставляя (16) в (14) и (15), обозначая р = 5 — Vt и считая малым смещение центра «масс» |р| << A / VA , получим уравнение

d 2р / dt2 +

совпадающее с (8) для возмущенного солитона системы уравнений Шредингера и Бк. Частота колебаний ВЧ и НЧ компонент волнового поля из (16) составляет

ю = 2W2/15 ,

которая совпадает с частотой колебаний (9), полученной для возмущенного солитона системы уравнений Шредингера и Бк. При начальном импульсе (16) с k0 Ф 0 и р(0)= 0 (скоростное смещение ВЧ компоненты) изменение параметров солитона во времени описывается соотношениями :

2

р = — ^sinЦ) a2 = b = - = 1 —3sVsinA k° sin2 ю A A2

л

ю —t

v 2 У

При начальном импульсе (16) с к0 = 0 и р(0)^ 0 (пространственное смещение ВЧ компоненты) изменение параметров солитона во времени описывается соотношениями:

р = р

А 2Я

Численные результаты

Рассмотрим начальную задачу динамики пакета

и

в рамках системы уравнений (10) при VSIN = 1, Р = у = е = 1/10. При к0 = 0 импульс (22) совпадает с солитоном (11) с Я = 1 в начальный момент времени. При к0 Ф 0 импульс (22) опи-

да

—да

сывает скоростное рассогласование компонент солитона

при к0 = 1/3 (скоростное рассогласование компонент). Начальное скоростное рассогласование компонент солитона приводит к периодическому колебанию параметров солитона с постоянным периодом Т « 8,6, близким к величине 2ж / ю, полученной из аналитического рассмотрения (соотношение (19)).

Рис. Пространственно-временные распределения ^

Аналогичные результаты наблюдаются и при пространственном рассогласовании компонент солитона в начальный момент времени. Т.е. ВЧ компонента солитона ДС осциллирует около центра масс НЧ компоненты с постоянной частотой, близко описываемой формулой (19) в аналитическом описании.

Заключение

В данной работе исследованы колебания высокочастотной компоненты векторного солитона Давыдова - Скотта (ДС), обусловленные их начальным пространственно-скоростным рассогласованием. Рассмотрение проведено в рамках Захаровского типа системы двух связанных уравнений для ВЧ и НЧ полей. В этой системе ВЧ поле описывается линейным уравнением Шредингера с переменным во времени и пространстве потенциалом, обусловленным присутствием НЧ компоненты. Последняя в этой системе описывается или уравнением Буссинеска, или Кортевега - де-Вриза с квадратичным (стрикционным) воздействием ВЧ поля. Показано, что при малом различии в пространственном положении ВЧ и НЧ компонент или при малом различии их скоростей компоненты колеблются друг относительно друга с постоянной частотой, т.е. имеет место динамическое равновесие этих компонент векторного солитона ДС. Аналитически с использованием метода моментов найдена частота пространственно-скоростных колебаний ВЧ компоненты векторного солитона ДС около НЧ компоненты. Аналитические результаты подтверждаются численным счетом.

Данная работа выполнена при поддержке Российского фонда фундаментальных исследований (проект РФФИ № 15-02-01919 а).

Библиографический список

1. Infeld, E. Nonlinear Waves, Solitons, and Chaos / E. Infeld, G. Rowlands. - Cambridge: Cambridge University Press, 2000.

2. Agrawal, G. P. Nonlinear Fiber Optic / G. P. Agrawal. - San Diego: Academic Press, 2001.

3. Yang, J. Solitons in Field Theory and Nonlinear Analysis / J. Yang. - New York: Springer, 2001.

4. Kivshar, Y. S. Optical Solitons: From Fibers to Photonic / Y. S. Crystals Kivshar, G. P. Agrawal. -San Diego: Academic, 2003.

5. Dickey, L. A. Soliton Equations and Hamiltonian Systems /L. A. Dickey. - New York:World Scientific, 2005.

6. Malomed, B. A. Soliton Management in Periodic Systems / B. A. Malomed. - New York: Springer, 2006.

7. Dauxois, T. Physics of Solitons / T. Dauxois, M. Peyrard. - Cambridge: Cambridge University Press, 2006.

8. Zakharov, V. E. Sov. Phys. JETP / V. E. Zakharov, A. B. Shabat. 1972. V 34. № 62.

9. Hasegawa, А. Appl. Phys. Lett / А. Hasegawa, F. Tappert. 1973. V 23. № 142.

10. Kartashov, Y. V. Spontaneous Symmetry Breaking, Self-Trapping, and Josephson Oscillations / Y. V. Kartashov [et al.]; edited by B. A. Malomed. - Berlin and Heidelberg: Springer, 2013.

11. Malomed, B.A. J. Opt. Soc. Am. 2014. B 31. Р. 2460.

12. Ostrovskaya, E.A. Phys Lett / E.A. Ostrovskaya [et al.]. 2001. A 282. № 157.

13. Janssen, P. The Interaction of Ocean Waves and Wind / P. Janssen. - Cambridge: Cambridge University Press, 2009.

14. Brunetti, M. Phys. Lett / M. Brunetti[et al.]. 2014. A 378. № 1025.

15. A.S. Davydov, A.S. Solitons in Molecular Systems / A.S. Davydov. - Dordrecht: Reidel, 1987.

Дата поступления в редакцию 10.05.2016

L.G. Blyahman, V.P. Morozov, V.V. Tyutin

DYNAMICS OF THE DAVYDOV-SCOTT SOLITONS AT SPACE-SPEED DISPLACEMENT OF THE SOLITON'S COMPONENT

National Research University Higher School of Economics

Purpose: The dynamics of the two-component (vector) Davydov-Scott (DS) solitons with a space-speed displacement of high-frequency (HF) and low-frequency (LF) component. Investigation was conducted in the framework of Zaharov type system of two coupled equations for HF and LF field. In this system, the HF field is described by the linear Schrodinger equation with variable in time and space potential, caused by LF component. LF component in this system is described by Boussinesq or Korteweg-de Vries equation with a quadratic (striction) action to HF fields. Approach: The investigation was considered as analytically as numerically.

Findings: Analytically, using the method of wave field integrals found frequency of spatially-speed oscillation HF components of the vector soliton DS. The analytical results are confirmed by a numerical simulations.

Key words: linear schrodinger equation, boussinesq equation, Korteweg-de Vries equation, soliton, dispersion, nonlinearity, numerical simulation.