Минимальное флегмовое число уравнение андервуда

Основные теоретические положения. Минимальное флегмовое число определяется по методу Андервуда

Минимальное флегмовое число определяется по методу Андервуда. При этом сначала находится промежуточная характеристическая величина w по уравнению

(10)

где a_i/v – относительная летучесть по высококипящему ключевому компоненту;

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

Уравнение (10) решается методом последовательных приближений. Вычислив w, находим минимальное флегмовое число из уравнения Андервуда

(11)

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

Примем давление на верху колонны на 0,01МПа ниже давления в зоне тарелки питания. Температуру Т_верха колонны определяем итерациями по уравнению:

(12)

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

Давление низа колонны выше давления в зоне тарелки питания на суммарное гидравлическое сопротивление тарелок. Если с выбором тарелок пока не определились, можно принять его на 0,03 – 0,05МПа выше давления на верху колонны. Температуру Т_низа низа колонны определяем итерациями по уравнению (13) при фиксированном составе куба и фиксированном давлении:

(13)

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

Для определения тарелки питания используют уравнение Керкбрайда:

(14)

где N₁ , N₂ – число тарелок укрепляющей и исчерпывающей частей колонны, соответственно; i – индекс тяжелого ключевого компонента,

k – индекс легкого ключевого компонента.

Еще один метод расчета минимального орошения — Львова-Серафимова.

В этом случае рассчитывают так называемое «псевдоминимальное» флегмовое число для каждой из разделяемых пар по уравнению

(15)

Определяющей паре будет соответствовать максимальное значение произведения (z + R). Далее расчет ведется по компонентам определяющей пары

Находят числа тарелок в укрепляющей N + и исчерпывающей N — частях колонны. Пусть определяющей, например, является пара c-d:

(16)

(17)

Далее находятся значения величин z + и z — на уровне тарелки питания по уравнениям:

(18)

(19)

Индекс а относится к наиболее легколетучему компоненту.

Находится истинное значение минимального флегмового числа, соответствующего заданному разделению рассматриваемой смеси:

(20)

Очевидно, что метод Львова – Серафимова позволяет определять не только минимальное флегмовое число, но и число тарелок, соответствующее заданному разделению смеси в условиях минимального орошения. Кроме того, использование метода позволяет сразу определить номер тарелки питания сложной ректификационной колонны. Однако результаты, полученные данным методом для ректификации непрерывных смесей, как правило, существенно отклоняются от экспериментальных показателей работы колонн. Число тарелок, рассчитанное этим методом, несколько завышено.

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

Пример расчета минимального флегмового числа методом Андервуда, выполненный в приложении Microsoft Excel, приведен в Таблицах 4- 6.

Ход выполнения работы

1. Определить Т_верха

2. Определить Т_низа

4. Если необходимо, рассчитать долю отгона питания

5. Определить коэффициенты относительной летучести

7. Определить R_min методом Андервуда

8. Определить место ввода сырья

9. Определить R_min , число тарелоки номер тарелки питания методом Львова – Серафимова

Минимальное и действительное флегмовое число

ЛЕКЦИЯ № 14

Ректификация

Принцип ректификации.Как отмечалось, достаточно высокая степень разделения однородных жидких смесей на компоненты может быть достигнута путем ректификации. Сущность процессов, из которых складывается ректификация, и получаемые при этом результаты можно проследить с помощью t – x – y диаграммы (рис. 13).

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

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

Материальный баланс ректификационной колонны

Пусть, согласно схеме на рис. 16, в колонну поступает F кмоль исходной смеси, состав которой x_F мол. долей НК. Сверху из колонны удаляется G кмоль паров, образующих после конденсации флегму и дистиллят. Количество получаемого дистиллята Р кмоль, его состав х_Р мол. долей НК. На орошение колонны возвращается флегма в количестве Ф кмоль, причем ее состав равен составу дистиллята (х_ф = х_Р мол. долей). Снизу из колонны удаляется W кмоль остатка состава x_w мол. долей НК.

Тогда уравнение материального баланса колонны будет:

Поскольку G = Р + Ф, то

Соответственно материальный баланс по НК:

