Расчет и проектирование двухкорпусной выпарной установки

Размещено на http://www.allbest.ru/

21

Размещено на http://www.allbest.ru/

Введение

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

Выпаривание применяется:

1. для концентрации растворов твердых веществ;

2. для концентрации растворов высококипящих жидкостей;

3. для выделения растворителя в чистом виде.

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

Выпаривание производят несколькими способами:

а) под вакуумом;

б) под высоким (избыточным) давлением;

в) под атмосферным давлением.

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

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

1. Описание технологической схемы

Исходный разбавленный раствор из промежуточной емкости 1 центробежным насосом 2 подается в теплообменник 3, где подогревается до температуры, близкой к температуре кипения, а затем - в первый корпус 4 выпарной установки. Предварительный подогрев раствора повышает интенсивность кипения в выпарном аппарате. Первый корпус обогревается свежим водяным паром. Вторичный пар, образующийся при концентрировании раствора квасного сусла в первом корпусе, направляется в качестве греющего во второй корпус 5. Сюда же поступает частично сконцентрированный раствор сусла из первого корпуса. Аналогично третий корпус 6 обогревается вторичным паром второго и в нем производится концетрирование раствора, поступившего из второго корпуса.

Самопроизвольный переток раствора и вторичного пара в следующие корпуса возможен благодаря общему перепаду давлений, возникающему в результате создания вакуума конденсацией вторичного пара последнего корпуса в барометрическом конденсаторе смешения 7, где заданное давление поддерживается подачей охлаждающей воды и отсосом неконденсирующихся газов вакуум - насосом 8. Смесь охлаждающей воды и конденсата выводится из конденсатора при помощи барометрической трубы с гидрозатвором 9. Образующийся во втором корпусе концентрированный раствор сусла с массовой долей 69% центробежным насосом 10 подается в промежуточную емкость упаренного раствора 11. Конденсат греющих паров из выпарных аппаратов выводиться с помощью конденсатоотводчиков 12.

2. Расчет выпарной установки

Рассчитать и спроектировать двухкорпусную выпарную установку непрерывного действия для сгущения квасного сусла производительностью (G_н) 3000 кг. Начальное сусло содержит (Х_н) 12% сухих веществ, а упаренное (Х_к) 69%. Сусло поступает на выпаривание подогретым до температуры кипения. Давление греющего пара - избыточное (Р_г1) 0,2 МПа, давление вторичного пара в барометрическом конденсаторе (Р_бк) - абсолютное 0,08МПа.

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

2.1 Определение поверхности теплопередачи выпарных аппаратов

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

F=Q/K Дt_п , (1)

насос вакуум конденсатор выпарной

где

Q - тепловая нагрузка, кВт

K - коэффициент теплопередачи, Вт/(м²*К)

Дt_п - полезная разность температур, ?С.

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

Первое приближение.

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

W = G_н(1-Х_н/Х_к) (2)

переведем кг/ч в кг/с: 3000/3600=0,83 кг/с,

Получим:

W=0,83* (1-12/69)=0,67 кг/с.

2.1.1 Расчет концентраций упариваемого раствора

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

w₁: w₂ = 1,0: 1,1

тогда:

w₁ = 1,0*W/1,0+1,1 = 1,0*0,67/2,1 =0,32 кг/с

w₂ = 1,1*W/1,0+1,1 = 1,1*0,67/2,1 = 0,35 кг/с

Далее рассчитывают концентрации растворов в корпусах:

Х₁ = G_н * x_н = 0,83*0,12 = 0,19 (19%)

G_н - w₁ 0,83-0,32

Х₂ = G_н * Х_н = 0,830 * 0,12 = 0,62 (62%)

G_н-w₁-w₂ 0,83-0,32-0,35

Концентрация раствора в последнем корпусе Х₂ соответствует заданной концентрации упаренного раствора Х_к.

2.1.2 Определение температур кипения растворов

Так как давление греющего избыточное в пересчете на нормальные условия

Р_г1 = 2 ат. + 1 ат. = 3 ат. или 0,3 МПа

Общий перепад давлений в установке равен:

ДР_об = Р_г1 - Р_бк = 0,3 - 0,08 = 0,22 МПа

В первом приближении общий перепад давлений распределяют между корпусами поровну. Тогда давления греющих паров в корпусах равны:

Р_г1 = 0,3

Р_г2 = Р_г1 - Д Р_об/2 = 0,3 - 0,22/2 = 0,19Па

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

Р_бк = Р_г2 - ДР_об/2 = 0,19-0,22/2 = 0,08 МПа,

что соответствует заданной величине Р_бк.

