ISSN 0236-235X (P)
ISSN 2311-2735 (E)

Journal influence

Higher Attestation Commission (VAK) - К1 quartile
Russian Science Citation Index (RSCI)

Bookmark

Next issue

4
Publication date:
09 December 2024

Software implementation of the algorithm for finding the optimal temperature condition of the catalytic process

Date of submission article: 29.09.2021
Date after edit article: 08.11.2021
UDC: 544.4:004.4
The article was published in issue no. № 1, 2022 [ pp. 106-112 ]
Abstract:The paper describes the software for solving problems of optimal control of the catalytic process in an ideal mixing reactor. The general form of the formulated optimal control problem is based on a mathe-matical model of the process in an ideal mixing reactor. The refrigerant temperature is considered a control parameter; its values are limited. In order to solve the problem numerically, the paper presents a step-by-step algorithm based on the evolutionary method of artificial immune systems. The method of artificial immune systems for solving optimal control problems makes it possible to obtain an ap-proximate solution to the problem in a time that is acceptable from a practical point of view. The formulated algorithm is a basis for the developed application for the catalytic reaction of α-methylstyrene dimerization in the presence of the NaHY zeolite catalyst, whose products (linear di-mers) are widely used in industrial production. To calculate the process in an ideal mixing reactor, the program implements a number of numerical methods. The software tool allows the user to customize both the catalytic process parameters and the parameters of the algorithm of the artificial immune system method. The user sets the optimization criterion during the program operation, which makes it possible to use it for various formulations of the tasks of the catalytic process of α-methylstyrene dimerization and to obtain a set of optimal concentrations of substances and optimal temperature conditions that correspond to the specified process indicators. The paper gives a solution to the problem of finding the optimal tem-perature condition for the considered catalytic process, the optimality criterion of which is the achievement of the maximum yield of linear dimers with a minimum yield of cyclic dimers and trimers.
Аннотация:В статье описано ПО для решения задач оптимального управления каталитическим процессом в реакторе идеального смешения. На основе математической модели процесса в реакторе идеального смешения сформулирована в общем виде задача оптимального управления. В качестве управляющего параметра рассматривается температура хладоагента, на значения которой наложены ограничения. Приведен пошаговый алгоритм численного решения поставленной задачи, в основу которого положен эволюционный метод искусственных иммунных систем. Применение данного метода дает возможность получить приближенное решение задачи за приемлемое с практической точки зрения время. На основе сформулированного алгоритма разработано приложение для каталитической реакции димеризации a-метилстирола в присутствии цеолитного катализатора NaHY, продукты которой (линейные димеры) широко применяются в промышленном производстве. Для расчета процесса в реакторе идеального смешения в программе реализован ряд численных методов. В программном средстве имеется возможность настройки пользователем параметров как каталитического процесса, так и алгоритма метода искусственных иммунных систем. Критерий оптимизации задается пользователем в ходе работы программы, что позволяет применять ее для различных постановок задач каталитического процесса димеризации α-метилстирола и получать набор оптимальных концентраций веществ и оптимальный температурный режим, соответствующие заданным показателям протекания процесса. Приведено решение задачи поиска оптимального температурного режима рассматриваемого каталитического процесса, критерием оптимальности которой является достижение максимального выхода линейных димеров при минимальном выходе циклических димеров и тримеров.
Authors: E.V. Antipina (stepashinaev@ya.ru) - Sterlitamak Branch of Bashkir State University (Junior Researcher), Sterlitamak, Russia, Ph.D, S.A. Mustafina (mustafina_sa@mail.ru) - Sterlitamak Branch of Bashkir State University (Professor), Sterlitamak, Russia, A.F. Antipin (andrejantipin@ya.ru) - Sterlitamak Branch of Bashkir State University (Associate Professor), Sterlitamak, Russia, Ph.D
Keywords: optimal temperature condition, ideal mixing reactor, linear dimers, artificial immune systems
Page views: 2895
PDF version article

