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

ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО ОБРАЗОВАНИЮ РФ

Государственное образовательное учреждение

высшего профессионального образования

СИБИРСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНОЛОГИЧЕСКИЙ УНИВЕРСИТЕТ

Кафедра " Процессов и аппаратов химических производств "

ПОЯСНИТЕЛЬНАЯ ЗАПИСКА

РАСЧЕТ НАСАДОЧНОЙ РЕКТИФИКАЦИОННОЙ КОЛОННЫ

Проверил:

Лис Е.В.

Выполнил:

студент группы 64-7

Мымликова Т.Н.

Авап

Содержание

Введение

1. Материальный баланс

2. Расчет флегмого числа

3. Определение скорости пара и диаметра колонны

4. Определение числа тарелок и высоты колонны

5. Тепловой расчет

6. Гидравлический расчет

7. Расчет тепловой изоляции колонны

8. Расчет вспомогательного оборудования

Заключение

Библиографический список

Приложение

Реферат

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

Курсовой проект содержит расчетно-пояснительную записку из 25 страниц текста, таблицы, рисунка, 2 приложений, 5 литературных источников и графическую часть из 2 листов формата А1 технологической схемы, чертежа общего вида.

Введение

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

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

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

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

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

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

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

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

I. Технологические расчёты

1.Материальный баланс

Концентрация исходной смеси, дистиллята и кубового остатка в массовых долях.

где х_А-мольная доля низкокипящего компонента;

М_А-молекулярная масса низкокипящего компонента, кг/кмоль;

М_В-молекулярная масса высококипящего компонента, кг/кмоль.

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

Материальный баланс для всей колонны:

G_f=G_d+G_w

G_f*x_f=G_d*x_d+G_w*x_w

G_d= (G_f*x_f- G_f* x_w)/ x_d - x_w

где G_f- массовый расход исходной смеси, кгс/с;

G_d- массовый расход дистиллята кгс/с;

G_w- массовый расход кубовой жидкости, кгс/с;

x_f, x_d, x_w- содержания низкокипящего компонента в массовых долях в исходной смеси, дистилляте, кубовой жидкости.

G_f= G_d (x_d - x_w)/( x_f - x_w)=0,7(0,95-0,05)/(0,25-0,05)=3,15 кг/с;

G_w = G_f + G_d=3,15-0,7=3,08 кг/с.

Относительный мольный расход питания:

2.24

2.Расчет флегмового числа

Минимальное флегмовое число можно рассчитать по формуле:

где y_f^*-мольная доля низкокипящего компонента в паре, равновесном с исходной смесью.

Оптимальное флегмовое число определим по формуле n_t(R+1), где n_t-число ступеней изменения концентрации.

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

· Задаемся рядом значений коэффициента избытка флегмы ? в пределах от 1,1 до 5,5;определяем рабочее флегмовое число R=?R_min и величину отрезка

;

Откладываем отрезок В на оси ординат и проводим линии АВ1, АВ2… . На оси абцисс откладываем содержание низкокипящего компонента в исходной смеси xf и в кубовом остатке xw. Из xf восстанавливаем перпендикуляр до пересечения с линиями АВ1, АВ2… . получаем точки D1, D2… . Линии AD1, AD2… являются линиями рабочих концентраций верхней части колонны. Соеденив точки D1, D2… с точкой С получаем линии рабочих концентраций нижней части колонны D1C, D2C… .

· Между равновесной и рабочими линиями в пределах концентраций x_w и x_d строим ступени, каждая из которых соответствует теоретической тарелке;

· При каждом значении ? определяем число теоретических тарелок n_t и величину произведения n_t(R+1). Результаты расчета сводим в таблицу 2.

· По данным таблицы 2 строим график зависимости n_t(R+1)=f(R) и находим минимальное значение величины n_t(R+1). Ему соответствует флегмовое число R=12,804.

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

Таблица-2

?	R	В	nt	nt(R+1)
1,10	11,737	0,075	33	420,321
1,20	12,804	0,069	28	386,512
2,00	21,34	0,043	25	558,5
2,80	29,876	0,031	23	710,148
3,60	38,412	0,024	21	827,652

3. Расчет насадочной колонны

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

Характеристика насадки:

размер кольца………………35х35х4мм

удельная поверхность,?.......140м²/м²

свободный объем, V_св……...0,78м³/м³

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

lg

где ?_п,?_ж- плотность пара и жидкости;

µ_ж- динамический коэффициент вязкости жидкости, Па*с;

µ- динамический коэффициент вязкости воды при 20⁰С, Па*с;

