Н.И. Перов

Введение

На основании 60-летних наблюдений астрономов Главной астрономической обсерватории РАН периодических изменений блеска системы Плутон-Харон, происходящих с периодом 7,8 лет [3], в работе [2] , с использованием модели движения частицы в рамках дважды осреднённой ограниченной задачи трёх тел, с учётом сжатия Плутона, было сделано предположение о существовании вблизи Плутона на расстоянии 2510-2520 км гипотетического кольца, состоящего из частиц с характерным радиусом Rr. Поскольку радиус Плутона составляет RPl =1137 км [4], то естественно ожидать, что атмосфера Плутона будет ограничивать космогонический срок жизни такого кольца. Для состава атмосферы Плутона и величины давления на его поверхности получены приближённые оценки. В частности, в работе [4] утверждается, что атмосфера Плутона состоит в основном из азота, окиси углерода и примеси метана, а давление (p0) на поверхности планеты составляет всего несколько микробар при температуре (T) 35-45 K.

Поставим задачу определения космогонического времени жизни (τ) спутникоподобного объекта радиуса Rr, плотности ρr, который в начальный момент времени находится на круговой орбите на расстоянии rr от центра Плутона и движется в атмосфере планеты. Примем, что плотность атмосферы Плутона (ρ), при неизменной температуре, изменяется по экспоненциальному закону, при изменении ускорения свободного падения (g) с высотой. Вследствие торможения радиус орбиты объекта будет уменьшаться, и по истечении времени (τ) этот объект упадёт на поверхность Плутона (поверхность Плутона на 70% состоит из скальных пород и на 30% - из льда).

Основные уравнения

Известно [1], что объект, движущийся со скоростью V относительно окружающей его атмосферы, подвержен воздействию аэродинамических сил, которые можно разложить на две компоненты: сопротивление D, действующее в направлении, противоположном V, и силу в плоскости, перпендикулярной V. Выражение для D записывают в аэродинамике [1] обычно в виде

ρVVSCD , (1)

где

ρ - плотность окружающей среды,

S - характерная площадь объекта (обычно - это площадь сечения объекта плоскостью, перпендикулярной направлению движения).

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

Значение CD примем равным 2,2 в соответствии с монографией [1] , где показано, что для искусственных спутников Земли самой разнообразной формы и находящихся на различных высотах над поверхностью Земли CD (2,1 ; 2,35).

Плотность атмосферы Плутона зависит от многих факторов (расстояния Плутона от Солнца, солнечной активности, температуры дневного и ночного полушарий Плутона). Учитывая изменение ускорения свободного падения (g) с высотой (h), плотность (ρ) изотермической атмосферы Плутона представим в виде

, (2)

где μ - молярная масса (азота), RPl - радиус Плутона, - универсальная газовая постоянная, G - гравитационная постоянная, MPl - масса Плутона.

Величина V в формуле (1) в общем случае является скоростью движения объекта относительно окружающей среды (атмосферы Плутона). Если обозначить через v скорость движения объекта относительно центра масс Плутона, а через VA - скорость движения атмосферы планеты (вследствие осевого вращения Плутона с периодом PPl), то

V=v-VA (3)

С учётом возмущений от сопротивления атмосферы Плутона дифференциальное уравнение, описывающее орбитальную эволюцию частицы массы mr гипотетического кольца планеты, примет следующий вид

d2r/dt2 = - rr + D/mr (4)

В уравнении (4) второе слагаемое - ускорение, обусловленное торможением в среде - должно быть существенно меньше первого члена - гравитационного ускорения. При малых значениях величины |D/mr| скорость |v| приближённо будет совпадать со скоростью движения объекта по невозмущённой орбите

(5)

Здесь a - большая полуось орбиты частицы-спутника Плутона.

Рассматривая круговые (rr=a), экваториальные орбиты объектов (лежащие в плоскости экватора Плутона), выражая h через a и RPl,

h=a-RPl, (6)

Выражая S и mr через Rr и ρr (для сферических частиц) ,

S =πRr2; mr=Rr3ρr, (7)

От уравнений (1) - (7) (дифференцируя (5) по времени (t) и используя (4)) , перейдём к уравнению изменения большой полуоси орбиты (a) с течением времени (t) - da/dt

(8)

В уравнении (8) знак "-" соответствует одинаковым направлениям скоростей объекта и атмосферы Плутона, а знак "+" - противоположным . Для больших полуосей орбит мелких спутников Плутона, удовлетворяющих условию a10 м) и плотные (ρr>1000 кг/м3 ) тела выпадают на поверхность планеты.

Разделяя в уравнении (8) переменные a и t, используя новую переменную ξ=1/ и интегрируя в пределах от ξн=6,299·10-4 м-1/2 до ξк=8,648·10-4 м-1/2, получим окончательную формулу для определения времени жизни τ объектов, образующих кольцо Плутона

(9)

Таким образом, при сделанных предположениях, время жизни спутникоподобных тел вблизи Плутона прямо пропорционально произведению радиусов этих тел и их плотности и обратно пропорционально давлению у поверхности планеты. Это время жизни также зависит сложным образом от других параметров, как атмосферы (μ, T), так и планеты (MPl, RPl).

ПРИМЕРЫ:

Для численных оценок определим время жизни: 1) "короткоживущих" и 2) "долгоживущих" частиц в атмосфере Плутона. Примем следующие значения параметров Плутона и его атмосферы [4]:

G=6,672·10-11 Нм2/кг2,

μ=28·10-3 кг/моль,

RPl=1137000 м,

MPl=1,27·1022 кг,

=8,31 Дж/(К·моль),

aн=2520000 м, ξн=6,299·10-4 м-1/2,

aк=1337000 м, ξк=8,648·10-4 м-1/2,

СВ=2,2.

Результаты вычислений для различных частиц и условий приведены в таблице (I - значение интеграла (9), вычисленного по методу Симпсона при 2m=200)

N	T К	p0 Па	ρr кг/м3	Rr м	I м1/2	Τ годы
1	45	1	100	0,1	3,47·10-10	3,31·105
2	35	0,1	1000	1	2,05·10-13	1,28·1012

Из вышеизложенного следует, что время жизни крупных (Rr>10 м) и плотных (ρr>1000 кг/м 3) частиц в атмосфере Плутона превосходит на несколько порядков время жизни (τ) мелких (Rr