Тепловое подобие космических объектов типовых конфигураций

6

ТЕХНИЧЕСКИЕ СИСТЕМЫ И ТЕХНОЛОГИИ ENGINEERING AND TECHNOLOGIES

УДК 520.224.2. 224.4

ТЕПЛОВОЕ ПОДОБИЕ КОСМИЧЕСКИХ ОБЪЕКТОВ ТИПОВЫХ

КОНФИГУРАЦИЙ А.М. Дзитоева, С.И. Ханкова

а Военно-космическая академия имени А.Ф. Можайского, Санкт-Петербург, Россия, dzitoi8@gmail.com

Тепловое подобие объектов различной конфигурации определяется равенством их стационарных среднеповерхност-ных температур в тени Земли, что эквивалентно равенству их эффективных коэффициентов облученности собственным тепловым излучением Земли. Среди типовых конфигураций выбраны конус, цилиндр и сфера. В отличие от двух последних фигур, расчет коэффициента облученности конического объекта является наиболее сложным и содержит ряд неопределенностей. Изложена методика расчета интегральных и эффективных коэффициентов облученности космического объекта конической формы, которая типична для фрагментов космических аппаратов. Интегральные коэффициенты облученности определяют средний тепловой баланс на боковой поверхности цилиндра и конуса, а также полный энергетический баланс на поверхности сферы. Эффективные коэффициенты облученности определяют полный падающий удельный поток излучения Земли на всю поверхность цилиндрического или конического объекта с учетом их оснований. По данным об эффективных коэффициентах облученности определяются средние стационарные температуры космических объектов в тени Земли, а также и на подсвеченном Солнцем участке траектории с учетом двух дополнительных компонент энергетического баланса - прямого и отраженного Землей солнечного излучения. Исследования проводились в диапазоне изменения высоты орбиты от 200 до 40000 км в зависимости от угла наклона оси цилиндра и конуса относительно линии зенит-надир. Условия подобия цилиндра и конуса определены при равных величинах отношения высоты фигуры к диаметру основания.

Ключевые слова: космический объект конической, цилиндрической и сферической формы, эффективный коэффициент облученности космического объекта, удельный тепловой поток излучения Земли, тепловой баланс космического объекта, тепловое подобие космических объектов.

3 Mozhaisky Military Space Academy, Saint Petersburg, Russia, dzitoi8@gmail.com

Thermal similarity of objects of various configuration is defined by equality of their stationary surface average temperatures in the Earth shadow that is equivalent to equality of their effective irradiance coefficients by own thermal radiation of the Earth. Cone, cylinder and sphere are chosen among standard configurations. Unlike two last figures, calculation of irradiance coefficient for conic object is the most difficult and contains a number of uncertainties. The method of calculation for integrated and effective irradiance coefficients of space object with a conic form is stated which is typical for fragments of spacecrafts. Integrated irradiance coefficients define the average thermal balance on a lateral surface of the cylinder and cone, and also full power balance on a sphere surface. Effective irradiance coefficients define a full falling specific stream of the Earth's radiation on the whole surface of cylindrical or conic object taking into account their bases. By data about effective irradiance coefficients, the average stationary temperatures of space objects in the Earth shadow are defined, as well as on the trajectory part illuminated by the Sun taking into account two additional components of power balance - direct sunlight and reflected by the Earth. Researches were conducted in the height change range for an orbit from 200 to 40000 km depending on a tilt angle of the cylinder and cone axis relative to zenith-nadir line. Similarity conditions for the cylinder and cone are defined at equal ratio sizes of the figure height to base diameter.

Keywords: space object of the conic, cylindrical and spherical form, effective irradiance coefficient of space object, specific thermal stream of the Earth's radiation, thermal balance of space object, thermal similarity of space objects.

Для расчетов теплового режима космического объекта (КО) при его движении по различным орбитам [1-4] необходимо определение поглощаемых поверхностью КО лучистых тепловых потоков, в том числе обусловленных собственным тепловым излучением Земли. Это требует определения интегрального коэффициента облученности (ИКО) [5-10] отдельных участков поверхностей КО земной подсветкой. Для сферического объекта коэффициент облученности описывается через ИКО достаточно просто [1, 6], что облегчает расчет средних температур сферического КО на любой орбите [6-10]. Ранее авторами было получено решение той же задачи для КО цилиндрической формы, причем для такого КО, имеющего три поверхности (боковую и две торцевых), используется эффективный коэффициент облученности (ЭКО) [5]. ИКО и ЭКО определяются через описание локального коэффициента облученности (ЛКО) для элементов боковой поверхности криволинейной формы или для всей поверхности плоских элементов. Актуальность данной задачи определяется тем, что космические аппараты содержат фрагменты различных конфигураций, в том числе конической формы.

