CC BY
80
УДК 520.224.2. 224.4
Ю. В. Баёва, Е. В. Лаповок, С. И. Ханков
МЕТОДИКА РАСЧЕТА НЕСТАЦИОНАРНЫХ ТЕМПЕРАТУР КОСМИЧЕСКОГО ОБЪЕКТА НА КРУГОВЫХ ОРБИТАХ
Получено аналитическое описание энергетического баланса изотермического сферического космического объекта с учетом высоты и угла наклона орбиты. Предложена методика, включающая определение интегральных коэффициентов облученности Землей и отраженным ею солнечным излучением. Разработана математическая модель для расчетов нестационарных температур таких объектов при попеременном заходе в тень Земли и выходе на освещенный Солнцем участок траектории.
Ключевые слова: космический объект, коэффициент облученности, тепловой режим, теплообмен излучением.
Введение. Определение нестационарных тепловых режимов космических объектов (КО) различного класса, например, космических аппаратов (КА), движущихся по разным траекториям вокруг Земли, является актуальной задачей. Расчет температурных полей реальных КА требует использования специальных программ расчетов. С точки зрения теории и практики очень важно исследовать общие закономерности формирования тепловых режимов идеализированных и обобщенных моделей КА, а именно КО простых типовых конфигураций.
Целями настоящей статьи являются разработка математической модели нестационарного теплового режима КО с учетом высоты и угла наклона орбиты, проведение расчетных исследований и получение на этой основе математической модели для тестирования более сложных расчетных программ. В качестве эталонной конфигурации принята сферическая форма КО. В математической модели особое внимание уделено циклограммам изменения интегральных коэффициентов облученности, которые ранее были введены без учета их изменения во времени на конкретных траекториях [ 1, 2].
Тепловая модель космического объекта, параметры его траектории. Для составления математической модели и получения аналитических формул для расчетов нестационарных температур КО введем ряд ограничений и допущений:
— рассматриваются круговые орбиты с различными углами наклона в, меньшими предельного значения во, т.е. орбиты, имеющие теневые участки;
— в расчетах Земля принимается изотермичным объектом с серой диффузной поверхностью; удельная мощность теплового излучения Земли принимается постоянной на ее теневой и освещенной Солнцем поверхностях, что соответствует общепринятой модели [1];
— в тепловом балансе КО не рассматриваются внутренние источники тепловыделений, тепловой режим КО формируется внешними тепловыми воздействиями: прямым солнечным излучением, собственным тепловым излучением Земли и отраженным ей солнечным излучением;
— поверхность КО, так же как и поверхность Земли, принимается серой и диффузной, при этом степени черноты и коэффициенты поглощения солнечного излучения не зависят ни от координат по поверхности, ни от температуры;
— пренебрегается рефракцией солнечного излучения при переходе КО через терминатор;
— рассчитывается среднеповерхностная температура КО в пренебрежении возможной температурной неравномерностью по поверхности и без учета его внутренней структуры, что соответствует либо сплошной структуре КО, либо оболочке с толстой стенкой из теплопроводного материала; такой подход соответствует принципу поэтапного моделирования,
согласно которому дальнейшая детализация температурных полей выполняется по результатам определения средних температур тел.
Схемы рассматриваемых траекторий КО представлены на рис. 1 (а — круговая траектория в направлении наблюдения от Солнца). Отрезок A'B' соответствует диаметру орбиты. Движение КО по круговой орбите АМВЫ, которая вследствие ее наклона видна в виде эллиптической проекции. Подсвеченный участок траектории находится в пределах части эллипса NAM. Теневой участок траектории заключен в пределах части эллипса NBM и угла 2ш. На рис. 1, б представлен вид сбоку, теневая часть Земли выделена темным цветом, орбита видна как проекция А'О'В', наклонение орбиты определяется углом р. Предельное значение угла наклона орбиты Ро соответствует величине угла наклона р, начиная с которого (при р > ро) КО никогда не заходит в тень Земли (рис. 1, в, г).
Рис. 1
Величину предельного угла нетрудно вычислить
ро = arcsin^/фО = 0о; Фо =
R
R + h
(1)
где Я — радиус Земли; И — высота орбиты над поверхностью Земли; фо — коэффициент об лученности Землей горизонтальной площадки, лежащей в плоскости местного горизонта (ло кальный коэффициент); 90 — половинный угол, под которым с КО видна Земля.
Период обращения КО вокруг Земли t0 зависит от высоты орбиты и описывается из вестным соотношением [3], которое для дальнейших расчетов удобно выразить через локаль
ный коэффициент облученности ф0 в виде
)3/2
t0 = 2л
(R + h)3
,1/2
2^4 = 5048ф0
Фо
3/4
c = 84,1фо3/4 мин
(2)
где g — ускорение свободного падения.
Для учета времени нахождения КО на солнечном и теневом участках траектории необходимо определить угол 2со (см. рис. 1). Из тригонометрических расчетов можно получить следующее соотношение:
(
2ш = 2 arcsin
Фо " sin р
Л
Фо
(1 - sin2 р)
при р < ро.
(3)
С учетом (1) из (3) следует: ш = 0 при Р = Р0. В диапазоне Р0 ^ Р ^ л/2 выполняется условие ш = 0- При Р > л/2 принимается схема орбиты, симметричная горизонтальной оси на рис. 1. В этом случае отсчет угла Р осуществляется зеркально схеме рис. 1, а.
При равномерном круговом движении КО время его нахождения в тени Земли t1 и на подсвеченном Солнцем участке траектории t2 может быть определено использованием схемы рис. 1 и соотношения (3):
t1 = 0,5(1 -t2 = 0,5(1 + n = 1 . (4)
\1 - sin Р
Оценим время нахождения КО в тени и на Солнце. При высоте орбиты h = 600 км (Р0 = 66°) и наклонении орбиты Р = 30° получим ф0 = 0,835, ш = 75,13°, n = 0,469 и t1 = 0,2655t0, t2 = 0,7345t0, при периоде обращения t0 = 96,3 мин t1 = 25,6 и t2 = 70,7 мин.
Тепловой баланс объекта в околоземном космическом пространстве. Из уравнений теплового баланса КО определяется его нестационарная среднеповерхностная температура, равная среднеобъемной для изотермического КО. В тепловом балансе должны быть учтены постоянная составляющая теплового потока от подсветки Земли, а также конкретная, зависящая от параметров орбиты, циклограмма изменения во времени прямой солнечной засветки и отраженного Землей солнечного излучения.
Уравнение теплового баланса КО (уравнение нестационарного теплообмена) может быть представлено в виде
Cs^ + S(1 -ф,) oT4 = QS, (5)
а т
C, = C0V/S = C0L; со = су,
Ql = s ф, Q+аАЕфь + а£ф/,
где T — среднеповерхностная температура КО; С, — поверхностная плотность полной теплоемкости КО; с0, с, у — соответственно объемная и удельная теплоемкость и плотность материала КО; s — степень черноты поверхности; V — объем КО, за исключением пустот в его теле, для сферической оболочки он равен объему самой оболочки; S — площадь наружной поверхности КО; L — определяющий размер КО, для сферической оболочки L = А, где А — толщина стенки КО [2]; ф, — интегральный коэффициент облученности сферического объек-
8 2 4
та Землей; o = 5,6710 Вт/м К — постоянная Стефана-Больцмана; Qi — сумма удельных тепловых потоков, поглощенных наружной поверхностью объекта в околоземном космическом пространстве; Q = 239 Вт/м — поверхностная плотность мощности, излучаемая Землей [4]; а, — коэффициент поглощения солнечного излучения покрытием на внешней поверхности наружного корпуса; А = 0,3 — альбедо Бонда Земли [4]; Е = 1366 Вт/м — солнечная постоянная [4]; фь — интегральный комбинированный коэффициент облученности внешней поверхности объекта отраженным Землей солнечным излучением; Ф — фактор формы, отношение площади миделя к общей площади поверхности КО (для КО сферической формы Ф = 0,25); f — функция, равная единице на освещенном Солнцем участке траектории и нулю — в тени Земли.
Интегральный коэффициент облученности объекта Землей известен [1, 2]:
2 (1 -тгф" ). (6)
Значение фь может быть рассчитано по формуле, выведенной с учетом данных [2]:
л ' t2 - 2(т-ti У
ф, =■ , ^ 2
фь = (1 -5)ф, cos Р cos
2
t2
5 = 0,25J-; H = h -10-3. (7)
Здесь t1 < т < t2. Моменты времени t1 и t2 соответствуют точкам траектории, при которых комбинированный коэффициент облученности ф^ уменьшается до нуля, а непосредственное воздействие на КО солнечного излучения остается постоянным до захода в тень Земли. В этом заключается различие в циклограммах воздействия на КО прямого и отраженного Землей солнечного излучения, что особенно проявляется при больших высотах, когда ^ >> t1, однако при этом условии значение фь пренебрежимо мало [2].
Опишем циклограмму воздействия на КО внешних тепловых потоков. В тени Земли поверхностью КО поглощается удельный тепловой поток
01 = в фs 0 при 0 < т < Ь. (8)
После прохождения теневого участка траектории суммарный поглощаемый поток равен & = 01 + Д0(т); Д0(т) = аАЕфь(т) + аДФ/ (9)
Поглощение потока прямого солнечного излучения продолжается в течение времени t1 < т < t2 (/= 1), начиная с момента времени т = t1 + ^ принимается / = 0.
Поглощение потока отраженного Землей солнечного излучения начинается после выхода КО из тени Земли с некоторым запаздыванием на отрезок времени Дт = (п/4)^. С этого момента величина комбинированного коэффициента облученности ф^ нарастает, достигая максимума при т = ^/2, а затем падает до малого значения к моменту времени т = ^ - (п/4)^. Далее циклы повторяются. Для расчетов t1, ^ и п используются формулы (2), (4).
Результаты расчетов. В расчетах использовались значения параметров: с = 922 Дж/кг К; у = 2640 кг/м3; с0 = 2,434-106 Дж/м3 К, что соответствует материалу КО (сплав АМГ6); в = 1; а, = 0,25 и 1. Толщина стенки сферического оболочечного КО А в расчетах варьировалась от 10-6 до 10-2 м. Расчеты проводились для орбиты с параметрами в = 30°, И = 600 км, ^ = 96,3 мин; t1 = 25,6 и ^ = 70,7 мин. Начальная температура КО принята равной 290 К.
Результаты расчетов периодических изменений поглощаемых тепловых потоков и температуры КО представлены на рис. 2, 3. На рис. 2 приведены циклограммы изменения во времени поглощаемых поверхностью КО тепловых потоков в пределах одного цикла. По оси абсцисс отсчитывается безразмерное время t = тЛ0, циклограмма построена в рамках одного цикла, т. е. полного периода обращения КО вокруг Земли. Линии соответствуют: 1 — поглощаемый удельный поток собственного теплового излучения Земли при в = 1; 2, 4 — к уровню линии 1 добавляется поглощаемый удельный поток от прямого солнечного излучения; 3, 5 — к уровню линий 2, 4 добавляется поглощаемый поток отраженного Землей солнечного излучения. Линии 2, 3 соответствуют коэффициенту поглощения поверхностью КО солнечного излучения а, = 1; 4, 5 — а, = 0,25. Вертикальные линии 6 определяют границы перехода через терминатор, расстояние между ними равно половине цикла.
На рис. 3 представлена динамика изменения температуры КО с различными теплоемко-стями в пределах одного цикла (а — а, = 1; б — 0,25). Графики построены для квазиустано-вившегося теплового режима КО, т.е. после прохождения участка переходного режима (линия 7 отсекает на оси абсцисс границу в пределах цикла между теневым и солнечным
-6.
участками траектории. Линии соответствуют толщине стенки сферического КО, равной 1 — 10 2 — 10-5; 3 — 10-4; 4 — 10-3; 5 — 10-2 м). Это означает, что на последующих циклах картина изменения температур полностью воспроизводится.
а)
Т, К
320
300
280
260
240
220
200
б)
Т, К
250
240
230
220
210
200
0 0,2 0,4 0,6 0,8 /*
Рис. 3
Из рис. 3 можно сделать общие выводы:
— безынерционному в тепловом отношении объекту соответствует толщина стенки сферического КО до 10-5 м (линии 1 и 2), температурная динамика таких КО полностью воспроизводит циклограмму изменения поглощаемых тепловых потоков, а амплитуда колебаний их температур максимальна — 132 К при а^ =1 и 44 К при а^ = 0,25. Нижний уровень температуры (в тени Земли) безынерционных КО одинаков (он зависит от в) и при в = 1 составляет - 206 К;
— при толщине стенки сферического КО 10-4 м объект нельзя считать безынерционным, хотя с некоторым запаздыванием (на ^/2 и ^/2) изменение температуры соответствует изменению во времени поглощаемых тепловых потоков. Такой КО можно считать малоинерционным;
-3
— при толщине стенки 10 м объект является термоинерционным, размах колебаний его температуры уменьшается в основном за счет повышения температурного уровня на теневом участке траектории (процесс охлаждения), на солнечном участке траектории динамика изменения температуры запаздывает относительно циклограммы изменения тепловых потоков и максимум смещен к концу цикла;
-2
— при толщине стенки 10 м термическая инерция КО велика, амплитуда колебаний температуры составляет всего 18 К при а^ =1 и 4 К при а^ = 0,25. Очевидно, что при дальнейшем увеличении толщины стенки КО можно будет отнести к классу термостабилизиро-ванных объектов, т.е. таких, у которых колебание температуры относительно среднего уровня можно считать практически незначимым.
Выводы. Разработана математическая модель для расчета нестационарных температур изотермического сферического КО с учетом высоты и наклона орбиты. Модель включает описание циклограмм поглощаемых поверхностью КО всех тепловых потоков в околоземном космическом пространстве. Проведенные расчеты позволили выявить закономерности формирования нестационарных тепловых режимов КО. Установлены критерии отнесения объектов к классу безынерционных или термостабилизированных (с большой тепловой инерцией). Разработанная методика расчета может применяться и для КО с другими конфигурациями (цилиндр, параллелепипед) после уточнения формул для коэффициентов облученности.
Модель нестационарного теплового режима изотермического сферического КО может использоваться как эталонная, она удобна для тестирования более сложных расчетных программ.
СПИСОК ЛИТЕРАТУРЫ
1. Моделирование тепловых режимов космического аппарата и окружающей его среды / Под ред. Г. И. Петрова. М.: Машиностроение, 1971. 382 с.
2. Каменев А. А., Лаповок Е. В., Ханков С. И. Аналитические методы расчета тепловых режимов и характеристик собственного теплового излучения объектов в околоземном космическом пространстве. 2006. 186 с.
3. Богомолов С. И. Введение в специальность „Радиосвязь, радиовещание и телевидение": Учеб. пособие. Томск: ТМЦДО, 2004. 162 с.
4. Trenberth K. E., Fasullo J. T., Keihl J. Earth's global energy budget // Bull. Amer. Meteor. Soc. 2009, Vol. 90, N 3. Р. 311—323.
Юлия Валерьевна Баёва
Евгений Владимирович Лаповок
Сергей Иванович Ханков
Рекомендована кафедрой компьютерной теплофизики и энергофизического мониторинга
Сведения об авторах аспирант; Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кафедра компьютерной теплофизики и энергофизического мониторинга; E-mail: [email protected]
канд. техн. наук; Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кафедра компьютерной теплофизики и энергофизического мониторинга; E-mail: [email protected]
д-р техн. наук; Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кафедра компьютерной теплофизики и энергофизического мониторинга
Поступила в редакцию 15.01.13 г.
