Научная статья на тему 'Теплообменная поверхность для подогрева кислорода в околокритическом состоянии в условиях невесомости'

Теплообменная поверхность для подогрева кислорода в околокритическом состоянии в условиях невесомости Текст научной статьи по специальности «Механика и машиностроение»

CC BY
0
0
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
теплообменник-подогреватель / теплоноситель / кислород в надкритическом состоянии / шаровая емкость / коэффициент случайной конвекции / число Фурье / число Нуссельта / коэффициент оребрения / heat exchanger-heater / heat carrier / oxygen in supercritical state / ball capacity / coefficient of random convection / Fourier number / Nusselt number / ribbing coefficient

Аннотация научной статьи по механике и машиностроению, автор научной работы — Светлов Юрий Валентинович, Пранцуз Оксана Сергеевна, Красниченко Александр Александрович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по механике и машиностроению , автор научной работы — Светлов Юрий Валентинович, Пранцуз Оксана Сергеевна, Красниченко Александр Александрович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Heat exchange surface for heating oxygen in the near-critical state in zero gravity

The problem of determining the large surface of a tubular finned heat exchanger located in a spherical vessel of a spacecraft (SC) filled with a single-phase liquid (oxygen) in a near-critical state under weightlessness conditions has been solved. The heat exchanger-heater is designed for constant pressure in the ball tank. The heat carrier circulating in the tube space of the heat exchanger is oxygen. The analysis of the influence of external disturbances on the heat transfer process in the heat exchanger-heater as well as on the propagation of heat in a homogeneous medium of supercritical state in zero gravity in a ball tank is carried out. The calculation was made for the most stressful conditions of the first hour operation of the system. The resulting heat exchange surface has a five-fold margin compared to the calculated one. In the future, in such cases it is possible to make a direct calculation of the required surface and find the convection coefficient taking into account the geometric structure of ribbing.

Текст научной работы на тему «Теплообменная поверхность для подогрева кислорода в околокритическом состоянии в условиях невесомости»

УДК 621.565.9

Теплообменная поверхность для подогрева кислорода в околокритическом состоянии в условиях невесомости

Д-р техн. наук Ю. В. СВЕТЛОВ1, канд. техн. наук О. С. ПРАНЦУЗ2, канд. техн. наук А. А. КРАСНИЧЕНКО3 Московский государственный университет технологий и управления имени К. Г. Разумовского 1yuvsvetlov@mail.ru, 2оksana.prantsuz@yandex.ru, 3a.krasnichenko@yandex.ru

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

Ключевые слова: теплообменник-подогреватель, теплоноситель, кислород в надкритическом состоянии, шаровая емкость, коэффициент случайной конвекции, число Фурье, число Нуссельта, коэффициент оребрения.

Информация о статье:

Поступила в редакцию 30.03.2023, одобрена после рецензирования 11.05.2023, принята к печати 22.05.2023 DOI: 10.17586/1606-4313-2023-22-3-29-36 Язык статьи — русский Для цитирования:

Светлов Ю. В., Пранцуз О. С., Красниченко А. А. Теплообменная поверхность для подогрева кислорода в околокритическом состоянии в условиях невесомости. // Вестник Международной академии холода. 2023. № 3. С. 29-36. DOI: 10.17586/1606-4313-2023-22-3-29-36

Heat exchange surface for heating oxygen in the near-critical state in zero gravity

D. Sc. Y. V. SVETLOV1, Ph. D. O. S. PRANTSUZ2, Ph. D. A. A. KRASNICHENKO3

K. G. Razumovsky Moscow State University of Technologies and Management 1yuvsvetlov@mail.ru, 2оksana.prantsuz@yandex.ru, 3a.krasnichenko@yandex.ru

