ПБС

[Старый]

Видно, что быстрее 0.06 км/день - почти только Deb. Не вся таблица, не влезла, менее 0.06 - 79 объектов.

Я же и говорю: спутники надувные! 

[cross-track]

OPS5712, он же NRL PL 160, он же Калсфер-4, был выведен на орбиту чуть повыше - 919х934 км. Сейчас его орбита 637х640 км и он потерял за год 22 км. На высоте 640 км, да

22 км/год - видимо от буйного Солнца. В среднем за 59 лет у него было 4.9 км/год. А сейчас за 10 дней 0.037 км/сут ( 13.6 км/год) , а мгновенное по текущему TLE 0.028 км/сут (10.3 км/год).

В таблице - объекты со средней высотой 620..640 км и разницей апогей-перигей не более 20 км.
Видно, что быстрее 0.06 км/день - почти только Deb. Не вся таблица, не влезла, менее 0.06 - 79 объектов.

Насколько я понимаю, в двух последних столбцах приведено вычисление темпа спада высоты двумя методами. Видимо, в левом столбце вы вычисляете dh/dt по прямым измерениям (т.е. берете отношение изменения высоты к интервалу времени, например, 12 дней), а в правом столбце находите  другим методом - с использованием значений n, dn/dt из ТЛЕ, причем тоже, видимо, усредняете эти величины по этому же интервалу времени, а затем пересчитываете эти усредненные зачения в dh/dt. При этом получаются близкие, но все же отличающиеся значения. Вы не могли бы для какого-то объекта из этой таблицы привести численные значения средних n, dn/dt, которые вы использовали как входные данные в расчетах?

[Олег]

с использованием значений n, dn/dt из ТЛЕ, причем тоже, видимо, усредняете эти величины по этому же интервалу времени, а затем пересчитываете эти усредненные зачения в dh/dt

Нет, не усреднял. Брал из второго TLE. Хотя логичнее наверное было взять среднее по двум TLE.
Код для расчета :

Спойлер

function dndtTodadt(dndt,n:Double):Double; //  dn/dt, n => da/dt =dh/dt
const KNA = 42241.1083269817;
begin
Result:=-(2/3)*KNA*Power(n,-5/3)*dndt;  // Power - n в степени минус 5/3
end;

[свернуть]

*  в TLE приведено значение (dn/dt) / 2
Примеры
0 CAPELLA-11 (ACADIA)
1 57693U 23126A  26121.45702955  .00005981  00000-0  78261-3 0  9993
2 57693  53.0094  27.2233 0000736  13.8846 346.2151 14.81562452145198
0 CAPELLA-11 (ACADIA)
1 57693U 23126A  26131.97901696  .00002440  00000-0  33263-3 0  9994
2 57693  53.0101 341.7351 0001340 123.2065 236.9041 14.82016278146750

n	dn/dt	-da/dt = -dh/dt
14.8156245	0.00011962	0.03769
14.8201628	4.88E-05	0.01537

0 OPS 5712 (P/L 160)
1  2826U 67053A  26121.53418036  .00007816  00000-0  11151-2 0  9994
2  2826  69.9154 152.5825 0002011 247.3763 112.7172 14.77069011 59346
0 OPS 5712 (P/L 160)
1 02826U 67053A  26132.30120145  .00004451  00000-0  63825-3 0  9990
2 02826  69.9150 126.2158 0001797 253.6468 106.4482 14.77195601 60937

n	dn/dt	-da/dt = -dh/dt
14.7706901	0.00015632	0.04951
14.771956	8.90E-05	0.02819

[cross-track]

с использованием значений n, dn/dt из ТЛЕ, причем тоже, видимо, усредняете эти величины по этому же интервалу времени, а затем пересчитываете эти усредненные зачения в dh/dt

Нет, не усреднял. Брал из второго TLE. Хотя логичнее наверное было взять среднее по двум TLE.
Код для расчета :

Спойлер

function dndtTodadt(dndt,n:Double):Double; //  dn/dt, n => da/dt =dh/dt
const KNA = 42241.1083269817;
begin
Result:=-(2/3)*KNA*Power(n,-5/3)*dndt;  // Power - n в степени минус 5/3
end;

[свернуть]

