Геометрия фазовых диаграмм четверных водно-солевых систем простого эвтонического типа

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

2

4

Пермский государственный национальный исследовательский университет

Россия, 614990, Пермь, ул. Букирева,

Геометрия фазовых диаграмм четверных водно-солевых систем простого эвтонического типа

М. С. Моисеенков,

С. А. Мазунин,

В. Л. Чечулин

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

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

растворимость солевая система кристаллогидрат

Введение© Моисеенков М.С., Мазунин С.А., Чечулин В.Л., 2012

Состояния многокомпонентной физико-химической системы описываются правилом фаз Гиббса, принципами соответствия и непрерывности, предложенными Курнаковым [1].

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

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

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

1. Фазовые диаграммы водно-солевых систем простого эвтоничсекого типа

1.1 Тройные водно-солевые системы простого эвтонического типа

Тройная водно-солевая система простого эвтнического типа согласно правилу фаз и принципам соответствия и непрерывности состоит из трех компонентов и включает восемь фазовых областей. Фазовая диаграмма ее изображается в виде равностороннего треугольника, образованного вершинами (S₁, S₂ и H₂O). В нем боковые ребра S₁-H₂O и S₂-H₂O являются оконтуривающими двойными системами. Точки R₁ и R₂ соответствуют растворимостям в этих системах при данной температуре, точка E - эвтоника - составу насыщенного при постоянном давлении и температуре раствора, находящегося в нонвариантном равновесии с двумя или более твердыми фазами [2].

Внутреннее пространство треугольника (рис. 2) делится на четыре фазовые области:

1) область совместной кристаллизации двух солей (рис. 3, а);

2) две области кристаллизации одной соли (рис. 3, б):

область кристаллизации соли S₁ - область S₁R₁E,

область кристаллизации соли S₂ - область S₂R₂E;

3) область ненасыщенных растворов (рис. 3, в).

Кривые R₁E и R₂E называются ветвями растворимости и являются линиями моновариантного равновесия солей S₁ и S₂ и соответствующих насыщенных растворов. В настоящее время не существует аналитического описания этих кривых, поэтому они строятся приближенно по имеющимся экспериментальным данным.



1.2 Четверные водно-солевые системы простого эвтонического типа

Изотерма и изобара четверной системы изображается уже в виде тетраэдра в трехмерном пространстве состояния (рис. 4, а).

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

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

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

1) область совместной кристаллизации всех трех солей - пирамида S₁S₂S₃E (рис. 4, б);

2) три области кристаллизации двух солей:

область совместной кристаллизации солей S₁ и S₂ (рис. 4, г),

Рис. 4. Фазовые области ФДР четверной водно-солевой системы простого эвтонического типа

область совместной кристаллизации солей S₁ и S₃ (рис. 4, д),

область совместной кристаллизации солей S₂ и S₃ (рис. 4, е);

3) три области кристаллизации одной соли:

область кристаллизации соли S₁ (рис.4,ж),

область кристаллизации соли S₂ (рис.4, з),

область кристаллизации соли S₃ (рис.4, и);

4) область ненасыщенных растворов - область H₂O-R₁E₁₂R₂E₂₃R₃E₁₃ (рис. 4, в).

Линии моновариантного равновесия E₁₂E, E₁₃E, E₂₃E являются кривыми, восстановленными по экспериментальным данным.

Поверхности кристаллизации R₁E₁₂EE₁₃, R₂E₂₃EE₁₂ и R₃E₁₃EE₂₃ солей S₁, S₂ и S₃ соответственно также не имеют аналитического описания и на практике строятся по данным экспериментов.

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

2. Структура фазовых диаграмм водно-солевых систем простого эвтонического типа

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

Для построения фазовых диаграмм растворимости простых эвтонических систем применяют метод Гиббса-Розебома [1], использующий барицентрические координаты.

Диаграмма k-компонентной водно-солевой системы изображается в виде d-мерного симплекса, каждая вершина которого соответствует компоненту системы, причем

. (1)

Любая точка внутри симплекса соответствует некоторой смеси компонентов, причем состав этой смеси определяет координаты этой точки. Пусть данной смеси X соответствуют массовые концентрации компонентов системы x₁, x₂,…, x_k, тогда сумма этих концентраций всегда составляет 100%. Поэтому, зная солевой состав, всегда можно вычислить долю воды:

.(2)

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

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

Пусть t - количество солей, находящихся в твердой фазе (t=0, 1, …, k-1). Тогда количество различных фазовых областей, содержащих t твердых фаз, равно числу сочетаний:

. (3)

Общее количество фазовых областей на фазовой диаграмме простой k-компонентной системы равно в изобарно-изотермических условиях:

. (4)

Описание структуры диаграммы сводится к определению ее фазовых областей. Рассмотрим фазовые области, в которых хотя бы одна соль находится в твердой фазе (t=1, 2, …, k-1). Каждая такая область содержит t точек состава компонентов, отвечающих солям, выпавшим в осадок, и точки, образующие геометрические образы фазовых равновесий различной вариантности.

Размерность этих образов можно определить с помощью правила фаз Гиббса [2], с учетом постоянного давления и температуры в системе:

. (5)

Здесь f - число фаз, v - число степеней свободы (или вариантность) [2]. Так как число фаз в данном случае равно t+1, число степеней свободы

. (6)

Составы фаз, участвующих в нонвариантном равновесии (v=0), отображается точками, в моновариантных равновесиях (v=1) - линиями, в дивариантных равновесиях (v=2) - поверхностями, в тривариантных и более (v?3) - гиперповерхностями. Поливариантные равновесия не имеют аналитического описания и строятся приближенно по наборам точек, которые находятся на этих поверхностях.

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

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

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

1) четырехфазная область кристаллизации трех солей (рис. 4, б);

2) три трехфазные области кристаллизации двух солей (рис. 4, г, д, е);

3) три двухфазные области кристаллизации одной соли (рис. 4, ж, з, и);

4) область ненасыщенных растворов (рис. 4, в).

На каждом шаге проверки фазовые области кристаллизации t солей (t=1, 2, …, k-2) разбиваются на совокупность непересекающихся d-мерных симплексов, и определяется положение точки F относительно каждого из них. Если фигуративная точка принадлежит хотя бы одному из симплексов разбиения, то делается вывод о принадлежности этой точки данной фазовой области.

2.1 Определение положения фигуративной точки относительно симплекса

Пусть фигуративная точка F имеет координаты в d-мерном пространстве:

.

Симплекс задан d+1 вершинами

.

Каждая грань симплекса S образована d его вершинами, определяющими d-мерную гиперплоскость G_i, (i=1, 2, …, d+1), описываемую уравнением

. (7)

Вычислив расстояния от грани G_i симплекса S до точки F и до вершины S_i этого симплекса, противолежащей данной грани, по формулам

,(8)

,(9)

запишем систему неравенств:

. (10)

Если все неравенства (10) выполняются строго, то фигуративная точка F, соответствующая ИСК, находится внутри данного симплекса. Если () неравенств системы (10) обращаются в равенство, то это означает, что фигуративная точка находится на грани d-мерного симплекса (при n = 1), на грани меньшей размерности, являющейся пересечением n граней симплекса (при 1 < n < d) или в одной из его вершин (при n=d).

2.2 Разбиение фазовых областей и вычисление равновесных фаз

В соответствии с предложенной последовательностью проверки областей кристаллизации фазовой диаграммы первой рассматривается четырехфазная область кристаллизации трех компонентов: S₁, S₂ и S₃.

Единственным равновесным составом в этой области является нонвариантный раствор E. Три другие твердые фазы, хотя и имеют постоянный состав, но могут присутствовать в донной фазе в самых разнообразных соотношениях. Соединив ее предельными нодами со всеми вершинами солевого основания (равновесными составами твердых фаз), получим четырехфазную область кристаллизации трех солей, геометрическим образом которой является трехмерный симплекс. Поэтому проверка принадлежности фигуративной точки F сводится к проверке принадлежности этой точки симплексу S₁S₂S₃E.

Если ИСК находится в четырехфазной области кристаллизации трех компонентов: S₁, S₂ и S₃ (рис. 5), то состав равновесной жидкой фазы является постоянным и соответствует тройному эвтоническому раствору (точка E), а состав донной фазы - переменным и отвечает точке F_S на солевом основании диаграммы. Соотношение твердых и жидких фаз определяется по правилу центра масс или путем последовательного применения правила рычага.

Для трех объемов совместной кристаллизации двух солей моделирование сводится к нахождению функциональной зависимости линии моновариантного равновесия, а также к определению принадлежности ИСК данной фазовой области. Линия моновариантного равновесия E₁₂E представляется ломаной, составленной из отрезков E₁₂M₁, M₁M₂, …, M_n-1M_n, M_nE, где M₁, M₂, …, M_n - экспериментальные моновариантные точки на данной линии или вычисленные интерполяционные значения.

Зафиксировав точки S₁ и S₂, и последовательно взяв по две ближайшие друг к другу моновариантные точки, получим разбиение трехфазной области на тетраэдры S₁S₂E₁₂M₁, S₁S₂M₁M₂, …, S₁S₂M_n-1M_n, S₁S₂M_nE (рис. 6).

Далее решается задача последовательной проверки принадлежности данной фигуративной точки каждому из (n+1) объемов. Если ИСК не принадлежит ни одному объему, то он не находится в данной фазовой области. Если же среди полученных объемов будет найден такой трехмерный симплекс S₁S₂M_iM_j, который включит точку F, то равновесная жидкая фаза является решением СЛАУ [2]:

(11)

где F(x₀, y₀, z₀) - фигуративная точка, F_L(x_L, y_L, z_L) - точка, соответствующая равновесной жидкой фазе, M_i(x_i, y_i, z_i) и M_j(x_j, y_j, z_j) - экспериментальные точки на линии моновариантного равновесия двояконасыщенных растворов.

По формулам (12) и (13) вычисляется равновесная донная фаза F_S(x_S, y_S, z_S):

, (12)

x_S + y_S = 100, z_S =0. (13)

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

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

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

Обозначим множество точек поверхности насыщенных растворов компонентом S₁ через M₁, M ₂, ..., M _i, ..., M _n-1, M _n. В первом приближении задача сводится к нахождению ближайших трех точек поверхности насыщенных растворов от точки состава перспективной проекции данной исходной смеси из вершины соответствующего солевого компонента на поверхность насыщенных растворов (рис. 7). В качестве такого критерия "близости" нами был выбран косинус угла между двумя прямыми в трехмерном пространстве (заданными в виде векторов с соответствующими координатами: a={x_a; y_a; z_a}, b_i={x_b; y_b, z_b}), определяемый по формуле

. (14)

При этом a=, (i=1, 2, …, n), b=. Чем меньше угол между прямыми, т.е. чем ближе значение к 1, тем ближе находится точка поверхности M_i к точке пересечения луча кристаллизации (прямая S₁F).



Рис. 7. Определение состава насыщенного раствора при использовании мозаичной модели поверхности кристаллизации

Вычислив значения всех n косинусов, выберем такие три ближайшие к F точки M_i, M_j и M_k, что

1) косинусы соответствующих им углов были наибольшими;

2) эти точки не лежали на одной прямой;

3) эти точки не лежали на одной линии моновариантного равновесия.

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

Равновесной твердой фазе соответствует вершина солевого компонента S₁, а равновесной жидкой фазе - точка F_L, точка пересечения прямой S₁F и плоскости M_iM_jM_k (рис.7).

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

3. Вычисление физико-химических характеристик в точках на границе ликвидуса

Пусть некоторая фигуративная точка F попала в фазовую область, в которой хотя бы одна из солей находится в твердой фазе. Тогда для нее могут быть вычислены координаты точки, отвечающей составу равновесной жидкой фазы F_L(x₁,x₂,…,x_k). Эта точка лежит на поверхности ликвидуса.

Если в точках ликвидуса известны физико-химические характеристики (p), такие как плотность или показатель преломления, то требуется вычислить их и для новой точки - F_L, как усредненного значения по ближайшим точкам.

Если точка равновесной жидкой фазы принадлежит v-вариантному геометрическому образу фазового равновесия, то существует такой минимальный v-мерный симплекс S в d_мерном пространстве, которому принадлежит эта точка:

. (15)

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

. (16)

Приравняв определитель этой матрицы к нулю и выразив p, получим формулу для вычисления заданной характеристики:

,(17)

где A_ij - алгебраическое дополнение элемента a_ij матрицы A.

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

Приведенное описание геометрической структуры фазовой диаграммы простых эвтонических водно-солевых систем без образования кристаллогидратов и ограниченных рядов твёрдых растворов позволяет построить объектную модель, программная реализация которой автоматизирует вычисления, связанные с обработкой информации о системе. При этом наглядное представление фазовых диаграмм с числом компонентов до 4 является важной прикладной задачей, поскольку визуальная обработка графической информации играет огромную роль при интерпретации результатов исследований. Поскольку в настоящее время далеко не все графические библиотеки поддерживают визуализацию в барицентрических координатах, возникла задача графического представления фазовых диаграмм четверных водно-солевых систем на основе декартовой системы координат. Координатным симплексом четверной системы служит трехмерный симплекс - тетраэдр. Его вершины имеют координаты H₂O(0;0;0), S₁(100;0;0), S₂(0;100;0), S₃(0;0;100). В декартовой системе координат этот симплекс не является правильным, поскольку прямые H₂O-S₁, H₂O-S₂, H₂O-S₃ взаимно ортогональны и записаны в ортонормированном базисе:

. (18)

Для того чтобы точки симплекса образовывали правильный тетраэдр, вершиной которого была бы точка H₂O, а основанием - правильный треугольник , применяется аффинное преобразование координат, с матрицей T перехода в новый базис:

. (19)

Используя приведенное описание геометрической структуры ФДР и преобразование координат (19) был разработан прототип модуля визуализации четверных ФДР в пакете Wolfram Mathematica 7. Модуль визуализации позволяет строить ФДР четверных водно-солевых систем простого эвтонического типа без образования кристаллогидратов и ограниченных рядов твердых растворов по входным экспериментальным данным (рис. 8).

5. Приложение визуализации

Одним из приложений визуализации диаграмм растворимости многокомпонентных водно-солевых систем является подтверждение и проверка открытого недавно коллигативного свойства плоскостности линий моновариантного равновесия [3_6].

Фундаментальное объяснение свойства плоскостности заключается в преимущественно физическом взаимодействии молекул в растворе [7].

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

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

Рис. 9. Плоскость, построенная по оконтуривающим эвтоникам

Заключение

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

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

1. Курнаков Н.С. Введение в физико-химический анализ. Л.; М.: Изд-во АН СССР, 1940.

2. Мазунин С.А. Основы физико-химического анализа: учеб. пособие по спецкурсу / Перм. ун-т. Пермь, 2000. Ч.2. Многокомпонентные водно-солевые системы. 212 с.

3. Чечулин В.Л., Мазунин С.А. О плоскостности координат точек моно- и нон-вариантных равновесий в 4-х и более компонентных водно-солевых системах // Изв. высш. учеб. заведений: Химия и хим. технология. 2010. Т.53, №3. С.152-154.

4. Чечулин В.Л., Мазунин С.А., Фролова С.И., Кистанова Н.С. Новое коллигативное свойство многократно насыщенных нон- и моновариантных водных растворов // IX Междунар. Курнаковское совещ. по физ.-хим. анализу: тез. докл. Пермь, 2010. С. 59.

5. Чечулин В.Л., Мазунин С.А. О плоскостности моно- и нонвариантных равновесий как коллигативном свойстве многократно насыщенных водных растворов // Журн. общей химии. 2012. Т. 82, №.2. С.202-204.

6. Chechulin V.L., Mazunin S.A. Planarity of Mono- and Non-Variant Equilibria as the Colligative Property of Multicomponent Saturated Water Solutions // Russian Journal of General Chemistry. 2012. Vol. 82, №2. P.199-201.

7. Чечулин В.Л. К вероятностному описанию многократно насыщенных растворов в многокомпонентных водно-солевых системах // IX Междунар. Курнаковское совещ. по физ.-хим. анализу: тез. докл. Пермь, 2010. С. 86.

8. Кистанова Н.С., Мазунин С.А., Фролова С.И. Прогнозирование составов многократно насыщенных нонвариантных растворов // Там же. С.52.

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

...