The problem of determining the large surface of a tubularfinned heat exchanger located in a spherical vessel of a spacecraft (SC) filled with a single-phase liquid (oxygen) in a near-critical state under weightlessness conditions has been solved. The heat exchanger-heater is designed for constant pressure in the ball tank. The heat carrier circulating in the tube space of the heat exchanger is oxygen. The analysis of the influence of external disturbances on the heat transfer process in the heat exchanger-heater as well as on the propagation of heat in a homogeneous medium of supercritical state in zero gravity in a ball tank is carried out. The calculation was made for the most stressful conditions of the first hour operation of the system. The resulting heat exchange surface has a five-fold margin compared to the calculated one. In the future, in such cases it is possible to make a direct calculation of the required surface and find the convection coefficient taking into account the geometric structure of ribbing.

Keywords: heat exchanger-heater, heat carrier, oxygen in supercritical state, ball capacity, coefficient of random convection, Fourier number, Nusselt number, ribbing coefficient.

Article info:

Received 30/03/2023, approved after reviewing 11/05/2023, accepted 22/05/2023 DOI: 10.17586/1606-4313-2023-22-3-29-36 Article in Russian

For citation:

Svetlov Y. V., Prantsuz O. S., Krasnichenko A. A. Heat exchange surface for heating oxygen in the near-critical state in zero gravity. Journal of International Academy of Refrigeration. 2023. No 3. p. 29-36. DOI: 10.17586/1606-4313-2023-22-3-29-36

Введение

В сложных условиях космического полета, при отсутствии воздуха, теплообменное оборудование космических аппаратов (КА) включает разнообразные системы терморегулирования с различными газожидкостными контурами, излучательными радиаторами, нагревателями [1]-[8]. Циркуляция жидкости в теплообменных аппаратах позволяет практически ликвидировать тепловую стратификацию. Обеспечение долговечности и эффективности оборудования в системе жизнеобеспечения КА требует решения задач по определению эффективной площади поверхности трубчатых оребренных теплообменников, оптимальной конструкции, компоновки и тепловых режимов работы.

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

Вместе с тем существует ряд практических соображений, которые с достаточной степенью вероятности могут облегчить решение задачи получения необходимой теплообменной поверхности [9, 10]. Каким бы ни был космический аппарат (КА) по своему назначению, форме и сложности, он при работе неизбежно оказывает влияние на происходящие внутри него процессы и явления. Определенные импульсы возникают в результате работы системы ориентации, ударов частиц, внутренних механических воздействий, перемещений космонавтов и др. Как указывает Б. Гербгард (Gebgard) [9], даже в случае, когда такие импульсы являются следствием относительно малого ускорения электрической ракетной двигательной системы их влияние, как и влияние неравномерной тяги двигательной системы, может быть существенным. Суммарным результатом воздействия подобных возмущений на движение аппарата с управляемой ориентацией должны быть последовательные небольшие флуктуации мгновенной скорости и ориентация в пространстве. Таким образом, создаются реальные предпосылки для получения эффектов случайной (беспорядочной) конвекции. Анализ имеющихся литературных данных [10]—[13] позволяет рассчитывать, что случайные импульсы в условиях, когда затруднительно обеспечить вынужденное движение, как например, в шаровой емкости, должны оказывать некоторое влияние на скорость процессов переноса. Более того, существует предположение, что такого рода возмущения имеют порядок величины, сравнимый с любыми имеющими место установившимися «массовыми» силами.

Задачи исследования

Все сказанное выше относится к нашему конкретному случаю, когда передача теплоты производится к однофазной жидкости, находящейся в шаровой емкости, от теплоносителя, циркулирующего в трубчатом оребренном теплообменнике (рис. 1, 2). Ребра предна-

значены для более равномерного и эффективного съема теплоты.

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

Для решения задачи необходимо знать:

— температуру греющей поверхности в любой точке, т. е. характер распределения температур, как по поверхности трубки, так и в ребрах;

— закон распространения теплоты в гомогенной среде надкритического состояния в условиях невесомости.

Метод для оценки влияния возмущений на процесс теплопереноса

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

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

Средний тепловой поток в единицу времени и на единицу площади можно выразить следующим образом:

q = i о"

(1)

Q

где--количество теплоты, передаваемой в единицу

т

времени;

ф(т) — вероятностный закон распределения времени.

Nu = hs,

X

(3)

где 90 — разность между температурами стенки и среды; X — коэффициент теплопроводности.