Font size:       Font:

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

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

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

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

В настоящее время метод искусственных иммунных систем используется при решении задач распознавания образов [4], в качестве инструмента интеллектуального анализа бизнес-процессов [5], а также находит практическое применение для настройки промышленных регуляторов [6] и в реализации диагностики электрических систем [7, 8].

Алгоритм поиска оптимального  температурного режима  каталитического процесса

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

Пусть каталитический процесс протекает в реакторе идеального смешения периодического действия. Математическая модель каталитического процесса в реакторе идеального смешения описывается уравнениями материального и теплового баланса:

                     (1)

         (2)

с начальными условиями

                                    (3)

                                              (4)

где x(t) = (x1(t), x2(t), …, xn(t))T – вектор переменных состояния процесса (концентрация веществ, переменный реакционный объем); t – время; Т – температура реакционной смеси;  f1(x(t), T(t), t) – функции, непрерывные вместе со своими частными производными; Cp – мольная теплоемкость реакционной смеси; Qj – теп- ловой эффект j-й стадии реакции; Wj – приведенная скорость j-й стадии реакции; ax – коэффициент теплопередачи; Sx – удельная поверхность теплосъема; C0 – начальная мольная плотность реакционной среды (начальная суммарная концентрация); T – температура реакционной смеси; Tx – температура хладоагента.

Пусть параметром управления является температура хладоагента Tx(t), на значения которой наложены ограничения:

                                        (5)

Допустимые значения температуры реакционной смеси задаются неравенством

Ta £ T(t) £ Tb.                                         (6)

Требуется определить оптимальный температурный режим Tx(t) каталитического процесса, описываемого системой дифференциальных уравнений (1), (3) с начальными условиями (2), (4) и с учетом ограничений (5), (6), доставляющий максимум критерию оптимальности:

G(x(tend), T, Tx) ® max.                           (7)

Критерий оптимальности выражает максимальный выход целевых продуктов реакции (для положительных значений коэффициентов μi) или минимальное содержание побочных продуктов реакции (для отрицательных значений коэффициентов μi) в конечный момент времени.

Сформулируем алгоритм поиска оптимального температурного режима каталитического процесса в реакторе идеального смешения, используя метод искусственных иммунных систем. Работа алгоритма заключается в переходе из поколения в поколение наиболее приспособленных иммунных клеток, представляющих собой набор векторов [9, 10]. Для определения приспособленности иммунной клетки необходимо определить значение фитнес-функции, в качестве которой выступает критерий оптимальности (7).

Шаг 1. Создание начальной популяции температуры хладоагента.

Задаются следующие параметры алгоритма: размер популяции rp, количество иммунных клеток для отбора cp, количество клеток для оператора клонирования cl, параметр оператора мутации mut, количество иммунных клеток с наименьшим значением фитнес-функции cb, количество точек разбиения временного интервала m, общее количество итераций алгоритма iter, счетчик итераций i_iter := 0.

Случайным образом (с учетом неравенства (5)) заполняются значения иммунных клеток начальной популяции температуры хладоагента:

  j = 1, …, rp.

Вычисляются значения фитнес-функции (7) для каждой иммунной клетки  путем решения системы дифференциальных уравнений (1), (2) с начальными условиями (3), (4).

Шаг 2. Клонирование иммунных клеток.

Из популяции i_iter выбираются cp иммунных клеток с наибольшим значением фитнес-функции. Для каждой выбранной клетки создаются cl копий (клонов).

Шаг 3. Мутация клеток-клонов.

С помощью генератора случайных чисел для каждого клона генерируются числа    Если l > 0,5, то   иначе    

Шаг 4. Вычисление значений фитнес-функции.

Для полученных клонов-мутантов вычисляется значение фитнес-функции (7).

Шаг 5. Обновление популяции.

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

Шаг 6. Переход к новой популяции.

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

