Реферат Расчеты схем теплонасосных и холодильных установок
Работа добавлена на сайт bukvasha.net: 2015-10-28Поможем написать учебную работу
Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
от 25%
договор
Федеральное агентство по образованию
Государственное образовательное учреждение высшего профессионального образования
«Южно-Уральский государственный университет»
Факультет Энергетический
Кафедра «Промышленная теплоэнергетика»
Пояснительная записка к семестровому заданию
по дисциплине «Теплоносители и холодильные установки»
Тема: “ Расчеты схем теплонасосных и холодильных установок”
Э-450.140104.012.00.ПЗ
Нормоконтролер: Руководитель:
Жиргалова Т.Б. _______ Заболотный В.В.
__________
“____” __________2010 г. “____” __________2010 г.
Выполнил:
Студент группы Э-450
________ Першукова В.И.
“____” __________2010 г.
Челябинск,
Аннотация
Заболотный В.В.– Челябинск: ЮУрГУ, Э, 2010, 19 с. Библиография литературы – 10 наим.
В данной семестровой работе проведен расчет компрессорной парожидкостной одноступенчатой холодильной установки с рекуперативным теплообменником. По полученным расчетам было выбрано основное оборудование. Проведен эксергетический баланс холодильной установки.
Содержание.
Введение
1.Расчет холодильного цикла
1.1 Принципиальная схема установки
1.2. Расчёт основных параметров цикла
2. Характерные величины цикла
3. Эксергетический анализ
4.Подбор оборудования
Заключение.
Литература
Приложения
Введение.
В настоящей работе приведён энергетический и эксергетический анализ холодильной установки, подобрано основное оборудование.
Приведённые материалы являются исходной базой для разработки проекта холодильной установки.
Разработанные мероприятия направлены на:
ü Повышение качества холодоснабжения;
ü Снижение капитальных затрат на создание системы холодоснабжения;
ü Снижение затрат на содержание системы холодоснабжения.
1.Расчет холодильного цикла
1.1. Принципиальная схема установки.
На рис.1. приведена схема для ХЛУ
Рис.1. Принципиальная схема парожидкостной холодильной установки: КМ - компрессор; К - конденсатор; ПО - переохладитель; РД - регулирующий вентиль; И - испаритель.
Принцип действия установки состоит в следующем. Компрессор отсасывает пары хладагента из испарителя при давлении Ро и сжимает их до давления Рк (процесс 1-2 на рис.2). Сжатый пар поступает в конденсатор, где он под воздействием ВИТ (вода) переходит в состояние насыщения, а затем конденсируется при температуре Тк (процесс 2-3 на рис.2). При этом ВИТ нагревается от температуры Т’в до температуры Т’’вк. Жидкий агент направляется в переохладитель ПО (рис.1.), где происходит понижение его температуры на величину DТ=Т3–Т4 (процесс 3-4 на рис.2). При дросселировании агента в РД происходит понижение его давления и температуры с переходом в состояние влажного пара. Жидкая фаза, поступая в испаритель, кипит при температуре Т0 за счет тепла, подводимого от НИТ
(процесс 5 – 1 на рис.2). Процесс на T,s-диаграмме парожидкостной холодильной установки приведен на рис.2.
Рис.2. Процесс в T,s-диаграмме парожидкостной холодильной установки.
1.2. Расчёт основных параметров цикла.
Расчетная температура испарения (рис.3), давление и энтальпия, энтропия в т.1
t0 = tн2 - Δtи = -10 - 5 = -15 оС
Рис.3. Расчетная температура испарения.
По T,s-диаграмме и таблицам(прил.1) определяем давление и энтальпию. энтропию в т.1: p0 = 0,22 МПа, h1 = 1675 кДж/кг, s1 = 9051 Дж/(кгК).
Расчетная температура (рис.4) и давление конденсации, энтальпия и энтропия в т.3:
tк = tв1 + Δtк = 24+7=31 оС
Рис.4. Расчетная температура конденсации.
По T,s-диаграмме и таблицам (прил.1) определяем давление и энтальпию. энтропию в т.3: pК = 1,203 МПа, h3 = 566 кДж/кг, s3 = 4683 Дж/(кгК).
Температура, энтальпия и энтропия в т.2:
Определяем энтальпию паров хладона после идеального компрессора т.2’
h2’ = 1922 кДж/кг.
Принимаю внутренний адиабатный КПД компрессора: ηi =0,8
Определяем энтальпию рабочего агента на выходе из компрессора т.2:
h2 = h1 +ав/ηi=h1+(h’2 - h1)/ηi = 1675+(1922-1675)/0,8 = 1983,75 кДж/кг.
По T,s-диаграмме (прил.1) определяем температуру энтропию в т.2:
t2 = 128 оС,
s2 = 9220 Дж/(кг К).
Температура, энтальпия и энтропия в т.4: Δt’ = 9 оС,
t4 = tкон + Δt’ = 19 оС
По T,s-диаграмме (прил.1) : h4 = 503 кДж/кг
s4 = 4410 Дж/(кг К).
Энтропия в т.5:
По T,s-диаграмме (прил.1) определяем энтропию в т.5: s5 = 4706,15 Дж/(кг К).
Проверочный расчет энтропии в т.5:
s5c = s’·(1–х) + s”·х=3927·(1-0,151)+ 9051·0,151 = 4700,724 Дж/(кг×К).
где s’=3927 Дж/(кг×К), s”=9051 Дж/(кг×К) по таблице прил.1,
х – степень сухости и равна x=0,151 (прил.1).
Погрешность расчета энтропии в т.5:
δ = 1 - s5c / s5 = 0,12%
Точность определения энтропии в т.5 достаточная для инженерных расчетов (0,12% < 1%)
Значения удельных эксергий хладона в характерных точках процесса определим по е, h -диаграмме и по формуле:
e = h - hо.с. - Tо.с.(s - sо.с) = h - Tо.с. s - (hо.с. - Tо.с. sо.с)
Найдем по T,s-диаграмме (прил.1) энтропию. энтальпию хладона R717 при параметрах окружающей среды ро.с. = 0,1 МПа, Tо.с. = 293 K:
hо.с.= 1750 кДж/кг, sо.с = 9765 Дж/(кг×К)
Значения температуры, давления, энтальпии, энтропии, эксергии точек 1, 2, 3, 4, 5 сведём в таблицу 1.
Таблица 1
| t,oC | T, K | р, Мпа | h, кДж/кг | v, м3/кг | s, Дж/(кг К) | е, кДж/кг |
1 | -15 | 258,16 | 0,22 | 1675 | 0,5 | 9051 | 134,202 |
2 | 128 | 401,16 | 1,203 | 1983,75 | | 9220 | 393,435 |
3 | 31 | 304,16 | 1,203 | 566 | | 4683 | 305,026 |
4 | 19 | 292,16 | 1,203 | 503 | | 4410 | 322,015 |
5 | -15 | 258,16 | 0,22 | 503 | | 4706,15 | 235,243 |
Строим процесс на T,s – и e,h - диаграмме (прил. 1)
2. Характерные величины цикла
Удельная нагрузка испарителя:
q0 = h1 - h5 = 1675 – 503 = 1172 кДж/кг
Массовый расход хладагента
G=Q0 / q0= (20,93*4,187)/(1172*3,6) = 0,02 кг/с ( 72 кг/ч)
Удельная внутренняя работа компрессора:
ав = h2 - h1= 1983,75 – 1675 = 308,75 кДж/кг
Удельная тепловая нагрузка конденсатора:
qК = h2 - h3 = 1983,75 – 566 = 1417,75 кДж/кг
Удельная тепловая нагрузка переохладителя:
qпо = h3 - h4 =566 - 503 = 63 кДж/кг
Полная тепловая нагрузка переохладителя:
Qпо = qпо G = 63*0,02 = 1,26 кВт
Проверка решения по I закону термодинамики:
q = aв + q0 = 308,75 + 1172 = 1480,75 кДж/кг
q = qк + qпо = 1417,75 +63 = 1480,75 кДж/кг
Объемная производительность компрессора:
V0 = G ν1 = 72 *0,5 = 36 м3/ч (0,01 м3/с)
Тепловая нагрузка конденсатора:
QК = G qК = 0,02*1417,75 = 28,36 кВт
Удельная работа, затраченная на компрессор с учетом электромеханического КПД, при ηэм = 0,9
аК = ав / ηэм = 308,75/0,9 =343,1 кДж/кг
Электрическая мощность компрессора:
Nэ = аК G = 343,1*0,02 = 6,862 кВт
Принимаем электрическую мощность компрессора в 7 кВт.
Холодильный коэффициент:
e
= q0 / аК = 1172/343,1 = 3,42
Удельный расход электроэнергии на единицу выработанного холода:
ex = аК / q0 = 343,1 / 1172 = 0,29
3. Эксергетический анализ цикла.
Составим эксергетический баланс для 1 кг/с расхода рабочего тела.
1. Удельное количество эксергии, подведенной к установке, по измерениям на зажимах электродвигателя компрессора:
eвх = Nэ / G = 6,862 / 0,02 =343,1 кДж/кг
2. Электромеханические потери эксергии:
dэм = eвх - eвх ηэм= eвх(1 - ηэм) = 343,1(1-0,9) = 34,31 кДж/кг
3. Удельная эксергия, подводимая к компрессору:
eв = eвх - dэм = 343,1 – 34,31 = 308,79 кДж/кг
4. Внутренние потери эксергии в компрессоре:
dкм = eв – (e2 – e1) = 308,79 – (393,435 – 134,202) = 49,56 кДж/кг
5. Удельная работа компрессора:
eк = e2 – e1 = 393,435 – 134,202 = 259,233 кДж/кг
6. Эксергетический коэффициент полезного действия компрессора:
ηк = ((eв - dкм) / eв ) 100 = (308,79 – 49,56)/ 308,79)100 = 83,95 %
7. Потери эксергии в конденсаторе:
Δ e2-3 = e2 – e3 = 393,435 – 305,026 = 88,409 кДж/кг
8. Средняя температура охлаждающей воды:
Tв.ср = (Tв1 – Tв2)/2 + 273 = (24-18) / 2 +273 = 276 К
9. Коэффициент работоспособности охлаждающей воды:
(τq)в = 1 - Tо.с / Tв.ср = 1 – 293/ 276 = -0,06
10. Эксергия, полученная водой:
Δ eв2-в3 = qк (τq)н = 1417,75*0,06 = 85,069 кДж/кг
11. Потери эксергии вследствие необратимого теплообмена:
dк.т = Δ e2-3 - Δ eв2-в3 = 88,409 – 85,065 = 3,34 кДж/кг
12. КПД конденсатора:
ηк = (Δ e в2-в3 / Δ e2-3) 100% = (85,065 / 88,409) 100 = 96,22 %
Так как эксергия охлаждающей воды после конденсаторов компрессионных установок обычно не используется, то суммарные потери эксергии в конденсаторе составят:
dк = dк.т + Δ e в2-в3 = 3,34 + 85,069 = 88,409 кДж/кг
ηк = 0%
13. Эксергия, полученная в переохладителе
W = Qох / Cp Δt - расход воды
Δ e3-4 = e4 – e3 = 322,015– 305,026 = 17 кДж/кг
Δe = W Δ e3-4 - Δ e3-4 = 66,308 – 17 = 49,308 кДж/кг
14. КПД переохладителя:
ηпо = (Δ e3/ Δ e4) 100% = (305,026 / 322,015)100 = 94,72 %
15. Потери эксергии в регулирующем клапане (дросселе):
dрв = e4 – e5= 322,015 – 235,243 = 86,77 кДж/кг
16. КПД дросселя:
ηдр =( e5 / e4) 100% = (235,243 / 322,015)100 = 73 %
17. Эксергия, отданная хладоном в испарителе:
Δ e1-5 = e1 – e5 = 134,202 –235,243 = 101,041 кДж/кг
18. Средняя температура воздуха:
Tн.ср = (Tн1 – Tн2)/2 + 273 = (-5 – 10)/2 +273 = 265,5 К
21. Коэффициент работоспособности воздуха:
(τq)н = 1 - Tо.с / Tн.ср = 1 – 293 / 265,5 = - 0,1
22. Эксергия, полученная воздухом:
е0 = Δ eн2-1 = eн2 - eн1=q0 (τq)срн = 1172 *0,1 = 117,2 кДж/кг
23. Потери эксергии вследствие необратимого теплообмена в испарителе:
dи = Δ e5-1 - Δ eн2-1 = 101,041 – 117,2 = 16,159 кДж/кг
24. КПД испарителя:
ηи = (Δ eн2-1 / Δ e5-1) 100% =(101,041 /117,2)100 = 86,21 %
Эксергетический баланс установки приведен в табл. 2.
Таблица 2
Коэффициент полезного действия теплообменно-дроссельной части установки:
ηт.д = (e0 /( Δe2-1)) 100% = (117,2 /259,233)100% = 45,21%
Полный КПД установки:
η' = e0 /( eвх) 100% = (117,2 /343,1)100% = 34,2%
4.Подбор оборудования.
Холодопроизводительность установки:
Q0 = 20,93 кВт
Подбор оборудования ведем по программе Select 7/39780
Основные преимущества хладагента аммиак (R717) обусловлены тем, что он:
- обладает термодинамическими и теплофизическими характеристиками, позволяющими получать высокий КПД в холодильных установках;
- химически нейтрален по отношению к большинству конструкционных материалов холодильных установок, за исключением меди и сплавов на ее основе;
- не растворяется в смазочных маслах, применяемых в конструкциях холодильных установок, не чувствителен к влаге и легко обнаруживается в случае утечки;
- не способствует созданию парникового эффекта;
- имеет невысокую стоимость (не более 2200 рублей за тонну) и легко доступен на рынке.
Вместе с тем у аммиака есть ряд серьезных недостатков. В частности, это вещество:
- обладает высокой токсичностью (считается, что предельно допустимая концентрация аммиака в рабочих помещениях должна быть не выше 20 мг/м3, однако даже при более слабой концентрации характерный запах аммиака в случае его появления вызывает сильную панику; при более высоких концентрациях появляются серьезные затруднения дыхания вплоть до удушья; смертельная концентрация аммиака - 30 г/м3);
- является взрывоопасным (при концентрации в воздухе 200-300 г/м3 возникает угроза самопроизвольного взрыва; температура самовоспламенения равна 650 °С);
- создает опасность ожогов при растворении в воде, поскольку этот процесс
сопровождается выделением значительного количества тепла;
- имеет высокую температуру нагнетания при сжатии в холодильных компрессорах.
Поэтому его могут заменять фреонами. В качестве замены для подбора оборудования применим R600a. Краткие физические свойства приведены в таблице 3.
Таблица 3
Расчет компрессоров заключается в определение необходимого действительного объема, описываемого поршнями в единицу времени
, м3/с ,
где V0 — теоретический объем, м3/с ,определяется из выражения(6.6);
l — коэффициент подачи компрессора, он рассчитывается:
l=lс l¶р lt lпл .
выбираем несколько однотипных компрессоров, включенных параллельно в схему установки и обеспечивающих суммарную подачу. Общее число компрессоров с учетом одного резервного определится:
К = Краб + 1 .
Объем, описываемый поршнями в единицу времени одним компрессором
, м3 /с.
Расход воды (воздуха) на конденсатор:
, м3/с,
где Cp — теплоёмкость воды (воздуха), кДж/(кг×К); — нагрев воды (воздуха) в конденсаторе.
Перепад давления Δp = 0,5 бар
Переохлаждение ΔT = 4 оС
Результаты расчета и схема установки с приборами измерения представлены в приложениях.
Заключение.
В ходе выполнения семестрового задания мной была рассчитана холодильная установка, определены основные энергетические параметры.
Рекомендации:
1. Заменить хладагент R717 на R600a
2. Разбить установку на несколько для увеличения надежности.
Литература
1. Везиришвили О.Ш., Меладзе Н.В. Энергосберегающие теплонасосные системы тепло – и хладоснабжения. — М.:Изд-во МЭИ, 1994. — 160с., ил.
2. Горбенко В.И., Юртаев М.А. Расчёт парожидкостных компрессионных одноступенчатых холодильных и теплонасосных установок. - Ч: Изд. ЮУрГУ, 2002. — 38 с., ил.
3. Горбенко В.И. Конспект лекций. - Ч: Изд. ЮУрГУ, 2002. — 38 с., ил.
4. Курылев Е.С., Герасимов Н.А. Холодильные установки. – Л.: Машиностроение, 1980. – 622 с., ил.
5. Мартынов А.В. Установки для трансформации тепла и охлаждения. – М.: Энергоиздат, 1989. – 320 с., ил.
6. Соколов Е.Я., Бродянский В.М. Энергетические основы трансформации тепла и процессов охлаждения. – М.: Энергоиздат, 1981. – 320 с., ил.
7. Теплофизические основы получения искусственного холода: Справочник серии «Холодильная техника». — М.: Пищевая промышленность, 1980. — 232с.
8. Холодильные компрессоры: Справочник серии «Холодильная техника». — М.: Легкая и пищевая промышленность, 1981. — 280с.
9. Холодильные машины: Справочник серии «Холодильная техника». — М.: Легкая и пищевая промышленность, 1982. — 223с.
10. Щербин Н.К., Гринберг Я.И. Холодильные станции и установки. — М.: Химия, 1979. — 376с.
dк = dк.т + Δ e в2-в3 = 3,34 + 85,069 = 88,409 кДж/кг
ηк = 0%
13. Эксергия, полученная в переохладителе
W = Qох / Cp Δt - расход воды
Δ e3-4 = e4 – e3 = 322,015– 305,026 = 17 кДж/кг
Δe = W Δ e3-4 - Δ e3-4 = 66,308 – 17 = 49,308 кДж/кг
14. КПД переохладителя:
ηпо = (Δ e3/ Δ e4) 100% = (305,026 / 322,015)100 = 94,72 %
15. Потери эксергии в регулирующем клапане (дросселе):
dрв = e4 – e5= 322,015 – 235,243 = 86,77 кДж/кг
16. КПД дросселя:
ηдр =( e5 / e4) 100% = (235,243 / 322,015)100 = 73 %
17. Эксергия, отданная хладоном в испарителе:
Δ e1-5 = e1 – e5 = 134,202 –235,243 = 101,041 кДж/кг
18. Средняя температура воздуха:
Tн.ср = (Tн1 – Tн2)/2 + 273 = (-5 – 10)/2 +273 = 265,5 К
21. Коэффициент работоспособности воздуха:
(τq)н = 1 - Tо.с / Tн.ср = 1 – 293 / 265,5 = - 0,1
22. Эксергия, полученная воздухом:
е0 = Δ eн2-1 = eн2 - eн1=q0 (τq)срн = 1172 *0,1 = 117,2 кДж/кг
23. Потери эксергии вследствие необратимого теплообмена в испарителе:
dи = Δ e5-1 - Δ eн2-1 = 101,041 – 117,2 = 16,159 кДж/кг
24. КПД испарителя:
ηи = (Δ eн2-1 / Δ e5-1) 100% =(101,041 /117,2)100 = 86,21 %
Эксергетический баланс установки приведен в табл. 2.
Таблица 2
Параметр | кДж/кг | % |
Подвод эксергии | ||
Электроэнергия на компрессор eвх | 343,1 | |
Артезианская вода | 49,308 | |
Всего | 392,408 | 100% |
Отвод эксергии | ||
Электромеханические потери в компрессоре dэм | 34,31 | 8,74% |
Внутренние потери в компрессоре dкм | 49,56 | 12,63% |
Потери эксергии в конденсаторе dк | 88,409 | 22,53% |
Потери эксергии в дросселе dдр | 86,77 | 22,11% |
Потери в испарителе dи | 16,159 | 4,12% |
Эксергетическая холодопроизводительность е0 | 117,2 | 29,87% |
Всего | 392,408 | 100% |
Коэффициент полезного действия теплообменно-дроссельной части установки:
ηт.д = (e0 /( Δe2-1)) 100% = (117,2 /259,233)100% = 45,21%
Полный КПД установки:
η' = e0 /( eвх) 100% = (117,2 /343,1)100% = 34,2%
4.Подбор оборудования.
Холодопроизводительность установки:
Q0 = 20,93 кВт
Подбор оборудования ведем по программе Select 7/39780
Основные преимущества хладагента аммиак (R717) обусловлены тем, что он:
- обладает термодинамическими и теплофизическими характеристиками, позволяющими получать высокий КПД в холодильных установках;
- химически нейтрален по отношению к большинству конструкционных материалов холодильных установок, за исключением меди и сплавов на ее основе;
- не растворяется в смазочных маслах, применяемых в конструкциях холодильных установок, не чувствителен к влаге и легко обнаруживается в случае утечки;
- не способствует созданию парникового эффекта;
- имеет невысокую стоимость (не более 2200 рублей за тонну) и легко доступен на рынке.
Вместе с тем у аммиака есть ряд серьезных недостатков. В частности, это вещество:
- обладает высокой токсичностью (считается, что предельно допустимая концентрация аммиака в рабочих помещениях должна быть не выше 20 мг/м3, однако даже при более слабой концентрации характерный запах аммиака в случае его появления вызывает сильную панику; при более высоких концентрациях появляются серьезные затруднения дыхания вплоть до удушья; смертельная концентрация аммиака - 30 г/м3);
- является взрывоопасным (при концентрации в воздухе 200-300 г/м3 возникает угроза самопроизвольного взрыва; температура самовоспламенения равна 650 °С);
- создает опасность ожогов при растворении в воде, поскольку этот процесс
сопровождается выделением значительного количества тепла;
- имеет высокую температуру нагнетания при сжатии в холодильных компрессорах.
Поэтому его могут заменять фреонами. В качестве замены для подбора оборудования применим R600a. Краткие физические свойства приведены в таблице 3.
Таблица 3
Основные физические свойства R134a в сравнении с R717 | ||
Параметр | R717 | R600a |
Нормальная температура кипения (p=0,1 МПа), оС | -33,4 | -12 |
Температура замерзания, оС | -77,7 | -159 |
Критическая температура, оС | 133 | 135 |
Критическое давление, МПа | 11,4 | 3,65 |
Давление всасывания при -15 оС, МПа | 0,273 | 0,089 |
Потенциал разрушения озона | 0 | 0 |
Расчет компрессоров заключается в определение необходимого действительного объема, описываемого поршнями в единицу времени
, м3/с ,
где V0 — теоретический объем, м3/с ,определяется из выражения(6.6);
l — коэффициент подачи компрессора, он рассчитывается:
l=lс l¶р lt lпл .
выбираем несколько однотипных компрессоров, включенных параллельно в схему установки и обеспечивающих суммарную подачу. Общее число компрессоров с учетом одного резервного определится:
К = Краб + 1 .
Объем, описываемый поршнями в единицу времени одним компрессором
, м3 /с.
Расход воды (воздуха) на конденсатор:
, м3/с,
где Cp — теплоёмкость воды (воздуха), кДж/(кг×К); — нагрев воды (воздуха) в конденсаторе.
Перепад давления Δp = 0,5 бар
Переохлаждение ΔT = 4 оС
Результаты расчета и схема установки с приборами измерения представлены в приложениях.
Заключение.
В ходе выполнения семестрового задания мной была рассчитана холодильная установка, определены основные энергетические параметры.
Рекомендации:
1. Заменить хладагент R717 на R600a
2. Разбить установку на несколько для увеличения надежности.
Литература
1. Везиришвили О.Ш., Меладзе Н.В. Энергосберегающие теплонасосные системы тепло – и хладоснабжения. — М.:Изд-во МЭИ, 1994. — 160с., ил.
2. Горбенко В.И., Юртаев М.А. Расчёт парожидкостных компрессионных одноступенчатых холодильных и теплонасосных установок. - Ч: Изд. ЮУрГУ, 2002. — 38 с., ил.
3. Горбенко В.И. Конспект лекций. - Ч: Изд. ЮУрГУ, 2002. — 38 с., ил.
4. Курылев Е.С., Герасимов Н.А. Холодильные установки. – Л.: Машиностроение, 1980. – 622 с., ил.
5. Мартынов А.В. Установки для трансформации тепла и охлаждения. – М.: Энергоиздат, 1989. – 320 с., ил.
6. Соколов Е.Я., Бродянский В.М. Энергетические основы трансформации тепла и процессов охлаждения. – М.: Энергоиздат, 1981. – 320 с., ил.
7. Теплофизические основы получения искусственного холода: Справочник серии «Холодильная техника». — М.: Пищевая промышленность, 1980. — 232с.
8. Холодильные компрессоры: Справочник серии «Холодильная техника». — М.: Легкая и пищевая промышленность, 1981. — 280с.
9. Холодильные машины: Справочник серии «Холодильная техника». — М.: Легкая и пищевая промышленность, 1982. — 223с.
10. Щербин Н.К., Гринберг Я.И. Холодильные станции и установки. — М.: Химия, 1979. — 376с.