Нетрудно видеть, что по физическому смыслу Л представляет собой коэффициент пропорциональности в уравнении Нуссельта а, а если учесть, что характерным размером для трубки будет ее диаметр, то получим известное выражение:

— ad Nu = —.

X

(4)

Рис. 1. Размещение трубчатой поверхности внутри шаровой емкости

Fig. 1. Placement of the tubular surface inside the ball tank

Разрез А—А i

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

В развернутом виде уравнение (3) запишется:

г д•и (у),

Л 0 у

Nu = -

dy,

(5)

где Бот — среднее или наиболее вероятное значение числа Фурье; у=¥о/¥от;

/ (у) — распределение вероятности ряда независимых переменных, в частности, представляется возможным использование экспоненциального закона / (у) = е-у или при условии, согласно которому общая вероятность равна 1.

/ (y) = V ^ -

(6)

Рис. 2. Разрез А-А шаровой емкости с размещенной в ней те-плообменной поверхностью Fig. 2. Section A-A of a ball tank with a heat exchange surface placed in it

Зависимость (1) выразим через число Фурье (Fo) при

наличии возмущении Fo = :

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

При п^-да все интервалы времени между возмущениями можно считать одинаковыми и характеризовать их величиной тт.

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

Опуская промежуточные рассуждения, укажем, что средние значения чисел Нуссельта должны быть рассчитаны по формуле (5), с учетом распределения вероятности (6).

Для теплопроводности в случае плоской поверхности получаем следующее выражение:

Nu =

2nn n ~ 0,5^. 0,5

я Fo0 nНо

! í

°yn~0,5e~'"dy = -C;

Fo0

(7)

— = 711 Fo ^ )d (Fo),

где Cn =

n0,5 (2n - l)(2n - 3) ...3 *1

(2)

n !2n

рутся из данных табл. 1.

или значения Cn бе-

где а — коэффициент температуропроводности (термической диффузии); 5 — характерный размер; ф (Бо) — вероятность распределения Бо.

Переходя к получению выражения для числа Нуссельта (^и), воспользуемся, введенным Гебгардом (Gebgard), понятием среднего коэффициента конвекции

и-зависимая константа

и-dependent constant

Таблица 1 Table 1

n 1 2 3 да

Cn 1 1,061 1,088 1,128

h = -

—— = —, таким образом

Для теплопроводности в случае цилиндрической поверхности используется выражение следующего вида:

- = П— Г(1

П! ^2

0,5™ 0,5 0,5

п , ^т У

)yne-nydy

или

Шп = - + -%-. п 2 Fo0,5

(9)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Преобразуя формулы (7) и (9) с использованием уравнения (3), получим выражения для случайного коэффициента конвекции (условного коэффициента теплоотдачи). Формулы Б. Гербгарда (Gebgard) [9] изначально приведены в оригинале, далее при расчетах для получения результатов в Международной системе (СИ), вводились коэффициенты перерасчета: — для плоской поверхности:

— I РСД

= С, ^ р

0,5

для цилиндрической поверхности:

^ = ^+С

РСРХ

0,5

(10)

(11)

К =

усл "1

-

13. Плотность р окружающей среды в сосуде в течение первого часа работы меняется от 1002 до 992 кг/м3 (при минимальном расходе р=581 кг/м3). Расчет производится для наиболее напряженных условий первого часа работы системы. Критериальные уравнения, а также общая методология теплового и конструктивного расчета теплообменных аппаратов, многократно проверенные для условий работы криогенных систем различного назначения [14, 15], остаются в основе настоящего расчета.

В качестве рабочих параметров принимается средняя температура потока в течение первого часа работы: Тср=(278 +130)/2=204 К.

Давление кислорода в трубках: Р1 = 58,84 бар. Удельный объем кислорода при Тср=204 К и Р1 =58,84 бар по -рР-диаграмме [17]: у=0,0072 м3/кг. Плотность равна: р! = 1/у= 139 кг/м3. Скорость газа в трубах w1 = G/ (р—). Принимаем диаметр трубки теплообменника й?н=8*1 мм, тогда площадь сечения по внутреннему диаметру трубки составит:

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