Шаг 7. Проверка условия окончания поиска.

Если i_iter < iter, то осуществляется переход к шагу 2, иначе поиск решения завершается. Решением задачи поиска оптимального температурного режима каталитического процесса является иммунная клетка из последней популяции, которой соответствует наибольшее значение фитнес-функции.

Описание ПО

Сформулированный алгоритм реализован  в виде программного средства в среде визуального программирования Delphi для каталитической реакции димеризации a-метилстирола в присутствии цеолитного катализатора NaHY.

Механизм реакции димеризации a-метил- стирола в присутствии цеолитного катали- затора NaHY описывается совокупностью стадий [11]:

(8)

где X1 – α-метилстирол; X2 – α-димер; X3 –  β-димер; X4 – циклический димер; X5 – тримеры.

Кинетические уравнения реакции имеют вид:

                   (9)

где Ci – концентрация i-го вещества  (кмоль/м3); kj  – константа j-й стадии  (м3/(кмоль×ч)), которая рассчитывается из уравнения Аррениуса.

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

 

где                               (10)

 

где                                    (11)

с начальными условиями:

                    (12)

где N = C/C0 – относительное изменение числа молей реакционной среды,  – начальная суммарная концентрация веществ (кмоль/м3); vij – матрица стехиометрических коэффициентов веществ ;  Wj = wj/C0 – приведенная скорость j-й стадии реакции (1/ч).

Уравнение теплового баланса рассматриваемого процесса имеет вид

              (13)

с начальным условием

T(0) = T0,                                              (14)

где Cp – мольная теплоемкость реакционной смеси (кДж/(кмоль∙К)); Qj – тепловой эффект  j-й стадии реакции ; Wj – приведенная скорость j-й стадии реакции ; ax – коэффициент теплопередачи (ккал/(м2∙К)); Sx – удельная поверхность теплосъема (1/м); T – температура протекания реакции (К); Tx – температура хладоагента (К).

Значения кинетических и технологических параметров процесса димеризации a-метилстирола приведены в [12].

Работа приложения состоит из двух этапов.

Этап 1. Настройка параметров процесса оптимизации.

На вкладке «Оптимизация процесса» за- даются значения времени протекания каталитического процесса и технологические ограничения температуры реакции и температуры хладоагента. Затем с помощью кнопок необходимо ввести выражение для критерия оптимальности, которое представляет собой математическое выражение вида  где mi = ±1. Знак коэффициента зависит от того, максимальное или минимальное значение компонента xi требуется найти в задаче оптимального управления. На рисунке 1 критерий опти- мальности

x2 + x3 – x4 – x5 ® max                           (15)

соответствует задаче поиска максимума выхода линейных димеров (x2 + x3) при минимальном выходе циклических димеров x4 и тримеров x5 в конечный момент времени протекания реакции. После ввода выражения критерия оптимальности следует нажать кнопку «Сформировать критерий».

Этап 2. Настройка параметров процесса алгоритма.

На вкладке «Параметры алгоритма» устанавливаются параметры алгоритма метода искусственных иммунных систем (рис. 2). Также для системы дифференциальных уравнений (10), (11), (13) необходимо задать начальные условия (12), (14) и выбрать численный метод ее решения. В программе реализованы следующие неявные методы решения системы дифференциальных уравнений: неявный метод Эйлера, метод прогноза и коррекции, неявный метод Рунге–Кутта.

После настройки всех параметров процесса и алгоритма следует нажать кнопку «Расчет». Полученное решение задачи поиска оптимального температурного режима для процесса димеризации a-метилстирола представляется в программе в численном и графическом видах.

На рисунке 2 показано решение задачи поиска оптимального температурного режима Tx(t) для процесса, описываемого системой дифференциальных уравнений (10), (11), (13) с начальными условиями (12), (14) и с учетом ограничений 303 K £ Tx(t) £ 403 K, 303 K £  £ T(t) £ 403 K, с целью достижения максимума критерием оптимальности (15).

