В статье предложен алгоритм расчета скорости адсорбции кислорода на поверхности кристалла с учетом глобальной стохастической модели. Построенная модель следует из микроскопической и учитывает стадию взаимодействия контактирующих фаз
алгоритм расчета, скорость адсорбции, глобальная стохастическая модель, флуктуации
Для прогресса в исследованиях адсорбции газов на поверхностях твердых тел приходится использовать, наряду с традиционными методами, теоретические исследования. При этом возможность строгих математических расчетов параметров взаимодействия на границе раздела фаз газ-твердое тело весьма ограничена. Подобные ограничения связаны с влиянием на процесс большого числа различных факторов: наличие предсостояний, характер адсорбционной связи, наличие латеральных взаимодействий и скорость поверхностной диффузии. Для неоднородных поверхностей учет влияния всевозможных факторов на скорость процесса возможен в рамках ограничения пространственной области реакционной системы. Еще одним из факторов, способных смещать непредсказуемым образом характеристики процесса, являются спонтанные флуктуации. Анализ воздействия флуктуаций на эволюцию реакционной системы можно оценить на пути сравнения результатов стохастических моделей с детерминистической. Изложение общего подхода в исследовании модельного примера предложено в работах [1-3].
При адсорбции газа на поверхности твердого тела все атомы и молекулы занимают определенные места – адсорбционные центры, которые определяются потенциалом взаимодействия частицы с поверхностью. Геометрия расположения адсорбционных центров определяется рельефом поверхности. Наличие атомарного слоя на поверхности может влиять на перестройку поверхности, но при условии, что заполнение адсорбционной поверхности газом достаточно высокое, то в первом приближении можно не учитывать перестройку поверхностного слоя твердого тела [4, 5]. В этой связи моделирование процесса адсорбции кислорода на поверхности кристалла выполнялось в двумерной пространственной области.
Описание кинетики процессов, происходящих на поверхности, составляет иерархическую последовательность математических моделей [6, 7]:
1) микроскопическая стохастическая модель, учитывающая влияние латеральных взаимодействий между адсорбированными частицами на скорости элементарных стадий и являющаяся наиболее детальной. Данная модель и результаты расчетов представлены в работах [1-3];
2) глобальная стохастическая модель, которая следует из микроскопической при наличии случайного перемешивания частиц на поверхности кристалла. Математическая модель расчета скорости процесса с использованием глобальной стохастической модели представлена в работе [8];
3) глобальная детерминистическая модель, следующая из стохастической модели, которую предстоит построить.
Цель данной работы заключается в разработке алгоритма расчета скорости адсорбции кислорода на неоднородной поверхности (моносульфида железа) с использованием глобальной стохастической модели.
Моделирование динамики взаимодействия с использованием глобальной стохастической модели выполнялось на основе решения численной задачи для системы дифференциальных уравнений:
1) скорость адсорбции 
–
,
2) скорость десорбции 
–
,
3) скорость образования 
–
, (1)
4) скорость образования 
,
, (2)
где 
– энергетический параметр взаимодействия, 
– номер стадии, 
– номер соседства 
, 
– сорт адсорбированной частицы 
; 
для одноузельных процессов, 
для двухузельных; 
– универсальная газовая постоянная; 
– предэкспоненциальный множитель и энергия активации стадии с номером 
[8, 9]. Процесс миграции выпадает из рассмотрения, так как не изменяет общее число частиц сорта на фрагменте решетки размером 
.
Начальное состояние системы на фрагменте поверхности 
либо полностью соответствует незанятой поверхности, либо формируется после случайного «разбрасывания» по узлам фрагмента заданного числа частиц кислорода. Состояние фрагмента кристаллической решетки 
в момент времени 
определяется совокупностью чисел заполнения 
всех его узлов 
. Каждое состояние принимает одно из значений:
(3)
[4].
Численные расчеты поставленной задачи для двумерной области можно получить с помощью динамического метода Монте-Карло [8]. Алгоритм реализации данного метода состоит из пяти этапов и представлен в [8].
Описанный стохастический алгоритм имитирует микроскопические процессы на решетке и учитывает спонтанные флуктуации, существенно влияющие на эволюцию реакционной системы с нелинейной кинетикой.
Для реализации моделирующей программы расчета скорости процесса с использованием глобальной стохастической модели адсобции кислорода на поверхности моносульфида железа построена блок-схема, представленная на рисунке 1. Каждый из блоков, обозначенных на рисунке 1, представляют собой алгоритмы расчета соответствующих скоростей элементарных стадий. Процедура расчета скоростей элементарных актов основана на решении следующих задач:
1) вероятность осуществления соответствующего элементарного акта. Для решения этой задачи используется математическая модель поверхности структуры кристалла, основанная на строении кристаллической решетки моносульфида железа [1] и распределения частиц адсорбата на неоднородной поверхности адсорбента;
2) непосредственный расчет скоростей элементарных актов на основе теории абсолютных скоростей реакций. Константы скоростей отражают процессы, происходящие на атомно-молекулярном уровне, и служат параметрами при построении кинетических уравнений [1]. В блоке 1 выполняется расчет скорости адсорбции 
для элемента 
, изменение параметров 
. В блоке 2 выполняется расчет скорости образования 
для элемента при 
, изменение параметров .
|
Начало |
|
Инициализация
|
|
Ввод данных:
|
|
|
|
|
|
|
|
|
|
да |
|
да |
|
да |
|
да |
|
|
|
|
|
|
|
1 |
|
|
|
да |
|
да |
|
Конец |
|
да |
|
нет |
|
нет |
|
нет |
|
нет |
|
|
|
нет |
|
нет |
|
нет |
|
нет |
|
|
|
|
|
Вывод
|
|
|
|
|
|
2 |
|
3 |
|
|
|
4 |
|
5 |
|
|
|
|
|
|
|
|
Рисунок 1 – Блок-схема алгоритма расчета скорости адсорбции кислорода на поверхности моносульфида железа (глобальная стохастическая модель).
В блоке 3 выполняется расчет скорости образования 
для при 
, изменение параметров . В блоке 5 выполняется расчет скорости образования для элемента при 
, изменение параметров . В блоке 5 выполняется расчет скорости образования 
для элемента при 
, изменение параметров .
Построенная блок-схема являются основной составной частью алгоритма расчета скорости адсорбции кислорода на поверхности кристалла с использованием глобальной стохастической модели [8]. Глобальная стохастическая модель представляет описание в терминах общего числа частиц кислорода, адсорбированных на поверхности, в отличие от микроскопической стохастической модели, описывающей процесс в терминах заполнения всех узлов решетки [10]. Сравнение результатов расчетов по микроскопической стохастической модели и глобальной стохастической модели позволит оценить влияние спонтанных флуктуаций на динамику реакционной системы.
1. Свердлова, О.Л. Автоматизация управления технологическими процессами разделения газов в промышленности: специальность 05.13.06: диссертация на соискание ученой степени кандидата технических наук защищена / Свердлова Ольга Леонидовна; Иркутский государственный университете путей сообщения. - Иркутск, 2014. - 117 с. - Библиогр.: с. 42-63. - Текст: непосредственный
2. Свердлова, О.Л. Математическая модель взаимодействия кислорода с поверхностью моносульфида железа / О.Л. Свердлова, С.В. Иванова, Н.М. Туркина. - Текст: непосредственный // Вестник АГТА. - 2014. - № 8. - С. 118-122.
3. Свердлова, О.Л. Общий подход к описанию нелинейной динамики взаимодействия на границе раздела фаз газ-твердое тело / О.Л. Свердлова, Л.М. Кондратьева. - Текст: непосредственный // Вестник АнГТУ. - 2019. - № 13. - С. 97-100.
4. Свердлова, О.Л. Математическое моделирование затухания автоколебаний адсорбции кислорода на неоднородной поверхности с использованием стохастического подхода / О.Л. Свердлова, Н.Н. Добрынина, Л.Г. Евсевлеева, Н.Н. Туркина. - Текст: не-посредственный // Современные технологии. Системный анализ. Моделирование. - Иркутск. - 2015. - № 4. - С. 88-91.
5. Товбин, Ю.К. Теория абсолютных скоростей реакций на границе раздела фаз газ-твердое тело / Ю.К. Товбин. - Текст: непосредственный // Электрохимия. - 2009. - Т. 45, № 9. - С. 1030-1045.
6. Макеев, А.Г. Автоколебания скорости гетерогенной каталитической реакции: сравнение детерминистического и стохастического подходов к моделированию / А.Г. Макеев, Н.Л. Семендяева. - Текст: непосредственный // Математическое моделирование. - 1996. - Т. 8, № 8. - С. 76-96.
7. Макеев, А.Г. Сравнение стохасти-ческого и детерминистического подходов к моделированию мономолекулярной термодесорбции / А.Г. Макеев, Н.Л. Семендяева. - Текст: непосредственный // Математическое моделирование. - 1995. - Т. 7, № 8. - С. 29-40.
8. Свердлова, О.Л. Глобальная сто-хастическая модель процесса адсорбции кислорода на поверхности кристалла / О.Л. Свердлова, Л.М. Кондратьева, С.В. Иванова. - Текст: непосредственный // Сборник научных трудов ангарского государственного технического университета. - Ангарск. - 2019. - С. 222-227.
9. Глесстон, С. Теория абсолютных скоростей реакций: [пер. с англ.] / С. Глес-стон, К. Лейдлер, Г. Эйринг. - М.: ИИЛ, 1948. - 583 с. - Библиогр.: С. 10-30. - Текст непосредственный
10. Свердлова, О.Л. Стохастическая модель процесса адсорбции кислорода / О.Л. Свердлова, Л.М. Кондратьева, С.В. Иванова. - Текст: непосредственный // Новые информационные технологии в исследовании сложных структур // Материалы XIII Международной конференции. - Томск. - 2020. - С. 106.
