Каждая технологическая схема химического производства подразумевает под собой комплекс упорядоченных действий (этапов), включающих тепло-, массообменную и реакционную аппаратуру, из которых особенно необходимым является этап химического взаимодействия, которое имеет неотъемлемое воздействие на расходные коэффициенты.
У экспертов, занимающихся конструированием химического производства, часто возникают определенные сложности, поскольку для создания проекта необходимы знания различных дисциплин химического профиля.
В некоторых источниках приведены таблицы стандартных термодинамических функций отдельных веществ. Для любых органических соединений, число которых слишком большое, невозможно найти такую информацию. К тому же реальный химический процесс в основном проходит в условиях, далеких от стандартных, поэтому в данной ситуации необходимо принимать к сведению переход от стандартных к реальным условиям. Данный случай требует расчетов с применением приближенных методов или экспериментальных формул.
Для решения уравнений, полученных при изучении структуры основных характеристик экономической эффективности химико-технологического процесса, необходимо численное представление зависимости выхода целевого продукта от степени превращения основного элемента. Так как в химической кинетике указанная зависимость обычно представлена в неявном виде, возникает потребность систематизации и вывода уравнений, группирующих избирательность и выход со степенью превращения.
Основными при расчетах являются материальный и тепловой балансы, на которые опираются все остальные расчеты: подбор и расчет оборудования, определение расходных коэффициентов, расчет себестоимости продукции и т.д. На самом деле, после составления материального баланса и нахождения всех материальных потоков возможны конструктивные расчеты оборудования [1, 2].
Исходя из вышеперечисленных причин, разработка программного комплекса по подборке константы скорости химических реакций и исследованию адекватности применяемой модели является актуальной задачей.
Требования к программному комплексу
Назначение разрабатываемого комплекса – подборка константы скорости химических реакций и исследование адекватности применяемой модели. Область применения – химическое производство.
Входные данные программы:
– количество параметров, среднее отклонение концентрации, количество опытов, данные, полученные в результате серии опытов – вводятся пользователем;
– порядок реакции – выбирается пользователем или программа считает автоматически.
Программный комплекс должен рассчитывать следующие параметры:
– константы скорости химических реакций;
– критерий Фишера;
– порядок реакции;
– коэффициент корреляции [3].
Все вышеописанные значения выводятся на экран в удобном упорядоченном виде.
Описание и структура программы
Для создания программного комплекса в наше время используется огромное количество языков программирования. Для реализации нашего проекта был выбран С#. Язык программирования С# – это один из распространенных средств разработки в рамках объектно-ориентированного подхода. Данный язык был создан в 1998–2001 гг. командой из инженеров под руководством Андерса Хейлсберга в компании Microsoft как основной язык создания приложений для платформы Microsoft NET Framework [4].
Самыми важными структурными компонентами программы являются модули, в которых находится программный код, и связанные с ними экранные формы. Модуль осуществляет обработку команд пользователя, поступающих с соответствующей формы. Программный код модуля состоит из взаимосвязанных подпрограмм (процедур), каждая из которых выполняет какую-либо конкретную задачу [5].
Программа состоит из следующих форм:
1. Форма Constants (окно расчета константы скорости химических реакций и исследование адекватности применяемой модели). Это основное окно, с которым будет работать пользователь. Здесь происходит ввод исходных данных, расчет всех параметров и вывод конечных результатов на экран.
2. Форма Кинетика. Позволяет увидеть кинетику химических реакций.
Интерфейс программы (рис. 1) представляет собой несколько областей, а его реализация является достаточно простой и интуитивно ясной для пользователей.
Рис. 1. Окно расчета константы скорости химических реакций и исследование адекватности применяемой модели
Рис. 2. Окно расчета константы скорости химических реакций и исследование адекватности применяемой модели. Ввод исходных данных
Основная часть программы разделена на две области – «Измерения» и «Результаты» (рис. 2). В области «Измерения» вводятся данные реакций, количество параметров, среднее отклонение концентрации, количество опытов. В зависимости от количества данных можно добавить / удалить строчки с помощью знака «+» / «–».
При вводе данных можно заметить в области «Результаты» автоматическое появление точек на координатной плоскости с соответствующими координатами, которые мы вводим. После нажатия кнопки «Вычислить» на экране появится график.
В области «Результаты» под графиком есть возможность выбора уравнения скорости реакции. Выбор уравнения скорости реакции можно осуществить, нажав на выпадающий список. Первое уравнения демонстрирует уравнение реакции первого порядка; второе и третье уравнение – второго порядка, два разных случая; четвёртое уравнение – уравнение реакции порядка n, где n не равно 1, n вводится рядом (рис. 3). При n = 2 получается то же самое уравнение, что и втором случае из списка. Если уравнение скорости реакции выбрано неправильно, тогда на экране появится красный крест, который сообщает о некорректности введённых данных для данного порядка. Если возникают такие затруднения, то можно воспользоваться возможностью нахождения порядка автоматически самой программой.
Рис. 3. Окно расчета константы скорости химических реакций и исследование адекватности применяемой модели. Определение порядка реакции
Рис. 4. Окно расчета константы скорости химических реакций и исследование адекватности применяемой модели и дополнительное окно с кинетическими кривыми
Рис. 5. Окно расчета константы скорости химических реакций и исследование адекватности применяемой модели с окончательными результатами
При проведении химической реакции получены кинетические кривые для компонентов. Данные кривые можно посмотреть в области «Результаты», нажав на кнопку Кинетика. На экране появится дополнительное окно с кинетическими кривыми (рис. 4).
В результате работы данного программного комплекса на экран выводятся следующие данные (рис. 5):
– график;
– константы скорости химических реакций;
– критерий Фишера;
– порядок реакции;
– коэффициент корреляции.
Адекватность модели оценивали по критерию Фишера. Сравнивали опытное значение критерия Фишера с табличным. Для нахождения констант уравнений в программу был заложен метод наименьших квадратов.
На рис. 6 изображено графическое представление алгоритма программы, описанной с помощью функциональных блоков данной схемы. Каждый блок соответствует выполнению одного либо нескольких действий.
Рис. 6. Блок-схема работы программы
Для вычисления математической основы использовали данные из источников [3, 6].
Заключение
На основе проведенных химиками-технологами различных экспериментов, выдвигаемых гипотез о механизме реакций и построении кинетических моделей возникает потребность исследования, состоящая в обработке экспериментальных данных в следующей последовательности:
1) находим численные значения констант скоростей химических реакций;
2) проводим проверку адекватности полученных уравнений эксперименту;
3) проверяем другие гипотезы о механизме реакций и выведенных из них моделей, производим дискриминацию других гипотез;
4) находим доверительные интервалы найденных параметров (констант) уравнений, адекватных эксперименту.
Для поиска констант уравнений используем метод наименьших квадратов (МНК). При обработке опытов по МНК требуется большая вычислительная работа. После нахождения оптимальных значений констант необходимо убедиться в адекватности описания эксперимента полученным кинетическим уравнениям (моделям). Для решения этих задач разработали программный комплекс для расчета константы скорости химической реакции и исследования адекватности применяемой модели.
Разработанный программный комплекс позволяет быстро и точно установить порядок реакции, произвести расчет константы скорости химических реакций и исследовать адекватность применяемой модели. Оценка адекватности модели произведена по критерию Фишера, сравнивая опытное значение с табличным.
Данный программный комплекс упростит работу специалистов при проектировании химического производства. За счет этого у сотрудников химического производства появляется больше времени для изучения и анализа возможных вариантов конструкции технологического аппарата.