Прямое численное моделирование некоторых физико-химических процессов и явлений В.Л. Ковалев, А.Н. Якунчиков Московский государственный университет им. М.В. Ломоносова

Московский государственный университет им. М.В.Ломоносова, механико-математический факультет ПРИМЕНЕНИЕ МЕТОДА МОНТЕ –КАРЛО ДЛЯ ИССЛЕДОВАНИЯ ГЕТЕРОГЕННОЙ РЕКОМБИНАЦИИ НА ТЕПЛОЗАЩИТНЫХ ПОКРЫТИЯХ МНОГОРАЗОВЫХ АППАРАТОВ В.Л. Ковалев, В.Ю. Сазонова, А.Н. Якунчиков

Введение

Химические реакции Моделируется течение реагирующей бинарной смеси вблизи каталитической поверхности Физическая адсорбция-десорбция Химическая адсорбция Диффузия физадсорбированных атомов Рекомбинация Eley-Rideal Рекомбинация Langmuir-Hinshelwood

Феноменологическая модель - доля заполнения активных мест физически адсорбированными атомами - доля заполнения активных мест химически адсорбированными атомами - константы скоростей реакций Решение (если пренебречь диффузией и рекомбинацией LH):

Модель метода Монте-Карло

Схема метода Монте-Карло данет Случайный выбор ячейки Рекомбинация ER ДесорбцияДиффузия Рекомбинация LH Адсорбция Обратима я? Ячейка свободна? Обратима я? Рек-ция ER ? Десорбция Диффузия Хим. адсор б. Физ. адсор б. да Рек-ция LH нет

Результаты График 1. Доля занятых мест для физической адсорбции от времени для одной компоненты (N). График 2. Доля занятых мест для химической адсорбции от времени для одной компоненты (N). Феноменологическая модель Метод Монте-Карло Феноменологическая модель Метод Монте-Карло

Результаты, Kim Y.C., Boudart M;, Marshall, T. C. J. Chem. Phys. 1962, 37, 2501;, Метод прямого численного моделирования Зависимость коэффициента рекомбинации азота от температуры Доля заполнения поверхности физически адсорбированными атомами азота и кислорода в зависимости от концентрации в газовой фазе. азот кислород

Исследование течения и теплообмена в микро и нано каналах методами молекулярной динамики. Москва 2007 Московский государственный университет им. М.В.Ломоносова, механико-математический факультет, кафедра газовой и волновой динамики. Ковалёв В. Л., Якунчиков А. Н.

Введение Тенденция в развитии электроники, которая проявилась а прошлое десятилетия, состоит в том, что количество энергии, рассеиваемое системами охлаждения электронных компонентов, неуклонно увеличивается. В будущем микроэлектронные компоненты будут только уменьшаться в размерах, поэтому вопрос об их охлаждении стоит достаточно остро. Предполагается, что системы охлаждения будут представлять из себя систему микро (или даже нано) каналов, пронизывающую электронный компонент. По этим каналам будет осуществляться циркуляция охлаждающей жидкости или газа (также возможны потоки с двумя фазами). Прототипы таких устройств уже появляются в исследовательских институтах США (Purdue University, Вашингтонский Университет). Поток газа или жидкости в них может создаваться за счёт действия электромагнитного поля на ионизированную среду, либо «микронасосом» - осциллирующей стенки канала. Эмиттер Коллектор

Постановка задачи LyLy TwTw TwTw T g, V o T g = 0.9 T w, T w = T o, T o – температура торможения

Макроскопические параметры течения вычисляются в некотором малом объёме V 0 пространства по распределению скоростей моделирующих частиц с помощью формул молекулярной динамики: Рассчитываются траектории и скорости каждой частицы, участвующей в движении. x y z Метод прямого численного моделирования

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

Взаимодействие с поверхностью Для описания взаимодействия газа со стенкой использовалась диффузная модель. При этом считалось, что скорости каждой из молекул после отражения не зависят от их индивидуальных скоростей падения, а распределяются согласно равновесной максвелловской функции распределения в том полупространстве скоростей, где вектор скорости молекул направлен от поверхности. Распределение соответствовало температуре стенки Tw. TwTw T = T w, V = 0

Модели взаимодействия Модель со стоком энергии Модель была предложена Бёрдом как модификация модели твёрдых сфер для многоатомных газов. С каждой молекулой связывается переменная, представляющая её внутреннюю энергию: Предполагая равнораспределение энергии, получим зависимость для внутренней энергии молекулы: Температура поступательного движения определяется в виде: Аналогично, определим температуру для внутренних степеней: Это условие проверяется при каждом столкновении. Если величины не равны, то часть разности между ними перераспределяется в направлении удовлетворения уравнения.

Расчётная область плоскость симметрии 0.5 L y A Б TwTw TgTg T T g Так как задача симметрична относительно плоскости равноудалённой от пластин, образующих канал, течение моделировалось по одну сторону от плоскости симметрии. Для этого на ней устанавливалось условие зеркального отражения частиц. изучаемая область зеркальное отражение

Расчётная область ПлотностьСкоростьТемпература Течение считалось стационарным. Для того, чтобы избежать флуктуаций в результатах, проводилось осреднение по времени. Б изучаемая область

Результаты 0.9 T w 1.0 T w V числовая плотность скорость температура

Результаты

Моделирование взаимодействия струи разреженного газа с преградой методами молекулярной динамики Ковалёв В.Л., Сазонова В.Ю., Якунчиков А.Н. Москва 2006 Московский государственный университет им. М.В.Ломоносова, механико-математический факультет