Уравнение рабочих линий. Для получения уравнений рабочих линий используем общее для всех массообменных процессов уравнение (11, а), выразив применительно к ректификации входящие в него концентрации в мольных долях:

. (11 а)

Уравнение (11 а) может быть представлено также в ином виде, если воспользоваться материальным балансом по распределяемому компоненту.

откуда следует, что

(11 в)

Делая подстановку в уравнение (11 а), получим

(11 г)

Укрепляющая часть колонны. Количество жидкости (флегмы), стекающей по этой части колонны

Здесь R = Ф/Р – флегмовое число, представляющее собой отношение количества флегмы к количеству дистиллята. Количество паров, поднимающихся по колонне

G=P + Ф = P + P R = P (R+1) (13)

Для верхнего конца укрепляющей части колонны состав паров у_G = у_P и, согласно принятому выше допущению, у_Р = х_Р. Следовательно, в данном случае у_к = х_Р.

В том же сечении колонны состав жидкости (флегмы), поступающей из дефлегматора, х_ф = х_Р, т. е. х_н = х_Р. Подставляя значения L, G, у_к и х_к в уравнение (11 а), получим

(14)

Зависимость (14) является уравнением рабочей линии укрепляющей части колонны. В этом уравнении – тангенс угла наклона рабочей линии к оси абсцисс, а – отрезок, отсекаемый рабочей линией на оси ординат диаграммы у—х (рис. 17).

Исчерпывающая часть колонны. Количество орошающей жидкости L’ в этой части колонны больше количества флегмы Ф, стекающей по укрепляющей части на количество исходной смеси, поступающей на питающую тарелку. Если обозначить количество питания, приходящегося на 1 кмоль дистиллята через f = F/P, то F = P f и количество жидкости, стекающей по исчерпывающей части колонны, составит:

L’ = Ф + F = P R + P f = P (R + f).

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

Для низа колонны состав удаляющейся жидкости (остатка) х’_к = x_w и, согласно допущению, состав поступающего сюда из кипятильника пара y’_н = y_W = x_W. Подставив значения L’, G’, х’_к и y’_н в уравнение (11 a), получим

.

После приведения к общему знаменателю и сокращения подобных членов находим:

. (15)

Зависимость (15) представляет собой уравнение рабочей линии исчерпывающей части колонны. В этом уравнении – тангенс угла наклона рабочей линии к оси ординат, а – отрезок, отсекаемый рабочей линией на оси абсцисс (см. рис. 17).

Умножив числитель и знаменатель выражений для А’ и А (для укрепляющей части колонны) на количество дистиллята Р, можно заметить, что они представляют собой отношения количеств жидкой и паровой фаз, или удельный расход жидкости, орошающей данную часть колонны.

Построение рабочих линий на диаграмме у – х. Для построения рабочих линий откладывают на оси абсцисс диаграммы (см. рис. 17) заданные составы жидкостей x_w, x_F и х_Р. Учитывая принятые допущения о равенстве составов пара и жидкости на концах колонны, из точки х_Р восстанавливают вертикаль до пересечения с диагональю диаграммы в точке а с координатами у_Р = х_р.

Величину R считаем известной. Откладывая на оси ординат отрезок [см. уравнение (14)], соединяют прямой конец отрезка (точку d) с точкой а. Из точки, отвечающей заданному составу x_F, проводят вертикаль до пересечения с линией ad в точке b. Прямая ab – рабочая линия укрепляющей части колонны. Согласно допущению y_w = x_w из точки, соответствующей составу x_w, восстанавливают вертикаль до пересечения с диагональю диаграммы и получают точку с – конечную точку рабочей линии исчерпывающей части колонны. Соединяют точку с прямой с точкой b, принадлежащей одновременно рабочим линиям укрепляющей и исчерпывающей частей колонны. Прямая bс представляет собой рабочую линию исчерпывающей части колонны.

Рабочие линии ab и bс в отличие от рабочих линий процесса абсорбции располагаются под линией равновесия. В данном случае, как уже отмечалось, НК переходит в паровую фазу, стремящуюся к равновесию с жидкой фазой, т. е. по существу десорбируется из жидкости.

Минимальное и действительное флегмовое число

