Интегрируемость и стохастичность в консервативных динамических системах тема диссертации и автореферата по ВАК РФ 01.04.02, кандидат физико-математических наук Симаков, Николай Николаевич

• Симаков, Николай Николаевич
• кандидат физико-математических наук
• 1999

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

Оглавление

1. Интегрируемые задачи вихревой динамики. 9 1.1. Движение двух точечных вихрей внутри цилиндрической области

1.1.1. Алгебраическая классификация

1.1.2. Бифуркационный анализ

1.1.3. Симплектические координаты

1.1.4. Бифуркационные диаграммы задачи

1.1.5. Возможность динамического коллапса

2. Неинтегрируемые задачи вихревой динамики

2.1. Движение двух точечных вихрей вне круговой области во внешнем стационарном течении

2.1.1. Общие уравнения движения

2.1.2. Движение двух точечных вихрей в набегающем с постоянной скоростью потоке

2.1.3. Стационарные конфигурации

2.2. Неинтегрируемость задачи о движении трех точечных вихрей в цилиндрической области

2.3. Стационарные конфигурации в задаче о движении трех вихрей внутри цилиндрической области

3. Неинтегрируемость и сценарии развития стохастичности в динамике твердого тела с одной неподвижной точкой

3.1. Бифуркации удвоения периода в динамике твердого тела

3.1.1. Состояние проблемы хаотизации периодических процессов в динамических системах через последовательность бифуркаций кратного увеличения периода

3.1.2. Методика построения последовательности бифуркаций удвоения периода периодических решений уравнения Эйлера-Пуассона

3.1.3. Результаты численных экспериментов

3.2. Численное доказательство неинтегрируемости уравнений Эйлера-

Пуассона

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

3.2.2. Построение сепаратрис отображения Пуанкаре

3.2.3. Результаты численного эксперимента

7 о

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

Введение диссертации (часть автореферата) на тему «Интегрируемость и стохастичность в консервативных динамических системах»

Введение.

Диссертационная работа посвящена исследованию динамических эффектов, наблюдающихся в консервативных динамических системах из различных разделов математической физики. Их исследование имеет значение, как для теории так и для приложений.

Первые две главы диссертации посвящены исследованию задач вихревой динамики. Теория вихревого движения имеет давнюю историю. На ее основе Декарт пытался объяснить движение планет, построить модель Вселенной [2]. Математическое описание процессов, связанных с движением завихренности в жидкости было предложено в 1858 г. Г. Гельмгольцем [44]. Первые исследования были в основном направлены на создание объясняющей инерцию и гравитацию "вихревой теории материи" (эфиродинамика, атом Кельвина и т.д.). В настоящее время интерес к вихревой динамике связан с изучением природных процессов в атмосфере и в океане, с исследованием процессов отрыва потока, сопротивления движению тел в жидкости и т.д. Задачи изучения движения точечных вихрей, в областях с твердыми неподвижными границами, возникли на раннем этапе исследования вихревых структур. Простейший случай движения одного точечного вихря в идеальной жидкости, ограниченной одной плоскостью, рассмотрен Г. Гельмгольцем [47]. При этом используется метод зеркальных отображений, в котором влияние границы заменяется взаимодействием с вихрем противоположной по знаку интенсивности, зеркально симметрично расположенным к исходному, что позволяет обеспечить выполнение граничных условий. Такой способ инверсии широко используется при решении задач гидродинамики.

Общая теория движения вихрей в произвольной области была заложена в работах Э. Рауса [45] и детально разработана Ц.Ц. Линем [46]. Согласно этой теории, движение N точечных вихрей с интенсивностями к\. ..к^ в односвя-занной области описывается гамильтоновой системой с N степенями свободы, в которой роль сопряженных переменных выполняют декартовые координаты особенностей.

1 дН . 1 дН .

К* ОУг кг ОХг

Интегрируемость этих систем зависит от числа вихрей и от вида области, в которой происходит движение.

На плоскости движение вихрей интегрируемо, если число вихрей не превышает трех. Задача о движении двух точечных вихрей на плоскости была полностью изучена Г. Гельмгольцем [34], который установил, что в общем случае два вихря совершают равномерное вращательное движение вокруг центра

