Выходные устройства центробежного компрессора. Потери в них и их математическое описание Текст научной статьи по специальности «Физика»

Ю. А. Фирсова, И. Г. Хисамеев

ВЫХОДНЫЕ УСТРОЙСТВА ЦЕНТРОБЕЖНОГО КОМПРЕССОРА.

ПОТЕРИ В НИХ И ИХ МАТЕМАТИЧЕСКОЕ ОПИСАНИЕ

Ключевые слова: математическое моделирование, центробежные компрессоры; концевая ступень, выходное устройство, потери в проточной части.

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

Keywords: mathematical modelling, centrifugal compressors, a trailer step, the target device,

losses in a flowing part.

The method of mathematical modelling is developed for calculation of losses in target devices of centrifugal compressors. Settlement formulas are resulted. Results of calculation of factors of losses in the target device on mathematical model and other techniques are compared.

Выходные устройства (ВУ) концевых ступеней (КС) центробежных компрессоров (ЦК) оказывают существенное влияние на КПД и зону устойчивой работы, а также надежность ЦК, поскольку в значительной степени определяют аэродинамические нагрузки, действующие на ротор, особенно на нерасчетных режимах работы. Влияние ВУ резко усиливается при использовании центробежных компрессоров в области высоких давлений и при сжатии газов из группы галоидозамещенных углеводородов, работа на которых происходит при больших числах Маха, а также при повышении уровня окружных скоростей. В связи с этим отработка высокоэффективных конструкций ВУ и создание методов их расчета и оптимального проектирования имеют большое значение [1].

Метод математического моделирования для расчета потерь в ВУ разработан на основе комплексной обработки большого числа экспериментальных данных, полученных при испытаниях ступеней с различными геометрическими и режимными параметрами [2]. В данной работе в качестве выходного устройства ЦК рассмотрена кольцевая сборная камера (КСК), которая является наиболее распространенной в современном машиностроении [3].

Основываясь на представлениях о характере течения потока, проточную часть (ПЧ) выходного устройства можно разбить на три участка (рис. 1):

• участок 4-5 - поворотный канал;

• участок 5-6 - торовая часть камеры;

• участок 6-К - выходной конический диффузорный патрубок.

На выделенных участках возникают потери, обусловленные радиальным подводом массы газа (hp.n); потери трения (hTp); потери, обусловленные наличием радиального градиента давления (hMH,q); потери в нагнетательном патрубке, обусловленные трением, изогнутостью канала и диффузорностью течения (hH.n), а также потери в меридиональном сечении (hMep). Тогда потери в КСК

hKCK _ hTp + hMep + hMHA + hH.n- (1)

Слагаемые выражения (1) можно представить в виде структурно идентичных формул, имеющих вид (например, для Итр)

°с

Чр " Сwтр рср 2“~т РКСК ’ (2)

где сwтр - коэффициент силы сопротивления; рСр - средняя плотность газа в проточной части КСК; ССр - средняя скорость продольного обтекания канала КСК в радиальной плоскости; т - массовый расход газа; Ркск - площадь камеры.

Рис. 1 - Схема сечения выходного устройства в меридиональной (а) и радиальной (б) плоскостях

Используя аналогичную структурную запись для остальных составляющих потерь, можно получить выражение для определения суммарного коэффициента сопротивления %к к в виде суммы составляющих коэффициентов Сщ . В качестве параметра, характеризующего эффективность КСК, выбрана доля потерь КПД

Л Ькск С3 Р £

П _ 1^1 "2 (1 + Ртр + Рпр) фт Ф С^СК ’ (3)

F Р c

Р = 1 КСК . £ = К Ср . — = Ср

П й2 / 4. р0 . и2 ’

где Ртр - коэффициент дискового трения; Рпр - коэффициент протечек в колесе;

фт- коэффициент теоретического напора; ; Ф - коэффициент расхода; й2 - диаметр колеса;

Ро - плотность газа на входе в ступень ЦК; и2 - окружная скорость колеса.

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

плексов параметров, характеризующих распределение скоростей идеального потока w и форму КСК Р :

