Быстроменяющиеся асимптотические решения некоторых нелинейных эволюционных уравнений в частных производных

Миненков, Дмитрий Сергеевич

Быстроменяющиеся асимптотические решения некоторых нелинейных эволюционных уравнений в частных производных тема автореферата и диссертации по математике, 01.01.03 ВАК РФ

Миненков, Дмитрий Сергеевич АВТОР

кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ

Москва МЕСТО ЗАЩИТЫ

2014 ГОД ЗАЩИТЫ

01.01.03 КОД ВАК РФ

Диссертация по математике на тему «Быстроменяющиеся асимптотические решения некоторых нелинейных эволюционных уравнений в частных производных»

Автореферат диссертации на тему "Быстроменяющиеся асимптотические решения некоторых нелинейных эволюционных уравнений в частных производных"

На правах рукописи

Миненков Дмитрий Сергеевич

БЫСТРОМЕНЯЮЩИЕСЯ АСИМПТОТИЧЕСКИЕ РЕШЕНИЯ НЕКОТОРЫХ НЕЛИНЕЙНЫХ ЭВОЛЮЦИОННЫХ УРАВНЕНИЙ В ЧАСТНЫХ ПРОИЗВОДНЫХ.

Специальность 01.01.03 - Математическая физика

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

3 П ЯНВ 2014

Москва, 2014

005544738

Работа выполнена в лаборатории механики природных катастроф Института проблем механики имени А.Ю. Ишлинского РАН

Научный руководитель доктор физико-математических наук,

профессор Сергей Юрьевич Доброхотов

Официальные оппоненты доктор физико-математических наук,

профессор Леонид Анатольевич Калякин

доктор физико-математических наук, профессор Владимир Григорьевич Данилов

Ведущая организация Санкт-Петербургское отделение Математического

института имени В.А. Стеклова РАН

Защита состоится 20 февраля 2014 г. в 15 ч. 30 мин. на заседании диссертационного совета Д 501.002.10 при Московском государственном университете имени М.В. Ломоносова по адресу: 119991, г. Москва, Ленинские горы, МГУ, дом 1, стр. 2, физический факультет, СФА.

С диссертацией можно ознакомиться в фундаментальной библиотеке МГУ имени М.В. Ломоносова.

Автореферат разослан 2014 г.

Ученый секретарь

диссертационного совета Д 501.002.10 доктор физико-математических наук профессор ' П.А. Поляков

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

С появлением и развитием таких программных пакетов, как МаЛетайса, МаШЬаЬ и им подобные, возникла возможность компьютерной реализации быстрых аналитико-численных алгоритмов для моделирования волновых процессов, что позволяет анализировать зависимость решения от параметров в режиме "онлайн". Однако, существующие формулы для асимптотических решений не всегда подходят для подобных целей, и есть потребность в модификации существующих и получении новых формул для асимптотических решений.

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

Общая методика исследования основана на сочетании адиабатического и квазиклассического приближений, методов осреднения и теории возмущений.

Научная новизна определяется следующими основными результатами:

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

2. Построены однофазовые (формальные) асимптотические решения для нелинейного уравнения Клейна-Гордона и уравнения Кортевега-де-Фриза, которые являются регулярными относительно перехода от "сильнонелинейного" случая к "слабонелинейному"; Доказано, что для этих уравнений фазовый сдвиг в анзаце Кузмака-Уизема определяется из задачи для системы уравнений Уизема с соответствующим образом подобранными начальными условиями;

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

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

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

Научная и практическая ценность. Работа носит теоретический характер. Асимптотические методы решения задач математической физики сами по себе представляют теоретический интерес. Асимптотические решения часто объясняют ключевые свойства точных решений, получаемых численно или экспериментально. Рассматриваемые задачи описывают несколько практически важных явлений, например, поведение жидкости в каналах или динамику волн цунами. Полученное представление, равномерное относительно перехода от "сильнонелинейного" случая к "слабонелинейному", позволяет описывать динамику цуга волн, когда на хвосте цуга амплитуда решения небольшая и ситуация "слабонелинейная", а у пика - "сильнонелинейная".

Личное участие автора. Задача построения асимптотических решений одномерной нелинейной системы уравнений мелкой воды с вырождающейся скоростью изучена автором самостоятельно. Некоторые результаты получены совместно с научным руководителем С.Ю. Доброхотовым (Институт Проблем Механики РАН, Московский Физико-Технический Институт) и С.Б. Медведевым (Институт вычислительных технологий СО РАН, Московский государственный университет). Здесь вклад автора заключается в проведении конкретных вычислений и доказательств, а также в реализации полученных асимптотических формул на компьютере.

Апробация работы. Основные результаты работы были представлены на Международных конференциях "Дни дифракции" (в 2010-2012 годах, Санкт-Петербург), на международной конференции "Mathematical analysis of Asymptotic and Applications" (в 2010 году, IIMAS, UNAM, Mexico City, Mexico), на научной конференции МФТИ (2012, Москва), на международной конференции "Нелинейный анализ и спектральные задачи" (2013, Уфа) и на международной конференции "Дифференциальные уравнения. Функциональные пространства. Теория приближений", посвященной 105-летию со дня рождения C.JI. Соболева (2013, Новосибирск). Также результаты били представлены и обсуждались на научных семинарах под руководством В.М. Бабича в ПОМИ РАН (2012, Санкт-Петербург) и под руководством Л.А. Калякина в Институте Математики с ВЦ УНЦ РАН (2013, Уфа).

Публикации. Основное содержание работы отражено в четырех публикациях, список которых приведен в конце автореферата. Работы опубликованы в изданиях, входящих в утвержденный ВАК перечень ведущих рецензируемых научных

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

Структура и объем диссертации. Диссертация состоит из введения, 2 глав и списка литературы. Материал диссертации изложен на 105 страницах. Список литературы содержит 68 наименования.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

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

Фазовый сдвиг в анзаце Кузмака-Уизема

В первой главе проводится обзор методов осреднения, разработанных Н. М. Крыловым, Н. Н. Боголюбовым, Ю. А. Митропольским и другими, и приводятся результаты, связанные с процедурой определения фазового сдвига в анзаце Кузмака-Уизема при построении быстроосциллирующих асимптотических решений в нелинейных уравнениях.

