Расчет и оценка эффективности холодильных центробежных компрессоров методом обобщенной политропы Текст научной статьи по специальности «Механика и машиностроение»

Раздел 1. ХОЛОД

УДК 621.565: 621.515

Расчет и оценка эффективности холодильных центробежных компрессоров методом обобщенной политропы

Д-р техн. наук Н.Н. БУХАРИН, А.А. ПОПОВ

СПбГУНиПТ

The article considers the most general and strict method of calculation of polytropic processes in uncooled refrigerating centrifugal compressors -the common method, when operating both on any real gas (or their mixture), or on ideal gas. The performed analysis gives the reasons to recommend in all the cases the determination of polytropic efficiency of compressors by the method of generalized polytrope, providing strict and stable results under all the initial and final conditions, not depending on the sphere of the constitution diagram, where the calculated processes take place. The method of generalized polytrope is useful as a basis of modern methods of calculations and evaluation of centrifugal compressors efficiency.

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

Применение уравнения процесса, справедливого для идеального газа,

/ю" = const, (1)

где р и и - текущие значения давления и удельного объема соответственно, при расчете процессов в реальных газах дает ощутимые погрешности даже при низких давлениях, что делает невозможным его прямое использование в таких случаях. Тем не менее относительная простота расчетов, основанных на использовании уравнения (1), оказалась настолько привлекательной, что до настоящего времени делаются попытки разработки уточненных методов расчета полит-ропных процессов в реальных газах, основанных на его использовании.

Одной из первых работ в этом направлении, оказывающей еще и сейчас значительное влияние на развитие методов расчета центробежных компрес-

соров, сжимающих реальные газы, была работа Дж. Шульца [4], который пришел к трем уравнениям с тремя разными показателями политропы в реальном газе вместо одного уравнения процесса (1) с одним показателем политропы в идеальном газе:

pv = const; = const;

Р_

Т

п-1 р"

(2)

- = const,

Z

где Т — текущее значение температуры.

Показатели степени в этих уравнениях определяются по зависимостям:

m : ТОЮ ,

h(pJpH) ’

_=М£к/£н).

Ш/ькУ m и —1_ lgOH/zK)

(3)

ЫРк/РнУ

где г — коэффициент сжимаемости;

индексы «н» и «к» соответствуют началу и концу процесса.

В реальном газе показатели (3) могут быть определены точно только в дифференциальной форме и в такой записи представляют собой осредненные значения. Расчет параметров трехступенчатого компрессора по результатам испытаний Дж. Шульц проводит по зависимостям:

1ё(р./ри).

/с 'еСРк/Рн).

•ё^н/Цс»)’

где к — цоказатель изоэнтропы;

к б — индекс конца изоэнтропного сжатия;

е

(4)

(5)

/ = -

^—[(Рк^кз-Рн^н)

к ^

и= |иф =/~т(РЛ - Р„Ч.).

^ П -1

н

где / — фактор политропного напора; /п0Л - политропная работа;

(6)

(7)

к

Ск ~ Сн

-|_______&___________Н.

Л,=Лпол= —

(8)

где — политропный КПД процесса;

/ — энтальпия рабочего вещества;

* — верхний индекс, обозначающий параметры торможения.

Удельные объемы ин, ик, ию и энтальпии I *, г*, 1к5 находят по результатам измерений с помощью диаграмм состояния, таблиц или непосредственно по уравнению состояния реального рабочего вещества и полученным на его основе зависимостям для определения калорических параметров. Введение в расчет «фактора политропного напора» / обусловлено расхождением значений изоэнтропной работы, определенной по формулам к

^ — £_ | ^Рч^КБ ~ Рн^н)» Ь ~ _ *Н >

и связано с тем, что значение показателя изоэнтропы к находилось по формуле (5). Возникающее при этом расхождение вызвано тем, что в реальном газе значение показателя изоэнтропы — переменная величина, а уравнение (5) дает некоторое ос-редненное значение к и, кроме того, это среднее значение относится к теоретическому изоэнтроп-ному процессу «н — кб», а не к действительному процессу «н — к», для которого, собственно, и определяется политропная работа (7).

В настоящее время в термогазодинамике турбомашин наиболее часто используется уравнение состояния идеального газа, которое лежит в основе почти всех теоретических и прикладных методов расчета. Поэтому желательно сохранить преемственность в использовании таких методов применительно к процессам, происходящим в реальных газах. Этим тре-