THERMAL SIMILARITY OF SPACE OBJECTS OF STANDARD CONFIGURATIONS A.M. Dzitoeva, S.I. Khankov3

.а

Введение

Частью расчета ЭКО является определение ИКО характерной боковой поверхности, определяющей тип конфигурации КО (цилиндрическая, коническая). Для цилиндрического и конического КО в эффективный коэффициент облученности входят также ЛКО торцов цилиндра или основания конуса [5, 6, 9]. Тепловое подобие определяется равенством температур КО различных конфигураций в тени Земли, что эквивалентно равенству их ЭКО [6], поскольку интегральный удельный тепловой поток, падающий на всю подсвеченную Землей часть поверхности КО, равен произведению ЭКО на удельную мощность теплового излучения Земли, которая на верхней границе атмосферы равна 239 Вт/м2 [11, 12].

Целью настоящей работы являлось исследование зависимости ИКО и ЭКО конического КО от высоты орбиты, угла наклона его оси относительно плоскости местного горизонта и половинного угла при вершине конуса, а также сопоставление ЭКО конуса, цилиндра и сферы с целью определения областей теплового подобия указанных фигур.

Модель подсветки конуса собственным тепловым излучением Земли

Для расчетов ИКО и ЭКО прямого кругового конуса тепловым излучением Земли необходимо задать все угловые координаты, определяющие ЛКО каждого элемента его поверхности. На рис. 1 представлены схемы различных вариантов ориентации конуса относительно линии зенит-надир ^Ы).

2 2 • щ г

Wt

ПМГ

Wb

Г

N

*

N б

N

Рис. 1. Варианты ориентаций оси конуса относительно линии зенит-надир () и плоскости местного горизонта (ПМГ): а - вертикальная ориентация оси конуса при направлении вершины конуса в зенит;

б - ориентация оси конуса под таким минимальным углом а, начиная с которого часть участка поверхности конуса не подсвечивается Землей; в - горизонтальная ориентация оси конуса.

Фа - угловая координата по конической поверхности; а - угол наклона оси конуса, отсчитываемой от направления в зенит; р - половинный угол при вершине конуса; фь - угол наклона торца конуса к ПМГ;

9о - половинный угол, под которым из любой точки конуса видна Земля. Линия АВ на образующей конической поверхности - наиболее удаленная от Земли

Необходимо отметить, что угловая координата отсчитывается от образующей конической поверхности, наиболее близкой к Земле. Исключение составляет вертикальное положение конуса, когда он обращен вершиной либо в зенит, либо в надир.

На рис. 1, а, представлен конический КО, вершина которого обращена в зенит, а основание - в надир и совпадает с плоскостью местного горизонта, при этом угол наклона оси конуса а = 0. Для описания модели конуса в данной ситуации достаточно двух углов, р и при этом ЛКО для любой образующей конуса не зависит от На рис. 1, б, приведены все угловые координаты, необходимые для расчетов ИКО боковой конической поверхности и ЛКО торцевой поверхности. На этом рисунке представлена ситуация, когда угол наклона оси конуса а достиг такой величины, что наиболее удаленная от Земли образующая конической поверхности АВ уже выходит из зоны подсветки Землей. При этом выполняется условие 60 = а+р. ЛКО основания конуса описывается через угол равный углу наклона диаметральной линии основания АС к ПМГ. На рис. 1, в, приведено положение конуса, при котором его ось лежит в ПМГ, а угол наклона оси по отношению к линии а = 90°. В этом положении основание конуса лежит в плоскости, перпендикулярной ПМГ и совпадающей с линией 2Ы, при этом = 90°.

Для вычислений ИКО и ЭКО конуса достаточно проводить расчеты в диапазоне изменения угла наклона оси 0 < а < 180°.

Локальные коэффициенты облученности элементов конической и всей торцевой поверхности конуса

Формулы для расчета ЛКО для каждого элемента конической поверхности подсветкой Земли в общем случае можно получить на основе преобразования известных соотношений [1, 6] для элементарной площадки.

а

в