Расчет минимального флегмового числа. При заданном составе дистиллята х_Р величина отрезка В (см. рис. 17), отсекаемого рабочей линией укрепляющей части колонны на оси ординат, зависит только от флегмового числа R, так как . С уменьшением R отрезок В увеличивается (В» > В) и рабочая линия как бы поворачивается вокруг точки а по часовой стрелке, занимая последовательно положения ab, ab» и т. д. Однако величину R можно уменьшать только до некоторого предела, определяемого движущей силой процесса массопередачи между жидкой и паровой фазами.

Движущая сила, выраженная в концентрациях паровой фазы, изображается на диаграмме у – х вертикальным отрезком между данной точкой на рабочей линии и линией равновесия. Например, при рабочей линии ab в точке ввода питания (x_F) движущая сила равна y * _F – у_F и изображается отрезком b»‘b. С уменьшением R точка b перемещается по вертикали, соответствующей абсциссе точки, которая отвечает составу x_F, и движущая сила снижается до тех пор, пока не обратится в нуль (точка b»‘). При этом рабочая линия ab«‘ отсекает на оси ординат максимальный отрезок В»‘ = В_max, которому при заданном х_Р соответствует минимальное флегмовое число R_mln:

.

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

С увеличением R отрезки В уменьшаются и рабочая линия поворачивается вокруг точки а против часовой стрелки. Очевидно, нижнее предельное положение рабочих линий должно соответствовать совпадению точки их пересечения с диагональю диаграммы (точка b’). При этом угол наклона рабочих линий к оси абсцисс равен 45°, А = А’ = 1 и В = В’ = 0, что возможно, как следует из выражений для В и В’, только при бесконечно большом флегмовом числе (R = ¥).

Действительное (рабочее) флегмовое число R_д, при котором работает колонна, должно находиться в пределах R_min и R = ¥. Исходной величиной для выбора действительного флегмового числа является R_min, значение которого можно найти расчетом.

Для определения R_min проведем из точки b» (см. рис. 17) горизонтальный отрезок b’е до пересечения с ординатой точки а. Тангенс угла наклона рабочей линии укрепляющей части колонны при R_min равен отношению катетов ае и b»e треугольника ab»‘e, причем катет ае = у_Р – y_f = х_Р – y’_F, а катет b»‘е = х_Р – x_F. Следовательно

. (А)

Вместе с тем, согласно уравнению (14), при минимальном флегмовом числе

. (Б)

Сопоставляя выражения (А) и (Б), получим

. (16)

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

Эксплуатационные расходы, определяемые расходом теплоносителя, возрастают прямо пропорционально величине R (рис. 18, кривая 1). Более сложной является зависимость капитальных затрат от величины флегмового числа. С увеличением R возрастает движущая сила процесса и уменьшается необходимое число теоретических и соответственно действительных ступеней. В итоге при некотором флегмовом числе рабочий объем колонны станет минимальным и, следовательно, минимальной будет ее стоимость. Поэтому зависимость капитальных затрат от флегмового числа имеет минимум (кривая 2). Отсюда следует, что суммарные затраты будет также иметь минимум, который не совпадает с минимумом капитальных затрат. Зависимость суммарных затрат З (в рублях) от флегмового числа изображается на рисунке кривой 3. Этому минимуму суммарных затрат соответствует оптимальное значение действительного флегмового числа (R_опт).

В связи со сложностью технико-экономического расчета R_oпm выбор действительного флегмового числа R_д часто производят приближенно. Так, при расчетах задаются отношением действительного флегмового числа к минимальному. Это отношение носит название коэффициента избытка флегмы:

.

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

Зависимость между флегмовым числом, высотой колонны и расходом теплоносителя (греющего пара).Рассмотрим, как связана величина флег-мового числа с рабочей высотой колонны и расходом тепла на ректификацию в двух предельных случаях: R = R_min и R = ¥. Рабочая высота колонны пропорциональна числу теоретических ступеней изменения концентрации, которое определяется построением «ступенек» между рабочими линиями и равновесной линией.

