Проектирование ректификационной установки для разделения бинарной бензол - уксусная кислота, производительностью по дистилляту 200 кг/ч

Температурный коэффициент  b определим по формуле:

          

где и принимают при температуре 20 °С

Плотность жидкой смеси при температуре 20 по Цельсию:

Верхняя часть колонны

         

 

 

 

Нижняя часть колонны:

         

Отсюда:

3.1.4 Скорость пара  и диаметр колонны

Выбор рабочей скорости паров  обусловлен многими факторами и  обычно осуществляется путем технико-экономического расчета для каждого конкретного  процесса. Для ректификационных колонн, работающих в пленочном режиме при атмосферном давление, рабочую скорость можно принять на 20-30 % ниже скорости захлёбывания. Предельную фиктивную скорость пара , при которой происходит захлёбывание насадочных колонн, определяют по уравнению:

где   - предельная скорость пара в критических точках, м/с; а – удельная поверхность насадки, м²/м³; ε – свободный объём насадки, м³/м³; μ_х – динамический коэффициент вязкости жидкости, мПа∙с; и - массовые расходы жидкой и паровой (газовой) фаз, кг/с; и - плотность пара и жидкости соответственно, кг/м³; А и В -  коэффициенты.

Выбираем кольца Рашига 50 мм:

Верх:

Низ:

Отсюда:

Рабочая скорость пара (газа) рассчитывается по соотношению:

Для определения  диаметра колонны пользуются формулой:

где - объёмный расход пара (газа) при рабочих условиях в колонне, м³/с.

Выбираем стандартный аппарат с  диаметром 0.4 м для обеих частей колонны и уточняем рабочую скорость пара и

Уточнённую рабочую скорость проверяют по графической зависимости  Эдулджи,  где комплексы Y и X  имеют следующий вид:

 

где – критерий Фруда рассчитывается по номинальному размеру насадки d (м); - критерий Рейнольдса (условный), также рассчитывается по номинальному размеру d; и - соответственно плотность воды и орошаемой жидкостью, кг/м³.

где и - объёмные расходы жидкости и пара (газа), м³/с; U – плотность орошения.

Расчет плотности орошения производится по формуле:

  

где U – плотность орошения, м³/(м²∙с); V_x – объёмный расход жидкости, м³/с; S – площадь поперечного сечения колонны, м².

Критерий Рейнольдса:

Критерий  Фруда:

Рассчитаем комплексыY и X:

Точки (X_в,Y_в) и (X_н, Y_н) лежат ниже линии захлебывания, которая соответствует неустойчивому режиму работы. Это означает, что работа колонны устойчива и выбор диаметра колонны правильный.

 

Определение гидравлического  сопротивления 1 м насадки:

где скорость пара (газа) в свободном сечении насадки (действительная);

 коэффициент сопротивления

  

  

                                  

   

В результате получаем:

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

Активная поверхность  насадки.

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

где U – плотность орошения, м³/(м²с)

a – удельная поверхность насадки, м²/м³

p и q – постоянные, зависящие от типа и размера насадки (кольца Рашига 50 мм: 1.019 10^-5 и 0.0086)

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

 

3.1.5  Расчет  высоты насадки и колонны.

Высоту насадки Н рассчитываем по модифицированному уравнению массопередачи.

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

                                      

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

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

 

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

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

 

где - высота единицы переноса по паровой фазе, м; ψ- коэффициент, определяемый по рис; - диффузионный критерий Прандтля; - массовая плотность орошения, ; - массовый расход жидкости, ; D – диаметр колонны, м; z – высота насадки одной секции (z не должна превышать – 3 м), м; , - динамический коэффициент вязкости жидкости, мПа∙с; ;

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

  

     

                  

    

    

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

                          

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

Выбранные коэффициенты из графика:

в:                              Ф =0,05        , с =1  , =200  ,

н:                              Ф =0,05   ,     с =1  , =215  ,

Рассчитаем высоту чисел переноса по жидкости:

                                 .

                                 .

Тогда высота чисел  переноса по пару равна:

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

Высота насадки для верха  и низа колонны:

 

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

        Общую высоту колонны определяют по формуле:

                    

где - высота насадки одной секции, м; - число секций ( );  - высота промежутков между секциями, м; и - высота сепарационного пространства над насадкой и расстояние между днищем колонны и насадкой, соответственно, м.

В соответствии с рекомендациями:

 

Диаметр  колонны, м	, м	, м
0,4 – 1,0	0,6	1,5

 

Величина  зависит от размеров распределительных тарелок (ТСН-3) и при проектировании м.

Тогда общая высота колонны  будет:

 

 

 

 

 

 

3.1.6  Расчет  гидравлического сопротивления  колонны

Гидравлическое сопротивление  насадки при пленочном режиме находят по приближенному соотношению:

где b – коэффициент;

      U – плотность орошения

   

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.1.7 Тепловой баланс  ректификационной колонны