В общем случае локальный коэффициент облученности ф в зависимости от угла наклона площадки у может быть описан одной из трех формул:

Ф = Ф0 х при 0 < у < л/2 - 60; ф = 0 при л/2 + 60 < у < л,

а в промежуточном диапазоне изменения у - при л/2 - 60 < у <л/2 + 60 используется формула

Ф0 х

ф= 2

2

l +—arcsin Y

+—arcsin,

n Vl - x2 n

-Ф0)-Z ,

л

где Y и Z - переменные, описываемые соотношениями

Y = — -1 ; |y| < 1; Z = ф0 - x2; Z > 0.

\Фо V1 - x2

В общем случае при произвольной ориентации оси конуса относительно плоскости местного горизонта можно записать выражения для х, а также для угла наклона произвольной элементарной площадки на конической поверхности к плоскости местного горизонта у, для коэффициента облученности ф0 плоской площадки, лежащей в плоскости местного горизонта, и половинного угла 0о, под которым из любой точки КО видна Земля, в виде

x = cosp- sina- cosy k - sinp- cosa; y = arccos x; ф0 =|—— | ; 80 = arcsin,

^ R + h )

где p - половинный угол при вершине конуса; a - угол наклона оси конуса к линии зенит-надир; ук - угловая координата в полярной системе координат, отсчитываемая по поверхности конуса, от ближней к Земле линии; R = 6371 км - радиус Земли; h - высота КО над поверхностью Земли.

В случае реализации условий |y| > 1 или Z < 0 необходимо принимать в расчетах по формуле

Y = +1 и Z = 0.

Интегральные коэффициенты облученности боковой поверхности и основания конуса

для характерных ориентаций

В ряде частных случаев выражение (4) для x упрощается. 1. При a = 0 и a = л координаты х и у не зависят от ук, а зависят только от p, и ИКО боковой поверхности конуса совпадает с ЛКО для каждой его образующей.

1. В случае a = 0 (вершина конуса направлена в зенит) x = x0 = -sinp при 0 < p < 00, расчет проводится по формуле (2).

2. В случае a = л (вершина конуса направлена в надир) x = xл = sinp при 0 < p < 00, расчет проводится по формуле (2), а при 60 < p < л/2 расчет проводится по формуле, вытекающей из (1):

ф = Ф0= .

1. При a = л/2 (ось конуса лежит в плоскости местного горизонта) x = xg = cospcosyk и ИКО конуса определяется интегрированием выражений для ЛКО.

1. В случае относительно малых значений половинного угла при вершине конуса p, задаваемого критерием малости 0 < p < л/2 - 00, вычисления проводятся интегрированием в трех диапазонах изменения у. Указанные три диапазона имеют свои пределы интегрирования, задаваемые областью значений у:

- для диапазона 0 < у < л/2 - 00 расчет проводится интегрированием формулы, вытекающей из (1),

Ф = Ф0 xg =Ф0^ p cos Ук;

- в пределах диапазона л/2 - 00 < у < л/2 + 00 расчет проводится интегрированием формулы (2);

- для диапазона л/2 + 00 < у < л ф = 0.

Два интеграла, определенные по формулам (2) и (5), суммируются и делятся на л.

2. В случае л/2 - 00 < p < л/2 расчет проводится интегрированием формулы (2) по координате у в пределах ее изменения 0 < у < 2л, полученное значение делится на 2л.

В случае вертикальной ориентации оси конуса (a = 0 и a = л) и при p = 0 конус вырождается в цилиндр, и значения ИКО боковой поверхности фя в точности равны коэффициенту облученности боковой поверхности цилиндра (его ЛКО), ось которого совпадает с направлением зенит-надир. В этом случае коэффициент облученности описывается соотношением [6]

Ф,

= i

77 \ /

Если вершина конуса направлена в надир и р = л/2, значения фя точно совпадают с величинами ко -эффициентов облученности плоской площадки, лежащей в плоскости местного горизонта фо, определяемыми формулой (4). При ориентации вершины конуса в зенит ЛКО его торца описывается формулой для Фо из (4).

При горизонтальной ориентации оси конуса (а = л/2) формулой (6) описывается ЛКО его торца.

Методики расчета ИКО и ЭКО в зависимости от ориентации оси конуса