*  в TLE приведено значение (dn/dt) / 2
Примеры
0 CAPELLA-11 (ACADIA)
1 57693U 23126A  26121.45702955  .00005981  00000-0  78261-3 0  9993
2 57693  53.0094  27.2233 0000736  13.8846 346.2151 14.81562452145198
0 CAPELLA-11 (ACADIA)
1 57693U 23126A  26131.97901696  .00002440  00000-0  33263-3 0  9994
2 57693  53.0101 341.7351 0001340 123.2065 236.9041 14.82016278146750

n	dn/dt	-da/dt = -dh/dt
14.8156245	0.00011962	0.03769
14.8201628	4.88E-05	0.01537

Для Капеллы-11 я сделал расчет (по второму нобору ТЛЕ), и получил ваш результат.
Но есть вопрос по усреднению. Когда вычисляется dh/dt по прямым измерениям (т.е. берется отношение изменения высоты к заданному интервалу времени), то при этом учитывается изменение высоты на всем этом интервале, к примеру, изменение высоты на протяжении 10 дней, а затем пересчитывается на изменение высоты за один день, т.е. мы фактически имеем интегральную характеристику. Если использовать расчет dh/dt через  значения n, dn/dt из ТЛЕ, и не усреднять эти величины по всему интервалу, а выбирать только один набор со всего интервала, то это может привести к большим погрешностям. Я взял данные, приведенные в вашей таблице, и посчитал отношения dh/dt, найденные по разным методикам - по интегральной (изменение высоты за 10 дней), и по локальной (значениям n, dn/dt в конечной точке). Вот результат для первых 15 строчек (включая Капеллу-11):Вы не можете просматривать это вложение.
Как видим, отношения этих скоростей изменения высоты в основном колеблются около значения 1.5, а для двух объектов (включая Капеллу-11), наблюдаются всплески, доходящие до 10. Почти не сомневаюсь, что если бы использовались усредненные значения dn/dt по всему интервалу времени, то согласие было бы ощутимо лучше.

[Олег]

Почти не сомневаюсь, что если бы использовались усредненные значения dn/dt по всему интервалу времени, то согласие было бы ощутимо лучше.

Да, при небольших изменениях высоты наверное так. При больших ( например близко к ПБС ) - не знаю. Наверное скажутся нелинейности. Хотя производная высоты dh/dt линейно зависит от производной среднего движения dn/dt, но при этом увеличивается само среднее движение n,  а от n - нелинейная зависимость.

[cross-track]

Почти не сомневаюсь, что если бы использовались усредненные значения dn/dt по всему интервалу времени, то согласие было бы ощутимо лучше.

Да, при небольших изменениях высоты наверное так. При больших ( например близко к ПБС ) - не знаю. Наверное скажутся нелинейности. Хотя производная высоты dh/dt линейно зависит от производной среднего движения dn/dt, но при этом увеличивается само среднее движение n,  а от n - нелинейная зависимость.

Сделал расчет для Капеллы-11. Усредненные значения dn/dt по всему интервалу времени вычисляются очень просто - как разность крайних значений n, поделенную на длительность интервала (10.52 дня). За этот промежуток времени полуось практически не меняется; в расчетах полагалось a=7000 км. Вот что получилось:
Вы не можете просматривать это вложение.
Как видим, значение dh/dt совпало с вашим, найденным по прямым вычислениям изменения высоты на этом интервале времени.

[Старый]

Надо брать неуправляемые спутники на круговых орбитах 600-650 км и смотреть за год. Жаль все Целины-Д уже ниже. 

[Олег]

значение dh/dt совпало с вашим, найденным по прямым вычислениям изменения высоты на этом интервале времени

Да, это получилось.