1 (12)

F =

= 28,210"6 м2;

0,275

=-«-

1 139 28,2 10"6-3600

= 0,0195 м/с.

а h

S

где -51 — коэффициент оребрения теплообменной поверхности;

£пл — поверхность ребер (плоская теплоотдающая поверхность), м2;

Sц — поверхность трубок (цилиндрическая теплоотдающая поверхность), м2.

Расчет теплообменника

Теплообменник — подогреватель располагается в шаровой емкости, заполненной однофазной жидкостью в надкритическом состоянии, и предназначается для поддержания постоянного давления в емкости. Исходные данные для расчета:

1. Характеристика емкости — объем 30 л, Д,н=385 мм.

2. Газ в теплообменнике — кислород.

3. Температура газа на входе в теплообменник — 278-298 К.

4. Температура газа на выходе из теплообменника 130 К.

5. Расход в течение первого часа работы — 0,275 кг/ч.

6. Минимальный расход — 0,076 кг/ч.

7. Средний расход через теплообменник — 0,132 кг/ч.

8. Расчетное давление в теплообменнике — 58,84 бар (60 кгс/см2).

9. Окружающая среда (содержимое емкости) — кислород.

10. Температура окружающей среды — 120-200 К.

11. Давление окружающей среды — 58,84 бар (60 кгс/см2).

12. Теплопритоки из окружающей среды — (1,4 ккал/ч) 1,63 Дж/с.

Физические константы кислорода при Тср=204 К и Р1 = 58,84 бар [14]-[16]:

— вязкость ц= 14710-7 Пас;

— теплопроводность Х=0,0233 Вт/ (мК). Число Рейнольдса:

Ке = = 0,0195-610~3 139 = 1106

ц 147 10"7 .

Для ламинарного режима течения потока в трубе (Яе < 2300):

D

1/3

№ = 1,61| Ре— | .

D

Поскольку Ре— < 12 , наступает режим тепловой

стабилизации потока и число Нуссельта принимает постоянное значение №=3,66. В этом случае

„ „„ X 3,66 0,0233 Вт а = 3,66-=-—-= 14,2—-—.

d„,

6 10"

Тепловая нагрузка на теплообменник: Я = 1,1629GAi = G (( - ^300 ) = = 1,1629 ■ 0,275 (93,16-17) = 24,42Вт.

Значение энтальпий определялись по диаграмме 1-1рР для кислорода [15].

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

т

т

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

В этих условиях можно считать, что число возмущений п^-да и тогда константа Сп= 1,128 (тем более, что значение Сп не слишком отличается при различных п). Очевидно, и время тт между отдельными возмущениями будет достаточно малым (если условно считать весь этот непрерывный процесс равномерно цикличным).

Принимаем тт=60 с=0,0167 ч. Физические константы кислорода при Т= 120 К и Р=58,84 бар:

— плотность р= 1000 кг/м3;

— теплоемкость Ср= 1843,6 Дж/ (кгК);

— теплопроводность Х=410,62 Дж/ (кгК).

Определяем коэффициент случайной конвекции:

а) для плоской поверхности по уравнению (10):

лпл = 1,128

1000 1843,6 ■ 410,62

0,5

= 238,83

кДж_

(11):

h = 1,189 + 238,83 = 240

кДж кг-К'

В связи с тем, что Лпл и кц практически одинаковы, в качестве расчетной осредненной величины примем

к = 239 КДЖ = 66,39- Вт

кг-К

кг ■ К

Таблица 2

Число рядов отверстий, в зависимости от диаметра пластин

Table 2

The number of rows of holes depending on the diameter of the plates

Диаметры пластин Число рядов отверстий

dj = 300 мм m1=d1/t=300/7=43

d2=360 мм m2=360/7=51

d3=380 мм m3=380/7=54

d4=380 мм m4=54

d5=360 мм m5 = 51

d6=300 мм m6=43

0,0167 ■ 41902 ) ' кг-К'

б) для цилиндрической поверхности по уравнению

