Энергетика тангенциального подводящего устройства микротурбины системы терморегулирования перспективного космического аппарата

Обложка

Цитировать

Полный текст

Аннотация

В настоящей работе представлен обзор современной технической проблемы, связанной с двухфазными системами терморегулирования космических аппаратов, и возможные технические приложения рекуперации тепловой энергии в органическом цикле Ренкина как составной части систем обеспечения теплового режима. Конструктивное решение подразумевает собой интегрирование паровой микротурбины за радиатором-испарителем. Микротурбина представляет собой тангенциальное подводящее устройство и радиально центростремительное рабочее колесо низкой быстроходности 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

 

Последовательности постов измерений соответствует:

– пост измерений перед входом в канальном подводе p2*, pвх*, Твх*=Твх;

– в канале подводящего устройства p0, p0*, Т0изм;

– на выходе из канального подводящего устройства, на входе в рабочее колесо p1, p1u*, Т1изм,

где p1u* – измеряется в окружном направлении;

– на выходе из рабочего колеса на радиусе R2, p2, p2u*, Т1изм.

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

Методика обработки результатов на участке /вх – о/.

Измеряемое pвх*, Твх*, p0, p0*Т0изм, М0, Т0, τf 0, С0.

Дополнительно рассчитывает массовый расход двумя способами для исключения промаха.

m·=ρ0C0F0=ρ0F0mqf0, (1)

где

m=k2k+1k+1k11R, (2)

qf0=k+12k+12k1M01+k12M0k+12k1. (3)

F0 = h0 b0 – площадь проходного сечения канала подводящего устройств.

Коэффициент потерь ζвх рассчитывается:

2С0kk1pвх*ρвхp0*ρ0=ξвх. (4)

Методика обработки результатов на участке /0 – 1/.

Измеренные: p1u*, p1, Т1изм.

Известны: p0, p0*, Т0изм.

Рассчитываются по выражениям для окружной составляющей С1u:

M1u, T1u, τf 1, С1u, pu, Т1u*. (5)

Параметры рассчитаны по окружной составляющей С1u, необходима коррекция:

С1=C1u2+C1R2, (6)

где C1R – радиальная составляющая абсолютной скорости на входе – определяется в первом приближении p1=p1u и

С1R=m˙ρ1F2=m˙ρ12πR1b1, (7)

где (1).

Корректируется полное давление:

p1*=p1+C122k1k. (8)

И совершается прогон пока:

С1i+1C1iС1i+10,01

При известных С1, С1u, С1R рассчитывается угол абсолютной скорости:

α1=arctgC1RC1u. (9)

Угол относительной скорости:

β1=arctgC1RC1uU1=arctgC1RC1uω1R1, (10)

W1=C1uU1cosβ.

Угол атаки:

i=β1лβ1, (11)

где β – конструктивный параметр.

Коэффициент потерь сопла (подводящего аппарата):

