Пред.
След.
Макеты страниц
Распознанный текст, спецсимволы и формулы могут содержать ошибки, поэтому с корректным вариантом рекомендуем ознакомиться на отсканированных изображениях учебника выше Также, советуем воспользоваться поиском по сайту, мы уверены, что вы сможете найти больше информации по нужной Вам тематике ДЛЯ СТУДЕНТОВ И ШКОЛЬНИКОВ ЕСТЬ
ZADANIA.TO
Знаменитый физик-ядерщик, лауреат Нобелевской премии по физике Энрико Ферми в середине 50 -х годов отошел от дел, связанных с американским атомным проектом. Продолжая работать в ядерном центре в Лос-Аламосе, он получил доступ к первым ЭВМ, только что разработанным и установленным там для сложных физических расчетов. Совместно с двумя своими сотрудниками, Дж. Паста и С. Уламом, он занялся, как в то время казалось, совершенно изолированной проблемой — исследованием поведения систем обыкновенных дифференциальных уравнений, которые первоначально были линейными, но в которые была привнесена нелинейность как возмущение. Если бы такого возмущения не было, то энергия каждой нормальной моды линейной системы, то есть колебаний с заданной частотой, была бы постоянной. Можно было надеяться, что нелинейные взаимодействия между модами приведут к тому, что энергия системы равномерно распределится между всеми модами — этот результат был бы в согласии с теоремой о равномерном распределении, которая предсказывает поведение нелинейных систем общего вида. Но полученные результаты противоречили этой идее. Рассмотрим динамическую систему, состоящую из $N$ идентичных частиц единичной массы на прямой с фиксированными конечными точками ( $N$ велико; ФПУ рассматривали $N=64$ ). Предположим, что между ближайшими соседями действуют упругие силы, как если бы они были связаны Для описания этого движения обозначим через $Q_{n}$ смещение из состояния равновесия $n$-й частицы. Тогда по второму закону Ньютона уравнение движения этой частицы можно записать Здесь $f(Q)$ — упругая сила взаимодействия между ближайшими соседями, которая включает в себя обычный линейный закон Гука для малых колебаний пружинки плюс некоторый малый нелинейный член. ФПУ рассматривали два случая: Здесь $\gamma$ — линейная константа цепочки, а константы $\alpha$ и $\beta$ выбираются таким образом, чтобы максимальное смещение $Q$, вызываемое нелинейным членом, было мало. Если использовать нелинейности такого вида и численно решить уравнение (2.1) с начальными условиями в виде, скажем, синусоидальной волны, то обнаружится, что энергия не распределяется между всеми нормальными модами, а остается в начальной моде и нескольких ближайших. К тому же плотность энергии в этих ближайших модах имеет почти периодический характер зависимости от времени. Рисунок взят из оригинальной статьи, написанной ФПУ, и представляет собой график зависимости энергии от времени для первых пяти мод с синусоидальными начальными данными. После большого числа колебаний энергия каждой нормальной моды выглядит как почти периодическая функция от времени, причем энергия не перераспределяется в более высокие моды с течением времени. Строгое объяснение этой периодичности стимулировало более глубокое изучение уравнений типа (2.1). При некоторых приближениях уравнение (2.1) можно преобразовать в уравнение КдФ в континуальном пределе, то есть при переходе к непрерывному случаю. Введем шаг решетки $l$ и определим переменную $x=n l$ как расстояние вдоль решетки. Для того чтобы заменить члены типа $Q_{n}$ непрерывными переменными, используем разложение в ряд Тейлора, полагая Рис. 2. Распределение энергии единичной частицы по Фурье-модам с течением времени Имеем Теперь уравнение (2.1) приобретает вид Если функция $f$ выбрана в виде первого выражения (2.2), то получаем $f^{\prime}(0)=\gamma, f^{\prime \prime}(0)=2 \alpha, f^{\prime \prime \prime}(0)=0$. Поэтому если положить $r=1$, то члены с квадратичной нелинейностью оказываются одного порядка с членом, содержащим четвертую производную. Полагая $u(x, t)=p_{x}(x, t)$, с точностью до $O\left(l^{4}\right)$ имеем Уравнение (2.5) известно в гидродинамике под названием уравнения Буссинеска. Если рассматривать это уравнение с точностью до $O\left(l^{2}\right)$, то оно превращается в линейное волновое уравнение и влияние нелинейности пропадает. Учет следующего приближения порядка $O\left(l^{4}\right)$ описывает уединенные волны, которые могут перемещаться как влево, так и вправо. Это видно из явной формулы для одиночного солитона ( на бесконечной прямой по $x$ ) где $\omega^{2}=\gamma a^{2} l^{2}+a^{4} l^{4} / 12$, константа $a$ произвольна. демонстрируют экспоненциально быстрое убывание коэффициентов Фурье Отсюда следует, что таким же свойством обладают коэффициенты Фурье для кинетической энергии $\left|Q_{n}^{\prime}(t)\right|^{2}$ единичной частицы, что находится в согласии с экспериментальными данными, изображенными на рисунке. При этом следует заметить, что солитоны вида (2.6) являются устойчивыми решениями, то есть мало изменяются при малом возмущении уравнения (2.5) и/или начальных условий. Этот факт будет доказан ниже для уравнения КдФ (см. Лекцию 12), хотя он справедлив для всех солитонных уравнений. Последнее означает, что выводы, сделанные относительно решения Покажем теперь, как с помощью еще одной автомодельной редукции уравнение (2.5) может быть сведено к уравнению КдФ. Для этого рассмотрим волны, распространяющиеся только в одном направлении вдоль цепочки. Используем масштабные прообразования переменных $x, t$ и $u$. Сначала введем малый параметр $\varepsilon$ и новые пространственную и временную переменные затем разложим $u$ по степеням параметра $\varepsilon$ : Для того чтобы все члены с четвертой производной, производные по времени, а также нелинейные члены были одного и того же порядка по $\varepsilon$, необходимо положить $c=\sqrt{\gamma} l, s=1 / 2$ и $q=3 / 2$. Полученное при этом уравнение для $U_{1}$ можно проинтегрировать один раз по $\xi$, при этом получится уравнение КдФ ЗАМЕЧАНИЕ. В некоторых случаях нелинейные цепочки вида (2.1) без всяких асимптотических приближений демонстрируют солитонное поведение. Таковой является, например, так называемая цепочка Toды Для кубически нелинейной цепочки с функцией $f(Q)$, определенной по второй формуле (2.2), имеем $f(0)=\gamma, f^{\prime \prime}(0)=0$ и $f^{\prime \prime \prime}(0)=6 \beta$. Теперь для того чтобы кубически нелинейные члены и члены с четвертой производной оказались одного порядка, придется выбрать $r=0$. Тогда с точностью до $O\left(l^{4}\right)$ получим Проделав такие же масштабные преобразования переменных, как и в предыдущем случае, мы сможем теперь привести интересующие нас члены к где $\omega=\frac{\sqrt{a^{3} l^{3}}}{24 \sqrt{\gamma}}$ и $a-$ произвольная постоянная.
|
1 |
Оглавление
|