L/G- отношение массовых расходов жидкости и пара, кг/кг;

g- ускорение свободного падения.

Средние массовые расходы по жидкости:

А) В верхней (укрепляющей) части колонны

L_жв=G_DR=8,6 кг/с

Б) Нижней (исчерпывающей) части колонны

L_жн=G_DRG_F11,5 кг/с

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

?^|_пв=2,24 м/с

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

?^|_пн=2,05 м/с

Примем рабочую скорость пара на 25% ниже предельной, тогда рабочая скорость:

в верхней части колонны ?^|_пв=2,24-0,75=1,49 м/с

в нижней части колонны ?^|_пн=2,05-0,75=1,3 м/с

Средние массовые потоки пара

а) в верхней части колонны

G_пв=G_D(R+1)=10,22 кг/с

б) в нижней части колонны

G_пн=G_D(R+1)= 11,03 кг/с

Диаметр верхней части колонны:

D_B= =2,7 м

Диаметр нижней части колонны:

D_н==2,5 м

Принимаем стандартный диаметр обечайки D_к=2800мм одинаковым для обеих частей колонны, при этом действительная скорость пара в колонне:

а) в верхней части колонны

?^|_пв=1,49 м/с

б) в нижней части колонны

?^|_пн=1,1 м/с

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

где ?- удельная поверхность насадки, м²/м³;

b- поправочный коэффициент, м³/мч; для процесса ректификации b=0,065

Фактическая плотность орошения для верхней части колонны

Отношение

Коэффициент смачиваемости насадки ?_В=1

Фактическая плотность орошения для нижней части колонны

Отношение

Коэффициент смачиваемости насадки ?_В=1

Равномерное смачивание насадки соблюдается при соотношении диаметра колонны и диаметра элемента насадки, если . Для нашего случая соблюдается и это условие:

дистиллят фаза холодильник

4. Высоты колонны

Высота слоя насадки в обычных насадочных колоннах, работающих при оптимальном режиме, может быть определена по уравнению:

H=n_оуh_оу

где n_оу- общее число единиц переноса по паровой фазе;

h_оу-общая высота единиц переноса, м.

n_oy=SM_xM_y

где S- площадь,ограниченная кривой, ординатами y_w и y_D и осью абцисс.

М_х, М_у- масштаб по оси х, у.

Для этого по диаграмме х - у снимают значения рабочих концентраций х, у и соответствующих им равновесных концентраций у^*.

Общее число единиц переноса верхней и нижней частях колонны

n_oy.в=10,2 n_оу.н=4,03

Определение высоты единиц переноса

Определяем общую высоту единиц переноса

где h_x, h_y- высота единиц переноса соответственно в жидкой и паровой фазах;

m- средний коэффициент распределения в условиях равновесия для соответствующих частей колонны.

Отношение нагрузок по пару и жидкости

а) для верхней части колонны

б) для нижней части колонны

где F-относительный (на один кмоль дистиллята) мольный расход питания

Относительный мольный расход питания:

тогда

Высота единиц переноса в жидкой фазе

h_x=0,258*Ф*с*Pr_ж^0,5*Z^0,15

где Ф- коэффициент, зависящий от массовой плотности орошения L/D, кг/м²ч;

с- коэффициент, зависящий от отношения рабочей скорости пара к предельной ?/?_п;

Pr_ж- критерий Прандтля для жидкостей при средней температуре в верхней и нижней частях колонны.

Расчет коэффициентов с, Ф для определения высоты единиц переноса.

Определение коэффициента с

а) для верхней части колонны

c=0,58

б) для нижней части колонны

с=0,61

Определение коэффициента Ф.

Рассчитываем часовую массовую плотность орошения:

а) для верхней части колонны

Коэффициент Ф_В=0,1

б) для нижней части колонны

Коэффициент Ф_В=0,12

Высота единиц переноса:

а) для верхней части колонны

h_XB=0,064м

б) для нижней части колонны

h_XН=0,095м

Высота единиц переноса в паровой фазе h_y

где L_s- массовая плотность орошения, кг/м²с;

D_k- диаметр колонны, м;

f₁, f₂, f₃ -коэффициенты, зависящие от физических свойств жидкости.

а) для верхней части колонны

f₁=µ_x^0,16(µ_x-в мПа с)=0,62

f₂=(1000/?_x)^1,25=1,28

f₃=(72,8*10^-3/?_ч)^0,8=2,82

б) для нижней части колонны

f₁=µ_x^0,16(µ_x-в мПа с)=0,58

f₂=(1000/?_x)^1,25=1,35

