Уединенная волна. Ударные волны

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

1. Ударная волна. Ударная волна (скачок уплотнения) - это распространяющаяся со сверхзвуковой скоростью тонкая переходная область, в которой происходит резкое увеличение плотности, давления и скорости вещества. Она возникает при взрывах, детонации, при сверхзвуковых движениях тел, при мощных электрических разрядах и т.д. Например, при взрыве образуются продукты взрыва, обладающие большой плотностью и находящиеся под большим давлением. Расширяющиеся продукты взрыва сжимают окружающий воздух, причем в каждый момент времени сжатым оказывается лишь воздух, находящийся в определенном объеме, вне этого объема воздух остается в невозмущенном состоянии. С течением времени объем сжатого воздуха возрастает. Поверхность, которая отделяет сжатый воздух от невозмущенного воздуха, и представляет собой ударную волну (или как говорят, фронт ударной волны). На рис. 6.27,а в качестве примера приведен график распределения плотности в ударной волне, распространяющейся в реальном газе ( – плотность газа перед фронтом волны).

При ускоренном движении тела ударная волна возникает не сразу. Сначала возникает волна сжатия с непрерывными распределениями плотности и давления. С течением времени крутизна передней части волны возрастает и в некоторый момент времени происходит резкий скачок всех гидродинамических величин, возникает ударная волна.

В случае движения тела со сверхзвуковой скоростью (
) звуковые волны охватывают лишь часть объема газа, лежащую позади движущегося тела и ограниченную некоторой поверхностью, называемой характеристической поверхностью, поверхностью слабого разрыва или фронтом ударной волны.

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

В случае электромагнитных волн аналогом ударной звуковой волны, возникающей при движении тел со скоростями, превышающими фазовые скорости упругих волн в данной среде, является излучение Вавилова – Черенкова (см. §7.4.4).

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

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

Солитон (от лат. solus – один) – структурно устойчивая уединенная волна в нелинейной диспергирующей среде. Структура солитона поддерживается стационарной за счет баланса между действием нелинейности среды и дисперсии.

Солитон впервые наблюдался на водяном канале в 1834 г., когда при резкой остановке баржи около ее носа образовался водяной выступ (водяной холм) и затем он стал самостоятельно распространяться по каналу, сохраняя на протяжении длительного времени свою структуру и скорость.

Рассмотрим возможность образования солитона на поверхности воды. Для волн, у которых длина волны значительно превышает глубину
водоема (
, мелкая вода) явление дисперсии отсутствует, они распространяются с фазовой скоростью
, где – ускорение свободного падения, а - смещение поверхности жидкости в вертикальном направлении в данной точке профиля волны (см. рис. 6.27,в). Из записанной формулы для фазовой скорости следует, что вершина водяного холма движется быстрее, чем точки вблизи его подножия. Это нелинейность среды приводит к тому, что крутизна фронта волны возрастает с течением времени, т.е. происходит пространственное сужение водяного холма (см. рис. 6.28,б).

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

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

Солитоны играют важную роль в теории конденсированного состояния вещества, в частности в квантовой статистике, теории фазовых переходов. Структуры в форме солитонов обнаружены во многих динамических системах – в плазме, радиосхемах, лазерах, нервных волокнах.

Учебное издание

Марс Гильманович Валишев

Александр Александрович Повзнер

Солитоны бывают различной природы:

Математическая модель

Уравнение Кортевега - де Фриза

Одной из простейших и наиболее известных моделей, допускающих существование солитонов в решении, является уравнение Кортевега - де Фриза:

u_t - 6 u u_x + u_{xxx} = 0

Одним из возможных решений данного уравнения является уединённый солитон:

$u(x,t)\; =\; -\; \backslash frac\{2\backslash varkappa^2\}\{\; \backslash mathrm\{ch\}^2\backslash ,\backslash varkappa(x-4\backslash varkappa^2\; t-\backslash varphi)\; \}$

где $2\backslash varkappa^2$ - амплитуда солитона, $\backslash varphi$ - фаза. Эффективная ширина основания солитона равна $\backslash varkappa^\{-1\}$. Такой солитон движется со скоростью $v\; =\; 4\backslash varkappa^2$. Видно, что солитоны с большой амплитудой оказываются более узкими и движутся быстрее .

В более общем случае можно показать, что существует класс многосолитонных решений, таких что асимптотически при $t\backslash to\; \backslash pm\backslash infty$ решение распадается на несколько удалённых одиночных солитонов, движущихся с попарно различными скоростями. Общее N-солитонное решение можно записать в виде

$u(x,t)\; =\; -2\; \backslash frac\{d^2\}\{dx^2\}\; \backslash ln\; \backslash det\; A(x,t)$

где матрица $A(x,t)$ даётся выражением

$A\_\{nm\}\; =\; \backslash delta\_\{nm\}\; +\; \backslash frac\{\backslash beta\_n\}\{\backslash varkappa\_n\; +\; \backslash varkappa\_m\}\backslash mathrm\{e\}^\{8\backslash varkappa\_n^3\; t\; -(\backslash varkappa\_n\; +\; \backslash varkappa\_m)x\}$

Здесь $\backslash beta\_n,\; n=1,\backslash dots,N$ и $\backslash varkappa\_n>0,\; n=1,\backslash dots,N$ - произвольные вещественные постоянные.

Замечательным свойством многосолитонных решений является безотражательность : при исследовании соответствующего одномерного уравнения Шрёдингера

$-\backslash partial^2\_x\backslash psi(x)\; +\; u(x)\backslash psi(x)\; =\; E\backslash psi(x)$

с потенциалом $u(x)$, убывающим на бесконечности быстрее чем $|x|^\{-1-\backslash varepsilon\}$, коэффициент отражения равен 0 тогда и только тогда, когда потенциал есть некоторое многосолитонное решение уравнения КдФ в некоторый момент времени $t$.

Интерпретация солитонов как некоторых упруго взаимодействующих квазичастиц основана на следующем свойстве решений уравнения КдФ. Пусть при $t\backslash to\; -\backslash infty$ решение имеет асимптотический вид $N$ солитонов, тогда при $t\backslash to\; +\backslash infty$ оно также имеет вид $N$ солитонов с теми же самыми скоростями, но другими фазами, причём многочастичные эффекты взаимодействия полностью отсутствуют. Это означает, что полный сдвиг фазы $k$-го солитона равен

$\backslash Delta\backslash varphi\_k\; =\; \backslash sum\_\{\backslash stackrel\{n=1\}\{n\backslash ne\; k\}\}^\{N\}\; \backslash Delta\backslash varphi\_\{nk\}$

Пусть $n$-ый солитон движется быстрее, чем $m$-ый, тогда

$\backslash Delta\backslash varphi^\{+\}\_\{n\}\; =\; \backslash Delta\backslash varphi\_\{kn\}\; =\; \backslash frac\{1\}\{\backslash varkappa\_n\}\backslash ln\backslash left|\; \backslash frac\{\backslash varkappa\_n+\backslash varkappa\_m\}\{\backslash varkappa\_n-\backslash varkappa\_m\}\; \backslash right|$ $\backslash Delta\backslash varphi^\{-\}\_\{k\}\; =\; \backslash Delta\backslash varphi\_\{nk\}\; =\; -\; \backslash frac\{1\}\{\backslash varkappa\_m\}\backslash ln\backslash left|\; \backslash frac\{\backslash varkappa\_n+\backslash varkappa\_m\}\{\backslash varkappa\_n-\backslash varkappa\_m\}\; \backslash right|$

то есть фаза более быстрого солитона при парном столкновении увеличивается на величину $\backslash Delta\backslash varphi^\{+\}\_\{n\}$, а фаза более медленного - уменьшается на $\backslash Delta\backslash varphi^\{-\}\_\{k\}$, причём полный сдвиг фазы солитона после взаимодействия равен сумме сдвигов фаз от попарного взаимодействия с каждым другим солитоном.

Нелинейное уравнение Шрёдингера

при значении параметра $\backslash nu\; >\; 0$ допустимы уединённые волны в виде:

$u\; \backslash left(x,t\; \backslash right)\; =\; \backslash left(\backslash sqrt\{\backslash frac\{2\; \backslash alpha\}\{\backslash nu\}\; \}\; \backslash right)\; \backslash mathrm\{ch\}^\{-1\}\; \backslash left(\backslash sqrt\{\backslash alpha\}(x\; -\; Ut)\; \backslash right)\; e^\{i(r\; x-st)\},$