завихренности с координатами

N N

12 кхХх kiУi

Х = Ц-= (0.2)

12 ^г 12 к{

¿=1 1=1

и частотой

п = ^ («■»)

где г расстояние между вихрями, являющееся интегралом движения. В случае вихревой пары центр завихренности находится на бесконечности, а два вихря движутся поступательно. Задача о движении трех вихрей значительно сложнее. Бе интегрируемость была отмечена Гребли [35] и Пуанкаре [36]. В работе [35] содержится анализ различных траекторий системы, в [36] указан полный набор некоммутативных интегралов. Классификацию пуассоновых структур, а так же бифуркационный анализ движения вихрей можно найти в работе [37]. Система четырех вихрей на плоскости уже неинтегрируема. Неинтегрируемость ограниченной постановки задачи была доказана в [38]. В работе [56] особое внимание уделено развитию методов численного анализа рассматриваемых систем, что позволяет получить результаты о поведении системы, недоступные для современных аналитических методов.

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

Задача о движение двух точечных вихрей в цилиндрической области интегрируема и изучалась в работах [40, 41]. В первой главе диссертации дается новый более полный анализ этой задачи, основанный на качественных методах исследования приведенной системы [37]. Излагаемый подход опирается на представлении уравнений движения на алгебре Ли и изучении регуляри-зованной системы. Данный метод позволяет исследовать проблему взаимного слияния (коллапса) вихрей, и дает более ясное геометрическое представление о структуре приведенного фазового пространства задачи.

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

В работе так же исследовано движение вихрей за цилиндром при наличии набегающего потока. Эта модель впервые была рассмотрена Фепплем [55] и ей уделялось значительное внимание, в виду ее математической простоты и связи с моделями вихревых дорожек. Наиболее простой является задача о движении

симметричной вихревой пары за цилиндром. Линии тока, такой модели хорошо согласуются с экспериментом для течений с числом Рейнольдса В.е = ^ в интервале 5 < Ее < 30.

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

Третья глава диссертации посвящена качественному исследованию динамических эффектов, сопровождающих развитие хаоса в неинтегрируемых гамиль-тоновых системах. На сегодняшний день известен ряд сценариев развития критических явлений (перехода к стохастичности). Эти сценарии не зависят от конкретного вида динамической системы и в этом смысле являются универсальными.

Одним из наиболее известных сценариев развития хаоса в системе является сценарий Фейгенбаума. Впервые универсальность Фейгенбаума была обнаружена в классе одномерных диссипативных отображений отрезка на себя [5]. Особенный интерес представляет наличие универсальности в классе двухмерных отображений сохраняющих площадь, которые представляют поведение гамиль-тоновых систем с двумя степенями свободы. Этот вопрос достаточно полно изучен на примере простейших модельных отображений в работах [7, 8, И]. Особый интерес вызывает изучение последовательности бифуркаций удвоения в реальных динамических системах, для которых фактически невозможно применение одних лишь аналитических методов исследования.

В третьей главе работы, рассмотрена классическая задача о движении тяжелого твердого тела с одной неподвижной точкой, которое описывается уравнениями Эйлера-Пуассона

М = М х АМ - д (г х 7), 7 = 7 х АМ, М, 7 е И3, (0.4)

где А - обратный тензор инерции; М, 7 - соответственно векторы кинетического момента и единичного орта вертикали в системе координат жестко связанной с телом; г = (х, у, г) - радиус-вектор центра масс в этой же системе координат; ¡л - вес тела. В работе исследован сценарий перехода к хаосу через последовательность бифуркаций удвоения при изменении параметров системы и полной энергии. Показано, что развитие локальной неустойчивости в системе и увеличение стохастического слоя сопровождается бесконечной последовательностью бифуркаций удвоения, обладающей универсальным асимптотическим поведением.

Уравнения Эйлера-Пуассона (3.7) могут быть сведены на уровне геометрического интеграла С = (т>7) и интеграла площадей К = (М, 7) к двухстепенной гамильтоновой системе с гамильтонианом

Н = 1/2 (АМ,М) -/х(г,7). (0.5)

Для таких систем возможна визуализация движения при помощи сечения Пуанкаре.

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

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

Проблема интегрируемости и неинтегрируемости в динамике твердого тела занимает одно из центральных мест. Особенно много усилий было потрачено на поиски интегрируемых случаев. Многое было сделано русскими учеными (C.B. Ковалевская, С.А. Чаплыгин, A.M. Ляпунов, Н.Е. Жуковский, В.А. Стеклов). Впоследствии в работах Пуанкаре, Кирхгофа было показано, что интегрируемые случаи являются изолированными в пространстве параметров. Строгое доказательство неинтегрируемости в при различных значениях параметров системы (0.4) содержится в работах [25, 27, 28, 30, 33]. В третьей главе диссертации приведено численное "доказательство" неинтегрируемости уравнений Эйлера-Пуассона, для тех значений параметров, при которых не известно аналитического доказательства. Это доказательство основано на явлении расщепления сепаратрис (поверхностей заполненных траекториями, асимптотически приближающимися к гиперболическим периодическим решениям при t ±оо). Расщепление сепаратрис является одним из динамических эффектов, препятствующих интегрируемости уравнений движения динамических систем и приводящим к стохастичности. В интегрируемых системах сепаратрисы, исходящие из двух различных гиперболических периодических решений, сдвоены или вообще не имеют общих точек. В общем случае эти сепаратрисы расщепляются и можно показать, что в случае трансверсального пересечения сепаратрис система не имеет полного набора аналитических интегралов [3].

В аналитическом подходе исследование расщепления сепаратрис сводится к вычислению интеграла Пуанкаре-Мельникова, который возникает при анализе величины расщепления в первом порядке теории возмущения вблизи интегрируемого случая. В некоторых случаях его удается вычислить с помощью вычетов [4]. Более сложной с аналитической точки зрения является ситуация, когда расщепление сепаратрис имеет более высокий порядок по малому параметру или экспоненциально мало. В этом случае для анализа существует мало конструктивных критериев. Численные методы позволяют изучать поведение

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

Заключение диссертации по теме «Теоретическая физика», Симаков, Николай Николаевич

Основные результаты диссертационной работы состоят в следующем:

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

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

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

4) Методом расщепления сепаратрис доказана неинтегрируемость ограниченной постановки задачи о движении трех вихрей внутри цилиндрической области.

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

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