При R = ¥ рабочие линии совпадают с диагональю диаграммы и движущая сила процесса или является наибольшей, а необходимое число теоретических ступеней – наименьшим (рис. 19, а). Количество действительных ступеней разделения пропорционально числу теоретических ступеней. Таким образом, при R = ¥ потребовалась бы наименьшая рабочая высота колонны. Однако флегмовое число R = Ф / Р может стать равным бесконечности только при Р = 0. Это означает, что при R = ¥ отбора дистиллята нет, и вся жидкость, полученная в результате полной конденсации паров в дефлегматоре, возвращается в колонну в виде флегмы. В данном случае колонна работает «на себя», без выдачи продукта, что в нормальных производственных условиях, естественно, исключается. Подобный режим работы колонны удобен только для исследовательских целей.

С увеличением R возрастает количество жидкости, которое необходимо испарить в кипятильнике. При R = ¥ требуется испарить максимально возможное количество жидкости. Следовательно, в этом случае расход греющего пара наибольший.

При R_min (рис. 19, б), когда рабочие линии пересекаются с линией равновесия, в точке пересечения движущая сила равна нулю. Значит, для того чтобы достигнуть концентраций фаз, соответствующих их составам на питающей тарелке, потребовалась бы бесконечно большая поверхность контакта фаз, т. е. бесконечно большое число «ступенек» – теоретических ступеней разделения. Таким образом, при R_min разделение возможно только в гипотетической ректификационной колонне бесконечно большой высоты. При этом расход греющего пара, которые при прочих равных условиях пропорционален флегмовому числу, т. к. G = Р (R + 1), будет наименьший.

На основе проведенного анализа можно заключить, что с увеличением флегмового числа высота аппарата уменьшается, а расход греющего пара возрастает. Вместе с тем с увеличением R возрастает количество орошающей жидкости и диаметр аппарата (при прочих равных условиях) увеличивается.

Не нашли, что искали? Воспользуйтесь поиском:

Расчёт оптимального флегмового числа

Расчёт ректификационной колонны

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

Для нахождения расхода дистиллята и флегмы составим материальный баланс для колонны:

F – молярный расход питания, [кмоль/с];

D – молярный расход дистиллята, [кмоль/с];

W – молярный расход кубового остатка, [кмоль/с].

x_F – молярная доля низкокипящего компонента (НК) в питании, [кг/кмоль];

x_D – молярная доля НК в дистилляте, [кг/кмоль];

x_W — молярная доля НК в кубовом остатке, [кг/кмоль].

Решая систему уравнений (1) и (2), можно найти расходы дистиллята и кубового остатка:

В задании даны массовые доли НК в питании, дистилляте и кубовом остатке. Чтобы найти расходы дистиллята и кубового остатка нужно преобразовать массовые доли в молярные. Дано:

– массовая доля НК в питании, [кг/кг];

– массовая доля НК в дистилляте, [кг/кг];

– массовая доля НК в кубовом остатке, [кг/кг].

M_НК – молярная масса НК, [кг/кмоль];

M_ВК – молярная масса ВК, [кг/кмоль].

По заданию в ректификационной колонне будет происходить разделение бинарной смеси хлорбензол-бензол. Температура кипения бензола – 80,1 °С, хлорбензола – 131 °С. Следовательно НК – бензол, ВК – хлорбензол. Молярная масса бензола – 78,11 г/моль, хлорбензола – 112,56 г/моль. Подставим значения в формулы.:

Найдём молярные расходы дистиллята и кубового остатка. Для этого найдём молярный расход питания, т.к. в задании дан массовый расход:

Связь массового и молярного расхода определяется по следующей зависимости:

M_F – молярная масса питания, [кг/кмоль].

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

Выразим молярный расход питания и подставим в выражение значение молярной массы питания:

Расчёт минимального флегмового числа

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

– молярная доля бензола в паре, равновесном с жидким питанием, [кмоль/кмоль].

Для того, чтобы найти необходима таблица со справочными данными:

Таблица 1. Данные по равновесию жидкость – пар для смеси хлорбензол-бензол при атмосферном давлении

Т.к. x_F известна, можно определить графическим методом y_F, но метод интерполяции в данном случае будет точнее:

Теперь по формуле (11) рассчитаем минимальное флегмовое число:

Для проверки расчёта R_min построим диаграмму x – y и обозначим там точки: А (x_D; y_D=x_D), B (x_W; y_W=x_W), C (x_F=y_F). Чтобы построить более точную диаграмму, воспользуемся пакетом MatLab. С помощью следующего кода у нас получится диаграмма с обозначенными на ней точками:

