ON DIFFUSION IN LATTICE GAS FOR LARGE CONCENTRATION GRADIENTS

Research article
DOI:
https://doi.org/10.23670/IRJ.2019.88.10.003
Issue: № 10 (88), 2019
Published:
2019/10/18
PDF

О ДИФФУЗИИ В РЕШЕТОЧНОМ ГАЗЕ ПРИ БОЛЬШИХ ГРАДИЕНТАХ КОНЦЕНТРАЦИЙ

Научная статья

Трощиев Ю.В. *

ORCID: 0000-0001-6805-7001,

Московский Государствнный Университет, Москва, Россия

* Корреспондирующий автор (yuvt[at]yandex.ru)

Аннотация

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

Ключевые слова: решеточный газ, коэффициенты диффузии, уравнения химической кинетики. 

ON DIFFUSION IN LATTICE GAS FOR LARGE CONCENTRATION GRADIENTS

Research article

Troshchiev Yu.V. *

ORCID: 0000-0001-6805-7001,

Moscow State University, Moscow, Russia

* Corresponding author (yuvt[at]yandex.ru)

Abstract

The article derives expressions for flows and diffusion coefficients in an imperfect lattice gas from equations of chemical kinetics convenient for study. Coefficients are linear operators depending on concentrations and their gradients. It is shown that the flow direction can build any angle with the direction of the gradient.

Keywords: lattice gas, diffusion coefficients, chemical kinetics equations. 

Введение

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

Наиболее точным способом вычисления более сложных коэффициентов диффузии является использование кинетических уравнений (при условии достаточной точности коэффициентов в самих кинетических уравнениях). Примерами применения такой методики являются вывод коэффициентов диффузии в [2] и кинетически согласованные разностные схемы [4].

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

Однако градиенты концентраций не всегда можно считать малыми. В [8] автором выведены выражения для потока мигрирующих частиц и для коэффициентов диффузии без предположения о малости градиентов. За основу взяты уравнения химической кинетики [5], [6], [7]. В [8] учитывается также несовпадение направлений градиентов и потоков, т.е. полученные коэффициенты диффузии являются матричными и зависят от концентраций адсорбированных веществ и от их градиентов. В [9], [10], [11] проводилось исследование в приближении среднего поля и квазихимическом (см. [5], [6]). Исследуем полученные коэффициенты в этих приближениях более подробно.

Диффузия в решеточном газе

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

Состояние решеточного газа задается числами заполнения ni, где i - индекс узла решетки. Значение ni равно номеру химической компоненты, которой занят узел (0 - для свободного узла). Числа заполнения изменяются за счет происходящих по определенным статистическим законам прыжков частиц в соседние узлы, обменов с газовой фазой и реакций между частицами, находящимися в соседних узлах.

Описание решеточного газа как сплошной среды осуществляется с помощью уравнений реакции-диффузии [7]:

29-10-2019 14-37-22     (1)

где N - число компонент, 29-10-2019 14-37-46 - плотность n-й компоненты, функция 29-10-2019 14-38-11 описывает обмен с газовой фазой и химические реакции. Поток 29-10-2019 14-38-19 определяется в соответствии с законом Фика:

29-10-2019 14-38-27

где 29-10-2019 14-38-38 - матрица диффузии.

Вывод выражений для потока и коэффициентов диффузии

При описании состояния решеточного газа на уровне уравнений химической кинетики используются следующие переменные [5], [6], [12]: 29-10-2019 14-48-23 - вероятность заполнения i-го узла n-ой компонентой, 29-10-2019 14-48-34 - вероятность одновременного заполнения i-го узла n-ой компонентой, а j-го узла m-ой компонентой для всех пар соседних узлов. Могут рассматриваться и более чем парные вероятности, но в дальнейшем это не учитывается, так как не вносит качественных изменений в схему вывода коэффициентов диффузии. В общем виде уравнения химической кинетики для решеточного газа являются обыкновенными дифференциальными уравнениями и имеют вид [7]

29-10-2019 14-50-04   (3)

