ТЕХНИЧЕСКИЕ СИСТЕМЫ И ТЕХНОЛОГИИ ENGINEERING AND TECHNOLOGIES
УДК 520.224.2. 224.4
МЕТОДИКА РАСЧЕТА КОЭФФИЦИЕНТОВ ОБЛУЧЕННОСТИ ЦИЛИНДРИЧЕСКОГО КОСМИЧЕСКОГО ОБЪЕКТА ПОДСВЕТКОЙ ЗЕМЛИ
А.М. Дзитоева, С.И. Ханкова
а Военно-космическая академия им. А.Ф. Можайского, Санкт-Петербург, Pocc^,[email protected] Разработана приближенная методика расчета эффективных коэффициентов облученности космического объекта цилиндрической формы, которая типична для наружных корпусов космических телескопов, а также фрагментов космических аппаратов. В настоящее время отсутствуют аналитические методики расчета интегральных и эффективных коэффициентов облученности для такой конфигурации космических объектов. Интегральные коэффициенты облученности определяют средний тепловой баланс на цилиндрической поверхности, а эффективные - на всей поверхности цилиндра с учетом его торцов. Расчеты эффективных коэффициентов облученности космических объектов подсветкой Земли необходимы для определения важнейшей компоненты энергетического баланса - удельной мощности собственного теплового излучения Земли. Такие расчеты используются для определения средних температур космических объектов. Методика основана на обобщении и аппроксимации результатов численных расчетов в диапазоне изменения высоты орбиты от 200 км до 40000 км, а также во всем возможном диапазоне изменения угла наклона оси цилиндра. Показаны неопределенности в математической модели при небольших отклонениях положения оси цилиндра от линии зенит-надир. Исследованы погрешности приближенного решения при определении интегрального коэффициента облученности боковой цилиндрической поверхности. Допустимая итоговая погрешность полученного приближенного решения оценивалась по отклонению температуры объекта в тени Земли от ее точного значения. При больших отношениях высоты цилиндра к диаметру погрешности расчета интегрального коэффициента облученности максимальны и могут приводить к погрешности определения температур в тени Земли до 1,5 К. Погрешности расчета температуры цилиндра уменьшаются с уменьшением отношения его высоты к диаметру Они минимальны для модели тонкого диска. Погрешности расчета температур значительно снижаются для случая нахождения цилиндрического объекта на подсвеченном Солнцем участке траектории, а еще более - при расчетах нестационарного теплового режима. Предложенная методика расчета достаточно точна, проста и удобна для расчетов.
Ключевые слова: космический объект цилиндрической формы, коэффициент облученности космического объекта, удельный тепловой поток излучения Земли, тепловой баланс космического объекта.
CALCULATION METHODS FOR IRRADIANCE COEFFICIENTS OF CYLINDRICAL SPACE OBJECT BY THE EARTH RADIATION
A. Dzitoevb, S. Khankovb
b Military Space Academy n.a. A.F. Mozhaisky, Saint Petersburg, Russia, [email protected]
The paper deals with approximate calculation methods for effective irradiance coefficients of cylindrical form space object typical for external cases of space telescopes and fragments of spacecrafts as well. There are no analytical calculation methods nowadays for integrated and effective irradiance coefficients for space objects of such configuration. Integrated irradi-ance coefficients define the average thermal balance on a cylindrical surface and effective ones - all over the cylinder surface taking into account its butt ends. Calculations of effective coefficients for space objects irradiance by the Earth radiation are necessary for definition of the major component of power balance - the specific power of the Earth inherent thermal radiation. Such calculations are used for determination of average temperatures of space objects. The technique is based on generalization and approximation of results for numerical calculations in the range of orbit height change from 200 km to 40000 km and also in all possible range of a tilt angle change of a cylinder axis. Uncertainties in mathematical model at small location deviations for cylinder axis from zenith-nadir line are shown. Errors of the approximate solution are analyzed during determination of integrated irradiance coefficient for a lateral cylindrical surface. The admissible total error of the received approximate solution was estimated at a temperature deviation of an object in the Earth shadow from its exact value. At the large ratios of the cylinder height to its diameter calculation, errors of integrated irradiance coefficient are maximum and can cause temperature error definitions in the Earth shadow up to 1,5 K. Errors of cylinder temperature calculation decrease with reduction of its height-to-diameter ratio. They have minimum values for a thin disk model. Temperatures calculation errors are considerably decreased for a case of cylindrical object location on the trajectory section illuminated by the Sun, and even more - at calculations of a non-stationary thermal mode. The offered calculation technique is rather accurate, simple and convenient for calculations.
Keywords: cylindrical form space object, irradiance coefficient of space object, specific thermal radiation flux of the Earth, space object thermal balance.
Введение
При расчетах нестационарных тепловых режимов космических объектов (КО) необходимо определение величин падающих и поглощаемых поверхностью КО тепловых потоков от прямого и отраженного Землей солнечного излучения, а также собственного теплового излучения Земли [1-4]. Первые две ком-
4
поненты теплового баланса КО вычисляются следующим образом. Прямая солнечная засветка поверхности КО определяется через солнечную постоянную с учетом площади его миделя и коэффициента поглощения солнечного излучения [1]. Эта компонента теплового баланса вычисляется без особых затруднений. Вторая компонента вычисляется через так называемый комбинированный коэффициент облученности, расчет которого связан с большими затруднениями [1] и, как правило, осуществляется численными методами. Третья компонента теплового баланса определяется через величину удельного теплового потока в полном спектре собственного теплового излучения Земли, а также через коэффициент облученности, который представляет собой отношение падающего на поверхность КО удельного теплового потока к тепловому потоку, излучаемому Землей с верхней границы атмосферы. В этом потоке учитывается тепловое излучение атмосферы и излучение подстилающей земной поверхности, проходящее через окна прозрачности атмосферы в спектре инфракрасного теплового излучения [1].
Известно, что спектр собственного теплового излучения Земли с учетом атмосферы имеет ярко выраженный селективный характер. Однако на практике принято пользоваться для расчетов теплового баланса интегральными удельными тепловыми потоками [1]. По этой причине значительные средства затрачиваются на экспериментальные определения глобальной (усредненной по всей поверхности Земли) величины излучаемой Землей в космическое пространство удельной мощности [5, 6]. Такие экспериментальные исследования, проводимые с космических аппаратов (КА), позволили уточнить эту величину с 235 Вт/м2 [5] до 239 Вт/м2 [6].
Для расчетов тепловых потоков, поглощаемых КО, необходимо определение интегрального коэффициента облученности (ИКО) [7-12] КО земной подсветкой в зависимости от его конфигурации. ИКО определяет удельный тепловой поток, который усредняется по всей поверхности КО и вычисляется интегрированием локальных коэффициентов облученности для каждого элемента поверхности. Для сферического объекта коэффициент облученности описывается просто [1, 12], что облегчает расчет нестационарных температур сферического КО на любой орбите [12-15]. Большой практический интерес представляет решение той же задачи для КО цилиндрической формы, которую имеют многие реальные объекты, в том числе корпуса телескопов, базируемых на борту КА [16]. Сами КА также содержат фрагменты цилиндрической формы. Это определяет актуальность данной задачи.
Для объектов с конфигурацией, составленной из разного рода поверхностей (параллелепипед, цилиндр, конус) используется эффективный коэффициент облученности (ЭКО), который описывает удельный тепловой поток, поглощенный всей поверхностью КО. Частью расчета ЭКО является определение ИКО цилиндрической боковой поверхности. Для цилиндрического КО (ЦКО) в эффективный коэффициент облученности входят также коэффициенты облученности Землей торцов цилиндра [14, 15].
Целью настоящей работы являлось получение приближенного аналитического описания зависимости эффективного коэффициента облученности ЦКО от высоты орбиты, угла наклона его оси относительно плоскости местного горизонта и отношения высоты цилиндра к его диаметру.
Разрабатываемая методика должна основываться на аналитических формулах, составляющих алгоритм расчетов, исключающий необходимость проведения интегрирования громоздких выражений для локального коэффициента облученности.
Физическая и математическая модели процесса подсветки цилиндра собственным
тепловым излучением Земли
Стационарная температура Т КО произвольной выпуклой формы в тени Земли может быть описана формулой [10-14]
т=т ёл; т=^, (1)
где Тг = 255 К - эффективная температура Земли [15, 16]; фе - ЭКО объекта; Q = 239 Вт/м2 - удельная мощность собственного теплового излучения Земли [15, 16]; ст = 5,67-10-8 Вт/(м2-К4) - постоянная Стефа-на-Больцмана.
Для ЦКО ЭКО определяется на основании данных [12, 13] в виде
т Фь + 2п(Рш . п L (2)
фс = , , 0-. п = ~7, (2)
1 + 2п d
где фь - коэффициент облученности торцевой поверхности цилиндра (в расчетных формулах учитывается наличие двух торцов); фш - интегральный коэффициент облученности боковой поверхности цилиндра; Ь - высота цилиндра; й - диаметр его основания; фс - эффективный коэффициент облученности.
На рис. 1 представлены варианты ориентации оси цилиндра относительно линии зенит-надир; а -угол наклона оси ЦКО относительно этой вертикальной линии.
При вертикальной ориентации цилиндра, т.е. при а = 0 (рис. 1, а), верхний торец цилиндра не подсвечен Землей, коэффициент облученности нижнего торца вычисляется как для площадки, лежащей в плоскости местного горизонта [1, 10]
/1 \-2 И Фо =(1+а) ; а = —, к
где И - высота площадки над поверхностью Земли; К = 6371 км - радиус Земли.
(3)
We
а б в
Рис. 1. Варианты ориентаций оси цилиндра относительно линии зенит-надир и плоскости местного горизонта (ПМГ): а - вертикальная ориентация оси цилиндра; б - ориентация оси цилиндра под углом а,
равным углу 0о; в - ориентация оси цилиндра под углом а, большим угла 0о и меньшим 90°. Углы соответствуют: фс - угловая координата по цилиндрической поверхности; а - угол наклона оси цилиндра, отсчитываемой от оси 0о - половинный угол, под которым из любой точки цилиндра видна Земля. Линия АВ на образующей цилиндрической поверхности - наиболее удаленная от Земли
Для боковой поверхности цилиндра в этом случае расчет проводится как для вертикально ориентированной площадки по формуле [1, 10]
1 (arcsin^-Тф0'V1 -Фо ) >
ТГ * '
(4)
где фъс - коэффициент облученности вертикальной площадки, одновременно локальный и интегральный; индексы Ь и с относятся к боковой поверхности цилиндра.
Для горизонтально ориентированного цилиндра также имеются расчетные формулы для определения ЭКО [15, 16]. Для промежуточных ориентаций оси 0 < а < п/2 (рис. 1, б, в) расчеты существенно усложняются. При наклоне оси цилиндра относительно оси ZN на угол а (рис. 1) угол наклона у элемента поверхности цилиндра, задаваемого угловой координатой на цилиндрической поверхности ус, отсчитываемой от образующей, наиболее близкой к Земле, может быть определен с использованием данных [1, 2] с помощью соотношения
у = arccos х ; х = sin а • cos у с
(5)
Важным параметром, определяющим характеристики облученности элемента поверхности КО любой формы, является половинный угол 90, под которым с этого элемента поверхности видна Земля, описываемый формулой [1, 12]
0о = агс8т(л/Ф0).
Зависимость локального коэффициента облученности плоского элемента поверхности собственным тепловым излучением Земли от угла наклона касательной к этому элементу, а также от высоты этого элемента над поверхностью Земли описывается известными соотношениями [1, 12], которые в наиболее компактной форме могут быть выражены через х в формуле (5). При этом в зависимости от величины угла наклона у коэффициент облученности описывается одной из трех формул:
_ п _ п _
Ф1 = ф0 • х при 0 < у < — - 00 ; ф3 = 0 при — + 00 < у < п ;
ф2 =
Ф1
1 2 . 1 + — arcsin п
-1 •
1
+1 arcsin п
Ф0 - х 1 - х2
(1-Ф0)'(Ф0 -х2)
(6) (7)
пп Формула (7) используется при — -00 < у < — + 00.
х
2
2
х
0
ИКО цилиндра определяется путем интегрирования (6) и (7) в указанных при этих формулах пределах. Затем полученные интегралы суммируются и делятся на л.
Для расчетов ЭКО ЦКО по формуле (2) необходимо определить фь и tyint.
При определении коэффициентов облученности торцов ЦКО введены обозначения: верхний торец, наиболее удаленный от поверхности Земли, характеризуется коэффициентом облученности фЬи, нижний торец, ближний к поверхности Земли, - коэффициентом облученности фы. Тогда в (2) вместо фь подставляется
ФА = 0,5 Ы +Фы ) . (8)
Для расчета коэффициентов облученности торцов ЦКО в формулы (6) и (7) подставляется значение х = cos у. Зависимость угловой координаты у от угла наклона оси ЦКО относительно плоскости местного горизонта для нижнего и верхнего торцов определяются соотношениями у = а ; уи = л - а . Использование формулы (8) обеспечивает предельные переходы и непрерывность ЭКО во всем диапазоне изменения угла наклона оси цилиндра а.
Методика расчета интегрального коэффициента облученности боковой поверхности цилиндра
При расчете ИКО ЦКО в пределах изменения угла наклона оси цилиндра 0 < а < 90 возникает принципиальная неопределенность в математической модели. При а = 0 (рис. 1, а) расчет фш в формуле (2) заменяется расчетом коэффициента облученности по формуле (4) - вариант модели 1 (ВМ 1). Эта модель предполагает использование локальных коэффициентов облученности ф! и ф2 из (7) при малых величинах х. Однако границы допустимых значений а и соответственно х неопределенны. При значительных величинах угла наклона оси, вплоть до а = 90, ВМ 1 неприменима. В то же время до а < 90 неправомерно использование уравнений (6) и (7) для определения ИКО на основе интегрирования ф! и ф2. Для а > 90 (рис. 1, в) ИКО определяется интегрированием по участкам координат ф! и ф2 - вариант модели 2 (ВМ 2).
Проведенные расчетные исследования подтвердили указанную неопределенность в выборе математической модели на участке 0 < а < 90. При использовании ВМ 1 и ВМ 2 вычисленные значения коэффициентов облученности боковой поверхности ЦКО на одних и тех же диапазонах изменения а недостаточно точно стыкуются друг с другом. Это обосновывает необходимость получения хотя бы приближенного единого описания ИКО во всем диапазоне изменения угла наклона 0 < а < л на основе получения аппроксимационных функций с определением погрешности получаемого конечного решения.
Зависимость коэффициента облученности фа боковой поверхности цилиндра от высоты и угла наклона оси аппроксимировалась формулами
Фа = ФЬ t1 - f (а)]+Фт • f (а) ; ФА =ф(а = 0) ; Фт =ф(а=л/2) , (9)
где фь и фт - функции только высоты; f (а) - функция только угла наклона оси ЦКО.
Зависимость f (а) была аппроксимирована функцией вида
f (а) = а1 • sin ( bl -а + c1 ) + а2 • sin (b2 - а- c2 ), (10)
которая изменяется от нуля при а = 0 до единицы при а = л/2.
Зависимости фЛН) аппроксимирована формулой
h
Фт = a3eXP(-b3H) + a4exP(-b4H); H = 1000 . (11)
Величина H соответствует высоте траектории, исчисляемой в тысячах километров.
Для обеспечения минимальных погрешностей аппроксимационных зависимостей (9)-(11) относительно исходных расчетных зависимостей ИКО, определяемых по формулам (4), (6) и (7), коэффициенты в аппроксимационных формулах определялись для пяти интервалов высот H. Значения коэффициентов аппроксимации представлены в таблице.
Погрешность аппроксимации в максимуме зависимости, при а = л/2, не превышает 1%. В то же время максимальная погрешность, реализуемая в области 20° < а < 60°, достигает 2,3%, 4%, 4,5%, -5,7%, - 4,7% при H равном 0,2-2, 2-5, 5-10, 10-20, 20-40 соответственно. Для вытянутого цилиндрического стержня (n > 10) из (2) следует: фс ~ фпг, поэтому для расчетов температур по формуле (1) достаточно использования ИКО. Погрешности определения стационарных температур ЦКО в тени Земли составляют для указанных диапазонов 1,2 К, 1,5 К, 1,4 К, -1,2 К и -0,8 К, что вполне удовлетворительно для практических расчетов. На подсвеченном Солнцем участке траектории погрешности расчетов температур будут еще меньше, а при расчетах нестационарных температур при движении ЦКО по траектории с попеременной сменой солнечных и теневых участков погрешности расчетов температур будут дополнительно сглаживаться за счет термической инерции ЦКО при его достаточно большой полной теплоемкости.
На рис. 2 представлены типовые зависимости фш и фс от угла наклона оси цилиндра а для разных высот ЦКО над поверхностью Земли в наиболее важном диапазоне изменения h. На рис. 2 представлены
вычисленные по разработанным формулам зависимости от угла наклона оси цилиндра для высот от 200 км до 2000 км. Сплошными линиями указаны зависимости ЭКО для отношения высоты цилиндра к диаметру его основания п = 3, штриховыми - ИКО боковой цилиндрической поверхности. Там же приведены для сравнения значения ИКО сферического КО.
№ Коэффициенты Диапазон высот H в тысячах километров
п/п аппроксимации 0,2-2 2-5 5-10 10-20 20-40
1 ai 0,7753 0,4995 0,4944 0,4063 0,3164
2 bi 0,01652 0,0001099 0,0003042 0,0001966 0,0009266
3 Ci 0,08263 1,591 1,863 1,491 1,879
4 0,2248 0,5114 0,5458 0,6072 0,7152
5 b2 0,04872 0,03407 0,03065 0,0269 0,02303
6 C2 2,822 1,502 1,201 0,8528 0,5199
7 аз 0,09587 0,1882 0,2149 0,1875 0,1162
8 Ьз 0,002238 0,0006294 0,1171 0,0002392 0,0001389
9 a4 0,3555 0,2039 0,1171 0,05989 0,02373
10 b4 0,0002796 0,0001494 0,0001 0,0001 0
Таблица. Значения коэффициентов аппроксимации в формулах (10) и (11)
Рис. 2. Зависимость интегрального (штрихпунктирные линии) и эффективного (сплошные линии) коэффициентов облученности от угла наклона оси цилиндра при высотах h, км, равных: 200 (1), 600 (2), 1000 (3), 1500 (4), 2000 (5). Горизонтальными сплошными линиями указаны интегральные коэффициенты
облученности сферического объекта
Как и ожидалось, погрешности расчета ЭКО меньше, чем вычисленных значений ИКО. Для вышеуказанных диапазонов изменения высоты ЦКО максимальные погрешности составляют 2%, -2,7%, 3,8%, -3,9% и 3,8%. Соответственно для цилиндра с умеренным отношением п погрешности определения среднеповерхностной температуры цилиндра будут меньше, чем в случае вытянутого цилиндра с большим п.
Обсуждение результатов
Предложенная методика расчета основывается на следующем алгоритме определения ЭКО для ЦКО. На первом этапе определяются значения функций f (а) и фт(Н) для заданных угла наклона оси цилиндра и высоты ЦКО над поверхностью Земли по формулам (10) и (11). При этом используется значения коэффициентов в формулах из таблицы. Затем определяется значение ИКО боковой поверхности цилиндра фа по формуле (9), в которую подставляются также формулы для щ = фЬс из (4) и фт = ф0 из (3). Затем полученное значение ИКО фа подставляется в (2) вместо значения ф/п/. В ту же формулу (2) подставляется значение фь из (8), а также конкретная величина п. В результате получается значение ЭКО, которое для определения температуры ЦКО в тени Земли подставляется в формулу (1). При любой высоте ЦКО над поверхностью Земли погрешность расчета ИКО по формулам (9)—(11) при горизонтальном положении оси цилиндра, соответствующем максимальному значению ИКО, не превышает 1%. Максимальные погрешности реализуются в диапазоне изменения углов 20°< а < 60°, в котором осуществлялось «сшивание» зависимостей, вычисленных по двум моделям - ВМ 1 и ВМ 2. В указанной области изменения углов погрешности определения стационарных температур по формуле (1) не превышают 1,5 К для случая длин-
ного стержня (при n > 10). С уменьшением n погрешности определения ЭКО по формуле (2) снижаются, а при n < 0,1 (модель тонкого диска) расчеты становятся точными. Погрешности значительно снижаются также на солнечной стороне траектории и при расчетах нестационарного теплового режима. При а, равном 60° и 120°, значения ИКО и ЭКО цилиндра при n = 3, а также ИКО сферы примерно равны, причем ИКО сферы занимает промежуточное положение между значениями ИКО цилиндра при а = 0 и а = л/2.
Заключение
Предложенная приближенная аналитическая методика расчета эффективного коэффициента облученности для космического объекта цилиндрической формы реализует единственный способ преодоления неопределенностей в физической и математической моделях в области изменения углов наклона оси цилиндра до а < 90. Методика обеспечивает высокую для практических расчетов точность. Она проста в использовании и существенно облегчает расчеты нестационарных температур цилиндрического космического объекта. Альтернативой может быть численный расчет, обеспечивающий определение распределений поглощаемых лучистых потоков, а также температур по поверхности цилиндрического космического объекта. Такой расчет оправдан для определения температурных полей в объеме цилиндрических оболо-чечных конструкций с учетом прямой солнечной подсветки и отраженного Землей солнечного излучения. Однако на начальной стадии теплового проектирования телескопов космического базирования аналитические методы расчета незаменимы, поскольку позволяют быстро и надежно определять требования к основным параметрам конструкции на основе расчетов интегрального теплового баланса цилиндрического космического объекта.
Литература
1. Моделирование тепловых режимов космического аппарата и окружающей его среды / Под ред. Г.П. Петрова. М.: Машиностроение, 1971. 382 с.
2. Дубошин Г.Н. Небесная механика. Основные задачи и методы. М.: Наука, 1968. 800 с.
3. Чеботарев Г. А. Аналитические и численные методы небесной механики. М.-Л.: Наука, 1965. 367 с.
4. Смарт УМ. Небесная механика. М.: Мир, 1965. 502 с.
5. Keihl J.T., Trenberth K.E. Earth's Annual Global Mean Energy Budget // Bull. Am. Meteorol. Soc. 1997. V. 78. N 2. P. 197-208.
6. Trenberth K.E., Fasullo J.T., Keihl J. Earth's global energy budget // Bull. Am. Meteorol. Soc. 2009. V. 90. N 3. P. 311-323.
7. Дульнев Г.Н., Семяшкин Э.М. Теплообмен в радиоэлектронных аппаратах. Л.: Энергия, 1968. 360 с.
8. Исаченко В.П., Осипова В.А., Сукомел А.С. Теплопередача. М.: Энергия, 1975. 488 с.
9. Уонг X. Основные формулы и данные по теплообмену для инженеров: Справочник: Пер. с англ. М.: Атомиздат, 1979. 216 с.
10. Зигель Р., Хауэлл Дж. Теплообмен излучением: Пер. с англ. М.: Мир, 1975. 934 с.
11. Абдусаматов Х.И., Лаповок Е.В., Ханков С.И. Методы обеспечения термостабильности космического телескопа - солнечного лимбографа. СПб: Изд-во Политехн. ун-та, 2008. 195 с.
12. Каменев А.А., Лаповок Е.В., Ханков С.И. Аналитические методы расчета тепловых режимов и характеристик собственного теплового излучения объектов в околоземном космическом пространстве. СПб: НТЦ им. Л.Т. Тучкова, 2006. 186 с.
13. Баёва Ю.В., Лаповок Е.В., Ханков С.И. Методика расчета нестационарных температур космического объекта, движущегося по эллиптической орбите // Научно-технический вестник информационных технологий, механики и оптики. 2013. № 6 (88). С. 67-72.
14. Баёва Ю.В., Лаповок Е.В., Ханков С.И. Аналитическая методика расчета тепловых потоков в околоземном космическом пространстве, формирующих тепловой режим космических телескопов // Оптический журнал. 2013. Т. 80. № 5. С. 30-37.
15. Баёва Ю.В., Лаповок Е.В., Ханков С.И. Метод поддержания заданного температурного диапазона космического аппарата, движущегося по круговой орбите с заходом в тень Земли // Изв. вузов. Приборостроение. 2013. Т. 56. № 7. С. 56-61.
16. Баёва Ю.В., Ханков С.И. Обеспечение термостабильности телескопа дистанционного зондирования Земли за счет выбора параметров бленды // Изв. вузов. Приборостроение. 2013. Т. 56. № 5. С. 76-79.
Дзитоев Азамат Миронович - адъюнкт, Военно-космическая академия им. А.Ф. Можайского, Россия, Санкт-
Петербург, [email protected] Ханков Сергей Иванович - доктор технических наук, ведущий научный сотрудник, Военно-космическая
академия им. А.Ф. Можайского, Санкт-Петербург, Россия, [email protected]
Azamat Dzitoev Postgraduate in a military academy, Military Space Academy n.a. A.F. Mozhaisky,
Saint Petersburg, Russia, [email protected] Sergei Khankov D.Sc., leading scientific researcher, Military Space Academy n.a. A.F. Mozhaisky,
Saint Petersburg, Russia, [email protected]