Явление буферности и хаос в нелинейных волновых уравнениях тема диссертации и автореферата по ВАК РФ 05.13.18, кандидат физико-математических наук Глызин, Дмитрий Сергеевич

• Глызин, Дмитрий Сергеевич
• кандидат физико-математических наук
• 2006

Оглавление диссертации кандидат физико-математических наук Глызин, Дмитрий Сергеевич

Введение

1. Двухмодовые циклы нелинейного телеграфного уравнения в случае резонанса 1:

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

1.2. Алгоритмическая часть.

1.3. Линейный анализ устойчивости.

1.4. Обоснование алгоритмической части.

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

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

2.2. Алгоритмическая часть.

2.3. Анализ квазинормальной формы.

2.4. Линейный анализ.

2.5. Основной результат.

2.6. Численный анализ квазинормальной формы.

3. Буферность в системе трех однонаправленно связанных телеграфных уравнений

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

3.2. Построение асимптотики цикла системы связанных телеграфных уравнений.

3.3. Анализ устойчивости цикла. Основной результат.

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

4. Программный комплекс ТгасегЗ 58 4.1. Компилятор формул.

4.1.1. Особенности работы с анализатором формул

4.2. ТгасегЗ, общая информация

4.2.1. Построение фазовых портретов.

4.2.2. Вычисление ляпуновских показателей

4.2.3. Метод Бенеттина

4.2.4. Метод динамической перенормировки

4.2.5. Метод динамической перенормировки для дифференциальных уравнений.

4.2.6. Эффективность и надежность метода динамической перенормировки

4.2.7. Зависимость старшего ляпуновского показателя от параметра

4.3. Работа с программой.

4.4. Подробное описание интерфейса.

Введение диссертации (часть автореферата) на тему «Явление буферности и хаос в нелинейных волновых уравнениях»

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

Одним из наиболее содержательных примеров привлечения к моделированию физических явлений нелинейных волновых уравнений является математическая модель RCLG автогенератора с отрезком длинной линии в цепи обратной связи. Эта задача была поставлена и частично решена A.A. Виттом в работе [1], она состоит из линейных уравнений в частных производных гиперболического типа, а нелинейность содержится в краевых условиях. При решении проблемы Витт использовал метод медленно меняющихся фаз и амплитуд и на этом пути получил асимптотические формулы для периодических решений изучаемой задачи. Вид асимптотических формул позволил, кроме того, заключить, что при подходящем выборе параметров данная краевая задача может обладать произвольным наперед заданным числом циклов. Такое явление приобрело название буферности, и под ним в настоящий момент понимают одновременное существование в фазовом пространстве динамической системы сколь угодно большого числа однотипных аттракторов (циклов, торов и т.д.) Следует отметить, что результаты, полученные Виттом, были в значительной степени эвристическими, поскольку отсутствовал необходимый математический аппарат для их обоснования.

В целом, нелинейные краевые задачи гиперболического типа довольно долго изучались лишь в квазилинейной постановке. Такой подход подробно изложен и обоснован, например, в книге Ю.А. Митропольского и Б.И. Мосеенкова [2]. Бифуркационные задачи для уравнений гиперболического типа начали изучаться лишь в конце прошлого века в работах Ю.С. Колесова, А.Ю. Колесова и С.А. Кащенко [3-11].

Основная трудность, с которой приходится сталкиваться в этой ситуации, состоит в том, что при критических значениях параметров счетное число точек спектра оператора линеаризованной задачи лежат на мнимой оси. Тем самым реализуется так называемое бесконечномерное вырождение. Одной из характерных задач, обладающих таким свойством, является телеграфное уравнение: с краевыми условиями Дирихле или Неймана на границах отрезка [0,1] изменения про

0.1) странственной переменной х. Здесь u{t,x) - скалярная функция, определенная при t > t0, О < х < 1, £ - положительный малый параметр, а > 0, f(u,v) — скалярная функция из С00, порядка малости в нуле выше первого.