где $r,\; s,\backslash alpha,U$ - некоторые постоянные, связанные соотношениями:

$U=2r$ $s=r^2-\backslash alpha$

См. также

Напишите отзыв о статье "Солитон"

Примечания

1. J.S.Russell «Report on Waves»: (Report of the fourteenth meeting of the British Association for the Advancement of Science, York, September 1844 (London 1845), pp 311-390, Plates XLVII-LVII)
2. J.S.Russell (1838), Report of the committee on waves, Report of the 7th Meeting of British Association for the Advancement of Science, John Murray, London, pp.417-496.
3. Абловиц М., Сигур Х. Солитоны и метод обратной задачи. М.: Мир, 1987, с.12.
4. N.J.Zabusky and M.D.Kruskal (1965), Interaction of solitons in a collisionless plasma and the recurrence of initial states, Phys.Rev.Lett., 15 pp. 240-243.
5. Дж. Л. Лэм. . - М .: Мир , 1983. - 294 с.
6. А. Т. Филиппов. Многоликий солитон. - С. 40-42.
7. А. Т. Филиппов. Многоликий солитон. - С. 227-23.
8. - статья из Физической энциклопедии
9. Vladimir Belinski, Enric Verdaguer. . - Cambridge University Press , 2001. - 258 с. - (Cambridge monographs on mathematical physics). - ISBN 0521805864 .
10. Н. Н. Розанов // Природа . - 2007. - № 6 .
11. А. Т. Филиппов. Многоликий солитон. - С. 241-246.
12. А. И. Маймистов // Квантовая электроника . - 2010. - Т. 40 , № 9 . - С. 756-781 .
13. Andrei I Maimistov (англ.) // Quantum Electronics . - 2010. - Vol. 40. - P. 756. - DOI :10.1070/QE2010v040n09ABEH014396 .
15. Сазонов С. В. Оптические солитоны в средах из двухуровневых атомов // Научно-технический вестник информационных технологий, механики и оптики. 2013. Т. 5. № 87. С. 1-22.

Литература

• Абловиц М., Сигур Х. Солитоны и метод обратной задачи. - М .: Мир, 1987. - 480 с.
• Додд Р., Эйлбек Дж., Гиббон Дж., Моррис Х. Солитоны и нелинейные волновые уравнения. - М .: Мир, 1988. - 696 с.
• Захаров В. Е., Манаков С. В., Новиков С. П., Питаевский Л. П. Теория солитонов: Метод обратной задачи. - М .: Наука, 1980. - 320 с.
• Инфельд Э., Роуландс Дж. Нелинейные волны, солитоны и хаос. - М .: Физматлит, 2006. - 480 с.
• Лэм Дж. Л. Введение в теорию солитонов. - М .: Мир, 1983. - 294 с.
• Ньюэлл А. Солитоны в математике и физике. - М .: Мир, 1989. - 328 с.
• Самарский А. А., Попов Ю. П. Разностные методы решения задач газовой динамики. - М .: URSS, 2004. - 424 с.
• Уизем Дж . Линейные и нелинейные волны. - М .: Мир, 1977. - 624 с.
• Филиппов А. Т. Многоликий солитон // Библиотечка «Квант». - Изд. 2, перераб. и доп.. - М .: Наука, 1990. - 288 с.
• Yaroslav V. Kartashov, Boris A. Malomed, Lluis Torner (англ.) // Reviews of Modern Physics . - 2011. - Vol. 83. - P. 247–306.
• (англ.) // Physics . - 2013. - Vol. 6. - P. 15. - DOI :10.1103/Physics.6.15 .

Ссылки

Отрывок, характеризующий Солитон

СОЛИТОН –это уединенная волна в средах различной физической природы, сохраняющая неизменной свою форму и скорость при распространении.От англ. solitary – уединенная (solitary wave – уединенная волна), «-он» – типичное окончание терминов такого рода (например, электрон, фотон, и т.д.), означающее подобие частицы.

Понятие солитон введено в 1965 американцами Норманом Забуски и Мартином Крускалом, но честь открытия солитона приписывают британскому инженеру Джону Скотту Расселу (1808–1882). В 1834 им впервые дано описание наблюдения солитона («большой уединенной волны»). В то время Рассел изучал пропускную способность канала Юнион близь Эдинбурга (Шотландия). Вот как сам автор открытия рассказывал о нем: «Я следил за движением баржи, которую быстро тянула по узкому каналу пара лошадей, когда баржа неожиданно остановилась; но масса воды, которую баржа привела в движение, не остановилась; вместо этого она собралась около носа судна в состоянии бешенного движения, затем неожиданно оставила его позади, катясь вперед с огромной скоростью и принимая форму большого одиночного возвышения, т.е. округлого, гладкого и четко выраженного водяного холма, который продолжал свой путь вдоль канала, нисколько не меняя своей формы и не снижая скорости. Я последовал за ним верхом, и когда я нагнал его, он по-прежнему катился вперед со скоростью приблизительно восемь или девять миль в час, сохранив свой первоначальный профиль возвышения длиной около тридцати футов и высотой от фута до фута с половиной. Его высота постепенно уменьшалась, и после одной или двух миль погони я потерял его в изгибах канала. Так в августе 1834 мне впервые довелось столкнуться с необычайным и красивым явлением, которое я назвал волной трансляции…».

Впоследствии Рассел экспериментальным путем, проведя ряд опытов, нашел зависимость скорости уединенной волны от ее высоты (максимальной высоты над уровнем свободной поверхности воды в канале).

Возможно, Рассел предвидел ту роль, которую играют солитоны в современной науке. В последние годы своей жизни он завершил книгу Волны трансляции в водном, воздушном и эфирном океанах , опубликованную посмертно в 1882. Эта книга содержит перепечатку Доклада о волнах – первое описание уединенной волны, и ряд догадок о строении материи. В частности, Рассел полагал, что звук есть уединенные волны (на самом деле это не так), иначе, по его мнению, распространение звука происходило бы с искажениями. Основываясь на этой гипотезе и используя найденную им зависимость скорости уединенной волны, Рассел нашел толщину атмосферы (5 миль). Более того, сделав предположение, что свет это тоже уединенные волны (что тоже не так), Рассел нашел и протяженность вселенной (5·10 17 миль).

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

Реакция на научное сообщение Рассела наиболее авторитетных в то время английских механиков Джорджа Байделя Эйри (1801–1892) (профессора астрономии в Кембридже с 1828 по 1835, астронома королевского двора с 1835 по 1881) и Джорджа Габриэля Стокса (1819–1903) (профессора математики в Кембридже с 1849 по 1903) была отрицательной. Много лет спустя солитон был переоткрыт при совсем иных обстоятельствах. Интересно, что и воспроизвести наблюдение Рассела оказалось не просто. Участникам конференции «Солитон-82», съехавшимся в Эдинбург на конференцию, приуроченную к столетию со дня смерти Рассела и пытавшимся получить уединенную волну на том самом месте, где ее наблюдал Рассел, ничего увидеть не удалось, при всем их опыте и обширных знаниях о солитонах.

В 1871–1872 были опубликованы результаты французского ученого Жозефа Валентена Буссинеска (1842–1929), посвященных теоретическим исследованиям уединенных волн в каналах (подобных уединенной волне Рассела). Буссинеск получил уравнение:

Описывающее такие волны (u – смещение свободной поверхности воды в канале, d – глубина канала, c 0 – скорость волны, t – время, x – пространственная переменная, индекс соответствует дифференцированию по соответствующей переменной), и определил их форму (гиперболический секанс, см . рис. 1) и скорость.

Исследуемые волны Буссинеск называл вспучиваниями и рассмотрел вспучивания положительной и отрицательной высоты. Буссинеск обосновал устойчивость положительных вспучиваний тем, что их малые возмущения, возникнув, быстро затухают. В случае отрицательного вспучивания образование устойчивой формы волны невозможно, как и для длинного и положительного очень короткого вспучиваний. Несколько позже, в 1876, опубликовал результаты своих исследований англичанин лорд Рэлей.

