CC BY
58
МАШИНЫ И АППАРАТЫ, ПРОЦЕССЫ ХОЛОДИЛЬНОЙ И КРИОГЕННОЙ ТЕХНИКИ, СИСТЕМ КОНДИЦИОНИРОВАНИЯ И ЖИЗНЕОБЕСПЕЧЕНИЯ
УДК 621.5
ИССЛЕДОВАНИЕ РАБОТЫ КАНАЛЬНОГО НАПРАВЛЯЮЩЕГО АППАРАТА В СОСТАВЕ РАДИАЛЬНО-ОСЕВОЙ ФРЕОНОВОЙ ТУРБИНЫ
А.А. Костенко1,2, Т.М. Розеноер1,2, В.П. Леонов1
1МГТУ им. Н.Э. Баумана, Москва, Российская Федерация e-mail: kosrk1@gmail.com
2НИИ механики МГУ им. М.В. Ломоносова, Москва, Российская Федерация e-mail: RozenoerTM@imec.msu.ru
Рассмотрены особенности работы канального направляющего аппарата в составе силовой турбины, работающей на фреоне. Приведены результаты численного моделирования течения в каналах проточной части турбомашины для режимов работы с большими степенями расширения в канальном направляющем аппарате. Выполнена оценка влияния шероховатости каналов на основные параметры эффективности работы направляющего аппарата при различных граничных условиях. На основе проведенного исследования сделан вывод о целесообразности применения канального направляющего аппарата для силовой фреоновой турбины.
Ключевые слова: канальный направляющий аппарат, фреоновая турбина, цикл Ренкина, утилизация теплоты, органическое рабочее вещество.
OPERATION ANALYSIS OF CHANNEL GUIDE VANE OF RADIAL-AXIAL FREON TURBINE
А.А. Kostenko1,2, Т.М. Rozenoer1,2, V.P. Leonov1
1Bauman Moscow State Technical University, Moscow, Russian Federation e-mail: kosrk1@gmail.com
2Research Institute of Mechanics of Lomonosov Moscow State University,
Moscow, Russian Federation e-mail: RozenoerTM@imec.msu.ru
Peculiarities of channel guide vane operation within a power freon turbine, are analyzed. The results of numerical modeling offlow in the channels of turbine wheel space are presented for operation modes with large expansion rates. Evaluation of the channel roughness influence on major efficiency parameters of the guide vane operation at various boundary conditions, is performed. Conclusion is drawn of expediency to apply channel guide vane for power turbomachines operating on freon.
Keywords: channel guide vane, freon turbine, Rankine cycle, heat utilization, organic working fluid.
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. № 1 101
Сегодня активно развиваются технологии энергосбережения, разрабатывается большое число установок для утилизации сбросной теплоты и теплоты из возобновляемых источников (солнечная энергия, геотермальная и т.д.). Часть таких установок основана на цикле Рен-кина (Organic Rankine Cycle, ORC), где в качестве рабочего тела применяются различные органические вещества, например фреоны R113, R245fa. Если в составе установки применяется турбомашина, то необходимо обеспечить ее эффективную работу (создать эффективное рабочее колесо, направляющий аппарат (НА) и т.д.). Цель настоящей работы — исследование особенностей применения канального направляющего аппарата в силовой фреоновой турбине.
Канальные направляющие аппараты (КНА) были разработаны в ОАО “Криогенмаш” и применяются во всем размерном ряду турбодетандеров (с диаметрами колес от 30 до 250 мм), их используют в ВРУ средней и большой холодопроизводительности. По сравнению с лопаточными НА они имеют малое число каналов z (z = 6... 10), т.е. увеличенную горловину сопла, позволяют обеспечить малый угол установки сопла ac (от 6°, оптимальные значения ac от 7° до 9°) и большую относительную высоту сопла (меньшие потери трения). Профилировка КНА способствует образованию устойчивого к отрыву пограничного слоя, что приводит к меньшим потерям при закритиче-ских режимах течения, чем в лопаточных НА [1].
В ОАО “Криогенмаш” были проведены исследования эффективности НА такого типа и получены приемлемые значения потерь [1, 2]. Однако все проведенные исследования говорят об успешном использовании КНА в криогенной области, а о возможности их применения в силовых турбинах, работающих на органических веществах, информации нет. Для теплоутилизирующих ORC установок характерен большой перепад давления в расширителе. В случае с турбомашиной весь перепад стараются реализовать в одной детандерной ступени, что приводит к высоким значениям приведенной скорости A1s рабочего вещества на выходе из направляющего аппарата. Для криогенных турбин КНА эффективен при значениях A1s до 1,4, но для фреоновой машины это значение может быть другим.
По сравнению с рабочими веществами низкотемпературных турбодетандеров (например, воздухом) фреоны имеют более высокую молярную массу, плотность, меньшую вязкость и показатель адиабаты к (что обеспечивает большие значения чисел Рейнольдса и благоприятно сказывается на эффективности проточной части), малую скорость изоэнтропийного расширения Cs (что позволяет создавать малоразмерные турбины с частотами вращения, приемлемыми для использования совместно с высокочастотными электрогенераторами). Использование фреонов позволяет создавать эффективные силовые малоразмерные фреоновые турбины.
102 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №1
Задание режимов работы КНА. Было выполнено исследование КНА, работающего в составе фреонового турбоэлектрогенератора малой мощности. Угол установки сопла в данном направляющем аппарате ас = 8°, число каналов z = 6. Канальный направляющий аппарат испытывался при различных режимах работы (разные давление Р0, температура газа Т0 на входе в каналы и давление Р1 - статическое давление перед рабочим колесом). В испытательном стенде был реализован замер Pi (отбор проводился из зазора между колесом и НА). Степень расширения п = Р0/Р2 (где Р2 — давление на выходе из колеса) в эксперименте менялась от 5 до 8,2 для всей ступени и пКНА = Р0/Р1 — от 2,5 до 3,7 для НА. Температура на входе в ступень Т0 менялась от 140 до 160 °С. Замер расхода рабочего вещества на входе в турбину проводился с использованием кориолисового расходомера.
Моделирование течения в проточной части турбины. На основе геометрии КНА и рабочего колеса была создана трехмерная модель канала проточной части для последующего моделирования течения в CFD (Computational fluid dynamics) программном пакете ANSYS CFX методом конечных элементов. По модели построены расчетные сетки (до 900 000 элементов для канала КНА) с высоким разрешением в пристеночной области и вокруг тонких кромок. Течение моделировали как для экспериментальных режимов, так и для режимов с более высокими степенями расширения (п до 13,3, пКНА до 5,3). В качестве граничных условий задавали полное давление Р0 и полную температуру Т0 на входе в КНА, статическое давление Р2 (среднее по площади на некотором расстоянии от выхода из колеса) и частоту вращения ротора. Утечки и перетечки между колесом и покрывным диском не учитывали. Расчет проводился для одного канала КНА (рис. 1, канал заштрихован) и одного канала рабочего колеса, на которые было наложено условие периодичности.
Рис. 1. Проточная часть канального направляющего аппарата
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. № 1 103
В первой серии расчетов канал проточной части рассчитывался как гидравлически гладкий, во второй серии расчетов для стенок канала была задана эквивалентная песочная шероховатость Ks. Поскольку значение песочной шероховатости зависит не только от высоты неровностей, но и от их формы и взаимного расположения [3], то определить его можно путем гидравлического эксперимента. Однако в работе [4] приведена методика приближенной оценки песочной шероховатости через значения технической шероховатости (Ra, Rz), приведены результаты сравнения оценочных данных и результатов экспериментов. Для малых значений технической шероховатости приближенные формулы можно использовать для оценки Ks.
По методике [4] имеем
для Ra = 2,5 мкм по формуле (1) Ks = 15 мкм.
При Ks = 15 мкм величина u*Ks/v > 70 для всех режимов (где и* — динамическая скорость, v — кинематическая вязкость) [3] означает, что в КНА течение происходит с полным проявлением шероховатости (все элементы шероховатости выступают из ламинарного подслоя) с квадратичной зависимостью сопротивления от скорости.
Оценка качества расчета в ANSYS CFX проводилась путем сравнения массового расхода G и осредненного по площади выхода из КНА статического давления Р1, полученных в CFX и в эксперименте. Погрешность расчета G не превышала 7 %, а для Pi — 5 % для всех режимов. В результате численного моделирования течения была получена трехмерная картина течения в НА и колесе. Показатель изоэнтропы к для рабочего вещества в процессе расширения в КНА изменялся в пределах от 1,03 до 1,3.
Обработка результатов. Для оценки общих характеристик эффективности, из-за большой неравномерности потока на выходе из канала, необходимо провести осреднение параметров в сечениях проточной части. По методике из [5] осреднение выходной скорости c1, скоростного коэффициента ф и угла выхода потока а1 выполняется по количеству движения:
Ks
s
5,863R,
a
(1)
S1
S1
0
0
S1
G=
Р1С1 sin a1ds, a1c = arctg(Im/Iu)
0
c1c =
vm+iu/G, ф1с = c1c/c1
1c =
s
104 ISSN 0236-3941. Вестник МГТУ им.Н.Э.Баумана. Сер. “Машиностроение”. 2015. №1
>
X < х1 L X 1 -— ( — -•
1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2
Рис. 2. Зависимость коэффициента расхода ^ от приведенной скорости Aie:
• — моделирование с гидравлически гладким каналом; ■ — моделирование с шероховатым каналом; х — результаты эксперимента
Рис. 3. Зависимость скоростного коэффициента ф от приведенной скорости A1s (обозначения см. рис. 2)
Индекс “1с” означает, что параметр в сечении 1-1 (см. рис. 1) осредняется по количеству движения, а индекс “1s” — параметр при изоэнтропийном течении в канале. Поскольку во всех рассчитываемых режимах течение в канале КНА сверхзвуковое, то теоретический расход Gs определяется как критический расход через сопло. Коэффициент расхода ^ = G/Gs. После обработки результатов были построены графики зависимостей коэффициента расхода ^ (рис. 2), скоростного коэффициента ф (рис. 3) и коэффициента потерь Z (рис. 4) в канале КНА от приведенной скорости Ais.
Характер изменения коэффициента потерь получился аналогичным представленному в работе [2] (где минимуму потерь соответствовало значение Z = 0,15), однако для фреона область минимума потерь оказалась смещена в сторону больших скоростей и соответственно больших перепадов давлений (минимальному значению потерь соот-
Рис. 4. Зависимость коэффициента потерь Z от приведенной скорости A1s:
• — моделирование с гидравлически гладким каналом; ■ — моделирование с шероховатым каналом; ▲ — результаты испытаний в ОАО “Криогенмаш”
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. № 1 105
ветствует приведенная скорость A1s = 1,98 и степень расширения в КНА пКНА = 4,3). Коэффициент расхода для всего рассчитанного диапазона режимов остается практически постоянным как для гладкого, так и для шероховатого каналов. Различие расходов, полученных в результате моделирования и эксперимента, может быть обусловлено отклонением расчетной геометрии от реальной (разностью размеров, шероховатостей).
Выводы. Эффективность КНА в составе фреоновой теплоутилизирующей силовой турбины сопоставима с эффективностью КНА криогенных детандеров, что расширяет область применения НА данного типа, причем максимум эффективности на фреоне смещен в сторону более высоких скоростей и степеней расширения в НА, что хорошо для турбомашин, работающих в цикле Ренкина (ORC).
Минимальное значение Z = 0,143 — для гидравлически гладкого канала, Z = 0,217 — для канала с шероховатостью Ks = 15 мкм, соответственно следует стремиться к уменьшению значения эквивалентной песочной шероховатости, так как именно от Ks напрямую зависит величина потерь.
ЛИТЕРАТУРА
1. Кулаков В.М., Миронов И.Ю., Кулаков В.В., Данилович В.И.Расчет потерь в ступени турбодетандера // Химическое и нефтегазовое машиностроение. 1997. №4. С. 51-55.
2. Машины низкотемпературной техники. Криогенные машины и инструменты / А.М. Архаров, И.А. Архаров, А.Н. Антонов и др. / под общ. ред. А.М. Архарова и И.К. Буткевича. М.: Изд-во МГТУ им. Н.Э. Баумана, 2011. 582 с.
3. Шлихтинг Г Теория пограничного слоя. М.: Наука, 1974. 712 с.
4. Adams T, Grant C., Watson H. A Simple Algorithm to Relate Measured Surface Roughness to Equivalent Sand-grain Roughness // International Journal of Mechanical Engineering and Mechatronics. 2012. 1. Р. 66-71. DOI: 10.11159 / ijmem. 2012.008
5. Самойлович Г.С. Гидрогазодинамика. М.: Машиностроение, 1990. 384 с.
6. Давыдов А.Б., Пересторонин Г.А., Стулов В.Л., Шерстюк А.Д.Центростреми-тельные турбодетандеры / под общ. ред. В.Н. Удута. М.: Колос-Пресс, 2002. 312 с.
REFERENCES
[1] Kulakov V.M., Mironov I.Yu., Kulakov V.V., Danilovich V.I. Calculation of losses in the turboexpander stage. Khimicheskoe i neftegazovoe mashinostroenie [Chemical and petroleum engineering], 1997, no. 4, pp. 51-55 (in Russ.).
[2] Arkharov A.M., Arkharov I.A., Antonov A.N., Butkevich I.K., et al. Mashiny nizkotemperaturnoy tekhniki. Kriogennye mashiny i instrumenty [Cryogenic machines and enginery]. Moscow, MGTU im. N.E. Baumana Publ., 2011. 582 p.
[3] Schlichting G. Teoriya pogranichnogo sloya [Boundary layer theory]. Moscow, Nauka Publ., 1974. 712 p.
[4] Adams T., Grant C., Watson H.A. Simple Algorithm to Relate Measured Surface Roughness to Equivalent Sand-grain Roughness. International J. of Mechanical Engineering and Mechatronics, 2012, no. 1, pp. 66-71. DOI: 10.11159/ijmem.2012.008
106 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №1
[5] Samoylovich G.S. Gidrogazodinamika [Fluid and gas dynamics]. Moscow, Mashinostroenie Publ., 1990. 384 p.
[6] Davydov A.B., Perestoronin G.A., Stulov V.L., Sherstyuk A.N., Udut V.N., ed. Tsentrostremitel'nye turbodetandery [Centripetal turboexpanders]. Moscow, Kolos-Press Publ., 2002. 312 p.
Статья поступила в редакцию 15.07.2014
Костенко Алексей Андреевич — аспирант кафедры “Холодильная и криогенная техника, системы кондиционирования и жизнеобеспечения” МГТУ им. Н.Э. Баумана, инженер НИИ механики МГУ им. М.В. Ломоносова. Автор трех научных работ в области холодильной техники и турбостроения.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
НИИ механики МГУ им. М.В. Ломоносова, 119192, Москва, Мичуринский пр-т, д. 1. Kostenko A.A. — post-graduate of “Refrigeration, Cryogenic Technology, Air Conditioning and Life Support Systems” department of the Bauman Moscow State Technical University, engineer of the Research Institute of Mechanics of Lomonosov Moscow State University. Author of 3 publications in the field of refrigeration technology and turbine designing.
Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
Research Institute of Mechanics of Lomonosov Moscow State University, Michurinsky prosp. 1, Moscow, 119192 Russian Federation.
Розеноер Татьяна Михайловна — канд. техн. наук, доцент кафедры “Холодильная и криогенная техника, системы кондиционирования и жизнеобеспечения” МГТУ им. Н.Э. Баумана, ведущий инженер НИИ механики МГУ им. М.В. Ломоносова. Автор более 30 научных работ в области криогенной техники, газодинамики и турбостроения.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
НИИ механики МГУ им. М.В. Ломоносова, Российская Федерация, 119192, Москва, Мичуринский пр-т, д. 1.
Rozenoer T.M. — Cand. Sci. (Eng.), assoc. professor of “Refrigeration, Cryogenic Technology, Air Conditioning and Life Support Systems” department of the Bauman Moscow State Technical University, principal engineer of the Research Institute of Mechanics of Lomonosov Moscow State University. Author of more than 30 publications in the field of cryogenics, gas dynamics, and turbine engineering.
Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
Research Institute of Mechanics of Lomonosov Moscow State University, Michurinsky prosp. 1, Moscow, 119192 Russian Federation.
Леонов Виктор Павлович — канд. техн. наук, доцент кафедры “Холодильная и криогенная техника, системы кондиционирования и жизнеобеспечения” МГТУ им. Н.Э. Баумана. Автор более 60 научных работ в области криогенной и холодильной техники.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
Leonov V.P. — Cand. Sci. (Eng.), assoc. professor of “Refrigeration, Cryogenic Technology, Air Conditioning and Life Support Systems” department of the Bauman Moscow State Technical University. Author of more than 60 publications in the field of cryogenic and refrigeration technologies.
Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. № 1 107