Как известно, цепочки и решетки связанных генераторов с сосредоточенными параметрами являются полезными физически содержательными моделями, позволяющими выяснить ряд закономерностей развития пространственно-временного хаоса в сплошных средах [12-14]. При этом, как правило, в качестве отдельно взятого звена цепочки (парциальной системы) рассматривается генератор, описывающийся системой обыкновенных дифференциальных уравнений с единственным устойчивым циклом. Например, в работах [12-14] бралась одна и та же парциальная система й — и - d\u\2u, d = 1 + гс0, с0 G К, (0.2) где и — комплекснозначная функция, но рассматривались различные отвечающие ей цепочки. А именно, в [12,13] изучалась цепочка однонаправленно связанных генераторов (0.2), т.е. система вида j + a(uj — Uj-i) = Uj - d\uj\2uj, j = 1,2,., a G €, (0.3) а в [14] — аналогичная цепочка диффузионно связанных генераторов йj = a(uj+1 — 2Uj + Uj-1) + Uj - d\uj\2uj, j = 1,2,., a G C, (0.4) где щ = i = Щ, Rea > 0. Было установлено, что в обоих случаях при достаточно большом числе звеньев в соответствующей системе может наблюдаться хаотическое поведение, обусловленное коллективным взаимодействием парциальных осцилляторов. Предположим теперь, что в цепочках (0,3), (0.4) или в какой-либо аналогичной цепочке каждое звено заменено генератором с распределенными параметрами вида (0.1), в итоге получаем систему, имеющую (при определенных дополнительных условиях) достаточно большое число сосуществующих аттракторов различной природы. Как будет показано ниже при рассмотрении конкретных примеров, для того, чтобы добиться требуемого эффекта, вовсе не обязательно брать цепочку из большого числа звеньев, как это обычно делается в случае сосредоточенных осцилляторов, или, как это было сделано в статье [15], в которой рассматривались цепочки диффузионно связанных обобщенных кубических уравнений Шредингера или нелинейных телеграфных уравнений. Достаточно ограничиться некоторым минимально допустимым их количеством. В работе [16] рассматривается система трех однонаправленно слабо связанных осцилляторов, задаваемых обыкновенными дифференциальными уравнениями.

В данной работе в качестве парциальных осцилляторов были выбраны нелинейные уравнения гиперболического типа. Таким образом, исследовались системы из трех одно-направленно связанных волновых уравнений вида d2Uj дщ 2d2Uj . du*, . 1 л „ + Щ + еащ-1 = о + ^ )' Э = 2'3' щ = с краевыми условиями Дирихле или Неймана.

Цель работы. Основной целью данной диссертационной работы является изучение динамических свойств нелинейных краевых задач (0.1) в различных постановках и с различными нелинейностями f(u,v), а также систем (0.5) связанных в кольцо осцилляторов такого типа, Особое внимание уделяется условиям возникновения явления буферности в этих задачах.

Методы исследования. Методом решения бифуркационных задач с бесконечномерным вырождением является метод квазинормальных форм, впервые предложенный в работе Ю.С. Колесова [3] для уравнений параболического типа с малой диффузией, а затем обоснованный и распространенный на уравнения гиперболического типа и на сингулярно возмущенные дифференциальные уравнения с запаздыванием (A.B. Васильева [4], С. А. Кащенко [5-8], А.Ю. Колесов [9-11], А.Ю. Колесов и Н.Х. Розов [17-19]). Однако, этот метод не позволяет автоматически распространять свойства решений квазинормальной формы на решения исходной задачи и как следствие говорить о их соответствии друг другу. Для обоснования такого соответствия необходимо оценивать невязку асимптотического приближения решения, полученного на основе грубых устойчивых режимов нормализованной системы (см., например, работы [17-19], где соответствующий анализ проделывается).

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

1) Доказано существование и устойчивость двухмодового цикла телеграфного уравнения с малым параметром при квадратичной нелинейности при наличии внутреннего резонанса 1:2.

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

3) Показано наличие буферности в системе трех однонаправленно связанных осцилляторов без квадратичной нелинейности.

4) Создан программный комплекс численного анализа динамических систем ТгасегЗ.

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

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

Апробация работы. Основные результаты работы были представлены на следующих научных конференциях:

1) XXVII Конференция молодых ученых механико-математического факультета Московского государственного университета имени М.В. Ломоносова, Москва, 2005;

2) XXVIII Конференция молодых ученых механико-математического факультета Московского государственного университета имени М.В. Ломоносова, Москва, 2006 (В рамках общеуниверситетской конференции молодых ученых «Ломоносов-2006».);

3) VIII Крымская международная математическая школа «Метод функций Ляпунова и его приложения» (МФЛ-2006), Крым, Алушта, 2006.

Кроме того, результаты диссертации докладывались на ряде семинаров кафедры математического моделирования Ярославского государственного университета им. П.Г. Демидова, а также обсуждались на семинаре «Моделирование и исследование нейронных сетей» кафедры компьютерных сетей Ярославского государственного университета им. П.Г. Демидова.

Публикации. По теме диссертации автором опубликовано 8 работ [20-27]: 5 статей и 3 тезисов докладов. Из работ, выполненных совместно, в диссертацию включены результаты, полученные автором.

Краткое содержание работы. Структурно диссертационная работа состоит из введения, четырех глав, заключения и списка литературы, содержащего 52 наименования. Диссертация содержит 13 рисунков и одно приложение, в котором приводится принципиальная часть разработанных автором вычислительных процедур, входящих в программный комплекс ТгасегЗ.

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Глызин, Дмитрий Сергеевич

