После расчетов и поиска аналогий эти ученые установили, что уравнение, которое использовали Ферми, Паста и Улам, при уменьшении расстояния между грузиками и при неограниченном росте их числа переходит в уравнение Кортевега-де Фриса. То есть по существу задача, предложенная Ферми, сводилась к численному решению уравнения Кортевега-де Фриса, предложенного в 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 + ии х + b и ххх = 0 (3.1)
Уравнение КдФ было использовано для описания магнитозвуковых волн, распространяющихся строго поперек магнитного поля или под углами, близкими к
.Основные предположения, которые делаются при выводе уравнения: 1) малая, но конечная амплитуда, 2) длина волны велика по сравнению с длиной дисперсии.
Компенсируя действие нелинейности, дисперсия дает возможность формироваться в дисперсионной среде стационарным волнам конечной амплитуды - уединенным и периодическим. Уединенные волны для уравнения КдФ после работы стали называться солитонами . Периодические волны носят название кноидальных волн. Соответствующие формулы для их описания даны в .
3.2. Постановка дифференциальной задачи.В работе исследуется численное решение задачи Коши для уравнения Кортевега-де Фриза с периодическими условиями по пространству в прямоугольнике Q T ={( t , x ):0< t < T , x Î [0, l ].
u t + ии х + b и ххх = 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))где С - окружность длины, равной периоду, на русском языке эти результаты представлены в книге .
Формат: doc
Дата создания: 31.05.2003
Размер: 125.1 KB
Скачать реферат 1. Введение | |
1.1. Волны в природе | |
2. Уравнение Кортевега - де Фриса | |
2.2. Групповой солитон | |
3. Постановка задачи | |
3.1. Описание модели | |
3.2. Постановка дифференциальной задачи. | |
4. Свойства уравнения Кортевега - де Фриза | |
4.1. Краткий обзор результатов по уравнению КдФ | |
4.2. Законы сохранения для уравнения КдФ | |
5. Разностные схемы для решения уравнения КдФ | |
5.1. Обозначения и постановка разностной задачи. | |
5.2. Явные разностные схемы (обзор) | |
5.3 Неявные разностные схемы (обзор). | |
6.Численное решение | |
7. Заключение | |
8. Литература | |
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 .
Определить, с какой (линейной или нелинейной) моделью имеет дело исследователь, не всегда легко, но когда математическая модель сформулирована, то решение этого вопроса упрощается и выполнение принципа суперпозиции решений можно проверить.
Возвращаясь к волнам на воде, заметим, что их можно анализировать используя хорошо известные уравнения гидродинамики, о которых известно, что они нелинейны. Поэтому и волны на воде в общем случае являются нелинейными. Только в предельном случае малых амплитуд эти волны могут считаться линейными.
Отметим, что и распространение звука не во всех случаях описывается линейным уравнением. Еще Рассел при обосновании своих наблюдений по уединенной волне отметил, что звук от выстрела пушки распространяется в воздухе быстрее, чем команда произвести этот выстрел. Это объясняется тем, что распространение мощного звука описывается уже не волновым уравнением, а уравнениями газовой динамики.
Уравнение Кортевега - де Фриса
Окончательная ясность в проблеме, которая возникла после опытов Рассела по уединенной волне, наступила после работы датских ученых Д.Д. Кортевега и Г. де Фриса, которые попытались разобраться в существе наблюдений Рассела. Обобщив метод Рэлея, эти ученые в 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)
Эту же схему можно представить в виде явной формулы
(5.7)Самой простой трехслойной схемой является следующая схема:
Эта схема была использована при получении первых численных решений КдФ . Эта схема аппроксимирует дифференциальную задачу с порядком О ( 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
), к-рое в каждый момент времени локализовано в конечной области пространства и относительно медленно изменяет свою структуру при распространении.
Примеры уединённых волн: а - стационарное возвышение (солитон) на мелкой воде; h - смещение поверхности жидкости; б - небольшой амплитуды в газе; р - изменение давления; в - возбуждения в аксоне нерва; и - мембраны. По оси абсцисс отложена переменная
Типичная У. в. имеет вид одиночного импульса или перепада (рис.), но У. в. может иметь и более сложную структуру.
В более узком смысле под У. в. понимают локализованную стационарную нелинейную волну, распространяющуюся без изменения формы с постоянной скоростью и описываемую ур-ниями в обыкновенных производных. В фазовом пространстве У. в. отвечает , соединяющая две различные точки равновесия или возвращающаяся в ту же самую точку. К У. в. относят, напр., такие типы нелинейных волн, как ударные волны в диссипативной среде, стационарные импульсные волны возбуждения в активных средах (напр., ) и в среде без потерь.
Физический энциклопедический словарь. - М.: Советская энциклопедия . . 1983 .
УЕДИНЁННАЯ ВОЛНА
Волновое движение (см. Волны),
к-рое в каждый момент времени локализовано в конечной области пространства и достаточно быстро убывает с удалением от этой области. Типичная У. в. имеет вид одиночного импульса или перепада (рис.), но У. в. может иметь и более сложную структуру.
В более узком смысле под У. в. понимают локализованную стационарную нелинейную волну, распространяющуюся без изменения формы с пост. скоростью и описываемую ур-ниями в обыкновенных производных. В фазовом пространстве У. в. отвечает траектория, соединяющая две разл. точки равновесия или возвращающаяся в ту же самую точку. К У. в. относят, напр., такие типы нелинейных волн, как ударные волны в диссипативной среде, стационарные импульсные волны возбуждения в активных средах (напр., нервный импульс) и солитон в среде без потерь. Лит. см. при ст. Солитон. Л. А. Островский.
Примеры уединённых волн: а - стационарное возвышение (соли-тон) на мелкой воде; h - смещение поверхности жидкости; б - ударная волна небольшой амплитуды в газе; p - изменение давления; в - импульс возбуждения в аксоне нерва; и - потенциал мембраны. По оси абсцисс отложена переменная где t - время, x -координата, u- скорость уединённой волны.
Физическая энциклопедия. В 5-ти томах. - М.: Советская энциклопедия . Главный редактор А. М. Прохоров . 1988 .
Смотреть что такое "УЕДИНЕННАЯ ВОЛНА" в других словарях:
- (уединенная волна), структурно устойчивая уединенная волна, которая, распространяясь, не расширяется и сохраняет свою форму и скорость. Солитоны ведут себя, как частицы. Они важны во многих областях МЕХАНИКИ ТЕКУЧИХ СРЕД, а также ФИЗИКИ ТВЕРДОГО… … Научно-технический энциклопедический словарь
Структурно устойчивая уединённая волна, распространяющаяся в нелинейной среде. Солитоны ведут себя подобно частицам (частицеподобная волна): при взаимодействии друг с другом или с некоторыми другими возмущениями они не разрушаются, а расходятся,… … Энциклопедический словарь
Структурно устойчивая уединенная волна, распространяющаяся в нелинейной среде. Солитоны ведут себя подобно частицам (частицеподобная волна): при взаимодействии друг с другом или с некоторыми другими возмущениями они не разрушаются, а расходятся,… … Большой Энциклопедический словарь
Солитон - структурно устойчивая уединенная волна, распространяющаяся в нелинейной среде, которая может характеризоваться как частицеподобная волна, частица … Начала современного естествознания
1) Л. т. в д е с к р и п т и в н о й теории множеств: топологич. отображение между двумя множествами в можно продолжить до гомеоморфизма нек рых содержащих их множеств типа Следствием этой Л. т. является топологич. инвариантность хаусдорфова типа … Математическая энциклопедия
Здесь описаны В.: а) водяные, б) воздушные звуковые, в) световые, г) электрические волны и д) математическая теория В. А) Волны в воде обыкновенно являются следствием косвенного удара ветра о воду. Поверхность воды от этого делается вогнутой, но… … Энциклопедический словарь Ф.А. Брокгауза и И.А. Ефрона
Рассмотрим среду без диссипации Пусть пока нелинейность в среде квадратична, т. е. тогда вместо (19.1) будем искать уравнение, полученное Кортевегом и де Вризом для волн на поверхности жидкости:
Решения этого уравнения сейчас изучены очень подробно, в том числе и нестационарные, но мы будем обсуждать только самые простые из них, дополнив обсуждение качественными соображениями. Прежде всего поразмышляем над тем, к чему может привести добавление к уравнению простой волны слагаемого, описывающего дисперсионное расплывание. Как мы уже знаем, дисперсионное расплывание может компенсировать процесс опрокидывания волны, и тогда ее профиль стабилизируется, т. е. возможно существование стационарных бегущих волн, профиль которых не меняется во времени. Такие волны определены во всем пространстве и бегут с постоянной скоростью V, т. е. все переменные в волне являются функцией бегущей координаты Для них т. е. стационарные волны уравнения (19.14) описываются уравнением в обыкновенных производных или после интегрирования,
Таким образом, стационарным волнам уравнения Кортевега-де Вриза соответствует уравнение консервативного нелинейного осциллятора. Постоянную будем считать равной нулю (это всегда можно сделать, введя полую переменную), тогда уравнение (19.15) представляется в виде где Потенциальная энергия стационарных волн и их фазовый портрет приведены на рис. 19.6.
Существуют различные классы решений уравнения Кортевега-де Вриза. Можно выделить два из них.
1. Квазисинусоидальные колебания с малыми амплитудами (фазовые траектории вблизи состояния центра); для них нелинейность почти не сказывается (рис. 19.7 а).
2. Движение вблизи сепаратрисы и по самой сепаратрисе. Именно эти сильно нелинейные волны и представляют для нас интерес. Периодические движения вблизи сепаратрисы (рис. 19.76) называются кноидальными волнами. Сепаратрисе соответствует локализованное в пространстве решение в виде одиночного возвышения или уединенной волны - солитона (рис. 19.7 в) с амплитудой Это решение аналитически записывается в виде
где - характерная ширина солитона. Справедливость решения легко проверить прямой подстановкой его в уравнение (19.15) при
Рис. 19.6. Потенциальная энергия и фазовый портрет стационарных волн. Состояние равновесия центр. Солитон соответствует сепаратрисе
Рис. 19.7. Различные классы решений уравнения Кортевега-де Вриза и их соответствие фазового портрету стационарных волн: а - квазисинусоидальные колебания малой амплитуды - вблизи состояния центра; - кноидальные волны (периодические солитонные решетки) - вблизи сепаратрисы; в - солитон (уединенная волна) - сепаратриса
Используя при подстановке тождество получаем
Отсюда можно найти . Тождество (19.16) выполняется при любых , следовательно, коэффициенты при одинаковых степенях должны быть равны, т. е.
Итак, мы получили: - чем выше солитон, тем он уже; - чем солитон шире, тем он медленнее бежит и тем меньше его амплитуда. Таким образом, ширина, скорость и амплитуда солитона, описываемого уравнением Кортевега-де Вриза, однозначно связаны, т. е. семейство решений в виде солитонов однопараметрическое - меняем, например, V, получаем разные солитоны.
Почему солитоны, т. е. частные виды стационарных волн, интересны? Фактически по тон же причине, что и другие стационарные волны:
нестационарные возмущения довольно широкого класса в процессе распространения асимптотически приближаются к солитону! Экспериментально этот факт был обнаружен давно; еще более ста лет назад Скотт-Рассел наблюдал солитон и поэтично описал его .
Новая жизнь солитона - одного из самых привлекательных объектов современной физики - в значительной степени связала с построением точных решении многих уравнений нелинейной теории волн. При их построении большую роль сыграл так называемый метод обратной задачи рассеяния . Этот метод берет начало от работы Гарднера, Грина, Крускала и Миуры , которые в 1967 г. установили связь между уравнениями Кортевега-де Вриза и Шредингера. Поясним кратко суть этой связи. Как известно , уравнение Шредингера в случае, когда потенциал положительно определен и спадает до пуля при имеет финитные решения, стремящиеся вместе со своими производными к нулю на бесконечности, а спектр собственных значений дискретен. Рассмотрим уравнение Шредингера
где зависит от времени как от параметра. Тогда и собственные значения, вообще говоря, будут зависеть от Покажем, что собственные значения не будут зависеть от если функция удовлетворяет уравнению Кортевега-де Вриза (точнее, если - любое положительно определенное решение уравнения Кортевега-де Вриза, спадающее на , то соответствующий ему спектр собственных значений остается неизменным). Из уравнения (19.17) находим
Подставим это выражение в уравнение (19.14). После вычислений получим
где штрихи означают соответствующие производные по х.
Проинтегрируем левую и правую части (19.18) по х от до При этом правая часть получившегося уравнения обратится в нуль,
поскольку собственные функции (вместе со своими производными) дискретного спектра уравнения Шредингера исчезают на бесконечности. Таким образом,
Поскольку в силу нормировки то Так как решение произвольно, спектр нам неизвестен. Покажем теперь, что если - солитон, то уравнение Шредингера имеет единственное собственное значение. Когда - солитон, уравнение (19.17) принимает вид
Здесь Дискретные собственные значения уравнения Шредингера даются формулой (см. , § 23, задача 4)
где причем должно быть Подставляя в выражение для выписанные выше значения и а, получим т. е. существует единственное собственное значение Итак, мы получили, что: а) спектр собственных значений не зависит от хотя изменяется со временем; б) каждому собственному значению соответствует солитон. Отсюда следует вывод: любое локализованное положительное возмущение представляет собой набор солитонов и, если достаточно долго подождать, эти солитоны сформируются и возмущение превратится в последовательность солитонов, выстроившихся по амплитуде (рис. 19.8 в). Поскольку «соли-тонный состав» - набор солитонов, из которых состоит возмущение - не зависит от времени, солитоны могут лишь меняться местами в пространстве. Число солитонов зависит от формы начального возмущения; вершины их лежат на одной прямой, так как расстояние, пройденное каждым солитоном, пропорционально его скорости, а последняя, как мы уже знаем, пропорциональна амплитуде.
Такой метод решения уравнения Кортевега-де Вриза называется методом обратной задачи рассеяния, поскольку мы решаем задачу на собственные значения для уравнения Шредингера с потенциалом где играет роль параметра. В квантовомеханическом Если падающая из бесконечности волна плоская с единичной амплитудой, то амплитуда отраженной волны называется коэффициентом отражения. Мы искали сам потенциал. Это и есть решение обратной задачи квантовой теории рассеяния: по известному При дисперсионные эффекты несущественны: основную роль играет нелинейность, приводящая к формированию коротких импульсов, и лишь потом сказывается дисперсия, уравновешивающая процесс (рис. 19.86). Именно так начальное возмущение большей амплитуды распадается на последовательность солитонов, вершины которых лежат на одной прямой (на рис. 19.8 в приведены результаты численных расчетов, взятые из работы ).
Морякам давно известны одиночные волны большой высоты, которые губят корабли. Долгое время считалось, что подобное встречается только в открытом океане. Однако последние данные говорят о том, что одиночные волны-убийцы (до 20-30 метров высотой), или солитоны (от английского solitary - «уединенный»), могут появляться и в прибрежных зонах. Происшествие с «Бирмингемом" Мы находились примерно в 100 милях к юго-западу от Дурбана на пути в Кейптаун. Крейсер шел быстро и почти без качки, встречая умеренную зыбь и ветровые волны, когда внезапно мы провалились в яму и понеслись вниз навстречу следующей волне, которая прокатилась через первые орудийные башни и обрушилась на наш открытый капитанский мостик. Я был сбит с ног и на высоте 10 метров над уровнем моря оказался в полуметровом слое воды. Корабль испытал такой удар, что многие решили, что нас торпедировали. Капитан сразу же уменьшил ход, но эта предосторожность оказалась напрасной, так как умеренные условия плавания восстановились и больше «ям» не попадалось. Это происшествие, случившееся ночью с затемненным кораблем. было одним из наиболее волнующих в море. Я охотно верю, что груженое судно при таких обстоятельствах может потонуть». Так описывает неожиданную встречу с одиночной катастрофической волной британский офицер с крейсера "Бирмингем-. Эта история произошла во время Второй мировой войны, поэтому понятна реакция экипажа, решившего, что крейсер торпедирован. Не столь удачно закончилось аналогичное происшествие с пароходом "Уарита" в 1909 году. На нем находились 211 пассажиров и команда. Погибли все. Такие одиночные неожиданно появляющиеся в океане волны, собственно, и получили название волн-убийц, или солитонов. Казалось бы. любой шторм можно назвать -убийцей.. Ведь действительно, сколько судов погибло во время бури и гибнет сейчас? Сколько моряков нашли свое последнее пристанище в пучинах бушующего моря? И все же волны. возникающие в результате морских штормов и даже ураганов, "убийцами" не называют. Считается, что встреча с солитоном наиболее вероятна у южного побережья Африки. Когда транспортные морские пути благодаря Суэцкому каналу изменились и суда перестали ходить вокруг Африки, количество встреч с волнами-убийцами уменьшилось. Тем не менее уже после Второй мировой войны с 1947 года примерно за 12 лет с солитонами повстречались весьма крупные корабли - "Босфонтейн". "Гиастеркерк", "Оринфонтейн" и "Яхерефонтейн", не считая более мелких местных судов. В период арабо-израильской войны Суэцкий канал был практически закрыт, и движение судов вокруг Африки снова стало интенсивным. От встречи с волной-убийцей в июне 1968 года погиб супертанкер «Уорлд Глори» водоизмещением более 28 тысяч тонн. Танкер получил штормовое предупреждение, и при подходе шторма все выполнялось по инструкции. Ничего плохого не предвиделось. Но среди обычных ветровых волн, которые серьезной опасности не представляли. неожиданно возникла огромная волна высотой около 20 метров с очень крутым фронтом. Она подняла танкер так, что его середина -опиралась» на волну, а носовая и кормовая части оказались в воздухе. Танкер был нагружен сырой нефтью и под своим весом разломился пополам. Эти половинки еще какое-то время сохраняли плавучесть, но через четыре часа танкер ушел на дно. Правда, большую часть экипажа удалось спасти. В 70-е годы «нападения» волн-убийц на корабли продолжались. В августе 1973 года судно "Нептун Сапфир", шедшее из Европы в Японию, в 15 милях от мыса Хермис при ветре около 20 метров в секунду испытало неожиданный удар неизвестно откуда взявшейся одиночной волны. Удар был такой силы, что носовая часть судна длиной примерно 60 метров отломилась от корпуса! Судно «Нептун Сапфир» имело самую совершенную конструкцию для тех лет. Тем не менее встреча с волной-убийцей оказалась для него роковой. Подобных случаев описано довольно много. В страшный перечень катастроф, естественно, попадают не только крупные суда, на которых существуют возможности спасения экипажа. Встреча с волнами-убийцами для малых судов чаще всего заканчивается намного трагичнее. Такие корабли не только испытывают сильнейший удар. способный их разрушить, но на крутом переднем фронте волны могут запросто опрокинуться. Это происходит столь быстро, что рассчитывать на спасение невозможно.Это не цунами Что же это такое - волны-убийцы? Первая мысль, которая приходит в голову осведомленному читателю, - это цунами. После катастрофического «набега» гравитационных волн на юго-восточные берега Азии многие представляют цунами как жуткую стену воды с крутым передним фронтом, обрушивающуюся на берег и смывающую дома и людей. Действительно, цунами способны на многое. После появления этой волны у северных Курил гидрографы, изучая последствия, обнаружили приличных размеров катер, переброшенный через прибрежные холмы в глубь острова. То есть энергия цунами просто поражает. Однако это все касается цунами, «нападающих» на берег. В переводе на русский язык термин "цунами" означает "большая волна в гавани". Ее очень трудно обнаружить в открытом океане. Там высота этой волны обычно не превышает одного метра, а средние, типичные размеры -десятки сантиметров. Да и уклон чрезвычайно маленький, ведь при такой высоте ее длина составляет несколько километров. Так что выявить цунами на фоне бегущих ветровых волн или зыби практически нереально. Почему же при «нападении» на берег цунами становятся такими устрашающими? Дело в том, что эта волна из-за своей большой длины приводит в движение воду по всей глубине океана. И, когда при распространении она достигает сравнительно мелководных районов, вся эта колоссальная масса воды из глубин поднимается вверх. Вот так «безобидная» в открытом океане волна становится разрушительной на побережье. Так что волны-убийцы - это не цунами. На самом дел солитоны - это необыкновенное и малоизученное явление. Их называют волнами, хотя на самом деле они нечто иное. Для возникновения солитонов, конечно, необходим некоторый изначальный импульс, удар, иначе откуда взяться энергии, но не только. В отличие от обычных волн солитоны распространяются на большие расстояния с очень малым рассеянием энергии. Это загадка, которая еще ждет изучения. Солитоны практически не взаимодействуют друг с другом. Как правило, они распространяются с разными скоростями. Конечно, может получиться так, что один солитон догонит другой, и тогда они суммируются по высоте, но потом все равно снова разбегаются по своим путям. Конечно, сложение солитонов - редкое событие. Но есть еще одна причина резкого возрастания у них крутизны и высоты. Причина эта - подводные уступы, через которые «пробегает» солитон. При этом в подводной части происходит отражение энергии, и волна как бы «выплескивается» вверх. Подобная ситуация изучалась на физических моделях международной научной группой. Опираясь на эти исследования, можно прокладывать более безопасные маршруты движения судов. Но загадок все же остается намного больше, чем изученных особенностей, и тайна волн-убийц по-прежнему ждет своих исследователей. Особенно загадочны солитоны внутри вод моря, на так называемом «слое скачка плотности». Эти солитоны могут приводить (или уже приводили) к катастрофам подводных лодок.