7) Методом численного анализа расщепления сепаратрис доказана неинтегрируемость уравнений Эйлера-Пуассона при значениях параметров отличных от известных интегрируемых случаев.

Период Параметр Координаты неподвижной точки S,

1 26.6488874 (3.876717147123, -0.000000000215)

2 31.3958859 (3.504207723971, -0.000000000043)

4 32.2921660 (3.398114068641, -0.000000000040) 5.296333702

8 32.4047973 (3.572686680801, -0.000000000101) 7.957646764

16 32.4178867 (3.566555047962,-0.000000000048) 8.604771800

32 32.4193897 (3.580298907896, -0.000000000019) 8.708848969

64 32.4195621 (3.577633696195, -0.000000000049) 8.718097448

I = diag{2,3,1), г = (1,1,0),Я = 5, Energy = 50,g{mod2тг) = тг/2,д > 0.

Заключение.

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

Литература

[1] Ю.А. Архангельский. Аналитическая динамика твердого тела. М., Наука, 1978.

[2] В. В. Козлов. Общая теория вихрей. Ижевск., Издательский дом "Удмуртский университет", 1998, 238 С.

[3] С. Л. Зиглин. Расщепление сепаратрис и несуществование первых интегралов в системах типа гамильтоновых с двумя степенями свободы. АН СССР, сер. Матем. 1987, т.51 с.1088-1103.

[4] В. В. Козлов. Методы качественного анализа в динамике твердого тела. М., МГУ, 1980.

[5] М. J. Feigenbaum. Quntative universality for class of non-linear transformation // J. Stat. Phys. 1978. 19. P. 25-52.

[6] A. J. Lichtenberg, M.A. Lieberman. Regular and stohastic motion. Springer. 1983.

[7] P. Collet, J.-P. Eckman, H. Koch. On universality for area preserving maps of the plane. //Physica 3D. 1981. P.457-467.