где индексы 29-10-2019 14-50-59 пробегают все значения, для которых введены соответствующие переменные. Хотя формально в правые части входят полные векторы  29-10-2019 14-51-09, реально в число аргументов функций  29-10-2019 14-51-21 входят лишь переменные, относящиеся к узлам решетки, находящимся вблизи узлов  29-10-2019 14-51-28.

Можно предположить, что локальные парные вероятности 29-10-2019 14-51-35 устанавливаются существенно быстрее, чем происходит процесс переноса массы по поверхности. Поэтому для вывода коэффициентов диффузии можно вместо соответствующего дифференциального уравнения рассматривать уравнение

29-10-2019 14-54-40   (4)

Будем считать, что уравнение (4) однозначно разрешимо относительно переменных g. Для большинства приближений это так, либо переменные g явно выражаются через переменные 29-10-2019 14-56-26. Поэтому определена функция 29-10-2019 14-56-34. Подстановка этой функции в первое уравнение системы (3) приводит к уравнению относительно плотностей заполнения:

29-10-2019 14-59-51   (5)

Функция 29-10-2019 15-00-24 имеет вид 29-10-2019 15-00-11    (6)

где 29-10-2019 15-03-01  - множество ближайших соседей узла i, 29-10-2019 15-03-10 - скорости миграции n-ой компоненты из узла i в узел j и из узла j в узел i, соответственно. Функция 29-10-2019 15-03-21 описывает обмен с газовой фазой и химические реакции; она также представляет собой сумму скоростей элементарных процессов.

Пусть в узле i заданы концентрации 29-10-2019 15-03-32 и их градиенты 29-10-2019 15-03-41. Тогда с помощью линейной экстраполяции могут быть найдены концентрации в любом узле j из окрестности этого узла:

29-10-2019 15-10-37   (7)

где 29-10-2019 15-11-08 - вектор из узла i в узел j; (,) – скалярное произведение. Подстановка (7) в (6) позволяет выразить29-10-2019 15-11-16 через  29-10-2019 15-11-28. Выражение для потока может быть получено суммированием всех потоков в узле i:

29-10-2019 15-11-50    (8)

Отсюда коэффициенты диффузии могут быть вычислены как производные от потока по градиенту, взятые с обратным знаком:

29-10-2019 15-15-59    (9)

Выражения (8) и (9) имеют отличия от классического случая (1), (2). Поток (8), в отличие от (2), зависит от градиента нелинейно и, вообще говоря, должен быть непосредственно подставлен в выражение (1). Для того чтобы выяснить качественные отличия между полученным выражением для потока (8) и выражением (2), удобнее рассматривать коэффициенты диффузии (9). В силу (8), выражения (9) зависят не только от концентраций, но и от их градиентов. Более того, коэффициенты диффузии являются линейными операторами, т.е. направление потока может отличаться от направления градиента концентрации даже в случае однокомпонентного газа.

Зависимость коэффициентов диффузии от величины градиента позволяет сделать вывод о том, существенна ли подстановка в уравнение (1) выражения (8), а не произведения коэффициентов диффузии (9) на градиенты концентраций. Существенных отличий от вариантов диффузии, рассмотренных в [6], [7], естественно ожидать при больших градиентах и сильных латеральных взаимодействиях.

Будем считать, что ось абсцисс параллельна линиям решетки. Направление потока, противоположное градиенту, назовем нормальным. Угол между направлением потока и положительным направлением оси абсцисс обозначим β, а угол между нормальным направлением и положительным направлением оси абсцисс - 29-10-2019 15-17-08. Угол между положительным направлением оси абсцисс и направлением градиента, соответственно, равен 29-10-2019 15-17-15. Угол между направлением градиента и направлением потока 29-10-2019 15-17-24.

Будем рассматривать поток не в узле кристаллической решетки, а в центре ячейки. При этом шаблон, который необходимо учитывать, состоит не из тринадцати точек, а из двенадцати, и он более компактный, что существенно при изучении больших градиентов концентраций (рис. 1). Концентрацию в центре шаблона обозначим A0. Четыре узла решетки, находящиеся в углах центральной ячейки этого шаблона, будем называть внутренними. Остальные восемь узлов – внешними. Будем считать, что в центре шаблона находится начало координат.

