4.3.2. Расчет плотности падающего на поверхность ка потока исходящего от планет излучения

Рассмотрим методический подход к расчету плотности поглощаемого поверхностью КА потока исходящего от планеты излучения на примере Земли. Обычно при расчете и физическом моделировании падающего на КА излучения от планет используются их модели, основанные на осреднении радиационных характеристик поверхности планеты по достаточно большой области и даже по всей ее поверхности [7,9,10]. При этом планета рассматривается как диффузное сферическое тело с эффективным радиусом , где - радиус планеты, а - высота верхней границы эффективно излучающего и отражающего слоя атмосферы, если таковая имеется. Для Земли обычно принимают =12 , а . Плотность потока собственного излучения Земли принимается постоянной по поверхности и определяется из теплового баланса Земли в целом:

,

где - среднее планетарное альбедо Земли. Чаще всего задают равным 0,35. Тогда . Предполагают также, что Земля излучает как абсолютно черное тело при соответствующей величине температуре. Если , то радиационная температура Земли составит 253 , т.е. -20 ^○ . При такой модели в каждой точке орбиты КА поле интенсивности собственного излучения планеты в пределах телесного угла ее обзора является изотропным. Условимся называть такую модель лучистого поля планеты изотропной моделью.

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

,

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

Условимся называть модель поля отраженного от планеты излучения анизотропной моделью. Кратко рассмотрим методический подход к определению поглощаемого каким-то элементом КА или его какой-то поверхностью потока исходящего от Земли излучения - собственного и отраженного солнечного.

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

Алгоритм расчета , разработанный Залетаевым В.М. [ 7 ], приводится ниже. В этом алгоритме угловой коэффициент представлен в виде функции двух параметров - углов и , где - угол между нормалью к -ому тепловоспринимающему элементу и местной вертикалью , - линейный угол, равный углу полураствора конуса обзора планеты (см. приведенный ниже рисунок 4.2.), причем , - высота орбиты.

В соответствии с [7], при

= ;

при

;

при .

На рисунке 4.3 представлены результаты расчета зависимости = .

Средняя плотность поглощаемого какой-то однородной поверхностью потока собственного излучения Земли ( определяется с помощью следующего выражения:

,

где - средний угловой коэффициент поверхности и планеты. .

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

Локальная ( ) и средняя ( плотность поглощаемого элементом или поверхностью КА потока отраженного от Земли солнечного излучения определяется с помощью следующих соотношений :

, ,

где и - соответственно локальный и средний обобщенные угловые коэффициенты [7].

Определение коэффициентов и является более сложной задачей, нежели определение коэффициентов и . Дело в том, что обобщенные угловые коэффициенты зависят от большего числа параметров, чем коэффициенты и . Методика расчета обобщенных угловых коэффициентов излагается в [7]. Хотя реализация этой методики при наличии соответствующей компьютерной программы не является трудоемкой задачей, однако, для околоземных орбит высотой менее 500 (а именно такие орбиты представляют наибольший интерес) можно рекомендовать упрощенную методику расчета величин и . Рекомендуется определять и с помощью следующих приближенных соотношений:

, .

Использование этих соотношений основано на предположении о равномерном распределении плотности потока излучения, отраженного от видимой с КА области планеты, равной плотности потока отраженного излучения, уходящего от того участка земной поверхности, над которым в данный момент находится КА. То есть фактически вводится предположение об изотропности поля интенсивности излучения в пределах телесного угла обзора планеты. Сказанное иллюстрируется на рисунке 4.4. На этом рисунке линия 1 характеризует фактическое распределение плотности потока отраженного солнечного излучения (разумеется, в соответствии с принятой моделью). Линия 2 на рисунке характеризует распределение плотности отраженного от Земли солнечного излучения при наличии предположения об изотропности в пределах телесного угла обзора планеты поля интенсивности излучения.

Рис.4.1.-4.4. Качественный характер распределения по поверхности планеты плотности потоков собственного и отраженного солнечного излучения.

Рисунки 4.1- 4.4

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

Рис. 4.5. Зависимость температурных погрешностей, обусловленных заменой анизотропного лучистого поля изотропным, от высоты орбиты при различных значениях зенитного угла солнца и различной ориентации тепловоспринимающих серых элементов.