Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора Текст научной статьи по специальности «Механика и машиностроение»

УДК 621.514

Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора

Д-р техн. наук, проф. Носков А.Н. dn.noskoff@rambler.ru Зверев Д.И kayff007@mail.ru, Тарасенков Д.С. d.tarasenkov@mail.ru Университет ИТМО Институт холода и биотехнологий 191002, Санкт-Петербург, ул. Ломоносова, 9

Приведена методика расчета процесса нагнетания винтового компрессора, учитывающая особенности холодильного маслозаполненного компрессора. Рассматривается два этапа процесса нагнетания при давлении внутреннего сжатия ниже давления нагнетания. Определены зависимости для расчета давления в парной полости винтового компрессора на этих этапах. Рассмотрена работа компрессора в составе среднетемпературной холодильной машины на хладагентах R717 и R404а, а в составе нижней ветви каскадной холодильной машины на хладагентах R23 и R744. Приведены зависимости индикаторного КПД компрессора, относительного сопротивления в процессе нагнетания, относительной потери работы компрессора в процессе нагнетания, относительной потери работы компрессора в процессе натекания хладагента в парную полость из камеры нагнетания от геометрической степени сжатия. В режиме среднетемпературной холодильной машины энергетическая эффективность работы винтового компрессора на хладагенте R717 выше, чем на R404а, а составе нижней ветви каскадной холодильной машины энергетическая эффективность работы на хладагенте R744 выше, чем на К23 при всех геометрических степенях сжатия.

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

Calculation of expulsion process of wet screw refrigeration compressor

Noskov A.N. dn.noskoff@rambler.ru Zverev D.I. kayff007@mail.ru, Tarasenkov D.S. d.tarasenkov@mail.ru University of ITMO.

Institute of Refrigeration and Biotechnologies 9, Lomonosov Street, St Petersburg, 191002

Calculation method for the pressurization of screw compressor, which considers features of oil-filled refrigeration compressor is given. Two stages of process when internal pressure is lower then discharge pressure is considered. Depending for calculation pressure in working space of these stages is identified. Considered compressor operation composed of medium-temperature chiller at refrigerant R717 and R404a, and composed of lower branch cascade chiller at refrigerant R23 and R744. Provides efficiency indicators, relative resistance in discharge process, relative energy loss in inflow process from discharge port depending on volume ratio. In medium-temperature chiller mode energy efficiency of screw compressor on refrigerant

R717 above then on R404a, and in lower branch cascade chiller on refrigerant R744 refrigerant R23 of all volume ratios.

Key words: wet screw refrigeration compressor, discharge process calculation, energy loss of compressor in discharge process

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

Процесс сжатия в холодильном винтовом маслозаполненном компрессоре (ВКМ) протекает в два этапа. Первый этап происходит с момента отсечения парной полости от окна всасывания до момента соединения этой полости с окном нагнетания. Второй этап происходит с момента соединения парной полости с окном нагнетания, до момента достижения давлением в этой полости давления в камере нагнетания. Далее следует процесс нагнетания (выталкивания), который продолжается до полного сокращения парной полости.

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

ри [1]. Режим работы ВКМ, при котором ра < ри схематично изображен на рис. 1.

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