29-10-2019 15-20-48

а)                                           б)

Рис. 1 – а) Шаблон для потока в узле кристаллической решетки; б) Шаблон для потока в центре ячейки кристаллической решетки

 

Будем считать, что 29-10-2019 15-22-23, а величина градиента находится в пределах 29-10-2019 15-22-33. В этом случае концентрации во внутренних узлах шаблона изменяются в интервале (0;1). В случае большого градиента формула (7) может выдавать выходящие за пределы физически допустимых значения концентраций во внешних узлах (рис. 2). Вероятности, выходящие за пределы интервала (0;1), можно, например, заменять на  или , но такое приближение не очень точное и не везде дифференцируемое. Другой вариант – заменять значения на максимальные и минимальные концентрации во внутренних точках шаблона. Реальные значения концентраций обычно находятся между этими двумя приближениями. Можно также построить гладкую функцию со значениями, заключенными в этом промежутке. В данной работе используется следующая функция:

29-10-2019 15-25-10   (10)

где i  − внутренний узел шаблона с максимальной концентрацией.   29-10-2019 15-26-25

Рис. 2 – Варианты вычисления концентрации во внешних узлах

 

Совпадение направления градиента и линий решетки

Будем рассматривать квадратную решетку, узлы которой некоторым образом пронумерованы. Пусть слой является однокомпонентным (компонента A), и направление градиента концентрации совпадает с направлением линий решетки. В этом случае поток в однокомпонентном слое из узла i в узел j в приближении среднего поля может быть записан следующим образом:

29-10-2019 15-28-11   (11)

где 29-10-2019 15-29-32 - вероятность заполнения узла i; 29-10-2019 15-29-38 - вектор, соединяющий узел i с узлом j; 29-10-2019 15-34-13 - множитель, учитывающий предэкспоненту и энергию активации; ε - параметр латерального взаимодействия; T - температура. Встречный поток из узла j в узел i равен

29-10-2019 15-34-24     (12)

Суммарный поток между узлами i и j равен разности 29-10-2019 15-34-34   (13) Если пренебречь производными 29-10-2019 15-37-54, получается точное соответствие с классической диффузией, и коэффициент равен

29-10-2019 15-35-14    (14)

Множитель  30-10-2019 10-49-31, которым мы пренебрегли, пропорционален разности концентраций 30-10-2019 10-49-39. Таким образом, пренебрегать этим множителем можно только в том случае, если

  30-10-2019 10-49-48    (15) Ясно, что это условие не выполняется, если плотность заполнения заметно меняется на нескольких (30-10-2019 10-51-20) шагах решетки (это соответствует границам между различными фазовыми состояниями.).

Условие (15) может также не выполняться, если 30-10-2019 10-51-28 достаточно мало, но испытывает достаточно большие относительные изменения от узла к узлу. В этом случае, если считать 30-10-2019 10-51-40, то значения экспонент в выражениях (11) и (12) можно считать приблизительно равными единице, а коэффициенты диффузии - равными 30-10-2019 10-51-59. Поэтому условие применимости выражения (14) принимает вид:

  30-10-2019 10-52-10     (16)

На рис. 3 приведена зависимость коэффициента диффузии от величины градиента для различных значений параметров. Излом в точке 30-10-2019 11-00-13 соответствует достижению в крайних точках шаблона значений концентрации 0 и 1. Этот излом особенно заметен при положительном параметре латерального взаимодействия.

30-10-2019 11-04-00

Рис. 3 – Зависимость коэффициента диффузии от величины градиента при совпадении направления градиента и направления линий решетки. 30-10-2019 11-05-58

Сплошная линия – 1-й способ вычисления концентраций во внешних узлах, пунктир – второй способ, длинный пунктир – третий способ