Постановка задачи Падение струи разряженного газа на стенку Газ двухатомный: Молярная масса: Эффективные размеры молекулы: Температура торможения струи: 9H9H H звуковое сопло стенка X

Результаты Распределение плотности 9H9H Распределение скорости Распределение числа Маха на оси симметрии 9H9H

Адсорбция водорода углеродными нанотрубками ЛАБОРАТОРИЯ МНОГОМАСШТАБНОГО МОДЕЛИРОВАНИЯ В.Л. Ковалев, А.Н. Якунчиков

Адсорбция водорода углеродными нанотрубками Введение Автомобиль на топливных элементах автомобильные топливные элементы Масса водорода при хранении его в баллонах составляет примерно 2 – 3% от массы баллона. При хранении водорода в жидком состоянии потери связаны с захолаживанием системы при заправке, а также испарением водорода во время хранения. экспериментальные исследования свидетельствуют о высоком массовом содержании водорода в УНТ Водород является самым энергоемким топливом, продуктом его сгорания является вода. Одно из основных препятствий широкого использования водорода в энергетике это отсутствие эффективных способов его хранения и транспортировки.

Графеновый лист n m (10,10) Нанотрубка (10,10) Углеродная нанотрубка (УНТ) Адсорбция водорода углеродными нанотрубками

Пучки нанотрубок Электронная микроскопия Y. Ye, C.C. Ahn, C. Witham, B. Fultz, J. Liu, A.G. Rinzler, D. Colbert, K.A. Smith, R.E. Smalley нанотрубка (10,10) пучок x1x1 x2x2 (n,n) 1.4 Ǻ графеновый лист Диаметр пучка 6-12 нм Пучок диаметром 10 нм содержит около 50 трубок Площадь внешней поверхности пучка приблизительно в 6 раз меньше площади внешней поверхности входящих в него трубок

H2H2 Взаимодействие атомов C между собой не рассматривалось. H2H2 Газовая фаза: молекулярный водород H 2 – H 2 C– H 2 Физическая адсорбция молекул H 2 на углеродной нанотрубке. Углеродная нанотрубка (УНТ) взаимодействие Адсорбция водорода углеродными нанотрубками С Обмен энергией между внутренними и внешними степенями свободы в молекуле водорода не рассматривался. Взаимодействия описываются потенциалом Леннарда- Джонса Система в термостате

Взаимодействие молекулы H 2 и УНТ ε/k, K σ, A H 2 – H C – H Потенциал Леннарда-Джонса E a / k = 428 K EaEa Адсорбция водорода углеродными нанотрубками x1 x2 x1 x2 r r x1 x2 r

нанотрубка поверхность потенциальной энергии H2H2 Зависимость потенциала между трубкой и молекулой водорода можно приблизить потенциалом Леннарда-Джонса 8-4, который зависит только от расстояния между ними: Адсорбция водорода углеродными нанотрубками

Прямое численное моделирование Уравнение движения для молекулы: - сила, действующая со стороны УНТ - сила, действующая на i -ю молекулу со стороны j -й молекулы Начальные условия: 1.Координаты молекул распределялись регулярно в пространстве 2.Скорости молекул распределялись согласно равновесной функции Максвелла в соответствии с температурой системы:, x y z

Адсорбция водорода углеродными нанотрубками H2H2 H2H2 H2H2 H2H2 H2H2 H2H2 Уравнения решались с постоянным шагом по времени, на котором для каждой молекулы один раз рассчитывалась правая часть. Учитывался только вклад молекул, находящихся в окрестности данной молекулы. Макроскопические параметры течения вычисляются по распределению координат и скоростей молекул:

Адсорбция водорода углеродными нанотрубками Начальное распределение Расчет до выхода на равновесие Нахождение макроскопических параметров осредненем по пространству и времени 1 2 3

Адсорбция водорода углеродными нанотрубками T=80K P=12atm T=80K P=40atm T=80K P=60atm T=298K P=90atm При низких температурах обнаружено образование второго слоя адсорбированных молекул rrr r При комнатной температуре образования второго слоя не наблюдалось Образование второго слоя адсорбции n nn n r УНТ

Адсорбция водорода углеродными нанотрубками Феноменологическая модель основана на теории идеального адсорбированного слоя Ленгмюра, в которой адсорбированные частицы связаны с определенными локализованными центрами на поверхности адсорбента. Предполагалось, что каждый центр может присоединить одну и только одну частицу, энергия адсорбированных частиц на всех центрах поверхности одинакова. Феноменологическая модель =>

Адсорбция водорода углеродными нанотрубками 1 - изотерма Ленгмюра, 2 - прямое численное моделирование для одиночных УНТ, 3 - прямое численное моделирование, пересчитанное для пучков УНТ, 4 - эксперимент Y. Ye, C.C. Ahn, C. Witham, B. Fultz, J. Liu, A.G. Rinzler, D. Colbert, K.A. Smith, R.E. Smalley. (1999) 5 - расчеты N. Hu, X. Sun, A. Hsu. (2005) 6 - эксперимент Lawrence J, Xu Gu (2004) Относительное массовое содержание водорода в УНТ T=80K T=298K

Адсорбция водорода углеродными нанотрубками Адсорбция в массиве трубок Адсорбция на внешней поверхности Адсорбция на одиночных трубках Адсорбция внутри массива Расчеты ведутся на СКИФ МГУ - "Чебышёв"

Спасибо за внимание!