Следующим важным этапом в развитии теории солитонов стала работа (1895) голландцев Дидерика Иоганна Кортевега (1848–1941) и его ученика Густава де Вриза (точные даты жизни не известны). По-видимому, ни Кортевег, ни де Вриз работ Буссинеска не читали. Ими было выведено уравнение для волн в достаточно широких каналах постоянного поперечного сечения, носящее ныне их имя – уравнение Кортевега-де Вриза (КдВ). Решение такого уравнения и описывает в свое время обнаруженную Расселом волну. Основные достижения этого исследования состояли в рассмотрении более простого уравнения, описывающего волны, бегущие в одном направлении, такие решения более наглядны. Из-за того, что в решение входит эллиптическая функция Якоби cn , эти решения были названы «кноидальными» волнами.

В нормальной форме уравнение КдВ для искомой функции и имеет вид:

Способность солитона сохранять при распространении свою форму неизменной объясняется тем, что поведение его определяется двумя действующими взаимно противоположно процессами. Во-первых, это, так называемое, нелинейное укручение (фронт волны достаточно большой амплитуды стремится опрокинуться на участках нарастания амплитуды, поскольку задние частицы, имеющие большую амплитуду, движутся быстрее впереди бегущих). Во-вторых, проявляется такой процесс как дисперсия (зависимость скорости волны от ее частоты, определяемая физическими и геометрическими свойствами среды; при дисперсии разные участки волны движутся с разными скоростями и волна расплывается). Таким образом, нелинейное укручение волны компенсируется ее расплыванием за счет дисперсии, что и обеспечивает сохранение формы такой волны при ее распространении.

Отсутствие вторичных волн при распространении солитона свидетельствует о том, что энергия волны не рассеивается по пространству, а сосредоточена в ограниченном пространстве (локализована). Локализация энергии есть отличительное качество частицы.

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

Долго считалось, что уединенные волны связаны только с волнами на воде и изучались они специалистами – гидродинамиками. В 1946 М.А.Лаврентьев (СССР), а в 1954 К.О.Фридрихс и Д.Г.Хайерс США опубликовали теоретические доказательства существования уединенных волн.

Современное развитие теории солитонов началось с 1955, когда была опубликована работа ученых из Лос Аламоса (США) – Энрико Ферми, Джона Пасты и Стена Улама, посвященная исследованию нелинейных дискретно нагруженных струн (такая модель использовалась для изучения теплопроводности твердых тел). Длинные волны, бегущие по таким струнам, оказались солитонами. Интересно, что методом исследования в этой работе стал численный эксперимент (расчеты на одной из первых созданных к этому времени ЭВМ).

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

нелинейное уравнение Шредингера (НУШ)

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

Уравнение sin-Гордона (СГ)

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

Солитонные решения имеют и так называемые, родственные КдВ уравнения. К таким уравнениям относятся,

модифицированное уравнение КдВ

уравнение Бенджамина, Бона и Магони (ББМ)

впервые появившееся при описании боры (волны на поверхности воды, возникающей при открывании ворот шлюзов, при «запирании» течения реки);

уравнение Бенджамина – Оно

полученное для волн внутри тонкого слоя неоднородной (стратифицированной) жидкости, расположенного внутри другой однородной жидкости. К уравнению Бенджамина – Оно приводит и исследованиее трансзвукового пограничного слоя.

К уравнениям с солитонными решениями относится и уравнение Борна – Инфельда

имеющее приложения в теории поля. Есть и другие уравнения с солитонными решениями.

Солитон, описываемый уравнением КдВ, однозначно характеризуется двумя параметрами: скоростью и положением максимума в фиксированный момент времени.

Солитон, описываемый уравнением Хироты

однозначно характеризуется четырьмя параметрами.

Начиная с 1960, на развитие теории солитонов повлиял ряд физических задач. Была предложена теория самоиндуцированной прозрачности и приведены экспериментальные результаты, ее подтверждающие.

В 1967 Крускалом и соавторами был найден метод получения точного решения уравнения КдВ – метод так называемой обратной задачи рассеяния. Суть метода обратной задачи рассеяния состоит в замене решаемого уравнения (например, уравнения КдВ) системой других, линейных уравнений, решение которых легко находится.

Этим же методом в 1971 советскими учеными В.Е.Захаровым и А.Б.Шабатом было решено НУШ.

Приложения солитонной теории в настоящее время находят применение при исследованиях линий передачи сигналов с нелинейными элементами (диоды, катушки сопротивления), пограничного слоя, атмосфер планет (Большое красное пятно Юпитера ), волн цунами, волновых процессов в плазме, в теории поля, физике твердого тела, теплофизике экстремальных состояний веществ, при изучении новых материалов (например, джозефсоновских контактов, состоящих из разделенных диэлектриком двух слоев сверхпроводящего металла), при создании моделей решеток кристаллов, в оптике, биологии и многих других. Высказано мнение, что бегущие по нервам импульсы – солитоны.

В настоящее время описаны разновидности солитонов и некоторые комбинаций из них, например:

антисолитон – солитон отрицательной амплитуды;

бризер (дублет) – пара солитон – антисолитон (рис. 2);

мультисолитон – несколько солитонов, движущихся как единое целое;

флюксон – квант магнитного потока, аналог солитона в распределенных джозефсоновских контактах;

кинк (монополь), от английского kink – перегиб.

Формально кинк можно ввести как решение уравнений КдВ, НУШ, СГ, описываемое гиперболическим тангенсом (рис. 3). Изменение знака решения типа «кинк» на противоположный дает «антикинк».

Кинки были обнаружены в 1962 англичанами Перрингом и Скирмом при численном (на ЭВМ) решении уравнения СГ. Таким образом, кинки были обнаружены раньше, чем появилось название солитон. Оказалось, что столкновение кинков не привело ни к их взаимному уничтожению, ни к последующему возникновению других волн: кинки, таким образом, проявили свойства солитонов, однако название кинк закрепилось за волнами такого рода.

Солитоны могут быть также двумерными и трехмерными. Изучение неодномерных солитонов осложнялось трудностями доказательства их устойчивости, однако в последнее время получены экспериментальные наблюдения неодномерных солитонов (например, подковообразные солитоны на пленке стекающей вязкой жидкости, изучавшиеся В.И.Петвиашвили и О.Ю.Цвелодубом). Двумерные солитонные решения имеет уравнение Кадомцева – Петвиашвили, используемое, например, для описания акустических (звуковых) волн:

Среди известных решений этого уравнения – нерасплывающиеся вихри или солитоны-вихри (вихревым является течение среды, при котором ее частицы имеют угловую скорость вращения относительно некоторой оси). Солитоны такого рода, найденные теоретически и смоделированные в лаборатории, могут самопроизвольно возникать в атмосферах планет. По своим свойствам и условиям существования солитон-вихрь подобен замечательной особенности атмосферы Юпитера – Большому Красному Пятну.

Солитоны являются существенно нелинейными образованиями и столь же фундаментальны, как линейные (слабые) волны (например, звук). Создание линейной теории, в значительной мере, трудами классиков Бернхарда Римана (1826–1866), Огюстена Коши (1789–1857), Жана Жозефа Фурье (1768–1830) позволило решить важные задачи, стоявшие перед естествознанием того времени. С помощью солитонов удается выяснить новые принципиальные вопросы при рассмотрении современных научных проблем.

Андрей Богданов

Формат: doc

Дата создания: 31.05.2003

Размер: 125.1 KB

1. Введение

Волны в природе

Из школьного курса физики хорошо извест­но, что если в какой-либо точке упругой среды (твердой, жидкой или газообразной) возбудить ко­лебания, то они будут передаваться в другие места. Эта передача возбуждений обусловлена тем, что близкие участки среды связаны друг с другом. При этом колебания, возбужденные в одном месте, рас­пространяются в пространстве с определенной ско­ростью. Волной принято называть процесс передачи возбуждений среды (в частности, колебательного процесса) от одной точки к другой.

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

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

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

Очень важным и интересным типом волн яв­ляются волны на поверхности воды. Это один из распространенных видов волн, который каждый наблюдал еще в детстве и который обычно демон­стрируется в рамках школьного курса физики. Од­нако, по выражению Ричарда Фейнмана , "более неудачного примера для демонстрации волн приду­мать трудно, ибо эти волны нисколько не похожи ни на звук, ни на свет; здесь собрались все труднос­ти, которые могут быть в волнах".

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

Если анализировать скорость распространения волны на воде, то оказывается, что она зависит от ее длины. Скорость длинных волн пропорциональна корню квадратному из ускорения свободного паде­ния, умноженному на длину волны. Причиной воз­никновения таких волн является сила тяжести.

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

1.2. Открытие уединенной волны

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