В результате выполненных расчетов с помощью программного средства вычислены оптимальные концентрации веществ и оптимальный температурный режим процесса димеризации a-метилстирола в реакторе идеального смешения продолжительностью 3 часа. Установлено, что для достижения максимального выхода линейных димеров при минимальном выходе циклических димеров и тримеров и начальной температуре реакционной смеси, равной 303 К, необходимо:

-     на протяжении 0,5 ч поддерживать температуру хладоагента на уровне 403 К;

-     резко снизить температуру хладоагента до 303 К и в течение 1,5 ч удерживать ее на этом уровне;

-     резко поднять температуру до 403 К и поддерживать ее до конца процесса (рис. 2).

При соблюдении рассчитанного температурного режима хладоагента в конце реакции будут получены максимальный выход линей- ных димеров, равный 0,83 мольных долей,  и минимальный выход циклических димеров и тримеров, равный 0,91 мольных долей.

Заключение

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

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

Исследование выполнено в рамках государственного задания Министерства науки и высшего образования Российской Федерации (код научной темы FZWU-2020-0027).

Литература

1.    Островский Г.М., Зиятдинов Н.Н., Емельянов И.И. Синтез оптимальных систем простых ректификационных колонн с рекуперацией тепла // Докл. Академии наук. 2015. Т. 461. № 2. С. 189–192. DOI: 10.7868/S0869565215080149.

2.    Santos L., Villas-Boas F., Oliveira A.R.L., Perin C. Optimized choice of parameters in interior-point methods for linear programming. Computational Optimization and Applications, 2019, vol. 73, pp. 535–574. DOI: 10.1007/S10589-019-00079-9.

3.    Biegler L.T. Integrated optimization strategies for dynamic process operations. Theoretical Foundations of Chemical Engineering, 2017, vol. 51, no. 6, pp. 910–927. DOI: 10.1134/S004057951706001X.

4.    Astachova I.F., Zolotukhin A.E., Kurklinskaya E.Y., Belyaeva N.V. The application of artificial immune system to solve recognition problems. J. Phys.: Conf. Ser., 2019, vol. 1202, art. 012003. DOI: 10.1088/1742-6596/1202/1/012003.

5.    Wei Y., Wang J. An artificial immune system approach to business process mining. Advanced Materials Research, 2012, vol. 472–475, pp. 35–38. DOI: 10.4028/www.scientific.net/AMR.472-475.35.

6.    Wang M., Feng S., He C., Li Z., Xue Y. An artificial immune system algorithm with social learning and its application in industrial PID controller design. Mathematical Problems in Engineering, 2017, vol. 2017,  pp. 1–13. DOI: 10.1155/2017/3959474.

7.    Lima F.P.A., Lopes M.L.M., Lotufo A.D.P., Minussi C.R. An artificial immune system with continuous-learning for voltage disturbance diagnosis in electrical distribution systems. Expert Systems with Applications, 2016, vol. 56, pp. 131–142. DOI: 10.1016/j.eswa.2016.03.010.

8.    Lima F.P.A., Lotufo A.D.P., Minussi C.R. Wavelet-artificial immune system algorithm applied to voltage disturbance diagnosis in electrical distribution systems. IET Generation, Transmission and Distribution, 2015, no. 9, pp. 1104–1111. DOI: 10.1049/IET-GTD.2014.1102.

9.    Антипина Е.В., Мустафина С.А., Антипин А.Ф. Численный алгоритм идентификации кинетической модели химической реакции // Вестн. технологич. ун-та. 2019. Т. 22. № 9. С. 13–17.

10. Колоколов А.А., Леванова Т.В., Поздняков Ю.С. Алгоритмы искусственной иммунной системы для вариантной задачи размещения телекоммуникационных центров // Изв. Иркутского гос. ун-та. Сер. Математика. 2013. Т. 6. № 1. С. 35–44.