Тепловой баланс ректификационной колонны выражается общим уравнением:

где Q_K – тепловая нагрузка куба; Q_D –количество теплоты, передаваемой от пара к воде; Q_пот – тепловые потери (5%); -теплоёмкости соответствующие дистилляту, кубовому остатку и исходной смеси; - температуры соответствующие дистилляту, кубовому остатку и исходной смеси (находим из диаграммы «Зависимость температуры от равновесных составов пара и жидкости»):

  , ,

Найдем удельную теплоту  конденсации паров дистиллята по аддитивной формуле:                                                                               Определим тепловую нагрузку дефлегматора:

 

Определим теплоёмкости смеси:

При температуре кипения кубовой  жидкости:

При температуре кипения дистиллята:

 

При температуре кипения  исходной смеси:

 

 

 

Тогда:  

 

 

 

 

   3.2 Приближённый  расчёт теплообменной аппаратуры 

В качестве греющего агента выберем  пар абсолютным давлением 3.5 кг·с/см², выбор обусловлен тем, что разница между температурой кипения кубовой жидкости и пара должна составлять 20-25 градусов

              

Примем, что вода имеет начальную  температуру 10°С.

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

 3.2.2 Куб-испаритель

Исходные данные:

Q=83430 Вт – тепловой баланс куба – испарителя

t_г.п.= 147°С – температура конденсации водяного пара

t_W=117.576 °С – температура кипения кубовой жидкости

Средняя движущая сила:

Для куба-испарителя коэффициент  теплопередачи К_ор=1000Вт/(м^2.К):

Определим поверхность  теплообмена по формуле:

 

Расход греющего пара:

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

 Диаметр кожуха 

Общее число труб (25×2 мм)

Число ходов 

Длина труб

Площадь поверхности теплообмена 

    

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.2.2  Дефлегматор

 

Исходные данные:

 

Q_D=74950 Вт

t_D=85.258°С – температура конденсации дистиллята

Примем, что воздух имеет  начальную температуру t_внач=20 °С и нагревается на 40 °С

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

Δt_б=t_D-t_внач=85.258-20=65.258 °С

Δt_м=t_D-t_вкон=85.258-40=45.258°С

 

 

 

Поверхность теплообмена (К_ор=50Вт/(м^2.К):

Теплоёмкость воздуха при его средней температуре:

Расход охлаждающей воды:

 

Диаметр кожуха

Общее число труб (20×2 мм)

Число ходов 

Длина труб

Площадь поверхности теплообмена 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

Исходные данные:

- расход дистиллята

t_D=85.258°С – начальная температура дистиллята

t_кон=20 °С – конечная температура дистиллята

t_внач=15 °С – начальная температура воды

t_вкон=30 °С – конечная температура воды

 

Определим среднюю движущую силу процесса:

Δt_Б=t_D-t_вкон=85.258-30=55.258 °С

Δt_М=t_кон-t_внач=20-15=5°С

   

Средняя температура и теплоёмкость воды:

Средняя температура дистиллята:

t_1ср=t_вср+ Δt_ср=22.5+20.918=43.418°С

Определим теплоемкость дистиллята при  t_1ср:

Тепловой баланс холодильника дистиллята:

Расход охлаждающей воды:

 

Поверхность теплообмена:

Коэффициент теплопередачи в холодильнике К_ор=500Вт/(м^2.К)

 

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

 Диаметр кожуховой трубы-

Диаметр теплообменной трубы-

_{Длина теплообменной  трубы - 6м}

_{Площадь теплообмена по одной трубе – 1.65м²}

_{Число труб -2}

_{Площадь теплообмена-3.3м²}

 

 

 

3.2.4 Подогреватель исходной смеси

 

 Исходные данные:

- расход исходной смеси

t_г.п.=147 °С – температура конденсации водяного пара

t_нач=25 °С – начальная температура исходной смеси

t_F=108.193°С – конечная температура исходной смеси

Определим среднюю движущую силу процесса:

 

Δt_м=t_гп-t_F=147-108.193=38.807°С

Δt_б=t_гп-t_нач=147-25=122°С

 

 

Средняя температура исходной смеси:

t_см=t_гп- Δt_ср=147-72.631=74.369°С

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

 

 

 

Тепловой баланс подогревателя  исходной смеси:

 

 

Ориентировочная площадь теплообмена  подогревателя:

 

Пусть К_ор=700Вт/(м^2.К), тогда

Расход греющего пара:

 

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

 

Диаметр кожуха

Общее число труб (25×2 мм)

Число ходов 

Длина труб

Площадь поверхности теплообмена 

 

 

 

 

 

 

 

 3.3 Подробный расчёт холодильника кубовой остатка

 

Исходные данные:

- расход кубовой жидкости

t_w=117.576°С – начальная температура кубовой жидкости

t_кон=20 °С – конечная температура кубовой жидкости

t_внач=15 °С – начальная температура воды