бованиям удовлетворяет метод условных температур, идея которого была высказана В. Траупелем, а впервые применена к расчетам холодильных центробежных машин Ф. М. Чистяковым [3]. В более завершенном виде метод сформулирован в работах [1,2].

Метод условных температур базируется на уравнении реального газа в виде ри = ?ЯТ, где произведение гТ, имеющее размерность температуры, рассматривается как некоторая условная температура Ту = zT. В этом случае уравнение состояния реального газа примет вид /ю = ЯТу, формально совпадающий с уравнением идеального газа. Главным достоинством метода является возможность выполнения всех расчетов процессов в реальных газах в том же виде, что и для идеального газа, включая поэлементные расчеты центробежных компрессов, что практически невозможно при использовании метода Шульца. Подробное описание метода дано в вышеупомянутых работах.

Общий недостаток обоих методов состоит в их ориентации на уравнение (1) политропного процесса, справедливого лишь для идеального газа.

Предполагаемый метод обобщенной политропы использует определение политропного процесса в реальном газе как процесса с постоянной теплоемкостью.

Количество теплоты, подведенной к 1 кг вещества в политропном процессе «н — к» (рис. 1),

?н-к=Спол(7’н-7’к), (9)

где спол - теплоемкость политропного процесса.

Приращение энтропии в политропном процессе

я,, - я = сл

,1п:

(10)

ЙГВ

и ^ у* ^ пол г^1

н к

Отсюда можно определить теплоемкость политропного процесса

5. -5,

С гтпп — '

г

1п^

(11)

и удельное количество теплоты, подведенной в процессе сжатия к рабочему веществу, равное потерянной работе 1т _ к:

<?н-*='гн-к=(5к(12)

1п^

Т*

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

*н-к='к-'н>0- (13)

Техническая работа политропного сжатия равна разности полной и потерянной работ:

14.

^пол.н-к = J = /нчс — /гннк > 0 ,

или

^пол.н-к ~ 1к *н (5к 5н)

Т -Т

т

1п——

Ти

(И)

(15)

Политропный КПД процесса сжатия определяется, как обычно:

(16)

Выражения (15) и (16), полученные в предположении постоянства теплоемкости политропного процесса, являются универсальными для любых рабочих веществ, включая идеальный газ как частный случаи. Эти формулы получены без каких-либо допущений о свойствах рабочего вещества, сжимаемого в компрессоре; все эти свойства учтены опосредованно через термические и калорические параметры, которые определяются для каждого рабочего вещества индивидуально в соответствии с уравнением состояния, описывающим его свойства. Это дает основание считать зависимости, полученные методом обобщенной политропы, наиболее общими, строгими и принять их за основу для сравнения при расчетах политропных процессов и оценке эффективности центробежных компрессоров.

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

Сравниваются группы уравнений обработки результатов поэлементного экспериментального ис-

следования и определение параметров рабочего колеса.

Измерены массовый расход С и статическое давление при выходе колесар2. Известны все параметры в предыдущих сечениях и геометрия колеса.

Метод обобщенной политропы:

*2 = *0 + ^0-2^

Р2 = ЯРг>Ь)\

ф1 = р/л:

1 + Рир + Ртр = Лф2г» <Р2«> Р2>;

ф2“ = «10+Рпр+Рч.) ;

п Ф2г+Ф2ц-Ф0

°-2_ го+Р^+Р,)*,’

•*2 =АРъЬ)>

^гО-2 =(52 — 5о)

т2-т0

1п туг/

1п Т2/Т0 Метод условных температур:

^у2 =^у0 + ^0-2 (^ук“^ун)5

Р2

Р2 =

ЯТу2

У2

Ф2г=- г

Р 2^2к«2

1 + Рпр + Ртр = Лф2г Ф2шР2>; (П С*к(Тук~Тун) Ш

ф2“ 0 + Рч> + Р,)и?'

П0_2 - 1-

„2 . ,л2 2

Ф 2г + Ф 2и Ф 0

2 О Рпр Р тр) Ф2и ’ 18Рг/Ро

а°-2 1§Гу2/Ту о’

1г0-2= ст0 г)^( Ту2 “ ^уо) I

о

0-2

О.