Сечение отверстий на пластинах:

^1=5'6=14,1 10-6 1765=0,0249 м2; 5,2=^5=14,110 6 2263=0,0319 м2; ^з=5,4=14,110-62455=0,0346 м2.

Общее сечение отверстий:

£отв=0,0249 2+0,03192+0,03462 м2.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Пластины пронизываются трубками в 108 местах (3 трубчатых змеевика по 36 контактов) и припаиваются по периметру, образуя ребра. Таким образом, поверхность пластин, занятая трубками, будет равна:

£Тр - ■

яС

-■ 2 108 -1,57 ■ 82 10"6 108 - 0,0108 м2.

Поверхность пластин, не участвующая в теплооб-

мене:

В формуле (12) для определения условного коэффициента теплопередачи фигурирует коэффициент оре-

S

брения , который отражает соотношение между

плоской и цилиндрической поверхностями. Для получения величины указанного коэффициента задаемся геометрической формой теплоотдающей поверхности (рис. 1, 2). Она строится по принципу наиболее равномерного размещения в объеме сферы и состоит из 3-х трубчатых змеевиков, расположенных под углом 60° относительно друг друга, и шести поперечных пластин (ребер) толщиной 1-1,5 мм, которые припаиваются к трубкам. С целью уменьшения массы конструкции, а также улучшения условий для конвекции жидкости пластины перфорируются отверстиями 03 мм с шагом 7 мм (табл. 2). Исходя из принятой формы теплообменника, определяем его фактическую теплоотдающую поверхность.

Общее число отверстий на пластинах при условии их размещения по периметрам правильных шестиугольников с учетом заполнения сегментов [16].

^=^=1765;

N2=N5=2263;

N3=N4=2455.

Сечение одного отверстия:

Kd2

S = —о™ .2N = 1,57.32 .ю-6N = 14,110"6N.

S' = S0TB + STp = 0,183 + 0,0108 = 0,194 м2 .

Общая поверхность пластин без учета перфорации и трубок

Sобщ ~ 2

nd? „ nd2 „ nd2, „^ 1 -2 + —2 ■ 2 + —1 ■ 2

= 3,14(df + d22 + d32) =

4 4 4

v у

= 3,14(0,32 + 0,342 + 0,352) = 3,14(0,09 + 0,116 + 0,1225) = = 3,14■ 0,3285 = 1,08 м2.

Плоская теплоотдающая поверхность:

«пл = «общ -& = 1,08 -0,194 = 0,836 м2.

Длина трубки ¿=2280 мм. Цилиндрическая (змее-виковая) теплоотдающая поверхность

£ц=3£лй?тр. м= 17210 3 м2=0,172 м2.

Коэффициент оребрения:

«пл = = 4,86.

0,172

Общая фактическая теплоотдающая поверхность теплообменника:

Г = япл + = 0,836 + 0,172 = 1,008 м2.

Коэффициент теплопередачи:

K =_L

усл 1

1

- = 2,8-

Вт

4,86 + -14,2 66,39

В связи с существенным изменением удельной теплоемкости кислорода в зависимости от температуры

AT, = -

XI at 1= 0'234;

10

X

^ 0,234 'AT

Расчетная теплоотдающая поверхность: Q 24,42

^ = *

At

усл инт

2,8 ■ 42,7

= 0,204 м2.

Таблица 3

Средний температурный напор в теплообменном устройстве

Table 3

Average temperature pressure in the heat exchange device

№ участка D, мм AT, K 1/AT

1 339 143,8 0,00695

2 289 116,3 0,00860

3 228 89,93 0,01112

4 174 68,95 0,01450

5 131 52,21 0,01915

6 108 43,02 0,02325

7 99 39,70 0,02519

8 91 35,93 0,02783

9 74 28,59 0,03498

10 46 15,96 0,06266

Z 0,234225

Рис. 3. Определение средне интегральной разности температур графоаналитическим методом Fig. 3. Determination of the average-integral temperature difference by graphoanalytic method

