Эволюция планетарных систем - (реферат)
Дата добавления: март 2006г.
Сабелев Г. И.
ЭВОЛЮЦИЯ ПЛАНЕТНЫХ СИСТЕМ
Москва, 1997
ББК. 527
Ш66
УДК 523. 84
АННОТАЦИЯ
Рассмотрена эволюция планетных систем на основе теории гравитационного захвата нейтральных тел. Сабелев Георгий Иванович. Эволюция планетных систем. -М. , ИВЦ Маркетинг, 1997. -стр. 79, ил. 11 г Информационно-внедренческий центр “Маркетинг”, Москва, 1997. Содержание
стр.
1. Введение ..........................................
4
2. Гравитационный захват нейтральных тел ......
6
3. Орбитальное движение тел......................
13
4. Уход тела с орбиты ...............................
31
5. Эволюция планет и комет ........................
43
6. Эволюция спутников планет .....................
65
7. Эволюция Земли ..................................
74
Литература .......................................
79
Любимой жене, другу и помощнику
Сабелевой Валентине Ивановне
с благодарностью посвящаю
1. Введение
Два тела, взаимодействующих в центральном поле, составляют систему тел. Однако с одним центральным телом может одновременно взаимодействовать любое количество тел, каждое из которых составляет с центральным телом свою cистему двух тел. В свою очередь центральное тело может быть периферийным в другой системе двух тел. Именно в этом плане рассматривается эволюция тел и систем, движущихся в центральных полях. Любая система материальных тел проходит во времени через этапы возникновения, развития и распада; в этом смысле принято говорить об эволюции системы. Планетная система представляет собой также систему материальных тел, взаимодействующих друг с другом и совершающих на разных этапах относительное движение под влиянием этих взаимодействий. Определяющим при этом является центральное взаимодействие двух тел. Планетная система состоит из центрального тела (Солнца), планет со спутниками, комет и метеоров. Общим для периферийных тел планетной системы является отсутствие значительных по величине собственных электрических и магнитных зарядов. Это позволяет считать периферийные тела планетной системы нейтральными пo отношению к электромагнитным взаимодействиям и, следовательно, рассматривать только гравитационное взаимодействие между периферийным и центральным телом. Теория центрального взаимодействия двух тел изложена в [ 1 ]. Там же рассмотрены некоторые вопросы гравитационного взаимодействия нейтральных тел и отмечены основные этапы эволюции планет, спутников и комет. Однако остались за пределами книги многие важные проблемы эволюции, имеющие непосредственное отношение к эволюции Земли и, следовательно, к эволюции жизни. Поэтому представляет интерес более детально рассмотреть эти проблемы, уделив особое внимание анализу факторов, непосредственно влияющих на условия существования жизни в планетной системе. 2. Гравитационный захват нейтральных тел В случае центрального взаимодействия нейтральных (незаряженных) тел возможно лишь их встречное движение под действием гравитационного притяжения. При этом может произойти гравитационный захват нейтрального тела меньшей массы телом большей массы: на некотором расстоянии между теламиroначнется их встречное движение. В дальнейшем прямолинейное встречное движение тел перейдет в движение по эллиптической орбите тела меньшей массы вокруг центрального тела. Рассмотрим вначале движение двух гравитационно взаимодействующих тел на этапе их прямолинейного сближения. Поскольку этапы прямолинейного и орбитального движения тел сопряжены в пространстве и во времени, то в качестве характерного масштаба скорости выберем среднюю орбитальную скоростьV, а в качестве характерного линейного масштаба - большую полуось орбитального эллипсаa. Тогда относительная скорость движения нейтральных тел на прямолинейном участке , (2. 1)
где , , ,
U - относительная скорость тел,
r - текущее расстояние между телами,
ro- расстояние между телами в момент начала захвата,
M1 , M2 - массы взаимодействующих тел ( M2 < M1 ).
Первая производная от (2. 1) по r
. (2. 2)
Приравняв производную нулю, найдем координату максимума скорости (2. 3) и, далее, подставив (2. 3) в (2. 1), найдем максимальную скорость в системе двух тел
. (2. 4)
Ускорение движения тел
, (2. 5)
.
При r = rm s = 0.
Свойства функций (2. 1), (2. 2) и (2. 5) рассмотрены в [1].
Координаты экстремумов функции s = f(r, N) найдем, дифференцируя (2. 5) по r и приравняв производную нулю. Получим . (2. 6)
Величину ускорения при r = rэ найдем из (2. 5)
. (2. 7)
Координаты экстремумов функции sэ = f( N )
,
.
Этим значениям N соответствует r э = 0, 5.
Величину s э при N1 и N2 найдем из (2. 7)
s э1 = ± 0, 84715; s э2= ± 1, 84715.
При r = r э скорость движения тел
. (2. 8)
При N << 1 0, 9428086.
Энергия движения тел
, (2. 9)
где .
Энергия покоя (при r = r m)
. (2. 10)
Пороговая энергия (при r = r э)
. (2. 11)
Следовательно
. (2. 12)
При N << 1 .
Зависимости r min = r э = f( N ) однозначно определяют границу сферы, на которой происходит переход от прямолинейного движения к орбитальному, а зависимостиr m = f( N ) - границу сферы максимума относительной скорости тел [ 2 ]. Выясним, при каких условиях может происходить орбитальный переход. Очевидно, на границе перехода rmin должно выполняться условие Ј , (2. 13)
где Umin - относительная скорость тел при r = rmin.
Но Umin = Vb min и, следовательно, в безразмерном виде условие (2. 13) будет . (2. 14)
Подставив (2. 8) в (2. 14), получим неравенство
. (2. 15)
Значения N < 0, 0314, при которых выполняется это неравенство, найдем, приравняв числитель и знаменатель (2. 15). Следовательно, условие перехода от прямолинейного движения к орбитальному состоит в том, что на орбиту переходят лишь тела, относительная масса которых N < 0, 0314. Тела большей относительной массы не подвержены орбитальному переходу. Время встречного движения при гравитационном захвате нейтральных тел на интервале[1, r min ] определяется зависимостью
. (2. 16)
Здесь , G = 6, 672Ч10-8 - гравитационная постоянная.
При r ® 0
. (2. 17)
Среднее время движения на интервале [1, rmin ] при N ® 0 будет . (2. 18) Время движения от точки захвата r = 1 до точки перехода одного из тел на эллиптическую орбиту вокруг другого тела r min, в соответствии с (2. 16), будет (2. 19) Cредняя скорость прямолинейного движения на интервале [1, rmin ] (2. 20)
При N<<1
(2. 21)
Энергия, затрачиваемая на интервале движения [1, rmin] при N<<1 будет . (2. 22) Зная ro , rmin и tmin , можно рассчитать среднеарифметическую скорость на интервале [ro , rmin] (2. 23) и затем определить скорость в точке rmin (точка перехода от прямолинейного движения к орбитальному) (2. 24) Из соотношения (2. 19) следует, что между расстоянием ro и временем движения tmin существует взаимосвязь. При N << 1 rmin << 1, поэтому справедливы следующие приближения
и зависимость (2. 19) можно представить в виде
(2. 25)
или . (2. 26)
Отметим, что зависимость (2. 26) аналогична третьему закону Кеплера для орбитального движения и отличается от него лишь постоянным множителем. Зависимости (2. 3), (2. 4), (2. 6), (2. 8) при N << 1 практически линейны относительно ln N. При N ® 0 эти зависимости будут
ln bm = - 0, 25 ln N - 0, 693147, (2. 27)
ln bmin = - 0, 25 ln N - 0, 752038, (2. 28)
ln rm = 0, 5 ln N + 0, 693247, (2. 29)
ln rmin = 0, 5 ln N + 1, 386294. (2. 30)
Приведенные основные характеристики гравитационного захвата нейтральных тел безразмерны. Для того, чтобы представит эти зависимости в реальном пространственно - временном измерении, необходимо знать величинуro , то есть начальное расстояние между телами. Оценку этой величины можно получить, сопоставляя характеристики движения тел на этапах гравитационного захвата и орбитального движения.
3. Орбитальное движение тел.
На границе происходит переход от прямолинейного движения к орбитальному: тело меньшей массы (N < 0, 0314) начинает двигаться по эллиптической орбите вокруг тела большей массы; при этом периферийное тело приобретает вращение вокруг собственной оси. Орбитальное движение совершается по эллипсу, параметр которого остается постоянным. При этом большая полуось эллипса , период обращения Т и средняя скорость движения по орбите V связаны соотношениями (3. 1)
или в безразмерном виде
(3. 2)
где .
В соответствии со схемой орбитального перехода, рассматриваемой в [2], большая полуось начального эллипса
, (3. 3)
а эксцентриситет начального эллипса
. (3. 4)
Так как параметр эллипса остается постоянным, то можно написать (3. 5) где индексы “н” и “к” означают начальное и конечное (современное) значение параметра, а индекс “о” означает значение параметра при е = 0.
Но и можно написать
откуда, с учетом (3. 3),
. (3. 6)
При N ® 0
(3. 7)
или , (3. 8)
. (3. 9)
Сопоставив зависимости (3. 1) и (2. 26), получим также
(3. 10)
(3. 11)
Зная параметры начальной и конечной (современной) орбит и , можно рассчитать начальную среднюю скорость VH и начальный период обращения , используя всегда верные соотношения (3. 1) и (3. 2): , (3. 12)
, (3. 13)
а также рассчитать значения параметров и V при эксцентриситете орбиты е = 0 (3. 14)
(3. 15)
. (3. 16)
При N ® 0
. (3. 17)
Далее определим скорость в перигелии
(3. 18)
усредненные параметры
(3. 19)
а также разности начальных и конечных параметров за время движения тела по орбите
(3. 20)
и отношения
При большом числе оборотов орбитального тела
(3. 21)
Но эти же отношения могут быть получены другим способом. Дифференцируя (3. 1) поV, получим
и, используя (3. 1)
(3. 22)
Дифференцируя (3. 1) по Т, получим
(3. 23)
Но и, следовательно,
(3. 24)
Подставив (3. 24) в (3. 23), получим
(3. 25)
Расчеты разностей и их отношений с использованием зависимостей гравитационного захвата и расчеты частных производных по зависимостям (3. 22), (3. 24), (3. 25), проведенные для планет и их спутников показали, что действительно соотношения (3. 21) выполняются. Теперь появляется возможность получить оценку числа оборотов, совершаемых нейтральным космическим телом за время орбитального движения и времени жизни на орбите.
Преобразуя (3. 24) к виду
, (3. 26)
интегрируем в начальных и конечных пределах
(3. 27)
где n - число оборотов, совершенных телом за время движения на орбите, - средний период обращения тела, - отношение среднего периода обращения Земли к единице времени 1 сек. (величина, равная , но безразмерная). Получим
(3. 28)
Аналогично, преобразуя (3. 25) к виду
(3. 29)
и учитывая, что
интегрируем(3. 29) в начальных и конечных пределах
. (3. 30)
Получим
. (3. 31)
Очевидно, зависимости (3. 31) и (3. 28) должны приводить к одинаковому результату, что и подтверждается расчетами. Другой подход к оценке числа оборотов, совершаемых орбитальным телом за время, соответствующее изменению эксцентриситета от до е=0, состоит в следующем. Угловое смещение перигелия орбитального тела постоянно от оборота к обороту. В самом деле, согласно А. Эйнштейну, угловое смещение перигелия орбитального тела
(3. 32)
где с - скорость света в вакууме.
Величины G, , с, входящие в это соотношение, есть константы, а произведение - параметр эллиптической орбиты, который за все время орбитального движения также остается постоянным.
Следовательно
(3. 33)
где
(3. 34)
где ,
, (3. 35)
где .
Дифференцируя эти соотношения по e, получим
, (3. 36)
, (3. 37)
. (3. 38)
Число оборотов, совершенных, например, от начала орбитального движения (эксцентриситет орбитыeн) до момента, когда эксцентриситет орбиты становится e=0 , (3. 39)
или (3. 40)
где - изменение большой полуоси орбитального эллипса, периода обращения и средней скорости движения по орбите за1 оборот. При большом числе оборотов ,
и можно написать
(3. 41)
или . (3. 42)
Из соотношений (3. 36), (3. 37), (3. 38) следует, что
, , . (3. 43)
Предположим, что эти соотношения справедливы и для средних значений, то есть
Тогда
или (3. 45)
Расчеты числа оборотов орбитальных тел (планет и спутников) по этим зависимостям и по зависимостям, аналогичным (3. 28) и (3. 31) приводят к одинаковым результатам. В соответствии с принятыми определениями
(3. 46)
Следовательно
(3. 47)
Но поэтому
(3. 48)
При N ®0 и
об.
Полученный результат означает, что любое орбитальное тело (планета, спутник, комета) от начала орбитального движения до моментае = 0 совершает примерно одинаковое число оборотов (разница числа оборотов обусловлена разницей значенийeн). Орбитальная скорость космического тела изменяется во времени. Изменение скорости происходит из-за превращения потенциальной энергии взаимодействия тел в кинетическую энергию орбитального движения. Очевидно, что скорость орбитального движения не может превысить некоторого предела, при котором произойдет переход от орбитального движения к движению по параболической траектории и освобождение орбитального тела от действия притяжения центрального тела.
Этим предельным условием будет
, (3. 49)
то есть скорость в перигелии станет равной максимальной скорости в системе двух тел.
Подставив (3. 18) в (3. 49), получим
(3. 50)
и так как , то
(3. 51)
При значении e = eкр произойдет переход от орбитального движения по эллипсу к движению по параболической траектории и освобождение орбитального тела от действия притяжения центрального тела. Зная величину eкр (критический эксцентриситет орбиты), а также величины , можно рассчитать величину , а затем последовательно значения остальных параметров, определяющих условия этого перехода.
Так . (3. 52)
Используя (3. 3), (3. 4), (3. 7), и (3. 51), получим при N << 1 , (3. 53)
. (3. 54)
Далее
, (3. 55)
. (3. 56)
При N << 1
, (3. 57)
. (3. 58)
Расстояние между телами в критическом перигелии , а в критическом афелии . Из полученных результатов следует, что , и . В то же время , и , то есть зависимости , , и должны иметь экстремум на интервале времени от начала орбитального движения до критического перехода (ухода с орбиты). Дважды дифференцируя поe (3. 33), (3. 34) и (3. 35), действительно убеждаемся, что функции и при e= 0 имеют минимум, а функция - максимум. Отметим, что несмотря на уменьшение средней орбитальной скорости на интервале , скорость в перигелии увеличивается, достигая критического значения при . На рис. 1 показан примерный вид этих зависимостей для Земли (аналогичный вид они имеют для других планет и спутников). Эти зависимости представляют собой параболы с вершиной прие = 0. Левая ветвь параболы начинается от линии , правая ветвь заканчивается при пересечении ею линии . Параболы частично симметричны относительно линии е = 0 (на интервале слева от (от до ) и на интервале справа от (от до ). Таким образом, в точках пересечения правой ветви параболы с прямой значения параметрических функций будут в точности равны их начальным значениям. Учитывая это свойство параметрических функций , и и интегрируя (3. 27) в пределах от до и от 0 до , получим зависимость (3. 59) или, интегрируя (3. 30) в пределах от до и от 0 до получим зависимость . (3. 60) Здесь: - число оборотов, совершенных телом за время движения от (правая ветвь параболы) до ; - средний период обращения за это же время; - отношение среднего периода обращения Земли за это же время к единице времени 1 сек (величина, равная , но безразмерная).
Рис. 1. Зависимости 1 - , 2 - и 3 - для Земли
Полное критическое число оборотов, совершаемых орбитальным телом будет . (3. 61) Здесь уместно сделать следующее замечание. Для расчета оценки числа оборотов на интервале значений эксцентриситета орбиты от до (левая ветвь параболы) и от до (правая ветвь параболы) необходимо знать не только значение эксцентриситета , но и знак градиента параметрических функций или . Если знак градиента этих функций отрицательный, то число оборотов рассчитывают по зависимости (3. 28). Если знак градиента этих функций положительный, то число оборотов будет, где - число оборотов, рассчитанное по зависимости (3. 28), но в пределах изменения параметрических функций от до и от до . Знак градиента функций и может быть определен только опытным путем. Время жизни космического тела на орбите
(сек), (3. 62)
или, современных земных лет
(лет), (3. 63)
где сек - современный период обращения Земли.
Среднее изменение периода обращения за год (на интервале)
или .
Но , следовательно
. (3. 64)
При и .
Это значение должно быть близким современным значениям планет и спутников, у которых значения начального эксцентриситета их орбит мало отличается от значения при , а современные значения эксцентриситета орбиты близки нулю. Строго говоря, значение Kв (3. 23), определяемое (3. 46), зависит от эксцентриситета орбиты Земли, но в пределах изменения эксцентриситета от до e = 0 может быть принято постоянным. При значительном изменении эксцентриситета
. (3. 65)
Учитывая, что
,
получим
, (3. 66)
Эта зависимость одинакова для всех орбитальных тел на интервале изменения эксцентриситета от до (правая ветвь параболы, рис. 1). Она же позволяет рассчитать изменение среднего периода обращения Земли на интервале от до , входящего в соотношения (3. 28), (3. 31). Как видно из рис. 1, градиент функции a= f(e), T= f (e), V= f (e) мал в области значений 0 < e < 0. 33 справа от e = 0). В области значений e >0. 33 справа от е = 0 градиент функций a= f(e) и T= f (e) возрастает, то есть при е > 0. 33 быстро увеличивается с каждым оборотом как большая полуось орбитального эллипса, так и период обращения периферийного тела. Из сравнения зависимостей (3. 36), (3. 37) и (3. 38) также следует, что наибольший темп изменения во времени имеют период обращения и большая полуось орбитального эллипса. Интенсивность полного излучения диполя, вращающегося в одной плоскости XY c постоянной угловой скоростью (средняя за период обращения) при V << c (3. 67) где - для одного заряда, движущегося по эллипсу с большой полуосью . При движении нейтрального тела по эллиптической орбите аналогом электрического зарядаq является величина , имеющая ту же размерность . Следовательно, интенсивность полного гравитационного излучения диполя, вращающегося с постоянной угловой скоростью при будет
(3. 68)
где сек/г - постоянная диполя.
Дипольное гравитационное излучение изменяется во времени, так как изменяются и T в соответствии с соотношениями (3. 33) и (3. 34). Подставив (3. 33) и (3. 34) в (3. 68), получим
(3. 69)
Максимальное значение I соответствует периоду обращения при е = 0, минимальное . Дипольное излучение поляризовано по эллипсу с отношением длин полуосей , где q - угол между направлением излучения и осью Z. В направлении оси Z дипольное излучение поляризовано по кругу. Энергия излучения гравитационного диполя
(3. 70)
В точке перехода от прямолинейного движения к орбитальному относительная скорость двух тел равнаи тело меньшей массы, кроме перехода к движению по эллиптической орбите, приобретает еще и вращение вокруг собственной оси. Поэтому уравнение энергетического баланса в точке перехода будет
(3. 71)
где М - масса орбитального тела,
- начальный момент инерции относительного фокуса
эллипса,
- начальные моменты инерции и угловая скорость
вращения тела вокруг собственной оси.
Уравнение эллипса в полярных координатах с центром в фокусе эллипса . (3. 72)
Следовательно
(3. 73)
Но
(3. 74)
В нашем случае n + 1 = 2, то есть n = 1, и (3. 74) вырождается в интеграл , (3, 75)
С учетом (3. 75)
(3. 76)
Известно, что
(3. 77)
Подставив (3. 76) и (3. 77) в (3. 71), получим
(3. 78)
и - (начальный период вращения). (3. 79)
Энергия, затрачиваемая на движение периферийного тела М по орбите, в соответствии с (3. 76) будет
где .
Энергия, затрачиваемая на вращение периферийного тела вокруг собственной оси , (3. 81)
где - средний радиус тела,
- период вращения.
Отметим, что орбитальный переход при гравитационном захвате нейтральных тел является своеобразным защитным барьером, препятствующим прямому проникновению свободно перемещающихся космических тел к центральному телу. Поверхности центрального тела могут достигнуть лишь периферийные тела, кинетическая энергия которых больше энергии орбитального перехода и направление вектора скорости которых по отношению к центральному телу пересекает его миделево сечение. 4. Уход тела с орбиты. Выясним, по какой параболической траектории будет двигаться тело после его ухода с орбиты. Уравнение этой траектории найдем из следующих соображений. Пусть тело совершает орбитальное движение по эллипсу. Начало координат разместим в центре эллипса (рис. 2). Уравнение траектории тела при орбитальном движении
(4. 1)
где и b - большая и малая полуоси эллипса.
Уравнение параболы с вершиной в перигелии эллипса и осью симметрии вдоль оси x будет , (4. 2)
где pў - параметр параболы.
Очевидно, что условие ухода тела с орбиты состоит в равенстве радиусов кривизны для эллипса и параболы в точке начала ухода, то есть в перигелии. Но радиус кривизны указанной параболы и эллипса в перигелии одинаков, то естьR = pў = p, где - параметр эллипса. Следовательно, уравнение траектории, по которой тело уходит с орбиты, будет
(4. 3)
В точке ухода тела с орбиты должен рассматриваться эллипс с критическим эксцентриситетом. Учитывая, что
(4. 4)
окончательно получим уравнение траектории ухода тела с орбиты (4. 5) При x = - y = 0; при x = 0 то есть на интервале траектория критического эллипса и траектория параболы совпадают (эллипс вписан в параболу). Таким образом уход тела с орбиты происходит в точке с координатамиx= 0, . При больших расстояниях от точки ухода с орбиты тело движется по траектории . (4. 6) Скорость ухода тела с орбиты, как это видно из рис. 2, равна Vкр .
Рис. 2. Схема ухода тела с орбиты.
В то же время скорость освобождения в точке ухода будет
(4. 7)
и направлена в сторону удаления от центрального тела вдоль линии, соединяющей точку ухода с фокусом эллипсаP. Из рис. 2 видно, что и, следовательно, скорость освобождения в точке ухода тела с орбиты (4. 8)
то есть (4. 9)
Полученный результат означает, что скорость ухода тела с орбиты по траектории (4. 5) меньше скорости освобождения. Следовательно, энергия освобождения периферийного тела, уходящего с орбиты по траектории (4. 5), в два раза меньше энергии освобождения того же тела, уходящего от тела большей массы по линии их центрального взаимодействия. Докажем, что в точке ухода тела с орбиты энергия вращения Esкр - малая величина. Если предположить, что сразу после точки ухода тело движется прямолинейно (а не по параболической траектории), то уносимая им кинетическая энергия (4. 10) так как в точке ухода (рис. 2) локальная скорость тела равна средней орбитальной скорости при движении тела по критическому эллипсу.
С другой стороны, потенциальная энергия в точке ухода
(4. 11)
где - расстояние между телами в момент ухода с орбиты тела массой M2. Сравнивая (4. 10) и (4. 11), получим то есть третий закон Кеплера для орбитального движения. Поэтому Esкр- малая величина, период вращения становится большим, то есть в момент ухода тела с орбиты оно практически не вращается вокруг собственной оси. Так как после ухода тела с орбиты гравитационный диполь перестает существовать, то интенсивность и энергия излучения диполя равны нулю. Выясним, как будет изменяться скорость тела после его ухода с орбиты. Дифференцируя (4. 5) поx, получим
(4. 12)
При x = 0 (4. 13)
Следовательно (4. 14)
Но в точке ухода V0 = Vкр и поэтому
(4. 15)
где Vx - скорость тела на расстоянии x от точки ухода.
Зная зависимость скорости тела от расстояния (после его ухода с орбиты), можно получить оценку среднего времени движения тела на этом же расстоянии. Обозначим Зависимость (4. 15) примет вид
где (4. 16)
Разделяя переменные в (4. 16) и интегрируя (при начальных условиях: t = 0 при r = 0), получим (4. 17)
При r >> 1 (4. 18)
Переходя к размерным величинам и учитывая, что получим
(4. 19)
то есть аналог третьего закона Кеплера.
Обратим внимание, что время свободного движения тела (после его ухода с орбиты) не зависит от массы тела, а определяется лишь пройденным расстоянием и массой центрального тела. Следовательно, время движения тел на одном и том же интервале движения неодинаково для разных центральных систем (принцип относительности времени движения периферийных тел в центральных системах). Среднее время движения на интервале [0, r ] при r >> 1 будет (4. 20)
или (4. 21)
Итак, мы знаем картину движения периферийного тела в начале орбитального перехода и после ухода с орбиты. Однако, не вполне ясно, по какой траектории будет двигаться тело при орбитальном движении. Как отмечалось, все параметры орбиты (большая полуось орбитального эллипса, период обращения, средняя орбитальная скорость, эксцентриситет) изменяются непрерывно от оборота к обороту, что отражают зависимости (3. 33 - 3. 38). Найдем средние значения этих функций на интервале изменения эксцентриситета орбиты от до = 0. Обозначим
(4. 22)
Тогда (4. 23)
При N ® 0 (среднеарифметическое значение на том же интервале). (4. 24)
При N ® 0
(4. 25)
При N ® 0
Далее обозначим
(4. 26)
(4. 27)
(4. 28)
Следовательно (4. 29)
(4. 30)
(4. 31)
При N ® 0
(4. 32)
Как видно, среднеинтегральные значения функций и на интервале от eн до e = 0мало отличаются от их среднеарифметических значений, что вполне оправдывает применение правил линейной интерполяции к этим функциям на этом интервале. Система из трех уравнений эллипса в пространственных полярных координатах r i, yi, zi (с центром в фокусе эллипса) (4. 33)
(где i = 1, 2, 3),
вместе с соотношениями
(4. 34)
и начальном условием ei = eнi при t= 0 определяет пространственную кривую, по которой движется орбитальное тело в системе координат с галактическим центром (здесь индекс 1 относится к движению Солнца, индекс 2 - к движению планет или комет, индекс 3 - к движению спутников планет). В этой системе пространственных полярных координат Солнце движется по эллипсу, расположенному в плоскости его эклиптики, планеты и кометы - по пространственной эллиптической спирали (элоиде) с изменяющимся эксцентриситетом каждого ее витка и переменным шагом между виткамиh = vTпл вдоль оси элоиды - траектории движения Солнца (орбитальная скорость движения Солнцаv , Tпл - период обращения планеты или кометы). Спутники планет движутся по аналогичным элоидам с переменным шагом между виткамиh = Vпл Tсп вдоль оси элоиды - траектории движения планеты (Vпл - орбитальная скорость планеты, Тсп - период обращения спутника). Описание такого вида траекторий движения космических тел достаточно сложно, поэтому ограничимся описанием траекторий орбитального движения планет вокруг Солнца (или, что то же, спутников вокруг планет). В этом случае траектория орбитального тела может быть приближенно представлена плоской эллиптической спиралью с изменяющимся от оборота к обороту эсцентриситетом и расположенной в плоскости эклиптики планеты (спутника). Уравнение каждого из эллипсов, последовательно составляющих эту спираль, в полярных координатахr, j(центр координат - в одном из фокусов эллипса, полярная ось направлена от фокуса к ближайшей вершине, рис. 3) будет
(4. 35)
где
Большая полуось эллипса
(4. 36)
Малая полуось
(4. 37)
Рис. 3.
Таким образом, уравнение (4. 35) описывает семейство эллипсов, практически переходящих один в другой и последовательно составляющих спираль (эллипсоиду) при изменении эксцентриситета эллипсов от до (рис. 3). Конечно, эксцентриситет орбиты изменяется и в течение одного оборота, но эти изменения столь незначительны, что ими в практических расчетах можно пренебречь. Рассмотрим основные характеристики эллипсоиды. Эллипсоида незамкнута и конечна, она состоит из двух ветвей (рис. 3): ветви, начинающейся в точке и заканчивающейся в точке и ветви, начинающейся в точке и заканчивающейся в точке В этой точке (или ей симметричной) тело уходит с орбиты и далее траектория его движения - парабола (4. 5). Тело может двигаться по эллипсоиде как по часовой стрелке, так и против нее в зависимости от того, какое направление оно приобретает в точке начала орбитального движения.
Характерные точки эллипсоиды:
вершины:
(4. 38)
,
сопряженные точки:
(4. 39)
Расстояние между фокусами каждого эллипса
(4. 40)
Расстояние между фокусом и ближайшей вершиной
(4. 41)
Площадь окружности r = 1 (при е= 0)
(4. 42)
Отношение площади каждого эллипса к площади окружности r = 1 (4. 43)
Длина окружности r = 1
(4. 44)
Отношение периметра каждого эллипса к длине окружности r = 1 (4. 45)
где - полный эллиптический интеграл 2 рода.
Директрисы каждого эллипса
(4. 46)
Угловые коэффициенты сопряженных диаметров каждого эллипса
(4. 47)
Радиус кривизны в любой точке каждого эллипса
(4. 48)
Радиус кривизны в характерных точках:
вершины А и В: R*=1
вершины С и D: (4. 49)
сопряженные точки и : . (4. 50)
Обратим внимание, что действительная траектория орбитального движения (элоида) никогда не проходит через одну и ту же точку пространства, а в сопряженных точках пересекаются лишь проекции витков элоиды (в эллиптических координатах с галактическим центром) на плоскость эклиптики планеты.
5. Эволюция планет и комет.
Применим полученные закономерности к анализу движения планет. В табл. 1 и 2 приведены характеристики движения планет. Признаем неоспоримо верными результаты астрономических наблюдений и вычислений (М, N, ak, Tk, Vk, ek) [3], а также величины, вычисленные из третьего закона Кеплера (aн, Tн, Vн, aо, Tо, Vо, aкр, Tкр, Vкр) с использованием соотношения (3. 4). Все остальные численные значения параметров, которые приведены в таблицах и тексте, назовем оценками, что, впрочем, не снижает их достоверность, а лишь указывает на ограниченную точность вычисления. Прежде всего отметим, что нейтральные тела относительной массы N< 0, 0314 подвержены гравитационному захвату. Эволюцию движения космических тел при гравитационном захвате можно представить в следующем виде. Первый этап прямолинейное движение тел навстречу друг другу отro до rmin (за время tmin). Второй этап - орбитальное движение тела меньшей массы вокруг тела большей массы. При этом эксцентриситет начальной орбиты для всех захватываемых тел определяется зависимостью (3. 4) и по мере движения уменьшается до e= 0. Третий этап - эксцентриситет орбиты изменяется в пределах 0 < e < eкр. Четвертый этап - эксцентриситет орбиты становится равным eкр, тело переходит от орбитального движения по эллипсу к движению по параболической траектории и освобождается от притяжения центрального тела. В соответствии с этой картиной эволюции движения нейтральных космических тел кометы - это планеты (или их остатки), прошедшие первые два этапа эволюции и находящиеся в настоящее время на ее третьем этапе (e ® eкр).
Табл. 2. Характеристики малых планет Солнечной системы
Малая планета
Параметр
Церера
Паллада
Юнона
Веста
Геба
Эвномия
Давида
М, гЧ1022
60
18
2
10
2
4
3
N, 10-10
3, 0166
0, 9050
0, 1005
0, 5027
0, 1005
0, 2010
0, 1508
ек
0, 0790
0, 2350
0, 2560
0, 0880
0, 2030
0, 1850
0, 1770
ак, смЧ1013
4, 1394
4, 1394
3, 9943
3, 5320
3, 6293
3, 9569
4, 7603
ro, смЧ1018
0, 8881
1, 5416
4, 4141
1, 8525
4, 1151
3, 1964
4, 4524
Tк, секЧ108
1, 4524
1, 4550
1, 3772
1, 1448
1, 1923
1, 3556
1, 7902
b m
119, 97
162, 11
230, 79
187, 77
280, 79
236, 14
253, 72
b min
113, 11
152, 84
264, 73
177, 03
264, 73
222, 63
239, 21
0, 0186
0, 0137
0, 0079
0, 0118
0, 0079
0, 0094
0, 0087
rm, 10-6
34, 3
19, 0
6, 40
14, 2
6, 4
9, 0
7, 8
rmin, 10-6
69, 6
38, 0
12, 8
28, 2
12, 8
18, 0
15, 6
ен
0, 3333
0, 3333
0, 3333
0, 3333
0, 3333
0, 3333
0, 3333
ан, смЧ1012
46, 2170
44, 4350
42, 375
38, 6520
39, 5050
43, 1510
52, 0930
tmin, секЧ1014
0, 5709
16, 3053
6, 3228
1, 7210
5, 6914
3, 8974
6, 4061
rmin, секЧ1012
61, 812
59, 580
56, 500
52, 269
52, 673
57, 535
69, 457
Umin, см. /сек Ч105
17, 59
18, 05
18, 50
19, 06
19, 14
18, 25
16, 62
Vк, , см. /сек Ч105
17, 90
17, 87
18, 22
19, 38
19, 12
18, 34
16, 70
Vн , см. /сек Ч105
16, 94
17, 25
17, 69
18, 53
18, 42
17, 56
15, 96
Tн, секЧ107
17, 134
16, 182
15, 049
13, 106
13, 540
25, 438
20, 494
DV, см. /сек Ч105
0. 96
0, 62
0, 53
0, 85
0, 70
0, 78
0, 74
Da, смЧ1012
4, 8230
3, 0410
2, 432
3, 332
3, 212
3, 592
4, 490
DT, секЧ107
2, 6110
1, 6320
1, 277
1, 658
1, 617
1, 882
2, 592
, секЧ107
15, 8290
15, 366
14, 410
12, 277
12, 731
14, 497
19, 198
n, обЧ106
5, 6538
3, 6360
3, 0359
4, 6277
4, 3547
4, 4554
4, 4485
n, секЧ1014
8, 9496
5, 5857
4, 3749
5, 6814
5, 5442
6, 4599
8, 5403
, лет 106
28, 3593
17, 6987
13, 863
18, 0033
17, 5683
20, 4669
27, 0623
екр
0, 999861
0, 999921
0, 999962
0, 999943
0, 999975
0, 999964
0, 99997
Vкр , см. /сек Ч105
0, 3046
0, 22781
0, 13434
0, 20726
0, 13910
0, 15783
0, 13372
Tкр, секЧ1010
2, 949410
7, 06356
34, 3938
9, 36427
30, 9731
21, 1772
34, 8435
, см. /сек Ч105
8, 6223
8, 7389
8, 9122
9, 3686
9, 2796
8, 8589
8, 0468
Малая планета
Параметр
Церера
Паллада
Юнона
Веста
Геба
Эвномия
Давида
, смЧ1016
7, 15315
12, 7868
36, 7637
15, 4477
34, 2944
26, 6388
37, 1059
, секЧ1013
1, 47470
3, 53178
17, 1969
4, 68214
15, 4865
10, 5886
17, 4217
акр, смЧ1016
14, 3017
25, 5693
73, 5233
30, 8917
68, 5850
53, 2733
74, 2067
nкр, обЧ108
34, 2756
34, 1065
34, 1521
34, 2752
33, 6764
34, 3221
34, 2946
nкр, секЧ1018
5, 05464
12, 04564
58, 7312
16, 0481
52, 1531
36, 3424
59, 7472
, лет 1012
16, 0160
38, 1698
186, 211
50, 8528
165, 261
115, 161
199, 325
ао, смЧ1012
41, 0818
39, 4978
37, 6667
34, 3573
35, 1156
38, 3564
46, 3049
Tо , секЧ107
14, 3592
13, 5614
12, 6119
10, 9835
11, 3472
12, 9379
17, 17504
Vо, см. /сек Ч105
17, 9676
18, 2964
18, 7631
19, 6540
19, 5374
18, 6252
16, 9281
Rпл , смЧ108
0, 350
0, 230
0, 110
0, 190
0, 110
0, 140
0, 130
Tsк , секЧ104
3, 270
-
2, 598
1, 9200
2, 622
2, 190
-
Tsн , секЧ104
0, 00941
0, 00600
0, 00280
0, 00473
0, 00276
0, 00363
0, 00369
0, 70643
0, 71108
0, 71078
0, 70109
0, 68260
0, 70716
0, 70797
cек/век
0, 1150
-
0, 1869
0, 1062
0, 1491
0, 1068
-
Отметим, что расстояние ro, при котором начинается захват периферийного тела центральным, не прямопропорционально относительной массе планет, а сложным образом зависит от их относительной массы и будущих параметров орбиты периферийного тела. Порядок этого расстояния для планет составляет 1015 - 1013 см, что намного меньше расстояния от Солнца до галактического центра, но соизмеримо с расстояниями от Солнца до ближайших звезд. Конечно, результаты расчета ro для планет и других нейтральных космических тел носят характер оценки, поскольку при постановке задачи о взаимодействии тел в центральном поле было принято [1], что в момент начала захвата скорость периферийного тела равна нулю. В действительности любое свободно движущееся космическое тело в этот момент имеет конечную скорость относительно центрального тела. На рис. 4 - 6 показаны зависимости (2. 27 - 2. 30) для планет. Сплошные линии - указанные предельные зависимости, точки на рисунках - результаты расчета этих же характеристик по данным астрономических наблюдений [3]. Видно, что каждая из планет располагается на рисунке в строгом соответствии с ее относительной массой. Как показано в [1] , зависимости (2. 27 - 2. 30) являются предельными и для других космических объектов (ближайшие галактики, звезды, спутники планет), подверженных гравитационному захвату. Таким образом, характеристики движения большого многообразия космических объектов удается упорядочить в виде единых зависимостей динамических (bm, bmin) и кинематических (rm, rmin) характеристик от относительной массы взаимодействующих тел (или их ассоциаций). Еще раз подчеркнем, что зависимости (2, 27) - (2. 30), также как (3. 8), (3. 9) и (2. 21) следуют из теории гравитационного захвата, а для сравнения с ними использованы опытные данные астрономических наблюдений и расчетов характеристик реальных космических объектов [3]. Эти. же зависимости могут быть использованы для идентификации искусственных или еще не захваченных: космических объектов. В этих случаях характеристики движения объектов не будут соответствовать зависимостям при гравитационном захвате тел. Заметим, что взаимное расположение планет на рис. 4 - 6 отличается от их расположения относительно Солнца. Это неудивительно, поскольку на рис. 4 - 6 используются координаты ln b = f(ln N) и ln r = f(ln N). При использовании реальных значений rm = rm ro, rmin=rmin ro, aн=ro и полярных координат эти начальные характеристики орбит планет соответствуют их современному расположению относительно Солнца. Рис. 4. Зависимости: I - ln bm = f( ln N) и II - ln bmin= f( ln l) для планет солнечной системы: М - Меркурий, Мс - Марс, В - Венера, П - Плутон, З - Земля, У - Уран, Н - Нептун, С - Сатурн, Ю - Юпитер. Номера присвоены планетам в порядке увеличения среднего расстояния от Солнца. Рис. 5. Зависимости: I - ln rmin = f( ln N) и II - ln rm= f( ln N) и III - ln = f( ln N) для планет солнечной системы; обозначения те же, что на рис. 4. Рис. 6. Зависимости: I - ln bm = f( ln N) и II - ln rmin= f( ln N) для малых планет солнечной системы: Г - Геба, Д - Давида, Э - Эвномия, В Веста, П - Паллада, Ц - Церера. Номера присвоены планетам в порядке увеличения среднего расстояния от Солнца. Начальные, современные и критические параметры, определяющие движение планет, представлены в табл. 1 и 2. Прежде всего отметим, что время жизни планет на орбитах и суммарное время с начала их гравитационного захвата для всех планет не превышает космогонического возраста Земли (4, 5109лет). Время от начало захвата планет до их перехода на орбиту намного меньше времени их орбитального движения. Полное время жизни планет на орбите составляет от 5, 31010 до 9, 11011 современных земных лет. Это время намного больше современной оценки времени существования Вселенной (24109 лет). Большие времена жизни планет на орбите вовсе не означают, что планеты (или иные космические тела) смогут существовать в течение этого времени. Рассчитанное время жизни планет (космических тел) лишь показывает возможность, ресурс их орбитального движения. Реальное время существования космического тела на орбите определяется многими космологическими обстоятельствами и, в первую очередь, временем жизни центрального тела орбитальной системы. Результаты расчетов показывают, что планетная система не образовалась в одно и тoже время, а постепенно формировалась по мере захвата Солнцем нейтральных тел. Оценки параметра Тi (время жизни планет от перехода на орбиту до современности) позволяют установить последовательность захвата планет во времени). Из табл. 1 видно, что этот порядок следующий: Нептун, Плутон, Уран, Сатурн, Юпитер, Марс, Земля, Венера, Меркурий. В промежутке времени между захватом Юпитера и захватом Марса были захвачены малые планеты. Таким образом, время от начала захвата планет составляет 1, 11109 лет. При анализе орбитального движения планет нас интересует прогноз: . в течение какого времени на планетах сохраняются примерно одинаковые космологические условия (расстояние от Солнца, период обращения)? Как видно из рис. 7 - 9 и табл. 1 и 2, достаточно близкие космологические условия сохраняются на планетах в течение времени эволюции их орбит от eн до e = 0 и далее, от e = 0 до e = к eн к . Действительно, на этом временном интервале как большие полуоси планет, так и их периоды обращения несущественно отличаются. Следовательно, этот временный интервал для любого космического тела может быть оценен как 2nо/Е. При последующем увеличении эксцентриситета орбиты планет космологические условия будут все более заметно отличаться от начальных. Рис. 7. Зависимости для планет солнечной системы: 1 - Меркурий, 2 - Марс, 3 - Венера, 4 - Плутон, 5 - Земля, 6 - Уран, 7 - Нептун, 8 - Сатурн, 9 - Юпитер. Рис. 8. Зависимости для планет. Обозначения те же, что на рис. 7. Рис. 9. Зависимости для планет. Обозначения те же, что на рис. 7. Поскольку современные эксцентриситеты орбит всех планет не превышают eн 0, 33, то можно с уверенностью полагать, что космологические условия на планетах от начала их орбитального движения до наших дней сохранялись близкими к начальным. Отметим, что космологические условия сохранялись постоянными и в дозахватный период существования тела во время его свободного движения в космическомпространстве. Однако эти условия по их физическим характеристикам существенно отличались от космологических условий при орбитальном движении тел. Зная разность Тн - Тo и число оборотов no, можно рассчитать среднее уменьшение периода обращения dТo за 1 оборот космического тела (планеты). Значения dТo, приведенные в табл. 1 и 2, должны быть близкими современным значениям dTo, поскольку на интервале Тн- Тo градиент функции Т = f (e) незначителен (рис. 8). В табл. 3 представлена интенсивность полного излучения гравитационного диполя планет. Видно, что дипольное излучение планет значительно по величине (10151022 ). Для сравнения укажем, что cовременное полное излучение электрического диполя Земли при ее заряде q = - 5, 710 к составляет I =1, 8 10 . С дипольным излучением периферийных тел непосредственно связана периодическая (со временем, равным периоду обращения) деформация объема центральных тел. Так дипольное излучение планет влияет на деформацию объема Солнца. Как видно из табл. 3, наибольшее влияние на периодическую деформацию объема Солнца оказывает дипольное излучение Венеры, Меркурия, Земли и Юпитера с цикличностью, соответствующей периодам обращения этих планет. Влияние остальных планет на деформацию объема Солнца мало, хотя и им, конечно, соответствуют определенные циклы деформации объема. Поскольку дипольное излучение поляризовано по эллипсу, то наибольшая деформация объема Солнца происходит в момент прохождения планетой ее перигелия. Та же самая картина соответствует движению спутников вокруг своих планет. Так, дипольное излучение Луны приводит к периодической деформации объема Земли, что проявляется как в деформации земной коры, так и в приливах и отливах океанов. Обратим внимание на малую величину среднего (за период обращения) дипольного излучения Луны, по сравнению, например, с дипольным излучением спутников Сатурна Мимаса, Энцелада, Тефии, Дионы, Реи и Титана и, тем более, со спутниками Юпитера. По-видимому, периодическая деформация объема Сатурна и Юпитера существенно превышает деформацию объема Земли. Величина Е = Мс есть энергия, запасенная одним телом. Энергия покоя Е = Мсb характеризует полную энергию каждого из двух тел, взаимодействующих и совершающих движение в центральном поле. Центральное тело в системе двух тел, в свою очередь, может быть периферийным в другой системе тел (например, Сoлнце является центральным телом нашей планетной системы и одновременно периферийным телом в системе с галактическим центром).
Величина
e = = N
характеризует соотношение между энергией покоя периферийных тел в каждой из взаимодействующих центральных систем. Для планет солнечной системы e 10 до 10 (табл. 3). Поскольку энергия покоя включает в себя все виды энергии, в том числе энергию движения и излучения, то представляет интерес рассмотреть соотношение между этими основными составляющими энергии на примере нашей планетной системы. Результаты расчета энергетических характеристик планет представлены в табл. 3. Как видно, энергетические составляющие движения и излучения периферийного тела намного меньше энергии покоя. Основную долю энергии составляет энергия орбитального движения, на вращение планет и их гравитационное излучение затрачивается значительно меньше энергии (за исключением планет - гигантов, для которыхэнергия, затрачиваемая на вращение, соизмерима с энергией их орбитального движения). Так как известны начальные (в момент перехода от прямолинейного к орбитальному движению) и критические (в момент ухода планеты с орбиты) характеристики планет, то можносопоставить и их энергетические характеристики. Результаты расчета начальных и критических энергетических характеристик планет приведены в табл. 3. Видно, что начальная энергия орбитального движения меньше современной (также как и критическая энергия орбитального движения). Следовательно, энергия орбитального движения имеет максимум (при е = 0), что соответствует изменению во времени средней скорости орбитального движения (рис. 9). Начальная энергия излучения гравитационного диполя планет меньше, чем современная (также как и критическая энергия излучения/. Следовательно, энергия излучения гравитационного диполя имеет максимум. Отношение и за время орбитального движения мало изменяется. Энергия, затачиваемая на вращение планет, наибольшая в начальный момент (то есть в момент перехода от прямолинейного движения к орбитальному), при этом отношение 2, 451, 67; в дальнейшем отношение уменьшается. Для планет - гигантов , то есть за период от момента перехода на орбиту до настоящего времени их энергия, затрачиваемая на вращение, мало изменилась. Перейдем к анализу движения комет. Как отмечалось, кометы - это планеты, находящиеся на третьем этапе эволюции (), то есть наиболее долгоживущие орбитальные тела планетной системы. Поэтому анализ движения комет может привести к оценке времени жизни нашей планетной системы. Масса комет намного меньше массы Солнца, поэтому N << 1 для оценки характеристик движения комет можно принять, что их начальный эксцентриситетeн 0, 333. Современные большие полуоси орбит комет и эксцентриситеты орбит известны [3]. Следовательно, можно рассчитать и остальные характеристики орбитального движения комет. Результаты расчета приведены в табл. 4. Для расчета числа оборотов n комет на интервале от eн до использовалось соотношение (3. 59), в котором определяется соотношением (3. 66). Эта же зависимость вместе с зависимостью показана на рис. 10. При зависимость плавно переходит в зависимость на интервале движения от до . Полное число оборотов кометы от начала ее перехода на орбиту по будет , где - число оборотов, совершенных кометой за время движения от (справа от , рис. 1) до , - число оборотов, совершенных кометой за время движения от eн до e = 0.
Рис. 10. Зависимость nk = f(e) для комет солнечной
системы (1) и зависимость Ke = f(e), (2).
Как видно, полное число оборотов, совершенных кометой с наибольшим эксцентриситетом орбиты (комета Галлея, = 0, 967) составляет . Отсюда следует оценка времени жизни кометы Галлея на орбите: лет. Она же является оценкой возраста нашей планетной системы. По порядку величины эта оценка согласуется с принятой в настоящее время оценкой возраста Вселенной (лет), но представляется более достоверной, так как получена из анализа движения достаточно близких к Земле и потому тщательно исследованных космических тел (комет). В целом этапы эволюции планет и комет, рассмотренные на примере Солнечной системы, характерны и для других нейтральных тел, захваченных массивной звездой. Таким образом, исходя из теории гравитационного захвата нейтральных тел, существенно увеличивается вероятность существования планетных систем воВселенной.
6. Эволюция спутников планет.
Поскольку движение спутников относительно планет подчиняется тем же закономерностям гравитационного захвата, что и движение планет относительно Солнца, то все полученные зависимости применимы и к анализу движения спутников планет. Рассмотрим основные особенности движения спутников планет. Как и для планет, для каждого спутника характерны те же этапы его эволюции: прямолинейное движение на интервале, орбитальное движение, уход с орбиты и движение по параболической траектории. Сравнение данных астрономических наблюдений с предельными зависимостями (2. 21), (2. 27 - 2. 30) и (3. 8) представлено на рис. 11. Видно, что на участке прямолинейного движения опытные данные согласуются с соответствующими зависимостями. (Исключение составляет Луна: при расчетах отмечены отклонения по ряду параметров от общих закономерностей гравитационного захвата нейтральных тел). Орбитальное движение спутников планет и уход с орбиты подчиняется тем же закономерностям, что и для планет. Результаты расчета начальных, современных и критических величин, определяющих движение спутников планет, приведены в табл. 4. Величина для спутников планет составляет cм, что соизмеримо с параметрами орбит планет. Современное время жизни на орбите спутников планет намного меньше времени жизни самих планет. То же самое относится к сравнению критического времени жизни спутников планет и самих планет. Поскольку критическое время жизни спутников планет намного меньше времени жизни самих планет, то завсе время жизни на орбитах возможен переход спутника от одной планеты к другой планете (или спутник может стать планетой). Рис. 11. Зависимости для спутников планет солнечной системы: I - ln = f(lnN), II - ln bmin = f( ln N ), III - ln bm = f( ln N), IV - ln rmin = f(ln N), V - ln = f( ln N), VI - ln rm = f( ln N), VII - ln eн = f( ln N). 1 - Мимас, 2 - Энцелад, 3 - Гиперион, 4 - Нереида, 5 - Тефия , 6 - Диона, 7 Рея , 8- Япет , 9 - Европа , 10 - Ио , 11 - Каллисто , 12 - Ганимед, 13 - Титан , 14 - Тритон , 15 - Луна; точки - расчетные значения характеристик, сплошные линии - зависимости (2. 21), (2, 27-2. 30), и (3. 8) и усредняющая кривая (зависимостьVII). Результаты расчетов показывают, что системы планетарных спутников не образовались в одно и то же время, а формировались по мере захвата планетой нейтральных тел. Это отчетливо видно на примере системы спутников Сатурна, а именно, последовательность захвата спутников во времени такова: Япет, Гиперион, Титан, Рея, Диона, Тефия, Энцелад, Мимас. Аналогичная закономерность отмечается и для спутников Юпитера. Как и для планет, на спутниках сохраняются достаточно близкие космологические условия (большая полуось орбиты, период обращения) в течение времени эволюции их орбит. При дальнейшем увеличении времени жизни спутника космологические условия для него будут все более заметно отличаться от начальных вследствие все более резкого изменения эксцентриситета орбиты. В табл. 6 представлена интенсивность полного излучения гравитационного диполя спутников планет. Видно, что дипольное излучение спутников планет составляети для спутников Юпитера Ио и Европа сравнимо по величине с дипольным излучением Земли. С дипольным излучением спутников планет связана периодическая (со временем, равным периоду обращения) деформация объема планет. Так, на деформацию объема Сатурна наибольшее влияние оказывают его спутники Тефия, Титан, Диона, Рея. и Мимас, а на деформацию объема Юпитера - его спутники Европа и Ганимед. Конечно, и остальным спутникам соответствуют определенные циклы деформации объема этих планет, причем наибольшая деформация объема планеты происходит в моментпрохождения спутником его перигелия. Результаты расчетов энергетических характеристик спутников планет представлены в табл. 6. Как видно, для спутников планет. Эта величина превышает соответствующие величины для планет. Полная энергия спутников планет существенно превышает их энергию орбитального движения и энергию излучения. Отношение энергии излучения гравитационного диполя к энергии, затрачиваемой на орбитальное движение, для спутников того же порядка, что и для планет. Табл. 6 Энергетика спутников планет
E
e
IK
EIK
EPK
Спутник
10-4
Луна
18, 376
2, 3634
0, 000441
0, 1125
17, 382
Мимас
3, 3078
0, 0430
4, 370
35, 58
1, 8566
Энцелад
5, 1353
0, 0631
2, 967
35, 13
2, 784
Тефия
13, 695
0, 1751
11, 68
190, 54
21, 77
Диона
17, 253
0, 2206
6, 871
162, 52
25, 87
Рея
26, 40
0, 3376
4, 030
157, 36
41, 43
Титан
204, 56
2, 6159
8, 280
1142, 4
1031
Гиперион
5, 622
0, 0719
0, 00308
0, 5673
0, 623
Япет
38, 12
0, 4875
0, 00409
2, 824
12, 75
Ио
280, 64
1, 911
12560
192062
5933
Европа
216, 32
1, 473
1185
36388
2252
Ганимед
390, 06
2, 656
583, 8
36147
45, 76
Каллисто
297, 93
2, 029
35, 59
5151
1505
Тритон
95, 173
2, 917
0, 00075
0, 0381
715, 4
Нереида
1, 5146
0, 046
Как отмечалось, некоторые расчетные характеристики гравитационного захвата Луны не соответствуют закономерностям гравитационного захвата нейтральных тел. Из рис. 11 видно, что такие характеристики гравитационного захвата Луны, как и незначительно отклоняются от общих закономерностей. Существенно отклонение от этих закономерностей таких характеристик, как. Расчет характеристики для Луны приводит к абсурдному результату: отношение , чего не может быть, так как отношение средней скорости к максимальной на интервале движения должно быть меньше единицы (что и подтверждается расчетами этой же характеристики для остальных спутников планет солнечной системы). Отличия характеристик гравитационного захвата Луны от общих закономерностей связаны со следующим обстоятельством. Внимательный читатель обратил внимание на приближенный характер зависимости (2. 13), определяющей верхнюю границу орбитального перехода при гравитационном захвате нейтральных тел (). Точное соотношение будет , (6. 1)
так как .
Для Луны N = 0, 0122954 (табл. 4) и, следовательно, = 0, 3911656. Решая неравенство (6. 1) при этом значении , получим N = 0, 0145. Таким образом, переход Луны на орбиту вокруг Земли совершился в условиях, близких к верхней границе, определяющей орбитальный переход. Этим, по-видимому, и обусловлены некоторые отклонения характеристик гравитационного захвата Луны от общих закономерностейгравитационного захвата нейтральных тел. 7. Эволюция Земли. При рассмотрении эволюции Земли главная особенность, по сравнению с другими планетами, состоит в том, что нам известен ее космогонический возрастлет [3]. Это позволяет более определенно судить об отдельных временных этапах ее эволюции, в том числе эволюции жизни на Земле. Согласно расчетам (табл. 1), расстояние от Солнца , на котором начался гравитационный захват Земли, составляло 3, 232см (2, 16 а. е. ). Время от начала захвата Земли до перехода на орбиту лет, а время жизни на орбите лет. Поскольку космогонический возраст Земли много больше суммы этих величин, то большую часть своей жизни планета Земля провела, свободно перемещаясь в космическом пространстве вне сфер притяжения массивных космических тел. Предположение о том, что планета Земля могла быть какое-то время в составе другой планетной системы, маловероятно, поскольку ее полное время существования на орбителет. Это время намного больше современной оценки времени жизни Вселенной. По-видимому, и для остальных планет и других космических тел, захваченных Солнцем, характерен длительный этап их свободного движения в космическом пространстве. Для решения этой проблемы необходимо знать космогонический возраст планет и иных космических тел. Теория гравитационного захвата нейтральных тел практически отвергает гипотезу о происхождении планетных систем вследствие выброса части материи из недр звезд. В самом деле, маловероятен выброс части материи из недр звезды на расстояния, намного превышающее ее поперечный размер. За время движения по орбите Земля совершила оборотов, а время ее жизни составило лет. Современный эксцентриситет орбиты Земли = 0, 01675, то есть близок к е = 0. Нас, как и ранее, интересует прогноз: сколько времени сохранялись и будут сохраняться на Земле космологические условия, близкие современным? Исходя из рис. 7 и рис. 8, можно утверждать, что космологические условия, близкие современным, сохранялись на Земле от начала ее орбитального движения до наших дней и будут сохраняться еще в течение лет. При дальнейшем увеличении эксцентриситета орбиты Земли космологические условия будут все более существенно отличаться от современных. Отличия эти будут тем значительнее, так как во время прохождения Землей афелия она будет намного дальше от Солила, чем в настоящее время, что существенно уменьшит интенсивность облучения Солнцем поверхности Земли. Кроме того, по мере увеличения эксцентриситета орбиты Земли, ее вращение вокруг собственной оси будет все более замедляться. Все эти изменения космологических условий приведут к существенным изменениям климата Земли. Заметим, что расчеты изменения космологических параметров во времени позволяют создавать математические модели изменения климата Земли в прошлом и будущем, учитывающее, конечно, изменения, происходящие на планете за это же время (состав атмосферы, радиоактивный нагрев и т. д. ). Что касается общего климатического прогноза на будущее, то он крайне неблагоприятен как для Земли, так и для планет солнечной системы с точки зрения обеспечения необходимой экологической жизненной ниши. Действительно, при увеличении эксцентриситета орбит планет период, в течение которого планета движется вблизи перигелия орбиты, будет все более уменьшаться и наоборот, период, в течение которого планета движется на значительно большем расстоянии от Солнца, чем в перигелии, будет все более увеличиваться. При этом, как показывают расчеты (табл. 1), расстояние от планеты до Солнца даже в критическом перигелии будет такого же порядка, как современные большие полуоси их орбит. В связи с тем, что климатические условия на Земле и планетах будут со временем становиться все более жесткими, центральной проблемой обеспечения необходимой экологической жизненной ниши станет создание искусственного климата с использованием искусственных и геотермальных источников энергии. Рассмотрим некоторые особенности строения Земли в древнейшие времена и в настоящее время. По космологическим понятиям Земля и ранее и в настоящее время представляла собой шар с весьма малым твердым ядром, окруженным жидким расплавом - мантией, на поверхности которой из-за охлаждения Земли в космосе образуется тонкая твердая оболочка (земная кора). Эта тонкая оболочка деформируется под влиянием притяжения Солнца. Главной особенностью строения Земли в древнейшие времена следует признать существование единого массива суши, который впоследствии расползался на отдельные материки (дрейф материков). Такое неравномерное распределение суши по поверхности Земного шара необъяснимо в рамках современных космологических гипотез, но оно легко объясняется, если придерживаться теории гравитационного захвата. В самом деле, на участке прямолинейного движения от до (в течение времени ) любое космическое тело, в том числе и Земля, не вращается, но испытывает притяжение центрального тела. Под влиянием этого притяжения та часть тонкой оболочки планеты (коры), которая обращена к центральному телу (Солнцу), прогнется в его сторону, так же как и противоположная ей часть оболочки. Благодаря этому на одной стороне шара образуется выпуклость поверхности (коры), в будущем - суша, а на другой стороне шара - вогнутость пoверхности; в будущем эта вогнутость заполнится сконденсировавшимся паром (на Земле - водой), образуя древнейший океан. Конечно, после перехода космического тела (Земли) к орбитальному движению с вращением вокруг собственной оси, притяжение центрального тела и кориолисово ускорение будут совсем по другому воздействовать на тонкую оболочку периферийного тела, что неизбежно приведет к изменению фигуры планеты и вида ее поверхности. Рассмотренная картина эволюции движения Земли не противоречит возможности возникновения живой материи. В самом деле, при свободном движении Земли или иного тела. в космосе в течение длительного времени тепло, требующееся для образования из неорганических элементов органических соединений и впоследствии живой материи, может поступать к ее поверхности как из недр Земли и как продукт радиоактивного распада элементов. Конечно, по мере остывания космического тела (уменьшения теплового потока из недр и от радиоактивного распада) все менее вероятным становится возникновение живой материи. В дальнейшей эволюции органических соединений и живой материи решающую роль играет переход космического тела на орбиту вокруг относительно горячей звезды, которая за счет излучения поддерживает температуру космического тела на уровне, достаточном для возникновения и существования живой материи. При этом определенно нельзя сказать, возникли ли какие-либо формы живой материи до перехода Земли на орбиту или уже при орбитальном движении. По своему влиянию на процессы возникновения и развития живой материи доорбитальный и орбитальный этапы движения имеют два существенных отличия: при орбитальном движении на всю массу орбитального тела действует кориолисово ускорение, а на поверхность орбитального тела воздействует излучение горячей звезды. Все остальные действующие факторы остаются такими же, как и при свободном движении тела в космосе. Конечно, на этапе прямолинейного сближения тел (от начала захвата тела до его перехода на орбиту) изменяются некоторые физические характеристики захватываемого тела (например, освещенность поверхности и количество тепла, поглощенного телом), но эти изменения действуют сравнительно недолго, так как время прямолинейного движения тел намного меньше времени их орбитального движения.
ЛИТЕРАТУРА
Сабелев Г. И. К движению двух тел в центральном поле. М. , ИВЦ “Маркетинг”, 1995. Ландау Л. Д. и Лифшиц Е. М. Механика. М. , ГИФМЛ, 1958. Таблицы физических величин. Справочник под ред. Кикоина И. Н. М. , Атомиздат, 1976.
|