давлений в парной полости и в окне нагнетания - (рв. Затем происходит выталкивание рабочего вещества из парной полости в камеру нагнетания. Этот процесс происходит на угле поворота ВЩ винта от <рв до (рс.

Рис. 1. Схема процесса сжатия и нагнетания

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

Метод определения изменения параметров рабочего вещества ВКМ на первом этапе процесса сжатия приведен в работе [2].

Рассмотрим второй этап процесса сжатия в ВКМ.

Масса пара, поступившего в парную полость из камеры нагнетания за рассматриваемый отрезок времени

сЮ = и ■с ■ Т7 (ф,) • р ■ Л ,

пат гиат яат он\т 1/ / // ?

где /ишп1 - коэффициент расхода при натекании рабочего вещества в парную полость (ПП) через нагнетательное окно; - текущее значение суммарной площади

цилиндрической и торцовой части окна нагнетания в функции угла поворота ВЩ винта; рн - плотность пара в камере нагнетания.

В качестве модельной принималась скорость потока пара при адиабатном течении, которая определяется выражением

с

^нат

Л

1-

к -1

г х р

к

УРн У

где к - показатель изоэнтропы, Я - газовая постоянная, Тн - температура пара на

нагнетании; р - величина текущего давления.

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

ц.нат= 0,146 -10"2 3 Яе+ 0,24-103 +0,1.

Расчет изменения параметров состояния хладагента на этом этапе производится по тем же зависимостям, что и на первом этапе процесса сжатия [2], но с учетом натекающего из окна нагнетания пара.

Рассмотрим процесс выталкивания компримируемой среды в-с (рис. 1). После того, как давление в парной полости достигнет давления в камере нагнетания, начнется собственно выталкивание пара через нагнетательное окно, при этом давление в парной полости будет выше давления нагнетания на величину потерь.

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

Р -с2

Ар =1 • " выт (1)

¿^н Ъвыт 2 ' ^ '

где рн - плотность пара хладагента по условиям нагнетания; свыт - средняя скорость

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

Действительный объем рабочего вещества, вытесненный в камеру нагнетания из одной парной полости

Уд,=

Рве

Р И

где А - коэффициент подачи компрессора; м>п - полезный объем парной полости; рес -

удельная плотность пара по условиям всасывания. С другой стороны

^д.н ^ выт ^он.ср ^выт'

где свыт - искомая средняя скорость выталкивания смеси в процессе в-с; 1:<т ср - средняя площадь окна нагнетания в этом процессе; 1вът- время процесса.

Среднюю площадь окна нагнетания процесса в-с определяем как

Фс /

рон.ср = ф).^ кфс-фв). а /

Время процесса выталкивания

, = Фс - фе

1выт

2-п-п

1

Тогда уравнение (2.27) принимает вид:

Ч'с

V =с *__П)

д. и выт г\ V /

2 • п ■ п

где п - частота вращения ВЩ винта. Сравнивая (2) и (3) получим:

2пп

выт * п т

р

'н