Согласно стандартной схеме, примененной Г. Е. Кузмаком (в 1959) для обыкновенных дифференциальных уравнений и Г. Б. Уиземом (1965) для уравнений в частных производных, главный член асимптотического разложения определяется из нелинейного уравнения, а поправки к нему из линейных уравнений. Главный член разложения может быть представлен в форме где фаза 5,

"медленно меняющийся" параметр или параметры I и так называемый фазовый сдвиг ф определяются из системы "осредненных" уравнений - условий совместности уравнений для поправок. Вопрос о вычислении фазового сдвига возник давно и обсуждался для уравнений в частных производных С. Ю. Доброхотовым и В. П. Масловым (1981), а также Р. Хаберманом (1988). Для обыкновенных дифференциальных уравнений этот вопрос обсуждался в работах Р. Хабермана и соавторов.

Основная сложность при определении фазового сдвига состоит в том, что в "сильнонелинейном" случае для вычисления фазового сдвига достаточно рассмотреть уравнение на первую поправку к главному члену разложения, а в линейном и "слабонелинейном" случае необходимо также исследовать уравнение для второй поправки и для этого необходимо определить первую поправку, что является нетривиальной задачей. Это препятствие имеет топологический характер: при переходе от "сильнонелинейного" случая к "слабонелинейному" размерность коядра соответствующего линейного оператора в вариациях уменьшается на 1. Практический результат, основанный на идеях работ С. Ю. Доброхотова и В. П. Маслова, и работ Р. Хабермана, заключается в том, что можно учесть фазовый сдвиг ф, если рассмотреть его как часть фазы в, и при этом соответствующим образом изменить начальные данные для уравнения на I. Построенные представления для главного члена асимптотического разложения являются инвариантными относительно перехода от "сильнонелинейного" случая к "слабонелинейному". Это позволяет, например, описывать динамику цуга волн, когда на хвосте цуга ситуация "слабонелинейная", а у пика - "сильнонелинейная".

В пункте 1.1 первой главы рассматривается одномерный нелинейный осциллятор с потенциалом V(x, ет), медленно зависящим от времени т, и малым неконсервативным членом ед(х,х,гт), описывающим трение или (и) "накачку", где е < 1 - малый параметр. Соответствующий математический объект - задача Коши для неизвестного положения х(т, е) и скорости (момента) р(т, е) = х =

р = -Vx(x,£t) -ед{р,х,гт), х = р, <=>• х + f(x,er) + ед(х,х,ет) = 0, (0.1) х\т=о = я0. р\т=о = ¿|т=о = Р°- (°-2)

Здесь / = Vx = |V{x,t), д - гладкие функции, д(р,х,ет) - нечетная функция р. Обычно переменную т называют "быстрым" временем в противоположность "медленному" времени t = ет.1 Задачу (0.1), (0.2) можно переписать, используя медленное время £, тогда она примет форму так называемой сингулярно возмущенной задачи:

E2~ + Hx,t) + eg(e^,x,t)= 0, x|í=0 = х°, e§|t-o=P°. (0.3)

Мы предполагаем, что по крайней мере локально, потенциал V имеет форму потенциальной ямы на интервале (x¡, хг) с минимумом в точке х = a;m¡„ 6 (x¡, xr). Потенциал V, функция / так же как и начальные значения х°,р° могут зависеть регулярным образом от некоторого действительного параметра (или параметров) ц 6 [0, /¿о], например, можно рассмотреть потенциал V = ui(t)x2/2 + ^,V1(x,t),V1(x,t) = 0(ху). Параметр /х может быть связан с параметром е, в простейшем случае ц = е. Обычно, ситуацию, когда V = w(t)x2/2 + eVl{x, t), называют слабонелинейным случаем в противоположность силънонелинейному случаю, когда 3jr|£=o ф 0. Хорошим примером слабонелинейного случая служит осцилятор Ван-дер-Поля: х + х — е ах (1 — х2) = 0. В качестве примера сильнонелинейного случая мы рассмотрим задачу Коши для физического маятника с медленноменя-ющейся частотой w(í) и нелинейным неконсервативным членом, описывающим трение при a(t) < 0 и накачку при a(t) > 0: х + w2(¿) sinx — ea(t) х cosa; = 0.

Упомянутое разделение на "слабо-" и "сильнонелинейные" случаи является несколько искусственным и зависит от начальных условий х°,р°. Действительно, предполагая, что V = + 7(*)т + 0(х6) и х° = s/ix°,p° = после замены переменных х = y/ix, р = л/ер мы получим слабонелинейное уравнение с V = u)(t) y + £7(í)t + О (г2). Для упрощения обозначений, мы не будем включать параметр ц в V, д, х°,р° и другие функции и параметры, но будем держать в уме такую возможность.

Понятно, что обычно точное решение задачи Коши (0.1), (0.2) неизвестно, и можно говорить только об асимптотическом решении при £ < 1. Интересно, что

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

Хир(г)Л(т)

Рис. 1: The potential V(x,t) and E for various t: t = 0.1, t = 0.2 (left). The example of a fast oscillating solution (right).

построение таких асимптотик возможно на болыпйх временах, по крайней мере г ~ j. В переменных медленного времени t = ет "большие времена" означают t G [0,0(1)], и желаемые решения являются быстроосциллирующими (см. рис. 1). В теории нелинейных колебаний существуют различные асимптотические методы, которые обычно называются методами осреднения. Существует большое число статей и монографий, посвященных методам осреднения, подробную библиографию можно найти в книгах В. И. Арнольда, В. В. Козлова и А. И. Нейштадта, Н. Н. Боголюбова и Ю. А. Митропольского, J. D. Cole, J. Kevorkian. Несмотря на то, что методы осреднения основаны на похожих идеях, их реализация различается. Можно условно разбить их на две группы: методы, связанные с "заменой переменных" и методы "прямого расчета". Так метод Крылова-Боголюбова, теория KAM и им подобные попадут в первую группу, а так называемый метод Кузмака и его обобщение на случай уравнений в частных производных - метод Уизема или нелинейный метод ВКБ следует отнести ко второй группе.

Заметим, что задачи, рассматриваемые в первой главе относятся к так называемым сингулярно возмущенным, которые допускают асимптотики с составляющей типа пограничного слоя. Исследованию последних посвящены, в частности, работы В. Ф. Бутузова, А. Б. Васильевой, М. И. Вишика, С. А. Ломова, JI. А. Люстерни-ка, H.H. Нефедова, А.Н. Тихонова, и других. В отличие от таких асимптотик, в диссертации рассматриваются асимптотики, имеющие быстрые осцилляции.