f₃=(72,8*10^-3/?_ч)^0,8=2,85

Критерий Прандля в паровой фазе для верхней части колонны

а) для верхней части колонны

Pr_пв=0,15

б) для нижней части колонны

Pr_пн=0,17

Массовая плотность орошения L_s

а) для верхней части колонны

L_sв=2,08 кг/м²с

б) для нижней части колонны

L_sн=2,8 кг/м²с

Высота единиц переноса в паравой фазе

а) для верхней части колонны

б) для нижней части колонны

По уравнению находим общую высоту единиц переноса

а) для верхней части колонны

б) для нижней части колонны

Высота насадки в верхней и нижней частях колонны, соответственно равны:

H_B=h_oyB*n_oyB=1,18*10,2=12,04м

H_Н= h_oyН*n_oyН=0,88*4,03=3,55м

Общая высота насадки в колонне равна:

Н=12,04+3,55=15,59м

С учетом того, что высота слоя в одной секции Z=3,25 м, общее число секций в колонне составляет n=4.

Общая высота ректификационной колонны:



где h_p- высота промежутков между секциями насадки

H_k=3,25*4+(4-1)0,5+1+2=17,5 м

5. Тепловой расчёт ректификационной колонны

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

где Q_D - расход теплоты, отнимаемой охлаждающей водой от конденсирующихся в дефлегматоре паров;

C_F, C_D, C_W - теплоёмкости исходной смеси, дистиллята кубовой жидкости, соответственно, определим по формуле

C_в., C_э.с, - теплоёмкости воды и уксусной кислоты при t_f=84⁰C, равны 4190 Дж/(кг*К) и 2284 Дж/(кг*К) соответственно (с. 562 [5]);

C_в., C_э.с, - теплоёмкости воды и уксусной кислоты при t_w=78.9⁰C, равны 4190 Дж/(кг*К) и 2263 Дж/(кг*К) соответственно (с. 562 [5]);

C_в., C_э.с, - теплоёмкости воды и уксусной кислоты при t_d=88.5⁰C, равны 4190 Дж/(кг*К) и 2288 Дж/(кг*К) соответственно (с. 562 [5]);

Расход теплоты, отнимаемой охлаждающей водой от конденсирующихся в дефлегматоре паров, Q_D, определим по формуле

где r_d - удельная теплота конденсации дистиллята, определим по формуле

где r_в и r_э.с. - удельные теплоты конденсации воды и уксусной кислоты, соответственно, при t_d=84⁰С, r_в=2308 кДж/кг, r_э.с.=388 кДж/кг (с. 541 [5]);

Тепловые потери колонны в окружающую среду, , определим по формуле

где - температура наружной стенки колонны, принимаем =40⁰С;

- температура воздуха в помещении, принимаем =20⁰С;

F_H - наружная поверхность изоляции колонны, определим по формуле

- суммарный коэффициент теплоотдачи конвекцией и излучением, определим по формуле

Подставляя значения в формулу, получим

Подставляя значения в формулу, получим

Расход греющего пара в кубе колонны определим по формуле

где r_гр.п. и - удельные теплоты парообразования при давлении P_гр.п.0,25 МПа, r_гр.п.=2227 кДж/кг (с. 549 [5]);

x - степень сухости, принимаем х=0,95

Расход тепла в подогревателе исходной смеси

где C_f. - теплоёмкость смеси определённая при средней температуре

Расход греющего пара в подогревателе равен



Общий расход греющего пара равен

Расход воды в дефлегматоре при нагревании её на 20⁰С определим по формуле

Расход воды в холодильнике дистиллята при нагревании её на 20⁰С определим по формуле

где C_d. - теплоёмкость смеси определённая при равна

Расход воды в холодильнике кубового остатка при нагревании её на 20⁰С:

где C_d. - теплоёмкость смеси определённая при равна

Общий расход воды в ректификационной установке

6. Гидравлический расчет

Определяем число Рейнольдса для верхней части колонны

Для турбулентного режима коэффициент сопротивления сухой насадки в виде беспорядочно рассыпанных колец Рашига

Потеря давления на 1м насадки

Па/м

Отношение массовых расходов жидкости и паров верхней части колонны

Сопротивление орошаемой насадки в верхней части колонны высотой 1 м

Па/м

Сопротивление слоя насадки верхней части колонны

Па

Определяем число Рейнольдса для нижней части колонны

Потеря давления на 1м насадки

Па/м

Отношение массовых расходов жидкости и пара в нижней части колонны

Сопротивление орошаемой насадки в верхней части колонны высотой 1 м