\Рон( Ф^Ф!

Фе

В результате обработки индикаторных диаграмм рабочего процесса ВКМ получена зависимость „ от Яе

~ выт

6,857 -106

зон Яе+\2\,1Л0Ъ

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

Были определены следующие величины, характеризующие работу компрессора, в зависимости от геометрической степени сжатия 8 г. Индикаторный КПД компрессора

Л/

где Ц - работа компрессора при изоэнтропном сжатии; Ц - индикаторная работа

компрессора, определенная по площади индикаторной диаграммы. Относительные потери давления в процессе нагнетания

£Фн=^н/(Рн-Рес)>

где Арн - величина сопротивления нагнетательного окна в процессе нагнетания; рн и рв - величины давлений нагнетания и всасывания.

Относительной потери работы компрессора в процессе нагнетания

АЬи = АЬи / ,

гг гг I

где АЬН определялась по величине соответствующей площади (рис. 1).

Потеря работы компрессора из-за процесса натекания в случае, когда

ра<рн, определялась как разница площади индикаторной диаграммы с учетом и без учета натечек из окна нагнетания (рис. 1).

Относительная потеря работы в этом процессе определялась по формуле

^выт = ^выт / Ц ■

Расчеты были произведены для винтового компрессора имеющего следующие основные геометрические характеристики. Соотношение числа зубьев ведущего и ведомого винтов 4/6; внешний диаметр ведущего винта 215 мм; относительная длина винтов 1,4; ход ведущего вита 360 мм; полезный объем парной полости 1653,1 см ; частота вращения ведущего винта 49 с-1; профили зубьев винтов выполнены по [4].

Расчеты были проведены при работе компрессора на озонобезопасных хладагентах Я717, Я404а, Я23 и Я744.

Были выбраны следующие режимы работы компрессора [5,6].

Режим работы компрессора в составе среднетемпературной холодильной машины на хладагентах Я717 и Я404а: температура кипения - 15 0С, температура конденсации +30 0С. Аммиак и хладон Я22 широко используются в промышленных холодильных установках. При работе ВКМ на этом температурном режиме был рассмотрен смесевый хладагент Я404а, так как в настоящее время он находит применение в качестве альтернативы Я22.

Режим работы компрессора в составе нижней ветви каскадной холодильной машины на хладагентах Я23 и Я744: температура кипения - 50 0С, температура конденсации - 10 0С.

При расчете термодинамических и теплофизических параметров хладагентов использовались работы [7.. .12].

Результаты расчетов приведены на рис. 2.5.

Выполненные расчеты показали, что в режиме среднетемпературной холодильной машины энергетическая эффективность работы винтового компрессора на хладагенте Я717 выше, чем на Я404а при всех геометрических степенях сжатия, прежде всего из-за большей потери давления при нагнетании хладагента Я404а вследствие большей его плотности в процессе выталкивания.

Рис. 2. Зависимость индикаторного КПД компрессора Т|г от геометрической

степени сжатия 8

Рис. 3. Зависимость относительного сопротивления в процессе нагнетания Арн от геометрической

степени сжатия 8 у

Рис. 4. Зависимость относительной потери работы компрессора в процессе нагнетания от

геометрической степени сжатия 8 у

АЬ,

на&,04 0,035 0,03 0,025 0,02 0,015 0,01 0,005 0

\

^

N \

\

* \

\ ч ч

Р404а

Р717

Р23

х Р744

2,0

2,5

3,0

3,5

Рис. 5. Зависимость относительной потери работы компрессора в процессе натекания АЬтт от

геометрической степени сжатия £ у

В режиме работы компрессора в составе нижней ветви каскадной холодильной машины энергетическая эффективность работы винтового компрессора на хладагенте Я744 выше, чем на Я23 при всех геометрических степенях сжатия. Это

происходит из-за большего относительного сопротивления в процессе нагнетания хладагента Я23 вследствие меньшего значения индикаторной работы компрессора и примерно равной абсолютной потери давления в процессе выталкивания.

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

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

1. Носков А.Н. Расчет процесса сжатия холодильного винтового компрессора // Компрессорная техника и пневматика. -1997. - № 1-2. - С. 35-39.

2. Носков А. Н., Зимков А.А. Расчет процесса сжатия маслозаполненного холодильного винтового компрессора // Холодильная техника и кондиционирование. - №1. - 2013.

3. Носков А. Н., Зимков А.А. Расчет процесса всасывания маслозаполненного холодильного винтового компрессора // Холодильная техника и кондиционирование. -№1. - 2012.

4. Пат. 2109170 России, МКИ Б 04 С 18/16. Зубчатое зацепление винтового компрессора / Носков А.Н. 1998. - Бюл. №11.

5. Бараненко А.В., Калюнов В.С., Румянцев Ю.Д. Практикум по холодильным установкам - СПб.: Профессия, 2001. -270 с.

6. Крупененков Н.Ф. К вопросу применения эффекта Ранка-Хильша (вихревая труба) на предприятиях по производству колбасных изделий // Процессы и аппараты пищевых производств. 2013. № 1. С. 23.

7. Бухарин Н.Н. Моделирование характеристик центробежных компрессоров. Л.: Машиностроение. 1983. -214 С.

8. Рыков А.В., Кудрявцева И.В., Рыков В.А Асимметричное масштабное уравнение состояния хладона Я23 // Вестник Международной академии холода. 2012. № 4. С. 26-28.

9. Кудрявцева И.В., Рыков В.А., Рыков С.В., Селина Е.Г., Курова Л.В. Методы расчета плотности и теплоты парообразования двуокиси углерода // Процессы и аппараты пищевых производств. 2013. № 1. С. 25.

10. Кудрявцева И.В., Камоцкий В.И., Рыков С.В., Рыков В.А. Расчет линии фазового равновесия диоксида углерода // Процессы и аппараты пищевых производств. 2013. № 2. С. 31.

11. Рыков С.В., Кудрявцева И.В., Демина Л.Ю. Единое уравнение состояния Я717, учитывающее особенности критической области // Вестник Международной академии холода. 2009. № 4. С. 29-32.

12. Рыков С.В., Самолетов В.А., Рыков В.А. Линия насыщения аммиака // Вестник Международной академии холода. 2008. № 4. С. 20-21.