Задача (0.1), (0.2) решена в статье F. J. Bourland and R. Haberman2 (см. также книгу J.D. Cole, J. Kevorkian и статью A.M. Ильина3), посредством метода Кузмака-Уизема, с использованием некоторых идей, предложенных в работе С.Ю. Доброхотова и В. П. Маслова4. Окончательный ответ представлен в двух разных формах: для слабонелинейного случая и для сильнонелинейного случая (при этом нет перехода от одного к другому). Основной результат данной главы - единая асимптотическая формула, работающая в обоих случаях. В этом разделе сформу-

2F. J. Bourland and R. Haberman, SIAM J. Appl. Math, 48, 1988, pp.737-748

3A.M. Ильин, ТМФ, 118:3, 1999, 383-389

4С.Ю. Доброхотов, В. П. Маслов, Итоги науки и техники. ВИНИТИ, т. 15, 1980

лирована и доказана следующая теорема [1]: Теорема 1

Пусть существуют положительные числа а, Ь, to, такие что для I £ [о, 6], í €

[ОМ-

(11) вектор-функция (Х{в, I, t), Р(в, I, t)) - 2тг-периодическая по в и является решением невозмущенной системы с "замороженным " временем t:

Р2

Р = —Hx(x,t), х = Нр, H = ^ + V(x,t), 6,1 - переменные "действие-угол" для невозмущенной системы:

/(«) = - [ + y/2(H{I,t)-V(x,t))dx,

4-ОС/.0, «аMmm-vbm-

и x±(t) - корни уравнения V(x, t) = E(I, t);

(12) для t € [0, ío] существует решение I(t, e) € [a, b], <í>(í, e) задачи Коши для "осредненных " уравнений

fr = Wl fr = ^ f = f =

1 Г2*

G(I,t) = — jf Xe(e,I,t)-g(P(e,Lt),X(e,I,t),t)de, с начальными условиями:

Z|t=o = Л*°,Р°) + Ф|г=о = еф°(х°,р0),

где ф°,1° определяются из условий Х(ф°,1°, 0) = х°, Р(ф°,1°, 0) = р°, и поправка IJ равна

1°(х°,Р°) = -

(П (Г

J^j J* (у,(Х(в, Л 0), 0) - ¿ jT Vt(X(<r, 0), О)dff) de+

jf [хе{в, Ia, 0)g{P{9, Io, 0), Х(в, Io, 0), 0) - G(I°, 0))м.

Тогда для главного члена асимптотического разложения х(т,е), р{т, е) для задачи Коши (0.1), (0.2) на временах т G [0, tQ/h] справедливо следующее пред-

ставленые

х(т, е) = Х(Ф(££Т'£),Х(£г,£);ег) + 0(е), р(г,е) = + 0(е).

Основная цель данного раздела в том, чтобы показать, что этот (и более сильный) результат может быть доказан при помощи теоремы А. И. Нейштадта5, которая развивает группу методов "замены переменных", а также сравнить эффективность обоих подходов: методов "замены переменных" и методов "прямого счета". По нашему мнению, методы осреднения, связанные с заменой переменных, лучше подходят для рассмотренной задачи Коши (0.1), (0.2), а привлекательность метода Кузмака-Уизема (или нелинейного метода ВКБ) состоит в том, что он естественным образом обобщается на случай уравнений в частных производных.

В пункте 1.2 первой главы рассмотрены однофазовые (формальные) асимптотические решения в форме Кузмака-Уизема для нелинейного уравнения Клейна-Гордона и уравнения Кортевега-де-Фриза. В этом случае главный член асимптотического решения также ищется в форме анзаца Кузмака-Уизема Х(3(х,1)/е + Ф(ге,¿)), /(х,ж,+ 0(е), где е < 1 - малый параметр, фаза и медлен-

но меняющиеся параметры 1{х,Ь) находятся из системы "осредненных" уравнений Уизема. Вообще говоря, уравнение для фазового сдвига Ф(х, £) получается из исследования второй поправки к главному члену. При этом соответствующая процедура нахождения фазового сдвига неравномерна относительно перехода к линейному (и слабонелинейному) случаю.

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

Именно, сформулированы и доказаны следующие теоремы. Рассмотрим задачу Коши для нелинейного уравнения Клейна-Гордона:

где У(и,х, с(х,Ь), и°(х) и г>°(х) е Я - гладкие функции; < € Я, х е Яп, 0 < е < 1 - малый параметр. От функции У(и, х, ¿) потребуем выполнения условия, обеспечивающего существование у уравнения (0.4) быстроосциллирую-щих вещественнозначных асимптотических решений: пусть существуют гладкие функции а(х, ¿) < Ъ(х, ¿), такие что при каждых фиксированных х,Ь на интер-

5А.И. Нейштадт, ПММ, 48:2, 1984, 197-204

е2{иа - С2(х, г)Аи) + Уи(и, X, Ь) = 0, и(х,Ь,£) |(=о = и°(х), ещ(х, ¿,е)|(=о = Ъ°(х)

(0.4) (0.5)

вале а < Е < Ь уравнение V(u, x,t) = Е имеет по крайней мере два решения Mmin(x,t, E),umax(x,t, Е),ит[П < ипшк, гладко зависящих от x,t,E и таких, что

KWrnin, X, t) ф 0, K(-iimax, X, t) ф О И V(u, X,t) < Е При Wmin < U < Mmax.

Теорема 2 Пусть существует положительное число to и гладкие функции a(x,t) < b(x,t), такие что для t € [0, ¿о]> I(x>t) 6 [a(x,t),b(x,t)\ выполнены условия:

(ii) Существует 2ж-периодичное по переменной в решение Х(0, х, t) = Х(в, I(x,t),x,t) уравнения нулевого приближения Г22(ж, t)Xgy + Vu(X,x,t) = 0. Это решение задается неявной формулой

rX(U,x,t) Jz

9 = ±Sl(x,t) .

Лт1„ y/2{E(I,x,t)-V(z,x,t))

Ги max Л?

П{х, t) = 7Г/ / ■

где для параметризации использована переменная типа действия: I{x,t) = ± y/2(E(I,x,t)-V(z,x,t))dz;

(¡2) Существует решение S(x, t, е), I(x, t, е) задачи Коши для системы Уизема:

Sf-c2(x,t)(VS)2 = n(I,x,t), д г I о / .,/„ I

dt№(I,x,t) J v ' ' f2(I,x,t)

(где операторы дифференцирования во втором уравнении действуют также на функцию I(x, t, е)), удовлетворяющее начальным условиям:

S(x,t,e) |<=о = £Ф°0г), ï(x,t,e)\t=0 = l°(x) + el°(x),

где поправка определяется формулой

i ГЩх)

^ = h Vtimx>0)) ~

Тогда главный член асимптотического разложения для решения исходной задачи Коши (0.4), (0.5) на временах t € [0, ¿о] определяется формулой с "нулевым фазовым сдвигом":

и(х, t, е) = Х(+ Ф(х, t), 1{х, t), X, t) + 0{е).

Рассмотрим теперь уравнение Кортевега-де-Фриза

щ — 6иих + е2иххх = 0, х £ R, t > 0 (0.6)

где 0<£<1- как и раньше, малый параметр. Некоторые формальные асимптотические решения этого уравнения представляются в виде анзаца Кузмака-Уизема. Дифференциальное уравнение, из которого определяется главный член асимптотического разложения Х(9, х, t), в этом случае имеет вид:

П{х, t)Xg - 6К(х, t)XXe + К3(х, t)Xm = 0,

где K(x,t) = Sx(x,t), Г1(х, t) = St{x,t) - волновое число и частота. Здесь также в - периодическая переменная, х, t - "замороженные" переменные. Для параметризации функция Х(в,х,t) удобно использовать корни Е0, Ei, Е2 многочлена R2(z,x,t) = z3 - §^z2 + Ax{x,t)z + A2(x,t) = (z - E0)(z - E{){z - E2) (см. работы G.B. Whitham (1965), Н.И. Ахиезера (1970), В.Б. Матвеева (1976), A.B. Гуревича и Л.П. Питаевского (1971), Б.А. Дубровина, В.Б. Матвеева и

C. П. Новикова (1976), А. Р. Итса и В. Б. Матвеева (1976), Н. Flaschka, М. Forest,

D. McLaughlin (1980), С. Ю. Доброхотова и В. П. Маслова (1980), И. М. Кричевера (1988), R. Haberman (1988) и других). Формула для функции Х(9,Е) найдена в работе А. Р. Итса и В. Б. Матвеева6, ее можно представить в терминах 0-функции Якоби, как в вышеуказанной статье С. Ю. Доброхотова и В. П. Маслова. Система уравнений Уизема при параметризации через концы зон Eq, Ei, Е2 получена в работе Н. Flaschka, М. Forest и D. McLaughlin7:

itE> = -A'WlTxE" Z = 0'1'2-

(Формулы для коэффициентов Л¡(Е) приведены в тексте диссертации.) Волновое число К и частота Г2 выражаются через концы зон, после чего фаза находится интегрированием уравнения dS = Cl(x,t)dt + k(x,t)dx. Чтобы получить уравнение на фазовый сдвиг Ф, необходимо исследовать уравнение второго приближения. В Теореме 3 рассматривается задача Коши для уравнения Кортевега-де-Фриза с начальными условиями из того же класса, что и строящиеся асимптотики. Показано, что, так же как и в случае волнового уравнения, уравнение на фазовый сдвиг можно не решать, если включить его в фазу S(x,t,e) = S(x,t) +еФ(а:, t) и одновременно скорректировать концы зон E(x,t,e) = E(x,t) + s SE(x,t). Приведены формулы для поправок 6Е° = 6E\t=о к начальным условиям для концов зон.

Полученное представление для главного члена асимптотического разложения (Теоремы 1, 2 и 3) является равномерным относительно перехода к "слабонелинейному" случаю, и поэтому его можно использовать для расчета динами-

'А.Р. Итс, В. Б. Матвеев, в сб. Пробл. мат. физ. № 8, Л., Ленингр. ун-т, 1976, 70-92

7Н. Flaschka, М. Forest и D. McLaughlin, Comm. Pure. Appl. Math 33 (1980), pp. 739-784

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

Системы уравнений мелкой воды

Во второй главе рассматривается система нелинейных уравнений для волн на воде в приближении мелкой воды. Одномерные уравнения мелкой воды являются одной из простейших моделей распространения жидкости в каналах (см. книгу J. J. Stoker, 1958). Кроме того, эти уравнения могут быть использованы для исследования цунами (см. работы E.H. Пелиновского и P.X. Мазовой). В случае ровного дна уравнения могут быть сведены к линейным уравнения с помощью преобразования годографа (см. книгу R. Courant, 1962). При постоянном наклоне дна уравнения также могут быть сведены к линейным уравнениям с помощью преобразования годографа, после приведения их к римановым инвариантам. Однако, возможны дальнейшие преобразования зависимых и независимых переменных, которые позволяют выразить параметрически решения исходных уравнений через два параметра и функцию от них (см. статью G. F. Carrier, H. P. Greenspan, 1958s). При этом функция является решением линейного волнового уравнения. Преимущество такого представления состоит, во-первых, в том, что решение нелинейной системы первого порядка сводится к решению линейного волнового уравнения, во-вторых, подвижная точка уреза (точка, в которой глубина жидкости равна нулю) переходит в нулевую точку для линейного уравнения.

Хотя рассматриваемая система является бездисперсионной, но вопрос о фазовом сдвиге все равно возникает. Для уравнений мелкой воды фазовый сдвиг связан со скачком индекса Маслова при прохождении особой точки специального типа, в которой вырождается скорость распространения волны. Это объясняет известную метаморфозу профиля волны, которая описывается преобразованием Гильберта и состоит например в том, что "размазанное" ¿-образное решение переходит в так называемую N-волну. В целом оказывается, что решение системы мелкой воды другого типа, чем рассматриваемые в Главе 1, и никаких новых эффектов, связанных с фазовым сдвигом, нелинейность не дает.

В пункте 2.1 второй главы приведены точечные преобразования, связывающие три системы уравнений первого порядка [4]: уравнения мелкой воды над ровным дном D(z) = О

du du dh g dh d[hu] ^ dt dz dz ' dt dz

8G.F. Carrier, H. P. Greenspan, J. Fluid Mcch., 4, 1958, 97-109

уравнения мелкой воды над дном постоянного уклона В(х) = с1{х) = х:

dv dv drj дг) д[(т) + x)v]

m+vdï + dï = 0> dt дх =

и линейной системы, которая получается формальной линеаризацией:

dU + dN=0 dN | д[уЦ] _ о дт ду ' дт ду

Все найденные преобразования сохраняют дивергентный вид уравнений. Также приведены преобразования, связывающие линеаризованную систему и линейное волновое уравнение, выведенное в работе G. F. Carrier, H. P. Greenspan. Кроме того, представлены формулы специального решения линеаризованной системы в виде волны, бегущей с переменной скоростью (полученного С. Ю. Доброхотовым и С. Я. Секерж-Зеньковичем), и показана связь полученного решения с решением двумерного волнового уравнения с постоянной скоростью звука.

В пункте 2.2 второй главы рассматривается система нелинейных уравнений мелкой воды в одномерном случае над неровным дном D(x) = с2 (х) (с(х) - скорость распространения волн):

771 + (D(x)u + 7]u)x = 0, щ + (Г) + U2/2)x = 0, x > xmia(t), t 6 R (0.7)

для уровня свободной поверхности т] = rj(x,t) и скорости и = u(x,t) (см. например, книги J.J. Stoker и E.H. Пелиновского). Считаем, что точка х = 0 соответствует берегу, а именно sign(Z?(x)) = sign(a:). Значения D(x) < 0 отвечают "отрицательной глубине", т.е. высоте берега над уровнем водоема. Задача состоит в построении некоторых ограниченных решений рассматриваемых уравнений в области с переменной левой границей х € [ïmm(£), оо), где xm\n(t) определяется из условия D^roin) + v{xmm, t) = 0.

Решения этой задачи в случае линейного дна, когда D(x) = -у2х изучались в работах J. J. Stoker, Е. Н. Пелиновского и P. X. Мазовой, G. F. Carrier, H. P. Greenspan, T. Vukasinac и P. Zhevandrov (2002), С.Ю. Доброхотова и Б. Тироцци (2010), причем, в последней получено семейство точных решений, представленных в виде алгебраических функций.

Мы будем рассматривать дно в виде D(x) = y2xf(ex), /(0) = 1, =

0(1), где характерный размер "неровностей" дна D(x) определяется малым параметром £ «С 1. Изменением масштаба (х —> x, t —> £/7, г] —> 72?], и —> 7и) можно сделать 7=1. Решаем задачу Коши с начальным условием, удаленным от берега на расстояние порядка длины волны:

rj\t=0 = rçV), тР(х) = иУ(х-1), w|(=o = 0. (0.8)

Рис. 2: Главный член М0(г, т) для линейной системы над линейным дном (штрих) и главный член щ(х, ¿) для нелинейной системы над дном И{х) = £(1—0.18т(27га-)) (сплошная) в окрестности фокальной точки: для времен т, £ = 1,2,3.

Здесь - гладкая финитная функция вирр V С [—1/2,1/2], < 1, пара-

метр V характеризует амплитуду в начальный момент времени. Будем считать, что амплитуда волны в начальный момент достаточно мала по сравнению с глубиной V <С £>(1), чтобы волна дошла до берега без опрокидывания. Это достаточно часто бывает для волн цунами (как показано в работах Е. Н. Пелиновского и Р. X. Мазо-вой).

Основной результат заключается в следующем (см. [3]). Вблизи берега существует замена переменных (х,£) —» (г,т), (т](х,Ь),и(х,1)) —» (Ы{2,т),и{г,т)), которая сведет исходную задачу к задаче Коши для нелинейной системы, представляющей собой слабовозмущенную линейную систему:

Мт + (ги)2 = еЬ(г,М,и)ги, ит + N. = е/2(г, И, ЩИ, (0.9) Щт=0 = №{г), ¿7|т=о = 0. (0.10)

Указанная замена переменных строится как композиция двух замен. Сначала проводится замена переменных для "выравнивания дна", которая переводит исходную систему в систему мелкой воды с дном постоянного уклона и с малой правой частью. К полученной системе применяются формулы линеаризации методом Кариера-Гринспена (см. также книгу Е. Н. Пелиновского) в форме, представленной в статье С. Ю. Доброхотова и Б. Тироцци9:

г = у + ч(у,{), т = г--ш(у,1), Ы{г,т) 1/(г,т) = ю.

Преобразование Кариера-Гринспена переводит задачу с переменным левым пределом утт : ут¡п + ? = 0 в задачу (0.9), (0.10) с постоянной левой границей 2,шП = 0, поскольку. Для этой задачи получено (формальное) асимптотическое решение методом регулярной теории возмущений.

Главный член асимптотики и поправки определяются из системы неоднород-

®С.Ю. Доброхотов, Б. Тироцци, УМН, 65:1(391), 2010, 185-186

ных линейных уравнений со скоростью с = ~/\fz (с дном постоянного уклона D = 72z). Для корректной постановки задачи Коши для системы линейных уравнений, необходимо наложить дополнительное условие на поведение решения при х —¥ 0, поскольку скорость с{х) вырождается в точке х = 0. Как было сформулировано T. Vukasinac и P. Zhevandrov (2002) и доказано С. Ю. Доброхотовым,

B.Е. Назайкинским и Б. Тироцци (2010), такая задача будет корректно поставлена, если потребовать с(х)г)х{х) G L2(R+), что в сущности эквивалентно тому, что энергия решения должна быть конечна. Для задачи Коши для линейной системы с локализованными начальными условиями известно точное решение в интегральной форме и достаточно эффективные асимптотические формулы для этого решения вне окрестности точки г = 0 (см. работы Е. Н. Пелиновского и P. X. Мазовой,

C.Ю. Доброхотова, В.Е. Назайкинского и Б. Тироцци10). Как показано С.Ю. Доброхотовым, В. Е. Назайкинским и Б. Тироцци, эта точка является фокальной точкой или точкой поворота.

В результате в пункте 2.2 главы 2 сформулирована и доказана следующая теорема (см. [3]).

Теорема 4 Пусть скорость с(х) (дно D(x) = с(х)2) достаточно мало отклоняется от \[х (от дна постоянного уклона х), и амплитуда начальной волны достаточно мала. Пусть кроме того начальные условия лежат в пространстве Шварца на полупрямой: №(z), U°(z) 6 <S[0, оо) (№, U° € С°°[0, оо) и убывают на бесконечности быстрее любой степени).

Тогда исходная задача сводится к задаче (0.9), (0.10) с границей z = 0, для которой можно построить формальное асимптотическое решение методом регулярной теории возмущений в виде ряда по малому параметру е: N(z, т, е) = ¿^о^Л^л.т), U(z,T,e) = Y^kLoehuk(z--T) - в классе 5[0,оо) по z, на временах т = 0(1). Обратная замена переводит его в формальное асимптотическое решение для исходной задачи (0.7), (0.8) с переменной границей xmia(t) - в классе (S[a;min(i), оо) по х, на временах t = 0(1).

С помощью построенных асимптотических формул и точных решений специального вида для линейных уравнений мелкой воды рассмотрено отражение волны от берега в случае нелинейной системы (см. Рис. 2). Как показано, для уравнений мелкой воды фазовый сдвиг связан со скачком индекса Маслова при прохождении особой точки специального типа, в которой вырождается скорость распространения волны. Это объясняет известную метаморфозу профиля волны (обнаруженную в работе E.H. Пелиновского и P.X. Мазовой (1992)), которая описывается преобразованием Гильберта и состоит в том, что "размазанное" ¿-образное падающее решение переходит в так называемую N-волну.

,0С. Ю. Доброхотов, В. Е. Назайкинский, Б. Тироцци, Алгебра и анализ 21 (6), 67-90 (2010)

ОСНОВНЫЕ РЕЗУЛЬТАТЫ ДИССЕРТАЦИИ

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

РАБОТЫ, ОПУБЛИКОВАННЫЕ ПО ТЕМЕ ДИССЕРТАЦИИ

1. S.Yu. Dobrokhotov, D.S. Minenkov, "On Various Averaging Methods for a Nonlinear Oscillator with Slow Time-dependent Potential and a Nonconservative Perturbation", Regular and Chaotic Dynamics, 2010, Vol. 15, No. 2-3, pp. 285-299.

2. С.Ю. Доброхотов, Д.С. Миненков, "О фазовом сдвиге в анзаце Кузмака-Уизема", ТМФ, 2011, 166 (3), 350-365.

3. Д.С. Миненков, "Асимптотические решения одномерной нелинейной системы уравнений мелкой воды с вырождающейся скоростью", Математические заметки, 2012, 92 (5), 721-730.

4. С.Ю. Доброхотов, С.Б. Медведев, Д.С. Миненков, "О заменах, приводящих одномерные системы уравнений мелкой воды к волновому уравнению со скоростью звука с2 = х", Математические заметки, 2013, 93 (5), 725-736.

БЛАГОДАРНОСТИ

Я очень признателен В. Е. Назайкинскому и А. И. Шафаревичу за моральную поддержку и полезные дискуссии во время подготовки работ. Особую благодарность я выражаю научному руководителю С.Ю. Доброхотову за помощь, оказанную автору за время обучения в аспирантуре. Я также благодарен С. Б. Медведеву и С. Я. Секерж-Зеньковичу за дискуссии и ценные советы. Результаты диссертации были получены в рамках проектов РФФИ 08-01-00726 и 11-01-00973.

Подписано в печать: 11.01.14

Объем: 1,0 п.л. Тираж: 120 экз. Заказ № 201 Отпечатано в типографии «Реглет» г. Москва, Ленинский проспект, д.2 (495) 978-66-63, www.reglet.ru

Текст научной работы диссертации и автореферата по математике, кандидата физико-математических наук, Миненков, Дмитрий Сергеевич, Москва

Российская академия наук Федеральное государственное бюджетное учреждение науки Институт проблем механики им. АЛО. Ишлинского Российской академии наук (ИПМех РАН)

0420145Б9У7

на правах рукописи

Миненков Дмитрий Сергеевич

БЫСТРОМЕНЯЮЩИЕСЯ АСИМПТОТИЧЕСКИЕ РЕШЕНИЯ НЕКОТОРЫХ НЕЛИНЕЙНЫХ ЭВОЛЮЦИОННЫХ УРАВНЕНИЙ

В ЧАСТНЫХ ПРОИЗВОДНЫХ.

Специальность: 01.01.03 "Математическая физика"

ДИССЕРТАЦИЯ на соискание ученой степени кандидата физико-математических наук

Научный руководитель доктор физико-математических наук профессор С.Ю. Доброхотов

Москва, 2013

Оглавление

Введение 4

Общая характеристика работы....................................4

Краткое содержание работы ......................................8

1 Фазовый сдвиг в анзаце Кузмака-Уизема 25

1.1 Сравнение различных методов осреднения, па примере ангармонического осциллятора с медленно зависящим от времени потенциалом и некопссрвативным возмущением .............................. 25

1.1.1 Постановка задачи и формулировка результатов . 25

1.1.2 Возмущение начальных данных и потенциала, "неустойчивость" асимптотик........... 35

1.1.3 Доказательство теоремы 1.............. 38

1.1.4 Пример........................ 44

1.1.5 Метод Кузмака-Уизема............... 46

1.2 Фазовый сдвиг в анзаце Кузмака-Уизема для нелинейного уравнения Клейна-Гордона и уравнения Кортевега-де-Фриза........................... 52

1.2.1 Нелинейное уравнение Клейна-Гордона..... 52

1.2.2 Однофазовое решение уравнения Кортевега-де-Фриза......................... 63

2 Системы уравнений мелкой воды 72

2.1 Замены, приводящие одномерные системы уравнений

мелкой воды к волновому уравнению со скоростью звука

с2 =х............................. 72

2.1.1 Системы нелинейных уравнений мелкой воды . . 72

2.1.2 Линеаризация систем уравнений мелкой воды, основанная на точечных преобразованиях .... 74

2.1.3 Преобразования Кариера-Гринспена и близкие исследования .................... 78

2.1.4 Точные решения линейной системы........ 80

2.1.5 Нули Якобианов и "дополнительные" решения . 84

2.2 Асимптотические решения одномерной нелинейной системы уравнений мелкой воды с вырождающейся скоростью общего вида...................... 87

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

2.2.2 Замена для выравнивания с2(.г) (дна)....... 89

2.2.3 Линеаризация, основанная па преобразовании Кариера-Гринспена................. 90

2.2.4 Построение асимптотики с помощью метода теории возмущений................... 92

2.2.5 Пример, иллюстрирующий отражение волны от берега в нелинейном случае............ 96

Список литературы 98

Введение

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

С появлением и развитием таких программных пакетов, как МаШетайса, МаЛЬаЬ и им подобные, возникла возможность компьютерной реализации быстрых аналитико-численных алгоритмов для моделирования волновых процессов, что позволяет анализировать зависимость решения от параметров в режиме "он-лайн". Однако, существующие формулы для асимптотических решений не всегда подходят для подобных целей, и есть потребность в модификации существующих и получении новых формул для асимптотических решений.

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

ствующей берегу). Построенные асимптотические решения эффективно реализуются на компьютере.

Общая методика исследования основана па сочетании адиабатического и квазиклассического приближений, методов осреднения и теории возмущений.

Научная новизна определяется следующими основными результатами:

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

2. Построены однофазовые (формальные) асимптотические решения для нелинейного уравнения Клейна-Гордона и уравнения Кортевега-де-Фриза, которые являются регулярными относительно перехода от "сильно-нелинейного" случая к "слабо-нелинейному"; Доказано, что для этих уравнений фазовый сдвиг в анзаце Кузмака-Уизема определяется из задачи для системы уравнений Уизема с соответствующим образом подобранными начальными условиями;

4. Получено представление в виде точечных преобразований для преобразований, связывающих три одномерные системы: уравнений

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

Научная и практическая ценность. Работа носит теоретический характер. Асимптотические методы решения задач математической физики сами по себе представляют теоретический интерес. Асимптотические решения часто объясняют ключевые свойства точных решений, получаемых численно или экспериментально. Рассматриваемые задачи описывают несколько практически важных явлений, например, поведение жидкости в каналах или динамику волн цунами. Полученное представление, равномерное относительно перехода от "сильнонелинейного" случая к "слабо-нелинейному", позволяет описывать динамику цуга волн, когда на хвосте цуга амплитуда решения небольшая и ситуация "слабо-нелинейная", а у пика - "сильно-нелинейная".

(2012, Москва), на международной конференции "Нелинейный анализ и спектральные задачи" (2013, Уфа) и на международной конференции "Дифференциальные уравнения. Функциональные пространства. Теория приближений", посвященной 105-летию со дня рождения С.Л. Соболева (2013, Новосибирск). Также результаты били представлены и обсуждались на научных семинарах под руководством В.М. Бабича в ПОМИ РАН (2012, Санкт-Петербург) и под руководством Л.А. Каляки-на в Институте Математики с ВЦ УНЦ РАН (2013, Уфа).

Публикации. Основное содержание работы отражено в четырех публикациях, список которых приведен в конце автореферата. Рабо1ы опубликованы в изданиях, входящих в утвержденный ВАК перечень ведущих рецензируемых научных изданий и журналов, в которых должны быть опубликованы основные научные результаты диссертации на соискание ученой степени киндидата и доктора наук.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

Фазовый сдвиг в анзаце Кузмака-Унзема

В первой главе проводится обзор методов осреднения, разработанных Н. М. Крыловым, Н. Н. Боголюбовым, Ю. Л. Митропольским [6,26] и другими, и приводятся результаты, связанные с процедурой определения фазового сдвига в анзаце Кузмака-Унзема при построении быстро-осциллирующих асимптотических решений в нелинейных уравнениях.

Согласно стандартной схеме, примененной Г. Е. Кузмаком (в 1959) [27] для обыкновенных дифференциальных уравнений и Г. Б. Уиземом (1965) [66,67] для уравнений в частных производных, главный член асимптотического разложения определяется из нелинейного уравнения, а поправки к нему из линейных уравнений. Главный член разложения может быть представлен в форме . 1{1)А), где фаза 3, "мед-