ζС0=p0*ρ0p1*ρ12C02kk1. (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вх – плечо момента количества движения, где:

lвх=Rвхh/2;

dг экв – диаметр горла эквивалентный, где

dг.экв=4bhπ;

– ε – степень парциальности, где

ε=αвх2π, αвх=arccosRвхhRвх;

Re0 – число Рейнольдса по скорости С0, где

Re0=C0dг.эквv

Reω – число Рейнольдса по окружной скорости, где

Reω=C0lвхν=ωlвх2ν;

lвх/R1k – относительное плечо момента количества движения, где R1k = 24 мм – радиус входа в рабочее колесо микротурбины.

Спиральное подводящее устройство формируется совмещением двух радиусов профилирования R1 и R2 со смещением центров профилирования на величину Δ.

Обмеряемые параметры спирального подводящего устройства микротурбин:

b = 3,5 мм – ширина канала;

D – внешний максимальный диаметр;

δ – расстояние от D до верхней (ближней от D) поверхности канала h;

L1 – максимальной расстояние по линии диаметра D;

h – ширина входного канала;

R1 – больший радиус профилирования;

R2 – меньший радиус профилирования;

Δ – смещение центров профилирования;

ε – степень парциальности;

lвх – плечо входного канала по средней линии;

h/R1k – относительная ширина входного канала;

lвх/R1k – относительное плечо входа.

dг.экв=4bhπ.

Выражения для обсчета спирального поводящего устройства по результатам обмера:

L1=R1+R2+;

Δ=0,5h+δ+0,5L10,5D,

R1=D0,5h0,5L1,

R2=0,5L1h,

lвх=0,5Dδ0,5h,

αвх=arccoslвх0,5hlвх+0,5h,

ε=αвх.

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

R1=L12,

lвх=0,5L10,5h,

или

lвх=R1h2.

Приняты следующие обозначения для сборников 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/R1= 0,083–0,25 – относительная ширина входного канала;

b = 3,5 мм – ширина канала;

R1k  = 24 мм – радиус входа в рабочее колесо микротурбины.

Методика обработки экспериментальных испытаний тангенциального подвода дозвуковой центростремительной турбины

Для разработки математической модели и алгоритма расчета центростремительной микротурбины, применяющейся для утилизации тепловой мощности КА, необходимо оценивать окружную составляющую абсолютной скорости на радиусе входа в рабочее колесо (РК) турбины. Теоретический анализ суммарных потерь потока как суперпозиция внезапного расширения с последующим вихревым течением до рабочего колеса представляет теоретически неопределенную задачу. Для предварительной оценки и формирования базы данных потерь на этом участке необходимо использовать данные экспериментальных продувок с регистрацией энергетических и скоростных параметров потока на участке: 0 – параметры во входном канале устройства; 1 – параметры на входном диаметре во входном канале устройства. Для обработки экспериментальных данных помимо коэффициента потерь энергии ζ01 удобно воспользоваться коэффициентом окружной составляющей скорости φu [13–15].

Обоснуем некоторые соображения, определяющие содержание коэффициента окружной составляющей скорости.

Определим массовый расход в подводящем канале:

m˙=ρV˙=С0h0b0, (13)

где ρ – плотность; C0 – расходная скорость в канале; h0 – высота канала; b0 – ширина.

Допустим, что C0 – const по высоте h0, плечо момента количества движения по h0 изменяется от l1 до l2, тогда момент количества движения во входном сечением запишется как интеграл:

M0=l1l2ρC02lb0dl=ρC0h0b0C0l22l1l2. (14)

Подставим пределы интегрирования l2 R0; l1 = R0 h0 и учтем (13):

M0=m˙C012l22l12,

M0=12m˙C02Rh00h02=m˙C0h0R0h02. (15)

Теоретический момент количества движения на текущем радиусе R:

MТ=m˙UТR, (16)

где – массовый расход определяется выражением:

m˙=2πRbCRρ. (17)

Для идеального случая, когда нет момента сил трения, моменты M0 и М равны:

M0 = МТ,

Подставив выражения для моментов (16), получаем:

m˙C0h0R0h02=m˙UТR. (18)

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

UT=С0h0RR012h0, (19)

или с учетом

R0=lвх+h02,

UT=С0h0Rlвх. (20)

Для радиуса входа в рабочее колесо R1  при отсутствии потерь

U1T=С0h0R1lвх. (21)

Коэффициент окружной составляющей скорости определяется

φu=U1U1T, (22)

где U1 – действительное (измеренное) значение; U – максимально теоретически возможное.

Действительное значение окружной скорости рассчитывается как среднее по результатам измерений полного и статического давлений и температуры по периферии рабочего колеса микротурбины в 4 точках окружности R1 (рис. 4).

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

 

Рис. 6. Общий вид сборки

Fig. 6. General view of the assembly

 

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

kk1p0ρ0+C022=kk1p1ρ1+C122+ζcC022, (23)

где ξc – коэффициент потерь полного напора в сопле; абсолютная скорость на входе в РК:

С12=U12+СR12, (24)

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

СR1=m˙ρ1F1, (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. В общем случае база данных представляет собой табличную функцию коэффициента потерь ζс0=Lвх/R1k,h/R1k,Re0; коэффициента скорости φu=Lвх/R1k,h/R1k,Re0; от трех переменных:

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. В качестве исходных данных для расчета используем значения R0=lвх+h/2=0,02 м; ΔR=R0Rk=0,002 м; δα**=0,00017 м; Сu0=U0R0=2,314 [м2/с]. Коэффициенты потерь рассчитываются по выражению, полученному из уравнения энергии:

kk1p0ρ0+U022=kk1p1ρ1+U122+ζрад.уч.U022,

откуда ζрад.уч:

ζрад.уч=2U02kk11ρp0pk+U02Uk22. (26)

 

Рис. 7. Поверхности уровней функции (базы данных). ζс0=(Lвх/R1k,h/R1k,Re0) : а – относительная ширина входного канала h/R1k=0,092; б – h/R1k=0,013

Fig. 7. Surfaces of function levels (databases). ζс0=(Lвх/R1k,h/R1k,Re0) : a – relative width of the input channel h/R1k=0,092; b – h/R1k=0,013

 

Следует отметить, что коэффициент потерь на участке радиального течения формируется за счет падения статического давления и окружной составляющей скорости. В одномерном потенциальном течении изменение полного и статического давлений совпадает, скорость не изменяется. Результаты расчетов представлены в табл. 4 [16–17]. Из результатов видно, что коэффициент потерь при изменении Рейнольдса изменяется практически в два раза и близок к единице. Коэффициент потерь формируется из двух приблизительно одинаковых членов: потерь статического давления, потерь окружной скорости.

 

Таблица 4. Результаты расчетов коэффициентов потерь на радиальном участке течения

С0i(U0) [м/с]

[Па]

[Па]

τ[Н/м2]

dCu/dR

[м/с]

Cu0

2/с]

ΔCu0

2/с]

