Россия
Работа посвящена теоретическому анализу моделирования процесса адсорбции газа на поверхности твердого тела. В основе математической модели реакционной системы лежит микроскопическая стохастическая модель. Показано влияние латеральных взаимодействий на скорость реакции. Установлена причина затухания усредненной скорости процесса по пространственной области
микроскопическая стохастическая модель, адсорбция, флуктуации, латеральные взаимодействия, скорость адсорбции, имитационное моделирование, глобальная стохастическая модель
При рассмотрении кинетики сложных физико-химических процессов, происходящих на границе раздела фаз, микроскопические стохастические модели являются наиболее полными моделями, описывающими динамику систем с вероятностной природой взаимодействия [1], а применение современных компьютерных технологий позволяет рассчитывать сложную нелинейную динамику таких взаимодействий.
Для более детального описания процесса адсорбции применяется модель решеточного газа с определенным множеством микросостояний. Она состоит из математической модели поверхности кристалла и модели элементарных стадий, заданных кинетической схемой реакции. Элементарные стадии образуют компактный механизм процессов, происходящих на поверхности, и устанавливаются в результате математического моделирования линейно независимых реакций [2, 3]. При этом необходимо учитывать, что сокращенная схема ранжирования элементарных стадий формируется при определенных условиях протекания реакции.
Для теоретического анализа гетерофазных реакций микроскопическая стохастическая модель наиболее детально учитывает влияние спонтанных флуктуаций на динамику процесса. Если поверхность кристалла неоднородна, то полноценный учет всех факторов возможен в рамках ограничения пространственной области реакционной системы [4].
Эволюция протекания реакции на границе раздела фаз рассматривается как случайный марковский процесс с непрерывным временем и дискретным множеством состояний для потока элементарных событий. Отдельные реализации рассматриваемого случайного процесса можно получить с помощью метода Монте-Карло [4]. Цель работы заключается в теоретическом анализе процесса адсорбции газа на неоднородной поверхности кристалла, проведенным на основе микроскопической стохастической модели. Реализация модели получена с помощью имитационного моделирования.
Адсорбция газа на поверхности твердого тела вызывается наличием адсорбционного силового поля, создаваемого за счет нескомпенсированности межмолекулярных сил в поверхностном слое. При физической адсорбции основной вклад в энергию взаимодействия вносят дисперсионные силы, которые не изменяют индивидуальных свойств молекул и носят обратимый характер. Например, при понижении давления или концентрации адсорбированного вещества [5]. Поглощение вещества может быть обусловлено образованием химических связей между молекулами адсорбата и поверхностным слоем адсорбента, что сопровождается перестройкой поверхностного слоя адсорбента. Химическая адсорбция необратима, но для некоторых технологических процессов данное обстоятельство имеет решающее значение.
При адсорбции атомы и молекулы газа занимают отдельные адсорбционные центры на поверхности твердого тела. Потенциал взаимодействия определяется геометрией поверхности и латеральными (межмолекулярными) взаимодействиями между адсорбированными частицами (адчастицы, адатомы), которые определяют взаимную ориентацию молекул в пространстве. Механизмы латерального взаимодействия зависят от физико-химических свойств адсорбционных систем и носят прямой и косвенный характер. Прямое взаимодействие возможно при отсутствии поверхности и обусловлено сохранением электронного состояния адсорбированных частиц. Поэтому состояние адсорбированного слоя может быть далеко от термодинамического равновесия и анализ воздействия вышеперечисленных факторов на реакционную систему позволяет оптимизировать проведение подобного рода реакций [6].
В работах [7, 8, 9] рассматривается пример построенной микроскопической стохастической модели, учитывающей стадию взаимодействия контактирующих фаз. Расчет эволюции реакционной системы был осуществлен методом Монте-Карло. В общей сложности было проведено около 200 испытаний. При проведении вычислительных экспериментов менялись размеры фрагмента и концентрация адсорбата. Значения внешних параметров модели 
и универсальная газовая постоянная оставались неизменными. Влияние латеральных взаимодействий было учтено с помощью внутренних параметров модели – энергетических параметров взаимодействия где
– номер стадии
– номер соседства
– сорт адсорбированной частицы
.
На рисунках 1 и 2 представлено графическое изображения колебания скорости процесса адсорбции для различной концентрации адсорбированного газа при изменении размеров фрагмента.
|
V – средняя скорость процесса |
|
Концентрация ʘ |
Рисунок 1 – Колебания скорости процесса адсорбции при изменении концентрации адсорбированного газа: 10%, 15%, 25%, 35%, 50%
|
V – средняя скорость процесса |
|
Концентрация ʘ |
Рисунок 2 – Колебания скорости процесса адсорбции при изменении концентрации адсорбированного газа: 65%, 75%, 80%, 90%, 100%
Проведенные имитационные исследования позволяют сделать качественный анализ результатов моделирования, а непосредственные численные расчеты показывают, что с увеличением размера поверхности и концентрации адсорбата степень влияния латеральных взаимодействий между адсорбированными частицами нивелируется. Наибольшее влияние латеральных взаимодействии на колебания скорости процесса достигается при 50% концентрации адсорбата на поверхности твердого тела.
По результатам исследования была изучена зависимость затухания средней скорости процесса при изменении концентрации адсорбированного газа. С этой целью была решена серия задач на фрагментах, содержащих от 1000 до 10000 узлов. Из полученных экспериментальных данных видно, что данная зависимость носит обратный линейный характер. Результаты вычислительных экспериментов представлены на рисунке 3. Уравнение регрессии имеет вид
(1)
где – средняя скорость процесса,
– концентрация адсорбата. Мерой качества описания процесса регрессионной моделью является коэффициент детерминации, значение которого
говорит о том, что изменение скорости процесса адсорбции полностью определяется изменением концентрации адсорбированного газа. Также этот коэффициент для линейной регрессии определяет высокую точность аппроксимации. На рисунке 3 представлены экспериментальные значения изменения скорости процесса адсорбции при изменении концентрации газа и график уравнения регрессии.
|
V – средняя скорость процесса |
|
Концентрация ʘ |
Рисунок 3 – График уравнения регрессии и экспериментальные точки
В результате проведенных имитационных исследований произведена оценка влияния латеральных взаимодействий на усредненную скорость процесса. Установлено, что флуктуации скорости реакции достигают своего максимума при .
Установлено, что зависимость усредненной скорости процесса от концентрации адсорбированного газа близка к линейной и монотонно убывает в первом приближении. При возрастании числа занятых адсорбционных центров увеличивается вероятность осуществления стадии хемосорбции и латеральные взаимодействий не оказывают существенного влияния на скорость процесса.
1. Куркина, Е.С. Математическое моделирование пространственно-временных структур, возникающих в гетерогенной каталитической системе / Е.С. Куркина, Н. Л. Семендяева – Текст: непосредственный // Прикладная математика и информатика: Труды факультета ВМК МГУ. М.: МАКС Пресс, 2011. – № 37. С.14−43.
2. Жданов, В.П. Модель решеточного газа для описания хемосорбции на поверхностях металлов / В.П. Жданов, К.И. Замараев // Успехи физических наук. – 1986. – Т.149. №4. – С. 636-669.
3. Кафаров, В.В. Математическое моделирование основных процессов химических производств: Учеб. пособие для вузов / В.В. Кафаров, М.Б Глебов. – М: Высшая школа, 1991. – 367 с. – ISBN: 5-06-002066-5. – Текст: непосредственный.
4. Макеев, А.Г. Автоколебания скорости гетерогенной каталитической реакции: сравнение детерминистического и стохастического подходов к моделированию / А.Г. Макеев, Н.Л. Семендяева. – Текст: непосредственный // Математическое моделирование. – 1996. – Т. 8, № 8. – С. 76-96.
5. Глесстон, С. Теория абсолютных скоростей реакций: [пер. с англ.] / С. Глес-стон, К. Лейдлер, Г. Эйринг. – М.: ИИЛ, 1948. – 583 с. – Библиогр.: С.10-30. – Текст непосредственный.
6. Товбин, Ю.К. Теория абсолютных скоростей реакций на границе раздела фаз газ-твердое тело / Ю.К. Товбин. – Текст: непосредственный // Электрохимия. – 2009. – Т. 45, № 9. – С. 1030-1045.
7. Свердлова, О.Л. Автоматизация управления технологическими процессами разделения газов в промышленности: специальность 05.13.06: диссертация на соискание ученой степени кандидата техниче-ских наук защищена / Свердлова Ольга Леонидовна; Иркутский государственный университете путей сообщения. – Иркутск, 2014. – 117 с. – Библиогр.: с. 42-63. – Текст: непосредственный.
8. Свердлова, О.Л. Математическая модель взаимодействия кислорода с по-верхностью моносульфида железа / О.Л. Свердлова, С.В. Иванова, Н.М. Туркина. – Текст: непосредственный // Вестник АГТА. – 2014. – № 8. – С. 118-122.
9. Свердлова, О.Л. Математическое моделирование затухания автоколебаний адсорбции кислорода на неоднородной поверхности с использованием стохастического подхода / О.Л. Свердлова, Н.Н. Добрынина, Л.Г. Евсевлеева, Н.Н.Туркина. – Текст: непосредственный // Современные технологии. Системный анализ. Модели-рование Иркутск. – 2015. – № 4. – С. 88-91.