в околокритической области, средний температурный напор в теплообменном устройстве определили как средне интегральное значение графоаналитическим методом [14, 15], показано на рис. 3. Результаты измерений приведены в табл. 3. Геометрические построения конструкций производились в масштабе М 1:1, и были разбиты на десять участков.

= 42,7К.

Фактически принятая поверхность имеет запас по сравнению с расчетной, равный = 1,008/0,204=4,84.

Заключение

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

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

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

Обозначения, применяемые в статье

Q — тепловой поток; Вт; Ц — средний тепловой поток; Q/т — количество теплоты, передаваемой в единицу времени; Тср — средняя температура потока, К; G — массовый расход, кг/ч; м> — скорость в трубке, м/с; i — энтальпия, кДж/кг; ф (т) — вероятностный закон распределения времени; Fo — число Фурье; Л — средний коэффициент конвекции, кДж/ (кгК); Л - коэффициент случайной конвекции для плоской поверхности, кДж/ (кгК); Лц — коэффициент случайной конвекции для цилиндрической поверхности, кДж/ (кг-К); d — диаметр пластин, м; т1 — число рядов пластин; а — коэффициент температуропроводности (термической диффузии), Вт/ (м2-К); 5 — характерный размер; ф (Бо) — вероятность распределения Бо; № — число Нуссельта; 00 — разность температур между температурами стенки и среды; X — коэффициент теплопроводности, Вт/ (мК); Ре — число Пекле; Re — число Рей-нольдса; N — число отверстий на пластинах; — площадь сечения одного отверстия, м2; 5общ — общая поверхность пластин, м2; 5" — поверхность пластин не участвующих в теплообмене; 5тр — поверхность пластин, занятая трубками, м2; 5пл — плоская теплоотдающая поверхность, м2; 5ц — цилиндрическая (змее-виковая) теплоотдающая поверхность, м2; — расчетная теплоотдающая поверхность, м2; ^ — фактическая теплоотдающая поверхность, м2; Кусл — коэффициент теплопередачи.

i=i

n

Литература

1. Igrickii V. A. Shaded spacecraft radiators to be used on the daytime surface of the mercury planet, the moon, and asteroids of the solar system inner part. // Aerospace Scientific Journal. 2016; 2 (6):69-93.

2. Актуальные вопросы проектирования автоматических космических аппаратов для фундаментальных и прикладных научных исследований / Сост. В. В. Ефанов. Химки: НПО им. С. А. Лавочкина, 2015. 350 с.

3. Финченко В. С., Котляров Е. Ю., Иванков А. А. Системы обеспечения тепловых режимов автоматических межпланетных станций / Под ред. В. В. Ефанова, В. С. Финченко. Химки: НПО им. С. А. Лавочкина, 2018. 400 с.

4. КокаревМ. А., Дзюбенко О. Л., Чмутин Е. В., Суязов Д. С. Перспективы применения новых кислородных систем в авиации // Современные научные исследования и инновации. 2016. № 4. с. 74-76.

5. Альтов В. В., Залетаев С. В., Копяткевич Р. М. Математическое моделирование тепловых режимов космического аппарата с помощью пакета программ ТЕРМ. // Космонавтика и ракетостроение. 2001. Вып. 23. с. 110-118.

6. Худяков С. А. Космические энергоустановки. М.: Знание, 1984. 64

7. Худяков С. А. Энергетические установки на основе топливных элементов для лунного орбитального аппарата и орбитального корабля «Буран». // Альтернативная энергетика и экология. 2020. (1-6). с. 122-124.

8. Jones B. G. Development for Space Use of BAe's Improved Single-Stage Stirling Cycle Cooler for Applications in the Range 50-80 K. Cryocoolers 8, New York, 1995. р. 1-11.

9. GebhardB. Random convection under conditions of weighless-ness. // AIAA Journal. 1963. Vol. 1. no 2. р. 380-383.

10. Козлов Л. В., Нусинов М. Д., Акишин А. И. и др. Моделирование тепловых режимов космического аппарата и окружающей его среды / Под ред. Г. И. Петрова. М: Машиностроение, 1971. 382 с.