Первые шесть уравнений обеих групп представляют собой замкнутые системы, остальные три параметра вычисляются после их решения. „ Нами проведена сравнительная оценка результатов определения термодинамических параметров холодильных турбоагрегатов на основе использования приближенных методов Шульца, условных

температур и метода обобщенной политропы.

Термические и калорические параметры рабочего вещества находились с помощью уравнения состояния Боголюбова — Майера с вириальны-ми коэффициентами, в размерном

_ .<?.

* = —= 1 + ТТь,Л

или безразмерном виде

О)'

к 1 1=1 _/=<> х Здесь z — коэффициент сжимаемости; р - давление, Па; и - удельный объем, м3/кг;

Я — газовая постоянная, Дж/(кг К);

Ьу — коэффициенты уравнения; т = Т/Тщ- приведенная температура;

Гкр — критическая температура, К;

(о = р/рКр— приведенная плотность;

р - плотность, г/см3;

рад - критическая плотность, г/см3.

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

температур и метода Шульца). Сравнение проводилось относительно метода обобщенной политропы.

Результаты обработки результатов экспериментальных исследований холодильного центробежного компрессора типа К60-82-1 приведены в таблице. Компрессор поддерживает температуры кипения этилена в трех испарительных системах (—90, —70, —56 °С) и подает этилен в конденсатор.

Для I секции отличие в значениях КПД, найденных методами Шульца, и условных темпера-тур лГж не превышает 0,5 %. Они меньше полученных методом обобщенной политропы Г|пол примерно на 1 %.

Для II секции компрессора все три метода дают близкие результаты — отличие сравниваемых величин не превышает 0,7 %.

Для III секции влияние метода обработки опытных данных невелико и не превышает 0,5 %.

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

Газодинамические характеристики компрессора К60-82-1

I секция

Режим 1 2 3 4 5 6 7 8

3,971 4,091 4,168 4,259 4,293 4,048 4,03 4,226

^пол 0,6005 0,6155 0,6271 0,6148 0,5883 0,6291 0,6055 0,6078

лусп ІПОЛ 0,5904 0,6054 0,6171 0,6041 0,5764 0,6194 0,5952 0,5968

шульц ІПОЛ 0,5929 0,608 0,6198 0,6069 0,5792 0,622 0,5978 0,5995

II секция

1,62 1,567 1,607 1,594 1,587 1,604 1,573 -

\и 0,7321 0,74 0,7325 0,7155 0,6964 0,6559 0,6067 -

цусл 'пол 0,7312 0,7392 0,7317 0,7147 0,6955 0,6549 0,6056 -

шульц ІПОЛ 0,7314 0,7393 0,7321 0,7151 0,6956 0,6553 0,6061 -

III секция

я к 2,424 2,348 2,585 2,555 2,598 2,580 2,572 -

^пол 0,7083 0,6863 0,6953 0,7172 0,7129 0,6997 0,7068 -

лусл шел 0,7056 0,6837 0,6921 0,7142 0,7097 0,6965 0,7036 -

„шульц ІПОЛ 0,7069 0,6847 0,6936 0,7155 0,7114 0,6979 0,7052 -

Рис. 2. Процессы сжатия в различных областях диаграммы состояния

дились значения политропных КПД процессов сжатия, проходящих в различных областях диаграмм состояния (рис. 2).

Расчеты проводились в областях низких, средних и высоких давлений при одном и том же отношении давлений. Начальные точки процессов сжатия выбирались в непосредственной близости от линии насыщенного пара, где свойства пара в наибольшей мере отличаются от свойств идеального газа. Исследовались процессы сжатия в рабочих веществах холодильной техники и химической и нефтегазовой промышленности: Ш2; 1122; аммиак (И.717); пропан (11290); этилен (Ю150); метан (Я50).

Для каждой точки 1 начала сжатия (см. рис. 2) рассчитывалось несколько процессов сжатия с постепенно возрастающими перепадами температур в точках, обозначенных 2^\ где у — текущий номер процесса. Это позволило, рассмотрев процессы сжатия с различными КПД (в пределах г|пол « » 0,8...0,5), установить, в какой мере влияет на величину погрешности отклонение рассчитываемого процесса от изоэнтропного.

Погрешности в определении удельных работ и КПД различными способами находятся относительно значений, вычисленных методом обобщенной политропы.