Любопытную волну на воде наблюдал шотланд­ский ученый Джон Скотт Рассел в 1834 году. Он за­нимался исследованием перемещения по каналу баржи, которую тянула пара лошадей. Неожиданно баржа остановилась, но масса воды, которую баржа привела в движение, не остановилась, а собралась у носа судна, а затем оторвалась от него. Далее эта масса воды покатилась по каналу с большой скоро­стью в виде уединенного возвышения, не меняя своей формы и не снижая скорости.

На протяжении всей жизни Рассел неоднократ­но возвращался к наблюдению за этой волной, по­скольку верил, что открытая им уединенная волна играет важную роль во многих явлениях в природе. Он установил некоторые свойства этой волны. Во-первых, заметил, что она движется с постоянной скоростью и без изменения формы . Во-вторых, нашел зависимость скорости С этой волны от глу­бины канала h и высоты волны а:

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

Работа Рассела, опубликованная в 1844 году как "Доклад о волнах", вызвала осторожную реакцию в среде ученых. На континенте ее не заметили сов­сем, а в самой Англии на нее обратили внимание Г.Р. Эйри и Дж.Г. Стоке. Эйри подверг критике ре­зультаты экспериментов, которые наблюдал Рассел. Он отмечал, что из теории длинных волн на мелкой воде выводы Рассела не получаются, и утверждал, что длинные волны не могут сохранять неизменную форму. И в конечном итоге подверг сомнению пра­вильность наблюдений Рассела. Один из основате­лей современной гидродинамики, Джордж Габриэль Стоке, также не согласился с результатами наблюде­ний, полученными Расселом, и критически отнесся к факту существования уединенной волны.

После столь негативного отношения к откры­тию уединенной волны долгое время о ней просто не вспоминали. Определенную ясность в наблюде­ния Рассела внесли Дж. Буссинеск (1872 год) и Дж.У. Рэлей (1876 год), которые независимо друг от друга нашли аналитическую формулу для возвыше­ния свободной поверхности на воде в виде квадрата гиперболического секанса и вычислили скорость распространения уединенной волны на воде.

Позже опыты Рассела были повторены другими исследователями и получили подтверждение.

1.3. Линейные и нелинейные волны

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

u tt = c 2 u xx (1.1)

Характеристика волны и в этом уравнении зависит от пространственной координаты х и времени t , а индексы у переменной и обозначают вторую произ­водную от и по времени (u tt ) и вторую производную от и по переменной x (u xx ). Уравнение (1) описывает плоскую одномерную волну, аналогом которой мо­жет служить волна в струне. В этом уравнении в ка­честве и можно принять плотность воздуха, если речь идет, например, о звуковой волне в воздухе. Ес­ли рассматривают электромагнитные волны, то под и следует понимать напряженность электрического или магнитного поля.

Решение волнового уравнения (1), которое впервые было получено Ж. Д"Аламбером в 1748 го­ду, имеет вид

u(x,t)=f(x-ct)+g(x+ct) (1.2)

Здесь функции f и g находят из начальных условий для и. Уравнение (1.1) содержит вторую производную от и по t , поэтому для него следует задавать два на­чальных условия: значение и при t = 0 и производ­ную и, при t = 0.

Волновое уравнение (1.1) имеет очень важное свойство, суть которого заключена в следующем. Оказалось, что если взять два любых решения этого уравнения, то их сумма снова будет решением этого же уравнения. Это свойство отражает принцип су­перпозиции решений уравнения (1.1) и соответствует линейности явления, которое оно описывает. Для нелинейных моделей это свойство не выполняется, что приводит к существенным отличиям протекания процессов в соответствующих моделях. В частности, из выражения для скорости уединенной волны, ко­торую наблюдал Рассел, следует, что ее значение за­висит от амплитуды, а для волны, описываемой уравнением (1.1), такой зависимости нет.

Непосредственной подстановкой в уравнение (1.1) можно убедиться, что зависимость

u(x,t)=a cos(kx-  t) (1.3)

в которой а, k и  - постоянные, при  =± k является решением уравнения (1). В этом решении а - амплитуда, k - волновое число, а  - частота. При­веденное решение представляет собой монохрома­тическую волну, переносимую в среде с фазовой скоростью

c p = (1.4)

На практике монохроматическую волну создать трудно, и обычно имеют дело с цугом (пакетом) волн, в котором каждая волна распространяется со своей скоростью, а скорость распространения паке­та характеризуется групповой скоростью

C g = , (1.5)

определяемой через производную от частоты  по волновому числу k .

Определить, с какой (линейной или нелиней­ной) моделью имеет дело исследователь, не всегда легко, но когда математическая модель сформули­рована, то решение этого вопроса упрощается и вы­полнение принципа суперпозиции решений можно проверить.

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

Отметим, что и распространение звука не во всех случаях описывается линейным уравнением. Еще Рассел при обосновании своих наблюдений по уе­диненной волне отметил, что звук от выстрела пуш­ки распространяется в воздухе быстрее, чем коман­да произвести этот выстрел. Это объясняется тем, что распространение мощного звука описывается уже не волновым уравнением, а уравнениями газо­вой динамики.

1. Уравнение Кортевега - де Фриса

Окончательная ясность в проблеме, которая воз­никла после опытов Рассела по уединенной волне, наступила после работы датских ученых Д.Д. Кортевега и Г. де Фриса, которые попытались разобраться в существе наблюдений Рассела. Обобщив метод Рэлея, эти ученые в 1895 году вывели уравнение для описания длинных волн на воде. Кортевег и де Фрис, используя уравнения гидродинамики, рас­смотрели отклонение и(х, t ) от положения равнове­сия поверхности воды при отсутствии вихрей и при постоянстве плотности воды. Сделанные ими на­чальные приближения были естественны. Они так­же предположили, что при распространении волны выполняются два условия для безразмерных пара­метров

= <<1,  = (2.1)

Здесь а - амплитуда волны, h - глубина бассейна, в котором рассматриваются волны, l - длина волны (рис. 1).

Суть приближений состояла в том, что амплиту­да рассматриваемых волн была много меньше, чем

Рис. 1. Уединенная волна, распространяющаяся по каналу, и ее параметры

глубина бассейна, но в то же время длина волны бы­ла много больше, чем глубина бассейна. Таким образом, Кортевег и де Фрис рассматривали длин­ные волны.

Уравнение, которое было ими получено, имеет вид

u t + 6uu x + u xxx = 0. (2.2)

Здесь u (x,t) - отклонение от положения равновесия поверхности воды (форма волны) - зависит от ко­ординаты x и времени t . Индексы у характеристики u означают соответствующие производные по t и по x . Это уравнение, как и (1), является уравнением в ча­стных производных. Изучаемая характеристика у него (в данном случае u ) зависит от пространствен­ной координаты x и времени t .

Решить уравнение такого типа - значит найти зависимость u от x и t, после подстановки которой в уравнение мы придем к тождеству.

Уравнение (2.2) имеет волновое решение, извест­ное с конца прошлого века. Оно выражается через специальную эллиптическую функцию, изученную Карлом Якоби, которая носит теперь его имя.

При некоторых условиях эллиптическая функ­ция Якоби переходит в гиперболический секанс и решение имеет вид

u(x,t)=2k 2 ch -2 {k(x-4k 2 t)+  0 } , (2.3)

где  0 - произвольная постоянная.

Решение (8) уравнения (7) является предельным случаем бесконечно большого периода волны. Именно этот предельный случай является уединен­ной волной, соответствующей наблюдению Рассела в 1834 году.

Решение (8) уравнения Кортевега- де Фриса яв­ляется бегущей волной. Это означает, что оно зави­сит от координаты x и времени t через переменную  = x - c 0 t . Эта переменная характеризует положение точки координат, движущейся со скоростью волны с0, то есть она обозначает положение наблюдателя, который постоянно находится на гребне волны. Та­ким образом, уравнение Кортевега- де Фриса в от­личие от решения Д"Аламбера (1.2) волнового реше­ния (1.1) имеет волну, распространяющуюся лишь в одном направлении. Однако оно учитывает прояв­ление более сложных эффектов вследствие дополнительных слагаемых uu x и u xxx .

В действительности это уравнение является так­же приближенным, поскольку при его выводе ис­пользованы малые параметры (2.1)  и  . Если прене­бречь влиянием этих параметров, устремляя их к нулю, мы получим одну из частей решения Д"Алам­бера.

