Энергетическое, металлургическое и химическое машиностроение
УДК 621.515 doi: 10.18698/0536-1044-2021-3-39-52
Первичное проектирование безлопаточных диффузоров центробежных компрессорных ступеней Методом универсального моделирования*
О.А. Соловьёва1, К.В. Солдатова2, Ю.Б. Галеркин1, А.Ф. Рекстин1
1 ФГАОУ ВО «Санкт-Петербургский политехнический университет Петра Великого»
2 Мельбурнский политехнический институт
Primary Design of Vaneless Diffusers of Centrifugal Compressor Stages by the Universal Modeling Method
O.A. Solovyeva1, K.V. Soldatova2, Yu.B. Galerkin1, A.F. Rekstin1
1 Federal State Autonomous Educational Institution of Higher Education — Peter the Great St. Petersburg Polytechnic University
2 Melbourne Institute of Technology
Безлопаточные диффузоры промышленных центробежных компрессоров чаще всего состоят из начального сужающегося участка и основной части с параллельными стенками. Предложено выбрать такую ширину основной части, при которой поток в диффузоре остается безотрывным на границе помпажа. Для реализации метода первичного проектирования получена эмпирическая формула для определения минимального безотрывного угла потока в зависимости от относительной ширины диффузора. Выполнено первичное проектирование восемнадцати ступеней, охватывающих практически важный диапазон параметров. Методом универсального моделирования рассчитаны соответствующие газодинамические характеристики, проанализированы размеры и углы потока. Предложенный метод первичного проектирования интегрирован в Метод универсального моделирования и используется в проектной практике.
Ключевые слова: ступень центробежного компрессора, первичное проектирование, безлопаточный диффузор, коэффициент расхода, КПД ступени, граница помпажа
Vaneless diffusers of industrial centrifugal compressors most often consist of a tapered inlet section and a parallel-walled main section. The study proposes to choose such a width of the main section, at which the flow in the diffuser remains unseparated at the surge limit. To implement the primary design method, an empirical formula was obtained to determine the minimum continuous flow angle depending on the relative width of the diffuser. The primary design of eighteen stages was completed, covering a practically important range of parameters. The corresponding gas-dynamic characteristics were calculated by the universal modeling method, the dimensions and angles of the flow were analyzed. The proposed primary design method is integrated into the universal modeling method and is used in design practice.
* Работа выполнена при поддержке гранта Президента Российской Федерации для молодых кандидатов наук МК-1893.2020.8
Keywords: centrifugal compressor stage, primary design, vaneless diffuser, flow coefficient, stage efficiency, surge limit
В СССР систему взглядов на рабочий процесс промышленных центробежных компрессоров и основы проектирования проточной части сформулировал главный конструктор Невского завода в Ленинграде (НЗЛ) В.Ф. Рис [1, 2]. Приоритет отдавался высокому коэффициенту полезного действия (КПД) в номинальном режиме и минимизации радиальных размеров, поэтому НЗЛ и другие производители выпускали компрессоры с лопаточными диффузорами. В то же время внимание исследователей привлекали безлопаточные диффузоры (БЛД).
В монографиях [3, 4] для анализа использованы две разные модели с противоречащими результатами. Аналогия с прямоосным плоским диффузором показывала, что минимальный коэффициент потерь соответствует некоторому углу потока в БЛД, при котором оптимален угол раскрытия эквивалентного плоского диффузора. В БЛД промышленных компрессоров угол потока незначительно меняется по радиусу, поэтому эквивалентный угол раскрытия V определяется из выражения
V
tg2 =
к _ Е
п У п
Г2 _ Г1
sin1,5 а + yfn'
sin —
0,2
ОД
0
\ \
\
\
\
_ \
\ \
Ч ,
1 1 1 1
0,03
0,04
0,05
0,06
0,07 Ф
Рис. 1. Экспериментальная ( ) и расчетная ( ) зависимости коэффициента потерь БЛД С модельной ступени 048-0,48-0,29 серии 20СЕ от условного коэффициента расхода Ф
полного давления на входе в БЛД на радиальном расстоянии от рабочего колеса (РК), равном 5 % его радиуса. Измеренная характеристика коэффициента потерь БЛД действительно имеет минимум. Результат испытания одной из модельных ступеней 048-0,48-0,29 серии 20СЕ (штриховая кривая), выполненного в Проблемной лаборатории компрессорострое-ния ЛПИ [5, 6], приведен на рис. 1. Там же показана аналогичная зависимость (сплошная кривая), полученная путем расчета Методом универсального моделирования.
Оптимальному для БЛД условному коэффициенту расхода Ф = 0,055 соответствует угол потока а = 30°. Однако расчеты того же автора [3, 4] методом теории пограничного слоя показали, что коэффициент потерь БЛД монотонно уменьшается с увеличением угла потока в БЛД вплоть до 90°.
Анализ характера течения потока и одномерные расчеты в работе [7] также показали, что с ростом угла потока снижаются потери трения и повышается способность потока сопротивляться отрыву. В БЛД с шириной b(r) = const и углом потока а = const коэффициент потерь
где /2, /1 — площади; ?2, г — радиусы; а — угол потока в БЛД.
В соответствии с рассмотренной концепцией оптимальный угол потока в БЛД меньше 90°, что подтверждалось результатами измерения
fi
4b sin a V D4
(1)
где Я — коэффициент трения; Ь — относительная ширина БЛД; В4 — относительный диаметр на выходе из БЛД, Б4 = Л4/02 (Б4 — диаметр на выходе из БЛД, Б2 — диаметр на входе в БЛД).
Одномерный расчет коэффициента потерь БЛД С в зависимости от угла потока а при коэффициенте трения Я = 0,015, относительном диаметре на выходе из БЛД В4 = 1,65 и относительной ширине Ь3 = Ь2 = 0,01, 0,02, 0,04, 0,08 [7] выявил, что оптимальный угол потока а = 90° (рис. 2).
Минимум потерь при угле потока а = 90° подтвердили испытания БЛД с созданием закрутки потока вращающейся системой сеток [8], а также исследования, проведенные в работе [9].
Несоответствие этому закону коэффициента потерь БЛД, измеренного в центробежной ступени (см. рис. 1), связано с некорректным измерением полного давления на границе РК и БЛД. Смешение следа и ядра потока происхо-
С
Рис. 2. Зависимость коэффициента потерь БЛД £ от угла потока а при коэффициенте трения
X = 0,015, относительном диаметре на выходе из БЛД = 1,65 и относительной ширине Ь3 = Ь2 = 0,01 (1), 0,02 (2), 0,04 (3), 0,08 (4)
дит вниз по потоку после измерения полного давления. Потери смешения РК происходят в БЛД, но это не потери БЛД. Разница между штриховой и сплошной линиями на рис. 1 соответствует потерям смешения в РК.
Расчеты по соотношениям, аналогичным формуле (1), показывают, что потери трения в БЛД с шириной, равной высоте лопаток РК, меньше, чем в диффузоре с уменьшающейся шириной. Хотя при Ь < Ь2 угол потока увеличивается, и траектория потока сокращается, уменьшение Ь повышает потери.
Однако было показано, что БЛД с шириной Ь > Ь2 можно применять лишь в малорасходных ступенях [10, 11]. Уменьшение угла потока ведет к отрыву потока, росту потерь и вращающемуся срыву. Это подтверждено результатами экспериментальных исследований [12].
Таким образом, оптимальным решением представляется классический БЛД с параллельными стенками и расстоянием между ними, равным высоте лопаток колеса. Однако еще в 1960-е годы [6] обширные эксперименты Проблемной лаборатории компрессорострое-ния показали определенные преимущества БЛД с начальным сужающимся участком и параллельными стенками в основной части.
Применение БЛД с относительной шириной Ь3 = Ь3 / Ь2 < 1 в составе начального суженного участка и основного участка с Ь3 = Ь4 стало общепринятым (рис. 3) [5, 7, 13, 14].
В работах [15-17] предложена математическая модель (ММ), описывающая уменьшение потерь смешения в сужающейся части БЛД.
Заметное влияние сужающейся части БЛД на РК фиксируют прямые измерения (рис. 4) [7].
Задача первичного проектирования (PD) — выбрать правильное соотношение высоты лопаток РК и ширины основного участка БЛД. Выбор формы сужающегося участка и радиального размера БЛД — задача последующих этапов проектирования проточной части.
Цель работы — продемонстрировать результаты исследований в области выбора оптимальной ширины БЛД. Данные о влиянии радиальной протяженности БЛД на его эффективность, полученные ранее, приведены в работе [18].
Существует мнение, что при излишней протяженности узкого БЛД давление перестает повышаться. Расчеты в трудах [19, 20] это не подтверждают. В целом же оптимальная радиальная протяженность БЛД — это предмет совместного исследования диффузора и следующего за ним обратно направляющего аппарата.
Рис. 3. Схема типичного БЛД с суженным начальным участком
т_
А
- °Л Л Л
о Л о
| 1 о Л о Л о 1 1
0,03 0,04 0,05 0,06 0,07 Ф
Рис. 4. Влияние начального участка БЛД на коэффициент напора РК ^ с выходным углом лопаток рл2 = 45° при относительной ширине БЛД Ь3 = 1,000 ( ) и 0,785 ( )
Принцип выбора относительной ширины БЛД при первичном проектировании. В качестве объекта первичного проектирования выбран БЛД ступени промышленного центробежного компрессора, в которой течение дозвуковое.
Если проектируют ступень компрессора, не предназначенного для работы при производительности меньше расчетной, то оптимальное техническое решение заключается в равенстве ширины диффузора БЛД Ьз и высоты лопаток РК Ь2. При снижении производительности угол потока в диффузоре уменьшается, и возникает возможность отрыва потока, что неприемлемо для работы компрессора. Необходимо уменьшить ширину БЛД. Это снизит КПД в расчетном режиме, но позволит работать при производительности меньше расчетной.
В статье [18] предложен такой принцип выбора относительной ширины БЛД: после первичного проектирования РК известен расчетный угол потока на выходе a2des. Расчетный угол потока на входе в БЛД a3des должен быть таким, чтобы при уменьшении расхода на границе помпажа угол потока все еще соответствовал безотрывному течению. Для выбора относительной ширины БЛД Ь3 необходимо решить следующие задачи:
• определить минимальный угол потока на входе в БЛД a3min, при котором нет отрыва потока, — границу безотрывного течения в БЛД — в зависимости от его относительной ширины Ь3 и критериев подобия;
• найти коэффициент запаса по помпажу для РК проектируемой ступени, равный отношению Фсг /Фdes, где Фсг — коэффициент расхода на границе помпажа; Ф^ — расчетный условный коэффициент расхода;
• вычислить угол потока на выходе из РК a2cr при коэффициенте расхода на границе помпажа Фсг ; _ _
• если при Ь3 = Ь угол потока соответствует безотрывному течению a2cr > a3min, то следует принять ширину БЛД, равной высоте лопаток РК на выходе;
• если при Ь3 = Ь2 угол потока не соответствует безотрывному течению a2cr <a3min, то необходимо уменьшить ширину БЛД относительно высоты лопаток РК на выходе, т. е. следует найти такое отношение Ь3 / Ь2 < 1, при котором выполняется условие a2cr =a3min.
Угол потока, при котором начинается его отрыв в БЛД. Обработка результатов расчет-
ного эксперимента с БЛД [19-21] подтвердила известные условия, при которых происходит отрыв потока в БЛД. На рис. 5 показано отрывное и безотрывное течение в БЛД по результатам СБЭ-расчетов при различных значениях относительной ширины Ь3, угла потока а3 и числа Рейнольдса ReЬ2, рассчитанного по высоте Ь2 [10, 19, 20].
Как и при экспериментах [7], СБЭ-расчет обнаруживает отрыв на одной из стенок. На рис. 6 показано изменение касательного напряжения по радиусу Я БЛД (по ограничивающим стенкам БЛД Wall1 (красная кривая) и Wall2 (синяя кривая) с относительной шириной Ь3 = 0,057 при Хс2 = 0,64, угле потока а3 = 10° и числе Рейнольдса ReЬ2 = 520 000.
Начало отрыва соответствует известному из теории пограничного слоя условию нулевого касательного напряжения на стенке. В зоне отрыва потока касательные напряжения отрицательные, так как он движется в противоположном направлении.
Зависимость минимального угла потока ат;п, при котором начинается отрыв потока,
а
Рис. 5. Линии тока в меридиональной плоскости БЛД с различными значениями относительной ширины, угла потока на входе и числа Рейнольдса, рассчитанного по высоте Ь2: а — Ь3= 0,100, а3 = 10°, ЯеЬ2 = 911 000; б — Ь3= 0,100, а3 = 20°, ЯеЬ2 = 911 000; в — Ь3 = 0,057, а3 = 10°, ЯеЬ2 = 521 000; г — Ь3 = 0,057, а3 = 20°, ЯеЬ2 = 521 000; д — Ь3 = 0,014, а3 = 10°, ЯеЬ2 = 130 000
«min. град
Рис. 6. Изменение касательного напряжения, Па, по радиусу Я, м, БЛД с относительной шириной Ь3 = 0,057 при скоростном коэффициенте Хс2 = 0,64, угле потока а3 = 10° и числе Рейнольдса Reb2 = 520 000
от скоростного коэффициента Хс2 при различных значениях относительной ширины БЛД Ь3 и числа Рейнольдса КеЬ2 приведена на рис. 7.
Очевидна сильная зависимость минимального угла потока атщ от относительной ширины БЛД. В узких диффузорах при Ь3 < 0,029 отрыв потока не обнаружен. Быстрое смыкание пограничного слоя усиливает касательные напряжения, препятствующие отрыву потока. В исследованном диапазоне влияние критерия Рейнольдса КеЬ2 на угол атщ не обнаружено. Влияние критерия сжимаемости проявляется в БЛД с относительной шириной Ь3 > 0,08.
Зависимость минимального угла потока аШщ от относительной ширины БЛД Ь3 при скоростном коэффициенте Хс3 = 0,39, 0,64 и 0,82 приведена на рис. 8.
Целью вычислений было не точное определение угла, при котором начинается отрыв потока, а построение ММ расчета коэффициента потерь. Для определения угла отрыва следовало бы рассчитывать течение с шагом угла потока, равным примерно 1°. С учетом большого шага углов потока, при которых проведены расчеты (см. рис. 8), значения а3тщ = / (Ь3) можно считать не зависящими от скоростного коэффициента Хс3 при относительной ширине БЛД Ьз < 0,08.
Рис. 7. Зависимость минимального угла потока атщ, при котором начинается отрыв, от скоростного коэффициента Хс2 при различных значениях относительной ширины БЛД и числа Рейнольдса: X — Ь3 = 0,100, ЯеЬ2 = 373 000...103 000; ■ — Ь3 = 0,057, ЯеЬ2 = 214 000.590 000; ▲ — Ь3 = 0,029, ЯеЬ2 = 106 000.294 000
а„
, град
25 л
20 1
15 "ж •
10 - •
| | 1
0,02
0,04
0,06
0,08
Рис. 8. Зависимости минимальных углов потока amin (точки) и a3min (сплошная линия) от относительной ширины БЛД b3 : — Хс3 = 0,39; ^ — Хс3 = 0,64; Л — Хс3 = 0,82
В узких БЛД отрыв не обнаружен при угле а3 > 10°, но даже на границе помпажа угол потока не должен быть меньше 5°. На этом основании предложена эмпирическая формула
«3min = arctg (0,0875 + 3,5b3 ). (2)
Зависимость a3min = f (b3) показана сплошной линией на рис. 8.
Оценка коэффициента запаса по помпажу.
Положение границы помпажа зависит от расчетного коэффициента теоретического напора РК des Чем больше коэффициент теоретического напора, тем более пологой будет напорная характеристика. При уменьшении расхода
рост потерь меньше компенсируется подводом механической работы. При дальнейшем снижении расхода конечное давление перестает увеличиваться и наступает помпаж.
Эта проблема подробно рассмотрена в публикации [22]. Предложен упрощенный способ оценки границы помпажа как режима с максимальным отношением давлений. В практике расчетов такое понимание границы помпажа общепринято [23, 24]. Далее способ оценки границы помпажа из работы [18] представлен в уточненном виде.
Граница помпажа соответствует коэффициенту расхода, при котором коэффициент по-литропного напора Нр /и| = (Нт - К, )/и\ (где Нр — политропный напор; и2 — окружная скорость; Нт и Н„ — теоретический и потерянный напоры) достигает максимума:
d\\t p йФ
- 0,
выходе из РК в относительном движении; с2 и — относительная скорость потока на выходе из РК при абсолютном и относительном движении; и2 — относительная скорость вращения РК на выходе.
В пределах практически используемой части размерной характеристики зависимость напора от расхода также линейная, что отмечено еще в монографии [1]. Поэтому при моделировании принято считать функцию ¥т = /(Ф) линейной, хотя связь между параметрами Ф и ф2 не вполне линейная.
При линейной зависимости коэффициента теоретического коэффициента напора от условного коэффициента расхода можно записать
Vт - Vт0 "(Vт0 т des )
Ф
Ф
des
где ¥ р — коэффициент политропного напора.
Линейная зависимость коэффициента теоретического напора ут от коэффициента расхода ф2 проанализирована в публикациях [25, 26]. Напорная характеристика и треугольник скоростей при линейной зависимости ¥ т = / (ф2) [26] показаны на рис. 9. Здесь введены следующие обозначения: ¥т0 — коэффициент теоретического напора при нулевом расходе; ф2 и ф2тах — коэффициент расхода на выходе из РК и его максимальное значение; ф2тах~ — максимальное значение коэффициента расхода на выходе из идеального РК с бесконечным числом лопаток; Рт — угол наклона напорной характеристики; р2 — угол потока на
Рл2>90° Рй = 90°
Коэффициент теоретического напора при нулевом расходе Vto зависит от многих факторов, но, главное, этот коэффициент больше у РК с большим коэффициентом напора в расчетном режиме. На основе экспериментальных данных, приведенных в работе [25], предложена приближенная формула, пригодная для целей первичного проектирования при коэффициенте теоретического напора РК vtdes - 0,40...0,70:
Vто - 0,84 + 0,27(vтdes "0,40).
Связь между политропным и теоретическим напорами можно описать следующим образом:
hp Йт hwdei
hw
h
(3)
wdes
где hwdes — расчетный потерянный напор.
Ф2тах Ф2тахоо
Рис. 9. Напорная характеристика (а) и треугольник скоростей при линейной зависимости коэффициента теоретического напора ¥т от коэффициента расхода (б)
В работе [25] для отношения Нц, рекомендовано следующее приближенное выражение:
К = 1 + 8,5 f 1 -Ф
hwdes V Фdes
Соотношение (3) представляется в виде
Vт0 - Vт des
V p =Vx0"
Ф
-Ф-
des
-¥т des (1 -^des )
1 + 8,51 1 -
Ф
Ф
des
(4)
Фс
Ф
= 1 -
des
-1
¥т des
17 (1 -^des )
(5)
ф2сг ~ 92des
Фс
Ф
(6)
des
¥т cr = ¥т0 ( ¥т0 V des )
Фс
Ф,
des
с Ф Л
^ сг
ф2d,
«2сг = arctg
Ф
des
Z
(7)
где
Z = 0,84 + 0,27(des -0,40)- [0,84 + 0,27 (V т des - 0,40 )) т des ]
Фс
Фdes
Приближенное значение расчетного коэф фициента расхода на выходе из РК определяет ся из уравнения неразрывности
Фdes Фdes
ф2des =■
(8)
4£2b2 4-1,18b2 где е2 — коэффициент сжимаемости на выходе из РК, е2 = 1,18.
Относительная высота лопаток на выходе из РК Ь2 при первичном проектировании РК задается эмпирическим соотношением.
Если при условии Ь3 = Ь2 угол а2сг >а3т;п, то применять сужающийся участок не следует. Если_а2сг <а3т;п, то надо определить отношение Ь3/Ь2, при котором эти углы будут равны, исходя из того, что tg а3т;п = tg а2сгЬ2/Ь3. В соответствии с формулой (2) решению подлежит уравнение
- -0,0875 + ^000766+14^02^
Ь3 =- *
7
(9)
где — расчетный КПД ступени.
После дифференцирования формулы (4) коэффициент запаса по помпажу определяется выражением
¥т0
Угол выхода потока из РК на границе помпа-жа. Выбор ширины основного участка диффузора. Угол выхода потока из РК на границе помпажа определяется расходной ф2сг и окружной сг составляющими скоростей:
где фмга — расчетный коэффициент расхода на выходе из РК.
Угол потока на границе помпажа
Безлопаточные диффузоры ступеней с разными параметрами проектирования. Методом первичного проектирования [27-29] определены размеры восемнадцати ступеней в диапазонах параметров проектирования Фdes = 0,150...0,015, Vтdes = 0,45...0,70. По размерам первичного проектирования Методом универсального моделирования рассчитаны газодинамические характеристики. Пример результатов расчета приведен на рис. 10. Здесь показано влияние условного коэффициента расхода Ф и расчетного коэффициента теоретического напора ¥тdes на отношение давлений П, КПД коэффициент внутреннего напора Vi и коэффициент политропного напора VР.
В процессе оптимизации размеры корректируются и КПД может быть повышен. Для анализа БЛД достаточно расчетов характеристик по первичному проектированию. Расчетные параметры, характеризующие работу БЛД, при условном числе Маха РК Мц = 0,70 приведены в таблице.
Значения ряда параметров, приведенных в таблице, рассчитаны по следующим выражениям:
• КПД ступени ^des — по упрощенной модели [27-29];
• относительная высота лопаток на выходе из РК b2 — по формуле первичного проектирования РК (в тексте не приведена);
• относительная ширина БЛД b3 — по формуле первичного проектирования БЛД (9);
• коэффициент расхода на границе помпажа Фсг (PD) — по уравнению первичного проектирования (5).
Остальные параметры, указанные в таблице, взяты из расчета газодинамических характеристик.
— П — ^ II
Рис. 10. Характеристики ступени с расчетными коэффициентами Фае5 = 0,15, ¿ж = 0,45 и размерами проточной части по первичному проектированию
Расчетные параметры БЛД ступеней при условном числе Маха РК М„ = 0,7
Фйев Ут ^ ЦйеБ (РЩ ьг ьъ а3йеэ Фсг (РЩ Фсг/Фаев (PD) И3сг (ММ)
0,15 0,450 0,8098 0,0852 0,0675 39,5 47,7 0,0752 0,501 17,9 18,90
- 0,575 0,8180 0,0852 0,0740 32,8 37,2 0,1037 0,691 19,1 23,80
- 0,700 0,8170 0,0852 0,0775 27,6 30,3 0,1232 0,821 19,7 24,50
0,10 0,450 0,8460 0,0776 0,0466 31,6 47,2 0,0446 0,446 14,1 15,50
- 0,575 0,8510 0,0776 0,0543 25,1 34,7 0,0659 0,659 15,5 20,00
- 0,700 0,8500 0,0776 0,0588 21,6 27,9 0,0802 0,802 16,3 21,90
0,070 0,450 0,8540 0,0723 0,0348 24,4 45,4 0,0302 0,431 11,8 14,50
- 0,575 0,8580 0,0723 0,0423 19,3 32,3 0,0455 0,650 13,3 18,30
- 0,700 0,8560 0,0723 0,0467 16,0 25,1 0,0559 0,799 14,1 18,80
0,045 0,450 0,8450 0,0597 0,0256 19,5 40,9 0,0201 0,467 10,0 12,50
- 0,575 0,8480 0,0597 0,0321 15,1 27,6 0,0298 0,662 11,3 15,50
- 0,700 0,8450 0,0597 0,0357 12,2 20,6 0,0362 0,782 12,0 15,80
0,030 0,450 0,8010 0,0400 0,0227 19,0 32,2 0,0154 0,513 9,5 11,60
- 0,575 0,8020 0,0440 0,0270 13,7 22,1 0,0211 0,703 10,3 13,20
- 0,700 0,7970 0,0482 0,0291 10,1 17,0 0,0249 0,830 10,7 13,40
0,015 0,450 0,6810 0,0260 0,0176 15,6 22,9 0,0092 0,613 8,5 9,97
- 0,575 0,6800 0,0240 0,0193 12,7 15,8 0,0115 0,767 8,8 10,40
0,700 0,6680 0,0230 0,0198 10,6 12,6 0,0130 0,867 8,9 10,40
Диапазон изменения скорости на входе в БЛД С3 = 0,81.0,49, что при условном числе Маха РК Ми = 0,7 определяет дозвуковой характер течения в БЛД.
Результаты расчетов из таблицы представлены в графическом виде на рис. 11, где относительная высота лопаток на выходе из РК Ь2 сопоставлена с относительной шириной БЛД Ь3.
Алгоритм расчета высоты лопаток РК не учитывает влияние коэффициента напора, и при окончательном проектировании значение относительной высоты лопаток Ь2 иногда корректируют. Это не влияет на относительную ширину БЛД Ь3. В диапазоне расхода Ф^ = 0,150.0,045 ширина БЛД практически пропорциональна расчетному коэффициенту расхода.
В ступенях с одинаковыми значениями условного коэффициента расхода при расчетном режиме и высоких коэффициентах напора угол потока в БЛД меньше. Но оказывается, что при больших значениях расчетного теоретического напора е БЛД должен быть шире. Причина в том, что коэффициент запаса по помпажу у высоконапорных РК ближе к единице. Угол потока на границе помпажа у них меньше отличается от расчетного значения.
У малорасходных ступеней отрыв потока в БЛД наступает при небольших углах потока, или не наступает вовсе. Требуется небольшое сужение БЛД. У высокорасходных ступеней углы потока на выходе из РК очень большие. Нет необходимости их значительно увеличивать в диффузоре. У среднерасходных ступеней потребность сужения БЛД максимальная. Хотя
a2des> а3 des
0,08 0,12 Ф^
Рис. 11. Зависимость параметров БЛД от расчетного коэффициента расхода Ф^ : — относительной высоты лопаток на выходе из РК Ь2; , ▲, х — относительной ширины БЛД Ь3 при ут ¡¡е$ = 0,450; 0,575; 0,700 соответственно
50 40 30 20 10
0
■ ■ ■
■ ■
А
А
■ ▲ ■ *
▲ X X
i А
■ ▲ X
■ i А
X ▲ X
А А X
X X 1 1 i
0,04
0,08
0,12
Ф,
des
Рис. 12. Зависимость расчетных углов потока на выходе из РК a2des ( , ▲, ) и на входе
в основной участок БЛД a3des ( , ▲, ) от расчетного коэффициента расхода Ф^ при различных значениях расчетного коэффициента теоретического напора Ут des :
■, ■ — Ут des = 0,450; ▲ , ▲ — Ут des = 0,575; X, X — Ут des = 0,700
при окончательном проектировании значения относительной ширины БЛД могут быть увеличены.
Зависимость расчетных углов потока на выходе из РК а2^ и на входе в основной участок БЛД от расчетного коэффициента расхода Ф^ при ^ = 0,450, 0,575 и 0,700 приведена на рис. 12.
У низконапорных ступеней при Ф^ > 0,025 углы потока в БЛД очень большие. Это связано с тем, что такие ступени имеют коэффициент запаса по помпажу Фсг /Фе = 0,50.0,45, [30], и необходимо обеспечить безотрывное течение в БЛД на границе помпажа. У ступеней со средними и большими коэффициентами напора граница помпажа ближе к расчетному режиму. Поэтому их БЛД могут быть шире, а углы потока в расчетном режиме — меньше.
При окончательном проектировании выбор ширины БЛД может быть существенно изменен. Например, первая ступень многоступенчатого компрессора работает в сравнительно узком диапазоне расходов Фсг /Фе < 1, определяемом запиранием и помпажом в последней ступени [5].
Зависимости минимального угла потока о3тщ, при котором не возникает отрыв в БЛД, и минимального угла а3сг, вычисленного по точной ММ, от расчетного условного коэффициента расхода Ф^ показаны на рис. 13.
®3min> сг
Рис. 13. Зависимость минимальных углов потока «3min ( , , ) и a3cr ( , ▲, ) от расчетного коэффициента расхода Ф^ : ■, ■ — ут des = 0,450; А, ▲ — Ут des = 0,575; х, X — Ут des = 0,700
В принципе эти углы должны были бы совпадать. Однако во всех случаях, приведенных на рис. 13, значения угла a3min =a2cr, вычисленные по формуле (7), меньше чем у значений угла a3cr, рассчитанных по ММ. Основная причина этого заключается в том, что в формулах, связывающих коэффициенты ф2 и Ф, сжимаемость потока учтена приближенно (8) или не учтена (6).
Расчеты характеристик, приведенные в таблице, выполнены при условном числе Маха РК Mu = 0,70. При других значениях Mu соотношения углов аз^ и а3т;п могут различаться.
Выводы
1. Предлагаемые первичные расчеты рациональной ширины БЛД основаны на упрощающих допущениях. В приведенных примерах точного расчета Методом универсального моделирования угла потока на выходе из РК показано, что рациональная ширина БЛД может быть больше, чем по первичному расчету. Также очевидно, что принятое условие является необязательным для всех случаев. Например, в многоступенчатом компрессоре первая ступень может иметь коэффициент запаса по помпажу значительно меньше единицы, но никогда не будет работать на большом удалении от расчетного режима из-за помпажа на последних ступенях.
2. Разработанный метод первичного проектирования интегрирован в Метод универсального моделирования и используется в проектной практике для расчета и проектирования промышленных дозвуковых центробежных компрессоров.
Литература
[1] Рис В.Ф. Центробежные компрессорные машины. Ленинград, Машиностроение, 1964.
336 с.
[2] Рис В.Ф. Центробежные компрессорные машины. Ленинград, Машиностроение, 1981.
351 с.
[3] Ден Г.Н. Механика потока в центробежных компрессорах. Ленинград, Машинострое-
ние, 1973. 269 с.
[4] Ден Г.Н. Проектирование проточной части центробежных компрессоров. Ленинград,
Машиностроение, Ленингр. отд-ние, 1980. 230 с.
[5] Галеркин Ю.Б. Турбокомпрессоры. Рабочий процесс, расчет и проектирование проточ-
ной части. Москва, Изд-во КХТ, 2010. 596 с.
[6] Галеркин Ю.Б., ред. Труды научной школы компрессоростроения СПбГТУ. Москва,
Изд-во КХТ, 2000. 443 с.
[7] Селезнев К.П., Галеркин Ю.Б. Центробежные компрессоры. Ленинград, Машинострое-
ние, 1982. 271 с.
[8] Jansen W. Steady fluid flow in a radial vaneless diffuser. Transactions ASME. Journal of basic
engineering, 1964, Sept, pp. 607-619.
[9] Johnston J.P., Dean R.C. Losses in vaneless diffusers of centrifugal compressors and pumps.
Transactions ASME, Journal of Engineering for Power, 1966, vol. 88, pp. 49-60.
[10] Соловьева О.А., Галеркин Ю.Б. Выбор оптимальных соотношений размеров безлопаточного диффузора малорасходной центробежной компрессорной ступени. Компрессорная техника и пневматика, 2017, № 5, с. 11-15.
[11] Соловьева О.А., Галеркин Ю.Б. Расчетный анализ безлопаточного диффузора малорасходной центробежной компрессорной ступени. Компрессорная техника и пневматика, 2017, № 3, с. 10-13.
[12] Лившиц С.П. Аэродинамика центробежных компрессорных машин. Москва, Ленинград, Машиностроение, 1966, с. 335-337.
[13] Кампсти Н. Аэродинамика компрессоров. Москва, Мир, 2000. 688 с.
[14] Aungier R.H. Centrifugal compressors: a stragedy for aerodynamic design and analysis. New York, ASME Press, 2000. 320 p.
[15] Нуждин А.С. Исследование течения в безлопаточных диффузорах. Автореф. дис. ... канд. техн. наук. Ленинград, ЛПИ, 1969. 19 с.
[16] Селезнев К.П., Галеркин Ю.Б., Нуждин А.С., Стрижак Л.Я., Суслина И.П. Повышение эффективности центробежных компрессорных ступеней с безлопаточными диффузорами. Компрессорные и вакуумные машины. Сб., ЦИНТИ-ХИМНЕФТЕМаш, 1968, № 2.
[17] Галеркин Ю.Б., Нуждин А.С., Селезнев К.П. Влияние форм профиля безлопаточного диффузора на эффективность работы центробежной компрессорной ступени. Исследования в области компрессорных машин. Тр. II конф. по компрессоростроению, Киев, Будивельник, 1970, с. 203-214.
[18] Галеркин Ю.Б., Рекстин А.Ф., Соловьева О.А. Выбор размеров безлопаточного диффузора центробежной компрессорной ступени на стадии первичного проектирования. Известия высших учебных заведений. Машиностроение, 2019, № 10, с. 43-57, doi: 10.18698/0536-1044-2019-10-43-57
[19] Галеркин Ю.Б., Соловьева О.А. Совершенствование методов расчета безлопаточных диффузоров центробежных компрессорных ступеней на основе вычислительных экспериментов. Часть 1. Компрессорная техника и пневматика, 2014, № 3, с. 35-41.
[20] Галеркин Ю.Б., Соловьева О.А. Совершенствование методов расчета безлопаточных диффузоров центробежных компрессорных ступеней на основе вычислительных экспериментов. Часть 2. Компрессорная техника и пневматика, 2014, № 4, с. 15-21.
[21] Rekstin A.F., Drozdov A.A., Solovyeva O.A., Galerkin Yu.B. Two mathematical models centrifugal compressor stage vaneless diffuser comparison. Oil and Gas Engineering, AIP Conf. Proc, 2007, Omsk, Russia, 2018, pp. 030035-1-030035-12, doi: https://doi.org/ 10.1063/1.5051896
[22] Галеркин Ю.Б. Вопросы проектирования проточной части центробежных компрессоров природного газа. Компрессорная техника и пневматика в XXI веке. Тр. Между-нар. науч.-техн. конф. по компрессоростроению, 15-17 сентября 2004, Украина, Сумы, Изд-во СумГУ, 2004, т. 2, с. 166-188.
[23] Sorokes J.M., Hutchinson B.R. The Practical Application of CFD in the Design of Industrial Centrifugal Compressors. Challenges and Goals in Pipeline Compressors, 2000, PID, vol. 5, pp. 47-54.
[24] Sorokes J.M., Nye D.A., D'Orsi N., Broberg R. Sidestream optimization through the use of computational fluid dynamics and model testing. Proceedings of the 29th Turbomachinery symposium, Texas, A&M, 2000, pp. 21-30, doi: 10.21423/R1QS8R
[25] Галеркин Ю.Б., Рекстин А.Ф., Солдатова К.В., Дроздов А.А. Альтернативный способ расчета характеристики коэффициента теоретического напора центробежного компрессорного колеса. Компрессорная техника и пневматика, 2016, № 6, с. 11-19.
[26] Галеркин Ю.Б., Солдатова К.В. Принципы моделирования напорной характеристики центробежного рабочего колеса. Компрессорная техника и пневматика, 2016, № 1, с. 24-34.
[27] Drozdov A., Rekstin A. Analysis of the velocity diagrams of impellers of centrifugal compressor stages after the preliminary design. International Scientific Conference on Energy, Environmental and Construction Engineering, MATEC Web Conf., 2018, vol. 245, doi: 10.1051/matecconf/201824504004
[28] Rekstin A.F., Galerkin Yu.B. The primary design method development of centrifugal compressor impellers based on the analysis of the geometrical parameters. Oil and Gas Engineering. AIP Conf. Proc., Omsk, Russia, 2019, vol. 2141, pp. 030052-030052-10, doi: https://doi.org/10.1063/L5122102
[29] Rekstin A.F., Galerkin Yu.B., Soldatova K.V. Computer programs application for development a primary design recommendations of low-flow rate centrifugal compressor stages. Oil and Gas Engineering. AIP Conf. Proc., Omsk, Russia, 2019, vol. 2141, pp. 030032030032-10, doi: https://doi.org/10.1063/L5122082
[30] Солдатова К.В. Создание новой математической модели проточной части центробежных компрессоров и базы данных модельных ступеней. Дис. ... д-ра техн. наук. Санкт-Петербург, СПбГПУ, 2017. 357 с.
References
[1] Ris V.F. Tsentrobezhnyye kompressornyye mashiny [Centrifugal compressor machines].
Leningrad, Mashinostroyeniye publ., 1964. 336 p.
[2] Ris V.F. Tsentrobezhnyye kompressornyye mashiny [Centrifugal compressor machines].
Leningrad, Mashinostroyeniye publ., 1981. 351 p.
[3] Den G.N. Mekhanika potoka v tsentrobezhnykh kompressorakh [The mechanics of the flow in
centrifugal compressors]. Leningrad, Mashinostroyeniye publ., 1973. 269 p.
[4] Den G.N. Proyektirovaniye protochnoy chasti tsentrobezhnykh kompressorov [Design of the
flow part of centrifugal compressors]. Leningrad, Mashinostroyeniye, Leningr. otd-niye publ., 1980. 230 p.
[5] Galerkin Yu.B. Turbokompressory. Rabochiy protsess, raschet i proyektirovaniye protochnoy
chasti [Turbochargers. Workflow, calculation and design of the flow section]. Moscow, KKhT publ., 2010. 596 p.
[6] Trudy nauchnoy shkoly kompressorostroyeniya SPbGTU [Proceedings of the scientific school
of compressor engineering SPbSTU]. Ed. Galerkin Yu.B. Moscow, KKhT publ., 2000. 443 p.
[7] Seleznev K.P., Galerkin Yu.B. Tsentrobezhnyye kompressory [Centrifugal compressor]. Lenin-
grad, Mashinostroyeniye publ., 1982. 271 p.
[8] Jansen W. Steady fluid flow in a radial vaneless diffuser. Transactions ASME. Journal of basic
engineering, 1964, Sept, pp. 607-619.
[9] Johnston J.P., Dean R.C. Losses in vaneless diffusers of centrifugal compressors and pumps.
Transactions ASME, Journal of Engineering for Power, 1966, vol. 88, pp. 49-60.
[10] Solov'yeva O.A., Galerkin Yu.B. The dimensions optimal ratios selection of low flowrate centrifugal compressor stage vaneless diffuser. Kompressornaya tekhnika i pnevmatika, 2017, no. 5, pp. 11-15 (in Russ.).
[11] Solov'yeva O.A., Galerkin Yu.B. Low flow rate centrifugal compressor stage vaneless diffuser analysis. Kompressornaya tekhnika i pnevmatika, 2017, no. 3, pp. 10-13 (in Russ.).
[12] Livshits S.P. Aerodinamika tsentrobezhnykh kompressornykh mashin [Aerodynamics of centrifugal compressor machines]. Moscow, Leningrad, Mashinostroyeniye publ., 1966, pp. 335-337.
[13] Cumpsty N. Compressor Aerodynamics. New York, Longman Scientific & Technical, 1990. 517 p. (Russ. Ed.: Kampsti N. Aerodinamika kompressorov. Moscow, Mir publ., 2000. 688 p.).
[14] Aungier R.H. Centrifugal compressors: a stragedy for aerodynamic design and analysis. New York, ASME Press, 2000. 320 p.
[15] Nuzhdin A.S. Issledovaniye techeniya v bezlopatochnykh diffuzorakh. Avtoref. Kand. Diss. [Investigation of flow in bladeless diffusers. Cand. Diss.]. Leningrad, 1969. 19 p.
[16] Seleznev K.P., Galerkin Yu.B., Nuzhdin A.S., Strizhak L.Ya., Suslina I.P. Povysheniye effek-tivnosti tsentrobezhnykh kompressornykh stupeney s bezlopatochnymi diffuzorami. Kom-pressornyye i vakuumnyye mashiny. Sb., TsINTI-KhIMNEFTEMash., 1968, no. 2.
[17] Galerkin Yu.B., Nuzhdin A.S., Seleznev K.P. Influence of the shape of the profile of a blade-less diffuser on the efficiency of the centrifugal compressor stage. Issledovaniya v oblasti kompressornykh mashin. Tr. II konf. po kompressorostroyeniyu [Research in the field of compressor machines. Proceedings of the II conference on compressor engineering]. Kiev, Budivel'nik publ., 1970, pp. 203-214.
[18] Galerkin Yu.B., Rekstin A.F., Solov'yeva O.A. Selecting the Dimensions of the Vaneless Diffuser of a Centrifugal Compressor Stage at the Primary Design Phase. Proceedings of Higher
Educational Institutions. Machine Building, 2019, no. 10, pp. 43-57 (in Russ.), doi: 10.18698/0536-1044-2019-10-43-57
[19] Galerkin Yu.B., Solov'yeva O.A. Improvement of vaneless diffuser calculations based on CFD experiment. Kompressornaya tekhnika ipnevmatika, 2014, no. 3, pp. 35-41 (in Russ.).
[20] Galerkin Yu.B., Solov'yeva O.A. Improvement of vaneless diffuser calculations based on CFD experiment. Part II. Kompressornaya tekhnika i pnevmatika, 2014, no. 4, pp. 15-21 (in Russ.).
[21] Rekstin A.F., Drozdov A.A., Solovyeva O.A., Galerkin Yu.B. Two mathematical models centrifugal compressor stage vaneless diffuser comparison. Oil and Gas Engineering, AIP Conf. Proc., 2007, Omsk, Russia, 2018, pp. 030035-1-030035-12, doi: https://doi.org/ 10.1063/1.5051896
[22] Galerkin Yu.B. Design issues of the flow part of natural gas centrifugal compressors. Kompressornaya tekhnika i pnevmatika v XXI veke. Tr. Mezhdunar. nauch.-tekhn. konf. po kom-pressorostroyeniyu [Compressor equipment and Pneumatics in the XXI century. Proceedings of the international scientific and technical conference on compressor engineering]. Ukraina, Sumy, 2004, vol. 2, pp. 166-188.
[23] Sorokes J.M., Hutchinson B.R. The Practical Application of CFD in the Design of Industrial Centrifugal Compressors. Challenges and Goals in Pipeline Compressors, 2000, PID, vol. 5, pp. 47-54.
[24] Sorokes J.M., Nye D.A., D'Orsi N., Broberg R. Sidestream optimization through the use of computational fluid dynamics and model testing. Proceedings of the 29th Turbomachinery symposium, Texas, A&M, 2000, pp. 21-30, doi: 10.21423/R1QS8R
[25] Galerkin Yu.B., Rekstin A.F., Soldatova K.V., Drozdov A.A. Alternative method of centrifugal compressor loading factor modeling. Kompressornaya tekhnika i pnevmatika, 2016, no. 6, pp. 11-19 (in Russ.).
[26] Galerkin Yu.B., Soldatova K.V. Loading factor performance of a centrifugal compressor impeller. Specific features and way of modeling. Kompressornaya tekhnika i pnevmatika, 2016, no. 1, pp. 24-34 (in Russ.).
[27] Drozdov A., Rekstin A. Analysis of the velocity diagrams of impellers of centrifugal compressor stages after the preliminary design. International Scientific Conference on Energy, Environmental and Construction Engineering, MATEC Web Conf., 2018, vol. 245, doi: 10.1051/matecconf/201824504004
[28] Rekstin A.F., Galerkin Yu.B. The primary design method development of centrifugal compressor impellers based on the analysis of the geometrical parameters. Oil and Gas Engineering. AIP Conf. Proc., Omsk, Russia, 2019, vol. 2141, pp. 030052-030052-10, doi: https://doi.org/10.1063/L5122102
[29] Rekstin A.F., Galerkin Yu.B., Soldatova K.V. Computer programs application for development a primary design recommendations of low-flow rate centrifugal compressor stages. Oil and Gas Engineering. AIP Conf. Proc., Omsk, Russia, 2019, vol. 2141, pp. 030032030032-10, doi: https://doi.org/10.1063/L5122082
[30] Soldatova K.V. Sozdaniye novoy matematicheskoy modeli protochnoy chasti tsentrobezhnykh kompressorov i bazy dannykh model'nykh stupeney. Dokt. Diss. [Creating a new mathematical model of the flow part of centrifugal compressors and a database of model stages. Doct. Diss.]. Sankt-Petersburg, SPbSTU publ., 2017. 357 p.
Статья поступила в редакцию 15.11.2020
Информация об авторах
СОЛОВЬЁВА Ольга Александровна — кандидат технических наук, старший преподаватель Высшей школы гидротехнического и энергетического строительства. Инженерно-строительный институт. ФГАОУ ВО «Санкт-Петербургский политехнический университет Петра Великого» (195251, Санкт-Петербург, Российская Федерация, ул. Политехническая, д. 29, e-mail: Solovyeva.OA@yandex.ru).
СОЛДАТОВА Кристина Валерьевна — доктор технических наук. Департамент информационных технологий и безопасности Мельбурнский политехнический институт (3072 Мельбурн, Австралия, Locked Bag 5 Preston Victoria, e-mail: buck02@list.ru).
ГАЛЕРКИН Юрий Борисович — доктор технических наук, профессор, главный инженер проекта Лидирующего исследовательского центра «Цифровое проектирование и моделирование (Smart Design)», Центр Национальной технологической инициативы Новые производственные технологии. ФГАОУ ВО «Санкт-Петербургский политехнический университет Петра Великого» (195251, Санкт-Петербург, Российская Федерация, ул. Политехническая, д. 29, e-mail: yuri_galerkin@mail.ru).
РЕКСТИН Алексей Феликсович — кандидат технических наук, ведущий научный сотрудник НИЛ «Газовая динамика турбомашин». Центр Национальной технологической инициативы «Новые производственные технологии». ФГАОУ ВО «Санкт-Петербургский политехнический университет Петра Великого» (195251, Санкт-Петербург, Российская Федерация, ул. Политехническая, д. 29, e-mail: rekstin2k7@mail.ru).
Information about the authors
SOLOVYEVA Olga Aleksandrovna — Candidate of Science (Eng.), Assistant Professor, Higher School of Hydraulic and Power Engineering Construction. Institute of Civil Engineering. Federal State Autonomous Educational Institution of Higher Education — Peter the Great St. Petersburg Polytechnic University (195251, Saint-Petersburg, Russian Federation, Politekhnicheskaya St., Bldg. 29, e-mail: Solovyeva.OA@yandex.ru).
SOLDATOVA Kristina Valerievna — Doctor of Science (Eng.). Melbourne Institute of Technology (3072 Melbourne, Australia, Locked Bag 5 Preston Victoria, e-mail: buck02@list.ru).
GALERKIN Yuriy Borisovich — Doctor of Science (Eng.), Professor, Chief Project Engineer, Leading Research Center "Digital Design and Modeling (Smart Design)", Center for the National Technology Initiative New Manufacturing Technologies. Federal State Autonomous Educational Institution of Higher Education — Peter the Great St. Petersburg Polytechnic University (195251, Saint-Petersburg, Russian Federation, Politekhnicheskaya St., Bldg. 29, e-mail: yuri_galerkin@mail.ru).
REKSTIN Aleksey Feliksovich — Candidate of Science (Eng.), Lead Research Fellow, Research Laboratory for Gas Dynamics of Turbomachines. Federal State Autonomous Educational Institution of Higher Education — Peter the Great St. Petersburg Polytechnic University (195251, Saint-Petersburg, Russian Federation, Politekhnicheskaya St., Bldg. 29, e-mail: rekstin2k7@mail.ru).
Просьба ссылаться на эту статью следующим образом:
Соловьёва О.А., Солдатова К.В., Галеркин Ю.Б., Рекстин А.Ф. Первичное проектирование безлопаточных диффузоров центробежных компрессорных ступеней Методом универсального моделирования. Известия высших учебных заведений. Машиностроение, 2021, № 3, с. 39-52, doi: 10.18698/0536-1044-2021-3-39-52
Please cite this article in English as: Solovyeva O.A., Soldatova K.V., Galerkin Yu.B., Rekstin A.F. Primary Design of Vaneless Diffusers of Centrifugal Compressor Stages by the Universal Modeling Method. BMSTU Journal of Mechanical Engineering, 2021, no. 3, pp. 39-52, doi: 10.18698/0536-1044-2021-3-39-52