11. Степашина Е.В., Байтимерова А.И., Мустафина С.А. Программный комплекс автоматизации процедуры уточнения механизма химической реакции на основе DRGEP-метода // Башкирский химический журнал. 2011. Т. 18. № 3. С. 112–115.

12. Степашина Е.В. Алгоритмы редукции кинетических схем сложных химических процессов: автореф. дис. Уфа: БашГУ, 2013. 22 c.

References

  1. Ostrovskiy G.M., Ziyatdinov N.N., Emelyanov I.I. Synthesis of optimal systems of simple distillation columns with heat recovery. Proc. of the Academy of Sciences, 2015, vol. 461, no. 2, pp. 189–192. DOI: 10.7868/S0869565215080149 (in Russ.).
  2. Santos L., Villas-Boas F., Oliveira A.R.L., Perin C. Optimized choice of parameters in interior-point methods for linear programming. Computational Optimization and Applications, 2019, vol. 73, pp. 535–574. DOI: 10.1007/S10589-019-00079-9.
  3. Biegler L.T. Integrated optimization strategies for dynamic process operations. Theoretical Foundations of Chemical Engineering, 2017, vol. 51, no. 6, pp. 910–927. DOI: 10.1134/S004057951706001X.
  4. Astachova I.F., Zolotukhin A.E., Kurklinskaya E.Y., Belyaeva N.V. The application of artificial immune system to solve recognition problems. J. Phys.: Conf. Ser., 2019, vol. 1202, art. 012003. DOI: 10.1088/1742-6596/1202/1/012003.
  5. Wei Y., Wang J. An artificial immune system approach to business process mining. Advanced Materials Research, 2012, vol. 472–475, pp. 35–38. DOI: 10.4028/www.scientific.net/AMR.472-475.35.
  6. Wang M., Feng S., He C., Li Z., Xue Y. An artificial immune system algorithm with social learning and its application in industrial PID controller design. Mathematical Problems in Engineering, 2017, vol. 2017, pp. 1–13. DOI: 10.1155/2017/3959474.
  7. Lima F.P.A., Lopes M.L.M., Lotufo A.D.P., Minussi C.R. An artificial immune system with continuous-learning for voltage disturbance diagnosis in electrical distribution systems. Expert Systems with Applications, 2016, vol. 56, pp. 131–142. DOI: 10.1016/j.eswa.2016.03.010.
  8. Lima F.P.A., Lotufo A.D.P., Minussi C.R. Wavelet-artificial immune system algorithm applied to voltage disturbance diagnosis in electrical distribution systems. IET Generation, Transmission and Distribution, 2015, no. 9, pp. 1104–1111. DOI: 10.1049/IET-GTD.2014.1102.
  9. Antipina E.V., Mustafina S.A., Antipin A.F. Numerical algorithm of identifying the kinetic model of a chemical reaction. Bull. of the Technological University, 2019, vol. 22, no. 9, pp. 13–17 (in Russ.).
  10. Kolokolov A.A., Levanova T.V., Pozdnyakov Yu.S. Artificial immune system algorithms for the multivariant problem of the telecommunications centers location. The Bull. of Irkutsk State University. Ser. Math., 2013, vol. 6, no. 1, pp. 35–44 (in Russ.).
  11. Stepashina E.V., Baytimerova A.I., Mustafina S.A. The software package for the computer aided procedure of qualified mechanism of chemical reactions based on DRGEP-method. Bashkir Chemical Journal, 2011, vol. 18, no. 3, pp. 112–115 (in Russ.).
  12. Stepashina E.V. Algorithms for Reducing Kinetic Schemes of Complex Chemical Processes. Ph.D. Thesis. Ufa, 2013, 22 p. (in Russ.).

Permanent link:
http://swsys.ru/index.php?page=article&id=4881&lang=en
Print version
The article was published in issue no. № 1, 2022 [ pp. 106-112 ]

Back to the list of articles