Конечно, при выводе уравнения для длинных волн на воде влияние параметров е и 6 может быть учтено более точно, но тогда получится уравнение, содержащее гораздо больше слагаемых, чем уравне­ние (2.2), и с производными более высокого порядка. Из сказанного следует, что решение уравнения Кортевега-де Фриса для описания волн справедливо только на определенном расстоянии от места обра­зования волны и на определенном промежутке вре­мени. На очень больших расстояниях нелинейные волны уже не будут описываться уравнением Кортевега-де Фриса, и для описания процесса потребует­ся более точная модель. Уравнение Кортевега-де Фриса в этом смысле следует рассматривать как не­которое приближение (математическую модель), со­ответствующее с определенной степенью точности реальному процессу распространения волн на воде.

Используя специальный подход, можно убе­диться, что принцип суперпозиции решений для уравнения Кортевега-де Фриса не выполняется, и поэтому это уравнение является нелинейным и описывает нелинейные волны.

2.1. Солитоны Кортевега - де Фриса

В настоящее время кажется странным, что от­крытие Рассела и его последующее подтверждение в работе Кортевега и де Фриса не получили замет­ного резонанса в науке. Эти работы оказались за­бытыми почти на 70 лет. Один из авторов уравне­ния, Д.Д. Кортевег, прожил долгую жизнь и был известным ученым. Но когда в 1945 году научная общественность отмечала его 100-летний юбилей, то в списке лучших публикаций работа, выполнен­ная им с де Фрисом, даже не значилась. Составите­ли списка сочли эту работу Кортевега не заслужива­ющей внимания. Только спустя еще четверть века именно эта работа стала считаться главным науч­ным достижением Кортевега.

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

Возвращение к открытию уединенной волны на воде произошло в какой-то степени случайно и вна­чале, казалось, не имело к нему никакого отноше­ния. Виновником этого события стал величайший физик нашего столетия Энрико Ферми. В 1952 году Ферми попросил двух молодых физиков С. Улама и Д. Паста решить одну из нелинейных задач на ЭВМ. Они должны были рассчитать колебания 64 гру­зиков, связанных друг с другом пружинками, ко­торые при отклонении от положения равновесия на  l приобретали возвращающуюся силу, равную k  l +a ( l ) 2 . Здесь k и a - постоянные коэффициен­ты. При этом нелинейная добавка предполагалась малой по сравнению с основной силой k  l . Созда­вая начальное колебание, исследователи хотели по­смотреть, как эта начальная мода будет распреде­ляться по всем другим модам. После проведения расчетов этой задачи на ЭВМ ожидаемого результа­та они не получили, но обнаружили, что перекачи­вание энергии в две или три моды на начальном этапе расчета действительно происходит, но затем наблюдается возврат к начальному состоянию. Об этом парадоксе, связанном с возвратом начального колебания, стало известно нескольким математи­кам и физикам. В частности, об этой задаче узнали американские физики М. Крускал и Н. Забуски, ко­торые решили продолжить вычислительные экспе­рименты с моделью, предложенной Ферми.

После расчетов и поиска аналогий эти ученые установили, что уравнение, которое использовали Ферми, Паста и Улам, при уменьшении расстояния между грузиками и при неограниченном росте их числа переходит в уравнение Кортевега-де Фриса. То есть по существу задача, предложенная Ферми, сводилась к численному решению уравнения Кор­тевега-де Фриса, предложенного в 1895 году для описания уединенной волны Рассела. Примерно в те же годы было показано, что для описания ионно-звуковых волн в плазме используется также уравне­ние Кортевега-де Фриса. Тогда стало ясно, что это уравнение встречается во многих областях физики и, следовательно, уединенная волна, которая опи­сывается этим уравнением, является широко рас­пространенным явлением.

Продолжая вычислительные эксперименты по моделированию распространения таких волн, Крус­кал и Забуски рассмотрели их столкновение. Оста­новимся подробнее на обсуждении этого замеча­тельного факта. Пусть имеются две уединенные волны, описываемые уравнением Кортевега-де Фриса, которые различаются амплитудами и дви­жутся друг за другом в одном направлении (рис. 2). Из формулы для уединенных волн (8) следует, что скорость движения таких волн тем выше, чем боль­ше их амплитуда, а ширина пика уменьшается с ростом амплитуды. Таким образом, высокие уеди­ненные волны движутся быстрее. Волна с большей амплитудой догонит движущуюся впереди волну с меньшей амплитудой. Далее в течение некоторого времени две волны будут двигаться вместе как еди­ное целое, взаимодействуя между собой, а затем они разъединятся. Замечательным свойством этих-волн является то, что после своего взаимодействия форма и

Рис. 2. Два солитона, описываемые уравнением Кортевега-де Фриса,

до взаимодействия (вверху) и после (внизу)

скорость этих волн восстанавливаются. Обе волны после столкновения лишь смещаются на не­которое расстояние по сравнению с тем, как если бы они двигались без взаимодействия.

Процесс, у которого после взаимодействия волн сохраняются форма и скорость, напоминает упру­гое столкновение двух частиц. Поэтому Крускал и Забуски такие уединенные волны назвали солитонами (от англ. solitary - уединенный). Это специ­альное название уединенных волн, созвучное элек­трону, протону и многим другим элементарным частицам, в настоящее время общепринято.

Уединенные волны, которые были открыты Рас­селом, и в самом деле ведут себя как частицы. Боль­шая волна не проходит через малую при их взаимо­действии. Когда уединенные волны соприкасаются, то большая волна замедляется и уменьшается, а волна, которая была малой, наоборот, ускоряется и подрастает. И когда малая волна дорастает до разме­ров большой, а большая уменьшается до размеров малой, солитоны разделяются и больший уходит вперед. Таким образом, солитоны ведут себя как уп­ругие теннисные мячи.

Дадим определение солитона . Солитоном на­зывается нелинейная уединенная волна, которая сохраняет свою форму и скорость при собственном движении и столкновении с себе подобными уеди­ненными волнами, то есть представляет собой ус­тойчивое образование. Единственным результатом взаимодействия солитонов может быть некоторый сдвиг фаз.

Открытия, связанные с уравнением Кортевега - де Фриса, не закончились открытием солитона. Следующим важным шагом, имеющим отношение к этому замечательному уравнению, было создание нового метода решения нелинейных уравнений в частных производных. Хорошо известно, что най­ти решения нелинейных уравнений очень сложно. До 60-х годов нашего столетия считалось, что такие уравнения могут иметь только некоторые частные решения, удовлетворяющие специально заданным начальным условиям. Однако уравнение Кортевега-де Фриса и в этом случае оказалось в исключи­тельном положении.

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

2.2. Групповой солитон

Выше мы говорили, что на практике волны, как правило, распространяются группами. Подобные группы волн на воде люди наблюдали с незапамят­ных времен. На вопрос о том, почему для волн на воде так типичны "стаи" волн, удалось ответить Т. Бенжамену и Дж. Фейеру только в 1967 году. Тео­ретическими расчетами они показали, что простая периодическая волна на глубокой воде неустойчива (теперь это явление называется неустойчивостью Бенжамена-Фейера), и поэтому волны на воде из-за неустойчивости разбиваются на группы. Уравнение, с помощью которого описывается распространение групп волн на воде, было получено В.Е. Захаровым в 1968 году. К тому времени это уравнение уже было известно в физике и носило название нелинейного уравнения Шрёдингера. В 1971 году В.Е. Захаров и А.Б. Шабат показали, что это нелинейное уравне­ние имеет решения также в виде солитонов, более того, нелинейное уравнение Шрёдингера, так же как и уравнение Кортевега-де Фриса, может быть проинтегрировано методом обратной задачи рассея­ния. Солитоны нелинейного уравнения Шрёдинге­ра отличаются от обсуждаемых выше солитонов Кортевега-де Фриса тем, что они соответствуют форме огибающей группы волн. Внешне они на­поминают модулированные радиоволны. Эти солитоны называются групповыми солитонами, а иногда солитонами огибающей. Это название от­ражает сохраняемость при взаимодействии огиба­ющей волнового пакета (аналог штриховой ли­нии, представленной на рис. 3), хотя сами волны под огибающей двигаются со скоростью, отличной от групповой. При этом форма огибающей описывается

Рис. 3. Пример группового солитона (штриховая линия)

зависимостью

a(x,t)=a 0 ch -1 (
)