x=[0 0.104 0.192 0.295 0.296 0.399 0.514 0.591 0.684 0.703 0.786 0.804 0.884 1]

y=[0 0.311 0.48 0.628 0.629 0.731 0.816 0.86 0.904 0.912 0.942 0.95 0.971 1]

plot(x, y, xr, yr, x0, y0, 0.32, 0.653, ‘O’, 0.04, 0.04, ‘O’, 0.99, 0.99, ‘O’), grid, xlabel ‘x, кмоль НК/кмоль смеси’, ylabel ‘y, кмоль НК/кмоль смеси’, gtext ‘y*=f(x)’, gtext ‘y=x’

text (0.32+0.01, 0.653-0.01, ‘C’)

text (0.04+0.01, 0.04-0.01, ‘A’)

text (0.99-0.01, 0.99-0.025, ‘B’)

legend(‘Равновесная линия’, ‘Рабочая линия’, ‘y=x’)

Рисунок 2.1. Диаграмма равновесия x – y для смеси хлорбензол-бензол при атмосферном давлении (P=0,1 Мпа)

Как видно по диаграмме, рабочая линия не пересекает равновесную, а касается её в точке С, следовательно, R_min рассчитано верно.

Расчёт оптимального флегмового числа

Оптимальное флегмовое число R_опт определяется путём построения и анализа графика функции N_T(R+1)=f(R). N_T – число теоретических тарелок, обеспечивающих заданное разделение смеси в колонне при данном значении R.

Примем, что флегмовое число изменяется в пределах от R¢ до R¢¢. Значения R¢ и R¢¢ рассчитываются по формулам:

В диапазоне значений изменения R установим n_R=4 промежуточных с интервалом DR:

После несложных расчётов получим следующие числа: R₁=1,352, R₂=1,644, R₃=1,936, R₄=2,228.

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

По уравнению (15) рассчитываем y_в, x_в=x_F=0,32 кмоль/кмоль, а рабочие линии исчерпывающей и укрепляющей колонны пересекаются в точке с координатами [x_в; y_в]. Строим диаграммы x – y и определяем количество теоретических тарелок для каждого флегмового числа. Опять-таки обратимся за помощью к пакету MatLab. Для R=R¢:

Код программы для вывода графиков:

x=[0 0.104 0.192 0.295 0.296 0.399 0.514 0.591 0.684 0.703 0.786 0.804 0.884 1]

y=[0 0.311 0.48 0.628 0.629 0.731 0.816 0.86 0.904 0.912 0.942 0.95 0.971 1]

plot(x, y, xr, yr, x0, y0, 0.32, 0.645, ‘O’, 0.04, 0.04, ‘Og’, 0.99, 0.99, ‘Og’), grid, xlabel ‘x, кмоль НК/кмоль смеси’, ylabel ‘y, кмоль НК/кмоль смеси’, gtext ‘y*=f(x)’, gtext ‘y=x’

text (0.32+0.01, 0.645-0.01, ‘C’)

text (0.04+0.01, 0.04-0.01, ‘A’)

text (0.99-0.01, 0.99-0.025, ‘B’)

legend(‘Равновесная линия’, ‘Рабочая линия’, ‘y=x’)

yr1=[0.04 0.605 0.99]

plot(x, y, xr1, yr1, x0, y0, 0.32, 0.605, ‘O’, 0.04, 0.04, ‘Og’, 0.99, 0.99, ‘Og’), grid, xlabel ‘x, кмоль НК/кмоль смеси’, ylabel ‘y, кмоль НК/кмоль смеси’, gtext ‘y*=f(x)’, gtext ‘y=x’

text (0.32+0.01, 0.605-0.01, ‘C1’)

text (0.04+0.01, 0.04-0.01, ‘A’)

text (0.99-0.01, 0.99-0.025, ‘B’)

legend(‘Равновесная линия’, ‘Рабочая линия’, ‘y=x’)

yr2=[0.04 0.573 0.99]

