Н. Батхина ВГИ ВолГУ

1. Моделирование регулярных и хаотических режимов в небесно-механических задачах (на примере модели Хилла) Н. БатхинаВГИ ВолГУ

2. Цель работы • Разработать пакет прикладных программ для комплексного исследования небесно-механических задач, использующий современные методы компьютерного моделирования гамильтоновых динамических систем. • Средствами пакета провести изучение основных регулярных и хаотических структур конкретной небесно-механической модели: плоской круговой модели Хилла.

3. СТРУКТУРА ПРОГРАММНОГО КОМПЛЕКСА Вспомогательные программы и модули ООБД FASTDB или STL Программа визуализации данных GNUPlot Библиотеки низкого уровня Высокоточная арифметика QD и ARPRES Процедуры регуляризации LIBREG Программы верхнего уровня Исследование хаотических режимов: Каскады бифуркаций CASCADE Показатели Ляпунова LYAPUNOV Расщепление сепаратрис SEPSPLIT Построение отображения Пуанкаре POINCAREMAP Построение траекторий ORBIT Матричные и численные алгоритмы, вычисление матриц монодромии NUMMETHODS Численное интегрирование ОДУ TAYLOR Продолжение неподвижных точек отображения DERPAR Продолжение семейств периодических решений CONTPER Алгебраические выражения GiNaC

4. Особенности алгоритмов • Большой объем вычислений и числовых данных • Изменение точности вычислений в процессе работы • Наличие особенностей в уравнениях движения и ударных периодических решений

5. Большой объем вычислений и числовых данных • Необходимо распараллеливание вычислительного процесса • Необходима организация хранения больших объемов числовых данных различной точности и эффективного доступа к ним

6. Изменение точности вычислений в процессе работы • Наличие явлений динамического хаоса (каскадов кратного увеличения периода, расщепления сепаратрис) требует применение библиотек высокоточной арифметики и эффективных алгоритмов численного интегрирования систем ОДУ

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

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

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

10. Периодическое решение однозначно задается вектором вектор невязок схема продолжения Метод Келлера (продолжение по длине дуги)

11. Продолжение семейств периодических решений • Алгоритм продолжения периодического решения является предикторно-корректорным методом продолжения по длине дуги семейства в пространстве параметров с адаптивным выбором шага смещения Δs. • Определяются дополнительные параметры – минимальный и максимальный шаги Δsmin Δsmax агрессивность смены шага a и максимальное число итераций Nmax.

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

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

14. Исследование каскадов кратного увеличения периода Наличие каскадов удвоения периода есть свидетельство резкого усложнения структуры фазового пространства. Каскады характеризуются набором универсальных скейлинговых констант: δ – постоянная Фейгенбаума, α и β – масштабные константы. Для их вычисления необходимо в процессе продолжения периодического решения определять последовательно бифуркационные значения параметра продолжения С

16. Характеристика каскада удвоения периода 1-2-4-8-16-32-64-128-256-512

17. Расщепление сепаратрис

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

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

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

21. Блок-схема работы программы

22. Технологические особенности программного комплекса • Использование высокоточной арифметики • Эффективный алгоритм интегрирования систем ОДУ методом Тейлора • Применение параллельного компьютинга • Использование Open Source проектов для переносимости комплекса

23. Высокоточная арифметика • Библиотека qd, реализующая арифметику вещественных чисел с плавающей точкой с 106 и/или 212 значащими битами (30 или 61 значащая десятичная цифра)в виде классов dd_real и qd_real. • Библиотека arpres, реализующая вычисления на произвольной разрядной сетке (классы mp_real, mp_complex). • Каждая из библиотек реализована на C++, имеет интерфейс для C и FORTRAN. Благодаря перегрузке основных арифметических операций, операций сравнения и вычисления элементарных функций библиотеки могут быть использованы в вычислительных алгоритмах без какого-либо изменения последних. • Доступны в исходных текстах по адресу http://www.nersc.gov/~dhbailey/mpdist/mpdist.html

24. Интегрирования систем ОДУ методом Тейлора • Для получения эффективной процедуры численного интегрирования систем ОДУ используется свободно распространяемая утилита Taylor. Она генерирует по формальному описанию системы набор процедур на языке C, реализующих одношаговый метод Тейлора с автоматическим выбором шага и порядка.

