Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Центробежные компрессоры

..pdf
Скачиваний:
14
Добавлен:
20.11.2023
Размер:
14.98 Mб
Скачать

няющих воедино Ф и ф , Попытки некоторых авторов представить такие коэффициенты как самостоятельные критерии, определяю­ щие свойства колес, следует оценить как нерациональное упро­ щение проблемы.

Широкое использование К ч или близкого по смыслу понятия удельной быстроходности как важного параметра колеса в практи­ ке насосостроения может быть оправдано постольку, поскольку у «насосных» колес величина ф обычно примерно одинакова, и тогда некоторое значение К п практически определяет важней­ ший параметр Ф.

Раздельное и очень сильное влияние Фр и фт. р на параметры колес, от которых зависит уровень потерь в самих колесах и не­ подвижных элементах, на форму напорной характеристики и размеры было показано выше. Колеса с одинаковыми /Сп, но разными Фр и фт. р могут оказаться попросту несопоставимыми.

По выбранному Фр с помощью соотношения (2.11) выбираются размеры входа в колесо. Втулочное отношение гвт и относитель­ ная толщина лопаток ô, влияющая на коэффициент стеснения ть задаются конструктивно из механических соображений. С аэро­ динамической точки зрения, чем они меньше, тем лучше, так как уменьшается величина w{. У колес с осевым входом на ло­

патки гвт не должно быть слишком малым,

соответствуя

rm/r0 ^

>

0,4-s-0,45, по аналогии с решетками осевых компрессоров.

S

Для РК ПЦК могут варьироваться величины KD>KF и

rlw»mm* Чем больше величина

тем

меньшими

потерями

сопровождается поворот потока из осевого направления в ра­ диальное. Но из-за увеличения ц>{ растет w{ и становится меньше удлинение лопаток Въ а значит, больше потери на ограничивающих поверхностях. Увеличение K D делает возможным^увеличить ра­ диус закругления покрывающего диска — в результате меньше потери на повороте, но создает и неблагоприятный ^эффект — увеличение rlw>min, а отсюда рост w[. При выбранных KF и K D принятие гх > rlw,mïn означает увеличение Вг (с одной стороны, это благоприятно, так как способствует росту удлинения лопатки; с другой — неблагоприятно, так как увеличивает неравномер­ ность потока во входном поворотном участке из-за увеличения размера канала по нормали). Одновременно уменьшается (возможно уменьшение тх), увеличивается угол наклона покры­ вающего диска 0, уменьшается фх.

При < rlüy'min параметры изменяются противоположно. Кон­ кретные рекомендации по выбору размеров входа для колес раз­ ных типов даны ниже, после обсуждения физической картины те­ чения и результатов экспериментов, так как сложность явления делает возможности аналитического исследования только каче­

ственными.

,

На выходе единственным геометрическим параметром в мери­

диональной плоскости является

В2. При заданном фт. р меньшие

В2 означают большие ср2, а 2, с2 (последнее означает большие потери

в НЭ при прочих равных условиях), меньшие гидравлические диа­ метры каналов НЭ, большие |32, большее отношение скоростей wf что улучшает условия течения в колесе. Приводимые ниже ре­ комендации по рациональным значениям w для колес и а 2 в за­ висимости от типа диффузора позволяют выбрать рациональные значения б2.

Заключительным этапом профилирования колеса с цилиндриче­ скими лопатками является определение формы лопаточной ре­ шетки на осесимметричной поверхности тока (ОПТ). За эту по­ верхность для упрощения считается допустимым принимать ра­ диальную плоскость. Ранее единственным и поныне широко рас­ пространенным приемом профилирования является чисто геоме­ трический: при определенных |3Л1 и |Зл2 средняя линия лопаток описывается дугами окружности, прямыми логарифмическими спиралями. Ниже обосновывается целесообразность применения аэродинамического подхода к профилированию, основанного на задании рационального распределения скоростей по поверхностям лопаток на ОПТ. Этот метод, разработанный на кафедре компрессоростроения Ленинградского политехнического института имени М. И. Калинина (ЛПИ), получил название метода ЛПИ. Ключе­ выми вопросами профилирования являются определение формы решетки, обеспечивающей заданный коэффициент теоретического напора, заданное положение оптимального режима, наивысший возможный к. п. д. и возможно более широкую зону устойчивой и экономичной работы колеса. Это осуществляется путем рациональ­ ного выбора густоты решетки, величины входных и выходных углов и формы средней линии лопатки в радиальной плоскости. Для конкретного колеса задача сводится к решению вопроса о применении одноили двухъярусной решетки, определению числа лопаток г2 и гъ зависимости рл = / (г) для первого и второго ярусов.

Для осерадиальных пространственных колес вместо профили­ рования в радиальной плоскости нужно спрофилировать решетку на некоторой средней (или периферийной) криволинейной кони­ ческой поверхности тока, а затем с принятием тех или иных законо­ мерностей построения квазитрехмерного потока и на ряде других поверхностей, что полностью описывает форму решетки. Вопросы профилирования колес типа «радиальная звезда» применительно к ПЦК и РК с пространственной решеткой изложены в п. 2.7.

2.3. Характер движения газа в рабочем колесе

Меридиональная плоскость. Для упрощения расчетов и ана­ лиза действительный трехмерный поток в РК часто представляют в виде квазитрехмерного (см. п. 1 .1), для чего вводится предпо­ ложение о бесконечности числа лопаток, средняя поверхность которых совпадает со средней поверхностью тока h — / в реаль­ ном межлопаточном канале. В силу осесимметричности поток

в таком РК становится двухмерным, не завися от угла Ô. Соот­ ветствующее распределение скоростей crn = / (г, z) в какой-то степени эквивалентно осредненному по 0 (т. е. по шагу t) меридио­ нальному потоку в реальном РК. Анализ осесимметричного ме­ ридионального потока полезен, так как способствует, по крайней мере, качественной оценке потерь и нахождению приемов рацио­ нального проектирования.

Запишем уравнение равновесия вдоль нормали в меридио­

нальной плоскости

при z ~

оо

[19; 96]

 

 

ÊJL

Pb +

- j - COSY- R

(2.24)

 

р дЪ

где Pb — проекция

на направление b условной силы

Р , отнесен­

ной к единице массы газа и заменяющей воздействие на поток разности давлений между передней и задней сторонами лопатки; R — радиус кривизны образующей ОПТ.

Сила Р определяется из условия равенства моментов dMA от­ носительно оси г, создаваемых разностью давлений Ар = рп — р3, и окружной составляющей силы Ри

dMz = zApdb dtmr = 2яг2т db dlmpPu = dfhd (cur),

(2.25)

где p — средняя

плотность в объеме dV =

2ягг db dlm.

const,

Принимая

в

пренебрежении

потерями,

р* = / (0) =

имеем Др/р =

0,5р3вд — 0,5рпо£,

считая средние плотности

в про­

цессах торможения р3 и рп от р до р* примерно равными р, по­ лучим

 

Ap/pæ Aw-w,

(2.26)

где Aw = w3 — доп; w = 0,5 (w3 +

wn).

 

Кроме того, приняв в (2.25) dm

— pcm2nrxdby получим другое

выражение (так как dr/dlm = sin у):

 

Ар/р =

ст (г |^ - + cHsinY^ .

(2.27)

Тогда из (2.25) величина

Ри может быть определена следующим

соотношением:

 

 

 

Ри = Ар/(р/т) = (kw -w yu •= ст

siny) =

= Cm ( 2(0 sin у + ^

- ctg p

ctg p sin y ) . (2.28)

Проекция силы P на нормаль b определяется через угол %ме­ жду b и касательной к поверхности лопатки в плоскости, перпен­ дикулярной к МП,

Рь = Р * Ч ъ

(2-29)

53

Так

как

 

p =

р* — 0,5рдо2,

р*

: p* +

0,5ри2,

при

hw — 0

w = Cm!sin р,

 

то p =

pi* +

0,5po)V -

0,5pCm/sin2 p и

 

 

 

1

dp

о

 

 

 

 

 

Ctn

àcm

(2 30)

 

-----== CoV COS V +

1 Sin2 P c t g P - f - -

sin2 p ~W

 

P

db

 

r

 

так как dr/db

 

— cos 7 .

постоянства

напора

по

высоте

лопатки

В частном

случае

hT = f (b) =

const

и при

culrt =

const должно выполняться ус­

ловие

сии =

(0сиг =

/ ф) = const,

тогда

из

p =

pï + р/гт —

— 0,5рс2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

др

си

 

 

дст

 

 

(2.31)

 

 

 

 

p

дЬ

— COS у Ст

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тригонометрическое соотношение между углами р и % в соот­ ветствии с [96] при принятом противоположном направлении b

 

+

<2 3 2 >

Приравнивая уравнения (2.24) и (2.30) и решая их относительно

dcjdb с учетом (2.32), получим

 

 

~ (sin2 p — cos2 р) -|- 2<о cos pcos у sin р —

 

- 2

cos2Рcos v - -gj. sin2p tg x - CmC° S t l i t

(2-33)

Выражение (2.33) с учетом (2.28) пригодно для численного оп­ ределения поля осредненных меридиональных скоростей в произ­ вольной вращающейся решетке. Оно также может служить для решения обратной задачи — получения необходимого распреде­ ления скоростей ст = f (b), например за счет выбора соответ­ ствующих значений угла %= f (b, ст). Одна из ОПТ при этом должна быть выбрана как исходная, т. е. для нее следует задать Р = / (1т)- Целесообразнее принимать за исходную ОПТ поверх­ ность покрывающего диска. В зависимости от формы РК в ме­ ридиональной плоскости, выбранных Р = / (1т) на исходной ОПТ и ст — f (b) получаются в соответствии с (2.32) вполне определен­

ные p = f (1т) для остальных ОПТ, т. е. те или иные W и распре­ деления напоров. Как будет показано ниже, для закрытых РК замедление потока на выпуклой ограничивающей поверхности не столь важно, как для неподвижного канала. Поэтому профили­ рование пространственных РК целесообразнее вести, задаваясь

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

ОПТ (например, задаваясь W),

а поле скоростей ст определять

по (2 .33).

Уравнение (2.33) несколько сложно для анализа, поэтому представим зависимость дст/дЬ в менее детализированной форме, используя те же соотношения (2.24) и (2 30)*

дст __

Сщ

|

1

à (rcu)

Ри

,

 

(2.34)

дЬ

R

^

г tg р " дЬ

ст 1&К

 

 

 

Из (2.34) следует,

что

и

при д (rcu)ldb =

0,

т.

е. при hT =

= f (b) = const, распределение

осредненных

по

0

скоростей ст

не соответствует распределению скоростей, имеющему место в ана­

логичном

осесимметричном

канале

без лопаток, где

dcjdb =

= сJ R .

Это следует из невозможности одновременно обеспечить

условия

hT = f (b) = const

и %=

0, так как углы р

(которые

определяют hr) и % связаны тригонометрическим соотношением (2.32).

Докажем сказанное анализом приведенных уравнений. При­

нимая

си =

и ст ctg р,

получаем

изменение

циркуляции сиг

вдоль

направления b

 

 

 

 

 

=

2»r cos V - С ctg PCOS т -

Г Ctg Р ^

^

% ■

При h, =

f (b) = const,

т. e. при

c j = / (b)

= const,

измене­

ние угла

потока по высоте лопатки в меридиональной плоскости

определяется так:

4f

- = i ^ ' ( rctgP'% L+ Стctg Р cos у — 2w cos Y) •

Если принять в некоторой точке отсутствие «навала» лопатки, т. е. %= 0, то из (2.34) в этой точке dcjdb = сJ R и, решая (2 .33) совместно с предыдущим уравнением относительно d % /d lm:

д% _

2(о cos у ■ , n

dim

Cm

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

ваются

осредненные по 0 меридиональные проекции скорости,

а в гл.

1 — действительные.

Анализ уравнения (2.34) применительно к РК ПЦК с радиально расположенными цилиндрическими лопатками показывает сле­ дующее. Последний член в уравнении (2.34) пренебрежимо мал, параллельные оси г образующие цилиндрических лопаток весьма близки к нормалям Ь, и углы %практически равны 0.

На расчетном режиме лопатка безударно обтекается на одной из средних ОПТ. В соответствии с рис. 2.4 скорости ст у пери­ ферийной ОПТ (кривые 2) значительно больше средней по b ве­ личины ст, а у втулки (кривые 1) меньше, т. е. соответственно там имеют место отрицательные и положительные углы атаки. В таком

случае периферийная часть лопатки обтекается в «турбинном» режиме, с уменьшением hT по длине траектории, пока углы р не станут примерно равными по всей высоте лопатки. Таким образом, на входном участке лопатки у периферии сиг меньше, чем у втулки, и д (cur)ldb < 0. Это значит, что по сравнению с каналом без лопаток в реальном колесе dcjdb должно быть мень­ ше, т. е. выравнивание меридиональных скоростей (cms cmh) (рис. 2.4) после входа потока в межлопаточные каналы должно происходить быстрее [см. уравнение (2.34)]. Впрочем, численный анализ показывает сравнительно слабую роль описанных эф­ фектов при обычных конфигурациях РК ПЦК.

Из сказанного можно сделать два вывода. Во-первых, харак­ тер течения на осесимметричных поверхностях тока практически не влияет на форму меридионального течения. Это позволяет изучать меридиональное течение изолированно от течения на осесимметричных поверхностях тока. Во-вторых, возможно теоре­ тическое и экспериментальное исследование меридионального те­ чения в РК ПЦК с помощью моделей в виде неподвижных кана­ лов. Однако при анализе результатов следует помнить о суще­ ственном влиянии вращения дисков колеса на состояние погранич­ ных слоев, меняющее картину течения вязкого газа в РК по сравне­ нию с неподвижным каналом. Эта проблема будет проанализиро­ вана ниже.

Рассмотрим картину невязкого меридионального течения. На рис. 2.4 показано распределение меридиональных скоростей по по­ верхностям основных 1 и покрывающих 2 дисков для трех вариан­ тов колес с различной формой покрывающего диска. В соответ­ ствии с известными законами движения газов в криволинейных каналах на выпуклой поверхности покрывающего диска в области поворота возникает пик скорости, величина которого при задан­ ном угле поворота зависит от кривизны канала. За безразмерный критерий кривизны по аналогии с криволинейными коленьями можно ориентировочно принять кривизну образующей ОПТ, средней между основным и покрывающим дисками:

RJb « (Rs + 0,5bj/bx = RJbx+ 0,5.

Ниже будет показано, что, поскольку в отличие от плоских коленьев речь идет об осесимметричных каналах, критерий кри­ визны, обобщающий количественные данные по распределению скоростей, должен учитывать и это обстоятельство.

Соответствующий «провал» скорости в области поворота на ос­ новном диске выражен значительно слабее. Пик скорости на по­ крывающем диске может отрицательно влиять на эффективность РК по следующим причинам.

1. Местная диффузорность после пика может привести к срыву потока и появлению соответствующих вихревых потерь.

2. Если входная кромка периферийного сечения лопатки рас­ положена в области, где меридиональная скорость значительно

превышает среднерасходную, £то может привести к снижению эф­ фективности на расчетном режиме из-за отрицательного угла атаки в периферийном сечении лопатки и ухудшить работу колеса

вправой ветви характеристики.

3.При достаточно высоких числах М на входе в решетку в периферийном сечении скорость может превзойти звуковую с соот­ ветствующим снижением к. п. д. и ухудшением зоны работы РК.

4.Различие в величине скоростей на основном и покрывающем

Рис. 2.4. Распределение меридиональных скоростей для

РК

ПЦК с различной формой покрывающего диска:

1 , 2

основной и покрывающий диски; 3 — средняя скорость

дисках увеличивает кинетическую энергию потока по сравнению со случаем, когда все частицы движутся со среднерасходной ско­ ростью. Это увеличивает все виды потерь, так как силы сопротив­ ления пропорциональны кинетической энергии.

Детальное изучение меридионального течения в каналах РК совместно с U-образными подводами различной формы было пред­ принято авторами совместно с В. И. Зыковым [64]. На рис. 2.5 показано распределение меридиональных скоростей на стенках одного из изученных каналов, полученное методом ЭГДА в сопо­ ставлении с измеренными скоростями при статической продувке этоцо канала. Видно, что на выпуклой поверхности канала воз­ никают зоны срыва, в доотрывной области расчетные и измерен-

нУе скорости хорошо совпадают. На рис. 2.6 показано поле ско­ ростей вязкого течения, демонстрирующее размеры и интенсивность срывной области. Коэффициент потерь, соответствующий мери­ диональному течению в неподвижном канале, оказался существен­ но зависящим от интенсивности срыва на выпуклой поверхности, которая определяется местной диффузорностью невязкого потока.

Рис. 2.5. Распределение скоростей на

 

 

 

 

стенках

переходного канала, имити­

в переходном канале

 

рующего

меридиональную форму РК

 

ПЦК и U-образного подвода:

 

 

 

 

— ------ методЭГДА,--------— — стати ie*

 

 

 

 

 

ские продукки

 

 

 

 

Экспериментальная зависимость £/£mln =

f (<У

на

рис. 2.7

под­

тверждает сказанное. Отметим, что при

cs >

0,6

срывные

зоны

в неподвижных каналах не возникали.

 

 

 

 

При изучении относительного потока в РК некоторые иссле­

дователи

наблюдали подобную зону срывов. Однако визуализация

течения с помощью напыления мелкодисперсного красителя в РК, имеющих современную меридиональную форму, не подтвердила наличие подобных срывов. На рис. 2.18 показана зона срыва на покрывающих дисках РК при разных режимах работы. Ни в од­ ном случае сплошной зоны срыва, характерной для неподвижных каналов, зарегистрировано не было. Срывные области возле по­ верхности лопаток связаны с ударным обтеканием периферийных сечений лопаточной решетки.

Метод визуализации напыления мелкодисперсного красителя был предложен В. П. Митрофановым и В. И. Зараевым. Порошко­ образный краситель вводится в проточную часть перед исследуе-

мым элементом. При соблюдении определенных условий краситель повторяет траекторию газовых частиц и прилипает к ограничива­ ющим каналы поверхностям в зонах малых касательных напряже­ ний, т. е. в областях срывов и утолщенных пограничных слоев. Наладочные опыты показали, что визуализационная картина устой­ чива. Она не зависит от количества введенного красителя (если оно не было слишком малым) и от времени работы ступени после его введения (при условии, что ступень работает на режиме и2

Рис. 2.7. Эксперименталь-

Рис. 2.8. Результаты

визуализаций

ная зависимость коэффи­

срывной зоны за

цилиндром

циента потерь

переходно­

 

 

го канала от местной диф-

 

 

фузорности

потока cs

 

 

= const, рк = const). На рис 2.8 показана картина обтекания цилиндра, полученная на пластине, перпендикулярной к его оси; она идентична иным получаемым другими способами визуализа­ ции. Условия подобия, соблюдение которых необходимо при ис­ пользовании метода (по В. Е. Маслову, В. Д. Лебедеву):

Fr = w2/gL ^ 50; St =

Ô^pT/(|xL)=^70;Re =

100,

где L — длина лопатки; ôT и рт — средний размер

и плотность

частиц красителя.

вращение дисков РК

препятствует

Как показывает анализ,

возникновению срывов на их поверхностях как в зоне осесимме­ тричного подвода, так и в самих межлопаточных каналах. На частицу газа, «прилипшую» к поверхности покрывающего диска,

в

направлении

вектора ст действует центробежная сила dPn =

=

cor2 cos у -dm,

препятствующая срыву. Сравним эту силу с си­

лой от градиента давления, возникающего из-за уменьшения ме­ ридиональной скорости:

йРц соУ COS у dm f cos у ое\

dPт

dV др/д1т

Сщ/СтоУо дст/д1щ

Здесь все величины, кроме относящихся к сечению 0—0 (рис. 1.1), следует рассматривать в вероятной точке срыва на поверхности покрывающего диска. Для типичных РК ПЦК в этом случае г *=» 0,6; cos у ^ 1; ф0 ^ 0,3, а при рекомендуемых Я8/Ьг + 0,5 ^ 1, как следует из расчетов невязкого течения, в окрестностях точки

отрыва

в

соответствующем

неподвижном канале ст/ст0 ^ 1 ,2 ;

dcjdlm

^

1,5—2. Отсюда следует, что соотношение между ста­

билизирующей центробежной

силой и ведущей к срыву силой от

встречного давления по (2.35) имеет порядок единицы. При мень­

ших R s/b1 + 0,5 градиент давления

получается большим, однако

и в этом случае стабилизирующая

роль центробежной силы ос-

Рис. 2.9. Визуализация течения в меридиональной плоскости вращающегося и неподвижного канала:

а — Re =

wblv = 370,

Ro = coblw =

0; б — Re =

= 370,

Ro = 0,105;

в — Re = 370,

Ro = 7,4

тается большой. Ее действие препятствует «вырождению» про­ филя скорости в пограничном слое (в МП). Остающиеся высокими касательные напряжения передают достаточную энергию от ча­ стиц ядра потока и препятствуют срыву.

Внутри межлопаточных каналов стабилизирующее действие на дисках могут создавать вторичные течения. Кориолисова сила, действующая на частицы пограничного слоя в направлении основ­ ного потока, равна 2 соwh dm. В результате по сравнению с течением в неподвижном канале профиль скорости в пограничном слое ока­ зывается более заполненным, возрастают касательные напряже­ ния на стенках, препятствующие срыву. На рис. 2.9 приведены экспериментальные данные Е. М. Смирнова и С. В. Юркина (кафедра гидроаэродинамики ЛПИ им. М. И. Калинина), демон­ стрирующие характер изменения скорости между ограничиваю­ щими стенками вращающегося и неподвижного канала в зави­ симости от Re и критерия Ro (Ro == сob/w). Через канал двига­ лась вода, профиль скорости получен фотографированием пу­ зырьков водорода, выделяющихся при подаче электрических импульсов на натянутую поперек сечения канала проволоку.

В межлопаточных каналах в принципе возможен срыв на огра­ ничивающих поверхностях, связанный с уменьшением относитель­ ной скорости. Однако такие срывы обычно не наблюдаются, что можно связать с описанным выше эффектом. Таким образом, вто­ ричные течения помимо рассмотренного ниже отрицательного влия*