Качественные методы исследования некоторых задач вихревой динамики тема автореферата и диссертации по физике, 01.04.02 ВАК РФ

С"/ >00-4 / -/.'О о

УДМУРТСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ физический факультет

На правах рукописи УДК 530.1

Лебедев Владимир Геннадьевич

КАЧЕСТВЕННЫЕ МЕТОДЫ ИССЛЕДОВАНИЯ НЕКОТОРЫХ ЗАДАЧ ВИХРЕВОЙ ДИНАМИКИ

01.04.02 — теоретическая физика Диссертация на соискание ученой степени кандидата физико-математических наук

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

Ижевск — 1999

Содержание

Введение 4

1 Движение точечных вихрей на плоскости и сфере.

Компактный случай 10

1.1 Вихревая алгебра на плоскости и сфере........................................12

1.2 Классификация алгебры скобок Ли-Пуассона -,.'...........................19

1.3 Бифуркационный анализ движения вихрей на плоскости....................22

1.4 Симплектические координаты для вихрей на плоскости....................27

1.5 Качественный анализ движения трех вихрей на сфере......................28

1.5.1 Канонические координаты.............................30

1.5.2 Бифуркационный анализ движения....................................37

2 Движение трех вихрей. Некомпактный случай. Проблемы

коллапса и рассеяния 44

2.1 Некомпактные движения на плоскости........................................44

2.2 Движения трех вихрей на сфере при условии некомпактности..............49

2.3 Условие коллапса вихрей на плоскости и сфере ..............................55

2.4 Рассеяние вихрей на плоскости ................................................61

2.5 Случай четырех вихрей с нулевой полной интенсивностью................62

3 Качественный анализ

совместной динамики вихревого пятна и точечных вихрей 64

3.1 Вихревые пятна................................ . 64

3.2 Взаимодействие вихревого пятна с точечным вихрем........................68

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

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

3.5 Неинтегрируемость взаимодействия двух точечных вихрей

с вихревым пятном..............................................................77

Заключение Библиография

Введение

Данная работа посвящена изучению нелинейных динамических систем, описывающих поведение вихревых структур в идеальной жидкости. Уравнения движения изучаемых систем могут быть представлены в гамильтоновой форме со скобкой Пуассона, определяемой некоторой алгеброй Ли (скобка Ли—Пуассона). Такая форма уравнений позволяет наиболее естественно учитывать геометрические и динамические симметрии задачи. Кроме того, она дает наиболее простой и наглядный способ установления аналогий между различными динамическими системами [14].

В первых двух главах диссертации исследуется динамика трех точечных вихрей на плоскости и сфере. Хотя исследование особенностей вихревого движения жидкости началось еще с работы Гельмгольца [18], значительный прогресс в понимании физических явлений, связанных с вихревой динамикой [52, 95, 102, 39], наблюдается лишь в последние десятилетия. Это обусловлено как появлением мощных компьютеров и эффективных численных методов, позволяющих моделировать трехмерные взаимодействия вихрей [81], так и совершенствованием экспериментального оборудования, позволяющего осуществлять более тонкие измерения. Экспериментальное открытие когерентных структур (крупномасштабных вихревых образований в свободных сдвиговых течениях — следах и струях, тонкой зоне течений на поверхностях раздела и в пограничных слоях [73]) заставило в значительной мере переосмыслить возможности классической статистической теории турбулентности и обратиться к детерминированным моделям переноса завихренности и энергии. Понимание природы вихревого движения является чрезвычайно важным для понимания природных процессов в в атмосфере и океане. Процессы отрыва, сопротивления движению и генерация шума в различных технических устройствах не могут быть описаны без использования тех или иных вихревых моделей. Динамика завихренности идеальной несжимаемой жидкости обеспечивает физически содержательные примеры нелинейных гамильтоновых систем бесконечной размерности, особенно интересными в связи с проблемами динамического хаоса [23].

Однако, несмотря на длительную историю развития вопроса, изучение модели точеч-

ных вихрей на плоскости сводилось в основном к численному построению траекторий движения [46, 48, 47, 49, 55, 53, 56, 57]. Движение вихрей на сфере было практически не исследовано вообще до недавнего времени [88, 89]. Актуальным на сегодняшний день становится качественное понимание процессов вихревого движения, хотя бы на таком уровне, который достигнут в динамике твердого тела. В данной работе для анализа движения были использованы такие подходы, как бифуркационный анализ, построение фазовых портретов, геометрических проекций, ставшие возможными после появления работы [62], в которой была найдена алгебраическая структура, определяющая особенности поведения вихрей на плоскости и сфере. Основываясь на результатах этой работы, в первых двух главах диссертации изучены компактные и некомпактные движения вихрей, рассмотрены такие явления как коллапс и рассеяние вихрей.

