CC BY
80
УДК 621.56
Описание расчетной модели установки для получения «бинарного льда» методом вакуумно-выпарной кристаллизации
Канд. техн. наук Круглов А.А. al-x-kru@yandex.ru Тазитдинов Р.Р. t.r_92@mail.ru Университет ИТМО 191002, Россия, Санкт-Петербург, ул. Ломоносова, 9
Рассмотрен вариант установки для получения бинарного льда методом вакуумно-выпарной кристаллизации. Выделены основные процессы в баке-кристаллизаторе, существенно влияющие на работу установки. Расчетная модель включает следующие этапы: составление тепловых и массовых балансов установки; расчёт системы откачки водяного пара и описание процесса распыления; расчёт процесса кристаллизации капель; расчет геометрических параметров бака-кристаллизатора. Составлена расчетная модель, позволяющая выполнять расчёты процессов, анализировать влияние технологических характеристик и определять геометрические параметры установки. Выполнен расчет технических характеристик установки для производства бинарного льда холодопроизводительностью 12 кВт. Ключевые слова: бинарный лёд, вакуумно-выпарная кристаллизация, тепловой баланс, массовый баланс, схема установки, бак-кристаллизатор.
DOI: 10.17586/2310-1148-2016-9-3-26-35
Description of analysis model of the plant for producing «ice slurry» method of vacuum-evaporating crystallization
Ph.D. Kruglov A.A. al-x-kru@yandex.ru Tazitdinov R.R. t.r_92@mail.ru ITMO University 191002, Russia, St. Petersburg, Lomonosov str., 9
The variant ofplant for producing ice slurry of method of vacuum-evaporation crystallization. The main processes highlighted in the tank-crystallizer, which significantly affect the performance of the plant. The analytic model includes the following stages: generation of heatbalance and material balanceof the system, calculation of the vacuum system water vapor and the description of the atomization process; a calculating process of crystallization drops; the calculation of the geometric parameters of the tank- crystallizer. Compiled computational model to calculations process, to analyze the impact of technological characteristics and to calculations the geometric parameters of the plant. The calculation of the technical characteristics of installations for the production of ice slurry cooling capacity of 12 kW.
Keywords: ice slurry, vacuum-evaporation crystallization, heat balance, mass balance, plant layout, tank-crystallizer.
Введение
Подписание Монреальского протокола по веществам, разрушающим озоновый слой Земли, и Киотского протокола о стабилизации выбросов парниковых газов в атмосферу оказали влияние на развитие холодильной техники в сторону применения природных хладагентов и снижения заправочной ёмкости систем с синтетическими хладагентами.
В 1972 году российский ученый В.А. Чижиков предложил использовать двухфазный хладоноситель «бинарный лёд». Бинарный лёд состоит из мелких кристаллов льда, размер которых не превышает 500 мкм, в смеси с водой или водным раствором [1]. Доля льда в смеси может варьироваться от 10 до 50 % [2].
Применение бинарного льда в качестве хладоносителя со стабильной температурой, меняющего своё агрегатное состояние в процессе охлаждения, наиболее экономично в случаях неравномерных тепловых нагрузок при околонулевых температурах. Бинарный лёд обладает высокой охлаждающей способностью, создаваемой многочисленными кристаллами, которые обволакивают продукт со всех сторон при контактном охлаждении, тем самым предотвращают его контакт с воздухом, обеспечивая при этом быстрое охлаждение и увеличивая срок
хранения охлаждаемой продукции [3]. Бинарная смесь с массовой долей льда 30 % потенциально может улучшить эффективность системы за счёт уменьшения объёма циркулирующей воды более чем на 80 %, или же увеличить охлаждающую способность в 6 раз при тех же тепловых нагрузках и мощностях на перекачивание хладоносителя. Опубликованные результаты исследований [4], свидетельствуют об увеличении коэффициента теплопередачи на 20 % для пластинчатых теплообменников и на 50 % для трубчато-ребристых теплообменников при использовании бинарного льда, что позволяет сократить расходы на теплообменное оборудование, трубопроводы, насосы и арматуру. Область применения и популярность использования бинарного льда с каждым годом растёт в странах Европы, в Америке и Японии [5, 6]. Это объясняется его теплофизическими и экологическими свойствами.
При производстве бинарного льда наиболее важным процессом является получение ледяных кристаллов с определёнными размером, формой и шероховатостью, которые влияют на свойства бинарной смеси в целом.
Процесс генерации ледяных кристаллов состоит из трёх основных этапов: переохлаждения жидкости; образования ядра и роста кристаллов. Этап переохлаждения жидкости представляет собой охлаждение ниже температуры кристаллизации без изменения фазового состояния. Затем, после переохлаждения жидкости, происходит гетерогенное образование ядер, которому способствуют чужеродные частицы и шероховатость поверхностей теплообменника. После образования ядра происходит рост кристаллов [7].
После этих трёх основных этапов ледяные кристаллы могут подвергнуться другим процессам, таким как истощение, агломерация и созревание [8]. Эти процессы влияют на форму и размеры кристаллов льда, а соответственно, и на свойства бинарной смеси в целом. При смешивании бинарного льда происходит истощение ледяных кристаллов за счёт напряжения, возникающего при столкновении кристаллов с твердыми поверхностями. Агломерация - это процесс слипания кристаллов друг с другом, с образованием при этом ледяных кристаллов с большим размером. Процесс агломерации может наблюдаться как при получении, так и при хранении (без перемешивания). При длительном хранении бинарного льда возникает процесс созревания, то есть мелкие кристаллы распадаются, и тогда основу бинарной смеси составляют большие частицы льда.
Как правило, мелкие кристаллы сферической формы легко перекачиваются по трубам и теплообменникам, поэтому процессы измельчения и сглаживания углов кристаллов занимают важную позицию в получении бинарной смеси. Исследования Kasza и Hayashi [9] показали, что полученные из пресной воды кристаллы имеют дендритную форму, то есть имеют неровную, шероховатую поверхность, это приводит к плохой текучести и быстрой агломерации бинарной смеси. Также им удалось обнаружить, что введение депрессантов в воду позволяет сглаживать поверхность кристаллов химическим путём.
Получение искусственного льда методом капельной кристаллизации является одним из высокопроизводительных способов. Высокая производительность данного способа достигается за счёт увеличения поверхности контакта воды при её распылении.
Применяемые в настоящее время парокомпрессионные генераторы бинарного льда стоят дороже водоохлаждающих холодильных машин. Новые способы получения бинарного льда должны составить конкуренцию существующим. Производство бинарного льда с помощью вакуумно-выпарных установок имеет ряд преимуществ: в качестве хладагента в вакуумных установках применяется вода; система находится под давлением ниже атмосферного; коэффициент эффективности вакуумно-выпарных установок выше, чем у традиционных способов получения. Таким образом, разработка новых установок, применение и исследование процессов в вакуумных установках для получения бинарного льда является актуальной темой [10].
Описание установки для получения бинарного льда
В качестве объекта рассматривается холодильная установка для получения бинарного льда методом вакуумно-выпарной кристаллизации.
Установка состоит из следующих основных элементов (см. рис. 1): вакуумного бака-кристаллизатора - 1, вакуумного насоса - 2, форвакуумного насоса - 3, конденсатора - 4, рециркуляционного насоса - 5, циркуляционного насоса - 6 и теплообменного аппарата - 7.
Процесс получения бинарного льда в вакуумно-кристаллизационной установке осуществляется в следующей последовательности. С помощью вспомогательного вакуумного насоса создаётся предварительное разряжение с остаточным давлением около 1100-2000 Па. После чего включается основной вакуумный насос, который создаёт рабочее давление в баке-кристаллизаторе на уровне 200-400 Па (абс.). Разрежение вспомогательным насосом необходимо для запуска основного вакуумного насоса, работа которого от атмосферного давления невозможна.
Рис. 1. Принципиальная схема вакуумно-выпарной установки для получения бинарного льда
Внутренний объем вакуумного бака-кристаллизатора для упрощения расчётов условно можно разделить на три основные зоны:
• зона распыления;
• зона кристаллизации;
• зона бинарной смеси.
В зоне распыления в верхней части бака-кристаллизатора происходит непрерывное распыление воды через форсунки. Распылительные устройства дают так называемый факел распыления, который представляет собой множество взаимодействующих между собой капель различного размера. Наиболее вероятный размер капель зависит от конструктивных особенностей форсунок и параметров процесса распыления.
Распыляемая жидкость обладает большой поверхностью испарения, вследствие чего, в процессе испарения при пониженном давлении часть жидкости в распыляемом факеле превращается в водяной пар. Самая мелкая часть капель в распыляемом факеле взвешивается потоками водяного пара, образуя при этом туман. Более крупные капли под действием силы тяжести выпадают из факела. Таким образом, в зоне распыления находятся водяной пар, туман и падающие капли.
Падающие из факела капли попадают в зону кристаллизации. В этой зоне капля попадает в псевдоожиженную среду, где на неё действует сила выталкивания, создаваемая потоками откачиваемого вакуумным насосом водяного пара, и сила тяжести. Эти две силы влияют на время нахождения капли в зоне кристаллизации и геометрические размеры бака-кристаллизатора. Во время нахождения капель в зоне кристаллизации, часть их переходит в паровую фазу, а часть кристаллизуется, образуя при этом падающие ледяные кристаллы льда.
Полученные кристаллы льда оседают в зоне бинарной смеси. В ней собираются неиспарившаяся жидкость и мелкодисперсные кристаллы льда. За счёт разности плотностей воды и льда, в верхней части этой зоны образуется смесь воды с высокой концентрацией кристаллов - бинарный лёд, откуда он забирается и подаётся циркуляционным насосом в пластинчатый теплообменник для охлаждения продукта. В нижней части скапливается вода с околонулевой температурой, которая непрерывно с помощью рециркуляционного насоса, обеспечивающего таким образом перемешивание смеси в объёме, что не даёт слипаться кристаллам, подаётся в форсунки (в верхнюю часть бака-кристаллизатора).
Откачиваемый вакуумным насосом водяной пар поступает в пластинчатый конденсатор 5, где он конденсируется. Конденсат возвращается в бак-кристаллизатор.
Таким образом, при описании расчётной модели вакуумно-выпарной установки можно выделить следующие этапы:
• составление тепловых и массовых балансов установки;
• расчёт системы откачки водяного пара и описание процесса распыления;
• расчёт процесса кристаллизации капель;
• расчет геометрических параметров бака-кристаллизатора.
Тепловые и массовые балансы
Уравнения тепловых и массовых балансов, являются важнейшими для описания математической модели холодильной установки. Они позволяют вести технологический расчёт и анализировать теплообменные процессы, влияющие на установку.
На рис. 2 изображена схема для составления тепловых и массовых балансов вакуумно-испарительной установки для получения бинарного льда.
Рис. 2. Схема расчёта балансов вакуумно-выпарной установки для получения бинарного льда
Для составления систем уравнений тепловых и массовых балансов, в схеме, изображённой на рис. 2, можно выделить четыре основных группы:
• балансы в баке-кристаллизаторе;
• балансы процесса кристаллизации капель;
• балансы в теплообменнике;
• балансы в конденсаторе.
Тепловые и массовые балансы в баке-кристаллизаторе (см. рис. 3) можно представить следующей системой уравнений.
Рис. 3. Схема расчёта тепловых и массовых балансов бака-кристаллизатора
6см • 1см - 6см • 1в + • 1п - • ¿к + бр • ¿р - вр • Ьр2 = 0;
ср'1р~ср' 1р2 + мр = 0; &см = (Сл + £в);
Ср = (Сл + Сп).
Принимаем:
С„ =
где Ссм - расход бинарной смеси, кг/с; 1СМ - энтальпия на входе в теплообменник, кДж/кг; /в - энтальпия бинарной смеси на выходе из теплообменника, кДж/кг; - расход водяного пара, кг/с; ¿„ - энтальпия водяного пара, кДж/кг; йк - расход воды после конденсатора, кг/с; ¿к - энтальпия воды после конденсатора, кДж/кг; Ср -расход воды перед рециркуляционным насосом, кг/с; 1р - энтальпия воды перед рециркуляционным насосом, кДж/кг; ¿р2 - энтальпия воды после рециркуляционного насоса, кДж/кг; Сл - расход льда в бинарной смеси, кг/с; Св - расход воды в бинарной смеси, кг/с; - мощность рециркуляционного насоса, кВт.
Тепловые и массовые балансы процесса кристаллизации капель (см. рис. 4) можно представить следующей системой уравнений.
^кр
ч>
Gn\ ¿п
<?л; ¿л
Рис. 4. Схема расчёта тепловых и массовых балансов процесса кристаллизации капель
Gp • iP2 — GjI-ij,-Gn-in + Qn — Q^, — Qm = 0; Qxp = (G„ + GjÏ) • ip2 + (Gn + Gj,) • ieQ + (Gn + Gj,) • iл) Qn = G„• r„,
где ¿л- энтальпия переохлажденного льда, кДж/кг; ie0- энтальпия воды при температуре фазового перехода, кДж/кг; гп - теплота парообразования, кДж/кг; Qn - тепловая нагрузка при испарении воды, кВт; QKp - тепловая нагрузка при кристаллизации воды, кВт; Qm - теплопритоки через стенки бака-кристаллизатора, кВт.
Тепловые и массовые балансы в теплообменном аппарате (см. рис. 5) можно представить следующей системой уравнений.
^пр' ^пр
ÇCM |
1 1 Çnp
гпр> ^пр2
Рис. 5. Схема расчёта тепловых и массовых балансов в теплообменнике
6Пр ' Iпр 6пр • 1Пр2 "Ь Носн + (¡см 0,
Gcм ' 1см Gcм • /в + Носн + 0,
0.см = Gл ' 1л — Сл • ¿в0 + Сл • + Ссм • ¿в0 — Ссм • ¿в;
0.Щ) &пр • Ьпр СПр ' 1Пр21 0 см <1пр>
где С„р - расход охлаждаемого продукта, кг/с; Iпр - энтальпия продукта на входе, кДж/кг; 1пр2 - энтальпия продукта на выходе, кДж/кг; (}см - тепловая нагрузка бинарной смеси, кВт; (¿пр - тепловая нагрузка от продукта, кВт.
Тепловые и массовые балансы в конденсаторе (см. рис. 6) можно представить следующей системой уравнений.
Рис. 6 Схема расчёта тепловых и массовых балансов в конденсаторе
^ ' К/ ~ • ¿ц,2 + Мвн + (}кд = 0; сп • ц - • ¿к + К» + Qw = 0;
Окд = ^л • гп> (¿ы = Qкд>
где С^,- расход охлаждающей воды в конденсаторе, кг/с; ¿ц,- энтальпия охлаждающей воды на входе в конденсатор, кДж/кг; - энтальпия охлаждающей воды на выходе из конденсатора, кДж/кг; (}кд - тепловая нагрузка от конденсации водяного пара, кВт; (¿у, - тепловая нагрузка, отводимая охлаждающей водой, кВт; Ывн - мощность вакуумного насоса, кВт.
Расчет системы откачки водяного пара
Испарившийся в баке-кристаллизаторе пар по вакуумной магистрали длиной I и диаметромс£ откачивается вакуумным насосом.
Эффективная скорость откачки установки определяется из следующего уравнения [11]:
_ Б'и
Яэф -
где 5 - номинальная скорость откачки вакуумного насоса, м3/с; и - пропускная способность вакуумной коммуникации, м3/с;
Пропускная способность для трубопровода круглого сечения при вязкостном режиме течения находится по формуле [12]:
7Г • С*4 • + Р2)
и =-—-—,
128-т? -I -2
где (I - диаметр трубопровода, м; Рг и Р2- давление в начале и в конце трубопровода, Па; Г] - динамическая вязкость откачиваемой среды, Нс/м2; I - длина трубопровода, м.
Удельная скорость откачки (эффективная скорость откачки на единицу площади сечения вакуумной полости), м3/(с-м2);
Рпф
где F„ф - площадь поверхности факела распыляемой воды, м2;
Fлф = 4лт2п
где п - количество капель в распыляемом факеле.
Количество капель в распыляемом факеле находится по формуле:
С.
п„ = —,
где тк- расход воды в одной капле, кг/с.
Расчет процесса кристаллизации капель в вакууме
Маринюком Б.Т. был изучен теплофизический процесс замерзания капли воды в вакуумной среде [13]. При описании процессов принимались следующие допущения:
• начальная температура капли воды равна нулю;
• теплофизические свойства воды и водного льда постоянны;
• температура фронта кристаллизации равна температуре стабильного равновесия между жидкой и твердой фазами воды;
• фронт кристаллизации представляет собой сферу с переменным радиусом;
• распределение температуры в капле подчиняется закону стационарного распределения;
с*
• удельная скорость откачки вакуумного насоса Ь постоянна;
Рис. 7 Модель промерзания капли в вакууме
Фронт промерзания капли в сферической системе координат (см. рис. 7) распространяется от поверхности к центру капли по координате г. Радиус границы фазового превращения обозначен символом ц; Я - радиус капли. Соответственно, толщина слоя льда в капле £ меняется со временем т.
Граничные условия «лёд-пар» можно записать в виде:
= 5* • Ь* • р,
дТ
г=Я
где Ял - коэффициент теплопроводности льда, Вт/(м-К); 5* - удельная скорость откачки (эффективная скорость откачки на единицу площади сечения вакуумной полости), м3/(с-м2); Ь* - теплота парообразования, Дж/кг; р - упругость насыщенных паров воды, Па.
На границе фазового превращения «лёд-вода» условия теплового взаимодействия можно записать в виде:
дТ
г дг>
= £ • — Рл. Г=7} ОТ
где Ь - теплота фазового перехода воды в лёд, Дж/кг; рл - плотность льда, кг/м3.
Согласно рекомендациям Л.С. Лейбензона, распределение температуры в сферическом слое твёрдой фазы примем как стационарное:
где Тф - температура фазового перехода воды, К; Т„ - начальная температура капли, К.
В результате преобразования перечисленных выражений Маринюком Б.Т. было получено уравнение [13], решение которого позволяет определить изменение времени замерзания т капли радиусом Я в вакууме:
Ьрл /V R2t) rj3\_S*L*nl_ 8940
Кт(Д -rj)\3R 6 2 J Ra \ т _ Ь£л_ (if я* _ д^ч >
ф Ллйт V 3 6 2 )'
где ц - молекулярная масса воды, г/моль; Яг - универсальная газовая постоянная.
Расчет геометрических параметров бака-кристаллизатора
Для изготовления бака-кристаллизатора требуется определить его размеры, которые зависят от процессов, происходящих в баке. Основные размеры бака определяются его диаметром и высотой зоны кристаллизации, представляющей из себя псевдоожиженный слой.
Уравнение скорости осаждения кристаллов в псевдоожиженном слое [14]:
Аг • £4'75 • т9„
Шг =
(18 + 0,59 • л/Аг • е4-75) • dK'
где Аг - критерий Архимеда; 8 - порозность; i9gn - кинематическая вязкость водяного пара, м2/с; dK - диаметр капли, м.
Критерий Архимеда можно вычислить по формуле:
9 ' d3 • (pm - Pep)
Аг =
• Pep
Уравнение скорости восходящих потоков водяного пара:
G„
где - расход водяного пара, кг/с; И - радиус сечения бака-кристаллизатора в зоне кристаллизации, м2; рп -плотность водяных паров, кг/м3;
Зависимость высоты зоны кристаллизации от скорости осаждения ледяных кристаллов можно описать с помощью уравнения для высоты псевдоожиженного слоя, м [15]:
н=
Т
-S
Применение расчетной модели
Расчетная модель позволяет определить основные характеристики вакуумно-выпарной установки для получения бинарного льда. Результаты расчетов на примере установки с холодопроизводительностью 12 кВт приведены в таблице 1.
Установка оснащена вакуумным насосом 2ДВН-1500 с номинальной скоростью откачки 1500 л/с, в качестве форвакуумного насоса первой ступени применяется насос R13ВЗ-180. С учётом пропускной способности вакуумной магистрали длиной 5 метров, диаметром 250 миллиметров и средним давлением 50 Па эффективная скорость вакуумного насоса составляет 1474 л/с.
Установка производит бинарный лёд с массовой долей льда 30%, и общим расходом смеси 401 кг/ч. За один час работы установка производит 120 килограмм чистого льда, с размером ледяных частиц 1,5 мм.
Таблица 1
Расчётные характеристики установки для производства бинарного льда
Характеристики вакуумно-испарительной установки
Холодопроизводительность 12 кВт
Производительность (по льду) 2,8 т/сут
Содержание кристаллов льда 30 %
Диаметр ледяных кристаллов 1,5 мм
Потребляемая мощность 9,5 кВт
Бак-кристаллизатор
Толщина стали 15 мм
Толщина теплоизоляции (пенополиуретан) 50 мм
Высота 2,57 м
Диаметр (основной) 1 м
Общая площадь поверхности 8,05 2 м
Вывод
Рассмотрена модель установки для получения бинарного льда методом вакуумно-выпарной
кристаллизации, позволяющая выполнять расчёты процессов, анализировать влияние технологических
характеристик и определять геометрическиепараметры установки.
Литература
1. Цветков О.Б., Лаптев Ю.А. Айс-сларри и однофазные хладоносители // Холодильная техника. 2004. N 3. С.2-4.
2. Маринюк Б.Т., Сусликов Д.В., Ермолаев А.Е. Экологически чистые методы получения водного льда // Холодильный бизнес. 2008. №2. С. 38-40.
3. Белозеров Г.А., Медникова Н.М., Пытченко В.П., Серова Е.Н. Исследование процесса охлаждения рыбы с использованием бинарного льда // Холодильная техника. 2012. N 6. С. 37-41.
4. Paul J. District cooling with binary ice. Design and execution of district cooling systems with a liquid, pumpable ice slurry «binary iceTechnicalUniversityofDenmark.
5. Цой А.П., Круглов А.А., Тазитдинов Р.Р. Основные типы аккумуляторов холода с использованием бинарного льда // VI Международная научно-техническая конференция "Казахстан-Холод 2016" (1-2 марта 2016 г.): сборник докладов конференции. 2016. С. 144-148.
6. Wang MJ; Goldstein V (2003) Ice slurry: advanced fish chilling and preservation technology [Текст]. American Fisheries Society Symposium 38:379-386
7. ДоронинЮ.П. Физика океана. - СПб.: изд. РГГМУ. 2000. - 340 с.
8. Миненков В.В., Зимин А.В., Хмельнюк М.Г. Влияние примесей на процессы получения и применения бинарного льда // Холодильная техника и технологии. 2013. № 4 (144). С. 24-30.
9. Hayashi K, Kasza KE. Ice slurry cooling research: micro scale study of ice particles characteristics, role of freezing point depressant, and influence on slurry fluidity. In: Proceedingsofthe ASHRAE meeting. Atlanta, Georgia; 2001.
10. Тазитдинов Р.Р., Круглов А.А. Анализ существующих методов получения бинарного льда // Альманах научных работ молодых ученых Университета ИТМО. 2016. Т. 5. С. 67-69
11. Кузнецов В.И. Механические вакуумные насосы. - М.: Под редакцией М.И. Меньшикова, 1959. - 279 с.
12. РозановЛ.Н. Вакуумные машины и установки. Л.: Машиностроение, 1975. - 302 с.
13. Маринюк Б.Т. Расчет теплообмена в аппаратах и системах низкотемпературной техники. - М.: Машиностроение, 2015. -272 с.
14. Айнштейн М.К., Захаров М.К., Носов Г.А., Захаренко В.В., Зиновкина Т.В., Таран А.Л., Костанан А.Е. Общий курс процессов и аппаратов химической технологии. - М.: Университетская книга; Лотос; Физматкнига, 2006. Кн. 1. - 912 с.:ил.
15. Айнштейн М.К., Захаров М.К., Носов Г.А., Захаренко В.В., Зиновкина Т.В., Таран А.Л., Костанан А.Е. Общий курс процессов и аппаратов химической технологии. - М.: Университетская книга; Лотос; Физматкнига, 2006. Кн. 2. - 872 с.:ил.
References
1. Tsvetkov O.B., Laptev Yu.A. Ais-slarri i odnofaznye khladonositeli // Refrigerating engineering. 2004. N 3. S.2-4.
2. Marinyuk B.T., Suslikov D.V., Ermolaev A.E. Ekologicheski chistye metody polucheniya vodnogo l'da // Kholodil'nyi biznes. 2008. №2. S. 38-40.
3. Belozerov G.A., Mednikova N.M., Pytchenko V.P., Serova E.N. Issledovanie protsessa okhlazhdeniya ryby s ispol'zovaniem binarnogo l'da // Refrigerating engineering. 2012. N 6. S. 37-41.
4. Paul J. District cooling with binary ice. Design and execution of district cooling systems with a liquid, pumpable ice slurry «binary iceTechnicalUniversityofDenmark.
5. Tsoi A.P., Kruglov A.A., Tazitdinov R.R. Osnovnye tipy akkumulyatorov kholoda s ispol'zovaniem binarnogo l'da // VI Mezhdunarodnaya nauchno-tekhnicheskaya konferentsiya "Kazakhstan-Kholod 2016" (1-2 marta 2016 g.): sbornik dokladov konferentsii. 2016. S. 144-148.
6. Wang MJ; Goldstein V (2003) Ice slurry: advanced fish chilling and preservation technology [Tekst]. American Fisheries Society Symposium 38:379-386
7. Doronin Yu.P. Fizika okeana. - SPb.: izd. RGGMU. 2000. - 340 s.
8. Minenkov V.V., Zimin A.V., Khmel'nyuk M.G. Vliyanie primesei na protsessy polucheniya i primeneniya binarnogo l'da // Refrigerating equipment and technologies. 2013. № 4 (144). S. 24-30.
9. Hayashi K, Kasza KE. Ice slurry cooling research: micro scale study of ice particles characteristics, role of freezing point depressant, and influence on slurry fluidity. In: Proceedingsofthe ASHRAE meeting. Atlanta, Georgia; 2001.
10. Tazitdinov R.R., Kruglov A.A. Analiz sushchestvuyushchikh metodov polucheniya binarnogo l'da // Al'manakh nauchnykh rabot molodykh uchenykh Universiteta ITMO. 2016. T. 5. S. 67-69
11. Kuznetsov V.I. Mekhanicheskie vakuumnye nasosy. - M.: Pod redaktsiei M.I. Men'shikova, 1959. - 279 s.
12. Rozanov L.N. Vakuumnye mashiny i ustanovki. L.: Mashinostroenie, 1975. - 302 s.
13. Marinyuk B.T. Raschet teploobmena v apparatakh i sistemakh nizkotemperaturnoi tekhniki. - M.: Mashinostroenie, 2015. -272 s.
14. Ainshtein M.K., Zakharov M.K., Nosov G.A., Zakharenko V.V., Zinovkina T.V., Taran A.L., Kostanan A.E. Obshchii kurs protsessov i apparatov khimicheskoi tekhnologii. - M.: Universitetskaya kniga; Lotos; Fizmatkniga, 2006. Kn. 1. - 912 s.:il.
15. Ainshtein M.K., Zakharov M.K., Nosov G.A., Zakharenko V.V., Zinovkina T.V., Taran A.L., Kostanan A.E. Obshchii kurs protsessov i apparatov khimicheskoi tekhnologii. - M.: Universitetskaya kniga; Lotos; Fizmatkniga, 2006. Kn. 2. - 872 s.:il.
Статья поступила в редакцию 29.09.2016 г.