Cuk

2/с]

Uk

[м/с]

Δp

[Па]

pk

[Па]

2kk1ρp0pkU021UkU02

ζ рад.уч.

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]:

ζм=0,95sin2δ2+2,05sin4δ2=0,986,

где δ = 90º – угол разворота, с учетом поправки коэффициент потерь на разворот ζпов=С1AζмC1 = 1 для симметричного течения; А=0,95+33,5/δ=1,32; тогда коэффициент разворота на  ζ90= 1,3. Следует обратить внимание, что предложенная энергетическая балансировка потерь на внезапное расширение, радиальный участок и разворот потока на 90º дает величину, размещенную в базе данных по коэффициенту потерь ζС0 в подводящем устройстве (табл. 2, 3).

Заключение

По результатам исследований видно, что подводящее тангенциальное устройство радиальной микротурбины топологически представляет собой граничные условия течения, сочетающие в себе последовательно гидравлические потери местных сопротивлений: внезапное расширение, разворот потока на 90º, поток на радиально-закрученном участке перед рабочим колесом. Формально гидравлически местные потери суммируются через гидравлические прямолинейные участки. В нашем случае потери на местное сопротивление объединены и взаимно влияют друг на друга, что исключает их суммирование в общем случае [14]. Поэтому с точки зрения адекватности результатов предпочтительно пользоваться базой данных по коэффициенту потерь в подводящем устройстве. Сложение потерь возможно только как приближенная оценка.

  1. Экспериментальные результаты по коэффициенту потерь и для кольцевого и спирального тангенциальных подводящих устройств представлены в виде базы данных. В общем случае база данных представляет собой табличную функцию от трех переменных:

Lвх/R1k – относительное плечо входа;

h/R1k – относительная ширина входного канала;

Re0 – число Рейнольдса по скорости С0(U0) в тангенциальном канале.

  1. Функция коэффициента потерь достаточно монотонна и не имеет ярких экстремумов, значение функции приблизительно в три раза превышает значение при внезапном расширении потока, что предполагает ее аддитивную структуру.
  2. По результатам сравнительного энергетического анализа в первом приближении предложена следующая структура энергетических потерь в тангенциальном подводящем устройстве микротурбины, практически согласующаяся с величиной экспериментальных потерь: ζс0=ζвн.расш.+ζ90+ζрад.уч., где ζвн.расш. – потери на внезапное расширение; ζ90– потери на разворот потока на 90º; ζрад.уч. – потери на трение на радиальном участке канала.
  3. Аналогичная база данных получена для коэффициента скорости φu.
×

Об авторах

Юлия Николаевна Шевченко

Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева

Автор, ответственный за переписку.
Email: gift_23j@mail.ru

аспирант

Россия, Красноярск

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

Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева

Email: spsp99@mail.ru

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

Россия, Красноярск

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

Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева

Email: dla2011@inbox.ru

доктор технических наук, профессор, заведующий кафедрой двигателей летательных аппаратов

Россия, Красноярск

Александр Викторович Делков

Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева

Email: delkov-mx01@mail.ru

кандидат технических наук, доцент кафедры холодильной, криогенной техники и кондиционирования

Россия, Красноярск