ленно меняющийся" параметр или параметры I и так называемый фазовый сдвиг ф определяются из системы "осрсдпснных" уравнений -условий совместности уравнений для поправок. Вопрос о вычислении фазового сдвига возник давно и обсуждался для уравнений в частных производных С. Ю. Доброхотовым и В. П. Масловым (1981) [13], а также Р. Хаберманом (1988) [59]. Для обыкновенных дифференциальных уравнений этот вопрос обсуждался в работах Р. Хабермана и соавт оров.

Основная сложность при определении фазового сдвига состоит в том, что в "сильно-нелинейном" случае для вычисления фазового сдвига достаточно рассмотреть уравнение на первую поправку к главному члену разложения, а в линейном и "слабо-иелинсйном" случае необходимо также исследовать уравнение для второй поправки и для этого необходимо определить первую поправку, что является нетривиальной задачей. Это препятствие имеет топологический характер: при пе-

реходе от "сильно-нелинейного" случая к "слабо-нелинейному" размерность коядра соответствующего линейного оператора в вариациях уменьшается на 1. Практический результат, основанный на идеях работ С.Ю. Доброхотова и В. П. Маслова [13], и работ Р. Хабермана [58,59], заключается в том, что можно учесть фазовый сдвиг ел если рассмотреть его как часть фазы 5, и при этом соответствующим образом изменить начальные данные для уравнения на /. Построенные представления для главного члена асимптотического разложения являются инвариантными относительно перехода от "силыю-нслинсйного" случая к "слабо-нелинейному". Это позволяет, например, описывать динамику цуга волн, когда на хвосте цуга ситуация "слабо-нелинейная", а у пика