[8] J. M. Green, R. S. MacKay, F. Vivaldi. M. J. Feigenbaum. Universal behavior in families of area - preseving maps. //Physica 3D. 1981. P. 468 - 486.

[9] Т. C. Bountis. Period doubling bifurcation and universality in conservative systems. //Physica 3D. 1981. P.577-589

[10] S. P. Kuznetsov, I. R. Sataev. Period-doubling for tow-dimensional non-invertible maps: Renormalization group analysis and quantative universality. Physica D. 101. 1997. P. 249-269.

[11] R. S. MacKay. Renormalisation in area-preserving maps. Wold Scientific. 1993.

[12] Я. Г. Синай. Современные проблемы эргодической теории. М.: Физматлит. 1995.

[13] J.-M.Mao, 1.1. Satija, B.Hu. Period doubling in four-dimensional sympletic maps. Physical Review A. 1986. V.34, N5. P.4325-4332.

[14] Jian-Min Mao, Robert H., G. He.lle.man. New Feigenbaum constant for four-dimensional volume-preserving sympletic maps. Physical Review A. 1987. V.35, N4. P.1847-1855.

[15] P. Collet, J.-P. Eckman, H. Koch . Period doubling bifurcation for families of maps in Rn J.Stat.Phys. 1981. V.25, N1. P.l-14.

[16] В. Ни, 1.1. Satija . A spectrum of universality classes in period doubling and period tripling. Phys.Lett. 1983. V.98A, N4. P. 143-146.

[17] А. И. Гольдберг, Я. Г. Синай, К. М. Ханин. Универсальные свойства для последовательности бифуркаций утроения. //УМН. 1983. Т.38, N1. С.159-160.

[18] В. В. Козлов. Методы качественного анализа в динамике твердого тела. М.,МГУ 1980.

[19] В. И. Арнольд. Математические методы классической механики. М.: Наука, Гл. ред. физ.-мат. лит., 1974. С.432.

[20] В. Н. Рубановский, В. А. Самсонов. Устойчивость стационарных движений в примерах и задачах. М.: Наука, Гл. ред. физ.-мат. лит., 1988. С.105-123.

[21] В. И. Гуляев, А. Л. Зубрицкая, В. Л. Кошкин . Построение последовательностей бифуркаций кратного увеличения периода решений нелинейных обыкновенных дифференциальных уравнений. //Дифф. уравнения. 1989. Т. 25. N6. С.929-933.

[22] Н. В. Емельянов . Методы составления алгоритмов и программ в задачах небесной механики. М.: Наука, Гл. ред. физ.-мат. лит., 1983. С.128.

[23] В. В. Козлов. Несуществование дополнительного аналитического интеграла в задаче о движении несимметричного тяжелого твердого тела вокруг неподвижной точки. //Вест. МГУ Сер. мат., мех. 1975. N1. С.105-110.

[24] В. В. Козлов. Расщепление сепаратрис возмущенной задачи Эйлера-Пуансо. //Вест. МГУ. Сер. мат., мех. 1976. N6. С.99-104.

[25] С. Л. Зиглин. Расщепление сепаратрис, ветвление решений и несуществование интеграла в динамике твердого тела. //Труды Моск. мат. об-ва. 1980. Т.41. С.287-303.

[26] Т. В. Сальникова. Неинтегрируемость возмущенной задачи Лагранжа. // Вестн. МГУ. Сер. мат., мех. 1984. N4. С.62-66.

[27] С. В. Болотин. Двоякоасимптотические траектории и условия интегрируемости гамильтоновых систем. // Вестн. МГУ. Сер. мат., мех. 1990. N1. С.55-63.

[28] С. А. Довбыш. Расщепление сепаратрис неустойивых равномерных вращений и неинтегрируемость возмущенной задачи Лагранжа. // Вестн. МГУ. Сер. мат., мех. 1990. N3. С.70-77.

[29] В. В. Козлов . О группах симметрии динамических систем. // ПММ.1998. Т.52, вып.4. С.531-541.

[30] С. В. Болотин . Вариационные принципы построенния хаотических движений в динамике твердого тела. // ПММ.1992. Т.56, вып.2. С.230-239.