В расчетах ИКО использовался подход, описанный для цилиндра в [5], с тем отличием, что пределы интегрирования выбирались с учетом угла наклона ближайшей к Земле образующей конуса, отличающегося от угла наклона оси конуса на угол р. Это приводит к необходимости при расчетах зависимостей ф1ИХа) ввести рассматриваемые диапазоны и пределы интегрирования с учетом поправки на половинный угол при вершине конуса. При расчете ИКО конической поверхности в пределах изменения угла наклона оси конуса 0 < а < 90 + р проводится интегрирование формулы (2) в пределах от нуля до 2л, и интеграл делится на 2л. При 90 + р < а < л - 90 - р проводится интегрирование локальных коэффициентов облученности ф! и ф2 из (1) и (2) в указанных пределах, и сумма интегралов делится на л. Для третьего участка изменения угла наклона оси конуса л - 90 - р < а < л осуществляется интегрирование формулы (2) в пределах от нуля до 2л, и полученный интеграл делится на 2л.

Для конуса ЭКО в случае равенства степеней черноты боковой поверхности и основания может быть определен из соотношения [6]

ф4 81Пр+ фш

Ф* =-

sinp +1

Каждый из коэффициентов облученности фЬ и <int зависит от высоты конического объекта над поверхностью Земли h и от угла наклона оси конуса относительно линии зенит-надир а, а от половинного угла при вершине конуса р зависит только <int.

Результаты расчетов ИКО и ЭКО конуса

На основе описанной методики расчетов <int, а также фЬ, определяемого из формул (1) или (2) как для локального коэффициента облученности (выбор формулы определяется с учетом величины угла наклона торца к плоскости местного горизонта), были получены зависимости ИКО от угла наклона оси конуса а для различных значений р и высоты конуса над поверхностью Земли h.

Исследования зависимостей ИКО <int от а проводились для диапазона изменения параметров 200 < h < 40000 км, 0° < р < 45°. На рис. 2 представлены зависимости ИКО боковой поверхности конуса от угла наклона его оси для разных значений h и р. При значениях угла наклона а = 0° и а = 180° реализуются предельные переходы к описанным выше частным случаям.

20 40 60 80 100 120 140 160 а°

Рис. 2. Зависимости ИКО конической поверхности от угла наклона оси конуса при р = 15° (сплошные линии), р = 30° (штриховые линии) и р = 45° (пунктирные линии).

Группы линий сверху вниз соответствуют высоте КО над поверхностью Земли 200, 2000 и 10000 км

Из рис. 2 можно сделать вывод о том, что с ростом угла р разрыв величины <int между точками а = 180° и а = 0° увеличивается. В области изменения величины а от 95° до 110° линии на каждой высоте пересекаются, т.е. величина ИКО не зависит от угла конусности р. Если при относительно малых углах р при больших высотах зависимости <int (а) имеют максимумы вблизи точек пересечения всех линий, то

при возрастании угла р зависимости приобретают все более выраженный монотонный характер. При угле Р, стремящемся к нулю, значения ИКО конуса приближаются к значениям ИКО цилиндра.

При исследовании ЭКО установлено, что все зависимости фк (а) симметричны относительно значения фк (а = 90°), при этом фк (а = 0°) = фк (а = 180°). На рис. 3 представлены зависимости фк (а) для высоты конического КО над поверхностью Земли h = 600 км при различных значениях угла конусности р.

Характер зависимостей ЭКО от угла наклона оси а воспроизводится для различных высот. Различия заключаются в количественных значениях фк, падающих с ростом высоты. Из данных рис. 3 можно сделать следующие выводы:

- в области значений а =57° и а =123°, что примерно соответствует а = 1 рад и а = (л - 1) рад, наблюдаются пересечения всех линий (своеобразные узлы);

- в указанных узлах реализуется тепловое подобие конических объектов с любым значением Р;

- значение фк в узловых точках очень близко к ИКО сферического КО при той же высоте;

- в точке а = 90° имеют место экстремумы зависимости (максимумы при малых значениях р и минимумы при больших);

- при р < 15° максимальные значения ЭКО реализуются в точках а = 90°, а при больших значениях р наибольшие значения фк имеют место при вертикальной ориентации оси конуса - в точках а = 0° и

а = 180°.

фк 0,35

0,34

0,33

0,32

0,31

0,30

0,29

0,28

0,27

0,26

20 40 60 80 100 120 140 160 а°

Рис. 3. Зависимости ЭКО конической поверхности при высоте над поверхностью Земли h = 600 км от угла наклона оси конуса при р = 5°

и р = 45°