- "сильно-нелинейная".

В пункте 1.1 первой главы рассматривается одномерный нелинейный осциллятор с потенциалом У(.г.ет), медленно зависящим от времени т, и малым неконсервативным членом £д(х. х. ет), описывающим трение или (и) "накачку", где с <С 1 - малый параметр. Соответствующий математический объект - задача Коши для неизвестного положения х(т, е) и скорости (момента) р(т. е) — х —

Р= -УХ(.Г,€Т) - £/1(р..г,ст). ¡=р. <(=> V + 1(1.ст) + СП{\, г.гт) = 0, (0.1)

.т|г=0 = х°, р\т~о ~ |г=0 = Р^■ (0.2)

Здесь / = Ух = ~{хЛ). У(хЛ). у - гладкие функции, д(р.х.Ет)

- нечетная функция р. Обычно переменную т называют "быстрым" временем в противоположность "медленному" времени i = ет} Задачу (0.1), (0.2) можно переписать, используя медленное время /, тогда она примет форму так называемой сингулярно возмущенной задачи:

С^Х* (¿Х (1 £

=0, х\,=0 = Л = (0.3)

Мы предполагаем, что по крайней мере локально, потенциал У имеет форму потенциальной ямы на интервале (х/. х,) с минимумом в точ-

'В классической теории осреднения обычно принято обошачать "быстрое" и "медленное" время соответственно через £ и т. Здесь же мы приняли дру! не обозначения, а именно 1С, коюрые используются в следующем разделе в случае уравнений в частых производных.

V(x) 2 г

Xup(r) \lr) I ■

1 0 -r

Рис. 1: The potential V(x,t) and E for various /: / = 0.1, t = 0.2 (left). The example of a fast oscillating solution (right).

КС ОС — тш il € {xi,xr). Потенциал V, функция f так же как и начальные значения х°,р° могут зависеть регулярным образом от некоторого действительного параметра (или параметров) ¡j G [0.р()], например, можно рассмотреть потенциал V = u(t)x'2/2 I- /i\/l(r. t). Vl( r. t) = 0(r3). Параметр р может быть связан с парамефом в простейшем случае fi = е. Обычно, ситуацию, когда V = u{i)r2 ¡2 f sVl(r.f), называют слабонелинейным случаем в противоположность сильнонелинейному случаю, когда ф 0. Хорошим примером слабонелинейного случая служит осцилятор Ван-дер-Поля: х I х — ¿a i (1 — х2) = 0. В качестве примера сильнонелинейного случая мы рассмотрим задачу Коши для физического маятника с медленноменяющейся частотой сo(t) и нелинейным неконсервативным членом, описывающим трение при a(t) < 0 и накачку при a (t) > 0: х + u2(t) мн.г - sa(í) х соьх = 0.

Упомянутое разделение на "слабо-" и "сильнонелинейные" случаи является несколько искусственным и зависит oí начальных условий

т°,р°. Действительно, предполагая, что \

£Х°,р° =

i С

Oil

- АПТ

и х" = = л/Ёр0, после замены переменных х = л/ех.р = \/ер

мы получим слабонелинейное уравнение с V = +

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

Понятно, что обычно точное решение задачи Коши (0.1), (0.2) неизвестно, и можно говорить только об асимптотическом решении при

£ « 1. Интересно, что построение таких асимптотик возможно на больших временах, по крайней мере т ~ В переменных медленного времени t = ет "большйе времена" означают t е [0.0(1)], и желаемые решения являются быстроосциллирующими (см. рис. 1). В теории нелинейных колебаний существуют различные асимптотические методы, которые обычно называются методами осреднения. Существует большое число статей и монографий, посвященных методам осреднения, подробную библиографию можно найти в книгах В. И. Арнольда, В. В. Козлова и А. И. Нейштадта [1,2], H.H. Боголюбова и Ю. А. Мит-ропольского [6], J.D. Cole, J. Kevorkian [46]. Несмотря па то, что методы осреднения основаны на похожих идеях, их реализация различается. Можно условно разбить их на две группы: методы, связанные с "заменой переменных" и методы "прямого расчета". Так метод Крылова-Боголюбова [6], теория KAM [1,2] и им подобные попадут в первую группу, а так называемый метод Кузмака [27,44,46] и его обобщение на случай уравнений в частных производных - метод Уизема [60, 66J или нелинейный метод ВКБ следует отнести ко второй группе.

Заметим, что задачи, рассматриваемые в первой главе относятся к так называемым сингулярно возмущенным, которые допускают асимптотики с составляющей типа пограничного слоя. Исследованию последних посвящены, в частности, работы В.Ф. Бутузова, А. Б. Васильевой, М. И. Вишика, С. А. Ломова, Л. А. Люстериика, Н. Н. Нефедова, А. Н. Тихонова, и других. В отличие от таких асимптотик, в диссертации рассматриваются асимптотики, имеющие быстрые осцилляции.

Задача (0.1), (0.2) решена в статье F. J. Boiiiiand and R. Haberman [44] (см. также книгу J.D. Cole, J. Kevorkian [46] и статью A.M. Ильина [22]), посредством метода Кузмака-Уизема, с использованием некоторых идей, предложенных в работе С.Ю. Доброхотова и В. П. Масло-ва [13]. Окончательный ответ представлен в двух разных формах: для слабонелинейного случая (см. также [22]) и для сильнонелинейпого случая (при этом нет перехода от одного к другому). Основной результат данной главы - единая асимптотическая формула, работающая в

обоих случаях. В этом разделе сформулирована и доказана следующая теорема [50]: Теорема 1

Пусть существуют положитетьиые чист а, Ь, /о, такие что дпя

X е [а, ь], г е [ол0}:

(11) вектор-функция (Х(0,1Л). Р(в.1Л)) - 2тг-периодическая по О и является решением невозмущенной системы с " заморолсенныи " временем t:

р=-Н,(хЛ), х = Н,,. II + У{хЛ)., 9,1 - переменные "действие-угол" для невозмущенной системы

1(1) = -1- Г х/2(Я(/.0 - У(гЛ))с1г. в = П{1Л), Щ1Л) = ///(/. £) - тг/ /

и x±{t) - корни уравнения У(х, t) =

(Í2) Элл £ £ [0, ¿o] существует решение Х(/. s) С [аЛ)}, Ф(/.£") задачи Коши для "осредненных" уравнений

<1Ф clX d Ф г/Х

^ = ^(х.í), ^ = О ^ ^ = 0- ^ = -ср. 0

£(/,£) = [ Х0{0.1Л) ■ С]{Г{0.1Л).Х{0.1Л)Л)С10. 2тг Уо

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

Х|т=0 = /VV) + Ф|т=() ^

где 0°, определяются из условий Х(ф°. 0) = Р(0°. 7°, 0) = р{\

и поправка /[' равна

Г,(Л'(0'./о,О).О)г/0') (10 Ь

Тогда для главного члена асимптотического разложения х(т,е), р{т.е) для задачи Коиш (0.1), (0.2) на временах т € [0./о//>] справедливо следующее представление

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