Журнал вычислительной математики и математической физики, 2020, T. 60, № 2, стр. 253-266
Асимптотики регулярных решений в задаче Камасса–Холма
С. А. Кащенко *
Ярославский гос. университет им. П.Г. Демидова
Национальный исследовательский ядерный университет “МИФИ”
150000 Ярославль, ул. Советская, 14, Россия
* E-mail: kasch@uniyar.ac.ru
Поступила в редакцию 16.07.2019
После доработки 16.07.2019
Принята к публикации 17.10.2019
Аннотация
Рассматривается периодическая краевая задача для модернизированного уравнения Камасса–Холма, которое отличается от известного классического уравнения наличием еще нескольких квадратичных слагаемых. Сформулированы три важных условия на коэффициенты уравнения, выполнение которых относит исходное уравнение к уравнениям типа Камасса–Холма. Исследуется вопрос о динамических свойствах так называемых регулярных решений в окрестностях всех состояний равновесия. Для определения “главных” составляющих решений строятся специальные нелинейные краевые задачи. Получены асимптотические формулы для множества периодических решений и конечномерных торов. Изучается вопрос о бесконечномерных торах. Показано, что компактная запись в виде уравнения в частных производных нормализованного уравнения в задаче о таких торах возможна лишь для классического уравнения Камасса–Холма. Приведен асимптотический анализ в случаях, когда один из коэффициентов линейной части уравнения является достаточно малым и когда значение периода в граничных условиях является достаточно большим. Библ. 16.
1. ВВЕДЕНИЕ
1. Описание краевой задачи. При $x \in ( - \infty ,\infty )$, $t \geqslant {{t}_{0}}$ рассматриваем уравнение
(1)
$\frac{{\partial {v}}}{{\partial t}} = {{a}_{1}}\frac{{{{\partial }^{3}}{v}}}{{\partial t\partial {{x}^{2}}}} + {{a}_{2}}\frac{{\partial {v}}}{{\partial x}} + {{a}_{3}}{v}\frac{{\partial {v}}}{{\partial x}} + {{a}_{4}}\frac{{\partial {v}}}{{\partial x}}\frac{{{{\partial }^{2}}{v}}}{{\partial {{x}^{2}}}} + {{a}_{5}}{v}\frac{{{{\partial }^{3}}{v}}}{{\partial {{x}^{3}}}}.$Уравнение (1) исследовалось многими авторами (см. [1]–[9]). Основные результаты посвящены вопросам интегрируемости и нахождению решений.
Вместе с уравнением (1) будем рассматривать $L$-периодические краевые условия
Сначала сформулируем три важных вывода о решениях краевой задачи (1), (3).
Вывод 1. Каждое стационарное значение функции ${v}$: ${v}(t,x) \equiv d$ $(d \in ( - \infty ,\infty ))$ является состоянием равновесия краевой задачи (1), (3).
Положим в (1)
и будем рассматривать все такие функции ${v}(t,x)$, которые при $t \geqslant {{t}_{0}}$, $x \in ( - \infty ,\infty )$ удовлевoряют условию (3) и соотношению(5)
$M({v}) \equiv d,\quad {\text{где}}\quad M(\varphi (x)) = \frac{1}{L}\int\limits_0^L {\varphi (x)dx} .$Вывод 2. Из условия $M({v}({{t}_{0}},x)) = d$ следует, что $M({v}(t,x)) \equiv d$ для всех таких значений $t > {{t}_{0}}$, для которых решение ${v}(t,x)$ существует.
Произведем в (1), (3), (5), замену (4) и перейдем к “движущейся” переменной
В результате получаем уравнение
(7)
$\frac{{\partial u}}{{\partial t}} = \alpha \frac{{{{\partial }^{3}}u}}{{\partial t\partial {{x}^{2}}}} + \kappa \frac{{\partial u}}{{\partial x}} + au\frac{{\partial u}}{{\partial x}} + b\frac{{\partial u}}{{\partial x}}\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}} + cu\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}}$(9)
$\frac{{\partial u}}{{\partial t}} = \alpha \frac{{{{\partial }^{3}}u}}{{\partial t\partial {{x}^{2}}}} + \delta \frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}} + au\frac{{\partial u}}{{\partial x}} + b\frac{{\partial u}}{{\partial x}}\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}} + cu\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}},$Поставим задачу исследования локального – в достаточно малой окрестности нулевого состояния равновесия – поведения всех решений краевой задачи (7), (8) при каждом $d \in ( - \infty ,\infty )$. Основу асимптотического исследования составляет построение уравнений для главных членов асимптотических представлений решений. Такое представление часто называют нормализацией.
Важную роль играет характер расположения корней ${{\lambda }_{k}}$ характеристического уравнения для линеаризованной в нуле краевой задачи (7), (8): ${{\lambda }_{k}} = i\lambda ({{z}_{k}}),$ где
Вывод 3. При всех значениях параметра $d$ все корни характеристического уравнения являются чисто мнимыми, а значит, в задаче об устойчивости нулевого решения краевой задачи (7), (8) при всех $d \in ( - \infty ,\infty )$ реализуется критический случай бесконечной размерности.
Обозначим через $F(u)$ квадратичное выражение, содержащие величины $u,\;\tfrac{{\partial u}}{{\partial x}},\;\tfrac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}},\;u\tfrac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}}$:
Рассмотрим более общую по сравнению с (7) краевую задачу
(10)
$\frac{{\partial u}}{{\partial t}} = \alpha \frac{{{{\partial }^{3}}u}}{{\partial t\partial {{x}^{2}}}} + \kappa \frac{{\partial u}}{{\partial x}} + F(u) - M(F(u)),$Согласно выводу 3 из условия $\operatorname{Re} {{\lambda }_{k}} = 0$, $k = \pm 1, \pm 2,...$, вытекает, что ${{f}_{1}} = {{f}_{3}} = 0$.
В ряде задач, например, в задачах об асимптотике периодических по времени решений полезно рассмотреть уравнение, которое отличается от (10) наличием линейных членов со второй и четвертой производными по пространственной переменной и с малым множителем $\mu $:
(12)
$\mu \left( {{{b}_{1}}\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}} + {{b}_{2}}\frac{{{{\partial }^{4}}u}}{{\partial {{x}^{4}}}}} \right),\quad {{b}_{2}} < 0,\quad 0 < \mu \ll 1.$Итак, после замены (4) и $x \to x + d{{f}_{4}}{{\alpha }^{{ - 1}}}t$ и с учетом выражения (12) приходим к краевой задаче
(13)
$\begin{gathered} \frac{{\partial u}}{{\partial t}} = \alpha \frac{{{{\partial }^{3}}u}}{{\partial t\partial {{x}^{2}}}} + {{\kappa }_{0}}\frac{{\partial u}}{{\partial x}} + {{g}_{1}}u\frac{{\partial u}}{{\partial x}} + {{g}_{2}}u\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}} + {{g}_{3}}\left( {\mathop {\left( {\frac{{\partial u}}{{\partial x}}} \right)}\nolimits^2 - M\left( {\mathop {\left( {\frac{{\partial u}}{{\partial x}}} \right)}\nolimits^2 } \right)} \right) + \\ + \;{{g}_{4}}\frac{{\partial u}}{{\partial x}}\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}} + {{g}_{5}}\frac{{\partial u}}{{\partial x}}\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}} + {{g}_{6}}\left( {\mathop {\left( {\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}}} \right)}\nolimits^2 - M\left( {\mathop {\left( {\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}}} \right)}\nolimits^2 } \right)ght) + {{g}_{7}}\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}}\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}} + } \right. \\ \, + {{g}_{8}}\left( {\mathop {\left( {\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}}} \right)}\nolimits^2 - M\left( {\mathop {\left( {\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}}} \right)}\nolimits^2 } \right)} \right) + \mu \left( {{{b}_{1}}\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}} + {{b}_{2}}\frac{{{{\partial }^{4}}u}}{{\partial {{x}^{4}}}}} \right). \\ \end{gathered} $Для этой краевой задачи можно предложить другую – эквивалентную запись. Сначала введем одно обозначение. Пусть $w(x)$ – есть ($2\pi {\text{/}}L$)-периодические функции с нулевым средним:
Продифференцируем (13) и (14) по $x$ и положим $w = \tfrac{{\partial u}}{{\partial x}}$, т.е. $u = J(w)$. Тогда приходим к уравнению, не содержащему оператор усреднения:
(16)
$\begin{gathered} \frac{{\partial w}}{{\partial t}} = \alpha \frac{{{{\partial }^{3}}w}}{{\partial t\partial {{x}^{2}}}} + {{\kappa }_{0}}\frac{{\partial w}}{{\partial x}} + {{g}_{1}}\left[ {{{w}^{2}} + J(w)\frac{{\partial w}}{{\partial x}}} \right] + {{g}_{2}}\left[ {w\frac{{{{\partial }^{2}}w}}{{\partial {{x}^{2}}}} + J(w)\frac{{{{\partial }^{3}}w}}{{\partial {{x}^{3}}}}} \right] + \\ \, + 2{{g}_{3}}w\frac{{\partial w}}{{\partial x}} + {{g}_{4}}\left[ {\mathop {\left( {\frac{{\partial w}}{{\partial x}}} \right)}\nolimits^2 + w\frac{{{{\partial }^{2}}w}}{{\partial {{x}^{2}}}}} \right] + \left[ {\frac{{\partial w}}{{\partial x}}\frac{{{{\partial }^{2}}w}}{{\partial {{x}^{2}}}} + w\frac{{{{\partial }^{3}}w}}{{\partial {{x}^{3}}}}} \right] + \\ \, + 2({{g}_{6}} + {{g}_{8}})\frac{{{{\partial }^{2}}w}}{{\partial {{x}^{2}}}}\frac{{{{\partial }^{3}}w}}{{\partial {{x}^{3}}}} + {{g}_{7}}\frac{{\partial w}}{{\partial x}}\frac{{{{\partial }^{3}}w}}{{\partial {{x}^{3}}}} + \kappa \left[ {{{b}_{1}}\frac{{{{\partial }^{2}}w}}{{\partial {{x}^{2}}}} + {{b}_{2}}\frac{{{{\partial }^{4}}w}}{{\partial {{x}^{4}}}}} \right]. \\ \end{gathered} $Важно отметить, что краевые задачи (13), (14) и (16), (17) зависят от параметра $d \in ( - \infty ,\infty )$ и от “внутреннего” параметра $L \in (0,\infty )$.
Ниже вводится еще один малый параметр $\varepsilon $, который учитывает варьирования значения периода $L$. Предполагаем, что значение ${{L}_{0}} > 0$ как-то фиксировано и для некоторого ${{L}_{1}}$ выполнено равенство
В работе рассматриваются так называемые регулярные решения (см. [12]–[14]). Регулярность означает, что эти решения и их производные являются ограниченными по $x \in [0,L]$ функциями при каждом фиксированном $t$. В частности, верно асимптотическое представление
Отметим, что пространственные производные нерегулярных решений могут достигать больших значений при стремлении к нулю малого параметра. Результаты исследования нерегулярных решений существенно отличаются от приводимых ниже. Им будет посвящена отдельная публикация.
Поставим задачу исследования поведения всех решений этих краевых задач из достаточно малой окрестности нулевого состояния равновесия. При условии (18) в разд. 2 изучается асимптотика периодических по $t$ решений и одномерных торов в рассматриваемых краевых задачах. Соответствующие построения можно применить и к рассмотрению асимптотики торов большей размерности. Для бесконечномерных торов методика оказывается неэффективной, поскольку приводит к бесконечной системе нелинейных обыкновенных дифференциальных уравнений, которую, вообще говоря, нельзя упростить.
В связи с этим возникает задача определения таких коэффициентов краевой задачи (13), (14) ((16), (17)), для которых соответствующую бесконечномерную систему удается записать в компактной форме в виде уравнения в частных производных. Решению этой задачи посвящен разд. 3.
Наконец, отметим, что без ограничения общности можно считать, что выполнено условие
Это равенство (при $\alpha > 0$) достигается нормировкой пространственной переменной и очевидными переобозначениями параметров краевой задачи. Возникает, конечно, необходимость асимптотического изучения рассматриваемых краевых задач отдельно при условиях и Кроме этого, весьма важно изучение решений в предположении, что значение параметра $L$ в краевых условиях достаточно велико:Этому случаю и случаю (20) посвящен разд. 4. Случай (21) будет изучен в отдельной публикации, посвященной нерегулярным решениям.
2. ПОСТРОЕНИЕ АСИМПТОТИКИ ПЕРИОДИЧЕСКИХ РЕШЕНИЙ И ТОРОВ
Во всех приведенных выше краевых задачах фигурирует “внутренний” параметр $L$, входящий в периодические краевые условия. Будем варьировать этот параметр и изучать зависимость от него асимптотики решений.
Фиксируем произвольно ${{L}_{0}} > 0$ и ${{L}_{1}}$, а через $\varepsilon $ обозначим малый положительный параметр. Положим
Еще один малый параметр $\mu $, фигурирующий в (13) и (16), удобно представить в виде
Остановимся для определенности на рассмотрении при условии (19) краевой задачи (13), (14). Для того чтобы период по пространственной переменной у решений не зависел от $\varepsilon $, в (13), (14) произведем замену
В результате приходим к краевой задаче
(26)
$\begin{gathered} \frac{{\partial u}}{{\partial t}} = \alpha (\varepsilon )\frac{{{{\partial }^{3}}u}}{{\partial t\partial {{x}^{2}}}} + \kappa (\varepsilon )\frac{{\partial u}}{{\partial x}} + {{g}_{1}}(\varepsilon )u\frac{{\partial u}}{{\partial x}} + {{g}_{2}}(\varepsilon )u\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}} + {{g}_{3}}(\varepsilon )\left[ {\mathop {\left( {\frac{{\partial u}}{{\partial x}}} \right)}\nolimits^2 - M\left( {\mathop {\left( {\frac{{\partial u}}{{\partial x}}} \right)}\nolimits^2 } \right)} \right] + \\ \, + {{g}_{4}}(\varepsilon )\frac{{\partial u}}{{\partial x}}\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}} + {{g}_{5}}(\varepsilon )\frac{{\partial u}}{{lx}}\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}} + {{g}_{6}}(\varepsilon )\left[ {\mathop {\left( {\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}}} \right)}\nolimits^2 - M\left( {\mathop {\left( {\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}}} \right)}\nolimits^2 } \right)} \right] + {{g}_{7}}(\varepsilon )\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}}\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}} + \\ \, + {{g}_{8}}\left[ {\mathop {\left( {\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}}} \right)}\nolimits^2 - M\left( {\mathop {\left( {\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}}} \right)}\nolimits^2 } \right)} \right] + \varepsilon \left[ {{{b}_{1}}(\varepsilon )\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}} + {{b}_{2}}(\varepsilon )\frac{{{{\partial }^{4}}u}}{{\partial {{x}^{4}}}}} \right], \\ \end{gathered} $2.1. Построение асимптотики периодических решений
Фиксируем произвольно целое ${{k}_{0}} > 0$ и положим ${{z}_{0}} = 2\pi {{k}_{0}}L_{0}^{{ - 1}}$. Отметим, что выполнены условия
(28)
$\lambda (2{{z}_{0}}) \ne z\lambda ({{z}_{0}})\quad {\text{и}}\quad \lambda (3{{z}_{0}}) \ne 3\lambda ({{z}_{0}}).$Линеаризованная в нуле краевая задача (26), (27) имеет, в частности, периодические решения
Введем в рассмотрение формальный ряд
(29)
$u = {{\varepsilon }^{{1/2}}}[\xi (\tau )exp[i{{z}_{0}}x + i\lambda ({{z}_{0}})t] + \overline \xi (\tau )exp( - i{{z}_{0}}x - i\lambda ({{z}_{0}})t)] + \varepsilon {{u}_{2}}(t,\tau ,x) + {{\varepsilon }^{{3/2}}}{{u}_{3}}(t,\tau ,x) + ...\;.$(30)
${{u}_{2}}(t,\tau ,x) = B{{\xi }^{2}}exp(2i{{z}_{0}}x + 2i\lambda ({{z}_{0}})t) + \overline B \overline \xi exp( - 2i{{z}_{0}}x - 2i\lambda ({{z}_{0}})t),$На следующем шаге, собирая в формальном тождестве коэффициенты при ${{\varepsilon }^{{3/2}}}$, получаем уравнение для ${{u}_{3}}(t,\tau ,x)$. Из условия его разрешимости в указанном классе функций приходим к уравнению для определения $\xi (\tau )$:
Для коэффициентов $\beta $ и $\delta $ имеем равенстваТеорема 1. Пусть ${{\xi }_{0}}(\tau )$ – ограниченное при $\tau \to \infty $ решение уравнения (29). Тогда краевая задача (26), (27) имеет асимптотическое по невязке с точностью до O(${{\varepsilon }^{{3/2}}}$) периодическое решение ${{u}_{0}}(t,\tau ,x)$, для которого
Сделаем ряд выводов. Сначала отметим, что при $\delta = {{g}_{3}} = {{g}_{5}} = {{g}_{6}} = {{g}_{8}} = 0$ коэффициенты в уравнении (31) являются чисто мнимыми. Поэтому все его решения периодические (негрубые). Это дает основание говорить о том, что в малой окрестности нуля краевая задача (26), (27) имеет континуальное семейство асимптотических по невязке решений. Их асимптотика определяется по решениям $\xi (\tau )$ уравнения (31) из формулы (29).
Уравнение (31) интегрируется в явном виде при любых значениях коэффициентов $\beta $ и $\delta $. Например, при условии $\operatorname{Re} \beta \cdot \operatorname{Re} \delta < 0$ это уравнение имеет устойчивый цикл
Обратим внимание, что произвольно фиксированный (континуальный) параметр ${{L}_{1}}$ и параметр ${{k}_{0}}$, принимающий целые значения, являются “внутренними” параметрами: само уравнение (26) от них не зависит. Поэтому формулы (29) дают представление для довольно богатого множества асимптотических по невязке решений в (26), (27).
2.2. Построение 2-мерных торов
Фиксируем произвольно два различных целых положительных значения ${{k}_{1}}$ и ${{k}_{2}}$. Рассмотрим формальный ряд
(32)
$\begin{gathered} u = {{\varepsilon }^{{1/2}}}[\xi (\tau )exp(i{{z}_{1}}x + i{{\lambda }_{1}}t) + \overline \xi (\tau )exp( - i{{z}_{1}}x - i{{\lambda }_{1}}) + \eta (\tau )exp(i{{z}_{2}}x + i{{\lambda }_{2}}t) + \\ \, + \bar {\eta }(\tau )exp( - i{{z}_{2}}x - i{{\lambda }_{2}}t)] + \varepsilon {{u}_{2}}(t,\tau ,x) + {{\varepsilon }^{{3/2}}}{{u}_{3}}(t,\tau ,x) + \ldots ,\quad \tau = \varepsilon t, \\ \end{gathered} $(33)
$\begin{gathered} \dot {\xi } = {{{\tilde {B}}}_{1}}\xi + \xi [{{D}_{1}}{{\xi }^{2}} + {{D}_{2}}{{\eta }^{2}}], \\ \dot {\eta } = {{{\tilde {B}}}_{2}}\eta + \eta [{{D}_{3}}{{\xi }^{2}} + {{D}_{4}}{{\eta }^{2}}]. \\ \end{gathered} $2.3. О более общих случаях
По аналогии с (32) можно рассмотреть вопрос об асимптотике 3, 4-мерных торов. В общем случае формальные выражения имеют вид
(34)
$u = {{\varepsilon }^{{1/2}}}\sum\limits_{k \ne 0\;k = - \infty }^\infty \,{{\xi }_{k}}(\tau )exp(i{{z}_{k}}x + i{{\lambda }_{k}}t) + \varepsilon {{u}_{2}}(t,\tau ,x) + {{\varepsilon }^{{3/2}}}{{u}_{3}}(t,\tau ,x) + \ldots ({{\xi }_{{ - k}}} = {{\bar {\xi }}_{k}}(\tau )),$3. АСИМПТОТИКА БЕСКОНЕЧНОМЕРНЫХ ТОРОВ
Остановимся подробнее на схеме построения асимптотики решений этой краевой задачи, базирующейся на формальном ряде (34).
3.1. Построение функций ${{u}_{2}}(t,\tau ,x)$
Ниже для простоты будем предполагать, что $\delta = 0$, значение $L$ не зависит от малого параметра $\varepsilon $, а варьируется только параметр $\kappa :\kappa = {{\kappa }_{0}} + \varepsilon {{\kappa }_{1}}$, $0 < \varepsilon \ll 1$. Таким образом, рассматривается краевая задача
(35)
$\begin{gathered} \frac{{\partial u}}{{\partial t}} = \alpha \frac{{{{\partial }^{3}}u}}{{\partial t\partial {{x}^{2}}}} + ({{\kappa }_{0}} + \varepsilon {{\kappa }_{1}})\frac{{\partial u}}{{\partial x}} + {{g}_{1}}u\frac{{\partial u}}{{\partial x}} + {{g}_{2}}u\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}} + {{g}_{3}}\left( {\mathop {\left( {\frac{{\partial u}}{{\partial x}}} \right)}\nolimits^2 - M\left( {\mathop {\left( {\frac{{\partial u}}{{\partial x}}} \right)}\nolimits^2 } \right)} \right) + {{g}_{4}}\frac{{\partial u}}{{\partial x}}\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}} + \\ \, + {{g}_{5}}\frac{{\partial u}}{{\partial x}}\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}} + {{g}_{6}}\left( {\mathop {\left( {\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}}} \right)}\nolimits^2 - M\left( {\mathop {\left( {\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}}} \right)}\nolimits^2 } \right)} \right) + {{g}_{7}}\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}}\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}} + {{g}_{8}}\left( {\mathop {\left( {\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}}} \right)}\nolimits^2 - M\left( {\mathop {\left( {\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}}} \right)}\nolimits^2 } \right)} \right), \\ \end{gathered} $Положим здесь
Отметим, для примера, что запись(38)
$\sum\limits_{k = - \infty ,k \ne 0}^\infty \,\sum\limits_{p = - \infty \;k + p \ne 0,p \ne 0}^\infty \,{{(k + p)}^{s}}{{k}^{m}}{{p}^{n}}{{\xi }_{k}}{{\xi }_{p}}$(39)
$\frac{{{{\partial }^{s}}}}{{\partial {{x}^{s}}}}\left( {\left( {\frac{{{{\partial }^{m}}\xi }}{{\partial {{x}^{m}}}}} \right)\left( {\frac{{{{\partial }^{n}}\xi }}{{\partial {{x}^{n}}}}} \right)} \right)$Подставим (34) в (26). Тогда для ${{u}_{2}}(t,\tau ,x)$ приходим к уравнению
(40)
$\frac{{\partial {{u}_{2}}}}{{\partial t}} = \frac{{{{\partial }^{3}}{{u}_{2}}}}{{\partial t\partial {{x}^{2}}}} + {{\kappa }_{0}}\frac{{\partial {{u}_{2}}}}{{\partial x}} + \sum\limits_{k = - \infty }^\infty \,\sum\limits_{p = - \infty }^\infty \,R({{z}_{k}},{{z}_{p}}){{\xi }_{k}}(\tau ){{\xi }_{p}}(\tau )exp[i({{z}_{k}} + {{z}_{p}})x + i(\lambda ({{z}_{k}}) + \lambda ({{z}_{p}}))t],$Функцию ${{u}_{2}}(t,\tau ,x)$ ищем в виде
(41)
${{u}_{2}} = \sum\limits_{k = - \infty }^\infty \,\sum\limits_{p = - \infty }^\infty \,{{u}_{{kp}}}{{\xi }_{k}}{{\xi }_{p}}exp[i({{z}_{k}} + {{z}_{p}})x + i(\lambda ({{z}_{k}}) + \lambda ({{z}_{p}}))].$(42)
$N(k,p) = i\kappa (k + p)kp[3 + ({{k}^{2}} + {{p}^{2}} + kp)]{{[(1 + {{k}^{2}})(1 + {{p}^{2}})]}^{{ - 1}}}.$Поставим задачу определения всех таких коэффициентов ${{g}_{i}}$, при которых функцию ${{u}_{2}}(t,\tau ,x)$ из (41) можно записать, не используя бесконечное суммирование, а используя только производные и интегралы от $\xi (\tau ,x)$ и от их квадратичных произведений. Ниже удобно нормировать функцию $u(t,x)$ так, чтобы выполнялось равенство
а коэффициент ${{g}_{6}}$ переобозначить через $g:{{g}_{6}} = g$.Сформулируем основной результат
Теорема 2. Функцию ${{u}_{2}}(t,x)$ можно выразить через $\xi (\tau ,x)$, не используя бесконечные суммирования, в том и только в том случае, когда выполнены равенства
Доказательство. Положим в (40) $kp = x$ и $(k + p) = y$. Выражение $N(k,p)$ можно тогда записать в виде
а выражение $R(k,p)$ имеет вид $P(x,y)$, гдеСогласно (36) и (37), условие представления ${{u}_{2}}$ через функции $\xi $ и через производные и “интегралы” от нее состоит в том, чтобы выражение $P(x,y)$ было представимо в виде
Тем самым при $x = 3 + {{y}^{2}}$ функция ${{P}_{0}}(y) = P(3 + {{y}^{2}},y)$ должна тождественно обращаться в нуль. Тем самым из равенства нулю тождественно многочлена 6-й степени относительно $y$В порядке обобщения теоремы 2 заметим, что при условии $g = 0$ на основании равенств (43)–(46) получаем классическое уравнение Камасса–Холма, а если бы вместо (43) принять условие ${{g}_{2}} = 0$, то имеем уравнение
(47)
$\frac{{\partial u}}{{\partial t}} = \frac{{{{\partial }^{3}}u}}{{\partial t\partial {{x}^{2}}}} + \kappa \frac{{\partial u}}{{\partial x}} - 3g\left[ {\mathop {\left( {\frac{{\partial u}}{{\partial x}}} \right)}\nolimits^2 - M\left( {\mathop {\left( {\frac{{\partial u}}{{\partial x}}} \right)}\nolimits^2 } \right)} \right] + 2g\frac{{\partial u}}{{\partial x}}\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}} + g\left[ {\mathop {\left( {\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}}} \right)}\nolimits^2 - M\left( {\mathop {\left( {\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}}} \right)}\nolimits^2 } \right)} \right].$В общем случае из теоремы 2 получаем уравнение (без учета слагаемого с множителем $\kappa $)
(48)
$\begin{gathered} \frac{{\partial u}}{{\partial t}} = \frac{{{{\partial }^{3}}u}}{{\partial t\partial {{x}^{2}}}} + \kappa \frac{{\partial u}}{{\partial x}} - 3u\frac{{\partial u}}{{\partial x}} + 2u\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}} - 3g\left[ {\mathop {\left( {\frac{{\partial u}}{{\partial x}}} \right)}\nolimits^2 - M\left( {\mathop {\left( {\frac{{\partial u}}{{\partial x}}} \right)}\nolimits^2 } \right)} \right] + \\ \, + 2\frac{{\partial u}}{{\partial x}}\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}} + 2g\frac{{\partial u}}{{\partial x}}\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}} + g\left[ {\mathop {\left( {\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}}} \right)}\nolimits^2 - M\left( {\mathop {\left( {\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}}} \right)}\nolimits^2 } \right)} \right]. \\ \end{gathered} $Приведем явный вид функции ${{u}_{2}}(t,x)$ при условиях (44)–(46). В случае ${{g}_{0}} = 0$, т.е. для уравнения Камасса–Холма, имеет равенство ${{u}_{2}}(t,\tau ,x) = U(t,\tau ,x)$, где
(49)
$U(t,\tau ,x) = - \frac{\kappa }{2}\left[ {\mathop {\left( {J(\xi )} \right)}\nolimits^2 - M(\mathop {(J(\xi ))}\nolimits^2 ) + \mathop {\left( {\frac{{\partial \xi }}{{\partial x}}} \right)}\nolimits^2 - M\left( {\mathop {\left( {\frac{{\partial \xi }}{{\partial x}}} \right)}\nolimits^2 } \right) + 2J(\xi )\frac{{\partial \xi }}{{\partial x}}} \right].$3.2. Построение нормализованного уравнения
На следующем шаге, подставляя ряд (34) в (26) соберем коэффициенты при ${{\varepsilon }^{{3/2}}}$. В результате с учетом (48) проходим к уравнению
(52)
$\frac{{\partial {{u}_{3}}}}{{\partial t}} - \frac{{{{\partial }^{3}}{{u}_{3}}}}{{\partial t\partial {{x}^{2}}}} - \kappa \frac{{\partial {{u}_{3}}}}{{\partial x}} = - \frac{1}{2}\frac{\partial }{{\partial x}}\left[ {3{{u}_{2}}\frac{{\partial \xi }}{{\partial x}} + 16\frac{{\partial \xi }}{{\partial x}}\frac{{\partial {{u}_{2}}}}{{lx}} - 2\frac{{{{\partial }^{2}}}}{{\partial {{x}^{2}}}}(\xi {{u}_{2}})} \right].$Критерий разрешимости (52) относительно ${{u}_{3}}(t,\tau ,x)$ в указанном классе функций состоит в отсутствии в правой части гармоник вида $exp(ikx + i{{\lambda }_{k}}t)$. Отсюда, учитывая формулы (49) и (37), для величин ${{\xi }_{k}}(\tau )$ приходим к бесконечной системе обыкновенных дифференциальных уравнений
(53)
$\frac{{d{{\xi }_{k}}}}{{d\tau }} = {{\Phi }_{k}}( \ldots ,{{\xi }_{{ - 2}}},{{\xi }_{{ - 1}}},{{\xi }_{1}},{{\xi }_{2}}, \ldots ),\quad l = \pm 1, \pm 2, \ldots .$(54)
$z(\tau ,x) = \sum\limits_{k = - \infty ,k \ne 0}^\infty \,{{\xi }_{k}}(\tau )exp(ikx),\quad {{\xi }_{{ - k}}} = {{\bar {\xi }}_{k}},$(55)
$\frac{{\partial z}}{{\partial \tau }} - \frac{{{{\partial }^{3}}z}}{{\partial \tau \partial {{x}^{2}}}} = {{\kappa }^{{ - 1}}}[{{\kappa }_{1}}(z) + {{\kappa }_{2}}(z)],$Будем говорить, что уравнение из (55) эквивалентно бесконечной системе разностных уравнений (53), если при подстановке в (55) ряда Фурье функции $z(\tau ,x)$ для коэффициентов Фурье получаем систему (53).
Из приведенных выше построений вытекает
Теорема 3. Уравнение (55) эквивалентно системе (53).
Сделаем важное замечание. Уравнение (55) можно существенно упростить, “убрав” слагаемые, содержащие функции $K(z)$. Для этого в системе (51) произведем замены
(57)
${{\xi }_{k}}(\tau ) = {{w}_{k}}(\tau )exp\left( {i{{{(\kappa (1 + {{k}^{2}}))}}^{{ - 1}}}C(k\int\limits_0^\tau {{{{\left| {{{\xi }_{k}}(\tau )} \right|}}^{2}}} d\tau } \right),$4. АСИМПТОТИКА РЕГУЛЯРНЫХ РЕШЕНИЙ ПРИ УСЛОВИЯХ $\alpha \ll 1$ И ПРИ $L \gg 1$
Обратим внимание, что с помощью функций $w(\tau ,x)$ определяются, согласно формуле (34), коэффициенты Фурье ${{\xi }_{k}}(\tau )$ функции $u(t,x,\varepsilon )$, а не сама эта функция.
В этом разделе коротко остановимся на построении нормализованных уравнений для двух случаев: когда $0 < \alpha \ll 1$ и когда $L \gg 1$.
4.1. Здесь предполагаем, что коэффициент $\alpha $ является малым параметром: $\alpha = \varepsilon $ и $0 < \varepsilon \ll 1$. Формально при $\varepsilon = 0$ главная часть уравнения (7) имеет вид
(60)
$\frac{{\partial u}}{{\partial t}} = \kappa \frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}} + au\frac{{\partial u}}{{\partial x}} + b\frac{{\partial u}}{{\partial x}}\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}} + c\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}}.$Учитывая условие регулярности и формулу (60), из (7) приходим к более точной, порядка $O(\varepsilon )$, записи уравнения (7)
(61)
$\frac{{\partial u}}{{\partial t}} = \kappa \frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}} + \ldots + \varepsilon \frac{{{{\partial }^{2}}}}{{\partial {{x}^{2}}}}\left( {\kappa \frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}} + \ldots } \right) + \ldots \;.$(62)
$\frac{{\partial u}}{{\partial t}} = \kappa \frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}} + au\frac{{\partial u}}{{\partial x}} + b\frac{{\partial u}}{{\partial x}}\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}} + cu\frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}} + \varepsilon \delta \left[ {{{b}_{1}}\frac{{{{\partial }^{2}}u}}{{\partial {{x}^{2}}}} + {{b}_{2}}\frac{{{{\partial }^{4}}u}}{{\partial {{x}^{4}}}}} \right] + \varepsilon \kappa \frac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}},$(64)
$\begin{gathered} u = {{\varepsilon }^{{1/2}}}[\xi (\tau )exp(2i{{k}_{0}}x - ik_{0}^{3}t) + \bar {\xi }(\tau )exp( - i{{k}_{0}}x + ik_{0}^{3}t)] + \\ + \;\varepsilon [Aexp(2i{{k}_{0}}x - 2ik_{0}^{3}t) + \bar {A}exp( - 2i{{k}_{0}}x + 2ik_{0}^{3}t)] + {{\varepsilon }^{{3/2}}}{{u}_{3}}(t,\tau ,x) + \ldots , \\ \end{gathered} $И здесь при $\delta = 0$ и ${{a}_{0}} = 0$ все коэффициенты в (65) являются чисто мнимыми, а значит все асимптотические по невязке решения (62), (63) согласно (64) тоже периодические. При $\delta \ne 0$ речь может идти об одном периодическом решении.
Как и выше, можно рассмотреть вопрос о построении асимптотики конечномерных асимптотических по невязке торов. Более интересен вопрос о бесконечномерных торах. В отличие от уравнения (7), здесь удается получить компактную запись системы уравнений, которая играет роль нормальной формы. Покажем это.
Рассмотрим формальный ряд
(66)
$u = {{\varepsilon }^{{1/2}}}\left[ {\sum\limits_{k = - \infty ,k \ne 0}^\infty \,{{\xi }_{k}}(\tau )exp(ikx - i{{k}^{3}}t)} \right] + \varepsilon {{u}_{2}}(t,\tau .x) + {{\varepsilon }^{{3/2}}}{{u}_{3}}(t,\tau ,x) + \ldots ,$(67)
$\frac{{\partial {{u}_{2}}}}{{\partial t}} = \kappa \frac{{{{\partial }^{3}}{{u}_{2}}}}{{\partial {{x}^{3}}}} + \frac{1}{2}a\frac{\partial }{{\partial x}}u_{1}^{2} + \frac{1}{2}b\frac{\partial }{{\partial x}}\mathop {\left( {\frac{{\partial {{u}_{1}}}}{{\partial x}}} \right)}\nolimits^2 + \frac{1}{2}c\frac{\partial }{{\partial x}}\left[ {\frac{{{{\partial }^{2}}}}{{\partial {{x}^{2}}}}(u_{1}^{2}) - 3\mathop {\left( {\frac{{\partial {{u}_{1}}}}{{\partial x}}} \right)}\nolimits^2 } \right].$Это уравнение разрешимо в указанном классе функций при выполнении условия равенства нулю коэффициентов при всех гармониках вида $exp(ikx - i\kappa {{k}^{3}}t)$ в его правой части:
На следующем шаге соберем коэффициенты при ${{\varepsilon }^{{3/2}}}$. В результате получим уравнение для определения ${{u}_{3}}(t,\tau ,x)$:
(68)
$\begin{gathered} \frac{{\partial {{u}_{3}}}}{{\partial t}} - \kappa \frac{{{{\partial }^{3}}{{u}_{3}}}}{{\partial t\partial {{x}^{2}}}} = - \frac{{\partial {{u}_{1}}}}{{\partial \tau }} + \delta \left[ {{{b}_{1}}\frac{{{{\partial }^{2}}{{u}_{1}}}}{{\partial {{x}^{2}}}} + {{b}_{2}}\frac{{{{\partial }^{4}}{{u}_{1}}}}{{\partial {{x}^{4}}}}} \right] + \kappa \frac{{{{\partial }^{5}}{{u}_{1}}}}{{\partial {{x}^{5}}}} + \Phi ({{u}_{1}},{{u}_{2}}), \\ {{u}_{3}}(t,\tau ,x + 2\pi ) \equiv {{u}_{1}}(t,\tau ,x), \\ \end{gathered} $Эта система допускает компактную запись подобно тому, как согласно теореме 3 систему (52) удалось записать в компактной форме (54). Соответствующее выражение для функции
не сложное, но достаточно громоздкое. Линейной части в (68) отвечает для $w(\tau ,x)$ выражение4.2. Асимптотика решений при условии достаточно больших значений параметра $L$.
Кратко остановимся на исследовании регулярных решений в случае, когда
Рассматривается модернизированное уравнение Камасса–Холма, которое отличается от формы записи в (10) тем, что вместо слагаемого $\kappa \tfrac{{\partial u}}{{\partial x}}$ имеем (как и в (9)) слагаемое $\kappa \tfrac{{{{\partial }^{3}}u}}{{\partial {{x}^{3}}}}$. Предполагается, конечно, что выполнены условия ${{f}_{1}} = {{f}_{3}} = 0$ на коэффициенты функции $F(u)$.
Произведем формализующие замены
В результате приходим к краевой задаче(70)
$\frac{{\partial {v}}}{{\partial \tau }} = \kappa \frac{{{{\partial }^{3}}{v}}}{{\partial {{x}^{3}}}} + {{f}_{2}}{v}\frac{{\partial {v}}}{{\partial x}} + O({{\varepsilon }^{2}}),$ВЫВОДЫ
Представлен класс модернизированных уравнений Камасса–Холма, который существенно обобщает классическое уравнение Камасса–Холма и сохраняет ряд его основополагающих свойств. Эти уравнения с периодическими ограниченными условиями имеют однопараметрическое семейство состояний равновесия, в окрестности каждого из которых реализуется критический случай бесконечной размерности в задаче об устойчивости стационаров.
Рассмотрено поведение так называемых регулярных решений в малых окрестностях состояний равновесия. Приведены асимптотические формулы для семейств периодических по времени решений и конечномерных торов.
Исследован вопрос о построении в окрестности состояний равновесия нормализованных уравнений, играющих роль нормальных форм. Эти нормальные формы представляют собой бесконечно мерные системы специальных нелинейных обыкновенных дифференциальных уравнений. Решена задача о компактной записи таких бесконечномерных систем в виде нелинейного уравнения в частных производных. В частности, показано, что такая запись возможна только для классического уравнения Камасса–Холма.
Изучены асимптотические свойства решений при условиях, когда коэффициент $\alpha $ в (7) является достаточно малым параметром или значение периода $L$ в граничных условиях (8) является достаточно большим. В обоих случаях исходная краевая задача сводится в главном к уравнению Кортевега–де Вриза.
Отметим, что построенные выше нормализованные уравнения, как правило, обладают тем свойством, что в окрестности их состояний равновесия тоже реализуется критический в задаче об устойчивости случай бесконечной размерности. Поэтому можно запустить процесс повторной нормализации.
Список литературы
Fokas A.S., Fuchssteiner B. On the structure of symplectic operators and hereditary symmetries // Lettere al nuvo cimento. 1980. № 28. P. 299–303.
Fokas A.S., Fuchssteiner B. Symplectic structure. Their Bäcklund transformations and hereditary symmetries // Phisica D. 1981. № 4. P. 47–66.
Camassa R., Holm D.D. An integrable shallow water equation with Peaked solitons // Phys, Rev. Lett. 1993. № 71. P. 1661–1664.
Camassa R., Hyman J., Holm D. A new integrable shallow water equation // Advances In Applied Mechanics. 1994. V. 31. P. 1–31.
Lakshmanan M. Integrable nonlinear wave equations and possible connections to tsunami dynamics // Tsunami and Nonlinear Waves. Springer. 2007.
Johnson R. Camassa-Holm, Korteweg-de Vries and related models for water waves // J. of Fluid Mechanics 2002. V. 455. P. 63–82.
Constantine A., Lannes D. The hydrodynamical relevance of the Camassa-Holm and degaperis procesi equation // Archive for Rational Mechanics and Analysis. 2009. V. 192. P. 165–186.
Dias F., Milewski P. On the fully-nonlinear shallow-water generalized serre equations // Physics Letters A. 2010 V. 374. P. 1049–1053.
Bhatt R., Mikhailov A.V. On the inconsistency of the Camassa-Holm model with the shallow water theory. arXiv:1010.1932v1.
Korteweg D.J., de Vries G. On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves // Phil. Mag. 1895. V. 5. № 39. P. 422–443.
Kudryashov N.A. On “new travelling wave solutions” of the KdV and KdV-Burgers equotions // Commun. Nonlinear Sci. Numer. Simul. 2009. V. 14. № 5. P. 1891–1900.
Кудряшов Н.А. Методы нелинейной математической физики. М.: “Интеллект”. 2010.
Кащенко С.А. Регулярные и нерегулярные решения в задаче о дислокациях в твердом теле // Теор. и матем. физ. 2018. Т. 195. № 3. С. 362–380.
Kashchenko S.A. Dynamics of a delay logistic equation with slowly varying coefficients // Comput. Math. and Math. Phys. 2018. V. 58. № 12. P. 1926–1936.
Кащенко С.А. Бифуркации в уравнении Курамото–Сивашинского // Теор. и матем. физ. 2017. Т. 192. № 1. С. 23–40.
Кащенко С.А., Преображенская М.М. Бифуркации в обобщенном уравнении Кортевега-де Вриза // Известия в. уч. з. Математика. 2018. № 2. С. 54–68.
Дополнительные материалы отсутствуют.
Инструменты
Журнал вычислительной математики и математической физики