APPROXIMATION OF SOME EQUILIBRIUM CURVES OF BINARY SYSTEMS BY ANALYTICAL DEPENDENCE

Research article
DOI:
https://doi.org/10.60797/IRJ.2025.155.83
Issue: № 5 (155), 2025
Suggested:
17.11.2024
Accepted:
25.04.2025
Published:
16.05.2025
112
4
XML
PDF

Abstract

The article examines the vapour-liquid equilibrium curves for a number of physicochemical systems. Binary systems play an important role in physicochemical processes such as absorption, rectification, distillation and others. The aim of the work is to select a formula that is approximating experimental data for a large set of physico-chemical vapour-liquid systems, and allows describing equilibrium curves for the noted processes in binary systems. The proposed formula is a mathematical model that indirectly takes into account the influence of physicochemical properties of components, temperature and pressure on the equilibrium distribution of components in the process. The task of the research was also the selection of the coefficients in the formula for some systems and the estimation of the error of correspondence between experimental and theoretical data.

The use of this approximating formula allows researchers and engineers to more accurately predict and optimise processes in binary systems. For example, for distillation processes, by optimising parameters such as pressure and temperature, it is possible to change the equilibrium concentrations of the components, which increases the yield of the target product.

The data on approximation of some binary systems at equilibrium are obtained, which are presented in the table, which is scientific novelty of the work, the evaluation of the error of this formula is given, confirming good accuracy. An original function of simple form for describing the process of simple distillation is found. It is of practical importance for the development of effective methods of mixture separation and increasing the energy efficiency of production processes.

1. Введение

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

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

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

2. Основные результаты

В данной работе используется двухпараметрическая формула:

img
(1)

В этом случае х — это мольная доля низкокипящего компонента (НК) в жидкой фазе, у*(х) — равновесная концентрация НК между паром и жидкостью, α и β — параметры, выбор которых позволяет аппроксимировать реальные системы. Соотношение (1) дает возможность аналитически находить определенные интегралы типа числа единиц переноса, встречающиеся при описании процессов перегонки, ректификации и возможно некоторых других. Данное соотношение было предложено и использовано для описания отмеченных процессов в насадочных колоннах в работах

,
Отметим также, что выражение (1) применено в работах
,
для описания ректификации в тарельчатой системе для одной конкретной системы.

Зависимость (1) удовлетворяет необходимым для равновесной функции требованиям: у*(0)=0, у*(1)=1. Отметим, что соотношение (1) не может аппроксимировать системы, имеющие азеотропную точку. Такая аппроксимация нецелесообразна, когда азеотропная точка имеет важное значение для описания конкретного процесса. В этом случае формулу (1) можно использовать в интервале от x = 0 до координаты азеотропной точки.

Заметим, что в частном случае α+β=1 формула (1) сводится к известной зависимости Рауля

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

Отметим также, что аппроксимационное соотношение (1) можно использовать и при расчете процессов (например абсорбции) когда кривая равновесия расположена ниже диагонали y = x (например кривая 3 рис. 1).

Равновесные кривые y*(x) определяемые формулой (2)

Рисунок 1 - Равновесные кривые y*(x) определяемые формулой (2)

Примечание: 1 – α = 1,5, β = 0,1; 2 – α = 1,9, β = 0,3; 3 – α = 2, β = 4; 4 – диагональ y = x (β = 1)

Проверка пригодности формулы (1) на конкретной системе
показала достаточную для практики точность. В связи с этим целью данной работы является распространение этой формулы на ряд других бинарных систем. Задачей исследования будет подбор коэффициентов α и β для некоторых систем и оценка погрешности соответствия экспериментальных и теоретических данных.

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

Таблица 1 - Значения параметров α и β для некоторых систем

Система

α

β

δ

Ацетон-бензол

1,586

0,281

8,177*10-3

Ацетон-вода

1,303

0,017

0,027

Ацетон-этиловый спирт

1,653

0,245

7,646*10-3

Бензол-толуол

0,568

0,411

1,351*10-3

Вода-уксусная кислота

0,697

0,547

2,73*10-3

Метиловый спирт-вода

1,323

0,118

5,863*10-3

Метиловый спирт-этиловый спирт

0,041

0,672

5,903*10-3

Муравьиная кислота-уксусная кислота

1,35

0,62

4,117*10-3

Сероуглерод-четырёххлористый углерод

0,859

0,326

4,456*10-3

Хлороформ-бензол

-0,26

0,679

9,004*10-3

Этилацетат-уксусная кислота

0,374

0,301

4,134*0-3

Бромистый этил-этиловый спирт

1,201

0,019

0,022

Бромистый этил-бензол

0,565

0,296

5,506*10-3

Бромистый этил-гептан

0,97

0,08

4,54*10-3

Акролеин-метилэтилкетон

0,283

0,543

0,017

Аммиак-вода

0,816

0,057

0,012

Там же приведены соответствующие данные об среднеквадратичном отклонении экспериментальных данных от рассматриваемой зависимости (1). Эта величина определяется так:

img
(2)

где N — это число экспериментальных точек, y*(xj) – значение эксперимен­тальной равновесной функции в точке xj. Поиск оптимальных значений параметров α и β проводится с помощью стандартной программы системы MathCad. При этом задавалось начальные приближенные значения этих параметров. В качестве наглядной иллюстрации точности апроксимации формулой (1) на рис. 2 представлена конкретная система.

Равновесная кривая y*(x) бинарной системы хлороформ-бензол

Рисунок 2 - Равновесная кривая y*(x) бинарной системы хлороформ-бензол

Примечание: 1 – равновесная кривая; 2 – диагональ

На рис. 2 можно увидеть аппроксимацию бинарной системы хлороформ-бензол: точки °– это экспериментальные данные, линия 1 – равновесная кривая, полученная с помощью формулы (1) и данных таблицы.

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

В качестве примера использования формулы (1) приведем расчетную зависимость для известного

,
уравнения простой перегонки:

img
(3)

где G — количество жидкости в кубе испарителе, GF — начальное количество перегоняемой смеси, x — состав НК в жидкости, xf — содержание НК в начальной смеси, y* — равновесная концентрация НК в паре. Для уравнения (3) после подстановки в него соотношения (1) и выполнения интегрирования, получаем выражение

img
(4)

Например, для системы бензол-толуол, взяв соответствующие значения параметров α = 0,568 и β = 0411 из таблицы, из формулы (4) следует выражение

img
(5)
Графики функции G/GF, определенные зависимостью (5)

Рисунок 3 - Графики функции G/GF, определенные зависимостью (5)

Примечание: 1 – xf = 0,5; 2 – xf = 0,7; 3 – xf = 0,9

Три кривые, в качестве примеров, приведены на рис. 3. Для кривой 3 пунк­тирная вертикальная линия, соответствующая x = xf = 0,9 иллюстрирует, что при этом G/GF = 1, что имеет место и при других xf, в частности показанных на рис 3 кривые 1 и 2, что вытекает из соотношения (4). При малых значениях состава x, формула (5) показывает, что в окрестности нуля функция G/GF степенная.

3. Заключение

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

Далее, найдена оригинальная (новая) функция простого вида (4) для описания процесса простой перегонки.

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

Article metrics

Views:112
Downloads:4
Views
Total:
Views:112