25. Особенности утилиты Taylor • Использование алгоритма автоматического вычисления коэффициентов Тейлора • Поддержка библиотек высокоточной арифметики • Применение адаптивных методов определения шага интегрирования и порядка разложения с учетом локальной погрешности и затрат на вычисление коэффициентов • Доступна в исходных текстах http://www.maia.ub.es/~angel/taylor/ • - положительная константа • –радиус сходимости • локальная погрешность • порядок разложения

26. Сравнение с другими интеграторами

27. Параллельный компьютинг • Для интенсификации вычислений используются кластерные вычисления. Распараллеливание применяется при построении отображения Пуанкаре для большого числа точек и для вычисления производных отображения на сетке. • Кластерные вычисления реализованы средствами MPI (только на платформе Linux)

28. Результаты тестирования параллельного алгоритма на гетерогенном кластере • Приведенный рейтинг узлов в гетерогенном кластере • 1. Pentium IV (Celeron) 1,7 Ghz, 256 Mb ram – 1.00 • 2. Atlhon 2.4 Ghz, 512 Mb ram – 2.64 • 3. Pentium III (Celeron) 900 Mhz, 256 Mb ram – 1.42 • 4. Pentium IV (Celeron) 1,7 Ghz, 256 Mb ram – 1.01

29. Технология Open Source • Комплекс поставляется в исходных текстах и использует только Open Source проекты для сборки и эксплуатации: gcc, g++, gmake, gdb, gnuplot • Комплекс апробирован на задаче Хилла, эксплуатировался на компьютерах с операционной системой Linux (ASPLinux 7.3, 9.2) и Windows (пакет Cygwin) • Модульная структура комплекса позволяет наращивать его функциональность

30. Общая схема работы с комплексом • Пользователь определяет набор вспомогательных процедур: вычисление функции Гамильтона задачи в физических и регулярных переменных, дает описание системы канонических и вариационных уравнений для утилиты Taylor. • С помощью сценариев сборки для утилиты make последовательно генерируются тексты процедур интегрирования систем ОДУ задачи (утилита Taylor) для разных классов высокоточной арифметики. • Затем происходит компиляция исходных текстов программ комплекса, использующих эти процедуры. • Пользователь последовательно использует программы и средства визуализации рядов данных (gnuplot)

31. Применение • Данный комплекс может быть применен к комплексному исследованию небесно-механических задач, автономных гамильтоновых систем • С его помощью были получены новые семейства периодических решений и изучены сценарии перехода к динамическому хаосу в плоской задаче Хилла

32. Уравнения движения в равномерно вращающейся системе координат Интеграл Якоби Гамильтониан Плоская задача Хилла

33. Регуляризация Регуляризация Леви-Чивитта Гамильтониан в регулярных переменных

34. Симметрии Канонические уравнения задачи Хилла инвариантны относительно преобразований

35. Вычисление матрицы монодромии M для симметричных периодических решений • Симметрия • Симметрия • Симметрия • J – симплектическая единица

36. Основные семейства периодических орбит • aи с – либрационные (ляпуновские) орбиты (неусточивые) • f – обратнооборотные спутниковые (хилловские) орбиты (устойчивые) • g, g’ – прямооборотные спутниковые (Хилловские) орбиты

37. Орбиты и индекс устойчивости семейства a

38. Орбиты и индекс устойчивости семейства f

39. Орбиты и индекс устойчивости семействаg

40. Орбиты и индекс устойчивости семействаg’

41. Характеристики семействa, c, f, g, g’

42. Результаты исследования семейств п.р. второго рода • Ветвление семейств от двоякосимметричной орбиты с соизмеримостью p/q, где либо p, либо qчетное, происходит с потерей одной из симметрий • Ветвление семейств от двоякосимметричной орбиты с соизмеримостью p/q, где p и qнечетные, происходит с сохранением симметрий

43. Орбиты семейства f_r(1/9)

44. Индекс устойчивости семейства fr(1/9)

45. Индексы устойчивости семейств с соизмеримостями (1/4), (1/5)

46. Орбиты семейств fr(1/4)Σ1и fr(1/4)Σ2 Устойчивая ветвь Неустойчивая ветвь

47. Орбиты семейства fr(1/5)Σ1Σ2 Неустойчивая ветвь Устойчивая ветвь

48. Характеристики некоторых семейств, ветвящихся от f

49. Орбиты семейств fl(1/4)Σ1и fl(1/4)Σ2

50. Структура сечения Пуанкаре в окрестности семейства g’ Индекс устойчивости семейства g’(1/4)