где а а - амплитуда, а l - половина размера солитона. Обычно под огибающей солитона находится от 14 до 20 волн, причем средняя волна самая большая. С этим связан хорошо известный факт, что самая вы­сокая волна в группе на воде находится между седь­мой и десятой (девятый вал). Если в группе волн об­разовалось большее количество волн, то произойдет ее распад на несколько групп.

Нелинейное уравнение Шрёдингера, как и урав­нение Кортевега- де Фриса, также имеет широкую распространенность при описании волн в различ­ных областях физики. Это уравнение было предло­жено в 1926 году выдающимся австрийским физи­ком Э. Шрёдингером для анализа фундаментальных свойств квантовых систем и первоначально ис­пользовано при описании взаимодействия внут­риатомных частиц. Обобщенное или нелинейное уравнение Шрёдингера описывает совокупность явлений в физике волновых процессов. Например, оно используется для описания эффекта самофоку­сировки при воздействии мощного лазерного луча на нелинейную диэлектрическую среду и для опи­сания распространения нелинейных волн в плазме.

3. Постановка задачи

3.1. Описание модели.В настоящее время наблюдается значи­тельно возрастающий интерес к исследованию нелинейных волно­вых процессов в различных областях физики (например, в оптике, физике плазмы, радиофизике, гидродинамике и т.д.). Для изучения волн малой, но конечной амплитуды в дисперсионных средах в каче­стве модельного уравнения часто используют уравнение Кортевега-де Фриза (КдФ):

u t + ии х +  и ххх = 0 (3.1)

Уравнение КдФ было использовано для описания магнитозвуковых волн, распространяющихся строго поперек магнитного поля или под углами, близкими к .

Основные предположения, которые делаются при выводе уравне­ния: 1) малая, но конечная амплитуда, 2) длина волны велика по сравнению с длиной дисперсии.

Компенсируя действие нелинейности, дисперсия дает возможность формироваться в дисперсионной среде стационарным волнам конеч­ной амплитуды - уединенным и периодическим. Уединенные волны для уравнения КдФ после работы стали называться солитонами . Периодические волны носят название кноидальных волн. Соот­ветствующие формулы для их описания даны в .

3.2. Постановка дифференциальной задачи.В работе иссле­дуется численное решение задачи Коши для уравнения Кортевега-де Фриза с периодическими условиями по пространству в прямоуголь­нике Q T ={(t , x ):0< t < T , x  [0, l ].

u t + ии х +  и ххх = 0 (3.2)

u(x,t)| x=0 =u(x,t)| x=l (3.3)

с начальным условием

u(x,t)| t=0 =u 0 (x) (3.4)

4. Свойства уравнения Кортевега - де Фриза

4.1. Краткий обзор результатов по уравнению КдФ.Задача Коши для уравнения КдФ при различных предположениях отно­сительно u 0 (х) рассматривалась во многих работах . Задача о существовании и единственности решения с условиями периодично­сти в качестве краевых условий была решена в работе с помощью метода конечных разностей. Позже, при менее сильных предположе­ниях, существование и единственность были доказана в статье в пространстве L  (0,T ,H s (R 1)), где s>3/2, а в случае периодической задачи - в пространстве L  (0,T ,H  (C ))где С - окружность дли­ны, равной периоду, на русском языке эти результаты представлены в книге .

Случай, когда не предполагается какая-либо гладкость началь­ной функции u 0  L 2 (R 1 ) , рассмотрен в работе . Там вводит­ся понятие обобщенного решения задачи (3.2),(3.4), устанавливает­ся существование обобщенного решения и(t ,х)  L  (0, T , L 2 (R 1 )) в случае произвольной начальной функции u 0  L 2 (R 1 ) ; при этом и(t ,х)  L 2 (0,Т;H -1 (- r , r )) для любого r>0 , и если для некото­рого  > 0 (x  u 0 2 (x ))  L 1 (0,+  ) , то

(4.1)

Используя обращение линейной части уравнения при помощи фун­даментального решения G (t,x) соответствующего линейного опера­тора
, вводится класс корректности задачи (3.2),(1.4) и уста­навливаются теоремы единственности и непрерывной зависимости решений этой задачи от начальных данных. Также исследуются во­просы регулярности обобщенных решений. Одним из основных ре­зультатов является достаточное условие существования непрерыв­ной по Гельдеру при t > 0 производной
в терминах существования моментов для начальной функции, для любых k и l .

Задача Коши для уравнения КдФ исследовалась также методом обратной задачи рассеяния, предложенном в работе . При по­мощи этого метода были получены результаты о существовании и гладкости решений при достаточно быстро убывающих начальных функциях, причем в установлен, в частности, результат о раз­решимости задачи (3.2),(3.4) в пространстве C  (О, Т; S(R 1 )) .

Наиболее полный обзор современных результатов по уравнению КдФ можно найти в .

4.2. Законы сохранения для уравнения КдФ. Как известно, для уравнения КдФ существует бесконечное число законов сохране­ ния. В работе приводится строгое доказательство этого факта. В работах , различные законы сохранения применялись для до­ казательства нелокальных теорем существования решения задачи (3.2),(3.4) из соответствующих пространств.

Продемонстрируем вывод первых трех законов сохранения для за­ дачи Коши на R 1 и периодической задачи.

Для получения первого закона сохранения достаточно проинте­ грировать уравнения (3.2) по пространственной переменной. Полу­ чим:

отсюда и следует первый закон сохранения:

Здесь в качестве a и b выступают +  и -  для задачи Коши и границы основного периода для периодической задачи. Поэтому второе и третье слагаемые обращаются в 0.

(4.2)

Для вывода второго закона сохранения следует умножить уравне­ ние (3.2) на 2 u (t,x) и проинтегрировать по пространственной пере­ менной. Тогда, используя формулу интегрирования по частям полу­ чим:

но в силу "краевых" условий все слагаемые кроме первого опять сокращаются

Таким образом второй интегральный закон сохранения имеет вид:

(4.3)

Для вывода третьего закона сохранения нужно умножить наше уравнение (3.2) на (и 2 + 2  и хх ), таким образом получим:

После применения несколько раз интегрирования по частям тре­тий и четвертый интегралы сокращаются. Второе и третье слагае­ мые исчезают из-за граничных условий. Таким образом из первого интеграла получаем:

что эквивалентно

А это и есть третий закон сохранения для уравнения (3.2). Под физическим смыслом первых двух интегральных законов со­ хранения в некоторых моделях можно понимать законы сохранения импульса и энергии, для третьего и последующих законов сохране­ния физический смысл охарактеризовать уже труднее, но с точки зрения математики эти законы дают дополнительную информацию о решении, которая используется потом для доказательств теорем существования и единственности решения, исследования его свойств и вывода априорных оценок.

5. Разностные схемы для решения уравнения КдФ

3.1. Обозначения и постановка разностной задачи. В области ={( x , t ):0  x  l ,0  t  T } обычным образом введем равномерные сетки, где

Введем линейное пространство  h сеточных функций, определен­ных на сетке
со значениями в узлах сетки y i = y h ( x i ). Пред­ полагается, что выполнены условия периодичности y 0 = y N . Кроме того, формально полагаем y i + N = y i для i  1 .

Введем скалярное произведение в пространстве  h

(5.1)

Снабдим линейное пространство П/г нормой:

Поскольку в пространство  h входят периодические функции, то это скалярное произведение эквивалентно скалярному произведе­ нию:

Будем строить разностные схемы для уравнения (3.2) на сетке с периодическими краевыми условиями. Нам потребуются обозна­чения разностных аппроксимаций. Введем их.

Используем стандартные обозначения для решения уравнения на очередном (n -м) временном слое, то есть

Введем обозначения для разностных аппроксимаций производных. Для первой производной по времени:

Аналогично для первой производной по пространству:

Теперь введем обозначения для вторых производных:

Третью пространственную производную будем аппроксимировать следующим образом:

Также нам потребуется аппроксимация у 2 , которую мы обозначим буквой Q и введем следующим образом:

(5.2)

Для записи уравнения на полу целых слоях будем использовать уравновешенную аппроксимацию, т.е.

за исключением аппроксимации у 2 на полу целом слое. Приведем одну из возможных аппроксимаций у 2 на полу целом слое:

Замечание 2. Стоит отметить, что для 1 выполняется равенство:

Определение 1. Следуя разностную схему для уравнения КдФ будем называть консервативной, если для нее имеет место сеточ­ ный аналог первого интегрального закона сохранения, справедливо­

Определение 2. Следуя разностную схему для уравнения КдФ будем называть L 2 -консервативной, если для нее имеет место сеточ­ ный аналог второго интегрального закона сохранения, справедливо­ го для дифференциальной задачи.

5.2. Явные разностные схемы (обзор). При построении раз­ ностных схем будем ориентироваться на простейшую разностную схему из работы для линеаризованного уравнения КдФ, кото­ рое сохраняет свойства самого уравнения КдФ в смысле двух первых законов сохранения.

(5.3)

Исследуем теперь схему (5.4) на свойства консервативности. Вы­ полнение первого закона сохранения очевидно. Достаточно просто умножить это уравнение скалярно на 1. Тогда второе и третье сла­ гаемые схемы (5.4) дадут 0, а от первого останется:

(5.4)

Это сеточный аналог первого закона сохранения.

Для вывода второго закона сохранения умножим скалярно урав­ нение (5.3) на 2  у. Приходим к энергетическому тождеству

(5.5)

Наличие отрицательного дисбаланса говорит не только о невыпол­ нении соответствующего закона сохранения, но и ставит под сомне­ние вопрос вообще об устойчивости схемы в наиболее слабой норме L 2 (). )- В работе показано, что схемы семейства (3.18) являются абсолютно неустойчивыми в норме L 2 ().