11. Ягодников Д. А., Антонов Ю. В., Новиков А. В., Стрижен-ко П. П., Быков Н. И. Исследование процесса течения кислорода в рубашке охлаждения камеры ЖРД.// Вестник МГТУ им. Н. Э. Баумана. Серия: Машиностроение. 2014. № 6. стр. 3-19.

12. Ross R. G. Chapter 11: Aerospace Coolers: a 50-Year Quest for Long-life Cryogenic Cooling in Space, Cryogenic Engineering: Fifty Years of Progress, Ed. by K. Timmerhaus and R. Reed, Springer Publishers, New York (2007), pp. 225-284.

13. Стриженко П. П. Особенности расчета теплового состояния камеры ЖРД с беззавесным охлаждением жидким кислородом // Молодежь, техника, космос: Труды II обще-росс. молодеж. науч.-техн. конф. / Балт. гос. техн. ун-т. СПб., 2010. С. 140-142.

14. Малков М. П., Данилов И. Б., Зельдович А. Г., Фрадков А. Б. Справочник по физико-техническим основам глубокого охлаждения. М.: Госэнергоиздат, 1963. 416 с.

15. Епифанова В. И., Аксельрод Л. С. Разделение воздуха методом глубокого охлаждения. Технология и оборудование. М.: Машиностроение, 1973. 468 с.

16. Термодинамические свойства кислорода: ГСССД. Серия мо-нографии/Авт.: В. В. Сычев, А. А. Вассерман, А. Д. Козлов, Г. А. Спиридонов, В. А. Цымарный. М.: Издательство стандартов, 1981. 304 с.

References

1. Igrickii V. A. Shaded spacecraft radiators to be used on the daytime surface of the mercury planet, the moon, and asteroids of the solar system inner part. Aerospace Scientific Journal. 2016; 2 (6):69-93.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

2. Actual issues of designing automatic spacecraft for fundamental and applied scientific research / Comp. V. V. Efanov. Khimki: S. A. Lavochkin NGO, 2015. 350 p. (in Russian)

3. Finchenko V. S., Kotlyarov E. Yu., Ivanov A. A. Systems for providing thermal modes of automatic interplanetary stations / Edited by V. V. Efanov, V. S. Finchenko. Khimki: NPO Lav-ochkina, 2018. 400 p. (in Russian)

4. Kokarev M. A., Dzyubenko O. L., Chmutin E. V., Suy-azov D. S. Prospects for the use of new oxygen systems in aviation. Modern scientific research and innovation. 2016. No. 4. pp. 74-76. (in Russian)

5. Altov V. V., Zaletaev S. V., Kopyatkevich R. M. Mathematical modeling of thermal modes of the spacecraft using the TERM software package. Cosmonautics and Rocket Science. 2001. issue 23. pp. 110-118. (in Russian)

6. Khudyakov S. A. Space power plants. M.: Znanie, 1984. 64 p. (in Russian)

7. Khudyakov S. A. Power plants based on fuel cells for the lunar orbiter and the Buran orbiter. Alternative Energy and Ecology. 2020. (1-6): p. 122-124. (in Russian)

8. Jones B. G. Development for Space Use of BAe's Improved Single-Stage Stirling Cycle Cooler for Applications in the Range 50-80 K. Cryocoolers 8, New York, 1995. pp. 1-11.

9. Gebhard B. Random convection under conditions of weighless-ness. AIAA Journal. 1963. Vol. 1. No 2, p. 380-383.

10. Kozlov L. V., Nusinov M. D., Akishin A. I., etc. Simulation of thermal regimes of the spacecraft and its environment / Ed. G. I. Petrova. M: Mechanical Engineering, 1971. 382 p. (in Russian)

11. Yagodnikov D. A., Antonov Yu. V., Novikov A. V., Strizhen-ko P. P., Bykov N. I. Investigation of the oxygen flow process in the cooling jacket of the LRE chamber. Bulletin of the Bau-man Moscow State Technical University. Ser. Mechanical Engineering. 2014. No. 6, pp. 3-19. (in Russian)