Излом в точке 30-10-2019 11-00-13 происходит в сторону графика, построенного по второму способу вычисления концентраций во внешних узлах шаблона. При 30-10-2019 11-12-01  этот способ дает меньшее значение коэффициента диффузии (а,б), а при 30-10-2019 11-12-16 – большее (в-е). Отметим также, что при стремлении градиента к нулю значения коэффициента диффузии для рассматриваемых способов вычисления концентраций во внешних узлах шаблона могут значительно различаться.

Коэффициент диффузии может считаться приблизительно постоянным лишь при значениях градиента <0.1. При значениях >0.1, а особенно при значениях ~1, зависимость коэффициента диффузии от величины градиента очень существенна. Коэффициент диффузии может, как возрастать, так и убывать с ростом величины градиента, однако характерно наличие участка заметного роста коэффициента, означающего что при очень большом градиенте система стремится градиент уменьшить. Это можно видеть из рисунков (г-е), где коэффициент переходит из отрицательной области в положительную. Отрицательное значение коэффициента соответствует тому, что решеточный газ будет собираться в капли. Но если граница капель станет слишком узкой, то она будет “размазываться” диффузией.

Отклонение направления градиента от линий решетки

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

30-10-2019 11-18-08    (17)

где 30-10-2019 11-18-36 - номера узлов кристаллической решетки, находящихся в углах рассматриваемой ячейки. Выражения 30-10-2019 11-18-46 принимают вид:

30-10-2019 11-25-08   (18)

где 30-10-2019 11-28-26  - градиент концентрации в узле i. Коэффициент диффузии в этом случае является матричным:

 30-10-2019 11-26-59    (19)

На рис. 4 приведены графики зависимости величины γ от угла α для различных значений параметров. Они являются нелинейными и немонотонными. Присутствуют зависимости двух типов. На рис. 3 (а-в) направление потока при 30-10-2019 11-27-20 совпадает с направлением градиента (коэффициент диффузии, соответственно, отрицателен). В промежутке между этими значениями α наблюдается отличие угла β от угла 30-10-2019 11-27-28. Можно видеть, также, что результаты, получаемые третьим способом вычисления концентраций во внешних узлах могут выходить за пределы, задаваемые первыми двумя способами.

30-10-2019 11-34-14

Рис. 4 – Угол между направлением потока и нормальным направлением в зависимости от угла между нормальным направлением и направлением линий решетки. 30-10-2019 11-36-34

На рис. 4 (г) направление потока изменяется от совпадающего с направлением градиента при 30-10-2019 11-58-39 до нормального при 30-10-2019 11-58-46. Из этого, в частности, следует возможность ортогонального направления потока по отношению к направлению градиента.

Обсудим это явление более подробно. Заметим, что фактически поток (17) является суммой двух ортогональных потоков. Далее, из формулы (18) следует, что на величину потока в данном направлении влияет лишь проекция градиента на это направление. Таким образом, поток расщепляется на два одномерных потока, каждый из которых подчиняется закону (11)-(13), где производную 30-10-2019 11-58-55 надо брать в направлении соответствующей оси координат.

Это позволяет объяснить возможность ортогонального направления потока по отношению к градиенту. Действительно, как следует из рис. 3, знак коэффициента диффузии может изменяться в зависимости от градиента. При вращении постоянного по величине вектора градиента от 30-10-2019 11-58-39 до 30-10-2019 11-58-46 величина проекции градиента на ось ординат изменяется от 30-10-2019 11-59-33, а величина проекции градиента на ось абсцисс изменяется от 30-10-2019 11-59-44. Если на этих интервалах коэффициент диффузии соответствующим образом меняет знак, то становятся возможными любые углы между направлениями потока и градиента.

В более точных приближениях такое расщепление, по-видимому, невозможно из-за наличия нелинейных связей между корреляционными функциями. Однако это не исключает возможности любых углов между направлениями потока и градиента. Действительно, если, например, при малых значениях градиента коэффициент диффузии отрицателен, то при больших значениях градиента он может стать положительным, так как система будет стремиться уменьшить слишком большой градиент. Таким образом, при некотором значении градиента коэффициент диффузии будет равен нулю. Причем эти значения могут различаться для 30-10-2019 11-58-39  и 30-10-2019 11-58-46. Выбрав значение градиента так, чтобы коэффициент обращался в нуль в интервале 30-10-2019 11-59-59, можно также получить любые углы между направлениями потока и градиента. На рис. 5 показан пример такой ситуации для одного из вариантов учета парных вероятностей (квазихимические вероятности [6]). В приведенном примере изменение угла 30-10-2019 12-00-20 произошло при значении градиента 0.325. И обусловлено это явление зависимостью коэффициента диффузии от величины и направления градиента.

  30-10-2019 12-17-20