Список литературы

  1. Органический цикл Ренкина в автономной теплоэнергетической системе : монография / А. А. Кишкин, О. В. Шилкин, А. В. Делков и др. ; СибГУ им. М. Ф. Решетнева. Красноярск, 2019. 234 с.
  2. Разработка установок-утилизаторов низкопотенциального тепла на основе органического цикла Ренкина / А. А. Кишкин, Д. В. Черненко, А. А. Ходенков и др. // Альтернативная энергетика и экология. 2014. № 3 (4). С. 35–36.
  3. Расчет и анализ тепловых технических систем, работающих по замкнутому циклу / А. А. Кишкин, Е. В. Черненко, Д. В. Черненко и др. // 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.
  4. Боровский Б. И. Энергетические параметры и характеристики высокооборотных лопастных насосов. М. : Машиностроение, 1989. 184 с.
  5. Теория пространственного пограничного слоя в гидродинамике турбомашин : монография / А. А. Кишкин, В. П. Назаров, Д. А. Жуйков, Д. В. Черненко ; Сиб. гос. аэрокосмич. ун-т. Красноярск, 2013. 250 с.
  6. Кишкин А. А., Зуев А. А., Делков А. В. Трехмерный температурный пограничный слой в теории конвективного теплообмена : монография / Сиб. гос. аэрокосмич. ун-т. Красноярск, 2015. 282 с.
  7. Смирнов М. В. Безлопаточные центробежные ступени для турбодетандоров малой мощности : дис. … канд. техн. наук. СПб. 2019. 154 с.
  8. Теплосиловые системы: Оптимизационные исследования / А. М. Клер, Н. П. Деканова, Э. А. Тюрина и др. Новосибирск : Наука, 2005. 236 с.
  9. Краев М. В., Лукин В. А., Овсянников Б. В. Малорасходные насосы авиациаонных и космических систем. М. : Машиностроение, 1985. 128 с.
  10. Краев М. В., Кишкин А. А., Сизых Д. Н. Гидродинамика малорасходных насосных агрегатов. Красноярск : САА, 1988. 157 с.
  11. Прототипирование микротурбогенератора и постановка задачи исследования / Ю. Н. Шевченко, О. В. Шилкин, А. А. Кишкин и др. // Испытания, диагностика, надежность. Теория и практика : сб. тр. Всеросс. науч.-практ. конф. Красноярск, 2023. С. 17–21.
  12. Моделирование и конструирование двухфазных систем термо- регулирования космических аппаратов : монография / О. В. Шилкин, А. А. Кишкин, А. В. Делков, и др. ; СибГУ им. М. Ф. Решетнева. Красноярск, 2022. 192 с.
  13. Коэффициент скорости тангенциального подвода дозвуковой центростремительной турбины / Ю. Н. Шевченко, А. А. Кишкин, А. В. Делков, М. У. Абдуллаев // Омский науч. вестник. Серия: Авиационно-ракетное и энергетическое машиностроение. 2022. Т. 6, № 2. С. 78–84.
  14. Идельчик И. Е. Справочник по гидравлическим сопротивлениям / под ред. М. О. Штейнберга. М. : Машиностроение, 1992. 672 с.
  15. 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.
  16. 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.
  17. Tn. Karman. Uber laminare und turbulente Reibung // ZAAM. 1921. No. 1. P. 233–252.

Дополнительные файлы

Доп. файлы
Действие
1. JATS XML
2. Рис. 1. Прототип микротурбогенератора: 1 – крышка; 2 – корпус микротурбины с тангенциальным подводящим устройством; 3 – рабочее колесо; 4 – генератор

Скачать (583KB)
3. Рис. 2. Вид на рабочее колесо со стороны крышки (а); крышка и корпус генератора (б); вид со стороны крышки генератора (в)

Скачать (373KB)
4. Рис. 3. Расчетная схема радиальной ступени

Скачать (260KB)
5. Рис. 4. Схема расположения датчиков в испытательной системе

Скачать (125KB)
6. Рис. 5. Экспериментальная установка с приемниками полного и статического давления

Скачать (311KB)
7. Рис. 6. Общий вид сборки

Скачать (405KB)
8. Рис. 7. Поверхности уровней функции (базы данных).  : а – относительная ширина входного канала ; б – 

Скачать (253KB)

© Шевченко Ю.Н., Кишкин А.А., Зуев А.А., Делков А.В., 2024

Creative Commons License
Эта статья доступна по лицензии Creative Commons Attribution 4.0 International License.