Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора
Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»
№ 1, 2014
УДК 621.514
Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора
Д-р техн. наук, проф. Носков А.Н. dn.noskoff@rambler.ru Зверев Д.И kayff007@mail.ru,
Тарасенков Д.С. d.tarasenkov@mail.ru
Университет ИТМО Институт холода и биотехнологий 191002, Санкт-Петербург, ул. Ломоносова, 9
Приведена методика расчета процесса нагнетания винтового компрессора, учитывающая особенности холодильного маслозаполненного компрессора. Рассматривается два этапа процесса нагнетания при давлении внутреннего сжатия ниже давления нагнетания. Определены зависимости для расчета давления в парной полости винтового компрессора на этих этапах. Рассмотрена работа компрессора в составе среднетемпературной холодильной машины на хладагентах R717 и R404а, а в составе нижней ветви каскадной холодильной машины на хладагентах R23 и R744. Приведены зависимости индикаторного КПД компрессора, относительного сопротивления в процессе нагнетания, относительной потери работы компрессора в процессе нагнетания, относительной потери работы компрессора в процессе натекания хладагента в парную полость из камеры нагнетания от геометрической степени сжатия. В режиме среднетемпературной холодильной машины энергетическая эффективность работы винтового компрессора на хладагенте R717 выше, чем на R404а, а составе нижней ветви каскадной холодильной машины энергетическая эффективность работы на хладагенте R744 выше, чем на R23 при всех геометрических степенях сжатия.
Ключевые слова: холодильный винтовой компрессор, процесс нагнетания, потери работы компрессора в процессе нагнетания
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
Носков А.Н. и др. Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора / А.Н. Носков, Д.И. Зверев, Д.С. Тарасенков // Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование», 2014. - № 1.[Электронный ресурс]: http://www.refrigeration.ihbt.ifmo.ru/
Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»
№ 1, 2014
R717 above then on R404а, 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с , до угла, на котором происходит выравнивание давлений в парной полости и в окне нагнетания − в . Затем происходит выталкивание рабочего вещества из парной полости в камеру нагнетания. Этот процесс происходит на угле поворота ВЩ винта от в до с .
Носков А.Н. и др. Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора / А.Н. Носков, Д.И. Зверев, Д.С. Тарасенков // Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование», 2014. - № 1.[Электронный ресурс]: http://www.refrigeration.ihbt.ifmo.ru/
Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»
№ 1, 2014
Рис. 1. Схема процесса сжатия и нагнетания
Процесс выталкивания сопровождается потерями давления, т.е. среднее по углу поворота давление в рассматриваемой парной полости будет больше давления нагнетания на величину рн .
Метод определения изменения параметров рабочего вещества ВКМ на первом этапе процесса сжатия приведен в работе [2].
Рассмотрим второй этап процесса сжатия в ВКМ. Масса пара, поступившего в парную полость из камеры нагнетания за рассматриваемый отрезок времени
dGнат
,нат снат Fон ( 1) н dt
где нат - коэффициент расхода при натекании рабочего вещества в парную полость
(ПП) через нагнетательное окно; Fон ( 1) - текущее значение суммарной площади
цилиндрической и торцовой части окна нагнетания в функции угла поворота ВЩ винта;
н - плотность пара в камере нагнетания. В качестве модельной принималась скорость потока пара при адиабатном течении,
которая определяется выражением
Носков А.Н. и др. Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора / А.Н. Носков, Д.И. Зверев, Д.С. Тарасенков // Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование», 2014. - № 1.[Электронный ресурс]: http://www.refrigeration.ihbt.ifmo.ru/
Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»
№ 1, 2014
k1
снат
2 k1
k R Тн
1
р рн
k. ,
где k - показатель изоэнтропы, R - газовая постоянная, Tн - температура пара на
нагнетании; p - величина текущего давления. Среднее значение коэффициента расхода при натекании пара в ПП из окна
нагнетания в зависимости от числа Рейнольдса определяется по формуле
нат 0,146 10 2 3 Re 0,24 103 0,1.
Расчет изменения параметров состояния хладагента на этом этапе производится по тем же зависимостям, что и на первом этапе процесса сжатия [2], но с учетом натекающего из окна нагнетания пара.
Рассмотрим процесс выталкивания компримируемой среды в с (рис. 1). После того, как давление в парной полости достигнет давления в камере нагнетания, начнется собственно выталкивание пара через нагнетательное окно, при этом давление в парной полости будет выше давления нагнетания на величину потерь. Поскольку величина потерь в этом процессе влияет, прежде всего, на площадь индикаторной диаграммы рабочего процесса холодильного ВКМ и через нее на энергетические коэффициенты компрессора, потеря давления в процессе выталкивания
оценивалась посредством коэффициента сопротивления н , определяемого по формуле
рн
выт
н св2ыт , 2
(1)
где н - плотность пара хладагента по условиям нагнетания; свыт - средняя скорость
выталкивания смеси из полости в рассматриваемом процессе.
Значение свыт может быть определено следующим образом.
Действительный объем рабочего вещества, вытесненный в камеру нагнетания из
одной парной полости
Vд.н wn вс ,
н
где - коэффициент подачи компрессора; wп - полезный объем парной полости;
удельная плотность пара по условиям всасывания. С другой стороны
вс -
Vд.н свыт Fон.ср tвыт ,
(2)
Носков А.Н. и др. Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора / А.Н. Носков, Д.И. Зверев, Д.С. Тарасенков // Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование», 2014. - № 1.[Электронный ресурс]: http://www.refrigeration.ihbt.ifmo.ru/
Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»
№ 1, 2014
где свыт - искомая средняя скорость выталкивания смеси в процессе в с ; Fон.ср - средняя
площадь окна нагнетания в этом процессе; tвыт - время процесса. Среднюю площадь окна нагнетания процесса в с определяем как
Fон.ср
Время процесса выталкивания
с
Fон (
)
1
d
1
в
(c
в) .
tвыт
с
2
в
п
.
1
Тогда уравнение (2.27) принимает вид:
где n1- частота вращения ВЩ винта. Сравнивая (2) и (3) получим:
cвыт
с
Fон (
)d
1
1
Vд.н свыт в 2
n 1
wп
вс н
2п 1 с
.
Fон (
)d
1
1
(3)
в
В результате обработки индикаторных диаграмм рабочего процесса ВКМ получена
зависимость выт от Re
6,857 106 он Re 121,7 103 .
По приведенному выше методу и методу представленному в работах [2,3] был произведен расчет рабочего процесса холодильного ВКМ, работающего на различных рабочих веществах, и построена расчетная индикаторная диаграмма винтового компрессора.
Были определены следующие величины, характеризующие работу компрессора, в
зависимости от геометрической степени сжатия Г . Индикаторный КПД компрессора
i Ls / Li , где Ls - работа компрессора при изоэнтропном сжатии; Li - индикаторная работа компрессора, определенная по площади индикаторной диаграммы.
Относительные потери давления в процессе нагнетания
pн pн /( pн pвс ),
Носков А.Н. и др. Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора / А.Н. Носков, Д.И. Зверев, Д.С. Тарасенков // Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование», 2014. - № 1.[Электронный ресурс]: http://www.refrigeration.ihbt.ifmo.ru/
Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»
№ 1, 2014
где pн - величина сопротивления нагнетательного окна в процессе нагнетания; pн и pв - величины давлений нагнетания и всасывания.
Относительной потери работы компрессора в процессе нагнетания
Lн Lн / Li , где Lн определялась по величине соответствующей площади (рис. 1).
Потеря работы компрессора из-за процесса натекания Lнат в случае, когда
ра рн , определялась как разница площади индикаторной диаграммы с учетом и без учета натечек из окна нагнетания (рис. 1).
Относительная потеря работы в этом процессе определялась по формуле
Lвыт Lвыт / Li . Расчеты были произведены для винтового компрессора имеющего следующие основные геометрические характеристики. Соотношение числа зубьев ведущего и ведомого винтов 4/6; внешний диаметр ведущего винта 215 мм; относительная длина винтов 1,4; ход ведущего вита 360 мм; полезный объем парной полости 1653,1 см3; частота вращения ведущего винта 49 с-1; профили зубьев винтов выполнены по [4]. Расчеты были проведены при работе компрессора на озонобезопасных хладагентах R717, R404а, R23 и R744. Были выбраны следующие режимы работы компрессора [5,6]. Режим работы компрессора в составе среднетемпературной холодильной машины на хладагентах R717 и R404а: температура кипения – 15 0С, температура конденсации +30 0С. Аммиак и хладон R22 широко используются в промышленных холодильных установках. При работе ВКМ на этом температурном режиме был рассмотрен смесевый хладагент R404а, так как в настоящее время он находит применение в качестве альтернативы R22. Режим работы компрессора в составе нижней ветви каскадной холодильной машины на хладагентах R23 и R744: температура кипения – 50 0С, температура конденсации – 10 0С. При расчете термодинамических и теплофизических параметров хладагентов использовались работы [7…12]. Результаты расчетов приведены на рис. 2…5. Выполненные расчеты показали, что в режиме среднетемпературной холодильной машины энергетическая эффективность работы винтового компрессора на хладагенте R717 выше, чем на R404а при всех геометрических степенях сжатия, прежде всего из-за большей потери давления при нагнетании хладагента R404а вследствие большей его плотности в процессе выталкивания.
Носков А.Н. и др. Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора / А.Н. Носков, Д.И. Зверев, Д.С. Тарасенков // Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование», 2014. - № 1.[Электронный ресурс]: http://www.refrigeration.ihbt.ifmo.ru/
Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»
№ 1, 2014
i 0,9
0,8
0,7
0,6
0,5
0,4 2,0 2,5 3,0 3,5 4,0 4,5 5,0 5,5
Г
R717 R404а R23 R744
Рис. 2. Зависимость индикаторного КПД компрессора i от геометрической степени сжатия Г
pн0,4
0,3 R404а R717
0,2 R23 R744
0,1
0 2,0 2,5 3,0 3,5 4,0 4,5 5,0 5,5 Г
Рис. 3. Зависимость относительного сопротивления в процессе нагнетания pн от геометрической
степени сжатия Г
Носков А.Н. и др. Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора / А.Н. Носков, Д.И. Зверев, Д.С. Тарасенков // Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование», 2014. - № 1.[Электронный ресурс]: http://www.refrigeration.ihbt.ifmo.ru/
Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»
№ 1, 2014
Рис. 4. Зависимость относительной потери работы компрессора в процессе нагнетания Lн от
геометрической степени сжатия Г
Lна0т,04
0,035
0,03 0,025
0,02 0,015
0,01
R404а R717 R23 R744
0,005
0 2,0 2,5 3,0 3,5 Г
Рис. 5. Зависимость относительной потери работы компрессора в процессе натекания Lнат от
геометрической степени сжатия Г
В режиме работы компрессора в составе нижней ветви каскадной холодильной машины энергетическая эффективность работы винтового компрессора на хладагенте R744 выше, чем на R23 при всех геометрических степенях сжатия. Это
Носков А.Н. и др. Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора / А.Н. Носков, Д.И. Зверев, Д.С. Тарасенков // Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование», 2014. - № 1.[Электронный ресурс]: http://www.refrigeration.ihbt.ifmo.ru/
Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»
№ 1, 2014
происходит из-за большего относительного сопротивления в процессе нагнетания хладагента R23 вследствие меньшего значения индикаторной работы компрессора и примерно равной абсолютной потери давления в процессе выталкивания.
Полученные данные позволяют более обосновано подойти к выбору озонобезопасных хладагентов при работе винтового компрессора на различных режимах.
Список литературы
1. Носков А.Н. Расчет процесса сжатия холодильного винтового компрессора // Компрессорная техника и пневматика. −1997. − № 1−2. − С. 35-39. 2. Носков А. Н., Зимков А.А. Расчет процесса сжатия маслозаполненного холодильного винтового компрессора // Холодильная техника и кондиционирование. − №1. − 2013. 3. Носков А. Н., Зимков А.А. Расчет процесса всасывания маслозаполненного холодильного винтового компрессора // Холодильная техника и кондиционирование. − №1. − 2012. 4. Пат. 2109170 России, МКИ F 04 С 18/16. Зубчатое зацепление винтового компрессора / Носков А.Н. 1998. − Бюл. №11. 5. Бараненко А.В., Калюнов В.С., Румянцев Ю.Д. Практикум по холодильным установкам − СПб.: Профессия, 2001. -270 с. 6. Крупененков Н.Ф. К вопросу применения эффекта Ранка-Хильша (вихревая труба) на предприятиях по производству колбасных изделий // Процессы и аппараты пищевых производств. 2013. № 1. С. 23. 7. Бухарин Н.Н. Моделирование характеристик центробежных компрессоров. Л.: Машиностроение. 1983. −214 С. 8. Рыков А.В., Кудрявцева И.В., Рыков В.А Асимметричное масштабное уравнение состояния хладона R23 // Вестник Международной академии холода. 2012. № 4. С. 26−28. 9. Кудрявцева И.В., Рыков В.А., Рыков С.В., Селина Е.Г., Курова Л.В. Методы расчета плотности и теплоты парообразования двуокиси углерода // Процессы и аппараты пищевых производств. 2013. № 1. С. 25. 10. Кудрявцева И.В., Камоцкий В.И., Рыков С.В., Рыков В.А. Расчет линии фазового равновесия диоксида углерода // Процессы и аппараты пищевых производств. 2013. № 2. С. 31. 11. Рыков С.В., Кудрявцева И.В., Демина Л.Ю. Единое уравнение состояния R717, учитывающее особенности критической области // Вестник Международной академии холода. 2009. № 4. С. 29−32. 12. Рыков С.В., Самолетов В.А., Рыков В.А. Линия насыщения аммиака // Вестник Международной академии холода. 2008. № 4. С. 20−21.
Носков А.Н. и др. Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора / А.Н. Носков, Д.И. Зверев, Д.С. Тарасенков // Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование», 2014. - № 1.[Электронный ресурс]: http://www.refrigeration.ihbt.ifmo.ru/
№ 1, 2014
УДК 621.514
Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора
Д-р техн. наук, проф. Носков А.Н. dn.noskoff@rambler.ru Зверев Д.И kayff007@mail.ru,
Тарасенков Д.С. d.tarasenkov@mail.ru
Университет ИТМО Институт холода и биотехнологий 191002, Санкт-Петербург, ул. Ломоносова, 9
Приведена методика расчета процесса нагнетания винтового компрессора, учитывающая особенности холодильного маслозаполненного компрессора. Рассматривается два этапа процесса нагнетания при давлении внутреннего сжатия ниже давления нагнетания. Определены зависимости для расчета давления в парной полости винтового компрессора на этих этапах. Рассмотрена работа компрессора в составе среднетемпературной холодильной машины на хладагентах R717 и R404а, а в составе нижней ветви каскадной холодильной машины на хладагентах R23 и R744. Приведены зависимости индикаторного КПД компрессора, относительного сопротивления в процессе нагнетания, относительной потери работы компрессора в процессе нагнетания, относительной потери работы компрессора в процессе натекания хладагента в парную полость из камеры нагнетания от геометрической степени сжатия. В режиме среднетемпературной холодильной машины энергетическая эффективность работы винтового компрессора на хладагенте R717 выше, чем на R404а, а составе нижней ветви каскадной холодильной машины энергетическая эффективность работы на хладагенте R744 выше, чем на R23 при всех геометрических степенях сжатия.
Ключевые слова: холодильный винтовой компрессор, процесс нагнетания, потери работы компрессора в процессе нагнетания
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
Носков А.Н. и др. Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора / А.Н. Носков, Д.И. Зверев, Д.С. Тарасенков // Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование», 2014. - № 1.[Электронный ресурс]: http://www.refrigeration.ihbt.ifmo.ru/
Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»
№ 1, 2014
R717 above then on R404а, 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с , до угла, на котором происходит выравнивание давлений в парной полости и в окне нагнетания − в . Затем происходит выталкивание рабочего вещества из парной полости в камеру нагнетания. Этот процесс происходит на угле поворота ВЩ винта от в до с .
Носков А.Н. и др. Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора / А.Н. Носков, Д.И. Зверев, Д.С. Тарасенков // Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование», 2014. - № 1.[Электронный ресурс]: http://www.refrigeration.ihbt.ifmo.ru/
Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»
№ 1, 2014
Рис. 1. Схема процесса сжатия и нагнетания
Процесс выталкивания сопровождается потерями давления, т.е. среднее по углу поворота давление в рассматриваемой парной полости будет больше давления нагнетания на величину рн .
Метод определения изменения параметров рабочего вещества ВКМ на первом этапе процесса сжатия приведен в работе [2].
Рассмотрим второй этап процесса сжатия в ВКМ. Масса пара, поступившего в парную полость из камеры нагнетания за рассматриваемый отрезок времени
dGнат
,нат снат Fон ( 1) н dt
где нат - коэффициент расхода при натекании рабочего вещества в парную полость
(ПП) через нагнетательное окно; Fон ( 1) - текущее значение суммарной площади
цилиндрической и торцовой части окна нагнетания в функции угла поворота ВЩ винта;
н - плотность пара в камере нагнетания. В качестве модельной принималась скорость потока пара при адиабатном течении,
которая определяется выражением
Носков А.Н. и др. Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора / А.Н. Носков, Д.И. Зверев, Д.С. Тарасенков // Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование», 2014. - № 1.[Электронный ресурс]: http://www.refrigeration.ihbt.ifmo.ru/
Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»
№ 1, 2014
k1
снат
2 k1
k R Тн
1
р рн
k. ,
где k - показатель изоэнтропы, R - газовая постоянная, Tн - температура пара на
нагнетании; p - величина текущего давления. Среднее значение коэффициента расхода при натекании пара в ПП из окна
нагнетания в зависимости от числа Рейнольдса определяется по формуле
нат 0,146 10 2 3 Re 0,24 103 0,1.
Расчет изменения параметров состояния хладагента на этом этапе производится по тем же зависимостям, что и на первом этапе процесса сжатия [2], но с учетом натекающего из окна нагнетания пара.
Рассмотрим процесс выталкивания компримируемой среды в с (рис. 1). После того, как давление в парной полости достигнет давления в камере нагнетания, начнется собственно выталкивание пара через нагнетательное окно, при этом давление в парной полости будет выше давления нагнетания на величину потерь. Поскольку величина потерь в этом процессе влияет, прежде всего, на площадь индикаторной диаграммы рабочего процесса холодильного ВКМ и через нее на энергетические коэффициенты компрессора, потеря давления в процессе выталкивания
оценивалась посредством коэффициента сопротивления н , определяемого по формуле
рн
выт
н св2ыт , 2
(1)
где н - плотность пара хладагента по условиям нагнетания; свыт - средняя скорость
выталкивания смеси из полости в рассматриваемом процессе.
Значение свыт может быть определено следующим образом.
Действительный объем рабочего вещества, вытесненный в камеру нагнетания из
одной парной полости
Vд.н wn вс ,
н
где - коэффициент подачи компрессора; wп - полезный объем парной полости;
удельная плотность пара по условиям всасывания. С другой стороны
вс -
Vд.н свыт Fон.ср tвыт ,
(2)
Носков А.Н. и др. Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора / А.Н. Носков, Д.И. Зверев, Д.С. Тарасенков // Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование», 2014. - № 1.[Электронный ресурс]: http://www.refrigeration.ihbt.ifmo.ru/
Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»
№ 1, 2014
где свыт - искомая средняя скорость выталкивания смеси в процессе в с ; Fон.ср - средняя
площадь окна нагнетания в этом процессе; tвыт - время процесса. Среднюю площадь окна нагнетания процесса в с определяем как
Fон.ср
Время процесса выталкивания
с
Fон (
)
1
d
1
в
(c
в) .
tвыт
с
2
в
п
.
1
Тогда уравнение (2.27) принимает вид:
где n1- частота вращения ВЩ винта. Сравнивая (2) и (3) получим:
cвыт
с
Fон (
)d
1
1
Vд.н свыт в 2
n 1
wп
вс н
2п 1 с
.
Fон (
)d
1
1
(3)
в
В результате обработки индикаторных диаграмм рабочего процесса ВКМ получена
зависимость выт от Re
6,857 106 он Re 121,7 103 .
По приведенному выше методу и методу представленному в работах [2,3] был произведен расчет рабочего процесса холодильного ВКМ, работающего на различных рабочих веществах, и построена расчетная индикаторная диаграмма винтового компрессора.
Были определены следующие величины, характеризующие работу компрессора, в
зависимости от геометрической степени сжатия Г . Индикаторный КПД компрессора
i Ls / Li , где Ls - работа компрессора при изоэнтропном сжатии; Li - индикаторная работа компрессора, определенная по площади индикаторной диаграммы.
Относительные потери давления в процессе нагнетания
pн pн /( pн pвс ),
Носков А.Н. и др. Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора / А.Н. Носков, Д.И. Зверев, Д.С. Тарасенков // Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование», 2014. - № 1.[Электронный ресурс]: http://www.refrigeration.ihbt.ifmo.ru/
Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»
№ 1, 2014
где pн - величина сопротивления нагнетательного окна в процессе нагнетания; pн и pв - величины давлений нагнетания и всасывания.
Относительной потери работы компрессора в процессе нагнетания
Lн Lн / Li , где Lн определялась по величине соответствующей площади (рис. 1).
Потеря работы компрессора из-за процесса натекания Lнат в случае, когда
ра рн , определялась как разница площади индикаторной диаграммы с учетом и без учета натечек из окна нагнетания (рис. 1).
Относительная потеря работы в этом процессе определялась по формуле
Lвыт Lвыт / Li . Расчеты были произведены для винтового компрессора имеющего следующие основные геометрические характеристики. Соотношение числа зубьев ведущего и ведомого винтов 4/6; внешний диаметр ведущего винта 215 мм; относительная длина винтов 1,4; ход ведущего вита 360 мм; полезный объем парной полости 1653,1 см3; частота вращения ведущего винта 49 с-1; профили зубьев винтов выполнены по [4]. Расчеты были проведены при работе компрессора на озонобезопасных хладагентах R717, R404а, R23 и R744. Были выбраны следующие режимы работы компрессора [5,6]. Режим работы компрессора в составе среднетемпературной холодильной машины на хладагентах R717 и R404а: температура кипения – 15 0С, температура конденсации +30 0С. Аммиак и хладон R22 широко используются в промышленных холодильных установках. При работе ВКМ на этом температурном режиме был рассмотрен смесевый хладагент R404а, так как в настоящее время он находит применение в качестве альтернативы R22. Режим работы компрессора в составе нижней ветви каскадной холодильной машины на хладагентах R23 и R744: температура кипения – 50 0С, температура конденсации – 10 0С. При расчете термодинамических и теплофизических параметров хладагентов использовались работы [7…12]. Результаты расчетов приведены на рис. 2…5. Выполненные расчеты показали, что в режиме среднетемпературной холодильной машины энергетическая эффективность работы винтового компрессора на хладагенте R717 выше, чем на R404а при всех геометрических степенях сжатия, прежде всего из-за большей потери давления при нагнетании хладагента R404а вследствие большей его плотности в процессе выталкивания.
Носков А.Н. и др. Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора / А.Н. Носков, Д.И. Зверев, Д.С. Тарасенков // Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование», 2014. - № 1.[Электронный ресурс]: http://www.refrigeration.ihbt.ifmo.ru/
Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»
№ 1, 2014
i 0,9
0,8
0,7
0,6
0,5
0,4 2,0 2,5 3,0 3,5 4,0 4,5 5,0 5,5
Г
R717 R404а R23 R744
Рис. 2. Зависимость индикаторного КПД компрессора i от геометрической степени сжатия Г
pн0,4
0,3 R404а R717
0,2 R23 R744
0,1
0 2,0 2,5 3,0 3,5 4,0 4,5 5,0 5,5 Г
Рис. 3. Зависимость относительного сопротивления в процессе нагнетания pн от геометрической
степени сжатия Г
Носков А.Н. и др. Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора / А.Н. Носков, Д.И. Зверев, Д.С. Тарасенков // Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование», 2014. - № 1.[Электронный ресурс]: http://www.refrigeration.ihbt.ifmo.ru/
Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»
№ 1, 2014
Рис. 4. Зависимость относительной потери работы компрессора в процессе нагнетания Lн от
геометрической степени сжатия Г
Lна0т,04
0,035
0,03 0,025
0,02 0,015
0,01
R404а R717 R23 R744
0,005
0 2,0 2,5 3,0 3,5 Г
Рис. 5. Зависимость относительной потери работы компрессора в процессе натекания Lнат от
геометрической степени сжатия Г
В режиме работы компрессора в составе нижней ветви каскадной холодильной машины энергетическая эффективность работы винтового компрессора на хладагенте R744 выше, чем на R23 при всех геометрических степенях сжатия. Это
Носков А.Н. и др. Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора / А.Н. Носков, Д.И. Зверев, Д.С. Тарасенков // Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование», 2014. - № 1.[Электронный ресурс]: http://www.refrigeration.ihbt.ifmo.ru/
Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»
№ 1, 2014
происходит из-за большего относительного сопротивления в процессе нагнетания хладагента R23 вследствие меньшего значения индикаторной работы компрессора и примерно равной абсолютной потери давления в процессе выталкивания.
Полученные данные позволяют более обосновано подойти к выбору озонобезопасных хладагентов при работе винтового компрессора на различных режимах.
Список литературы
1. Носков А.Н. Расчет процесса сжатия холодильного винтового компрессора // Компрессорная техника и пневматика. −1997. − № 1−2. − С. 35-39. 2. Носков А. Н., Зимков А.А. Расчет процесса сжатия маслозаполненного холодильного винтового компрессора // Холодильная техника и кондиционирование. − №1. − 2013. 3. Носков А. Н., Зимков А.А. Расчет процесса всасывания маслозаполненного холодильного винтового компрессора // Холодильная техника и кондиционирование. − №1. − 2012. 4. Пат. 2109170 России, МКИ F 04 С 18/16. Зубчатое зацепление винтового компрессора / Носков А.Н. 1998. − Бюл. №11. 5. Бараненко А.В., Калюнов В.С., Румянцев Ю.Д. Практикум по холодильным установкам − СПб.: Профессия, 2001. -270 с. 6. Крупененков Н.Ф. К вопросу применения эффекта Ранка-Хильша (вихревая труба) на предприятиях по производству колбасных изделий // Процессы и аппараты пищевых производств. 2013. № 1. С. 23. 7. Бухарин Н.Н. Моделирование характеристик центробежных компрессоров. Л.: Машиностроение. 1983. −214 С. 8. Рыков А.В., Кудрявцева И.В., Рыков В.А Асимметричное масштабное уравнение состояния хладона R23 // Вестник Международной академии холода. 2012. № 4. С. 26−28. 9. Кудрявцева И.В., Рыков В.А., Рыков С.В., Селина Е.Г., Курова Л.В. Методы расчета плотности и теплоты парообразования двуокиси углерода // Процессы и аппараты пищевых производств. 2013. № 1. С. 25. 10. Кудрявцева И.В., Камоцкий В.И., Рыков С.В., Рыков В.А. Расчет линии фазового равновесия диоксида углерода // Процессы и аппараты пищевых производств. 2013. № 2. С. 31. 11. Рыков С.В., Кудрявцева И.В., Демина Л.Ю. Единое уравнение состояния R717, учитывающее особенности критической области // Вестник Международной академии холода. 2009. № 4. С. 29−32. 12. Рыков С.В., Самолетов В.А., Рыков В.А. Линия насыщения аммиака // Вестник Международной академии холода. 2008. № 4. С. 20−21.
Носков А.Н. и др. Расчет процесса нагнетания маслозаполненного холодильного винтового компрессора / А.Н. Носков, Д.И. Зверев, Д.С. Тарасенков // Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование», 2014. - № 1.[Электронный ресурс]: http://www.refrigeration.ihbt.ifmo.ru/