В последней главе диссертации изучается взаимодействие точечных вихрей с так называемым вихрем Кирхгофа [28]. Хотя вихри Кирхгофа были обнаружены еще в прошлом веке, но их использование как модели реальных вихрей началось сравнительно недавно, после того как подобные конфигурации стали обнаруживать в натурных и численных экспериментах [72, 84, 91]. Модель эллиптических вихрей Кирхгофа является более сложной и более нелинейной моделью, поэтому ее практическое использование связано с компьютерным моделированием. Численные расчеты и сопоставление с экспериментом [85] убедительно показывают, что модель адекватно описывает, например, многие явления связанные с перемешиванием жидкостей. Однако понимания свойств модели вне рамок компьютерного моделирования отсутствовало до настоящего времени. Целью исследования является понимание данной модели с точки зрения качественного анализа ее внутренних свойств и связанных с ней математических структур. Такой качественный анализ проведен до конца для случая взаимодействия одного точечного вихря с вихрем Кирхгофа. Для случая взаимодействия двух точечных вихрей с вихрем Кирхгофа показано, что задача является неинтегрируемой. С этой целью исследовано расщепление сепаратрис для случая ограниченной задачи (задачи о движении частицы жидкости — точечного вихря малой интенсивности, не влияющего на движение остальных вихрей, в поле скоростей, являющемся суперпозицией скоростей, созданных другим точечным вихрем и вихрем Кирхгофа).

Подробная структура диссертации — следующая.

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

Во втором параграфе первой главы рассмотрена алгебраическая классификация движений скобок Ли-Пуассона трех вихрей на плоскости. Показано, что возможность вещественного разложения вихревой алгебры в прямую сумму алгебр зависит от соотношений между интенсивностями вихрей. Всего существует две различных возможности 5о(3) © И и зо(2,1) © К Еще один случай соответствует четырехмерной разрешимой алгебре. В каждом из случаев найден базис образующих алгебры, приводящий алгебру вихрей к каноническому виду, указан симплектический лист.

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

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

Пятый параграф посвящен вопросам качественного изучения вихрей на сфере. Благодаря траекторному изоморфизму задачи вихрей с задачей Лотки-Вольтерра, обнаруженному в [62], качественный анализ движений на сфере может опираться на соответствующий анализ движения вихрей на плоскости. Отличие между случаем плоскости и сферы проявляется лишь в изменении характерной области, допускающей физически приемлемые решения. Таким образом, изоморфизм с задачей Лотки-Вольтерра дает практическую возможность исследования задачи трех вихрей на сфере. Это особенно необходимо потому, что задача вихрей на сфере значительно сложнее соответствую-

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

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

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

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

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

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

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

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

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

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

Во втором параграфе проведен анализ взаимодействия точечного вихря с вихрем Кирхгофа, введены относительные координаты, получена алгебра скобок Ли-Пуассона для относительных переменных. Показано, что в зависимости от соотношения между интенсивностями алгебра скобок сводится к прямой сумме либо зо(2,1)©#, либо е(2)ф Я. Указан канонический базис и найдены центральные функции алгебры.

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

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

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

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

Глава 1

Движение точечных вихрей на плоскости и сфере. Компактный случай

Модель точечных вихрей является простейшей моделью, описывающей вихревое движение в жидкостях и газах, восходящей еще к работам Гельмгольца, Кирхгофа, Пуанкаре [18, 28, 93]. Современные обзоры состояния в этой области могут быть найдены в [15, 31, 39, 50, 52]. По сути дела, модель представляет из себя особые ¿-образные решения уравнений Эйлера [39, 31, 15]. Такие решения возникают в том случае, если размерами вихрей можно пренебречь (например, по сравнению с расстоянием между различными вихрями). В этом случае можно считать, что вся завихренность отдельного вихря Г (интенсивность вихря) сосредоточена в одной точке и наличие этой завихренности приводит к вращательному движению жидкости вокруг вихря. Результирующее поле скоростей в жидкости является суперпозицией поля скоростей от каждого точечного вихря. Взаимодействие вихрей связано с тем обстоятельством, что каждый вихрь движется в поле скоростей, создаваемых остальными вихрями.

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

где Г* — интенсивности точечных вихрей, г — расстояние между вихрями. Центр завихренности при этом движется равномерно и прямолинейно. При условии Г\ = —Г2 центр завих