CC BY
78
2015 ВЕСТНИК САНКТ-ПЕТЕРБУРГСКОГО УНИВЕРСИТЕТА Сер. 1. Том 2(60). Вып. 1
АСТРОНОМИЯ
УДК 523.44, 52-17
ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ВОЗДЕЙСТВИЯ СВЕТОВОГО ДАВЛЕНИЯ НА ДВИЖЕНИЕ АСТЕРОИДОВ, В ТОМ ЧИСЛЕ СБЛИЖАЮЩИХСЯ С ЗЕМЛЕЙ
А. А. Мартюшева1, Н. А. Петров2, Е. Н. Поляхова2
1 Главная Астрономическая Обсерватория РАН,
Российская Федерация, 196140, Санкт-Петербург, Пулковское шоссе, 65
2 С.-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7-9
В работе исследуется влияние светового давления Солнца на эволюцию движения нескольких групп относительно небольших (сотни метров), различных по своим физическим и орбитальным параметрам, модельных астероидов. Исследования проводились для трех областей Главного пояса астероидов — в окрестностях резонансных зон 2:1, 3:2, 1:1 с Юпитером, а также в резонансной зоне 1:1 с Землей, на различных интервалах времени. Начальные условия для модельной системы трех тел были взяты близкими к реальным соотношениям в Солнечной системе. Численным интегрированием уравнений движения получены количественные оценки, характеризующие отклонения в движении астероида, вызванные световым давлением. Максимальные изменения в орбитальном движении достигаются вблизи так называемой «зоны скачка», связанной с резонансом. Порядок изменений величины полного смещения астероида в случае резонансов с Юпитером составил 103 км за 100 лет; в случае резонанса с Землей — (1-4)-103 км. Для четырех резонансных областей при различных начальных значениях большой полуоси и эксцентриситета, вычислено время, за которое величина полного смещения астероида достигает заданного отклонения (радиуса Земли), порядок величины составляет 10—10000 лет. Библиогр. 7 назв. Ил. 8. Табл. 3.
Ключевые слова: астероиды, световое давление, орбитальные возмущения, резонансы.
Введение. Вопрос о влиянии давления солнечной радиации на движение малоразмерных астероидов, сближающихся с Землей (АСЗ), широко обсуждается последние десятилетия в связи с осознанием мировой и научной общественностью проблемы астероидно-кометной опасности. Имеются в виду АСЗ диаметром порядка сотен метров. Хотя радиационные эффекты в движении астероидов, как известно, малы, тем не менее на относительно больших интервалах времени они могут приводить к заметным изменениям параметров орбитального движения и могут быть одним из механизмов, участвующих в формировании группы астероидов, сближающихся с Землей, и дальнейшей эволюции их орбит. В данной работе исследовано влияние светового давления Солнца на движение модельных астероидов в окрестностях резонансных зон 2:1, 3:2, 1:1 с Юпитером, а также в резонансной зоне 1:1 с Землей на
нескольких временных интервалах с целью оценки максимальных значений влияния этого эффекта на орбиты, находящиеся вблизи областей неустойчивого движения.
Модели движения и начальная конфигурация тел. Для решения этой задачи выбрана система четырех тел (Солнце —Юпитер —астероид—астероид), которая распадается на две независимые модели:
• ограниченная плоская круговая задача трех тел (Солнце — Юпитер — астероид (точечный)), где астероид рассматривается как материальная точка, на которую действует только гравитация;
• фотогравитационная ограниченная плоская круговая задача трех тел (Солнце — Юпитер — астероид (протяженный)), где астероид рассматривается как сферическое тело, имеющее конечный радиус, заданную плотность и определенный оптический коэффициент, на который действует как гравитация, так и световое давление
Таким образом, сравнение двух моделей позволяет вычислить отклонения положений одного астероида относительно другого, вызываемые действием сил светового давления. Кроме того, система пяти тел, включающая Землю в качестве пятого тела (Солнце — Юпитер — астероид — астероид — Земля), распадается на четыре независимые модели трех тел разной конфигурации и позволяет выявить влияние Земли на движение астероида, что актуально в связи с проблемой астероидно-кометной опасности. Движение тел задается в гелиоцентрической прямоугольной системе координат. Движение Юпитера принимается круговым. Начальные орбиты астероидов в общем случае эллиптические. Начальные значения долготы перигелия и средней аномалии для обоих астероидов равны нулю.
Оптический коэффициент. К фотогравитационным задачам небесной механики принято относить задачи о движении малых небесных тел естественного или искусственного происхождения, испытывающих в силу своих физических свойств существенное влияние давления солнечных лучей. Для учета светового давления необходим переход от традиционной модели гравитационного поля в космическом пространстве к модели фотогравитационного поля, в котором наряду с гравитационным притяжением существует и радиационное отталкивание. Изменение орбитальной энергии тела приводит к эволюции его орбиты, темп и интенсивность которой определяется уровнем воздействия сил светового давления на это тело. Как известно, сила светового давления определяется формулой [1]
Е(т)А £©А
— к е — к — ~ е,
с 4 пг2 с
где Е(г) = у — освещенность на расстоянии г от Солнца, — светимость Солнца, А — площадь миделева сечения, с — скорость света, к — оптический коэффициент астероида, е — единичный вектор радиального гелиоцентрического направления. Поскольку сила тяготения
ге — ^ г2 е'
где О — гравитационная постоянная, т© — масса Солнца, т — масса притягиваемого тела, сумма сил может быть представлена в виде
П = ^ + ^ = к^-ё= -с^а - /3)е = -С^е,
где ¡3 = к-ф ^^ — фотогравитационная редукция массы Солнца, <3© = — вспомогательная константа, Шд = ш©(1 — в) —редуцированная масса Солнца, к — оптический коэффициент астероида. В общем случае оптическая модель поверхности астероида исходит из предположения о преобладании диффузного характера переизлучения света, причем условие непрерывности потока световой энергии есть а + р + 6 = 1, где а — коэффициент поглощения, р + 6 — коэффициент отражения, складывающийся из коэффициентов зеркального отражения р и диффузного отражения 6. Для диффузного переизлучения оптический коэффициент имеет вид [2] к = а + р + Для природного тела принимаем, что зеркальное отражение отсутствует (р = 0), а коэффициент диффузного отражения 6 есть геометрическое альбедо астероида. Итак, для неточечного сферического диффузно отражающего протяженного тела а = 1 — 3, к = а-\- = 1 + В настоящих расчетах альбедо поверхности астероида принимали в пределах от нуля до единицы 6 € [0; 1]. Таким образом, для любого альбедо можно найти пропорциональный оптический коэффициент к € [1; 1.44]. Если протяженный астероид считать черным телом, то его оптический коэффициент к = а = 1, то есть в первом приближении можно считать, что все излучение поглощается и не приводит к нагреву астероида, а вся поглощенная энергия переходит в орбитальное движение. Для точечного астероида оптический коэффициент к = 0, возмущения нет.
Уравнения движения и начальные условия. Для записи двух систем уравнений ограниченной задачи трех тел вводим индексы: 1 —Солнце, 2 — Юпитер, 3 — астероид. Уравнения движения в гелиоцентрической прямоугольной системе координат для двух систем, состоящих из трех тел, имеют следующий вид:
'¿з = ~Стф + + Ст.р^,
,¿3 = 0.
Здесь фотогравитационная редукция массы Солнца /3 = к®р°, к — оптический коэффициент поверхности астероида, Ро — световое давление на орбите Земли на 1 м2, о-дг — гравитационное ускорение на расстоянии 1 а.е. от Солнца (м/с2), адг = Щ^2- — ускорение (м/с2), го — среднее расстояние от Земли до Солнца (м), в = ^ = -щ; — баллистический коэффициент (парусность) (м/кг2), А — площадь миделева сечения протяженного астероида (м2), ш — масса астероида (кг), Д — радиус астероида (м), р — плотность астероида (кг/м3). В первом случае принимаем, что на точечный астероид действует только гравитация (к = 0), на протяженный — гравитация и световое давление (6 = 0, к = 1). Во втором случае рассматриваем два протяженных астероида, на которые действует световое давление. При этом альбедо первого астероида считаем нулевым 6 = 0, к = 1 (черное тело), альбедо второго равно единице 6 = 1 , к = 1.44 (белое тело). В результате исследований по двум этапам предполагалось найти максимальные отклонения положений астероида с учетом светового давления относительно положений астероида без учета данного эффекта при различных начальных условиях. Были приняты следующие начальные условия.
• Солнце: Ш1 = 1.0; в гелиоцентрической системе координат координаты и скорости Солнца равны нулю.
• Юпитер: т2 = 0.000955т1, а2 = 5.202603 а.е. — большая полуось [3], е2 = 0.0 — эксцентриситет, ^2 = 0.0° —долгота перигелия, М = 0.0° —средняя аномалия. Поскольку орбита Юпитера круговая, имеем Х2 = 5.202603 а.е., у2 = 0.0 а.е., х2 = 0.0 а.е./зв. год, у2 = 2.755984 а.е./зв. год.
• Астероиды: Д3 = 0, тз = 0, причем тз ^ т2 < т-1. Параметры орбит астероидов (большая полуось и эксцентриситет) варьировались в процессе исследования. В начальный момент параметры орбит астероидов принимались одинаковыми.
При выборе физических параметров астероида нужно учитывать, что влияние светового давления ослабляется с уменьшением баллистического коэффициента. В то же время потенциальную опасность представляют только астероиды средних и больших радиусов. Поэтому для учета влияния светового давления были назначены следующие физические параметры протяженного астероида: Д =100 м — радиус, т = 4.2 • 109 кг = 2.1 • 10-21т1 — масса, р = 1000 кг/м3 —плотность, в = 7.5 • 10-6 м2/кг — баллистический коэффициент (парусность), адг = 5.9 • 10-3 м/с2 - гравитационное ускорение на расстоянии 1 а.е. от Солнца. Рассмотрены два случая. В первом случае рассматривались точечный астероид (к = 0) и протяженный астероид (к = 1), где к — оптический коэффициент поверхности астероида; в = 0 и в = 5.7 • 10-9 —фотогравитационные редукции массы Солнца соответственно. Во втором случае рассмотрены два протяженных астероида: к = 1, альбедо 6 = 0, в = 5.7 • 10-9; к = 1.44, альбедо 6 = 1, в = 8.5 • 10-9. В нашей задаче уравнения движения системы тел численно интегрировались методом Эверхарта 19 порядка [4]. Для вычислений приняты параметры модельного астероида. Как следует из правых частей уравнений движения (*), малые возмущения от светового давления приводят к малым отклонениям одного астероида относительно другого. Эти малые возмущения прямо пропорциональны оптическому коэффициенту к и обратно пропорциональны радиусу астероида Д и его плотности р. Поэтому, во сколько раз величины к, Д и р для истинного астероида отличаются от тех же величин модельного астероида, во столько раз и малые отклонения истинной пары астероидов отличаются от малых отклонений модельной пары астероидов.
Исследуемые области движения астероидов. Проводились исследования движения астероидов, орбиты которых находились в трех областях Главного пояса астероидов —в окрестности резонансных зон 2:1, 3:2, 1:1 с Юпитером (рис. 1) и в резонансной зоне 1:1 с Землей. Орбиты астероидов имели большую полуось и эксцентриситеты соответственно из резонансной зоны 2:1 с Юпитером: аз„ = 3.277435 а.е., е3 € [0.0; 0.9]; из резонансной зоны 3:2 с Юпитером: а3п = 3.970329 а.е., е3 € [0.0; 0.8]; из резонансной зоны 1:1 с Юпитером: а3п = 5.202603 а.е., е3 € [0.2; 0.9]; из резонансной зоны 1:1 с Землей: а3„ = 1.000000 а.е., е3 € [0.0; 0.9]. Вычислялись следующие три величины: Дг = г3 — г 4 — смещение астероида вдоль гелиоцентрического радиус-вектора (км), А1 = /3 — /4 —смещение астероида по долготе (км), Ас1 = VДг2 + АР— полное смещение астероида (км).
Для того чтобы понять общий характер изменения смещений Дг, Д/, Дй по времени, были проведены вычисления для начальных данных, соответствующих узлам решетки (Да3,е3), в каждой из трех окрестностей резонансных зон Главного пояса астероидов 2 : 1, 3 : 2, 1 : 1 с Юпитером, а также в резонансной зоне 1 : 1 с Землей с достаточно большими шагами по большой полуоси и эксцентриситету. На рис. 2-4 представлены области устойчивости движения астероида по Лагранжу с на-
2.00 2.50 3.00 3.50 4.00 4.50 5.00
2/1 3/2 1/1
а, а.е
Рис. 1. Области устойчивости и неустойчивости (заштрихованные области) движения астероида. Пунктирными вертикальными линиями показаны точные невозмущенные резонансы для резонансных областей с Юпитером.
несенными решетками узлов (Даз, ез), где Даз = аз — а3„ была введена для удобства представления многоразрядных величин, первые цифры которых остаются постоянными. К примеру, для резонанса 2 : 1 значение а = 3.077435 а.е. можно записать как Да = —0.2. Представление информации в таком виде значительно упрощает восприятие большого объема данных. Оси абсцисс на рис. 2-4 представлены в относительных координатах. Каждый рисунок можно условно разделить на три области: область устойчивости (линии), неустойчивости движения (пустоты, отдельные точки) и переходная область. Область устойчивости движения характерна тем, что за время интегрирования траектория астероида находится в кольце, границы которого определяются с большой точностью и остаются неизменными со временем. Область неустойчивости движения характеризуется тем, что в процессе интегрирования проявляется хаотическое изменение минимальных и максимальных значений параметров орбиты вследствие тесных сближений с возмущающим телом. Результатом таких изменений может быть уход из Солнечной системы, падение на Юпитер или падение на Солнце. Можно предположить, что практически все точки из зоны неустойчивости имеют подобный финал, но некоторые за более длительный период, чем время интегрирования в данной работе. В переходной зоне устойчивые и неустойчивые траектории перемешиваются настолько сильно, что разделить их не удается, поэтому трудно выявить четкие границы устойчивости, так как, по-видимому, проявляются неустойчивости, вызванные влияниями других резонансов.
Результаты вычислений и их анализ. Первоначально вычисления проводились на двух интервалах: 100 и 1000 лет. Исходя из предположения о линейности изме-
а.е
Рис. 2. Область устойчивости движения астероида для резонанса 2:1 с Юпитером (заштрихованная область). Вычисления выполнены в 36 узлах решетки (Ааз, ез), а3п = 3.277435 a.e.
Да,, а.е
Рис. 3. Область устойчивости движения астероида для резонанса 3:2 с Юпитером (заштрихованная область). Вычисления выполнены в 33 узлах решетки (Ааз, ез), азп = 3.970329 a.e.
нения величин по времени (вследствие малости возмущения от светового давления), при увеличении времени интегрирования в 10 раз ожидали получить десятикратное увеличение смещений Дг, Д/, Дй. Результаты показали нелинейность характера изменения Дг, Д/, Дй по времени (табл.1): их изменение для некоторых начальных значений эксцентриситета оказалось существенно непропорционально времени. Максимальные значения по Д/ достигаются тогда, когда изменения по Дг минимальны (см. развертку взаимных изменений величин Д/ и Дг (рис.5)), поэтому максималь-
0001141921530248892927010200
Да3, а.е
Рис. 4. Область устойчивости движения астероида для резонанса 1:1 с Юпитером (светлая область). Вычисления выполнены в 38 узлах решетки (Ааз, ез), азп = 5.202603 a.e.
ные значения по Д/, как правило, близки к значениям Дй. Для того чтобы выяснить причину полученных нелинейных отклонений были предприняты более подробные исследования, а также сделаны графические представления отклонений в виде разверток по времени и взаимных изменений величин Дг и Д/.
Дг, км
Д/, км
Т, лет
ДI, км
Рис. 5. Смещения Аг, А1, А<1 (км) астероида на интервале времени 1000 лет. Резонансная зона 2:1 с Юпитером, а = 3.377435 а.е. (Аа = 0.1 а.е.), е = 0.2.
Таблица 1. Значения величин Дг, Д1, Д<1 в узлах решетки в окрестности резонанса 2:1 с ^Юпитером на интервалах времени 100 и 1000 лет для некоторых начальных значений эксцентриситета
е а (а.е.) Да (а.е.) Дг (км) А1 (км) Ас1 (км)
100 лет 1000 лет 100 лет 1000 лет 100 лет 1000 лет
0.2 3.077435 -0.200 148 1832 979 10965 979 10965
3.177435 -0.100 181 3214 1226 18197 1227 18198
3.277435 0.000 123 223 713 1112 713 1112
3.377435 0.100 169 889 648 3926 648 3926
3.477435 0.200 107 624 648 2804 649 2804
0.4 3.077435 -0.200 247 1565 865 5416 866 5418
3.177435 -0.100 212 2356 761 8570 761 8570
3.277435 0.000 355 395 1225 1346 1227 1346
3.377435 0.100 292 4047 919 12649 919 12653
0.6 3.177435 -0.100 580 7835 1551 21742 1555 21771
3.277435 0.000 930 1045 2461 2744 2462 2748
3.377435 0.100 1163 9224 2908 23479 2924 23479
0.8 3.177435 -0.100 2144 24194 4619 54285 4718 54580
3.277435 0.000 2886 7885 6427 14309 6441 17595
Анализ графического материала показал, что на малом временном промежутке хорошо видны короткопериодические колебания всех трех величин, обусловленные движением астероида вокруг Солнца, и долгопериодические колебания, обусловленные либрационным движением астероида. Полученные оценки средних значений периодов короткопериодических и долгопериодических колебаний согласуются с результатами по исследованию ограниченной плоской круговой задачи трех тел [5]. Детальное исследование изменений величин отклонений позволило сделать вывод, что причиной нелинейности является либрационное движение, которое в свою очередь существенно зависит от начальных значений большой полуоси и эксцентриситета.
В ходе вычислений для резонансной зоны 2:1 с Юпитером было установлено, что величины смещений Дг, Д1, Д! интенсивно возрастают (почти экспоненциально) по мере приближения с двух сторон к узкой области, проходящей через точки (Да = 0.05883082; е = 0.0), (Да = 0.08714965; е = 0.1), (Да = 0.0152; е = 0.2), см. рис. 2. При удалении от так называемой «зоны скачка» в обе стороны значения величин заметно уменьшаются, что можно проследить и по величине либрационно-го периода Р, который вычислялся для контроля. В связи с этим, было проведено специальное исследование для этой особой области, на интервале времени 5000 лет. По мере приближения к «зоне скачка» с двух сторон, исследуемые области (слева и справа) разбивались на 10 частей. Из них выбирались две соседние наиболее близкие к «зоне скачка» точки, интервал между которыми также разбивался на 10 частей, и так далее. На каждом этапе разбиений вычислялись величины Дг, Д1, Д! Таким образом, для значений эксцентриситета 0.0 и 0.1 таких этапов разбиений получилось всего семь, то есть расстояние между двумя соседними точками слева и справа от «зоны скачка» удалось уменьшить до 10-8 а.е. « 1.5 км, а значит, удалось очень близко подойти с двух сторон к «зоне скачка». При дальнейшей попытке разбиения окрестности (уменьшения интервала Д£) вычисления привели к выбросам, что заметно по данным табл. 2, поэтому, дальнейшее разбиение прекращалось. По этой же причине для е = 0.2 таких этапов вычислений получилось всего три. Вычисления при различных значениях эксцентриситета привели к заключению: ширина «зоны скачка» увеличивается с ростом эксцентриситета.
Таблица 2. Полное смещение астероида А<1 за время 5000 лет и либрационный период Р. Астероид находится в окрестности «зоны скачка» для резонанса 2:1 с ^Юпитером
AS (а.е.) Ad (км) Р (лет)
ю-'2 6.729 • 104 1.431 • 10Б 556 385
5.260 • 10ь 6.610 • 10Б 833 500
10"4 9.259 • 10ь 1.427- 106 1250 556
10~ь 2.211 • 10» 8.335 • 106 1667 714
кгь 2.211 • 10й 6.336 • 108 1667 1000
10-' 1.000 • 10у 1.059 • 109 1000 1000
10-» 1.026 • 10у 1.034 • 109 1000 1000
е = 0.0, а = 3.33626582 (Аа = 0.05883082) а.е., А5 — расстояние между двумя соседними точками с двух сторон от «зоны скачка»
После полученных численным интегрированием результатов предпринята попытка их обобщения. Для этого выполнено исследование некоторых функций в плоскости переменных (а, e). Вначале исследовалось поведение максимального отклонения как функции начальных переменных /i: Ad = /i(p4, Д4, T, а4, e4), затем при фиксированных Д4, ^4, р4 как функция /2: Ad = /2(T, а4, e4) и для фиксированных а4, e4 как функция /3: Ad = /3 (T). Были введены функции максимального смещения астероида dmax (км) и минимального времени Tmin (лет) достижения заданного отклонения следующим образом:
dmax(T) = max /з(4), 0 < t < T, Tmin = {min T : dmax(T )= Де } ;
dmax(T) < Де, T < Tmin-
Экстремумы этих функций представляют собой двумерные поверхности величин максимального смещения астероида dmax (км) и минимального времени Tmin (лет) достижения некоторого заданного отклонения, обусловленного световым давлением. В качестве такого отклонения было принято значение радиуса Земли Де . Исследование проводилось методом линейных сечений этих поверхностей вдоль большой полуоси с мелким шагом da = 0.001 (а.е.) при фиксированных значениях эксцентриситета для каждой из трех резонансных областей с Юпитером и для резонансной зоны 1 : 1 с Землей. Использовались результаты численного интегрирования, полученные на интервалах времени 5000 и 500 лет. Рисунок 6 иллюстрирует схему вычисления минимального времени Tmin (лет) достижения критического расстояния Де из функции максимального отклонения, зависящей только от времени. Линейные сечения позволяют выяснить, как выглядит поверхность dmax и Tmin и какие у неё особенности.
Для каждого значения эксцентриситета были получены сложные картины линейных сечений с проявлениями «всплесков». Эти «всплески» сильно искажают поведение кривых линейного сечения в области устойчивости, за счет большого количества пиков, которые увеличиваются и учащаются по мере приближения к границам
Дс/, км
8000
Рис. 6. Схема вычисления минимального времени Тт;п достижения Не по функции максимального отклонения ётах, зависящей только от времени. Случай резонанса 2 : 1 с Юпитером, а = 3.377435 a.e. (Да = 0.1 а.е.), е = 0.2.
неустойчивости движения. Поэтому для выявления экстремумов гладкого поведения этих кривых подобные «всплески» были сглажены. По-видимому, наблюдаемое проявление «всплесков» обусловлено проявлениями резонансов более высоких порядков. Для примера в крупном масштабе на рис. 7 показаны линейные сечения вдоль большой полуоси а двумерной поверхности !тах и Тт¡п.
Да, а.е Да, а.е
Рис. 7. Вычисление линейных сечений двумерной поверхности (а, е) вдоль большой полуоси с шагом ёа = 0.001 (а.е.) для случая, когда эксцентриситет е = 0.0. Резонансная зона 2:1 с Юпитером.
Для каждого линейного сечения были получены три значения, характеризующие это сечение — максимум слева, максимум справа (относительно точного невозмущенного резонанса а3„ для каждой области) и минимум. Затем были получены значения максимумов и минимумов по всем линейным сечениям, то есть была получена зависимость !тах и Тт\п от эксцентриситета. И, как следствие, для каждой из
резонансных областей найдены глобальные экстремумы по всей поверхности для величины максимального смещения астероида dmax (км) и минимального времени Tm¡n (лет) достижения заданного отклонения Re (табл.3). Для трех резонансных областей с Юпитером глобальные экстремумы были достигнуты при различных значениях эксцентриситета. Для резонансной области 1:1 с Землей величины минимумов и максимумов dmax (км) и Tmin (лет) уменьшаются и увеличиваются в зависимости от эксцентриситета гораздо равномернее, что согласуется с наиболее линейной картиной сечений по сравнению с другими резонансами. Между графиками зависимости экстремумов максимального смещения астероида dmax (км) и минимального времени Tmin (лет) достижения радиуса Земли RE от эксцентриситета для каждой из четырех резонансных областей заметно хорошее инверсивное соответствие.
Рис. 8. Графики зависимости экстремумов максимального отклонения астероида ^тах и минимального времени Тт;п достижения значения Не от эксцентриситета е на интервале времени 5000 лет для всей резонансной области 2:1 с Юпитером. Шаг по большой полуоси <1а = 0.001 а.е., е € [0.0; 0.8], Да € [-0.2; 0.2] а.е.
Таблица 3. Глобальные экстремумы по всем резонансным областям
Резонанс 2:1 Резонанс 3 : 2 Резонанс 1 : 1 Резонанс 1: 1
с Юпитером с Юпитером с Юпитером с Землей
924 км 842 км 641 км 4953 км
min <imax (а = 3.227435, (а = 3.958329, (а = 5.255603, (а = 1.199,
е = 0.2) е = 0.4) е = 0.5) е = 0.0)
2924364 км 98586 км 429114 км 128065 км
Шах C?max (а = 3.164435, (а = 3.900329, (а = 4.909603, (а = 0.8,
е = 0.3) е = 0.6) е = 0.8) е = 0.8)
78 лет 388 лет 94 года 25 лет
min Tmin (а = 3.356435, (а = 3.909329, (а = 4.909603, (а = 0.807,
е = 0.8) е = 0.1) е = 0.8) е = 0.8)
34475 лет 37851 год 49711 лет 643 года
шах 7,, i, (а = 3.226435, (а = 3.958329, (а = 5.255603, (а = 1.2,
е = 0.2) е = 0.4) е = 0.5) е = 0.0)
Оценка светового давления для астероида Апофис. Были приняты следующие начальные физические условия для астероида Апофис (см. http://ssd.jpl.nasa.gov): Я = 130 м, р = 3000 кг/м3, тогда т = 2.75 • 1010 кг = 1.39 • 10-20т1; в = 1.93 • 10-6 м2/кг, р = 0.33, отсюда к = 1.15, и, следовательно в = 1.70 • 10-9. Для орбитальных параметров было принято: е = 0.191, а = 0.922
а.е. Численное интегрирование (t = 0 и t = 20 лет — начальный и конечный моменты интегрирования соответственно) показало, что на этом интервале максимальное отклонение астероида Апофис, обусловленное световым давлением, составляет вдоль гелиоцентрического расстояния |Дг| = 16 км, а вдоль орбиты— |Д1| = 106 км и полное смещение Д< = 107 км, что соответствует в среднем изменению за год |Дг| « 0.8 км, |Д1| « 5.3 км, Д< « 5.4 км соответственно, что согласуется с опубликованными данными (см. [6], возмущения, которые необходимо учитывать при уточнении орбиты и прогнозе движения: гл. 7, пункт 7.3, с. 195).
Заключение. По результатам вычислений изменений радиус-векторов модельных астероидов получены количественные оценки влияния светового давления на орбитальную эволюцию в зависимости от начальных значений большой полуоси и эксцентриситета. Для исследования были выбраны следующие параметры сферического астероида: радиус R = 100 м, плотность р = 1000 кг/м3, оптический коэффициент k = 1. Для резонансных зон 2:1, 3:2, 1:1 с Юпитером и резонанса 1 : 1 с Землей наибольшие значения величины максимального смещения астероида dmax были достигнуты на границах областей устойчивости и неустойчивости движений. Отклонения растут с увеличением начального эксцентриситета. Максимальные изменения величины полного смещения астероида Д< были достигнуты вблизи зоны скачка, связанной с резонансом 2 : 1. Для каждой из областей вычислено минимальное время, за которое отклонение астероида достигает значения, равного радиусу Земли. Выполненное исследование показывает, что световое давление может быть важным фактором, который необходимо учитывать при изучении движения астероидов на значительных интервалах времени, особенно тех астероидов, орбиты которых находятся вблизи резонансных зон с Юпитером и Землей.
Литература
1. Радзиевский В. В. Фотогравитационная небесная механика. Н.Новгород.: Изд. Николаев Ю. А., 2003. 196 с.
2. Поляхова Е. Н., Шмыров А. С. Физическая модель сил давления световой радиации на плоскость и сферу // Вестн. С.-Петерб. ун-та. Сер. 1. 1994. Вып. 2 (№8). C.87-104.
3. Bretagnon P., Francou G. Planetary solutions VSOP87. 1988.
4. Бордовицына Т. В. Современные численные методы в задачах небесной механики. М: Наука, 1984. C. 67-70.
5. Петров Н.А., Орлов С. А. Исследование с помощью программы TRIO свойств движения в резонансной зоне 2:1 в ограниченной задаче трех тел // Всероссийская астрономическая конференция 6-12 августа 2001 г. СПб. Тез. докл. СПб: НИИХ СПбГУ, 2001. C. 142.
6. Астероидно-кометная опасность: вчера, сегодня, завтра / под ред. Б. М. Шустова, Л. В. Рыхловой. М.: ФИЗМАТЛИТ, 2010. 385 с.
Статья поступила в редакцию 23 октября 2014 г.
Сведения об авторах
А.лекса,ндра, Андреевна Мартюшева — младший научный сотрудник; alex.mart13@gmail.com
Никита Александрович Петров — научный сотрудник; petrov@astro.spbu.ru
Елена Николавна Поляхова —доцент; pol@astro.spbu.ru
SOLAR RADIATION PRESSURE EFFECTS ON ASTEROID MOTIONS, INCLUDING NEAR-EARTH OBJECTS
Alexandra A. Martyusheva1, Nikita A. Petrov2, Elena N. Polyakhova2
1 Pulkovo Observatory RAS, Pulkovskoe shaussee, 65, St. Petersburg, 196140, Russian Federation; alex.mart13@gmail.com
2 St.Petersburg State University, Universitetskaya nab., 7-9, St.Petersburg, 199034, Russian Federation; petrov@astro.spbu.ru, pol@astro.spbu.ru
The orbital evolution of several groups of relatively small asteroids (of about several hundreds meters in size) is investigated by numerical integration, with solar radiation pressure effects taken into account. Basic equations are written in the frame of the three-body problem (Sun-planet-asteroid) together with the Solar System planetary perturbations. Numerical integration was carried out to model objects with diverse physical properties and orbital parameters. Three resonance regions with Jupiter (2:1, 3:2, 1:1) and one region with the Earth (1 : 1) are considered. Initial conditions for integration were chosen to correspond to real relationships in asteroidal belts. It was discovered that solar radiation pressure effects increase significantly with orbital eccentricity for a given asteroid orbit, no matter how strong all other perturbing actions are. Maximal variations of orbital elements occur near the above-mentioned resonance regions with planets. For example, the derived estimates show that the total variation of the distance between an asteroid's non-perturbed orbit and the orbit as perturbed by solar radiation pressure force is as high as 103 km over 100 years for Jupiter and (1—4)x103 km for the Earth (where solar radiation is more intensive). In terms of the actual asteroid hazard for the Earth, the estimated time until a hazardous close approach (of about the Earth one radius) is 10—10000 years when solar radiation pressure is taken into account. Refs 7. Figs 8. Tables 3.
Keywords: asteroids, solar radiation pressure, orbital perturbation, resonances.
References
1. Radzievsky V. V., The photogravitational celestial mechanic (publisher Nikolaev Yu.A., N. Novgorod, 2003) [in Russian].
2. Polyakhova E. N., Shmyrov A. S., "The physical model of the solar radiation pressure on a sphere and a plane", Vestnik of St.Petersburg University, Series 1 Issue 2(8), 87—104 (1994) [in Russian].
3. Bretagnon P., Francou G., Planetary solutions VSOP87 (1988).
4. Bordovitsyna T. V., The modern numerical methods in the celestial mechanic problems 67—70 (Nauka, Moscow, 1984) [in Russian].
5. Petrov N. A., Orlov S. A., "The investigation of the motion properties in the 2:1 resonance zone in the restricted three-body problem by means of the TRIO program", Proceedings of the All-Russian Astronomy Conference, St. Petersburg, August 6-12, 2010 (Institute of Chemistry SPbU, St.Petersburg, p. 142) [in Russian].
6. The asteroid-comet hazard: yesterday, today, tomorrow (Shustov B.M., Rykhlova L. V. (ed.); FIZMATLIT, Moscow, 2010).
