CC BY
18
DOI: 10.24143/2073-1574-2021-1-74-80 УДК 621.515.001.5
ТРЕХМЕРНОЕ ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ МАЛОРАСХОДНОЙ ЦЕНТРОСТРЕМИТЕЛЬНОЙ ТУРБИНЫ С ЧАСТИЧНЫМ ОБЛОПАЧИВАНИЕМ РАБОЧЕГО КОЛЕСА
А. А. Крюков, С. В. Чехранов
Дальневосточный государственный технический рыбохозяйственный университет, Владивосток, Российская Федерация
Центростремительные малорасходные турбины - это надежные, производительные, малогабаритные приводы для различных агрегатов и устройств. Практика производства и использования центростремительных ступеней турбин выдвигает цель по повышению эффективности, упрощению и удешевлению технологии изготовления ступеней. Повышение эффективности центростремительных турбин требует решения задач по совершенствованию проточной части и расчетного определения оптимальной геометрии и режимов работы потока рабочего колеса. Значительная часть центростремительных малорасходных турбин относится к парциальным, что приводит к дополнительным потерям. Использование турбинных ступеней безвентиляционного типа позволяет уменьшить потери, вызванные парциальностью, посредством совершенствования конструктивных параметров. К подобным ступеням относятся турбины с частичным облопачиванием рабочего колеса, которые наделены некоторыми отличительными особенностями по сравнению с обычными парциальными турбинами. Определение потери кинетической энергии в сопловых и рабочих каналах ступени с помощью пакета ANSYS CFX является инновационным путем повышения эффективности проектирования данных турбин. Геометрическая модель построена с использованием ANSYS Design Modeler, выбрана сетка, заданы граничные условия. Приведены зависимости коэффициентов скорости соплового аппарата и рабочего колеса турбины с частичным облопачиванием рабочего колеса в диапазоне изменения степени расширения ступени п = 1,5^2,5 и степени парциальности в диапазоне 0,059^1. В результате аппроксимации получены эмпирические зависимости, учитывающие влияние степени парциальности и степени расширения ступени на коэффициенты скорости соплового аппарата и рабочего колеса. Полученные итоговые эмпирические зависимости незаменимы при моделировании центростремительных турбин.
Ключевые слова: малорасходная турбина, сопловой аппарат, рабочее колесо, полином, численный эксперимент, степень парциальности.
Для цитирования: Крюков А. А., Чехранов С. В. Трехмерное численное моделирование малорасходной центростремительной турбины с частичным облопачиванием рабочего колеса // Вестник Астраханского государственного технического университета. Серия: Морская техника и технология. 2021. № 1. С. 74-80. DOI: 10.24143/2073-1574-2021-1-74-80.
Введение
В современных условиях имеется значительный объем прикладного материала, который позволяет рассмотреть общие принципы проектирования парциальных центростремительных турбин. Турбинные ступени с частичным облопачиванием рабочего колеса (РК) занимают особое место в решении проблем с потерями от парциальности [1]. Эти ступени отличаются от традиционных парциальных не только конструктивно, но и в плане протекающих в них газодинамических процессов, что побуждает необходимость в их моделировании в проточной части ступени этого типа.
В такой ступени рабочее тело (в настоящем исследовании в качестве рабочего тела применяется воздух) подается по всей длине окружности соплового аппарата (СА). В сопловых каналах, расположенных напротив активной дуги РК, течение потока происходит аналогично потоку в СА традиционных парциальных турбин. В сопловых каналах, расположенных напротив неподвижного диска РК, поток ограничен осевым и радиальным зазорами. Перемещение рабочего тела в этих соплах будет объясняться оттоком его через зазоры.
При большей площади выходного участка сопла и поперечного кольцевого зазора между венцами увеличивается интенсивность перетекания рабочего тела из сопел в зазоры даже на неактивной дуге РК. Это явление необходимо учитывать для определения значения перетечки рабочего тела из неактивной области в активную дугу.
Цель работы
По аналогии с ранее проведенными численными экспериментами [2, 3] создается модель для определения коэффициентов скорости СА и РК турбинной ступени с частичным облопачива-нием РК путем аппроксимирования экспериментальных данных, полученных на имитационном стенде с использованием программной системы конечно-элементного анализа ANSYS CFX [4]. Коэффициенты ф и у определяются методом численного моделирования, что позволяет рассматривать каждый отдельный сопловой и рабочий каналы и определять их газодинамические характеристики. Суммарные коэффициенты скорости ф и у могут быть определены как значения ф и у в каждом сопловом и рабочем каналах, осредненные по активной дуге. Эмпирическая зависимость, учитывающая влияние степени парциальности s и расширения п, позволит ее использовать при моделировании в многорежимной оптимизации малорасходных турбин.
Постановка задачи
В работе [5] исследуются модели, в которых коэффициенты скорости СА и РК приняты по результатам обобщений экспериментальных исследований традиционных парциальных турбин. Используя численный эксперимент, возможно определить отличительные особенности характера течения потока от традиционных парциальных турбин и принять решение о необходимости полуэкспериментального исследования характеристик потока в СА и РК с помощью имитационного стенда.
Имитационный стенд создан на основе программной системы ANSYS CFX. С помощью CAD-моделирования в системе ANSYS Design Modeler построены геометрические модели. На рис. 1 приведены ступени с различной степенью парциальности в порядке ее увеличения от 0,059 до 1,00.
Рис. 1. Турбинные ступени с частичным облопачиванием РК c различной степенью парциальности в диапазоне е = 0,059^1
Для выполнения газодинамических расчетов по средствам сеткогенератора ANSYS Meshing создается расчетная сетка. Расчетная сетка в моделируемой турбинной ступени удовлетворяет основным требованиям, обеспечивающим необходимую достоверность получаемых данных.
Граничные условия, такие как параметры торможения, рабочее тело, параметры на выходе и частота вращения, схожи с ранее проведенными исследованиями [2, 3].
При нестационарном взаимодействии СА и РК шаг по времени соотносится с шагом по пространству с использованием критерия Куранта-Фридрихса-Леви (CFL). Критерий CFL (число Куранта) позволяет произвести расчет по десяти характерным положениям РК относительно СА.
ISSN2073-1574. Вестник АГТУ. Сер.: Морская техника и технология. 2021. № 1 Результаты исследования
Используя результаты предыдущего исследования [3] по определению значения коэффициента ф в каждом сопловом канале и суммарного коэффициента скорости СА, осредненного по дуге, подобным образом определяем коэффициент скорости РК у. Зависимости коэффициентов скорости СА ступени с частичным облопачиванием РК в диапазоне изменения пТ =1,5^2,5 и степенью парциальности в диапазоне 0,059^1 приведены на рис. 2.
б
Рис. 2. Двухмерная зависимость коэффициента скорости СА: а - от степени расширения пТ; б - от степени парциальности е
Необходимость в получении полиномной зависимости обусловлена отличием течения потока между обычной парциальной турбиной и турбиной с частичным облопачиванием РК. Полиномы, полученные из графиков зависимости коэффициента ф, могут быть выражены функциями ф = /(п) для каждой степени парциальности в диапазоне 0,059^1,0 и ф = /(е) для соответствующей степени расширения в диапазоне 1,5^2,5:
ф(пТ = 1,5) = -0,1923в4 + 0,5862в3 - 0,6984в2 + 0,3973в + 0,7874, ф( пТ = 2,0) = -0,0657в4 + 0,2709в3 - 0,4311в2 + 0,3106в + 0,8013, ф(пТ = 2,5) = -0,2503в4 + 0,8130в3 - 0,0097в2 + 0,5754в + 0,7757
(1)
ф(в = 0,059) = 0,0149п3Т - 0,126Ц + 0,3247пТ + 0,5569, ф(в = 0,118) = 0,0185п3Т -0,1639п2т + 0,4241пТ + 0,4905, ф(в = 0,206) = -0,0675п3Т + 0,4241п2т - 0,8834пТ +1,4642, ф(в = 0,412) = -0,0500п3Т + 0,3267п2Т -0,7172пТ +1,3982, ф(в = 1,00) = 0,0314п3Т -0,163Ц + 0,2428пТ + 0,8015
(2)
а
Общая математическая модель, которая позволила бы определить коэффициент скорости ф в зависимости от двух определяющих факторов (е и пТ), может выражаться зависимостью ф = /(е, пТ). Представленная на рис. 3 двухпараметрическая зависимость ф = /е, пТ) аппроксимируется кубическим полиномом:
ф = /(в,пТ) = 0,7686- 0,03433пТ + 0,87в + 0,04796пТ -0,3557впТ - 0,8089в2 -- 0,01233пТ + 0,06528пТв + 0,06176пТв2 + 0,3565в3.
Рис. 3. Трехмерная зависимость коэффициента скорости СА от степени расширения пТ и степени парциальности е
На графике зависимости коэффициента скорости СА от степени расширения наблюдается максимальное значение ф для турбинных ступеней со степенью парциальности е = 0,059^0,118 при пТ = 2. Для остальных ступеней со степенью парциальности е = 0,206^1,00 максимум ф наблюдается при пТ = 1,5. Это объясняется тем, что в сопловой дозвуковой решетке оптимальное значение числа Маха (М = 0,55^0,6) приходится на соответствующие степени расширения пТ. Подобным образом определяется эмпирическая зависимость для коэффициента скорости РК (рис. 4).
б
Рис. 4. Двухмерная зависимость коэффициента скорости РК: а - от степени расширения пТ; б - от степени парциальности е
Для коэффициента скорости РК у зависимости у = длт) и у = де) также для каждой степени парциальности и соответствующей степени расширения могут быть представлены кубическими полиномами:
у (пТ = 1,5) = 0,2086в3 - 0,5203в2 + 0,4148в + 0,5862, у (пТ = 2,0) = 0,1015в3 - 0,2918в2 + 0,2639в + 0,6325, у(пТ = 2,5) = 0,2970в3 - 0,5830в2 + 0,3590в + 0,6679
(4)
у(в = 0,059) = -0,066п3Т + 0,389п2Т -0,672пТ + 0,9625, у(в = 0,118) = 0,0794п3Т - 0,418п2т + 0,7797пТ + 0,1387, у(в = 0,206) = 0,0713п3Т -0,380Ц + 0,7286пТ + 0,1639, у(в = 0,412) = -0,025п3Т + 0,169Ц -0,3293пТ + 0,8860, у(в = 1,00) = 0,0518п3Т -0,2875п2т + 0,5640пТ + 0,3175
(5)
Представленная на рис. 5 двухпараметрическая зависимость у = /(е, лт) аппроксимируется кубическим полиномом:
ш = /"(в,пТ) = 0,4515 + 0,156пТ + 0,55в-0,069пТ -0,083впТ -0,0617в2 + Т 7 Т Т Т Т ' (6) + 0,017пТ - 0,007пТв + 0,07пТв2 + 0,2099в3.
v
0.74 -0 72 -0.7 -0 68 0.66 -0 64 -0.62 -
2.5
Рис. 5. Трехмерная зависимость коэффициента скорости РК от степени расширения лт и степени парциальности е
Ввиду того, что рабочие каналы постоянно открыты и не работают в режиме «открытие -закрытие», графики зависимости отличаются. На графике зависимости у от пт наблюдается максимальное значение у для всех турбинных ступеней при л = 2,5. Это объясняется тем, что для рабочей трансзвуковой решетки оптимальное значение числа Маха М = 0,90^0,95, поэтому по мере увеличения пт будет увеличиваться значение у, до достижения оптимума числа Маха. При уменьшении значения л в рабочие каналы поступает газ с более низким давлением, т. к. основной перепад давлений выработан в СА, что и приводит к более низкому значению коэффициента у.
Сравнение данных, полученных в ходе численного эксперимента, с расчетом по формулам (1)-(6) проиллюстрировало, что приведенные модели в настоящей работе для определения коэффициентов скорости СА и РК адекватны, т. к. расхождение между расчетом и численным экспериментом не превышает 2,5 % в широком диапазоне степени парциальности и степени расширения. Поэтому созданные двухпараметрические зависимости ф = /(е, лт) и у = /(е, лт) можно считать адекватными и рекомендовать к прикладному использованию при моделировании центростремительных турбин.
Заключение
По результатам анализа трехмерного моделирования можно сделать следующие выводы:
- получены однопараметрические зависимости для определения коэффициентов скорости СА ф и РК у для каждой степени парциальности ф = /(е), у = /(е) и степени расширения ф = /пт), у = /пт);
- полученные значения скоростных коэффициентов СА и РК при различных значениях е в диапазоне изменения л от 1,5 до 2,5 могут быть аппроксимированы кубическими полиномами (3) и (6);
- при сравнении опытных данных с расчетом по формулам (1)-(6) делаем вывод о целесообразности использования этих математических зависимостей, т. к. погрешность вычисления не превышает значения численного эксперимента более чем на 2,3 %.
СПИСОК ЛИТЕРА ТУРЫ
1. Рассохин В. А. Турбины конструкции ЛПИ: преимущества, характеристики, опыт разработки и применение // Тр. Санкт-Петербург. политехн. ун-та. 2004. № 491. С. 152-161.
2. Крюков А. А., Чехранов С. В. Численное исследование течения потока в ступени центростремительной турбины с частичным облопачиванием рабочего колеса // Мор. интеллектуал. технологии. 2020. № 4 (50). Т. 1. С. 114-120.
3. Kryukov A. A. et al. 2G2G IOP Conf. Ser.: Mater. Sci. Eng. 986 G12G45 IOP Conf. Series: Materials Science and Engineering. 986 (2G2G) G12G45. URL: https://iopscience.iop.org/article/1G.1G88/1757-899X/986/1/G12G45 (дата обращения: 1G.G1.2G21).
4. Епифанов А. А., Кириллов А. И., Рассохин В. А. Pасчет трехмерного течения в ступенях малорасходных турбин // Науч.-техн. вед. Cанкт-Петербург. политехн. ун-та. Наука и образование. 2G12. № 1. C. 65-7G.
5. Чехранов С. В., Симашов P. P. Математическая модель радиальной малорасходной турбины с частичным облопачиванием рабочего колеса // Трансп. дело Pоссии. 2G15. № 6. C. 16G-164.
Отатья поступила в редакцию 19.G1.2G21
ИНФОРМАЦИЯ ОБ АВТОРАХ
Крюков Алексей Алексеевич - Россия, 690087, Владивосток; Дальневосточный государственный технический рыбохозяйственный университет; аспирант кафедры судовых энергетических установок; dsl918@ya.ru.
Чехранов Сергей Валентинович - Россия, 690087, Владивосток; Дальневосточный государственный технический рыбохозяйственный университет; д-р техн. наук, профессор; профессор кафедры судовых энергетических установок; turboroom@yandex.ru.
THREE-DIMENSIONAL NUMERICAL SIMULATION OF LOW-CONSUMPTION INFLOW TURBINE WITH PARTIAL BLADING OF RUNNER
A. A. Kriukov, S. V. Chekhranov
Far Eastern State Technical Fisheries University, Vladivostok, Russian Federation
Abstract. The article describes the inflow turbines as reliable, productive, small-sized drives for various units and devices. The practice of manufacturing and using inflow turbine stages sets out the goal of increasing efficiency, simplifying technology and reducing cost of manufacturing stages. Increasing the efficiency of inflow turbines requires solving the problems of improving the flow part and calculating the optimal geometry and operating modes of the impeller flow. A significant part of low-consumption inflow turbines is partial, which leads to additional losses. The use of fanless turbine stages reduces the losses caused by partiality by means of improving the design parameters. Such stages include turbines with partial blading of the runner, which have the distinctive features compared to conventional partial turbines. An innovative way to improve the design efficiency of these turbines is to determine the loss of kinetic energy in the nozzle and runner using the ANSYS CFX package. The geometric model is constructed using ANSYS Design Modeler, a grid is selected, and boundary conditions are set. There are shown the dependences of speed coefficients of a nozzle diaphragm and a runner of the turbine with partial blading in the range of expansion stage nt =1,5^2,5, and the admission intensity in the range 0,059^1. As a result of approximation, the empirical dependences are obtained that take into account the influence of the degree of partiality and the degree of expansion of the stage on the coefficients of speed of the nozzle diaphragm and the runner. The resulting empirical dependences are indispensable for modeling an inflow turbine.
Key words: low-consumption turbine, nozzle diaphragm, rotor wheel, polynomial, numerical experiment, admission intensity.
For citation: Kriukov A. A., Chekhranov S. V. Three-dimensional numerical simulation of low-consumption inflow turbine with partial blading of runner. Vestnik of Astrakhan State Technical University. Series: Marine Engineering and Technologies. 2021;1:74-80. (In Russ.) DOI: 10.24143/2073-1574-2021-1-74-80.
REFERENCES
1. Rassokhin V. A. Turbiny konstruktsii LPI: preimushchestva, kharakteristiki, opyt razrabotki i prime-nenie [Turbines of LPI design: advantages, characteristics, development experience and application]. Trudy Sankt-Peterburgskogo politekhnicheskogo universiteta, 2004, no. 491, pp. 152-161.
2. Kriukov A. A., Chekhranov S. V. Chislennoe issledovanie techeniia potoka v stupeni tsentrostremi-tel'noi turbiny s chastichnym oblopachivaniem rabochego kolesa [Numerical study of flow in centripetal turbine stage with partial blading of impeller]. Morskie intellektual'nye tekhnologii, 2020, no. 4 (50), vol. 1, pp. 114-120.
3. Kryukov A. A. et al. 2020 IOP Conf. Ser.: Mater. Sci. Eng. 986 012045 IOP Conf. Series: Materials Science and Engineering. 986 (2020) 012045. Available at: https://iopscience.iop.org/article/10.1088/1757-899X/986/1/012045 (accessed: 10.01.2021).
4. Epifanov A. A., Kirillov A. I., Rassokhin V. A. Raschet trekhmernogo techeniia v stupeniakh maloraskhodnykh turbin [Calculating three-dimensional flow in stages of low-flow turbines]. Nauchno-tekhnicheskie vedomosti Sankt-Peterburgskogo politekhnicheskogo universiteta. Nauka i obrazovanie, 2012, no. 1, pp. 65-70.
5. Chekhranov S. V., Simashov P. P. Matematicheskaia model' radial'noi maloraskhodnoi turbiny s chastichnym oblopachivaniem rabochego kolesa [Mathematical model of radial low-flow turbine with partial blading of impeller]. Transportnoe delo Rossii, 2015, no. 6, pp. 160-164.
The article submitted to the editors 19.01.2021
INFORMATION ABOUT THE AUTHORS
Kriukov Aleksei Alekseevich - Russia, 690087, Vladivostok; Far Eastern State Technical Fisheries University; Postgraduate Student of the Department of Marine Power System; dsl918@ya.ru.
Chekhranov Sergei Valentinovich - Russia, 690087, Vladivostok; Far Eastern State Technical Fisheries University; Doctor of Technical Sciences, Professor; Professor of the Department of Marine Power System; turboroom@yandex.ru.