Рис. 5 – Угол между направлением потока и нормальным направлением в зависимости от угла между нормальным направлением и линиями решетки.  30-10-2019 12-27-32 третий способ вычисления концентраций

Ортогональность потока и градиента может, например, приводить к образованию капель поверхностной жидкости различной формы. Если при угле 30-10-2019 12-27-52 поток и градиент ортогональны, то они ортогональны также при 30-10-2019 12-28-01. Таким образом, если на поверхности имеется круглая капля, то на ее границе есть восемь точек, в которых вещество диффундирует вдоль границы. И, следовательно, капля будет деформироваться.

Заключение

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

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

Конфликт интересов Не указан. Conflict of Interest None declared.

Список литературы / References

  1. Ландау Л.Д. Статистическая физика / Ландау Л.Д., Лифшиц Е.М. – М: Наука. – 1964. – 567 c.
  2. Левич В.Г. Курс теоретической физики. т. 1 / Левич В.Г. – М: Физматгиз. – 1962. – 695 c.
  3. Хилл Т. Статистическая механика / Хилл Т. – М: ИЛ. – 1960. – 486 c.
  4. Елизарова Т.Г. Кинетические алгоритмы для расчета гидродинамических течений / Елизарова Т.Г., Четверушкин Б.Н. // ЖВМ и МФ. – 1985. – Т. 25. – N – С. 1526.
  5. Товбин Ю.К. Теория физико-химических процессов на границе газ-твердое тело / Товбин Ю.К. – М.: Наука. – 1988. – 285 c.
  6. Жданов В.П. Элементарные физико-химические процессы на поверхности / Жданов В.П.– Новосибирск: “Наука” Сиб. отд. – 1988. – 319 c.
  7. Еленин Г.Г. Математическое моделирование гетерогенных каталитических реакций на гранях монокристаллов благородных металлов. Част I. Сверхструктуры и фазовые переходы / Еленин Г.Г. // Российский химический журнал. – 1996. – N  – С. 19-50.
  8. Трощиев Ю.В. О коэффициентах диффузии в неидеальном решеточном газе при больших градиентах концентраций / Трощиев Ю.В. – Труды факультета ВМК, Прикладная математика и информатика. – М: Диалог-МГУ. – 1999. – N – С. 51-57.
  9. ТрощиевЮ.В. Исследование матричных коэффициентов диффузии в решеточном газе при больших градиентах концентрации в приближении среднего поля / ТрощиевЮ.В.//  Препринт. – М: МАКС Пресс. – 2003. – 10 c.
  10. ТрощиевЮ.В. Матричные коэффициенты диффузии на квадратной решетке / ТрощиевЮ.В. // Ломоносовские чтения. – Москва, МГУ, ф-т ВМК. – Апрель 2006.
  11. Трощиев Ю.В. Матричные коэффициенты диффузии в решеточном газе при больших градиентах концентраций (квазихимическое приближение) / Трощиев Ю.В.// Тихоновские чтения, Секция Математическое моделирование и вычислительные методы. – Москва, МГУ, ф-т ВМК. – 29 октября – 2 ноября 2007.
  12. Еленин Г.Г. Математическое моделирование процесса образования островковых структур на поверхности монокристалла / Еленин Г.Г., Макеев А.Г. // Математическое моделирование. –1991. – N – С. 29.

