Актуальность темы
В настоящее время достижения в области математического моделирования и вычислительного эксперимента как информационной технологии получений новых знаний об окружающем нас мире приобретают большое значение для различных областей наук, в том числе и для современной астрономии.
В связи с возросшим объемом информации об элементах орбит малых тел, принадлежащих Солнечной системе, возрос интерес к проблеме астероидной опасности. Разработка математической модели и программного обеспечения для численных теорий движения малых тел Солнечной системы является одним из составных этапов при решении этой проблемы. Наибольшую опасность для Земли, наряду с короткопериодическими кометами и крупными фрагментами в метеорных потоках, представляют астероиды (в частности астероиды групп Аполлона, Амура и Атона). Исследование их происхождения, устойчивости движения, оценка вероятности столкновения и предотвращение катастрофических последствий является лишь неполным перечнем проблем, требующих решения. Точность получаемых прогнозов диктуется адекватностью математической модели и адекватностью получаемых на её основе выводов и зависит от ряда факторов: учета в физической модели полного спектра действующих силвыбора численного метода и его сходимостиустойчивости и погрешности самой математической модели.
Изучение движения малых тел Солнечной системы на основе решения стандартной задачи и-тел с учетом лишь гравитационных эффектов, не позволяют получать долгосрочные адекватные прогнозы. Учет релятивистских эффектов необходим для создания высокоточных численных теорий движения небесных объектов.
Бурное развитие численных методов решения обыкновенных дифференциальных уравнений движения небесных объектов, произошло с созданием вычислительных машин. Д. К. Куликовым [58] были построены алгоритмы Коуэлла 8−10 порядков. Ш. Когеном и Э. Хаббартом [111] алгоритмы типа Адамса-Мултона-Коуэлла до 16 порядка. В. Ф. Мячин и О. А. Сизова [66], а также А. Ф. Заусаев [48, 49] разработали метод тейлоровских разложений для задачи я-тел. Гибридные и неявные методы разработали М. С. Яров-Яровой [99, 100], Э. Эверхарт [114] и другие.
Недостатками методов Коуэлла, Адамса-Мултона, гибридных методов является то, что они основаны на использовании разностных схем, что отражается на точности получаемых результатов. Методы, разработанные В. Ф. Мячиным, О. А. Сизовой и А. Ф. Заусаевым ориентированны на решение задачи «-тел и не допускают введения в правые части дифференциальных уравнений движения небесных объектов действия дополнительных возмущений.
Учитывая потребность в повышении точности применяемых методов, а также необходимость учета дополнительных действующих сил, устранение выше перечисленных недостатков является актуальной задачей, определяющей направление исследования диссертационной работы.
Цель работы.
1. Разработка универсального метода решения обыкновенных дифференциальных уравнений на основе метода тейлоровских разложений, позволяющего решать обыкновенные дифференциальные уравнения, правые части которых являются степенными функциями.
2. Создание программного обеспечения для реализации универсального метода на основе метода тейлоровских разложений.
3. Разработка математической модели движения небесных тел (больших планет, астероидов) с учетом в дифференциальных уравнениях как гравитационных сил, так и релятивистских эффектов.
4. Исследование сходимости, устойчивости, погрешности математической модели и обоснование оптимального выбора параметров численного эксперимента.
5. Составление на основе разработанного метода объектно-ориентированных программ численного интегрирования уравнений движения и исследование движения астероидов, сближающихся с Землей-, выявление наиболее опасных объектов среди астероидов, представляющих потенциальную угрозу для Земли.
Научная новизна.
1. Модифицирован метод численного интегрирования дифференциальных уравнений на основе метода тейлоровских разложений за счет получения универсальных формул нахождения точных производных любого порядка от правых частей, содержащих степенные функции.
2. Проведено исследование сходимости, устойчивости и погрешности метода Тейлора на примере решения модельных задач, описывающих движение различных небесных тел (больших планет, астероидов) для различных порядков аппроксимации и шагов интегрирования.
3. Разработан комплекс программного обеспечения для реализации метода тейлоровских разложений и его применения при математическом моделировании движения небесных объектов.
4. На основе математического и программного обеспечения проведено исследование движения 1230 астероидов (на примере групп Аполлона, Амура, Атона) на длительные интервалы времени, выполнена их систематизация и классификация, выявлены 49 потенциально опасных для Земли астероидов и создан каталог их орбитальной эволюции.
На защиту выносятся следующие положения.
1. Модификация метод численного интегрирования дифференциальных уравнений на основе метода тейлоровских разложений за счет получения универсальных формул нахождения точных производных любого порядка от правых частей, содержащих степенные функции.
2. Метод вычислительного эксперимента при анализе математической модели движения небесных тел с использованием банка данных координат и скоростей планет, позволяющий на порядок понизить размерность системы уравнений и повысить вычислительную эффективность модели.
3. Алгоритмы и комплекс проблемно-ориентированных программ для решения уравнений движения небесных тел на основе метода тейлоровских разложений.
4. Новые качественные и количественные результаты исследования орбитальной эволюции малых тел (на примере астероидов групп Аполлона, Амура и Атона), сближающихся с большими планетами и их классификация и систематизация.
Практическая значимость работы. Работа имеет теоретический характер. Разработанный метод может применяться для решения любых обыкновенных дифференциальных уравнений, содержащих в правых частях степенные функции. Разработано математическое и программное обеспечение, имеющее универсальный характер, для исследования малых тел Солнечной системы: комет, крупных фрагментов в метеорных потоках и т. д. Результатом исследования явилось создание каталога орбитальной эволюции 49 потенциально опасных для Земли астероидов. Результаты прогноза относительно потенциально опасных астероидов могут быть использованы при решении проблем, связанных с астероидной опасностью.
Обоснованность выносимых на защиту научных положений, выводов и рекомендаций, а также достоверность полученных результатов исследований проверена путем всестороннего исследования численного метода, основанного на представлении решения в виде тейлоровских разложений, на устойчивость, сходимость, а также согласованностью координат и скоростей больших планет на различные моменты времени с результатами, полученными высокоточными численными методами другими авторами.
Связь диссертационной работы с планами научных исследований.
Работа выполнялась в рамках плана НИР Самарского государственного технического университета (тема «Математическое моделирование движения небесных объектов, разработка высокоточных численных методов интегрирования уравнений движения небесных тел и их программного обеспечения»).
Апробация работы. Основные результаты диссертации докладывались на третьей, четвертой и пятой Международных конференциях молодых" ученых «Актуальные проблемы современной науки» (Самара, 2002, 2003, 2004 гг.), на 8 двенадцатой Межвузовской конференции «Математическое моделирование и краевые задачи» (Самара, 2002), на Всероссийской научной конференции «Математическое моделирование и краевые задачи» (Самара, 2004), на Всероссийской астрономической конференции ВАК 2004 «Горизонты Вселенной» (Москва, 2004), на Международной научной конференции «Актуальные проблемы математики и механики» (Казань, 2004), на шестой Международной Петрозаводской конференции «Вероятностные методы в дискретной математике» (Кисловодск, 2004), на пятом Всероссийском симпозиуме по прикладной и промышленной математике (Сочи, 2004), на шестом всероссийском симпозиуме по прикладной и промышленной математике (Санкт-Петербург, 2005), на Международном симпозиуме «Астрономия — 2005: состояние и перспективы развития» (Москва, 2005), на Второй Всероссийской научной конференции «Математическое моделирование и краевые задачи» (Самара, 2005), на Всероссийской конференции «Дифференциальные уравнения и их приложения» (Самара, 2005), на первом Международном форуме молодых ученых «Актуальные проблемы современной науки» (Самара, 2005), на научных семинарах «Механика и прикладная математика» Самарского государственного технического университета (рук. проф. Радченко В. П., 2003;2005 гг.).
Публикации. По теме диссертации опубликовано 17 работ, список которых приведен в конце автореферата.
Личный вклад автора. Имеется одна работа, выполненная в соавторстве, в которой автору в равной степени принадлежит как постановка задачи, так и результат выполненного исследования. Все остальные работы выполнены автором самостоятельно.
Структура работы. Диссертация состоит из введения, четырех глав, общих выводов, списка литературы и трех приложений, в которых приведен каталог эволюции 49 потенциально опасных астероидов и листинги разработанных программ. Общий объем диссертации 218 страниц, основной текст изложен на 144 страницах. Диссертация содержит 35 рисунков и 11 таблиц.
Список литературы
включает 127 наименований.
Основные результаты выполненных в диссертационной работе исследований состоят в следующем:
1. Модифицирован метод численного интегрирования дифференциальных уравнений на основе метода тейлоровских разложений, за счет получения универсальных формул нахождения точных производных любого порядка от правых частей, содержащих степенные функции.
2. Разработано математическое и программное обеспечение для исследования эволюции орбит больших планет и малых тел Солнечной системы с учетом, как гравитационных сил, так и релятивистских эффектов на основе представления решения отрезком ряда Тейлора.
3. Проведено исследование сходимости, устойчивости и погрешности метода Тейлора на примере решения модельных задач, описывающих движение различных небесных тел (больших планет, астероидов) для различных порядков аппроксимации и шагов интегрирования.
4. На основе численного эксперимента:
— исследована эффективность применения метода разложения в ряд Тейлора для решения задачи п-тел;
— найдены эффективные параметры использования предложенного высокоточного метода Тейлора.
5. На основе тейлоровских разложений разработана новая программно-информационная среда, в рамках которой создан банк данных координат и скоростей больших планет (Меркурий-Плутон) на интервале времени с 3030 года до н. э. по 6970 год н. э., при этом:
— выполнена проверка адекватности полученного банка данных с результатами исследования, проведенных методом Эверхарта и данными теории DE405;
— получены формулы вычисления координат и скоростей больших планет при использовании банка данных путем уточнения значений, полученных методом экстраполяции от ближайших точных значений.
6. Проанализирован полный информационный банк данных элементов орбит, известных на 10 мая 2005 года астероидов, при этом:
— выделены астероиды групп Аполлона, Амура, Атона, как потенциально опасные для столкновения с Землей;
— проведен статистический анализ распределения количества астероидов групп Аполлона, Амура и Атона как функций эксцентриситета, наклонения и большой полуоси;
— проведена классификация астероидов групп Аполлона, Амура и Атона по величине критерия Тиссерана, характеризующего общность происхождения малых тел;
— выявлено 1230 потенциально опасных на предмет столкновения с Землей астероидов.
7. С помощью специально разработанной программы численного интегрирования уравнений движения малых тел Солнечной системы, на основе высокоточного метода Тейлора проведено исследование эволюции орбит потенциально опасных астероидов на 1000 летнем интервале времени с 2005 по 3005 года. Получена новая информационная среда в виде банка данных эволюций орбит потенциально опасных для Земли астероидов.
8. На основании анализа полученного расчетным путем банка данных эволюций орбит астероидов:
— исследовано сохранение принадлежности астероидов к группам;
— проверена устойчивость движения астероидов при нахождении в резонансе с большими планетами.
9. Выявлены 49 астероидов, которые в процессе своей эволюции на интервале времени с 2005 по 3005 года пройдут сквозь сферу действия Земли. Составлен каталог эволюций орбит этих астероидов.
ЗАКЛЮЧЕНИЕ
.