Для Ш 2 на режиме низких давлений (рн = 0,4 бар, рк = 4 бар) (рис.З) погрешности определения по-литропного КПД обоими методами достигают наименьших значений при наибольшем КПД и возрастают примерно на 0,8 % при наименьшем КПД. Так, при расчетах методом условных температур погрешности возрастают с 1 до 1,89 %. При расчетах методом Шульца погрешности меньше на ~ 1 % во всем диапазоне и также возрастают от 0,06 до

0,95 % с уменьшением КПД. На режиме средних давлений (рп = 1 бар, рк = 10 бар) погрешности определения политропного КПД изменили свой ха-

Рис. 3. Зависимость относительной погрешности от т| (Я12):

1 — метод относительных температур; 2 — метод Шульца

рактер. При расчетах методом условных температур они положительны и возрастают с уменьшением КПД от 1,38 до 1,63 %. При расчете методом Шульца погрешности стали отрицательными, т.е. значения КПД, рассчитанные методом Шульца, стали меньше значений, полученных по базовым зависимостям обобщенной политропы. При этом характер кривой сохранился таким же, как и в случае низких давлений, но из-за изменения знака наибольшая погрешность (—0,36 %) наблюдается при максимальном КПД и постепенно уменьшается до —0,11 % при минимальном. На режиме рн = 4 бар, рк = 40 бар погрешности определения КПД возросли до 1,62...3,07 % для метода условных температур и до -1,24...-2,76 % для метода Шульца.

Для Л22 на режиме низких давлений (рн = 0,5 бар, рк = 5 бар) относительные погрешности определения политропного КПД обоими методами не превышают 0,8% (рис. 4), достигая наименьших значений при низком КПД, и возрастают с его увели-

Рис. 4. Зависимость относительной погрешности от г| (1122): 1 — метод условных температур; 2 — метод Шульца

0,45 0,5 0,55 0,6 0,65 0,7 0,75 0,8 0,85 7

Рис. 5. Зависимость относительной погрешности от л (111150): 1 - метод условных температур; 2 — метод Шульца

чением. При расчетах методом условных температур погрешности слабо зависят от характера процессов сжатия и изменяются в пределах 0,6...0,8 %. При расчетах методом Шульца погрешности практически такие же во всем диапазоне. На режиме рн= 5 бар, рк = 50 бар погрешности определения политропного КПД практически такие же, как и на режиме низких давлений. На режиме высоких давлений (рн = 20 бар, рк = 200 бар) погрешности определения политропного КПД выше, чем в первых двух случаях. При расчетах методом условных температур погрешности составляют 1,3... 1,5 %. При расчетах методом Шульца погрешности изменяют знак (-0,6...-0,3 %).

Для этилена К1150 на режиме низких давлений (рн = 0,5 бар, рк = 5 бар) погрешности определения политропного КПД при расчетах методом условных температур равномерно возрастают с повышением перепада температур и составляют 1,53...3,47 % (рис. 5). При расчетах методом Шульца погрешности меньше на ~ 1 % во всем диапазоне и также возрастают от 0,62 до 2,6 %. На режиме рн = 5 бар,

рк = 50 бар погрешности определения политропного КПД при расчетах методом условных температур составляют 3,1...4,8 %. При расчете методом Шульца они равны —1,5...—3,2 %. На режиме рн -= 50 бар, рк = 400 бар погрешности определения политропного КПД при расчетах методом условных температур составляют —8...—15 %, уменьшаясь по модулю с ростом перепада температур и падением г|пол. При расчетах методом Шульца погрешности меньше и возрастают в пределах 1,9...6 %, увеличиваясь с падением г|пол.

Для пропана (11290) на режиме низких давлений (рн = 0,5 бар, рк = 5 бар) погрешности определения политропного КПД (рис. 6) при расчетах методом условных температур равномерно возрастают с повышением перепада температур и составляют 1,41...3,21 %. При расчетах методом Шульца погрешности меньше на - 1 % во всем диапазоне и также возрастают от 0,22 до 2%. На режиме рн = 5 бар, рк = 50 бар погрешности определения политропного КПД при расчетах методом условных температур составляют 4,7...3,27 %. При расчете методом Шульца погрешности отрицательны и равны —1,03...—2,56 %. На режиме рн = 22 бар, рк = 220 бар погрешности определения политропного КПД выше, чем на первых двух режимах. При расчетах методом условных температур они составляют —15,35...-8,06 %, увеличиваясь по модулю с ростом г|пол. При расчетах методом Шульца погрешности меньше, возрастают в пределах 1,22...7,47 %, увеличиваясь с падением т^пол.

Для аммиака Я717 на режиме низких давлений (рн = 0,5 бар, рк = 5 бар) погрешности определения политропного КПД (рис.7) при расчетах методом условных температур составляют 1,06...1,82 %. При

6

1,8

1.6

1,4

1,2

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

0,65 0,67 0,69 0,71 0,73 0,75 0,77 0,79 0,81 г\

Рис. 6. Зависимость относительной погрешности от г| (1(290): Рис. 7. Зависимость относительной погрешности от т| (1(717):

1 — метод условных температур; 2 — метод Шульца 1 — метод условных температур; 2 — метод Шульца

Рис. Я. Зависимость относительной погрешности от т) (А50):

1 — метод условных температур; 2 — метод Шульца

расчетах методом Шульца погрешности меньше на -0,5 % во всем диапазоне и также возрастают от

0,5 до 1,26 %. На режиме р„ = 5 бар, рк = 50 бар погрешности определения политропного КПД при расчетах методом условных температур составляют 1,05...1,67 %. При расчете методом Шульца погрешности меньше (0,35...0,98 %). На режиме рИ= 10 бар,рк = 100 бар погрешности определения политропного КПД при расчетах методом условных температур составляют 0,73...1,28 %, увеличиваясь с ростом перепада температур и пацением "Ппол- При расчетах методом Шульца погрешности меньше, возрастают в пределах 0,22...0,77 % и также увеличиваются с падением г|пол.

Для метана (1150) на режиме низких давлений (рИ = 0,7 бар, рк = 7 бар) погрешности определения политропного КПД (рис. 8) обоими рассматриваемыми методами практически одинаковы. При расчетах методом условных температур они равномерно возрастают с повышением перепада температур и составляют 0,02...0,43 %. При расчетах методом Шульца погрешности также возрастают от 0 до 0,41 %. На режиме р„ = 7 бар, рк = 70 бар погрешности определения политропного КПД возросли, но незначительно. При расчетах методом условных температур они составляют -0,74...+0,46 %. Смена знака погрешности означает ее переход через ноль, когда погрешность вовсе отсутствует. При расчете методом Шульца погрешности положительны и составляют 0,31...1,49 %. На режиме рн = 70 бар, рк = 700 бар погрешности определения политропного КПД при расчетах методом условных температур составляют —3,14...—1,13 %, уменьшаясь по модулю с ростом перепада температур и падением Ппол- При расчетах методом Шульца погрешности возрастают в пределах 2,05...3,96 %, увеличиваясь с падением г|пол.

Обобщая результаты сопоставления погрешностей в определении КПД процессов сжатия, прохо-

дящих в разных областях диаграмм состояния, можно сделать следующие выводы.

1. Практически для всех изученных процессов политропный КПД, вычисленный приближенными методами Шульца и условных температур, отличается от точного значения, найденного методом обобщенной политропы, в большей или меньшей мере.

2. По мере увеличения давлений погрешности возрастают, но для разных веществ в разной мере. При очень высоких сверхкритических давлениях в конце процесса сжатия (200...400 бар) у пропана Я290, этилена Ю150 и Я22 погрешности для метода Шульца при высоких КПД составляют 1,24...2,05 %,апри низких КПД равны 2,76...7,47 %. При этом метод Шульца в ряде случаев дает более точные результаты, чем метод условных температур.

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

4. Применение приближенных методов следует ограничить: метода условных температур — низкими давлениями, а метода Шульца - низкими и средними давлениями, и в обоих случаях по возможности невысокими отношениями давлений. Правомерность использования любого из этих методов всегда необходимо проверять выполнением контрольных расчетов методом обобщенной политропы.

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

Список литературы

1. Епифанова В. И. Низкотемпературные радиальные турбодетандеры. — М.: Машиностроение, 1974.

2. Холодильные машины / А.В. Бараненко, Н.Н. Бухарин, В.И. Пекарев, Л.С. Тимофеевский; под ред. Л.С. Тимофеевского.- СПб.: Политехника, 2006.

3. Чистяков Ф.М. Холодильные турбоагрегаты,— М.: Машиностроение, 1967.

4. ШулъцДж. Подтропический анализ центробежного компрессора// Тр. Американского общества инженеров-механиков. Энергетическое машиностроение, 1962, № 1.