с^ =ф|„Ре„*„Я)[ (4)

Тогда общий вид аналитической зависимости потерь в выходном устройстве имеет

вид

1 £ Р С3 Г

Ап = ,—~^ ~---------------:т (км кРе сш (w, р) +

1 1 + р+ р 2 ф ф мтр тр ттр

+ ЧЕр К„еИЕр с 1«мЕ. км,„ КИецнд с Иинд ((',Р) +

+ КМн,п, К«е„д СИ„Д (^р)} (5)

где км, кре - поправочные коэффициенты, учитывающие влияние критериев соответственно М и Ре; Си - коэффициенты силы сопротивления соответствующих потерь.

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

1. коэффициент силы сопротивления трению СИ является функцией относительного

замедления потока - Ас/сСр;

2. оценка влияния пространственности течения газа в канале КСК производится путем введения соответствующих потерь Имер и ИИнД в выражение общего вида потерь;

3. коэффициент силы сопротивления СИинд является функцией поперечного градиента

скорости и относительной ширины канала;

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

Составляющие потерь, входящие в выражение (1), на основании принятых допущений определяются следующими выражениями:

- потери на трение [4]

£ Р с3.............................................Ас

'ср

итр = ----------ф------- х, ( 1 + Х2 ( — ) )Хз; (6)

потери в меридиональном сечении

1\ер = Х4*91,25 £ ( 1 - Ь4 Р4 )2 (Сг24) (7)

ср ср

индуктивные потери

и _ £ Р с 3 х . (сасР). . (8)

Инд ф 5 ( /I ). ( )

\т ср ' 1 ср /

5

'ср ' 'ср

- потери в нагнетательном патрубке

Ин„ =|св2к [(1 - 1/п2) ^э/2) + Хв 1дХд(/2) (1 - 1/п)2 ], (9)

где £ - отношение средней плотности воздуха сечения к плотности воздуха на входе ступень; Х1 — Хд - искомые параметры; Уэ - эквивалентный угол раскрытия; Ь4 - ширина се-

чения на входе в ВУ; Г4 - радиус сечения на входе в ВУ; ЬСр - средняя ширина сечения; ЇЧСр

- средний радиус сечения; Сг4 - среднерасходная радиальная составляющая скорости на входе в ВУ; СаСр - усредненный поперечный градиент скорости; 1Ср - усредненная длина ПЧ ВУ; С0к - скорость в сечении 9 = 9п; 9п - угол расположения нагнетательного патрубка;

П - коэффициент геометрической диффузорности.

В выражении (9) первая составляющая указывает потери на трение, вторая - потери от расширения.

На основе выражения (5) с учетом уравнений (6) - (9) получена в общем виде математическая модель для расчета потерь в КСК

1 1

ЛП = 1 + в в (Хіо МХіі) (і + Х^еХіз)

1 + Ртр + Рпр 2 Фт

X

£ F C X! ( 1 + X2 ^c) + X4 tg125 V, ( 1 ^i)2 (cr24>

Ф 14 2 ccp ' 2 ' bcp Rcp

+£ F c3 Х5 ( (cacp+ c2 (1- 1/n2) . /X\„\ + X8 tgX9(3/2) (1 1/n)2 , (10)

Ф 5 (bcp/lcp)Xe sin(v3/2) 8 a Va ’

Для нахождения конкретного вида ММ запишем уравнение (10) в виде

АП = (х, А), (11)

где АП - расчетное значение функции, в данном случае - доля потерь КПД, характеризующая эффективность КСК; x = (X1, X2, X3...XI, Xn) - вектор искомых параметров;

A = (a1,a2a3...aI,am) - вектор известных параметров, т.е. множество значений чисел М, Re, F и др.

На невязки, т.е. разницу между вычисленным и известным опытным значением Ап, накладывается условие, которое требует минимизации суммы квадратов невязок:

к 2

G = 1Кээк - (Хм Аi)] , (12)

1=1

где Ап экс - данные, полученные экспериментальным путем; (хI, AI) - значение функции,

полученное из расчетов по ММ; i - номер варианта или режима экспериментального испытания.

Задача решается перебором всех значений вектора х с некоторым шагом, обеспечивающим практическую точность ММ, и выбором из этого множества значений конкретного значения х *. В данной работе рабочим вариантом был принят метод случайного поиска с переменным шагом [2].

При определении неизвестных коэффициентов Х| методом случайного поиска границы исходных параметров варьировались в широком диапазоне. Конкретизация ММ потерь с использованием данных по 135 вариантам концевых ступеней производилась поэтапно с помощью программы, составленной на языке Visual Basic. Каждый этап поиска заканчивался анализом полученного решения, по результатам которого вносились необходимые коррективы в частные функциональные зависимости и в границы поиска вектора неизвестных коэффициентов. Среднее отклонение значений АПрасч от Апэкс было уменьшено с 18 до 5,7%.

По результатам проведенной идентификации ММ для определения Ап имеет окончательный вид:

Ап = п 1 п 7Г1— ( +15 М11) ( 1+58 )

1 + РТР + Рпр 2 уТ Яе

є Б с3 Ас

0,6

0,007( 1 + 0,03 — ) +

Ф сср

1,25 Уэ , л Ь4 Г4 \2 (-2 ) , Є Б С П ПППО , ч~асР

+ 0,12 tg1,25-2 ( 1 )2 (с24)+є ф с 0,0002

2 Ьср Кср ф

(саср )

(Ьср/іср)5,

1+

2 1 0,0018 / ч 1 2 /10Ч

Се2к ( 1 -2 ) Л . 2) + 8 tg(vэ /2) ( 1 - )2 . (13)

П 81ЩV э / 2 п

+Е*- (1 - ) ^(72)+8 ^э/2 (1 - )2

Диапазон применимости математической модели для расчета потерь в выходном устройстве: по числам Ре = 2,Ы05...4,2^105 и М = 0,3...0,4; по относительной площади камер Р = 1,5...4,4.

На рис. 2 представлены зависимости между составляющими потерь в ПЧ КСК, вычисленные по выражению (13). Как видно, потери трения составляют 10-15% суммарных потерь. Потери, обусловленные особенностями меридионального течения в КСК, составляют 12-17% от общих потерь; индуктивные потери, обусловленные наличием поперечного градиента давления, - 20-35%. Сложный характер картины течения в начальной зоне нагнетательного патрубка обусловливает сравнительно высокие значения потерь (40-50% суммарных потерь).

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

Рис. 2 - Соотношение между коэффициентами сопротивления составляющих потерь

Сопоставление расчетных и экспериментальных значений величин Дп для трех концевых ступеней (КС-52, КС-54, КС-57) (рис. 3) показало, что в пределах области применимости ММ, установленной в процессе экспериментального исследования, и идентификации модели, расхождение между теоретическими и экспериментальными значениями не превышает экспериментальной погрешности определения величины и составляет 9%.

■’ААг' >>,■

1,5

І Ч--------1--------1--------1-------1

0,6 0,8 I 1,2 1,4

ф/фр™.

♦ КС-52 ■ КС-54 * КС-57 Рис. 3 - Сравнение расчетных и экспериментальных значений Дп

Сопоставление характеристик, рассчитанных с использованием выражения (13) с характеристиками, полученными по рекомендациям других авторов (рис. 4) ([5], [6]), показывает, что данные, полученные расчетом по разработанной ММ [7], имеют большую точность.. Следовательно, математическая модель пригодна для анализа влияния определяющих параметров на потери в выходном устройстве центробежного компрессора.

1 2

• 3

0,01 0,02 0,03 0,04 - 0,05

Рис. 4 - Сравнение результатов расчета коэффициента потерь в КСК с использованием разработанной ММ и других методик: 1 - разработанная ММ; 2 - [5]; 3 - [6]

Литература

1. Мифтахов, А.А. Аэродинамика выходных устройств турбокомпрессоров/ А.А. Мифтахов. - М.: Машиностроение, 1999. - 480с.

2. 2 .Селезнев, К.П. Оценка эффективности двухзвенных ступеней центробежных компрессоров на основе статистической обработки результатов эксперимента / К.П. Селезнев, Ю.Б. Галеркин., А.Г. Никифоров // Труды ЛПИ. - СПб.: 1977. - С. 300-305

3. Максимов, В.А. Компрессорное и холодильное машиностроение на современном этапе/ В.А. Максимов, А.А. Мифтахов, И.Г. Хисамеев // Вестник Казан. технол. ун-та. - 1998. - №1. - С. 104113.

4. Идельчик, И.Е. Справочник по гидравлическим сопротивлениям /И.Е. Идельчик..- М.: Машиностроение, 1992. - 560с.

5. Никитин, А.А. Расчет потерь в выходном устройстве центробежного компрессора/ А.А.. Ники-тин,С.В. Цукерман // Энергомашиностроение. - 1979. - № 6. - С. 17-19.

6. Тарабарин, О.И. Исследование кольцевых сборных камер с различным расположением и числом нагнетательных патрубков / О.И. Тарабарин, А.А. Мифтахов // Тезисы докладов VI Всесоюзной научно-технической конференции по компрессоростроению. - Л., 1981.

7. Фирсова, Ю.А.. Совершенствование методов расчета кольцевых сборных камер центробежных компрессоров: дис. ... канд. техн. наук: 05.04.06: защищена 29.05.2009: утв. 18.09.2009 / Юлия Александровна Фирсова. - Казань, - 173 с.

© Ю. А. Фирсова - канд. техн. наук, доц. каф. холодильной техники и технологии КГТУ, firsovay@mail.ru; И. Г. Хисамеев - д-р техн. наук, проф., зав. каф. холодильной техники и технологии КГТУ.