По давлениям паров находим их температуры и энтальпии [1]:

Давление, МПа Температура, ?С Энтальпия, кДж/кг

Р_г1 = 0,3 t ₁ = 132,9 I₁ = 2730

Р_г2 = 0,19 t₂ = 118,0 I₂ = 2710

Р_бк = 0,08 t_бк = 93,0 I_бк = 2663

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

Изменение температуры кипения по высоте кипятильных труб происходит в следствие изменения гидростатического давления столба жидкости. Температуру кипения раствора в корпусе принимают в соответствующей температуре кипения в среднем слое жидкости. Таким образом, температура кипения раствора в корпусе, отличается от температуры греющего пара в последующем на сумму температурных потерь (УД) от температурной (Д^/), гидростатической (Д^//) и гидродинамической (Д^///) депрессий.

Гидродинамическая депрессия обусловлена потерей давления пара на преодоление гидравлических сопротивлений трубопроводов при переходе из корпуса в корпус. Обычно в расчетах принимают Д^/// = 1,0 - 1,5 град на корпус. Примем для каждого корпуса Д^/// = 1 град. Тогда температуры вторичных паров в корпусах равны:

t_вп1 = t_г2 + Д₁^/// = 118,0 + 1,0 = 119,0 ?С

t_вп2 = t_бк + Д₂^/// = 93,0 + 1,0 = 94,0 ? С

Сумма гидродинамический депрессий

УД^/// = Д₁^/// + Д₂^/// = 1,0 +1,0 = 2,0 ?С

По температурам вторичных паров определим их давления:

Температура, ?С Давление, МПа

t_вп1 =119.0 Р_вп1 = 0,14

t_вп2 = 94.0 Р_вп2 = 0,08

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

Р_ср = Р_вп + pgH (3)

(1-е),

где:

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

p - плотность кипящего раствора, кг/м³;

е - паронаполнение (объёмная доля пара в кипящем растворе), м³/м³.

Для выбора значения Н необходимо ориентировочно оценить поверхность теплопередачи выпарного аппарата (F_ор.) При кипении водных растворов можно принять удельную тепловую нагрузку аппаратов с естественной циркуляцией q = 20000-50000 Вт/м². примем q = 20 000 Вт/м². Тогда поверхность теплопередачи 1-го корпуса ориентировочно равна:

F_ор = Q/q = w₁r₁ / q = 0,32 * 2068 * 10³ / 20000 = 154,01 м²,

где

r - теплота парообразования вторичного пара [1].

Выпарные аппараты с восходящей пленкой с естественной циркуляцией состоят из кипятильных труб высотой 5 и 7 м при диаметре d_н = 38 мм и толщине стенки д_ст= 2 мм. Примем высоту кипятильных труб Н = 5 м.

При пузырьковом (ядерном) режиме кипения паронаполнение составляет е = = 0,4 - 0,6. Примем е = 0,5. Плотность водных растворов, в том числе раствора квасного сусла [3] при температуре 15?С и соответствующих концентрациях равна:

Концентрация, % Плотность, кг/м³

Х₁ = 19 p₁ = 1997

Х₂ = 62 p₂ =2318

При определении плотности растворов в корпусах пренебрегаем изменением ее с повышением температуры от 15 ?С до температуры кипения ввиду малого значения коэффициента объемного расширения и ориентировочно принятого значения е.

Давления в среднем слое кипятильных труб корпусов равны:

Р_1ср = Р_вп1 + p₁gH * (1 - е) = 19,0*10⁴ +1997 * 9,8 * 5 * (1-0,5) = 21,4 * 10⁴ Па

Р_2ср = р_вп2 + p₂gH * (1 - е) = 8,0 *10⁴ + 2318 * 9,8 * 5 * (1 - 0,5) = 10,8 * 10⁴ Па

Этим давлениям соответствуют следующие температуры кипения и теплоты испарения растворителя [1]:

Давление, МПа Температура, ?С Теплота испарения, кДж/кг

Р_1ср = 0,214 t_1cр = 120,9 r_1ср = 2204

Р_2ср = 0,108 t_2ср= 110,5 r_2ср= 2252

Определяем гидростатическую депрессию по корпусам

Д₁^// = t_1ср - t_вп1 = 120,9 - 109,7 = 11,2 ?С

Д₂^// = t_2ср - t_вп2 = 110,5 - 94,0 = 6,5 ?С

Сумма гидростатических депрессий равна:

УД^// = Д₁^// + Д₂^// = 11,2+6,5 = 18,7 ?С

Температурная депрессия Д^/ определяется по уравнению

Д^/ = 1,62 * 10^-2 *(Т² / r_вп) * Д ^/_атм, (4)

где

Т - температура паров в среднем слое кипятильных труб, ?С;

Д^/_атм - температурная депрессия при атмосферном давлении, ?С [2]

Находим значение Д^/ по корпусам:

Д₁^/ = 1,62 * 10^-2 * (120,9 + 273)² * 6,5 = 7,4 ?С

2204

Д₂^/ = 1,62 * 10^-2 * (110,5 + 273)² * 17,4 = 18,4 ?С

2252

Сумма температурных депрессий равна:

УД^/ =Д₁^/ + Д₂^/ = 7,4 + 18,4 = 25,8 ?С

Температуры кипения растворов в корпусах равны

t_к1 = t_г2 + Д₁^/ + Д₁^// + Д₁^///

t_к2 = t_бк + Д₂^/ + Д₂^// + Д₂^///

При расчете температуры кипения в пленочных выпарных аппаратах гидростатическую депрессию (Д^//) не учитывают. Тогда:

t_к1 = 118,0 + 7,4 +1 =126,4 ?С

t_к2 = 93,0 + 18,4 + 1 = 112,4 ?С

2.1.3 Расчет полезной разности температур

Общая полезная разность температур равна:

УДt_п = Дt_п1 + Дt_п2

Полезные разности температур по корпусам

Дt_п1 = t_г1 - t_к1 = 132,9 - 126,4 = 6,5 ?С

Дt_п2 = t_г2 - t_к2 = 118,0 - 112,4 = 5,6 ?С

УДt_п = 6,5 + 5,6 = 12,1 ?С

Проверим полезную разность температур:

УДt_п = t_г1 - t_бк - (УД^/ + УД^///) = 132,9 - 93,0 - (25,8 + 2) = 12,1 ?С

2.1.4 Определение тепловых нагрузок

Расход греющего пара в 1-й корпус, производительность каждого корпуса по выпаренной воде и тепловые нагрузки по корпусам определим путем совместного решения уравнений тепловых балансов по корпусам и уравнения баланса по воде для всей установки:

Q₁ = D(I_г1 - i₁₎ = 1,03 [G_нс_н(t_к1 - t_н) + w₁(I_вп1 - с_вt_к1) + Q_конц] (5)

Q₂ = w₁(I_г2 - i₂₎ = 1,03 [(G_н - w₁) с₁ (t_к2 - t_к1) + w₂(I_вп2 - с_вt_к2) + Q_2конц] (6),

где

1,03 - коэффициент, учитывающий 3% потерь тепла в окружающую среду;

с_н, с₁, с₂ - теплоемкости растворов соответственно исходного (начальной концентрации), в первом и во втором корпусе, Дж / кг * К [3];

t_н - температура кипения исходного раствора при давлении в 1-м корпусе, ?С;

Q_конц, Q_2конц - теплота концентрирования по корпусам, кВт;

При решении уравнений (5) - (6) можно принять I_вп1 ? I_г2; I_вп2 ? I_бк.

Величины Q_конц, Q_2конц настолько малы, что при расчетах ими можно пренебречь.

Тогда получим:

Q₁ = D (2730 - 502,4) = 1,03 [0,83*3,7 (126,4 - 110,5)+ w₁(2693 - 4,19*126,4)

Q₂ = w₁(2693 - 456,3) = 1,03 [(0,83 - w₁) 3,5 (112,4 - 126,4)+

+ w₂(2663 - 4,19*112,4)]

W = w₁ + w₂ = 0,67

Решение системы уравнений дает следующие результаты:

D = 0,36 кг/с; w₁ = 0,34 кг/с; w₁ = 0,33 кг/с;

Q₁ = 807,8 кВт; Q₂ = 760,5 кВт.

Полученные величины сводим в таблицу I.

Наибольшее отклонение вычисленных нагрузок по испаряемой воде в каждом корпусе предварительно принятых (w₁ = 0,32 кг/с, w₂ = 0,35 кг/с) не превышает 3%, поэтому не будем пересчитывать концентрации и температуры кипения растворов по корпусам.

Таблица I. Параметры растворов и паров по корпусам

Параметр	Корпус
	1	2
Производительность по испаряемой воде w, кг/с	0,34	0,33
Концентрация растворов х, %	19	62
Давление греющих паров Рг, МПа	0,3	0,14
Температура греющих паров tг, ?С	132,9	118,0
Температурные потери УД, ?С	19,6	25,9
Температура кипения раствора tк, ?С	126,4	112,4
Полезная разность температур Дtп, ?С	6,5	5,6

2.1.5 Выбор конструкционного материала

Выбираем конструкционный материал, стойкий в среде кипящего раствора квасного сусла в интервале от 12% до 69% [4]. В этих условиях химически стойкой является сталь марки Х17. Скорость коррозии ее менее 0,1 мм/год, коэффициент теплопроводности л_ст = 25,1 Вт/(м*К).

Физические свойства кипящего раствора квасного сусла

Параметр	Корпус	Литера-тура
	1	2
Теплопроводность раствора л Вт/(м*К)	0,61	0,62	[3]
Плотность раствора p, кг/м3	1062	1104	[3]
Теплоемкость раствора с, Дж/(кг*К)	3771	3561	[3]
Вязкость раствора м, Па*с	0,1 * 10-3	0,29*10-3	[3]
Теплота парообразования rв, Дж/кг	594*103	618*103	[1]
Плотность пара р, кг/м3	3,75	2,0	[1]

б₂ = 18,76 (б₁ /Дt₁₎^0,6 = 18,76 (1544)^0,6 = 1100 Вт/(м²*К)

Проверим правильность первого приближения по равенству удельных тепловых нагрузок:

q^/ = a₁ * Дt₁ = 1544 * 2 = 3088 Вт/м²

q^// = a₂* Дt₂ = 1100 * 2,47 = 2717 Вт/м²

Как видим q^/ ? q^//

Для второго приближения примем Дt₁ = 3,0?С. Пренебрегая применением физических свойств конденсата при изменении температуры на 1,0?С, рассчитываем б₁ по соотношению:

б₁ = 3536 ⁴v 2/3 = 1377 Вт/(м²*К)

Тогда получим:

Дt_ст = 1377 * 3 * 2,87 * 10^-4 = 1,18 ?С

Дt₂ = 6,5- 3 - 1,18 = 2,32 ?С

б₂ = 18,76 * (1337 * 3)^0,6 = 1429 Вт/(м² * К)

q^/ = 1377 * 3 = 4011 Вт/м²

q^// = 1429 * 2,32 = 3315,3 Вт/м²

Очевидно, что q^/ ? q^//

Для расчета в третьем приближении строим графическую зависимость удельной тепловой нагрузки q от разности температур между паром и стенкой в первом корпусе и определяем Дt.

Проверим общую полезную разность температур установки:

УДt_п = Дt_п1 + Дt_п2 = 4,5 + 7,6 = 12,1?C.

Найденные значения мало отличаются от ориентировочно определенной ранее поверхности F_ор. Поэтому в следующих приближениях нет необходимости вносить коррективы на изменение конструктивных размеров аппаратов (высоты, диаметра и числа труб). Сравнение распределенных из условия равенства поверхностей теплопередачи и предварительно рассчитанных значений полезных разностей температур представлено ниже:

	1-й корпус	2-й корпус
Распределенные в I-м приближении значения Дtп, ?C	4,5	7,6
Предварительно рассчитанные значения Дtп,?C	6,5	5,6

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

2.1.6 Уточненный расчет поверхности теплопередачи

Второе приближение.

В связи с тем, что существенное изменение давлений по сравнению с рассчитанным в первом приближении происходит только 1-м корпусе, где суммарные температурные потери незначительны, во втором приближении принимаем такие же значения Д^/, Д^//, Д^/// для каждого корпуса, как в первом приближении. Полученные после перераспределения температур (давлений) параметры растворов и паров по корпусам сводим в таблицу III:

Параметры растворов и паров по корпусам после перераспределения температур

Параметры	Корпус
	1	2
Производительность по испаряемой воде w, кг/с	0,28	0,31
Концентрации растворов х, %	19	62
Температура греющего пара в 1-й корпус tг1, ?C	132,9	-
Полезная разность температур Дtп,?C	6,5	5,6
Температура кипения раствора tк =tг1 -Дtп, ?C	126,4	127,3
Температура вторичного пара tвп = tк - (Д/ + Д//), ?C	107,8	102,4
Давление вторичного пара Рвп, МПа	0,13	0,12
Температура греющего пара tг = tвп - Д ///,?C	106,8	101,4
Удельная энтальпия вторичного пара rвп, при давлениях Рвп	2690	2686

Рассчитаем тепловые нагрузки (в кВт)

Q₁ = 1,03 * [0,83 * 3,7 * (126,4 - 119,0) + 0,34 * (2690 - 4,19 * 126,4)] = 779,9 кВт

Q₂ = 1,03 * [0,49 * 3,5 * (127,3 - 126,4) + 0,33 * (2686 - 4,19 * 127,3)] = 733,3 кВт

Расчет коэффициентов теплопередачи, выполненный описанным выше методом, приводит к следующим результатам:

К₁ = 843 Вт/(м²*К)

К₂ = 658 Вт/(м²*К)

Проверка суммарной полезной разности температур:

УДt_п = 5,5 + 6,6 = 12,1 ?C

Сравнение полезных разностей температур, полученных во втором и первом приближениях сводим в таблицу IV:

Таблица IV

	Корпус
	1	2
Значения Дtп во втором приближении, град	5,5	6,6
Значения Дtп во первом приближении, град	4,5	7,6

Различия между полезными разностями температур по корпусам в 1-м и 2-м приближениях не превышают 5%.

По ГОСТ 11987 - 81 [6] выбираем выпарной аппарат со следующими характеристиками:

Номинальная поверхность теплообмена F_н ………………160 м²

Диаметр труб d………………………………………………38 х 2

Высота труб Н……………………………………………….5000 мм

Диаметр греющей камеры d_к……………………………….1000 мм

Диаметр сепаратора d_c………………………………………2400 мм

Общая высота аппарата Н_а………………………………….12000 мм

Масса аппарата М_а…………………………………………..10000 кг

2.2 Определение толщины тепловой изоляции

Толщину тепловой изоляции д_и находят из равенства удельных тепловых потоков через слой изоляции от поверхности изоляции в окружающую среду

б_в * (t_ст2 - t_в) = (д_и / л_и) * (t_ст1 - t_ст2), (11)

где

б_в - коэффициент теплоотдачи от внешней поверхности изоляционного материала в окружающую среду, Вт/ (м² * К) [7]:

б_в = 9,3 + 0,058 * t_ст2

t_ст2 - температура изоляции со стороны окружающей среды (воздуха); для аппаратов, работающих в закрытом помещении, выбирается в интервале 35 - 45 ?С, а для аппаратов, работающих на открытом воздухе в зимнее время - в интервале 0 - 10 ?С;

t_ст1 - температура изоляции со стороны аппарата; ввиду незначительного термического сопротивления стенки аппарата по сравнению с термическим сопротивлением слоя изоляции t_ст1 принимают равной температуре греющего пара t_г1;

t_в - температура окружающей среды (воздуха), град;

л_и - коэффициент теплопроводности изоляционного материала, Вт/(м *К)

Выберем в качестве материала для тепловой изоляции совелит (85% магнезии + 15% асбеста) [8], имеющий коэффициент теплопроводности л_и = 0,09 Вт/(м * К)

б_в = 9,3 + 0,058 * 40 = 11,6 Вт/(м * К)

Принимаем толщину тепловой изоляции 0,03 м и для других корпусов.

2.3 Расчет барометрического конденсатора

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

Необходимо рассчитать расход воды, основные размеры (диаметр, высоту) барометрического конденсатора и барометрической трубы, производительность вакуум-насоса.

2.4 Расчет производительности вакуум - насоса

Производительность вакуум-насоса

G_возд = 2,5*10^-5(w₂ + G_в) + 0,01 w₂ (15),

где

2,5*10^-5 - количество газа, выделяющегося из 1 кг воды;

0,01 - количество газа, подсасываемого в конденсатор через неплотности на 1 кг паров.

Тогда:

G_возд = 2,5*10^-5 (0,33+64,63) + 0,01*0,31 = 474 кг/с

Температуру воздуха рассчитывают по уравнению

t_возд = t_н + 4 + 0,1 (t_к - t_н) = 20+4+0,1 (90-20)=31 ?С

Давление воздуха равно:

Р_возд = Р_бк - Р_п

где

Р_п - давление сухого насыщенного пара (Па) при t_возд = 31 ?С

Зная объемную производительность V_возд и остаточное давление Р_бк, по ГОСТ 1867-57 [8] подбираем вакуум-насос типа ВВН-25 мощностью на валу N=48 кВт

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

1. К.Ф. Павлов, П.Г. Романков, А.А. Носков «Примеры и задачи по курсу процессов и аппаратов химической технологии», Ленинград, 1976 г.

2. Ю.И. Дытнерский «Основные процессы и аппараты химической технологии», Москва, 2007 г.

3. Гинзбург М.П., Громов С.И. «Теплофизические характеристики пищевых производств», Москва, 1985 г.

4. Воробьева Т.Я. «Коррозионная стойкость материалов в агрессивных средах химических производств», Москва, 1994 г.

5 Касаткин А.Г. «Основные процессы и аппараты химической технологии», Москва, «Химия», 1973 г.

7. ГОСТ 11987-81, «Аппараты выпарные»

8. ГОСТ 1867-57, «Вакуум-насосы»

Размещено на Allbest.ru