В данный момент наиболее популярными схемами для уравнения КдФ считаются трехслойные схемы ввиду их простоты, точности и удобства реализации.

(5.6)

Эту же схему можно представить в виде явной формулы

Самой простой трехслойной схемой является следующая схема:

Эта схема была использована при получении первых численных решений КдФ . Эта схема аппроксимирует дифференциальную задачу с порядком О ( 2 + h 2 ). Согласно , схема является устой­ чивой при выполнении условия (при малых Ь):

Приведем еще несколько схем. Трехслойная явная схема с поряд­ ком аппроксимации O (  2 + h 4 ) :

Третья производная по пространству аппроксимируется на семи­ точечном шаблоне, а первая строится по пяти точкам. Согласно , эта схема устойчива при выполнении условия (при малых h ):

Легко видеть, что для этой схемы с более высоким порядком ап­проксимации условие устойчивости является более жестким.

В работе предлагается следующая явная разностная схема с порядком аппроксимации О( 2 + h 2 ) :

(5.8)

Так как разностное уравнение (5.8) можно записать в дивергент­ ном виде

то, скалярно умножив уравнение (5.9) на 1, получим

следовательно, выполняется соотношение:

которое можно считать сеточным аналогом первого закона сохране­ ния. Таким образом, схема (5.8) является консервативной. В доказано, что схема (5.8) является L 2 -консервативной и ее решение удовлетворяет сеточному аналогу интегрального закона сохранения

5.3. Неявные разностные схемы (обзор). В этом параграфе мы рассмотрим неявные разностные схемы для уравнения Кортевега-де Фриза.

Вариант двухслойной схемы - неявная абсолютно устойчивая схе­ ма с порядком аппроксимации О ( 2 , h 4 ) :

Решение разностной схемы (3.29) вычисляется с помощью семи д иагональной циклической прогонки . Вопрос о консервативности этой схемы не исследовался.

В работе предлагается неявная трехслойная схема с весами:

(5.10)

Разностная схемы (5.10) с периодическими по пространству реше­ниями, консервативна, L 2 -консервативна при  =1/2 и  =1/4 для ее решения имеют место сеточные аналоги интегральных законов сохранения.

6. Численное решение

Численное решение для (3.2), (3.3), (3.4) было проделано с использованием явной схемы

Решалась начально-краевая задача на отрезке . В качестве начальных условий бралась функция

u 0 (x)=sin (x).

Явным образом было получено решение.

Программа для расчетов была написана на языке Turbo Pascal 7.0. Текст основных частей программы прилагается.

Расчеты велись на вычислительной машине с процессором AMD -K 6-2 300 МГц с технологией 3DNOW !, размер оперативной памяти 32 Мб.

7. Заключение

Настоящая работа посвящена исследованию уравнения Кортевега – де Фриза. Проведен обширный литературный обзор по теме исследования. Изучены различные разностные схемы для уравнения КдФ. Выполнен практический счет с использованием явной пяти точечной разносной схемы

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

8. Литература

1. Ландсберг Г.С. Элементарный учебник физики. М.: Наука, 1964. Т. 3.

2. Фейнман Р., Лейтон Р., Сэндс М. Фейнмановские лекции по физике. М.: Мир, 1965. Вып.4.

3. Филиппов А. Г Многоликий солитон. М.: Наука, 1986. (Б-чка "Квант"; Вып. 48).

4. Рубанков В.Н. Солитоны, новое в жизни, науке, тех­нике. М.: Знание, 1983. (Физика; Вып. 12).

5. Korteweg D.J., de Vries G. On the change form of long waves advancing in a rectangular channel and on new type of long stationary waves.//Phyl.May. 1895. e5. P . 422-443.

6. Сагдеев Р.З. Коллективные процессы и ударные волны в разре­женной плазме.-В кн.: Вопросы теории плазмы, Вып.4. М.: Атомиз-дат, 1964, с.20-80.

7. Березин Ю.А., Карпман В.И. К теории нестационарных волн конечной амплитуды в разреженной плазме. // ЖЭТФ, 1964, т.46, вып.5, с. 1880-1890.

8. Zabusky N.J., Kruskal M.D. Interactions of "solitons"in a collisionless plasma and the reccurence of initial states // Phys.Rev.Lett. 1965. V .15. еб. Р.240-243.

9. Буллаф Р., Кодри Ф. Солитоны. М.: Мир; 1983

10. Sjoberg A. On the Korteweg-de Vries equation, existence and uniqueness, Uppsala University, Department of Computers, 1967

11. Temam R. Sur un probleme non lineare // J.Math.Pures Anal. 1969, V.48, 2, P. 159-172.

12. Лионе Ж.-Л. Некоторые методы решения нелинейных краевых задач. М.: Мир, 1972.

13. Кружков С.Н. Фаминский А.В. Обобщенные решения для урав­нения Кортевега-де Фриза.// Матем. сборник, 1983, т. 120(162), еЗ, с.396-445

14.. Gardner C.S., Green J.M., Kruskal M.D., Miura R.M. Method for solving the Korteweg-de Vries equation // Phys.Rev.Lett. 1967. V . 19. P . 1095-1097.

15. Шабат А.Б. Об уравнении Кортевега-де Фриза // ДАН СССР, 1973, т.211, еб, с.1310-1313.

16. Фаминский А.В. Граничные задачи для уравнения Кортевега-де Фриза и его обобщений: Дисс.... докт. физ.-матем. наук,М:РУДН,2001

17. Miura R.M., Gardner C.S., Kruscal M.D. Korteweg-de Vries equation and generlization. II. Existence of conservation laws and constants of motion. // J.Math.Phys. 1968. V .9. P . 1204-1209.

18. Амосов А.А., Злотник А.А. Разностная схема для уравнений движений газа.

19. Самарский А.А., Мажукин В.И., Матус П.П., Михайлик И.А. Z /2-консервативные схемы для уравнения Кортевега-де Фриса.// ДАН, 1997, т.357, е4, с.458-461

20. Березин Ю.А. Моделирование нелинейных волновых процес­сов. Новосибирск: Наука. 1982.

21. Березин Ю.А., О численных решениях уравнения Кортевега-де Вриза.// Численные методы механики сплошной среды. Новоси­бирск, 1973, т.4, е2, с.20-31

22. Самарский А.А., Николаев Методы решения сеточных уравнений. М: Наука, 1978

23. Самарский А.А., Гулин А.В. Численные методы. М: Наука, 1989

24. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. М: Наука, 1987

На теперешнем курсе семинары стали заключаються не в решении задач, а докладах на различную тематику. Думаю, будет верным оставлять их здесь в более или менее популярном виде.

Слово «солитон» происходит от английского solitary wave и означает именно уединенную волну (или говоря языком физики некоторое возбуждение).

Солитон возле острова Молокаи (Гавайский архипелаг)

Цунами - тоже солитон, но значительно более крупный. Уединенность не означает, что волна будет одна единственная на весь мир. Солитоны иногда встречаются группами, как возле Бирмы.

Солитоны в Андаманском море, омывающем берега Бирмы, Бенгалии и Тайланда.