Заключение

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

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

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

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

Список литературы диссертационного исследования кандидат физико-математических наук Глызин, Дмитрий Сергеевич, 2006 год

1. Витт, A.A. Распределенные автоколебательные системы / A.A. Витт // Журн. технич. физики. - 1934. - Т.4, № 1. - С. 144 - 157.

2. Митрополъский, Ю.А. Асимптотические решения уравнений в частных производных / Ю.А. Митрополъский, Б.И. Мосеенков. — Киев: Вища школа, 1976.

3. Колесов, Ю.С. Метод квазилинейных форм в задаче об установившихся режимах параболических систем с малой диффузией /Ю.С. Колесов // Укр. матем. журн. — 1987. Т. 39, № 1. - С. 28- 34.

4. Бифуркация автоколебаний нелинейных параболических уравнений с малой диффузией / А.Б. Васильева и др. // Математический сборник. — 1986. — 130(172), №4(8). С. 488 499.

5. Кащенко, С. А. Применение метода нормализации к изучению динамики дифференциально-разностных уравнений с малым множителем при производной /С. А. Кащенко // Диф. уравнения. 1989. - Т.25, № 2. - С. 262 - 270.

6. Кащенко, С. А. О нормализации в окрестности цикла систем параболических уравнений с малой диффузией ¡С. А. Кащенко // Укр. матем. журн. — 1991. — Т. 43, № 9. С. 1155 -1161.

7. Кащенко, С. А. Построение нормализованных систем для исследования динамики гибридных и гиперболических уравнений ¡С. А. Кащенко // Журн. вычисл. матем. и матем. физ. 1994. Т. 34. № 4. С. 564 - 576.

8. Кащенко, С. А. Уравнения Гинзбурга-Ландау — нормальная форма для дифференциал ьно-разностного уравнения второго порядка с большим запаздыванием /С. А. Кащенко // Журн. вычисл. матем. и матем. физ. — 1998. — Т.38, № 3. — С. 457 465.

9. Колесов, А.Ю. Бифуркация автоколебаний сингулярно возмущенного волнового уравнения /А.Ю. Колесов, Ю.С. Колесов // ДАН СССР. 1990. - Т. 315, № 2. С. 281 -284.

10. Колесов, А.Ю. Устойчивость автоколебаний телеграфного уравнения, бифурцирую-щих из состояния равновесия /А.Ю. Колесов // Матем. заметки. — 1992, — Т. 51. — Вып. 2. С. 59 - 65.

11. Колесов, А.Ю. Параметрические колебания решений телеграфного уравнения с умеренно малой диффузией /А.Ю. Колесов // Сиб. мат. журн. — 1992. — Т. 33, N2 6. — С. 79 86.

12. Апищенко, B.C. Сложные колебания в простых системах /B.C. Анищенко. — М.: Наука, 1990.

13. Гапонов-Грехов A.B., Рабинович М.И., Старобипец И.М. // Письма в ЖЭТФ. 1984. Т. 39, № 12. С. 561-564.

14. Нелинейные волны: Структуры и бифуркации / Под ред. A.B. Гапонова-Грехова, М.И. Рабиновича. — М: Наука, 1987.

15. Глызин, С. Д. Хаотическая буферность в цепочках связанных осцилляторов / С.Д. Глызин, А.Ю. Колесов, Н.Х. Розов // Дифференциальные уравнения. — 2005. — Т. 41, № 1. С. 41 - 49.

16. Глызин, С. Д. О явлениях хаоса в кольце из трех однонаправленно связанных генераторов / С.Д. Глызин, А.Ю. Колесов, Н.Х. Розов // Журнал вычислительной математики и математической физики. — 2006. — Т. 46, № 10. — С. 1809 1821.

17. Колесов, А.Ю. Явление буферности в RCLG-автогенераторе: теоретический анализ и результаты эксперимента / А.Ю. Колесов, Н.Х. Розов // Тр. МИАН. — 233. М.:Наука, 2001. С. 153 207.

18. Колесов, А.Ю. Двухчастотные автоволновые процессы в комплексном уравнении Гинзбурга-Ландау / А.Ю. Колесов, Н.Х. Розов // ТМФ. 2003. - Т. 134, № 3. - С. 353-373.

19. Колесов, А. Ю. Инвариантные торы нелинейных волновых уравнений / А. Ю. Колесов, Н.Х. Розов. М., 2004.

20. Метод динамической перенормировки для нахождения максимального ляпуновского показателя хаотического аттрактора /Д. С. Глызин. и др. // Дифференциальные уравнения. 2005. - Т. 41, № 2. - С. 268-273.