Тепловое подобие КО различных конфигураций

Тепловое подобие КО различных конфигураций определяется совпадением или весьма малым различием значений ЭКО, когда стационарные температуры таких КО в тени Земли равны.

На рис. 4 представлены зависимости ЭКО конуса и цилиндра от угла наклона их осей, а также ИКО сферы для малых величин угла конусности р. Расчеты проводились для случая, когда конус вписан в цилиндр, т.е. когда отношения высот обеих фигур к диаметру основания равны, что описывается соотношением

n = 0,5^р .

По данным расчетных исследований, частный результат которых представлен на рис. 4, можно сделать выводы, которые действительны и для других высот:

- при р < 11° ЭКО конуса с высокой точностью (погрешность менее 3%) совпадают с ЭКО цилиндра при условии, описываемом соотношением (8);

- в диапазоне изменения углов 18° < р < 21° зависимость фк (а) имеет пять экстремальных значений (при а, равном 0°, 40°, 90°, 140° и 180°), при этом колебания значений фк в диапазоне изменения а от 0° до 180° происходят относительно линии 4, соответствующей постоянному значению ИКО сферы;

- в результате можно сделать вывод о тепловом подобии конического и цилиндрического КО при р < 11° (что соответствует диапазону изменения n, примерно равному 2,8 < n < да), а также о тепловом подобии конического и сферического КО при 18° < р < 21°;

- в узловых точках а = 57° и а = 123° при всех высотах и углах при вершине конуса реализуется тепловое подобие конуса, цилиндра и сферы.

Рис. 4. Зависимости от угла наклона осей ЭКО конического (сплошные линии) и цилиндрического КО (штриховые линии), а также ИКО сферического КО (линия 4) при высоте орбиты 600 км.

Линии соответствуют: 1 - р = 5°; 2 - р = 10°; 3 - р = 19°; 4 - ИКО сферы

Подобие ЭКО конического КО сферическому создает большие удобства, поскольку ИКО сферы вычисляется по простой формуле [6]:

ф, = 0,5 (1 -^Тф).

Отмеченное тепловое подобие выполняется в указанных диапазонах с высокой точностью, это подобие сохраняется при увеличении высот КО над поверхностью Земли и свыше 2000 км, однако погрешность замены конуса цилиндром или сферой возрастает, а диапазоны изменения углов р, в пределах кото -рых реализуется подобие, сужаются.

Обсуждение результатов

Расчет ИКО конического КО базируется на использовании формул для ЛКО элементарной площадки для каждого участка боковой поверхности конуса (1)-(4), те же формулы используются для расчета ЛКО основания. В первом случае проводится интегрирование по угловой координате уравнений (1) и (2) в пределах, отличающихся от пределов, используемых при расчетах цилиндра. Затем определяется ЭКО конуса по формуле (7). Очень важно следить за выполнением условий (3), поскольку в некоторых случаях эти условия нарушаются, что приводит к разрывам зависимостей ИКО и ЭКО от угла ориентации оси а. Во избежание этой ситуации достаточно использовать в расчетах предельные значения У = +1 и г = 0.

На основе проведенных исследований установлено, что при равенстве отношения высоты к диаметру конуса и цилиндра (модель вписанного в цилиндр конуса) в пределах изменения половинного угла при вершине конуса 0 < р < 11° и высот 200 < к < 2000 км реализуется тепловое подобие конуса и цилиндра с погрешностью, не превышающей 3%. При изменении половинного угла при вершине конуса в пределах 18° < р < 21° на тех же высотах имеет место тепловое подобие конуса и сферы с той же погрешностью. При углах наклона оси конуса и цилиндра к линии зенит-надир примерно а = 57° и а = 123° цилиндр, конус и сфера с большой точностью подобны друг другу в тепловом отношении.

Заключение

Полученные результаты позволяют для ряда характерных случаев существенно упростить расчет эффективного коэффициента облученности конуса, заменяя его расчетом цилиндра, для которого ранее были получены простые в использовании аппроксимационные зависимости [5], а также расчетом интегрального коэффициента облученности сферы, определяемого по весьма простой формуле.

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

В отличие от цилиндрического космического объекта, для космического объекта конической формы весьма сложно подобрать достаточно простые и наглядные аппроксимационные зависимости эффективного коэффициента облученности от исходных параметров, которые позволили бы упростить расчеты. В определенной степени это объясняется введением дополнительного параметра - половинного угла при вершине конуса.