Но я еще использовал dn/dt имеено из TLE.
CAPELLA-11 (ACADIA)
1 57693U 23126A  26121.45702955 .00005981  00000-0  78261-3 0  9993
2 57693  53.0094  27.2233 0000736  13.8846 346.2151 14.81562452145198
Зеленое - это половина dn/dt  . Если по ним ( и еще n) считать, то вот именно у этой Капеллы-Акадии я обнаружил проблему.
Я когда выводил формулу dh/dt = f(dn/dt,n) , я проверял по разным объектам . Сравнивал 2 способа
1) по формуле dh/dt = f(dn/dt,n) , где dn/dt - из TLE
2) по высотам dh/dt = Δh/Δt, где высоты h - расчет из n.  (  Δt - между эпохами ).
Проверил по нескольким объектам с разными высотами. Причем не два TLE через 10 дней, все что есть на неск. дней, между эпохами обычно нес. часов. Более-менее совпадало разница порядка 10%.
А вот Капеллы-Акадии - по 2 способу в неск. раз. больше.

Сейчас буду искать еще объекты с таким несовпадением.

[Старый]

Ну вот Целина-Д №68, она же Космос-2221, она же  1992-080A она же №22236
Сейчас её орбита 569х587 км. Не 630х630 а на 50 км ниже. И что? За год она потеряла 5 км. Пять всего, хотя у неё очень растопырчатая конструкция.

[cross-track]

Сейчас буду искать еще объекты с таким несовпадением.

Я в таблице с 15-ю строками привел еще один объект - STARK. У него отличие еще больше - в 10.6 раз.

[Олег]

Надо брать неуправляемые спутники на круговых орбитах 600-650 км и смотреть за год.

Ступени например. На 600 .. 650 сейчас 20 штук R/B.

[cross-track]

Но я еще использовал dn/dt имеено из TLE.

Можно усреднить все значения dn/dt по всему интервалу, но вряд ли имеет смысл это делать. гораздо проще взять разность крайних (на интервале) значений n, и поделить на этот интервал времени. Собственно, мы же так и делаем, когда прямо вычисляем по высотам dh/dt = Δh/Δt, где высоты h - по данным ТЛЕ, а Δt - между эпохами.

[Олег]

Я в таблице с 15-ю строками привел еще один объект - STARK. У него отличие еще больше - в 10.6 раз.

Да, при этом у его корешей из 185 запуска - не очень большие различия. Макс. у E - BASIL, в 3.1 раза.

[Старый]

Здравствуй, Индия!
IRS-P5, он же Картосат-1, он же 2005-017A, он же №28649.
Сейчас его орбита 594х596 км, что на 20 км выше Целины и уже ближе к 600-м.
За год он потерял 4.3 км. Тоже, чай, не чугунный шар.

[Старый]

О, вот нашлось 630, спасибо, Индия! 
Картосат-2А, 2008-021A, №32783
Правда опять нескладуха, последнюю коррекцию он сделал 16.10.25 г, поэтому данные о естественном снижении только за полгода. Так вот за полгода он потерял... ровно километр. С 632.9 до 631.9 км. То есть за год это было бы ровно два км. 

[Олег]

на круговых орбитах 600-650 км

Для ступеней с перигеем более 600 и апогеем менее 650. Также с промежутком в дне1 10..12.
В пересчете на год - 0.6 .. 4 км/год снижение.

Фиолетовый столб - разница ( отношение ) двух вариантов расчета.

[Старый]

В пересчете на год - 0.6 .. 4 км/год снижение.

Замечательно. И где тут 20 км/год? 

[Старый]

Вот ещё Уосат-12. 1999-021А. №25693
Орбита на сейчас 595х637 км
Потеря высоты за год - 2 км - с 617.9 до 616.0 км. 

[Старый]

Итого показано что потеря высоты стандартным спутником на высоте 580 км составляет 4 км/год, на высоте 630 км - 2 км/год. 
 Терять на высоте 630 км по 20 км/год это аномалия. 

[Старый]

Ну и раз пошла такая пьянка то надо проверить остальные похожие запуски с участием Поппи. И найти. 
Запуск 1965-016А.  Восемь спутников были выведены на круговые орбиты периодом 103.50 мин. Перигеи в районе 900-910 км, апогеи - 938-946 км. 
 Все объекты этого запуска естественно так там и болтаются, самый низкий на сегодня объект F - 851х878 км. Потерял 50 км высоты за 60 лет. 
 Все кроме одного. Объект G  прекратил своё существование 27.3.81 г. Сойдя с 900-км орбиты за 16 лет (5862 суток)! 
 Вот вам и ещё один надувной шарик. И надо же - в каталогах он тоже числится как Сёркал - Сёркал-4.