12. Ross R. G. Chapter 11: Aerospace Coolers: a 50-Year Quest for Long-life Cryogenic Cooling in Space,» Cryogenic Engineering: Fifty Years of Progress, Ed. by K. Timmerhaus and R. Reed, Springer Publishers, New York (2007), pp. 225-284.

13. Strizhenko P. P. Features of calculating the thermal state of the chamber of a liquid rocket engine with self-cooling with liquid oxygen. Youth, technology, space: Proceedings of the II All-Russian Youth. sci.-tech. conf. / Baltic State Technical University. SPb., 2010. pp. 140-142. (in Russian)

14. Malkov M. P., Danilov I. B., Zeldovich A. G., Fradkov A. B. Handbook on the physical and technical foundations of deep cooling. M.: Gosenergoizdat, 1963. 416 p. (in Russian)

15. Epifanova V. I., Axelrod L. S. Air separation by deep cooling method. Technology and equipment M.: mechanical engineering, 1973. 468 p. (in Russian)

16. Thermodynamic properties of oxygen: GSSSD. Monograph series/Author: V. V. Sychev, A. A. Wasserman, A. D. Kozlov, G. A. Spiridonov, V. A. Tsymarny. M.: Publishing House of Standards, 1981. 304 p. (in Russian)

Сведения об авторах Светлов Юрий Валентинович

Д. т. н., профессор-консультант, кафедра пожарной безопасности и низкотемпературных систем, Московский государственный университет технологий и управления имени К. Г. Разумовского (ПКУ), 109004, Москва, ул. Земляной Вал, 73, yuvsvetlov@mail.ru

Information about authors

Svetlov Yuriy V.

D. Sc., Consultant Professor, Department of Fire Safety and Low-Temperature Systems, K. G. Razumovsky Moscow State University of Technologies and Management, 73 Zemlyanoi Val str., Moscow, 109004, yuvsvetlov@mail.ru

Пранцуз Оксана Сергеевна

К. т. н., доцент кафедры пожарной безопасности и низкотемпературных систем, Московский государственный университет технологий и управления имени К. Г. Разумовского (ПКУ), 109004, Москва, ул. Земляной Вал, 73, oksana.prantsuz@yandex.ru

Prantsuz Oksana S.

Ph. D., Associate Professor, Department of Fire Safety and Low-Temperature Systems, K. G. Razumovsky Moscow State University of Technologies and Management, 73 Zemlyanoi Val str., Moscow, 109004, oksana.prantsuz@yandex.ru

Красниченко Александр Александрович

К. т. н., доцент кафедры пожарной безопасности и низкотемпературных систем, Московский государственный университет технологий и управления имени К. Г. Разумовского (ПКУ), 109004, Москва, ул. Земляной Вал, 73, a.krasnichenko@yandex.ru

Krasnichenko Alexander A.

Ph. D., Associate Professor, Department of Fire Safety and Low-Temperature Systems, K. G. Razumovsky Moscow State University of Technologies and Management, 73 Zemlyanoi Val str., Moscow, 109004, a.krasnichenko@yandex.ru

@0®

Статья доступна по лицензии

Creative Commons «Attribution-NonCommercial»

Образовательный центр «Энергоэффективные инженерные системы»

Университета ИТМО

Приглашает принять участие в XI Международной научно-технической конференции

Искусственный холод в XXI веке 15-17 ноября 2023 г.

Секции конференции:

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

- Криогенная техника, водородные технологии и технологии СПГ.

- Системы жизнеобеспечения.

- Пищевые системы и консервирование холодом.

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

Регистрация участников - до 30 октября 2023 г., прием тезисов - до 15 июня 2023 г. Участие в конференции - бесплатно.

По итогам конференции будет издан Сборник тезисов конференции, индексируемый в РИНЦ.Материалы конференции, рекомендованные организационным комитетом, будут опубликованы в виде научных статей в изданиях, индексируемых в международной базе данных Scopus. Условия публикации будут сообщены дополнительно, в зависимости от уровня издания и его публикационной политики.

Информация о конференции на сайте: https://rft21.itmo.ru

i Надоели баннеры? Вы всегда можете отключить рекламу.