Список литературы на английском языке / References in English

  1. Landau L.D. Statisticheskaya fizika [Statistical Physics] / L.D. Landau, E.M. Lifshitz. // Butterworth-Heinemann. – 1980. – V. 5. ISBN 978-0-7506-3372-7.
  2. Levich B.G. Kurs teoreticheskoj fiziki [Theoretical Physics]. V. 1 / B.G. Levich.– Мoscow: PhysMathGIz. – 1962. – 695 p. [in Russian]
  3. Hill T.L. Statisticheskaya mekhanika [Statistical mechanics]. / T.L. Hill. – New York a. o., McHaw-Hill. – 1956.
  4. Elizarova T.G. Kineticheskie algoritmy dlya rascheta gidrodinamicheskih techenij [Kinetic algorithms for calculating gas dynamic flows]. / T.G. Elizarova, B.N. Chetverushkin. // USSR Computational Mathematics and Mathematical Physics. – V. 25. – Issue 5. – 1985. – P. 164–169. [in Russian]
  5. Tovbin Yu K. Teoriya fiziko-himicheskih processov na granice gaz-tverdoe telo [Theory of Physical Chemistry Processes at a Gas-Solid Interface] (Advances in Science and Technology in the USSR) / Yu K. Tovbin – Crc Pr I Llc. – 1991. – 349 p. [in Russian]
  6. Zhdanov V.P. Elementarnye fiziko-himicheskie processy na poverhnosti [Elementary Physicochemical Processes on Solid Surfaces] / V.P. Zhdanov. – Springer US. – 1991. DOI: 10.1007/978-1-4899-2373-8. [in Russian]
  7. Elenin G.G. Matematicheskoe modelirovanie geterogennyh kataliticheskih reakcij na granyah monokristallov blagorodnyh metallov. CHast I. Sverhstruktury i fazovye perekhody [Mathematical modeling of heterogeneous catalytic reactions on the faces of crystals of precious metals. Part I. Superstructures and phase transitions] / G.G. Elenin. // Russian Chemical Journal. – 1996. – N 2. – P. 19-50. [in Russian]
  8. Troshchiev Yu.V. O koefficientah diffuzii v neideal'nom reshetochnom gaze pri bol'shih gradientah koncentracij [Diffusion coefficients in a nonideal lattice gas with high concentration gradients]. – Computational Mathematics and Modelling / Yu.V. Troshchiev. – 2000. – V. 11. – N 1. – P. 40-45. DOI: 10.1007/BF02359062 [in Russian]
  9. Troshchiev Yu.V. Issledovanie matrichnyh koefficientov diffuzii v reshetochnom gaze pri bol'shih gradientah koncentracii v priblizhenii srednego polya [Investigation of Matrix Diffusion Coefficients in Lattice Gas under the Large Concentration Gradients in Mean Field Approximation]. – Preprint / Yu.V. Troshchiev. – Moscow: MAKS Press. – 2003. – 10 p. [in Russian]
  10. Troshchiev Yu.V. Matrichnye koefficienty diffuzii na kvadratnoj reshetke [The Matrix Diffusion Coefficients at Square Lattice] / Yu.V. Troshchiev. // Sci. Conf. Lomonosov Readings. – Moscow, MSU, Dept. CMC. – April 2006. [in Russian]
  11. Troshchiev Yu.V. Matrichnye koefficienty diffuzii v reshetochnom gaze pri bol'shih gradientah koncentracij (kvazihimicheskoe priblizhenie) [Matrix Diffusion Coefficients in Lattice Gas under Large Concentration Gradients (Quasi-chemical Approximation)] / Yu.V. Troshchiev // Tihonovskie chteniya, Sekciya Matematicheskoe modelirovanie i vychislitel'nye metody [Sci. Conf. Tikhonov Readings, Section Mathematical Modeling and Numerical Methods] – Moscow, MSU, Dept. CMC. – October 29– November 2 2007. [in Russian]
  12. Elenin G.G. Matematicheskoe modelirovanie processa obrazovaniya ostrovkovyh struktur na poverhnosti monokristalla [Mathematical Modeling of the Process of Island Structure Formation on the Single Cristal Surface] / G.G. Elenin, A.G. Makeev. – Matematicheskoe modelirovanie [Mathematical Modeling]. – 1991. – N 7. – p. 29. [in Russian]