Па/м

Сопротивление слоя насадки верхней части колонны

Па

Общее сопротивление колонны

Па

7. Расчет тепловой изоляции колонны

В качестве изоляции берем асбест (?_из=0,151 Вт/(м*К)). Исходя из упрощенного соотношения (для плоской стенки) имеем

где ?_из-толщина изоляции, м; t_из.в-температура внутренней поверхности изоляции, принимаем её ориентировочно на 10-20⁰С ниже средней температуры в колонне - t_из.в=110⁰С.

Определяем толщину изоляции

мм

⁰С

расхождение: 110,34-110=0,34⁰С<1,0⁰C

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

8. Расчет вспомогательного оборудования

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

Кипятильник (куб-испаритель)

Температурные условия процесса.

Кубовой остаток кипит при 106,8 ⁰С. Согласно заданию температура конденсации греющего пара равна 120⁰С ( = 1 атм). Следовательно, средняя разность температур, определяется по формуле

Подставив значения в формулу (78) получим

⁰С.

Определяем максимальную величину площади теплообмена , м². По [2, с 172] для данного случая теплообмена (от конденсации водяного пара к кипящей жидкости) принимаем значение минимального коэффициента теплопередачи . Тогда максимальная поверхность теплообмена определяется по формуле

м²

Определение коэффициента теплопередачи.

Принимаем среднее значение тепловой проводимости загрязнении стенок : со стороны конденсирующего водяного пара Вт/м² К, со стороны кубового остатка Вт/м² К.

Теплопроводность стали ?=46,5 Вт/м К.

Таким образом, термическое сопротивление стенки и ее загрязнений

, , определяется по формуле

Подставим значения в формулу (80) получим

Коэффициент теплоотдачи со стороны конденсирующегося пара , Вт/м²•К , определяем по формуле

Где - коэффициент теплопроводности конденсата Вт/м•К;

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

- плотность конденсата, кг/м³;

- удельная теплота конденсации греющего пара, Дж/кг;

- динамический коэффициент вязкости конденсата, Па•с

Значение коэффициентов теплопроводности, динамической вязкости, плотности определяется [2, с. 537].

Принимаем = 0,54 Вт/м•К, = 0,54•10⁶ Па•с, = 1121 кг/м³.

Вт/м²•К

Коэффициент теплоотдачи для кипящей муравьиной кислоты , Вт/м²•К , определяем по формуле



Где - коэффициент теплопроводности метанола Вт/м•К;

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

- динамический коэффициент вязкости кипящей муравьиной кислоты,Па•с

- поверхностное натяжение метанола.

Коэффициент ,определяется по формуле

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

формуле

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

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

Т - температура кипения муравьиной кислоты, К;

Р,Р₀ - давление в кипятильнике при рабочих и нормальных условиях

Плотность паров муравьиной кислоты при нормальных условиях , кг/м³, определяется по формуле



где М - мольная масса уксусной кислоты,

кг/м³

кг/м³

Значение коэффициентов теплопроводности, динамической вязкости, плотности определяется из таблиц [2].

Принимаем = 0,20 Вт/м•К, = 0,22•10^-3 Па•с, = 717,3 кг/м³,

= 15,99•10³ Н/м.

Коэффициент теплопередачи определяется по формуле :

Удельная тепловая нагрузка определяется по формуле



Откуда

Это уравнение решаем графически, задаваясь значениями q (рис), где y - левая часть уравнения.

При q = 2000 Вт/м² y = 4,42

При q = 1000 Вт/м² y = -1,42

При y = 0 находим q = 1250 Вт/м²

Коэффициент теплопередачи , Вт/м²•К, определяется по формуле

Вт/м²•К

Коэффициент поверхности теплообмена , м², определяется по формуле

м²

С запасом 15-20% принимаем по каталогу [2,с.51] 4 - ходовой теплообменник.

Характеристика теплообменника:

поверхность теплообмена - 397 м²

диаметр кожуха - 1200 мм

диаметр труб - 20х2 мм

длина труб - 4м

количество труб - 1580

Дефлегматор

В дефлегматоре конденсируется вода с небольшим количеством уксусной кислоты. Температура конденсации пара дистиллятора = 101,75⁰С. Температуру воды на входе в теплообменник примем 18 ⁰С, на выходе 38 ⁰С.

Составляем температурную схему процесса и определяем движущую силу процесса теплопередачи:

101,75>101,75

38>18

⁰С

По таблице [2.с,172] коэффициент теплопередачи от конденсирующегося пара органических веществ к воде находиться в пределах 240-870 Вт/м²•К. Принимаем наименьший коэффициент теплопередачи К=340 Вт/м²•К