В математическом смысле солитон является решением нелинейного уравнения в частных производных. Означает это следующее. Решать линейные уравнения что обыкновенные из школы, что дифференциальные человечество уже умеет достаточно давно. Но стоит возникнуть квадрату, кубу или еще более хитрой зависимости в дифференциальном уравнении от неизвестной величины и наработанный за все века математический аппарат терпит фиаско - человек пока не научился их решать и решения чаще всего угадываются или подбираются из различных соображений. Но Природу описывают именно они. Так нелинейные зависимости рождают практически все явления, чарующие глаз, да и позволяющие существовать жизни тоже. Радуга в своей математической глубине описывается функцией Ейри (правда, говорящая фамилия для ученого, чье исследование рассказывает о радуге?)

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

Функции, объекты или явления, называющиеся солитонами, имеют два важных свойства: они стабильны во времени и сохраняют свою форму. Конечно, в жизни никто и ничто бесконечно долго им удовлетворять не будет, поэтому нужно сравнивать с аналогичными явлениями. Вернувшись к морской глади, рябь на её поверхности возникает и исчезает за доли секунды, большие волны, вздымаемые ветром взлетают и рассыпаются брызгами. Но цунами движется глухой стеной на сотни километров не теряя заметно в высоте волны и силе.

Есть несколько типов уравнений, приводящих к солитонам. Прежде всего, это задача Штурма-Лиувилля

В квантовой теории это уравнение известно под названием нелинейного уравнения Шредингера (Schrödinger) если функция имеет произвольный вид. В этой записи число называют собственным. Оно такое особенное, что его тоже находят при решении задачи, потому как не каждое его значение может дать решение. Роль собственных чисел в физике очень велика. Например, энергия является собственным числом в квантовой механике, переходы между различными системами координат так же не обходятся без них. Если потребовать, чтобы изменение параметра t в не изменяли собственные числа (а t может быть временем, например, или каким-то внешним влиянием на физическую систему), то придем к уравнению Кортевега-де Фриза (Korteweg-de Vries):

Есть и иные уравнения, но сейчас они не так важны.

В оптике фундаментальную роль играет явление дисперсии - зависимость частоты волны от её длины , а точнее так называемого волнового числа :

В простейшем случае она может быть линейна (, где - скорость света). В жизни ж часто получаем квадрат волнового числа, а то и что-то более хитрое. На практике, дисперсия ограничивает пропускную возможность оптоволокна, по которому только что бежали эти слова к вашему интернет-провайдеру с серверов WordPress’а. Но так же она позволяет пропускать по одному оптоволокну не один луч, а несколько. И в терминах оптики приведенные выше уравнения рассматривают простейшие случаи дисперсии.

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

В оптике так же говорят про временные и пространственные солитоны. Из названия становится ясно, будем мы наблюдать солитон как некую волну в пространстве, или же это будет всплеск во времени. Временные возникают из-за балансировки нелинейных эффектов дифракцией - отклонения лучей от прямолинейного распространения. Например, посветили лазером в стекло (оптоволокно), и внутри лазерного луча показатель преломления стал зависеть от мощности лазера. Пространственные солитоны возникают из-за балансировки нелинейностей дисперсией.

Фундаментальный солитон

Как уже говорилось, широкополосность (то есть возможность передать много частот, а значит и полезной информации) волоконно-оптических линий связи ограничивается нелинейными эффектами и дисперсией, меняющими амплитуду сигналов и их частоту. Но с другой стороны, те же самые нелинейность и дисперсия могут привести к созданию солитонов, которые сохраняют свою форму и иные параметры существенно дольше чем все остальное. Естественным выводом отсюда является желание использовать сам солитон в качестве информационного сигнала (есть вспышка-солитон на конце волокна - передали единичку, нет - передали нолик).

Пример с лазером, изменяющим коэффициент преломления внутри оптоволокна по мере своего распространения достаточно жизненный, особенно если «запихнуть» в волокно тоньше человеческого волоса импульс в несколько ватт. Для сравнения много это или нет, типичная энергосберегающая лампочка мощностью в 9 Вт освещает письменный стол, но при этом размером с ладонь. В общем, мы не отойдем далеко от действительности предположив, что зависимость коэффициента преломления от мощности импульса внутри волокна будет выглядеть так:

После физических размышлений и математических преобразований различной сложности на амплитуду электрического поля внутри волокна можно получить уравнение вида

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

Распределение интенсивности в поперечном сечении лазерного луча в форме фундаментального солитона.

Именно это решение и называется фундаментальным солитоном. Мнимая экспонента определяет распространение солитона вдоль оси волокна. На практике это все означает, что посветив на стенку мы увидели б яркое пятно в центре, интенсивность которого быстро спадала бы на краях.

Фундаментальный солитон как и все солитоны, возникающие с использованием лазеров, имеет определенные особенности. Во-первых, если мощность лазера окажется недостаточной, он не появится. Во-вторых, даже если где-то слесарь излишне перегнет волокно, капнет на него маслом или сделает иную пакость, солитон проходя сквозь поврежденную область возмутится (в физическом и переносном смыслах), но быстро вернется к своим изначальным параметрам. Люди и иные живые существа так же попадают под определение автосолитона и это умение возвращаться в спокойное состояние очень важно в жизни 😉

Потоки энергии внутри фундаментального солитона выглядят так:

Направление потоков энергии внутри фундаментального солитона.

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

На практике можно получить несколько солитонов, если лазер имеет несколько каналов генерации, параллельных его оси. Тогда взаимодействие солитонов будет определяться степенью перекрытия их «юбок». Если рассеяние энергии не очень велико, можно считать, что потоки энергии внутри каждого солитона сохраняются во времени. Тогда солитоны начинают кружиться и сцепляться вместе. На следующем рисунке приведено моделирование столкновения двух троек солитонов.

Моделирование столкновения солитонов. На сером фоне изображены амплитуды (как рельеф), а на черном - распределение фазы.

Группы солитонов встречаются, цепляются и образуя Z-подобную структуру начинают вращаться. Еще более интересные результаты можно получить нарушением симметрии. Если расставить лазерные солитоны в шахматном порядке и выбросить один, структура начнет вращаться.

Нарушение симметрии в группе солитонов приводит к вращению центра инерции структуры в направлении стрелки на рис. справа и вращению вокруг мгновенного положения центра инерции

Вращений будет два. Центр инерции будет обращаться против часовой стрелки, а так же сама структура будет крутиться вокруг его положения в каждый момент времени. При чем периоды вращений будут равны, например, как у Земли и Луны, которая повернута к нашей планете лишь одной стороной.

Эксперименты

Столь необычные свойства солитонов обращают на себя внимание и заставляют задуматься о практическом применении уже около 40 лет. Сразу можно сказать, что солитоны можно использовать для сжатия импульсов. На сегодняшний день так можно получить длительность импульса до 6 фемтосекунд ( сек или дважды брать от секунды одну миллионную и результат поделить на тысячу). Отдельный интерес представляют солитонные линии связи, разработка которых идет уже довольно давно. Так Хасегавой было предложено следующую схему еще в 1983 году.

Солитонная линия связи.

Линия связи формируется из секций длиной около 50 км. Всего длина линии составляла 600 км. Каждая секция состоит из приемника с лазером передающих в следующий волновод усиленный сигнал, что позволило достичь скорости 160 Гбит/сек.

Презентация

Литература

1. Дж. Лем. Введение в теорию солитонов. Пер. с англ. М.: Мир, - 1983. -294 с.
2. Дж. Уизем Линейные и нелинейные волны. - М.: Мир, 1977. - 624 с.
3. И. Р. Шен. Принципы нелинейной оптики: Пер. с англ./Под ред. С. А. Ахманова. - М.: Наука., 1989. - 560 с.
4. С. А. Булгакова, А. Л. Дмитриев. Нелинейно-оптические устройства обработки информации// Учебное пособие. - СПб: СПбГУИТМО, 2009. - 56 с.
5. Werner Alpers et. al. Observation of Internal Waves in the Andaman Sea by ERS SAR // Earthnet Online
6. А. И. Латкин, А. В. Якасов. Автосолитонные режимы распространения импульса в волоконно-оптической линии связи с нелинейными кольцевыми зеркалами // Автометрия, 4 (2004), т.40.
7. Н. Н. Розанов. Мир лазерных солитонов // Природа, 6 (2006). С. 51-60.
8. О. А. Татаркина. Некоторые аспекты проектирования солитонных волоконно-оптических систем передачи // Фундаментальные исследования, 1 (2006), С. 83-84.

P. S. О диаграммах в .