21. Глызии, Д. С. Вычисление старшего ляпуновского показателя отображений усовершенствованным методом / Д. С. Глызин // XXVI Конференция молодых ученых механико-математического факультета МГУ им. М.В. Ломоносова. Тезисы докладов. М.: МГУ, 2004. С. 37 - 38.

22. Глызин, Д. С. Пространственно неоднородные циклы одной краевой задачи в критическом случае / Д. С. Глызин // Труды XXVI Конференции молодых ученых механико-математического факультета МГУ им. М.В. Ломоносова. — М.: МГУ, 2005.- С. 34.

23. Д. Гукенхеймер, Ф. Холмс. Нелинейные колебания, динамические системы и бифуркации векторных полей. Москва-Ижевск: Ин-т компьютерных исследований, 2002.29. http://tracer3.narod.ru

24. Ахо, А. Теория синтаксического анализа, перевода и компиляции. Т.1 / А. Ахо, Дою. Ульман. — М.: Наука, 1978.

25. Лебедев, В.Н. Введение в системы программирования / В.Н. Лебедев. — М.: Статистика, 1975.

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

27. Dormand, J. R. A family of embedded Runge-Kutta formulae / J. R. Dormand and B. J. Brince // J. Сотр. Appl. Math. 1980. - Vol. 6. - P. 19 - 26.

28. Бахвалов, H.C. Численные методы : учеб. пособие для физ.-мат. специальностей вузов / Н.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков; под общ. ред. Н.И. Тихонова.- 2-е изд. — М.: Физматлит, 2002.

29. Кузнецов, С. П. Динамический хаос (курс лекций) / С. П. Кузнецов. — М.: Физ.-мат. лит., 2001.

30. Малинецкий, Г. Г. Современные проблемы нелинейной динамики. Изд. 2-е, исправл. и доп. / Г. Г. Малинецкий, А. Б. Потапов. — М.: Едиториал УРСС, 2002.

31. Benettin, G. Kolmogorov entropy and numerical experiments / G. Benettin, L. Galgani, J.-M. Strelcyn // Phys. Rev. 1976. - V. A14. - P. 2338 - 2345.

32. Демидович, В. П. Лекции по математической теории устойчивости / Б. П. Демидо-вич. — М.: Наука, 1967.

33. Оселедец, В. И. Мультипликативная эргодическая теорема. Характеристические показатели Ляпунова динамических систем / В.И. Оселедец. Тр. Моск. мат. об-ва. — 1968. Т. 19 - С. 179 - 210.

34. Milnor, J. On the concept of attractor: Correction and remarks //. Milnor // Commun. Math. Phys. 1985. V.99, № 2. - P. 177 - 196.

35. Devaney, R. An introduction to chaotic dynamical systems / R. Devaney. — Addison-Wesley: Reading, MA, 1989.

36. On the Devaney's definition of chaos / J. Banks and others. // Amer. Math. Monthly.- 1992. V.99, № 4. - P. 332 - 334.

37. Афраймович, B.C. О притягивающих негрубых предельных множествах типа аттрактора Лоренца / B.C. Афраймович, В.В. Быков, Л.П. Шилъников // Труды Московского мат. общества. — 1982. Т. 44. — С. 150 212.

38. Плыкин, Р.В. Странные аттракторы / Р.В. Плыкин, Е.А. Сатаев, С.В. Шлячков // Динамические системы с гиперболическим поведением. Итоги науки и техники. Совр. проблемы математики. М.:ВИНИТИ, 1991. - Т.66. - Гл. 2. - С. 100 - 148.

39. Robinson, С. Homoclinic bifurcation to a transitive attractor of Lorenz type // Nonlinearity.- 1989. V.2, № 4. P. 495 - 518.

40. Синай, Я.Г. Стохастичность динамических систем / Я.Г. Синай. // Нелинейные волны. М.: Наука, 1979. - С. 192 - 212.

41. Бланк, М.Л. Устойчивость и локализация в хаотической динамике / М.Л. Бланк. — М.: МЦНМО, 2001.

42. Wolf, A. Determining Lyapunov exponents from a time series / A. Wolf, J.B. Swift, H.L. Swinney, J.A. Vastano. Physica D. 1985. - V. 16. - C. 285 - 317.

43. Глызин, С. Д. Динамические свойства простейших конечноразностных аппроксимаций краевой задачи «реакция-диффузия» / С. Д. Глызин // Дифференциальные уравнения. 1997. - Т.ЗЗ, № 6. - С. 805 - 811.

44. Рубаник, В.Н. Колебания квазилинейных систем с запаздыванием / В. Н. Рубаник.- М.: Наука, 1969.

45. Колесов, Ю. С. Проблема адекватности экологических уравнений / Ю. С. Колесов.- Ярославль, 1985. Деп. в ВИНИТИ 1985, №1901-85.