CC BY
151
ПЕЛЕВИН Федор Викторович (МГТУ им. Н.Э. Баумана)
PELEVIN Fedor Viktorovich
(Moscow, Russian Federation, Bauman Moscow State Technical University)
ПОНОМАРЕВ Алексей Викторович (МГТУ им. Н.Э. Баумана)
PONOMAREV Aleksey Viktorovich
(Moscow, Russian Federation, Bauman Moscow State Technical University)
СЕМЕНОВ Павел Юрьевич (МГТУ им. Н.Э. Баумана)
SEMENOV Pavel Yur'evich
(Moscow, Russian Federation, Bauman Moscow State Technical University)
УДК 621.454.2: 536.24
К вопросу о применении пористых металлов в регенеративной системе охлаждения жидкостных ракетных двигателей
Ф.В. Пелевин, А.В. Пономарев, П.Ю. Семенов
Требования по повышению надежности, ресурса, увеличению удельного импульса тяги жидкостного ракетного двигателя (ЖРД) обосновывают необходимость перехода к новым конструкциям и технологиям изготовления регенеративной системы охлаждения двигателя. В работе рассмотрена новая высокоэффективная система регенеративного охлаждения на основе принципа межканальной транспирации теплоносителя сквозь пористый сетчатый металл (ПСМ). Исследовано влияние числа каналов на гидравлические потери и эффективность теплообмена. Установлена высокая интенсификация теплообмена в пористом тракте при малых гидравлических потерях. Впервые определены оптимальные параметры регенеративной системы охлаждения ЖРД с межканальной транспирацией теплоносителя. На основе экспериментальных зависимостей по гидравлическому сопротивлению и теплообмену расчетно-экспери-ментальным путем показана высокая эффективность охлаждения ЖРД с использованием межканальной транспирации теплоносителя сквозь ПСМ. Представлена математическая модель и результаты расчетного исследования полей давления, температуры теплоносителя и пористого каркаса стенки камеры ЖРД.
Ключевые слова: регенеративная система охлаждения, пористый сетчатый металл, межканальная транспирация теплоносителя, интенсификация теплообмена, гидравлические потери, эффективность теплообмена.
On the application of porous metals in regenerative cooling systems of liquid propellant engines
F.V. Pelevin, A.V. Ponomarev, P.Yu. Semenov
To improve the reliability, service life, and specific impulse of a liquid propellant engine, new designs and manufacturing technologies of regenerative cooling systems are required. The paper deals with a new high-performance regenerative cooling system based on the principle of inter-channel coolant transpiration through a porous metal. The influence of the number of channels on the hydraulic losses and heat transfer efficiency is studied. The heat transfer intensification is established in a porous tract at low hydraulic losses. The optimum parameters of the regenerative cooling system of a rocket engine with inter-channel coolant transpiration are first determined. Using experimental dependences for hydraulic resistance and heat transfer, the high efficiency of liquid propellant engine
МАШИНОСТРОЕНИ
cooling with inter-channel coolant transpiration through a porous metal is demonstrated experimentally and by calculation. A mathematical model and the results of numerical analysis of pressure fields, coolant temperature, and porous engine chamber walls are presented. Porous heat exchange tracts with inter-channel coolant transpiration should be used as a regenerative cooling system in a liquid propellant engine and regenerative heat exchangers of various purposes.
Keywords: regenerative cooling system, porous metal, inter-channel coolant transpiration, heat transfer intensification, hydraulic losses, heat transfer efficiency.
ростом давления в камере сгорания и повышением коэффициента массового соотношения окислителя и горючего кт увеличивается не только удельный импульс двигателя, но и удельный тепловой поток q в стенку камеры жидкостного ракетного двигателя (ЖРД). Поэтому создание новых высокоэкономичных двигателей во многом зависит от эффективности системы регенеративного охлаждения камеры ЖРД.
В настоящее время в подавляющем большинстве эксплуатирующихся камер ЖРД используется оребренный тракт охлаждения, который обеспечивает интенсификацию теплообмена по сравнению с гладким каналом Nu/Nura в 1,5—3 раза [1]. Причем эффективность оребрения п уменьшается с ростом числа Рейнольдса. Применение искусственной шероховатости на дне канала в трактах охлаждения ЖРД также не дает высокого эффекта. Наибольшее значение теплогидравлической эффективности E достигается в области низких значений чисел Рейнольдса Re « 104, что для ЖРД с оребренными трактами не характерно. Такой уровень интенсификации теплообмена зачастую не обеспечивает удовлетворительного температурного состояния конструкции при ограниченных гидравлических потерях в регенеративной системе охлаждения. Поэтому для уменьшения удельных тепловых потоков в стенку двигателя практически всегда используется внутреннее (завесное) охлаждение, которое
обусловливает потери удельного импульса и усложняет конструкцию двигателя.
Необходимы новые конструкции и технологии изготовления теплообменных трактов, обеспечивающих надежное высокоэкономичное регенеративное охлаждение камеры ЖРД.
Теоретическое обоснование необходимости перехода от одномерного (продольно-канального) к двумерному (межканальному) движению теплоносителя сквозь пористый металл. Один из перспективных и эффективных методов интенсификации теплообмена заключается в использовании высокотеплопроводных пористых сетчатых металлов (ПСМ) в теплообменных устройствах. Такие металлы давно используются в ЖРД в качестве фильтров, капиллярно-заборных устройств и при изготовлении проницаемых огневых днищ смесительных головок. Однако в качестве теплообменного элемента (рис. 1) в проточной системе охлаждения ЖРД ПСМ не применяется из-за сверхвысоких гидравлических потерь, связанных с большими продольными габаритами камеры ЖРД, высокими коэффициентами гидравлического сопротивления и скоростями движения теплоносителя сквозь ПСМ, хотя, заполнение теплообмен-ного тракта пористым высокотеплопроводным металлом с малым термическим сопротивлением между стенкой и пористым металлом максимально интенсифицирует теплообмен [2—9]. Для уменьшения потерь давления следует снижать скорости движения теплоносителя в ПСМ за счет увеличения проходного сечения тракта, что приводит к уменьшению интенсификации теплообмена и повышению массы и габаритов тракта. Но есть и другие пути уменьшения потерь давления. Анализ модифицированного уравнения Дарси
Рис. 1. Регенеративная система охлаждения ЖРД с пористым наполнителем [2]:
1 — внутренняя стенка камеры; 2 — пористый металл
—йр / = аум + Рр^2 (1)
(р — давление; Z — координата; w — скорость движения; р — плотность; ц — коэффициент динамический вязкости; а, в — вязкостный и инерционный коэффициенты сопротивления ПСМ) показывает, что для снижения потерь давления в пористом теплообменном тракте, следует снижать скорость движения теплоносителя (на порядок и более), сократить путь движения теплоносителя сквозь ПСМ и уменьшить коэффициенты сопротивления ПСМ.
Для увеличения скорости движения теплоносителя через ПСМ и, соответственно, теплоотдачи, при заданных потерях давления, надо также уменьшать путь движения и коэффициенты сопротивления ПСМ.
Уменьшить потери давления в системе охлаждения без изменения габаритных размеров те-плообменного тракта можно, если перейти от общеизвестного продольно-канального к продольно-поперечному (межканальному) движению теплоносителя сквозь ПСМ, изготовленный методом диффузионной сварки в вакууме металлических тканых сеток. Будем называть этот способ движения теплоносителя межканальная транспирация теплоносителя (МКТТ) сквозь ПСМ [2].
Организовать МКТТ сквозь ПСМ можно следующим образом (рис. 2). Теплоноситель из
А-А
Рис. 2. Схема течения теплоносителя в тракте с МКТТ сквозь ПСМ:
1 — подводящий коллектор; 2 — наружная оболочка; 3 — сборный коллектор; 4 — внутренняя оболочка; 5 — подводящие каналы; 6 — отводящие каналы;
7 — ПСМ
подводящего коллектора 1 поступает в продольно расположенные в осесимметричной наружной оболочке 2 тракта подводящие каналы 5, заполняет их и под действием перепада давления движется сквозь ПСМ 7 в соседние отводящие каналы 6. Внутренняя оболочка 4 и торцы каналов непроницаемые для теплоносителя. Из отводящих каналов теплоноситель поступает в сборный коллектор 3, а в ЖРД — в смесительную головку. Подводящие и отводящие каналы чередуются и расположены симметрично относительно друг друга. Меняя число каналов, при заданном перепаде давления можно добиться требуемых скоростей движения теплоносителя сквозь ПСМ, не увеличивая толщину пористой вставки. Расстояние между подводящими и отводящими каналами может быть соизмеримо с толщиной пористого металла (3...5 мм).
При турбулентном режиме течения теплоносителя сквозь ПСМ, характерном для регенеративной системы охлаждения ЖРД, уравнение (1) принимает вид
Ар // = Ррг2.
Здесь / — путь движения теплоносителя сквозь ПСМ из подводящего в отводящий канал; V — скорость движения теплоносителя через ПСМ.
Выигрыш в скорости движения теплоносителя в тракте с МКТТ по сравнению с продольно-канальным движением при одинаковых значениях Ар, в, Р составляет
^мктт/Уп.к = (Ь //Г,
где Ь — продольный размер теплообменного тракта; — скорость продольно-канального движения теплоносителя.
Теплообменный тракт с межканальной транспирацией теплоносителя особенно эффективен при большой протяженности тракта охлаждения (Ь>>/), т. е. в области, где кольцевой тракт с ПСМ становится малоэффективным.
Увеличение теплоотдачи при одинаковых числах Яе и Рг составляет ^мктт/^и^к = (Ь//)0'4.
Важным преимуществом тракта с МКТТ является то, что площадь проходного сечения проходного тракта возрастает не за счет увеличения поперечных габаритов тракта, а за счет
большого числа подводящих и отводящих каналов N. Площадь проходного сечения тракта с МКТТ
Рпр = т&.
Здесь 5 — толщина пористого металла.
Для рассмотрения влияния числа каналов N на потери давления в тракте с МКТТ воспользуемся формулой (1). Для двух сравниваемых пористых теплообменных трактов с МКТТ секундный массовый расход теплоносителя т, толщину ПСМ 5, габариты тракта Б и Ь, коэффициенты сопротивления ПСМ а, в и температуру теплоносителя Тбудем считать постоянными, а число каналов N, расстояние между каналами /, площадь проходного сечения — переменными. Число каналов и расстояние между каналами связаны зависимостью
N1/1 = N2/2.
Для несжимаемой жидкости на ламинарном режиме течения применим закон Дар-си. Выражая скорость фильтрации через секундный массовый расход и учитывая, что = = NL5, получим
Ар =
ацт /
р5Ж'
Потери давления в пористых трактах с МКТТ
АР1 / Ар2 = N 22/ N12.
Для турбулентного режима течения несжимаемой жидкости
т 2р/
Ар =
р52 N2 Ь2'
Отношение потерь давления
АР1 / Ар2 = N 23/ N3.
Для расчета потерь давления при течении сжимаемой жидкости воспользуемся уравнением состояния идеального газа
Р:
р
ср
ЯТ
ср
где рср — осредненное давление жидкости в пористом металле.
Тогда для ламинарного режима
Р2ср Т1
АР1 / Ар 2 =
Р1срТ 2 N2
2
для турбулентного режима
N2 Р2срТ 1
АР1 / Ар 2 =
N1Р1 ср Т 2 '
Изменяя число кантов N, можно добиться требуемых скоростей фильтрации теплоносителя сквозь ПСМ и, следовательно, необходимых потерь давления, не увеличивая толщину пористого металла 5.
Из приведенных выше зависимостей следует, что для уменьшения затрат мощности на прокачку теплоносителя необходимо использовать тракты с максимально большим числом каналов.
Правильный выбор пористого материала для теплообменных трактов повышает их теп-логидравлическую эффективность.
Для эффективной работы теплообменного тракта с МКТТ пористый металл должен обладать следующими свойствами:
1) минимальным гидравлическим сопротивлением в направлении движения теплоносителя;
2) высокой теплопроводностью в направлении теплового потока;
3) высокой удельной прочностью, необходимой при изготовлении тонкостенных (2...5 мм) протяженных осесимметричных оболочек;
3) равномерной стабильной проницаемостью;
4) высокой технологичностью.
Всем этим требованиям отвечают ПСМ, изготовленные методом диффузионной сварки в вакууме металлических тканых сеток [3]. Эти металлы характеризуются прогнозируемой структурой, стабильной проницаемостью, высокой теплопроводностью в направлении слоев сеток, высокой прочностью и технологичностью [3]. Среди всех пористых металлов ПСМ имеет самые лучшие прочностные свойства.
Пористый сетчатый металл — это анизотропный материал. При одинаковой пористости гидравлическое сопротивление, теплопроводность, прочность вдоль слоев сеток и поперек плоскости сеток различны. В тракте с МКТТ
анизотропия свойств ПСМ полностью учитывается, поэтому ПСМ можно применять в регенеративной системе охлаждения ЖРД и в рекуперативных теплообменных аппаратах.
Регенеративная система охлаждения камеры ЖРД с МКТТ представлена на рис. 3.
5
4
2" 1
Рис. 3. Камера ЖРД с МКТТ:
1 — коллектор подвода охладителя; 2 — внутренняя
огневая стенка камеры; 3 — наружная силовая стенка камеры; 4 — ПСМ; 5 — подводящие и отводящие каналы
Результаты экспериментальных исследований гидравлического сопротивления и теплообмена в ПСМ. Анализируя результаты гидравлического сопротивления ПСМ при одномерной и двумерной фильтрации теплоносителя, следует отметить, что ПСМ, изготовленные из тканых фильтровых сеток, как говорилось выше, являются анизотропными материалами [2, 3]. Например, для ПСМ из тканой полотняной сетки П60 с пористостью П = 0,48 отношение вязкостных и инерционных коэффициентов ПСМ а±/ац и Р^/Рц при фильтрации теплоносителя ортогонально к плоскости сварки сеток (±) и при межсеточной фильтрации (||) составило а±/а || = 2,7 и Р^/Рц = 2,4.
Для каждого типа сетки (полотняного, саржевого переплетения и др.) получены свои экс-п е р и м е н тал ьн ы е з ав и с и м о с ти в я з ко с тн ы х и инерционных коэффициентов сопротивления ПСМ от пористости. Для ПСМ из полотняной сетки П60 получены следующие зависимости для вязкостного и инерционного коэффициентов сопротивления с восстановленной проницаемостью после механической обработки при межсеточной фильтрации теплоносителя:
а = (—5,148-П4 + 7,917-П3 - 4,409-П2 + + 1,034-П - 0,0824)-1013;
Р = (3,345-П4 - 5,232 П3 + 3,022 П2-
— 0,7718-П + 0,07531)-108. (2)
Здесь пористость ПСМ П = 0,291...0,476.
По сравнению с пористыми порошковыми металлами коэффициенты сопротивления ПСМ а и Р из тканой фильтровой полотняной сетки П24-П60 при одинаковой пористости и межсеточной фильтрации теплоносителя меньше в 10 и более раз [2, 3]. Это можно объяснить низкой шероховатостью (1-10-6 м) холоднокатан-ной проволоки, меньшей извилистостью поро-вых каналов и регулярностью структуры.
Такое значительное уменьшение коэффициентов сопротивления резко увеличивает эффективность теплообмена и выдвигает пористый теплообменный тракт с МКТТ и межсеточной фильтрацией теплоносителя на первый план.
На основании экспериментальных данных получено обобщающее критериальное уравнен и е п о в е р х н о с т н о й т е п л о о т д а ч и в т р акт е с МКТТ и двумерным движением теплоносителя в диапазоне изменения чисел Яе = = 20...2-104, чисел Рг = 0,7...7, пористости ПСМ П = 0,2...0,61, относительного пути движения теплоносителя //5 = 1,4...11,3. Критериальное уравнение теплоотдачи имеет вид
St = 0,57Яе"0'2Рг "°'7 г / . (3)
Здесь г / — коэффициент, учитывающий теплообмен на начальном участке; гх — коэффициент, учитывающий влияние теплопроводности материала сетки.
Из критериального уравнения (3) следует, что теплоотдача в пористом тракте с МКТТ в десятки раз больше, чем в гладком канале (трубе). Интенсификация теплообмена зависит от коэффициентов теплопроводности теплоносителя и ПСМ, а также от коэффициента г/, который учитывает увеличение теплоотдачи на начальном участке теплообмена в ПСМ. Как известно, интенсификация теплообмена происходит на начальном участке и в гладком канале.
Высокая интенсификация теплоотдачи в пористых трактах с МКТТ хорошо согласуется с данными работ [4—9].
Из характера зависимости (3) можно предположить, что зависимость будет действитель-
на и при больших числах Рейнольдса, так как развитый турбулентный режим движения полностью стабилизировался.
Влияние коэффициента теплопроводности теплоносителя на теплоотдачу учитывается с помощью числа Рг, а влияние коэффициента теплопроводности материала сетки X на теплоотдачу —
/ \ 0,4
с помощью коэффициента вх =ук/ Х12Х18Н10Т j ,
где Х12Х18Н10Т — коэффициент теплопроводности хромоникелевой стали 12Х18Н10Т.
Увеличение теплоотдачи с уменьшением отношения //Ъ учитывается введением поправочного коэффициента е/, который вычисляется по экспериментально полученной формуле [10]
в / = 0,5945 + 4,279(Ъ / /) + + 3,86(Ъ / / )2 — 4,995(Ъ / / )3,
справедливой в диапазоне толщин пористого металла Ъ = (2...5,6) -10-3 м и расстояний между каналами / = (7,97...33,77) -10-3 м.
Результаты моделирования температурного состояния камеры ЖРД с МКТТ. Особенность расчета ЖРД с МКТТ заключается в том, что расход охлаждающего компонента по длине регенеративного тракта охлаждения камеры переменен. При постоянном по длине камеры перепаде давления охлаждающего компонента между подводящими и отводящими каналами и неизменном числе каналов из-за переменного расстояния между каналами в критическом сечении камеры расход охлаждающего компонента будет наибольший. С увеличением расстояния между каналами расход компонента уменьшается, что положительно сказывается на равномерности температуры стенки камеры по длине ЖРД. Изменяя толщину ПСМ, пористость, тип сетки ПСМ по длине камеры, можно регулировать расход (скорость) компонента, а следовательно и температуру огневой стенки камеры. Варьируя этими параметрами можно добиться практически постоянной температуры стенки по длине камеры и тем самым уменьшить термическое напряжение в огневой стенке камеры.
На основании проведенных экспериментальных исследований был разработан метод поверочного расчета температурного состоя-
ния внутренней стенки охлаждающего тракта ЖРД с использованием МКТТ сквозь ПСМ.
Отличие от общеизвестного алгоритма расчета наружного охлаждения ЖРД заключается в нахождении количества компонента т,, протекающего через ,-й участок камеры.
По известной методике поверочного расчета охлаждающего тракта ЖРД [1], используя зависимость (3), определяется температура стенки камеры со стороны охлаждающего компонента (жидкости) Тст,ж;-, подогрев компонента на каждом участке разбиения А Т, суммарный подогрев компонента в тракте охлаждения ЖРД.
На основании методики расчета температурного состояния камеры ЖРД с МКТТ разработана компьютерная программа. Из примера расчета камеры разгонного двигателя, приведенного в таблице, следует, что при заданном постоянном перепаде давления 1,5 МПа компонента по длине камеры, числе подводящих и отводящих каналов N = 22, расход компонента О2ж переменен по длине камеры и зависит от расстояния между каналами и от значения коэффициентов сопротивления ПСМ а, в на участках разбиения камеры вдоль продольной оси (при одинаковых участках разбиения). При увеличении коэффициентов сопротивления ПСМ (пористость ПСМ уменьшилась с 0,29 до 0,23), снизился расход на 12-м участке и, следовательно, скорость компонента, а температура стенки камеры со стороны охлаждающего компонента возросла. Таким образом, можно получить почти постоянную температуру стенки камеры без резкого пика температуры в районе критического сечения.
Суммарный секундный расход компонента (жидкого кислорода) Ът = 10,0 кг/с; ТвхО2 =
= 108 К; ЪАТ = 151,14 К; ЪТвых = 182,86 К; средняя плотность, вязкость, теплоемкость жидкого кислорода: р = 752 кг/м3, ц = 0,4784-10-04 кг/(м-с), Ср = 2 074 Дж/(кг-К); ПСМ из полотняной сетки П60; материал сетки — хромистая бронза БРХ-08, П = 0,29; толщина пористой оболочки Ъ = 0,00289 м.
На основании результатов численного моделирования установлено, что минимальная температура огневой стенки при заданных потерях
Исходные данные и результаты расчета температурного состояния стенки камеры со стороны охлаждающего компонента
Номер участка АХ, м е-10-7, Вт/м2 а-10-10, —2 м р-10-5, м-1 м, м/с Яе-Ю-4 т, кг/с АТ, К Т.*, К
1 0,105 3,95 3,31 1,7 0,795 7,85 3,988 160,8 433,2
2 0,015 4,07 3,31 1,7 0,800 7,95 0,576 161,8 438,0
3 0,015 4,30 3,31 1,7 0,811 8,15 0,584 163,3 447,2
4 0,0155 4,77 3,31 1,7 0,829 8,55 0,637 166,9 465,9
5 0,012 5,47 3,31 1,7 0,866 8,77 0,538 164,2 483,9
6 0,010 6,51 3,31 1,7 0,910 9,22 0,486 164,4 514,9
7 0,010 7,44 3,31 1,7 0,950 9,41 0,473 161,5 537,9
8 0,010 7,74 3,31 1,7 0,971 9,27 0,466 156,4 539,4
9 0,010 5,64 3,31 1,7 0,953 7,40 0,490 130,4 448,4
10 0,0125 2,91 3,31 1,7 0,883 5,19 0,675 97,8 327,1
11 0,0175 1,69 3,31 1,7 0,790 4,17 0,870 85,1 271,9
12 0,016 1,01 150 15 0,226 3,19 0,216 227,4 456,7
давления в регенеративном тракте охлаждения устанавливается при расстоянии между каналами в районе критического сечения около 8...10 мм и высокотеплопроводном ПСМ на основе хромистой бронзы. Толщина пористого сетчатого металла соответствует оптимальной высоте оребренного тракта (3...4 мм). Для достижения оптимальной толщины ПСМ при заданном перепаде давления и числе каналов пористость ПСМ по длине камеры должна составлять 0,35...0,23.
Варьирование пористостью ПСМ по длине камеры ЖРД способствует выравниванию температуры огневой стенки по длине камеры и уменьшению термической деформации камеры. Низкие значения пористости ПСМ наблюдаются в области малых тепловых потоков. Низким значениям пористости ПСМ соответствует более высокая прочность пористого материала.
Математическая модель и результаты численных расчетов полей давления и температуры в пористом тракте с МКТТ. Для качественного понимания структуры двумерного движения теплоносителя и распределения температуры в ПСМ рассчитаны поля давления и температуры для ячейки тракта с МКТТ.
Когда отношение //5 мало течение сжимаемой жидкости в элементарной ячейке тракта с МКТТ носит двумерный характер. При расчете течения
и теплообмена в тракте с МКТТ необходимо использовать численное моделирование.
Рассмотрим пористое тело произвольной формы, частично заключенное в непроницаемую оболочку и имеющее включения из непроницаемого однородного материала. Часть внешней границы пористого тела нагревается. Охлаждение тела осуществляется прокачкой вязкого сжимаемого газа через пористую среду, движение газа происходит вследствие разности давления газа на входной и выходной границах. Пористая среда имеет анизотропную теплопроводность и проницаемость.
При построении математической модели процесса тепло- и массообмена в пористой среде для ее последующего численного решения методом контрольных объемов необходимо учесть следующие особенности:
1) вектор массовой скорости pw является фильтрационной скоростью и определяется как предел отношения секундного расхода теплоносителя через площадку, нормальную к направлению максимального расхода, к величине этой площадки при стремлении последней к нулю [10]. Поэтому при определении граничных условий на непроницаемых границах следует задаваться условием, что (pw,п) = 0, а не условием прилипания = 0;
МАШИНОСТРОЕНИ
2) при задании граничных условий для уравнения неразрывности на проницаемых участках внешней границы вместо задания расхода можно задавать распределение давления газа. Это вызвано тем, что дискретные аналоги для определения поля давления получаются интегрированием уравнения неразрывности;
3) граничные условия для уравнения энергии теплоносителя для соблюдения интегрального теплового баланса при численном решении уравнения энергии и уравнения теплопро-в о д н о с т и д л я кар кас а д о л жн ы з а д а в ат ь с я соотношением
((рс^ wTG — X0%га&Г0 ,л)| = /0(г,т).
Для расчета температурного состояния пористого элемента при двумерном движении теплоносителя необходимо решить уравнения [11]:
1 д
и = —-
ца1 + рР2|w| дх
1 др ца2 + рР2|w| ду
др 1 д(ри) _ д(рг) ду
— +
дт дх
+ Su; + Sv;
= 0;
(4)
(5)
(6)
(рс )G 8 Т + ¿1у((рс )G wTG — X G gгadTG ) =
= ку (тs —TG );
(7)
(рс) S (1 — г)Т + diy
Л дт
/
X1 0 0 X,
gгad Ts
= Ну (Тв — Т5 ) + ду
(8)
где w = и! + vj; SuSv — источниковые члены в уравнениях движения для и и V; и, v,TG ,Т8, р — неизвестные величины. Вязкость охлаждающего газа задана, плотность определяется из уравнения состояния. Уравнение (4) может рассматриваться, как уравнение для определения и — проекции вектора w на ось х, уравнение (5) — для определения V — проекции вектора w на ось у, уравнение (7) — для определения температуры охладителя TG, уравнение (8) — для определения температуры каркаса ПМ Т8. Для прямого определения давления
р уравнения нет, оно определяется косвенно через уравнение неразрывности (6).
В работе [11] приведено описание математической модели гидродинамики и теплообмена в ячейке тракта с МКТТ. Подводящие и отводящие каналы выполнены в ПСМ. Расчетная схема ячейки тракта с МКТТ представлена на рис. 4.
Рис. 4. Расчетная схема ячейки тракта с каналами в теле ПСМ:
р0 — давление на входе; р1 — давление на выходе
При моделировании процесса теплообмена при МКТТ сквозь ПСМ были заданы следующие параметры ПСМ:
— вязкостной коэффициент сопротивления а =3,97-109 м-2;
— инерционный коэффициент сопротивления Р =3,03104 м-1;
— теплопроводность ПСМ вдоль плоскости сеток X = 8,19 Вт/(мК).
Массовый расход газообразного кислорода через пористую вставку составил 0,620 кг/с. Поля давлений и температур при перепаде давления Ар = 1-105 Па изображены на рис. 5, 6.
Рис. 5. Поле давлений в ячейке тракта с МКТТ
На рисунке 5 видно, что застойные зоны наблюдаются вблизи теплоотдающей поверхности (огневой стенки) в зонах подвода и отвода теплоносителя и в теле ПСМ при непроницае-
мых боковых стенках подводящего и отводящего каналов.
На проницаемом дне подводящего канала температура теплоносителя минимальная (рис. 6, а). По мере приближения к отводящему каналу она постепенно растет. На дне отводящего канала наблюдаются потоки теплоносителя с разной температурой. У непроницаемой боковой стенки отводящего канала температура теплоносителя минимальная, а у центра канала — почти максимальная. Максимальная температура теплоносителя наблюдается в застойных зонах у теплоотдающей стенки напротив центров подводящего и отводящего каналов.
На рисунке 6, б видно, что поле температуры каркаса ПСМ в ячейке тракта качественно соответствует полю температуры теплоносителя. Отличие наблюдается только в численных значениях.
Рис. 6. Поле температуры в ячейке тракта с МКТТ; а — теплоносителя; б — каркаса ПСМ
Для снижения температуры и выравнивания ее полей необходимо использовать высокотеплопроводный ПСМ на основе сеток из хромистой бронзы.
Выводы
1. Теоретически обоснован переход от традиционного продольно-канального движения к продольно-поперечному движению охладите-
ля (межканальному движению) сквозь ПСМ с целью достижения максимально возможной интенсификации теплообмена в теплообменном тракте при минимальных потерях давления в нем.
2. Расчетно-экспериментальным путем определено температурное состояние камеры.
3. Установлены параметры тракта охлаждения ЖРД с МКТТ.
4. Переход от оребренного тракта охлаждения камеры ЖРД к тракту с МКТТ на основе ПСМ из хромистой бронзы позволяет уменьшить потери давления в тракте охлаждения ЖРД в 3 раза.
Литература
[1] Добровольский М.В. Жидкостные ракетные двигатели. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2005. 488 с.
[2] Леонтьев А.И, Пилюгин Н.Н., Полежаев Ю.В., Поля-ев В.М., ред. Научные основы технологий XXI века. Москва, УНПЦ ЭНЕРГОМАШ, 2000. 136 с.
[3] Резник С.В., ред. Материалы и покрытия в экстремальных условиях. Взгляд в будущее. В 3 т. Т. 2: Передовые технологии производства. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2002. 296 с.
[4] Зейгарник Ю.А., Поляев В.М. Теплообмен и гидродинамика двухфазных сред в условиях вынужденного движения в пористых структурах. Инженерно-физический журнал, 2000, т. 73, № 6, с. 1125-1134.
[5] Зейгарник Ю.А., Иванов Ф.П. К оценке теплогидрав-лических характеристик пористых структур. Тр. 5 Рос. нац. конф. по теплообмену. В 8 т. 2010, т. 5, с. 172-175.
[6] Кирсанов Ю.А., Назипов РА., Башкирцев Г.В. Теплообмен и сопротивление при течении однофазного теплоносителя в высокопористой вставке. Тр. 5 Рос. нац. конф. по теплообмену. В 8 т. 2010, т. 5, с. 176-179.
[7] Зейгарник Ю.А., Иванов Ф.П. Обобщение опытных данных по внутреннему теплообмену в пористых структурах. ТВТ, 2010, т. 48, № 3, с. 402-408.
[8] Кирсанов Ю.А., Назипов Р.А., Данилов В.А. Геометрические и теплофизические характеристики высокопористой структуры. Известия вузов. Авиационная техника, 2010, № 2, с. 49-52.
[9] Коротеев А.С., ред. Ядерные ракетные двигатели. Москва, ООО Норма-Информ, 2001. 416 с.
[10] Пелевин Ф.В., Лозовецкий В.В., Семенов П.Ю. Теплообмен и гидродинамика в пористых металлах при двумерном течении теплоносителя. XIV Минский международный форум по тепло-и массообмену. Минск, 10-13 сентября 2012. Тез. докл. и сообщ., т. 1, ч. 2, с. 577-579.
[11] Пелевин Ф.В., Лозовецкий В.В., Пономарев А.В., Семенов П.Ю. Моделирование гидродинамики и теплообмена в пористой среде при двумерном течении теплоносителя. 17-я Школа-семинар под руков. Академика А.И. Леонтьева. Проблемы газодинамики и тепломассообмена в аэрокосмических технологиях. Жуковский, 25-29 мая 2009 г., т. 1, с. 221-224.
References
[1] Dobrovol'skii M.V. Zhidkostnye raketnye dvigateli [Liquid rocket engine]. Moscow, Bauman Press, 2005. 488 p.
[2] Nauchnye osnovy tekhnologii 21 veka [Scientific bases of technology of the 21 century]. Ed. Leont'ev A.I., Piliugin N.N.,
Polezhaev Iu.V., Poliaev V.M. Moscow, UNPTs Energomash publ., 2000. 136 p.
[3] Materialy i pokrytiia v ekstremal'nykh usloviiakh. Vzgliad v budushchee [Materials and coatings under extreme conditions. Prospection]. Peredovye tekhnologii proizvodstva [Advanced production technology]. Ed. Reznik S.V. Vol. 2, Moscow, Bauman Press, 2002. 296 p.
[4] Zeigarnik Yu.A., Polyaev V.M. Heat Exchange and Hydrodynamics of Two-Phase Media under Conditions of Forced Motion in Porous Structures. Journal of Engineering Physics and Thermophysics. 2000, vol. 73, no. 6. pp. 1093—1101.
[5] Zeigarnik Iu.A., Ivanov F.P Kotsenke teplogidravlicheskikh kharakteristik poristykh struktur [Estimating the thermal and hydraulic characteristics of porous structures]. Trudy 5 Rossiiskoi natsional'noi konferentsii po teploobmenu [5 Proceedings of the Russian National Conference on Heat]. 2010, vol. 5, pp. 172—175.
[6] Kirsanov Iu.A., Nazipov R.A., Bashkirtsev G.V. Teploobmen i soprotivlenie pri techenii odnofaznogo teplonositelia v vysokoporistoi vstavke [Heat transfer and flow resistance in a single-phase coolant in highly porous insert]. Trudy 5 Rossiiskoi natsional'noi konferentsii po teploobmenu [5 Proceedings of the Russian National Conference on Heat]. 2010, vol. 5, pp. 176—179.
[7] Zeigarnik Yu.A., Ivanov F.P. Generalization of experimental data on internal heat transfer in porous structures. High Temperature, 2010, vol. 48, no. 3, pp. 382—387.
[8] Kirsanov Yu.A., Nazipov R.A., Danilov V.A. Geometric and thermophysical characteristics of a highly-porous structure. Russian Aeronautics, 2010, vol. 53, no. 2, pp. 191—197.
[9] Demianko Iu. G., Koniukhov G. V., Koroteev A. S., Kuz'min E. P., Pavel'ev A. A. Iadernye raketnye dvigateli [Nuclear rocket engines]. Moscow, OOO Norma-Informpubl., 2001. 416 p.
[10] PelevinF.V., Lozovetskii V.V., Semenov P.Iu. Teploobmen i gidrodinamika v poristykh metallakh pri dvumernom techenii teplonositelia [Heat transfer and hydrodynamics in porous metals with two-dimensional flow of the coolant]. 15 Minskii mezhdunarodnyi forum po teplo-i massoobmenu [15 Minsk International Forum on Heat and Mass Transfer]. Minsk, 10—13 September 2012, vol. 1, pt. 2, pp. 577—579.
[11] PelevinFV, Lozovetskii V.V, PonomarevAV., Semenov PIu. Modelirovanie gidrodinamiki i teploobmena v poristoi srede pri dvumernom techenii teplonositelia [Modeling hydrodynamics and heat transfer in porous media with two-dimensional flow of the coolant]. 17-ia Shkola-seminarpod rukovodstvom Akademika A.I. Leont'eva Problemy gazodinamiki i teplomassoobmena v aerokosmicheskikh tekhnologiiakh [17th School-workshop led by Academician Leontiev A.I. problems of gas dynamics and heat and mass transfer in aerospace technology]. Zhukovskii, 25—29 May 2009, vol. 1, pp. 221—224.
Статья поступила в редакцию 16.01.2014
Информация об авторах
ПЕЛЕВИН Федор Викторович (Москва) — доктор технических наук, профессор кафедры «Ракетные двигатели». МГТУ им. Н.Э. Баумана (105005, Москва, Российская Федерация, 2-я Бауманская ул., д. 5, стр. 1, e-mail: pelfv@rambler.ru).
ПОНОМАРЕВ Алексей Викторович (Москва) — аспирант кафедры «Ракетные двигатели». МГТУ им. Н.Э. Баумана (105005, Москва, Российская Федерация, 2-я Бауманская ул., д. 5, стр. 1, e-mail: 9722525@mail.ru).
СЕМЕНОВ Павел Юрьевич (Москва) — аспирант кафедры «Ракетные двигатели». МГТУ им. Н.Э. Баумана (105005, Москва, Российская Федерация, 2-я Бауманская ул., д. 5, стр. 1, e-mail: pavel.y.semenov@rambler.ru).
Information about the authors
PELEVIN Fedor Viktorovich (Moscow) — Dr. Sc. (Eng.), Professor of «Rocket Engines» Department. Bauman Moscow State Technical University (BMSTU, building 1, 2-nd Baumanskaya str., 5, 105005, Moscow, Russian Federation, e-mail: pelfv@rambler.ru).
PONOMAREV Aleksey Viktorovich (Moscow) — Post-Graduate of «Rocket Engines» Department. Bauman Moscow State Technical University (BMSTU, building 1, 2-nd Baumanskaya str., 5, 105005, Moscow, Russian Federation, e-mail: 9722525@mail.ru).
SEMENOV Pavel Yur'evich (Moscow) — Post-Graduate of «Rocket Engines» Department. Bauman Moscow State Technical University (BMSTU, building 1, 2-nd Baumanskaya str., 5, 105005, Moscow, Russian Federation, e-mail: pavel.y.semenov@rambler.ru).