[31] С. JI. Зиглин. Ветвление решений и несуществование первых интегралов в гамильтоновой механике. // Функц. анализ и его прил. 1982. Т.16, N3. С.30-41. 1983. Т.17, N1. С.8-23.

[32] С. JI. Зиглин. Об отсутствии вещественно-аналитического первого интеграла в некоторых задачах динамики. // Функц. анализ и его прил. 1997. Т.31, вып. 1. С.3-11.

[33] В. В. Козлов, Д. В. Трещев. Неинтегрируемость общей задачи о вращении динамически симметричного тяжелого твердого тела с неподвижной точкой. // Вестн. МГУ. Сер. мат., мех. 1985. N6. С.73-81. 1986. N1. С.39-44.

[34] В. В. Мелешко, М. Ю. Константинов. Динамика вихревых структур. Киев, Наукова думка, 1993, 277С.

[35] W. Gröbli. Specialle Probleme über die Bewegung geredliniger paralleler Wirbelfäden. Vierteljahrsch. d. Naturforsch. Geselsch. Zürich, v.22, 1877, s.37-81, 129-165.

[36] H. Poincaré. Théorie des tourbillous. Paris, 1893, 205p.

[37] A. V. Borisov, V. G. Lebe.de.v. Dynamics of three vortices on a plane and a sphere - II. //Reg. & Chaot. Dyn., 1998, v.3, N2.

[38] C. JI. Зиглин. Неинтегрируемость задачи о движении четырех точечных вихрей. //ДАН СССР. N9 1979. с.1296-1300.

[39] L. Zannetti, P. Franze.ce.. The non-integrability of the restricted problem of tow vortices in closed domain. //Phisica D. 1994, v.76, p.99-109.

[40] Y. Kimura. Motion of two point vortices in circular domain. //Ibid. 1998, 57, N5. p.1641-1649

[41] J. С. Hardin, J. P. Mason. Periodic motion of two and four vortices in a cylinder pipe. //Phys.Fluids v.27, N7, 1984, P.1583-1589

[42] A. V. Borisov, A. E. Pavlov. Dynamics and statics of vortices on a plane and a sphere - I. //Reg. & Chaot. Dyn., 1998, v.3, N1, p.28-38.

[43] Б. А. Дубровин, С. П. Новиков, А. Т. Фоменко. Современная геометрия: Методы и приложения. Том I. М.: Эдиториал УРСС. 1998.

[44] Г. Гелъмголъц. Два исследования по гидродинамике. М., 1902, С 5-51.

[45] Е. J. Routh. Some application of conjugate function. //Proc. Lond. Math. Soc.-1881.-12, N 170/171. P.73-89.

[46] C. C.Lin. On the motion of vortices in tow dimensions. //Proc. Nat. Acad. Sci. USA. -1941. -27, N 12. P.570-577.

[47] Г. Билля. Теория вихрей. Ленинград,ОНТИ,1936, 266 с.

[48] Г. Биркгоф Гидродинамика. Методы. Факты. Подобие. М.: Изд-во иностр. лит., 1963

[49] В. В. Козлов. Симметрия топология и резонансы в гамильтоновой механике. Ижевск, Изд-во УдГУ, 1995, 432с.

[50] К. М. Khanin Quasi-periodic motion of vortex systems . Physica 4D, 1982, 261-269 p.

[51] H. Poincare.. Les Methodes Nouvelles de la Mechanique Celeste, vols. 1-3 // Gauthier-Villars, Paris, 1892; (Dover, New York, 1957), reprint.

[52] Т. M. Danby. The evolution of periodic orbits close to hetroclinic points. //Cel. Mech. 1984. v.33. N3. P.261-270.

[53] А. П. Маркеев Теоретическая механика. М.:Наука. Гл.ред. физ.-мат. лит., 1990. 416 с.

[54] К. Магнус. Гироскоп. Теория и применение. М: Изд-во "Мир", 1974, 526 С.

[55] Г. Ламб. Гидродинамика. М: ОГИЗ-Гостехиздат, 1947, 280 С.

[56] Я. Aref, N. Pomphrey. Integrable and chaotic motions of four vortixes. Proc. R.Soc., Lond., A380, 1982, P.359-387.