Энергетика тангенциального подводящего устройства микротурбины системы терморегулирования перспективного космического аппарата
- Авторы: Шевченко Ю.Н.1, Кишкин А.А.1, Зуев А.А.1, Делков А.В.1
-
Учреждения:
- Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева
- Выпуск: Том 25, № 1 (2024)
- Страницы: 126-142
- Раздел: Раздел 2. Авиационная и ракетно-космическая техника
- Статья опубликована: 29.07.2024
- URL: https://journals.eco-vector.com/2712-8970/article/view/634627
- DOI: https://doi.org/10.31772/2712-8970-2024-25-1-126-142
- ID: 634627
Цитировать
Аннотация
В настоящей работе представлен обзор современной технической проблемы, связанной с двухфазными системами терморегулирования космических аппаратов, и возможные технические приложения рекуперации тепловой энергии в органическом цикле Ренкина как составной части систем обеспечения теплового режима. Конструктивное решение подразумевает собой интегрирование паровой микротурбины за радиатором-испарителем. Микротурбина представляет собой тангенциальное подводящее устройство и радиально центростремительное рабочее колесо низкой быстроходности nst<40. В этой области не существует достоверных данных по проектированию и энергетики как подводящего устройства, так и рабочего колеса. Энергетика (потери энтальпии) подводящего устройства определяет в основном транспорт закрученного потока к рабочему колесу и, как следствие, окружную работу на турбине. Разработан и представлен прототип радиальной микротурбины с целью оценки конструктивного исполнения проточной части как подводящего устройства, так и рабочего колеса. В результате анализа выделены основные определяющие гидродинамические участки, необходимые для гидродинамического анализа и математической проработки алгоритма расчета течений с оценкой энергетических потерь: течение закрученного потока радиально-кольцевой щели; осе-кольцевой щели и тангенциального подводящего устройства. Первые два алгоритма предполагают расчетное моделирование. Модель энергетических потерь в тангенциальном подводящем устройстве не поддается аналитическому моделированию, поскольку включает в себя последовательность (или совместность) течения в граничных условиях, определяемых как «местные сопротивления»: внезапное расширение, разворот потока, совместно с участком радиально окружного течения. Взаимовлияние этих граничных условий предполагает только экспериментальную оценку энергетических потерь в тангенциальном подводящем устройстве через коэффициент потерь местного сопротивления в диапазоне изменения геометрических и режимных параметров.
В результате экспериментальных исследований предложена база данных по коэффициенту потерь тангенциальных подводящих устройств микротурбины в области практического диапазона существования режимных и конструктивных параметров.
Полный текст
Введение
Анализ информации о современных зарубежных разработках в области использования двухфазных контуров (ДФК) на космическом аппарате (КА), имеющих высокое энергопотребление, показывает техническую перспективность использования их в системах терморегулирования. В числе КА с использованием фазового перехода в системах терморегулирования упоминаются «Эврика», «Колумб» NASA, на которых были успешно проведены натурные испытания высокоэффективных двухфазных модулей. На КА COMET (NASA) обеспечение теплового режима индивидуальных теплонапряженных блоков полезной нагрузки осуществляется с помощью систем терморегулирования (СТР) на капиллярных насосах. На КА STENTOR (NASA) в СТР также используются капиллярные насосы в раскрывающейся конструкции излучательного радиатора. Проведенные натурные испытания на этих КА показали высокую эффективность СТР на базе двухфазных контуров.
Применение генерации с использованием контуров на низкокипящих рабочих телах на кораблях для дальних космических полетов для их энергообеспечения было предложено в СССР еще в 70-е гг. прошлого века. В качестве источника тепла для силового паротурбинного цикла использовались унифицированные модули радиоизотопного источника (Pu-238) [1]. Результатом исследований стала установка номинальной мощности 1,3 кВт, массы 205 кг и оснащенная холодильником-излучателем площадью 10,8 м2. На стендовых испытаниях массовый расход через контур достигал 0,0107 кг/с. Установка проработала 20 ч. За это время существенных отклонений параметров не было выявлено.
Для решения задач по терморегулированию перспективных крупногабаритных КА и станций с увеличенным энерговооружением при необходимости снижения относительных массы и размеров СТР КА наиболее перспективный путь – это создание базовых элементов интегрированных СТР КА с ДФК, как наиболее эффективной по сравнению с существующими в настоящее время в отечественной и в зарубежной практике средствами терморегулирования. При существующих системах необходимо значительно увеличивать площади наружных радиаторов-излучателей, что приводит к увеличению их массы и габаритов. Совмещение системы терморегулирования с паровой турбиной позволяет отводить часть энергии через паровую турбину в систему энергоснабжения, что уменьшает тепловую нагрузку на радиатор-излучатель. Проблема радиаторов-излучателей возникает при ограниченных возможностях размещения их под обтекателем ракеты-носителя на участке выведения КА на орбиту.
Основное техническое предложение
Содержание технического решения определяется в соответствии с источниками [2–4]:
– прямой цикл на низкокипящих телах с выработкой механической энергии. В качестве источника энергии используется внутренний источник тепла КА – тепло, выделяемое работающими приборами и установками. В традиционной схеме это тепло отводится в космос с помощью холодильников-излучателей. Предлагаемый вариант позволяет сэкономить энергоресурсы и вырабатывать механическую энергию, а также сократить площади холодильников-излучателей;
– прямой цикл на низкокипящих телах с выработкой механической энергии с использованием в качестве источника солнечного тепла, снимаемого с солнечных панелей. В связи со сравнительно низкими равновесными температурами панелей имеет низкую эффективность, однако позволяет использовать солнечное излучение в более широком инфракрасном (тепловом) диапазоне. Установки, изготовленные по этому способу, не подвержены старению и деградации характеристик (по сравнению с солнечными батареями);
– то же, что и в предыдущем пункте, с применением солнечных концентраторов параболического типа. Позволяет получить высокие температуры, а значит и высокую эффективность установки. Однако в этом случае необходимо решить проблему высоких перепадов давления рабочего тела в системе.
Конструктивная схема микротурбины
В теории турбостроения в основном используется соотношение для расчета и проектирования ступени турбин высокой быстроходности как осевой, так и радиально-центростремительной схемы [2; 5–9]. Этой теме посвящены большая часть литературы, касающейся различных отраслей турбостроения: турбины большой электрогенерации (электростанции), турбины транспортных устройств (летательные аппараты, наземный, железнодорожный транспорт, различные устройства специального назначения и т. п.). Большая часть турбин проектируется как ступени активного типа большой мощности, больше 100 кВт. В области распределенной энергетики и практического использования низкопотенциального бросового тепла используются турбины радиально-центростремительного типа мощностью менее 100 кВт с коэффициентом быстроходности [4; 6; 10]. В основном применяется соотношение газовой динамики на сверхкритических (сверхзвуковых) перепадах давлений и температур. Тема расчетов конструкций и оптимизации параметров изложена достаточно объемно. Для диапазона низкой быстроходности nst < 60 тема недостаточна раскрыта в части инженерных расчетных приложений, что вызывает очевидные трудности при моделировании подобных объектов мощностью меньше 100 кВт. На сегодняшний момент для утилизации энергии газотранспортных систем и систем добычи углеводородов используются турбины различных специальных конструкций: в том числе вихревые, безлопаточные – центробежные и т. п. [7].
В нашем случае турбины мощностью менее 1 кВт можно идентифицировать (классифицировать) как микротурбины малой быстроходности, размерности и мощности [11]. Вопрос оптимального проектирования и выбора типа турбины остается открытым из-за низкого КПД тихоходных турбомашин, имеющего сходные значения для микротурбин лопаточно-осевого и радиального типа, лабиринтно-вихревых, дисковых турбин, турбин центробежных и центростремительного типа и т. п. Отдать предпочтение какому-то типу на данном этапе невозможно. Стоит отметить, что для турбин как активного, так и реактивного типа важнейшим элементом, формирующим окружное направление потока, обеспечивающего окружную работу рабочего колеса, является сопловой или направляющий аппарат для быстроходных турбин, выполненных в виде сопловых решеток (в лопаточном венце). Для тихоходных и малорасходных машин выполняется одно сопло (сопловой тангенциальный канал) [2]. Диапазон выходной мощности по техническому заданию заказчика составляет от 100 до 1000 Вт на дозвуковом перепаде рабочего тела (пара), при перепаде температур не более 60 ºС с ограничением оборотов до 5000 об/мин из-за требований высокого ресурса. При таких данных не стоит ожидать высоких КПД, сравнимых с КПД турбин большой энергетики. Задача оптимального проектирования с целью наивысшего достижения КПД является актуальной. Для постановки задач исследования спроектирован и изготовлен прототип микротурбины с генератором с характерными элементами конструкций на базе центростремительного рабочего колеса детандера (рис. 1).
Рис. 1. Прототип микротурбогенератора: 1 – крышка; 2 – корпус микротурбины с тангенциальным подводящим устройством; 3 – рабочее колесо; 4 – генератор
Fig. 1. Prototype of a microturbo generator: 1 – cover; 2 – microturbine housing with tangential supply device; 3 – impeller; 4 – generator
Из конструкции видно, что тангенцильный канал является основным элементом, формирующим окружное течение потока. Конструктивно следует выделить область радиально-окружного потока RU и осе-окружного потока LU. Общий вид турбогенератора показан на рис. 2.
Рис. 2. Вид на рабочее колесо со стороны крышки (а); крышка и корпус генератора (б); вид со стороны крышки генератора (в)
Fig. 2. View of the impeller from the cover side (a); сover and housing of the generator (b); view from the side of the generator cover (c)
Цель и постановка задачи исследования
Существующие на сегодняшний день конструкции СТР автономных КА с длительным сроком активного существования (до 15 лет) используют практически весь спектр известных конструкторских решений: естественную теплопроводность от источника до радиатора-излучателя в черный космос; сложную конструкцию системы транспортных тепловых труб; жидкостно-контурную передачу тепла с южной на северную панель, ДФК с динамической (насосной) системой подачи в контуре до контурной тепловой трубы. Указанный перечень в последовательности обеспечивает рост массоэнергетической характеристики системы СТР [12]. Следует отметить, что габариты северной и южной панелей определяются площадью испарителя и конденсатора СТР, интегрированных в сотопанели, и определяют, в основном, габариты и силовую схему КА в целом. Следующий логический шаг – рекуперация электрической энергии в системе двухфазной СТР на основе органического цикла Ренкина с целью снижения тепловой нагрузки на конденсатор-излучатель. Самым важным узлом такой системы является турбопривод, состоящий из подводящего устройства, обеспечивающего поле скоростей и давлений на входе в последующее рабочее колесо турбины. Для турбин малой быстроходности с малорасходным парциальным подводом характерна значительная несимметричность полей термодинамических параметров, приводящая к значительным отклонениям в результатах расчета по методикам полноразмерных агрегатов. В соответствии с изложенным, следуя поставленной цели повышения массоэнергетических характеристик двухфазных СТР за счет применения турбогенератора в цикле Ренкина, необходимо решить следующие задачи:
– разработать и изготовить комплект типоразмеров подводящих устройств с целью получения энергетических характеристик в диапазоне изменения конструктивных и режимных параметров;
– провести экспериментальные исследования энергетических характеристик подводящих устройств с оформлением базы данных по основным конструктивным и режимным параметрам.
Методика проведения энергетических и частичных балансовых испытаний микротурбин
Энергетические испытания проводятся с элементами балансовых испытаний, возможных на материальной части физического макета с дополнительной установкой измерительных (станций) постов: p* – полного давления; р – статического; Tизм – измеренной равновесной температуры.
Расположение постов (измерительных станций), с целью коррекции расчетного алгоритма, полностью совпадает с принятой расчетной схемой (рис. 3).
Рис. 3. Расчетная схема радиальной ступени
Fig. 3. Design scheme of the radial stage
Последовательности постов измерений соответствует:
– пост измерений перед входом в канальном подводе , , ;
– в канале подводящего устройства p0, , Т0изм;
– на выходе из канального подводящего устройства, на входе в рабочее колесо p1, , Т1изм,
где – измеряется в окружном направлении;
– на выходе из рабочего колеса на радиусе R2, p2, , Т1изм.
Поскольку технически сложно разместить приемник полного давления на радиусе R1, продувка канальных подводящих устройств проводится в специальном приспособлении, без рабочего колеса, на радиусе входа в рабочее колесо. При необходимости эти испытания согласуются по p1 – статическому давлению при энергетических испытаниях.
Методика обработки результатов на участке /вх – о/.
Измеряемое , , p0, , Т0изм, М0, Т0, τf 0, С0.
Дополнительно рассчитывает массовый расход двумя способами для исключения промаха.
, (1)
где
, (2)
. (3)
F0 = h0 ∙ b0 – площадь проходного сечения канала подводящего устройств.
Коэффициент потерь ζвх рассчитывается:
. (4)
Методика обработки результатов на участке /0 – 1/.
Измеренные: , p1, Т1изм.
Известны: p0, , Т0изм.
Рассчитываются по выражениям для окружной составляющей С1u:
M1u, T1u, τf 1, С1u, pu, . (5)
Параметры рассчитаны по окружной составляющей С1u, необходима коррекция:
, (6)
где C1R – радиальная составляющая абсолютной скорости на входе – определяется в первом приближении и
, (7)
где ṁ (1).
Корректируется полное давление:
. (8)
И совершается прогон пока:
При известных С1, С1u, С1R рассчитывается угол абсолютной скорости:
. (9)
Угол относительной скорости:
, (10)
.
Угол атаки:
, (11)
где β1л – конструктивный параметр.
Коэффициент потерь сопла (подводящего аппарата):
. (12)
Описание и конструктивные особенности исследуемых канальных подводящих устройств
Экспериментальная продувка канальных подводящих устройств производилась в специальном приспособлении, имитирующем энергетическую установку и позволяющем производить измерения по шести измерительным постам, согласно схеме (рис. 4).
Рис. 4. Схема расположения датчиков в испытательной системе
Fig. 4. The layout of the sensors in the test system
Измерительный пост на входе pвх, Tвх, измерительный пост в канале подводящего устройства p0 T0 и четырех измерительных постах на радиусе рабочего колеса p1T1, p2T2, p3T3, p4T4. На посту p2T2 измерения не производились, пост резервный. На каждом посту измерялись полное давление (p*) приемником полного давления и статическое давление (pст) приемником статического давления, равновесная измеренная температура (t) – термопарой.
Экспериментальная установка с размещенными приемниками полного и статического давления показана на рис. 5.
Рис. 5. Экспериментальная установка с приемниками полного и статического давления
Fig. 5. Experimental installation with full and static pressure receivers
Для геометрических параметров используются следующие определения и расчетные соотношения:
– lвх – плечо момента количества движения, где:
;
– dг экв – диаметр горла эквивалентный, где
;
– ε – степень парциальности, где
;
Re0 – число Рейнольдса по скорости С0, где
Reω – число Рейнольдса по окружной скорости, где
;
– lвх/R1k – относительное плечо момента количества движения, где R1k = 24 мм – радиус входа в рабочее колесо микротурбины.
Спиральное подводящее устройство формируется совмещением двух радиусов профилирования R1 и R2 со смещением центров профилирования на величину Δ.
Обмеряемые параметры спирального подводящего устройства микротурбин:
b = 3,5 мм – ширина канала;
D – внешний максимальный диаметр;
δ – расстояние от D до верхней (ближней от D) поверхности канала h;
L1 – максимальной расстояние по линии диаметра D;
h – ширина входного канала;
R1 – больший радиус профилирования;
R2 – меньший радиус профилирования;
Δ – смещение центров профилирования;
ε – степень парциальности;
lвх – плечо входного канала по средней линии;
h/R1k – относительная ширина входного канала;
lвх/R1k – относительное плечо входа.
.
Выражения для обсчета спирального поводящего устройства по результатам обмера:
;
,
,
,
,
,
.
Выражения для обсчета кольцевого поводящего устройства:
,
,
или
.
Приняты следующие обозначения для сборников C – спиральный, первая цифра в обозначении: высота h в мм. Вторая цифра в обозначении: lвх плечо момента количества движения в мм (округлено). Например: C6–39 – спиральный h = 6 мм, lвх ≈ 39 мм; К6–32 – кольцевой h = 6 мм, lвх ≈ 32 мм.
Согласно обозначениям, количественные значения параметров 19 спиральных и 13 кольцевых подводящих устройств представлены в табл. 1.
Таблица 1. Диапазон конструктивных параметров спиральных и кольцевых подводящих устройств
Обозначение | h (мм) | lвх (мм) | drэкв (мм) | lвх/R1k | h/R1k |
С2…6–28…39 | 2–6 | 28–38,9 | 3–5,2 | 1,16–1,62 | 0,083–0,25 |
К2...6–25...32 | 2–6 | 25–32 | 3–5,2 | 1,04–1,33 | 0,083–0,25 |
Диапазон изменения геометрических параметров в безразмерном (относительном) виде (22 конструктивные сборки):
Re0 = 10000–60000 = 1 – 6∙104 – число Рейнольдса по скорости С0;
lвх/R1k = 0,97–1,62 – относительное плечо входа;
h/R1k = 0,083–0,25 – относительная ширина входного канала;
b = 3,5 мм – ширина канала;
R1k = 24 мм – радиус входа в рабочее колесо микротурбины.
Методика обработки экспериментальных испытаний тангенциального подвода дозвуковой центростремительной турбины
Для разработки математической модели и алгоритма расчета центростремительной микротурбины, применяющейся для утилизации тепловой мощности КА, необходимо оценивать окружную составляющую абсолютной скорости на радиусе входа в рабочее колесо (РК) турбины. Теоретический анализ суммарных потерь потока как суперпозиция внезапного расширения с последующим вихревым течением до рабочего колеса представляет теоретически неопределенную задачу. Для предварительной оценки и формирования базы данных потерь на этом участке необходимо использовать данные экспериментальных продувок с регистрацией энергетических и скоростных параметров потока на участке: 0 – параметры во входном канале устройства; 1 – параметры на входном диаметре во входном канале устройства. Для обработки экспериментальных данных помимо коэффициента потерь энергии ζ01 удобно воспользоваться коэффициентом окружной составляющей скорости φu [13–15].
Обоснуем некоторые соображения, определяющие содержание коэффициента окружной составляющей скорости.
Определим массовый расход в подводящем канале:
, (13)
где ρ – плотность; C0 – расходная скорость в канале; h0 – высота канала; b0 – ширина.
Допустим, что C0 – const по высоте h0, плечо момента количества движения по h0 изменяется от l1 до l2, тогда момент количества движения во входном сечением запишется как интеграл:
. (14)
Подставим пределы интегрирования l2 R0; l1 = R0 – h0 и учтем (13):
,
. (15)
Теоретический момент количества движения на текущем радиусе R:
, (16)
где ṁ – массовый расход определяется выражением:
. (17)
Для идеального случая, когда нет момента сил трения, моменты M0 и М равны:
M0 = МТ,
Подставив выражения для моментов (16), получаем:
. (18)
Продолжим преобразование, выразим значения окружной составляющей на текущем радиусе через параметры во входном канале:
, (19)
или с учетом
,
. (20)
Для радиуса входа в рабочее колесо R1 при отсутствии потерь
. (21)
Коэффициент окружной составляющей скорости определяется
, (22)
где – действительное (измеренное) значение; U1т – максимально теоретически возможное.
Действительное значение окружной скорости рассчитывается как среднее по результатам измерений полного и статического давлений и температуры по периферии рабочего колеса микротурбины в 4 точках окружности R1 (рис. 4).
Общий вид сборки представлен на рис. 6, где показана подводка к термисторам, приемники давления не подключены к измерительным трубкам датчиков давлений.
Рис. 6. Общий вид сборки
Fig. 6. General view of the assembly
Для предварительной оценки окружной скорости на входе в РК микротурбины необходимо использовать уравнения энергии по периферии РК:
, (23)
где ξc – коэффициент потерь полного напора в сопле; абсолютная скорость на входе в РК:
, (24)
расходная составляющая определяется только расчетно:
, (25)
где ṁ – массовый расход рабочего тела; ρ1 – плотность на входе в РК; F1 – площадь проходного сечения в РК.
Коэффициент скорости представляет собой параметр, рассчитанный непосредственно по результатам измерений без дополнительных допущений. По результатам эксперимента в диапазоне изменений геометрических и режимных параметров подводящего устройства микротурбины (h0 – высота канала; lвх – среднее плечо момента количества движения; R1 – радиус входа в рабочее колесо; Re – число Рейнольдса на входе) реализуется возможность формирования базы данных для φu с целью использования при математическом моделировании и проектировании микротурбины.
Использование баз данных по коэффициенту окружной составляющей скорости φu (22) и коэффициенту потерь полного давления ζс (23) позволяет однозначно определить полное давление на входе в РК микротурбины при математическом моделировании и проектировании, а также проводить оптимизацию в диапазоне изменений режимных и геометрических параметров микротурбины.
Алгоритм обработки экспериментальных данных испытаний тангенциальных, спиральных и кольцевых подводов реализован в виде программного алгоритма, позволяющего по замерам в 0 (канальном подводе) и 1 сечении (по периферии рабочего колеса), значений давлений, температур и расчетных скоростей определять энергетические параметры подводящего устройства: ζс – коэффициент потерь (12) и φu – коэффициент скорости (21).
Согласно конструктивным схемам материальной части по спиральным и кольцевым подводящим устройствам, проведены экспериментальные замеры коэффициента потерь ζс0 и φu в диапазоне изменения геометрии подводящих устройств (табл. 1, 2) и режимного параметра Re0, рассчитанного по значению скорости C0 (U0) в тангенциальном подводе. Данные представлены в виде базы данных в табл. 2, 3. В общем случае база данных представляет собой табличную функцию коэффициента потерь ; коэффициента скорости ; от трех переменных:
Lвх/R1k – относительное плечо входа;
h/R1k – относительная ширина входного канала;
Re0 – число Рейнольдса по скорости С0(U0) в тангенциальном канале.
Таблица 2. Коэффициенты потерь ζС0 и скорости φu для кольцевого подводящего устройства
К2-25 | К6-32 | ||||||||
ζ С0 | Lвх/R1k | h/R1k | Re0 | φu | ζ С0 | Lвх/R1k | h/R1k | Re0 | φu |
3,0209067 | 1,04 | 0,083 | 17520,65 | 0,2889737 | 2,6914906 | 1,33 | 0,25 | 28199,46 | 0,350395 |
2,81332385 | 1,04 | 0,083 | 21327,45 | 0,4220766 | 1,3009716 | 1,33 | 0,25 | 24833,99 | 0,5772179 |
2,7729739 | 1,04 | 0,083 | 26275,56 | 0,4445341 | 1,2320646 | 1,33 | 0,25 | 31322,3 | 0,5997001 |
3,1156121 | 1,04 | 0,083 | 29045,19 | 0,3552896 | 1,732688 | 1,33 | 0,25 | 39639,36 | 0,5403541 |
2,8642894 | 1,04 | 0,083 | 30264,05 | 0,4435942 | 2,0409651 | 1,33 | 0,25 | 44025,96 | 0,5012578 |
2,880883 | 1,04 | 0,083 | 32048,86 | 0,4319966 | 1,7411898 | 1,33 | 0,25 | 47618,09 | 0,5481886 |
2,7712986 | 1,04 | 0,083 | 32068,46 | 0,4661496 | 1,8086491 | 1,33 | 0,25 | 50644,91 | 0,541099 |
2,8155358 | 1,04 | 0,083 | 30754,69 | 0,4596947 | 1,7884947 | 1,33 | 0,25 | 50892,56 | 0,5436699 |
Таблица 3. Коэффициенты потерь ζС0 и скорости φu для спирального подводящего устройства
С2-28 | С6-39 | ||||||||
ζ С0 | Lвх/R1k | h/R1k | Re0 | φu | ζ С0 | Lвх/R1k | h/R1k | Re0 | φu |
3,3439123 | 1,16 | 0,083 | 30933,17 | 0,2453211 | 1,7325754 | 1,62 | 0,25 | 19870,44 | 0,4304552 |
3,01471724 | 1,16 | 0,083 | 38247,2 | 0,3147172 | 2,6571707 | 1,62 | 0,25 | 30967,94 | 0,3141631 |
3,1387546 | 1,16 | 0,083 | 45548,23 | 0,3981814 | 1,9631256 | 1,62 | 0,25 | 32378,77 | 0,4201886 |
2,98048791 | 1,16 | 0,083 | 55432,55 | 0,3672183 | 2,1385318 | 1,62 | 0,25 | 40561,44 | 0,4038012 |
2,9605327 | 1,16 | 0,083 | 60051,83 | 0,3609476 | 2,0742571 | 1,62 | 0,25 | 44835,47 | 0,4082927 |
2,9455596 | 1,16 | 0,083 | 64664,08 | 0,3731699 | 1,8497709 | 1,62 | 0,25 | 52284,86 | 0,4391877 |
3,12578929 | 1,16 | 0,083 | 67022,12 | 0,3387341 | 1,801266 | 1,62 | 0,25 | 54236,36 | 0,445411 |
3,3556589 | 1,16 | 0,083 | 69246,83 | 0,2580974 | 1,6874583 | 1,62 | 0,25 | 58712,84 | 0,4565281 |
По табличным данным визуализация поверхности уровня функция достаточно монотонна, не имеет резких экстремумов и достаточно легко апроксимируется даже линейными сплайнами (рис. 7). Значение функции в 2–3 раза превышает коэффициент потерь на внезапное расширение. Очевидно, в структуре энергетических потерь должен присутствовать дополнительный элемент. В качестве такого элемента следует рассмотреть потери энергии потока при течении потока от радиуса R0 до радиуса R1k, что следует из конструктивной геометрии подводящего устройства. Для оценки коэффициента потерь на радиально окружном участке течения воспользуемся численной оценкой по данным испытания кольцевого подводящего устройства К2-25. В качестве исходных данных для расчета используем значения 0,02 м; 0,002 м; 0,00017 м; 2,314 [м2/с]. Коэффициенты потерь рассчитываются по выражению, полученному из уравнения энергии:
,
откуда ζрад.уч:
. (26)
Рис. 7. Поверхности уровней функции (базы данных). : а – относительная ширина входного канала ; б –
Fig. 7. Surfaces of function levels (databases). : a – relative width of the input channel ; b –
Следует отметить, что коэффициент потерь на участке радиального течения формируется за счет падения статического давления и окружной составляющей скорости. В одномерном потенциальном течении изменение полного и статического давлений совпадает, скорость не изменяется. Результаты расчетов представлены в табл. 4 [16–17]. Из результатов видно, что коэффициент потерь при изменении Рейнольдса изменяется практически в два раза и близок к единице. Коэффициент потерь формируется из двух приблизительно одинаковых членов: потерь статического давления, потерь окружной скорости.
Таблица 4. Результаты расчетов коэффициентов потерь на радиальном участке течения
№ | С0i(U0) [м/с] | [Па] | [Па] | τ0α [Н/м2] | dCu/dR [м/с] | Cu0 [м2/с] | ΔCu0 [м2/с] | Cuk [м2/с] | Uk [м/с] | Δp [Па] | pk [Па] | ζ рад.уч. | Re0 | ||
1 | 89 | 103750 | 98955 | 31,4 | 408,2 | 2,314 | 0,816 | 1,5 | 62,4 | 662 | 98923 | 0,487 | 0,508 | 0,995 | 17500 |
2 | 107 | 111260 | 103710 | 43,32 | 458,4 | 2,782 | 0,916 | 1,87 | 77,8 | 961 | 102749 | 0,518 | 0,471 | 0,983 | 21300 |
3 | 133 | 123277 | 110844 | 63,4 | 544,3 | 3,458 | 1,09 | 2,37 | 98,7 | 1493 | 109351 | 0,492 | 0,449 | 0,941 | 26286 |
4 | 147 | 139803 | 122733 | 75,5 | 582,2 | 3,828 | 1,16 | 2,662 | 111 | 1828 | 120905 | 0,493 | 0,43 | 0,923 | 29045 |
5 | 153 | 157831 | 13700 | 81,0 | 594,4 | 3,98 | 1,19 | 2,8 | 116 | 1982 | 135018 | 0,4936 | 0,425 | 0,919 | 30264 |
6 | 162 | 174356 | 148890 | 89,23 | 610,5 | 4,2 | 1,221 | 2,98 | 124 | 2216 | 146674 | 0,494 | 0,412 | 0,906 | 32048 |
7 | 162 | 193866 | 165535 | 89,4 | 598,2 | 4,21 | 1,196 | 2,014 | 125,6 | 2223 | 163312 | 0,4944 | 0,398 | 0,8924 | 32068 |
С учетом данных коэффициента на внезапное расширение и коэффициента потерь, на радиальном участке остается приблизительно 1/3 от общего коэффициента потерь, т. е. величина приблизительно равная единице. В качестве гипотезы предлагается эту величину оценить коэффициентом потерь на разворот потока на 90º по формуле Вейсбаха [14]:
,
где δ = 90º – угол разворота, с учетом поправки коэффициент потерь на разворот ; C1 = 1 для симметричного течения; ; тогда коэффициент разворота на = 1,3. Следует обратить внимание, что предложенная энергетическая балансировка потерь на внезапное расширение, радиальный участок и разворот потока на 90º дает величину, размещенную в базе данных по коэффициенту потерь ζС0 в подводящем устройстве (табл. 2, 3).
Заключение
По результатам исследований видно, что подводящее тангенциальное устройство радиальной микротурбины топологически представляет собой граничные условия течения, сочетающие в себе последовательно гидравлические потери местных сопротивлений: внезапное расширение, разворот потока на 90º, поток на радиально-закрученном участке перед рабочим колесом. Формально гидравлически местные потери суммируются через гидравлические прямолинейные участки. В нашем случае потери на местное сопротивление объединены и взаимно влияют друг на друга, что исключает их суммирование в общем случае [14]. Поэтому с точки зрения адекватности результатов предпочтительно пользоваться базой данных по коэффициенту потерь в подводящем устройстве. Сложение потерь возможно только как приближенная оценка.
- Экспериментальные результаты по коэффициенту потерь и для кольцевого и спирального тангенциальных подводящих устройств представлены в виде базы данных. В общем случае база данных представляет собой табличную функцию от трех переменных:
Lвх/R1k – относительное плечо входа;
h/R1k – относительная ширина входного канала;
Re0 – число Рейнольдса по скорости С0(U0) в тангенциальном канале.
- Функция коэффициента потерь достаточно монотонна и не имеет ярких экстремумов, значение функции приблизительно в три раза превышает значение при внезапном расширении потока, что предполагает ее аддитивную структуру.
- По результатам сравнительного энергетического анализа в первом приближении предложена следующая структура энергетических потерь в тангенциальном подводящем устройстве микротурбины, практически согласующаяся с величиной экспериментальных потерь: , где – потери на внезапное расширение; – потери на разворот потока на 90º; – потери на трение на радиальном участке канала.
- Аналогичная база данных получена для коэффициента скорости φu.
Об авторах
Юлия Николаевна Шевченко
Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева
Автор, ответственный за переписку.
Email: gift_23j@mail.ru
аспирант
Россия, КрасноярскАлександр Анатольевич Кишкин
Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева
Email: spsp99@mail.ru
доктор технических наук, профессор, заведующий кафедрой холодильной, криогенной техники и кондиционирования
Россия, КрасноярскАлександр Александрович Зуев
Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева
Email: dla2011@inbox.ru
доктор технических наук, профессор, заведующий кафедрой двигателей летательных аппаратов
Россия, КрасноярскАлександр Викторович Делков
Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева
Email: delkov-mx01@mail.ru
кандидат технических наук, доцент кафедры холодильной, криогенной техники и кондиционирования
Россия, КрасноярскСписок литературы
- Органический цикл Ренкина в автономной теплоэнергетической системе : монография / А. А. Кишкин, О. В. Шилкин, А. В. Делков и др. ; СибГУ им. М. Ф. Решетнева. Красноярск, 2019. 234 с.
- Разработка установок-утилизаторов низкопотенциального тепла на основе органического цикла Ренкина / А. А. Кишкин, Д. В. Черненко, А. А. Ходенков и др. // Альтернативная энергетика и экология. 2014. № 3 (4). С. 35–36.
- Расчет и анализ тепловых технических систем, работающих по замкнутому циклу / А. А. Кишкин, Е. В. Черненко, Д. В. Черненко и др. // Materiály VIII mezinárodní vědecko – praktická conference “Dny vědy – 2012”. Díl 91. Technické vědy: Praha. Publishing House Education and Science s.r.o, 2012.
- Боровский Б. И. Энергетические параметры и характеристики высокооборотных лопастных насосов. М. : Машиностроение, 1989. 184 с.
- Теория пространственного пограничного слоя в гидродинамике турбомашин : монография / А. А. Кишкин, В. П. Назаров, Д. А. Жуйков, Д. В. Черненко ; Сиб. гос. аэрокосмич. ун-т. Красноярск, 2013. 250 с.
- Кишкин А. А., Зуев А. А., Делков А. В. Трехмерный температурный пограничный слой в теории конвективного теплообмена : монография / Сиб. гос. аэрокосмич. ун-т. Красноярск, 2015. 282 с.
- Смирнов М. В. Безлопаточные центробежные ступени для турбодетандоров малой мощности : дис. … канд. техн. наук. СПб. 2019. 154 с.
- Теплосиловые системы: Оптимизационные исследования / А. М. Клер, Н. П. Деканова, Э. А. Тюрина и др. Новосибирск : Наука, 2005. 236 с.
- Краев М. В., Лукин В. А., Овсянников Б. В. Малорасходные насосы авиациаонных и космических систем. М. : Машиностроение, 1985. 128 с.
- Краев М. В., Кишкин А. А., Сизых Д. Н. Гидродинамика малорасходных насосных агрегатов. Красноярск : САА, 1988. 157 с.
- Прототипирование микротурбогенератора и постановка задачи исследования / Ю. Н. Шевченко, О. В. Шилкин, А. А. Кишкин и др. // Испытания, диагностика, надежность. Теория и практика : сб. тр. Всеросс. науч.-практ. конф. Красноярск, 2023. С. 17–21.
- Моделирование и конструирование двухфазных систем термо- регулирования космических аппаратов : монография / О. В. Шилкин, А. А. Кишкин, А. В. Делков, и др. ; СибГУ им. М. Ф. Решетнева. Красноярск, 2022. 192 с.
- Коэффициент скорости тангенциального подвода дозвуковой центростремительной турбины / Ю. Н. Шевченко, А. А. Кишкин, А. В. Делков, М. У. Абдуллаев // Омский науч. вестник. Серия: Авиационно-ракетное и энергетическое машиностроение. 2022. Т. 6, № 2. С. 78–84.
- Идельчик И. Е. Справочник по гидравлическим сопротивлениям / под ред. М. О. Штейнберга. М. : Машиностроение, 1992. 672 с.
- Kishkin A. A. , Shevchenko Yu. N., Delkov A. V. Analysis of the key design features of low-power turbines for electricity generation // IOP Conference Series: Materials Science and Engineering. 2020. Vol. 919. doi: 10.1088/1757-899X/919/6/062030.
- Energy equations for the temperature three-dimensional boundary layer for the flow within boundary conditions of turbo machinery / A. A. Zuev, A. A. Kishkin, D. A. Zhuikov et al. // IOP Conference Series: Materials Science and Engineering. 2019. Vol. 537. P. 22008. doi: 10.1088/1757-899X/537/2/022008.
- Tn. Karman. Uber laminare und turbulente Reibung // ZAAM. 1921. No. 1. P. 233–252.
Дополнительные файлы