Reference

1. Modelirovanie teplovykh rezhimov kosmicheskogo apparata i okruzhayushchei ego sredy [Modeling of the thermal regimes of the spacecraft and its environment]. Ed. G.P. Petrov. Moscow, Mashinostroenie Publ, 1971, 382 p.

2. Duboshin G.N. Nebesnaya mekhanika. Osnovnye zadachi i metody [Celestial mechanics. The main tasks and methods]. Moscow, Nauka Publ., 1968, 800 p.

3. Chebotarev G.A. Analiticheskie i chislennye metody nebesnoi mekhaniki [Analytical and numerical methods of celestial mechanics]. Moscow-Leningrad, Nauka Publ., 1965, 367 p.

4. Smart W.M. Celestial Mechanics. NY, John Wiley, 1961, 381 p.

5. Dzitoev A.M., Hankov S.I. Metodika rascheta koeffitsientov obluchennosti tsilindricheskogo kosmicheskogo ob'ekta podsvetkoi Zemli [Calculation methods for irradiance coefficients of cylindrical space object by the Earth radiation]. Scientific and Technical Journal of Information Technologies, Mechanics and Optics, 2014, no. 1 (89), pp. 145-150.

6. Kamenev A.A., Lapovok Ye.V., Khankov S.I. Analiticheskie metody rascheta teplovykh rezhimov i kharakteristik sobstvennogo teplovogo izlucheniya ob"ektov v okolozemnom kosmicheskom prostranstve [Analytical methods for calculating the thermal conditions and the characteristics of the intrinsic thermal radiation of objects in near-Earth space]. St. Petersburg, L.T. Tuchkov Scientific Technical Center, 2006, 186 p.

7. Bayova Yu.V., Lapovok Ye.V., Hankov S.I. Metodika rascheta nestatsionarnykh temperatur kosmicheskogo ob'ekta na krugovykh orbitakh [A method for calculation of non-stationary temperature of space object on a circular orbit]. Izv. vuzov. Priborostroenie, 2013, vol. 56, no. №12, pp. 51-56.

8. Bayova Yu.V., Lapovok Ye.V., Hankov S.I. Metodika rascheta nestatsionarnykh temperatur kosmicheskogo ob"ekta, dvizhushchegosya po ellipticheskoi orbite [Calculation method of the transient temperatures for moving space object on the elliptical orbit]. Scientific and Technical Journal of Information Technologies, Mechanics and Optics, 2013, no. 6 (88), pp. 67-72.

9. Bayova Yu.V., Lapovok Ye.V., Hankov S.I. Analytical technique for calculating the heat flexes in near-Earth space that form the thermal regime of space telescopes. Journal of Optical Technology, 2013, vol. 80, no 5, pp. 283-288. doi: 10.1364/J0T.80.000283

10. Bayova Yu.V., Lapovok Ye.V., Khankov S.I. Metod podderzhaniya zadannogo temperaturnogo diapazona kosmicheskogo apparata, dvizhushchegosya po krugovoi orbite s zakhodom v ten' Zemli [Method of maintaining the specified temperature range of the spacecraft moving in a circular orbit with entering the Earth's shadow]. Izv. vuzov. Priborostroenie, 2013, vol. 56, no. 7, pp. 56-61.

11. Keihl J.T., Trenberth K.E. Earth's annual global mean energy budget. Bulletin of the American Meteorological Society, 1997, vol. 78, no. 2, pp. 197-208.

12. Trenberth K.E., Fasullo J.T., Keihl J. Earth's global energy budget. Bulletin of the American Meteorological Society, 2009, vol. 90, no. 3, pp. 311-323.

Дзитоев Азамат Миронович

Ханков Сергей Иванович

Azamat M. Dzitoev Sergey I. Khankov

Принято к печати 10.02.14 Accepted 10.02.14

адъюнкт, Военно-космическая академия имени А.Ф. Можайского, Санкт-Петербург, Россия, dzitoi8@gmail.com

доктор технических наук, ведущий научный сотрудник, Военно-космическая академия имени А.Ф. Можайского, Санкт-Петербург, Россия, Leva0007@rambler.ru

postgraduate, Mozhaisky Military Space Academy, Saint Petersburg, Russia, dzitoi8@gmail.com

D.Sc., chief staff scientist, Mozhaisky Military Space Academy, Saint Petersburg, Russia, Leva0007@rambler.ru