Количества тепла, отнимаемого охлаждающей водой от конденсирующегося в дефлегматоре пара: = 410905 Вт.

Поверхность дефлегматора , м², находим из основного уравнения теплопередачи .

м²

С запасом принимаем по каталогу [2,с.51] 2 - ходовой теплообменник с поверхностью.

Характеристика теплообменника:

поверхность теплообмена - 38 м²

диаметр кожуха - 600 мм

диаметр труб - 25х2 мм

длина труб - 2 м

количество труб - 240

Холодильник для дистиллята

В холодильнике происходит охлаждение дистиллята от температуры конденсации до 30 ⁰С. Составляем температурную схему процесса и определяем движущую силу процесса теплопередачи:

101,75>30

38>18

⁰С

Количества тепла, отнимаемого охлаждающей водой от дистиллята в дефлегматоре, определяется по формуле

Вт

Поверхность дефлегматора , м², находим из основного уравнения теплопередачи.

Принимаем коэффициент теплопередачи К = 300 Вт/м²•К

м²

С запасом принимаем по каталогу [2,с.51] 1 - ходовой теплообменник с поверхностью.

Характеристика теплообменника:

поверхность теплообмена - 25 м²

диаметр кожуха - 325 мм

диаметр труб - 20х2 мм

длина труб - 4 м

количество труб - 400

Холодильник для кубового остатка.

В холодильнике кубового остатка происходит охлаждение кубовой жидкости от температуры кипения до 30 ⁰С. Составляем температурную схему процесса и определяем движущую силу процесса теплопередачи:

100,8>30

38>18

⁰С

Количества тепла, отнимаемого охлаждающей водой от кубовой жидкости, определяется по формуле

Вт

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

Принимаем коэффициент теплопередачи К = 250 Вт/м²•К

м²

С запасом принимаем по каталогу [2,с.51] 4 - ходовой теплообменник с поверхностью.

Характеристика теплообменника:

поверхность теплообмена - 95 м²

диаметр кожуха -800 мм

диаметр труб - 25х2 мм

длина труб - 3 м

количество труб - 404

Подогреватель

Служит для подогрева исходной смеси 20 ⁰С до температуры = 88 ⁰C. Исходная смесь подогревается водяным насыщенным паром с температурой 120 ⁰С. Составляем температурную схему процесса и определяем движущую силу процесса теплопередачи:

160>160

20>111,75

⁰С

Принимаем коэффициент теплопередачи К = 120 Вт/м²•К

м²

С запасом принимаем по каталогу [2,с.51] 2 - ходовой теплообменник с поверхностью.

Характеристика теплообменника:

поверхность теплообмена - 113 м²

диаметр кожуха -600 мм

диаметр труб - 25х2 мм

длина труб - 6 м

количество труб - 240

Заключение

В данной курсовой работе было определено количество дистиллята и кубового остатка; диаметр и высота колонны Н=17,5 м., D=3 м.; число тарелок n=4; параметры вспомогательного оборудования.

Библиографический списрк

1. Ченцова,Л.И. Процессы и аппараты химической технологии:учебное пособие к самостоятельной работе/ Л.И. Ченцова, М.К. Шайхудинова, В.М. Ушанова.- Красноярск: СибГТУ,2006.-267с.2. Павлов, К.Ф., Романков П.Г., Носков А.А. Примеры и задачи по курсу процессов и аппаратов х

2. мической технологии: Учебное пособие для вузов. Перепечатка с изд. 1987г.- М.:ООО «РусМедиаКонсалт», 2006.-576с.

3. Шайхудинова М.К., Ченцова Л.И., Борисова Т.В. Процессы и аппараты химической технологии. Расчет выпарной установки: учебное пособие к выполнению курсового проекта.-Красноярск: СибГТУ, 2007.- 80с.

4. Левин Б.Д., Ченцова Л.И., Шайхутдинова М.Н., Ушанова В.М. процессы и аппараты химических и биологических технолгий. Учеб. пособие для студентов химических специальностей вузов / Под общ. ред. д-ра. хим. Наук С.М. Репяха. - Красноярск: Сибирский государственный технологический университет, 2006. - 430с.

5. Основные процессы и аппараты химической технологии: Пособие по проектированию/ Г.С. Борисов, В.П. Брыков, Ю.И. Дытнерский и др. Под ред. Ю.И. Дытнерского, 2-е изд., перераб. и дополн. М.: Химия, 2008. - 496 с.

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