plot(x, y, xr2, yr2, x0, y0, 0.32, 0.573, ‘O’, 0.04, 0.04, ‘Og’, 0.99, 0.99, ‘Og’), grid, xlabel ‘x, кмоль НК/кмоль смеси’, ylabel ‘y, кмоль НК/кмоль смеси’, gtext ‘y*=f(x)’, gtext ‘y=x’

text (0.32+0.01, 0.573-0.01, ‘C2’)

text (0.04+0.01, 0.04-0.01, ‘A’)

text (0.99-0.01, 0.99-0.025, ‘B’)

legend(‘Равновесная линия’, ‘Рабочая линия’, ‘y=x’)

yr3=[0.04 0.548 0.99]

plot(x, y, xr3, yr3, x0, y0, 0.32, 0.548, ‘O’, 0.04, 0.04, ‘Og’, 0.99, 0.99, ‘Og’), grid, xlabel ‘x, кмоль НК/кмоль смеси’, ylabel ‘y, кмоль НК/кмоль смеси’, gtext ‘y*=f(x)’, gtext ‘y=x’

text (0.32+0.01, 0.548-0.01, ‘C3’)

text (0.04+0.01, 0.04-0.01, ‘A’)

text (0.99-0.01, 0.99-0.025, ‘B’)

legend(‘Равновесная линия’, ‘Рабочая линия’, ‘y=x’)

yr4=[0.04 0.526 0.99]

plot(x, y, xr4, yr4, x0, y0, 0.32, 0.526, ‘O’, 0.04, 0.04, ‘Og’, 0.99, 0.99, ‘Og’), grid, xlabel ‘x, кмоль НК/кмоль смеси’, ylabel ‘y, кмоль НК/кмоль смеси’, gtext ‘y*=f(x)’, gtext ‘y=x’

text (0.32+0.01, 0.526-0.01, ‘C4’)

text (0.04+0.01, 0.04-0.01, ‘A’)

text (0.99-0.01, 0.99-0.025, ‘B’)

legend(‘Равновесная линия’, ‘Рабочая линия’, ‘y=x’)

plot(x, y, xrm, yrm, x0, y0, 0.32, 0.51, ‘O’, 0.04, 0.04, ‘Og’, 0.99, 0.99, ‘Og’), grid, xlabel ‘x, кмоль НК/кмоль смеси’, ylabel ‘y, кмоль НК/кмоль смеси’, gtext ‘y*=f(x)’, gtext ‘y=x’

text (0.32+0.01, 0.51-0.01, ‘Cm’)

text (0.04+0.01, 0.04-0.01, ‘A’)

text (0.99-0.01, 0.99-0.025, ‘B’)

legend(‘Равновесная линия’, ‘Рабочая линия’, ‘y=x’)

Теоретическое количество тарелок для ректификационной колонны определим по графическому методу Мак-Кэба и Тиле.

Рисунок 2.2. Диаграмма x – y при R=R¢. N_T=15.

Рисунок 2.3. Диаграмма x – y при R=R₁. N_T=11.

Рисунок 2.4. Диаграмма x – y при R=R₂. N_T=10.

Рисунок 2.5. Диаграмма x – y при R=R₃. N_T=9,5.

Рисунок 2.6. Диаграмма x – y при R=R₄. N_T=9.

Рисунок 2.7. Диаграмма x – y при R=R¢¢. N_T=8.

Для более удобного построения графика зависимости N_T(R+1)=f(R) составим следующую таблицу:

Таблица 2. Данные для построения графика зависимости N_T(R+1)=f(R)

Для построения графика N_T(R+1)=f(R) воспользуемся уже хорошо нам знакомым пакетом MatLab:

R=[1.06 1.352 1.644 1.936 2.228 2.52]

N=[30.9 25.872 26.44 27.833 29.052 28.16]

plot(R, N, ‘or’, x,y), grid, xlabel ‘R’, ylabel ‘N(R+1)’

Как видно из графика R_опт=1,43. Для дальнейших расчётов примем рабочее флегмовое число R=R_опт.

Уравнения рабочих линий

Уравнения рабочих линий укрепляющей и исчерпывающей частей колонны принимают вид:

f – относительный молярный расход питания.

Уравнения рабочих линий укрепляющей и исчерпывающей частей колонны для данной установки: