Спектр применения фазовых диаграмм растворимости достаточно широк: химическое производство, обслуживание водных объектов, моделирование процессов добычи полезных ископаемых и т.д.
Математически фазовую диаграмму растворимости для пятирной водно-солевой системы можно определить следующим образом: диаграмма представляется в виде четырёхмерного симплекса; каждая вершина этого симплекса соответствует 100 % доле компонента этой системы; координаты любой точки внутри симплекса (фигуративной точки) определяют состав исходной реакционной смеси (СИРС) в процентном соотношении.
В ходе химических экспериментов для определённых пятирных водно-солевых растворов были найдены экспериментальные точки, являющиеся границами твёрдой и жидкой фаз, входящих в него солей. В симплексе эти точки задают искривлённые поверхности.
Таким образом, сформулированы следующие задачи:
1) построение непрерывных приближений поверхностей, заданных набором экспериментальных точек;
2) определение положения произвольных фигуративных точек относительно областей, ограниченных поверхностями внутри симплекса.
Решение обеих задач требует их автоматизации. Поэтому появилась необходимость разработать программный продукт, позволяющий работать с фазовыми диаграммами растворимости для пятирных водно-солевых систем простого эвтонического типа.
Работа была разбита на две подзадачи:
1) разработать программную реализацию выбранной математической модели;
2) по полученной диаграмме построить прогноз фаз для произвольных точек.
Попытки реализации уже предпринимались на кафедре неорганической химии химического факультета ПГУ, однако существующий программный продукт обладает рядом недостатков: с одной стороны, это недостаточный объём реализованной функциональности, а с другой - отсутствие анализа реализованных алгоритмов.
По этой причине было проведено теоретическое исследование. Ранее было установлено, что для построения математических моделей изотрем растворимости водно-солевых систем целесообразно использовать «мозаичные» поверхности. Поэтому для решения первой задачи за основу взят метод, основанный на построении триангуляции Делоне. Вычислительная сложность этого этапа составит в худшем случае O(n2) от числа экспериментальных точек.
За счёт разбиения объёма симплекса на канонические фигуры удаётся свести вторую задачу к определению положения точки по отношению к симплексу. Вычислительная сложность этого этапа составит O(d3 c), где d = 4 - размерность пространства, c - количество граней поверхности, полученных после триангуляции.
В результате проведенного исследования предложены алгоритмы для решения поставленных задач и получены теоретические оценки вычислительной сложности этих алгоритмов и времени работы программы.