НаномеханикаНИР

Nanomechanics

Источник финансирования НИР

госбюджет, раздел 0110 (для тем по госзаданию)

Этапы НИР

# Сроки Название
1 1 января 2016 г.-31 декабря 2016 г. Наномеханика
Результаты этапа: 1) Проведено моделирование свободномолекулярного течения газа в микроканале с рядом осциллирующих перегородок. Получено, что вероятность прохождения через канал тем или иным газом существенно зависит от безразмерного соотношения тепловой скорости молекул, расстояния между перегородками и частоты осцилляций. При этом наблюдается аналогия с эффектом резонанса, пик проводимости достигается, когда это соотношение (которое можно рассматривать как отношение некоторых характерных времен или обратных величин - частот) имеет значение порядка единицы. Данный эффект можно использоваться для разделения газов с разными молекулярными массами, т.к при одинаковой температуре они будут иметь разные тепловые скорости. Данный эффект исследован для различных безразмерных геометрических параметров канала, кол-ва перегородок, и для разных моделей рассеяния молекул на стенке канала.Разработан эффективный многомасштабный метод для моделирования течений бинарных газовых смесей при малых числах Маха в современных устройствах, сочетающих макро- и микрокомпоненты. Данный метод объединяет решение уравнения Больцмана в модели МакКормака (McCormack model) для моделирования течений в микроструктурах, и численное решение уравнений газовой динамики в приближении малых чисел Маха для повышения эффективности расчета. Данный метод был применен для моделирования работы прототипа нового газоразделительного устройства для различных бинарных смесей газов. Основным элементом функциональным элементом устройства являются мембраны, в которых за счет малого характерного размера пор под действием противоположно направленных градиентов температуры и давления возникают новые эффекты, приводящие к разделению смеси. Была проведен анализ эффективности устройства от всех безразмерных параметров его работы (геометрия, перепад температур, давлений, характерная скорость потока, число Кнудсена в мембране). Был выдвинут ряд предложений по улучшению прототипа. 2) Исследовано влияние примесей меди на скорость зарождения дислокаций в сплавах алюминия с медью. Зависимости скорости гомогенного и гетерогенного зарождения в сплавах алюминия с медью найдены из результатов молекулярно-динамического моделирования и аппроксимированы в аррениусовской форме. Исследован детальный механизм зарождения дислокационной петли в металлах с гранецентрированной кристаллической решеткой. Показано, что образованию петли частичной дислокации предшествует зарождения наноразмерного структурного перехода ГЦК-ГПУ типа, из которого, затем, формируется дислокационная петля. Показано, что имеются два различных механизма зарождения дислокаций вблизи наноразмерного медного диска в алюминии, зоны Гинье-Престона (ГП): с краев зоны ГП, расположенных вдоль плоскости (111), а также из ее центральной зоны. Показано, что эти механизмы реализуются при температурах до 400 К и выше 400 К, соответственно, и что им соответствует разная активационная энтропия в соответствующих им интервалах температур. Изучены механизмы распространения дислокационной петли в алюминии. Показано, что изначально созданная петля частичной дислокации может либо схлопнуться, либо трансформироваться в полную дислокацию, либо внутри нее образовывается двойник, который заполняет всю ее плоскость. Получены скорости движения краевого и винтового сегмента петли частичной дислокации как функции сдвигового напряжения и температуры. Построена модель, позволяющая моделировать релаксацию сдвигового напряжения в сплавах алюминия с медью, которая учитывает зарождения и распространение дислокационных петель. Методом функционала плотности показано, что образование карбонат-анионов в процессе адсорбции углекислого газа на поверхности цеолитов может приводить к снижению модулей упругости. С помощью методов переходного состояния изучены различные механизмы химической адсорбции углекислого газа на поверхности цеолитов типа X. Получено, что величина барьера активации образования гидрокарбоната из углекислого газа и воды снижается от LiX к NaX, KX и CsX. Показано, что в результате деалюминирования цеолита HY уменьшаются значения его модулей упругости. Методом функционала плотности проведены расчеты средних упругих свойств различных цеолитов и показано, что при наличие в них воды, удается получать результаты, очень близкие к экспериментальным. 3)Исследовано взаимодействие диссоциированного воздуха с поверхностью меди, которая для реальных условий экспериментов в плазматроне моделировалась закисью меди Cu2O. Моделирование поверхности кристаллического оксида меди проводилось на основе кластерного подхода. Кластер Cu2O5 использовался для моделирования процессов адсорбции/десорбции атомов кислорода и азота и ударной рекомбинации, а Cu4O8 и Cu5O8 - для процессов миграции адсорбированных атомов кислорода и азота по поверхности. Численное моделирование основных энергетических,структурных и колебательных характеристик элементарных химических реакций проводилось методами квантовой механики. Для расчетов использовался программный пакет GAUSSIAN. Уравнение Шредингера решалось методом теории функционала электронной плотности с гибридным трехпараметрическим обменным функционалом Беке, дополненный электронной корреляцией Ли, Янга и Пара (B3LYP), и использованием валентно-расщепленного базиса Попла 6-31G*¬ . Для атомов и молекул, участвующих в рассматриваемой задаче (O,N,Cu), указанное выше приближение, позволяет с хорошей точностью описать исследуемые процессы. В рамках теории переходного состояния Эйринга-Поляни были получены конкретные выражения для констант элементарных стадий гетерогенного катализа диссоциированного воздуха на поверхности оксида меди , включающие адсорбцию/десорбцию атомов кислорода и азота, поверхностную диффузию адсорбированных атомов, ударную и ассоциативную рекомбинацию с образованием молекул кислорода, азота и оксида азота. На поверхности рассматривались центры адсорбции с различными энергетическими характеристикаи. Обнаружено, что один из центров адсорбции характеризуется практически безбарьерным характером протекания реакции ударной рекомбинации. Поверхностная миграция атомов моделировалась двумя путями между тремя различными (относительно взаимного расположения) местами адсорбции. Процесс ассоциативной рекомбинации определялся суперпозицией поверхностной миграции и ударной рекомбинации. Полученные эффективные коэффициенты рекомбинации атомов кислорода и азота на поверхности оксида меди для различных парциальных давлений хорошо согласуются с экспериментальными данными. Определена вероятность гетерогенной рекомбинации атомов кислорода и азота в зависимости от парциального давления атомов азота и температуры поверхности в рабочем диапазоне 330-2000°K. Показано, что для атомов кислорода ударная рекомбинация (ИР) преобладает при Т < 1000°K , а ассоциативная (ЛХ) дает вклад в увеличение коэффициента рекомбинации при T > 850°К. Высокая степень каталитичности рассматриваемой поверхности относительно атомов кислорода при высоких и средних температурах (Т > 500°K) вызвана безбарьерным характером ударной рекомбинацией на втором центре адсорбции и является частной особенностью взаимодействия именно данной системы Cu2O+O. Показано, что ударная рекомбинация при малых значениях парциального давления атомов азота (PN=0.1 Па) дает основной вклад в вероятность гетерогенной рекомбинации при Т < 800°K, а механизм ассоциативной рекомбинации является определяющим при высоких температурах (T > 800°К). При увеличении парциального давления (PN=105 Па) значение температуры в точке смены приоритетов механизмов гетерогенной рекомбинации повышается до 1400°К. Относительный вклад ассоциативной рекомбинации в вероятность гетерогенной рекомбинации растет с уменьшением парциального давления.
2 1 января 2017 г.-31 декабря 2017 г. Наномеханика
Результаты этапа: Проведено моделирование свободномолекулярного течения газа в микроканале с границей, движущейся по волновому закону. Получено, что вероятность прохождения через канал тем или иным газом существенно зависит от безразмерного соотношения характерной тепловой скорости молекул, и волновой скорости движения границы. Дано теоретическое объяснение данному эффекту, и показано, как он может быть использован как для разделения газов, так и для создания газовых микронасосов. Показано, что эффект разделения наблюдается лишь при определенном соотношении амплитуды волн к диаметру канала, и эффективность разделения повышается с увеличением длины канала, а также с увеличением длины волны. Также показано, что наибольшая эффективность достигается , когда нижняя и верхняя граница канала движутся в одной фазе. Проведено моделирование течения двухкомпонентной смеси в неизотермической постановке в новом газоразделительном устройстве, использующем мембрану с приложенным градиентом температуры. Определено влияние на эффективность разделения таких параметров как скорость потока, толщина и длина мембраны, диаметр канала внешней части устройства, а также перепад температур на сторонах мембраны. Была выведена модель МакКормака уравнения Больцмана для течений бинарных смесей в прямом канале под влиянием перепадов температуры, давления и концентрации компонент. Ведется работа над расширением модели на случай трех и более компонент, а также над реализацией численного метода для решения систем уравнения для модели бинарной смеси. Построена модель релаксации сдвиговых напряжений в сплавах алюминия с медью, учитывает зарождения и распространение дислокационных петель. Показано, что при релаксации сдвигового напряжения скорость пластической деформации является функцией не только времени, но и величины пластической деформации. Установлено, что зависимости плотности дислокаций от скорости пластической деформации практически не зависят от температуры, что подтверждается экспериментальными данными. Показано, что время релаксации сдвигового напряжения, не превышающего 0.85 ГПа, увеличивается с ростом температуры, а при напряжениях более 0.85 ГПа, снижается. Показано, что зарождение дислокаций является основным механизмом пластической деформации при скоростях нагружения выше 107 с-1 для сплавов алюминия с начальной плотностью дислокаций 1010-1011 м-2. Методом функционала плотности проведены расчеты упругих свойств различных цеолитов с учетом влияния воды. Показано, что вода снижается упругие характеристики цеолитов с одновалентными катионами, и, наоборот, повышает их у цеолитов с двухвалентными катионами. Рассчитанные осредненные по Фойгту-Ройссу-Хиллу объемные модули и модули Юнга цеолитов согласуются с экспериментальными данными. Методами квантовой молекулярной динамики проведены расчеты динамики углекислого газа в поре цеолита NaKA в присутствии молекул воды, карбонат анионов и дикарбонатов. Показано, что и два карбонат-аниона, и дикарбонат могут оттягивать из 8-членного окна катион К так, что молекулы СО2 проникают это окно в другую пору. При моделировании динамики углекислого газа и воды в поре цеолита NaKA было установлен механизм образования гидрокарбонат аниона при реакции H2O и CO2, который был ранее предложен на основе статических расчетов методом переходного состояния. Разработан метод определения полей напряжений в теле конечного размера, которое содержит краевые дислокации, лежащие на разные плоскостях скольжения. Искомое поле напряжений находится суперпозицией напряжений, которые создают дислокации в бесконечном пространстве, и напряжений “изображения”. Напряжения “изображения” определяются при решении краевой задачи теории упругости, в которой из исходных граничных условий вычитаются упругие поля (перемещения и напряжений) от всех дислокаций на границе тела. Исследовано взаимодействие диссоциированного воздуха с поверхностью меди, которая для реальных условий экспериментов в плазматроне моделировалась закисью меди Cu2O. Моделирование поверхности кристаллического оксида меди проводилось на основе кластерного подхода. Кластер Cu2O5 использовался для моделирования процессов адсорбции/десорбции атомов кислорода и азота и ударной рекомбинации, а Cu4O8 и Cu5O8 - для процессов миграции адсорбированных атомов кислорода и азота по поверхности. Численное моделирование основных энергетических,структурных и колебательных характеристик элементарных химических реакций проводилось методами квантовой механики. Для расчетов использовался программный пакет GAUSSIAN. Уравнение Шредингера решалось методом теории функционала электронной плотности с гибридным трехпараметрическим обменным функционалом Беке, дополненный электронной корреляцией Ли, Янга и Пара (B3LYP), и использованием валентно-расщепленного базиса Попла 6-31G*¬ . Для атомов и молекул, участвующих в рассматриваемой задаче (O,N,Cu), указанное выше приближение, позволяет с хорошей точностью описать исследуемые процессы. В рамках теории переходного состояния Эйринга-Поляни были получены конкретные выражения для констант элементарных стадий гетерогенного катализа диссоциированного воздуха на поверхности оксида меди , включающие адсорбцию/десорбцию атомов кислорода и азота, поверхностную диффузию адсорбированных атомов, ударную и ассоциативную рекомбинацию с образованием молекул кислорода, азота и оксида азота. На поверхности рассматривались центры адсорбции с различными энергетическими характеристикаи. Обнаружено, что один из центров адсорбции характеризуется практически безбарьерным характером протекания реакции ударной рекомбинации. Поверхностная миграция атомов моделировалась двумя путями между тремя различными (относительно взаимного расположения) местами адсорбции. Процесс ассоциативной рекомбинации определялся суперпозицией поверхностной миграции и ударной рекомбинации. Полученные эффективные коэффициенты рекомбинации атомов кислорода и азота на поверхности оксида меди для различных парциальных давлений хорошо согласуются с экспериментальными данными. Определена вероятность гетерогенной рекомбинации атомов кислорода и азота в зависимости от парциального давления атомов азота и температуры поверхности в рабочем диапазоне 330-2000°K. Показано, что для атомов кислорода ударная рекомбинация (ИР) преобладает при Т < 1000°K , а ассоциативная (ЛХ) дает вклад в увеличение коэффициента рекомбинации при T > 850°К. Высокая степень каталитичности рассматриваемой поверхности относительно атомов кислорода при высоких и средних температурах (Т > 500°K) вызвана безбарьерным характером ударной рекомбинацией на втором центре адсорбции и является частной особенностью взаимодействия именно данной системы Cu2O+O. Показано, что ударная рекомбинация при малых значениях парциального давления атомов азота (PN=0.1 Па) дает основной вклад в вероятность гетерогенной рекомбинации при Т < 800°K, а механизм ассоциативной рекомбинации является определяющим при высоких температурах (T > 800°К). При увеличении парциального давления (PN=105 Па) значение температуры в точке смены приоритетов механизмов гетерогенной рекомбинации повышается до 1400°К. Относительный вклад ассоциативной рекомбинации в вероятность гетерогенной рекомбинации растет с уменьшением парциального давления.
3 1 января 2018 г.-31 декабря 2018 г. Наномеханика
Результаты этапа: 1.Разработан алгоритм построения атомных моделей нанотрубок с цилиндрической анизотропией упругих свойств различной хиральности. Реализована программная реализация алгоритма. Рассчитаны коэффициенты Пуассона нанотрубок из железа и меди, свернутых из плоскости (010) под различными углами хиральности, Показано, что с увеличением угла хиральности снижаются значения коэффициента Пуассона, и при угле хиральности выше 33 градусов он становится отрицательным. Установлено, что нанотрубки из Fe и Cu проявляют линейный эффект Пойнтинга, то есть закручиваются при растяжении. Максимальный эффект достигается при углах хиральности около 25 градусов. Реализован алгоритм метода дискретных дислокаций, который позволяет рассчитывать характеристики пластической деформации тел, имеющие нано- и микроразмеры. Показано, что барьер образования дикарбоната находится в пределах 0.3 эВ, что является достаточным для активации такой реакции при комнатной температуре. Получено, что катион K существенно выходит из плоскости 8-членного окна при координации к нему дикарбоната, карбоната или гидрокарбоната, позволяя молекулам газа проходит через это окно. Этот результат подтверждает, что именно процессы химической адсорбции имеют существенно влияние на характер разделение газовых смесей, включающий углекислый газ, в цеолитах. 2.Получены новые результаты о течении смеси газов в микроканале, стенки которого вынужденно изгибаются по волновому закону. Решена задача о медленном течении многокомпонентного газа в переходном и свободномолекулярном режимах в области со сложной границей и передами температуры на стенках. Показано, что сложная геометрия в сочетании с градиентом температуры приводит к эффекту прокачки газа и разделения смеси. Созданный в 2017 гибридный численный метод для моделирования течения бинарных смесей в устройствах с разными режимами течения (континуальный, переходный, свободномолекулярный) был развит на случай течения смесей из трех компонент. В настоящий момент с помощью данного метода ведется расчет задачи о течении смеси из трех в новом газоразделительном устройстве, представляющем собой U-образный канал с помещенной между внутренними сторонами мембраной или рядом нагретых нитей (филаментов). Показано, что при течении в такой системе обнаруживается эффект разделения смеси. Также метод событийного молекулярно-динамического моделирования был обобщен на случай наличия границ, движущихся с постоянными, но различными скоростями. 3.Проведены исследования обменных процессов между угарным газом и оксидами азота CO+N2O  СO2+N2, CO+NO2  СO2+NO, приводящим к образованию СO2. Методами квантовой механики получены переходные состояния, частоты колебаний и пути реакций, а также их ключевые энергетические характеристики. В рамках теории переходного состояния были рассчитаны константы скоростей прямых и обратных реакций. Проведено их сравнение с имеющимися в литературе экспериментальными данными. Представлены их аппроксимации в виде Аррениуса в широком диапазоне температур 300-2500 K. Проведены расчеты константы скорости диссоциации N2+N → 3N для различных уровней описания с использованием вычислительного комплекса MD Trajectory и интернет каталога теоретических моделей.Проведено сравнение результатов траекторных расчетов с теоретической моделью Мэрроуна-Тринора по уровневому фактору. Получено существенное отличие вплоть до 4 порядков при невысоких поступательных температурах для высоких колебательных уровней молекулы N2. Для описания траекторных расчетов предложена модификация модели Мэрроуна-Тринора с введением температурной зависимости параметра q(T). Проведено сравнение результатов траекторных расчетов и теоретических моделей с имеющимися экспериментальными данными для процесса диссоциации N2 + N → 3N. Показано, что наиболее адекватно экспериментальные данные описывает двухтемпературная модель Мэрроуна-Тринора.
4 1 января 2019 г.-31 декабря 2019 г. Наномеханика
Результаты этапа: 1) В настоящее время актуален поиск материалов и функциональных компонент, которые имеют необычные механические свойства на микро- и наномасштабных уровнях. Изучены механические свойства хиральных нанотрубок из металлов кубической симметрии (Fe, Cu, Co, Al), свернутых из плоскости (010) под различными углами хиральности. Разработан алгоритм построения атомных моделей нанотрубок с цилиндрической анизотропией упругих свойств различной хиральности. Установлено, что при определенных значениях углов хиральности трубки из Fe, Cu и Сo могут иметь ауксетические свойства, а также связь между одноосной и крутильной деформацией. Показано, что наличие ауксетических свойств у нанотрубок из кобальта является следствием влияния свойств поверхности, так как сами кристаллы кобальта ауксетиками не являются. Исследованы зависимости механических свойств нанотрубок от их геометрических параметров. Показано, что с уменьшением толщины медных нанотрубок с 2.5 нм до 0.8 нм коэффициент связи между одноосной и крутильной деформацией увеличивается в три раза. 2) Эффективное разделение газовых смесей СО2 является одной из центральных проблем улавливания и хранения углерода. Катионные формы цеолитов, содержащие в своей структуре восьмичленные окна 8R или двойные призмы из них (D8R), имеют исключительно высокую селективность по СО2. В работе изучается динамика катионов при физической и химической адсорбции СО2 с образованием карбонатов методом квантовой молекулярной динамики. Образование гидрокарбоната HCO3- было подтверждено низким активационным барьером 0.15 эВ, полученным при моделировании реакции между СО2 и водой на поверхности цеолита NaKA методом NEB (nudged elastic band). В результате расчетов определялись траектории молекул СО2, анионов HCO3-, CO32-, и катионов Na и K. Показано, что карбонаты образуют комплексы с катионами за счет кулоновского притяжения и вытягивают их окон цеолита, открывая проход для молекул газа. При этом гидрокарбонаты HCO3- вытягивают катион К из 8R окна на большее расстояние, чем CO32-: смещения К в направлении нормали к 8R окну составляют 1.2 А для CO32- и 1.6 А для HCO3-. Показано, что наличие в ячейке цеолита молекул углекислого газа не влияет на взаимодействие между HCO3- и катионом К, в результате которого катион вытягивается из окна, открывая его для прохождения молекул СО2. Обнаружено, что наличие одного лишь углекислого газа в ячейке цеолита не позволяет изменить среднее положение катиона К, который колеблется вблизи равновесного положения в 8R окне. По результатам моделирования сделан вывод о том, что эксперименты по измерению индивидуальной адсорбции компонент газовых смесей, содержащих СО2, не позволяют судить о способности цеолита к разделению этих компонент. 3)Исследована плоская задача о течении смеси из 5 инертных газов в разных режимах (континуальный, переходный, свободномолекулярный) в системе, состоящей из нескольких рядов натянутых нитей (филаментов), имеющих различную температуру. Получены распределения температуры, давления, плотности и концентраций компонент в системе. Показано, что под действием перепада температур возникает эффект разделения смеси. Числовые концентрации компонент увеличиваются в направлении холодной стороны камеры между ступенями, а в самих ступенях снижаются скачком. Показано, что разница в молярных концентрациях вырабатывается именно в камерах между ступенями, а в самих ступенях нитей состав практически не меняется. Получено, что наибольший эффект разделения достигается при условии, когда режим течения в ступени приближается к свободномолекулярному, а в камере между ступенями ‒ наоборот, к континуальному. Предложен и исследован новый принцип работы разделительного устройства с подвижной границей, коэффициент разделения которого существенно выше, чем у диффузионного метода. Принцип работы данного устройства основан на организации периодически прерываемого истечения смеси газов в вакуум. Показано, что за счет нестационарности истечения можно получить коэффициент разделения существенно выше, чем у диффузионного метода. Причем это возможно уже при скоростях подвижных элементов устройства порядка нескольких сотых от характерной тепловой скорости молекул. Такие скорости значительно легче достичь в микротехнике (по сравнению со скоростями, характерными для турбомолекулярных насосов). Исследована задача о течении газа в сужающемся микроканале в свободномолекулярном и переходном режимах под действием градиента температуры на стенках. Показано, как наклон стенок влияет на величину эффекта температурной транспирации. Получено, что без приложения градиента температуры на стенках (равномерно нагретые стенки), только за счет их наклона, поток газа создать не удается. Изучен процесс столкновения и перераспределения энергии двухатомных молекул азота и кислорода в широком диапазоне прицельных параметров (поступательной и вращательной энергии, прицельного расстояния, ориентации молекул, фаз вращения) при температурах до 2000 K. Для решения задачи проведена серия из большого числа (~〖10〗^13) молекулярно-траекторных расчетов парных столкновений молекул. Для моделирования использовались классические потенциалы Леннарда-Джонса (между атомами разных молекул) и упругий потенциал (между атомами одной молекулы). В результате анализа полученных результатов выявлены свойства, которые не отмечаются в существующих классических моделях. Так, показано, что угол разлета молекул (определяющий сечение столкновений), и коэффициент обмена энергией в столкновении зависят не только от поступательной энергии относительного движения молекул (как в существующих моделях), но и от вращательных энергий обеих молекул. При этом соотношения вращательных энергии молекул до столкновения также играет важную роль, как при рассеянии, так и при обмене энергией. С применением теории функционала плотности (DFT) были рассчитаны поверхности потенциальной энергии (ППЭ) для пар N_2-N_2,O_2-O_2,N_2-O_2. Расчеты проводились при равновесном значении длины межатомной связи, но с высокой детализацией по углам взаимной ориентации молекул (10 градусов) и межмолекулярного расстояния (0.5 Ангстрем), для детального описания анизотропии взаимодействия двухатомных молекул. В общей сложности, для каждой из трех пар ППЭ была рассчитана в 18810 точках. Показано, что использование произвольных парных потенциалов не позволяет аппроксимировать результаты DFT расчетов с точностью выше 5-6%. В связи с этим, был предложен анизотропный многочастичный потенциал, приближающий полученные DFT результаты с точность 0.75%. В дальнейшем планируется сравнить результаты MD моделирования для полученного потенциала с результатами при использовании парного потенциала Леннарда-Джонса. 4) Механизм обменных реакций CO+N<->CN+O, CO+N<->NO+C, CN+O <->NO+C, протекающих в условиях входа в атмосферу Марса, был исследован на основе DFT моделирования. В результате расчета поверхности потенциальной энергии для рассматриваемой системы были найдены пути реакций, переходные состояния, промежуточные комплексы и соответствующие им частоты колебаний и энергетические характеристики. Показано, что механизмы всех реакций являются многоступенчатыми, и идут через образование нескольких последовательных промежуточных комплексов. Рассчитаны константы скорости исследуемых реакций в широком диапазоне температур и проведено их сравнение с имеющимися экспериментальными данными. Проведено исследование обменного процесса между угарным газом и оксидом азота CO+N2O<->СO2+N2, приводящим к образованию СO2. Методами квантовой механики рассчитаны переходные состояния, частоты колебаний и пути реакции. Рассмотрен как прямой перенос атома кислорода, так и многостадийный процес, идущий через образование нескольких промежуточных комплексов. В рамках теории переходного состояния рассчитаны константы скорости прямой и обратной реакции. Проведен критический анализ имеющихся в литературе экспериментальных данных на основе сравнения с результатами проведенных расчетов. Рекомендуемые константы скорости представлены в обобщенной форме Аррениуса в широком диапазоне температур 300-2500 K. Разработана кинетическая схема образования молекул NO в основном и возбужденных состояниях при взаимодействии диссоциированного воздуха с поверхностью бета-кристобалита (SiO2). Построена кластерная модель бета-кристобалита, учитывающая различные активные центры физической и химической адсорбции. Проведены квантово-механические расчеты поверхности потенциальной энергии для моделирования различных процессов на поверхности SiO2- адсорбция/десорбция атомов и молекул (N, O, NO, N2, O2), ударная и ассоциативная рекомбинация и получены константы скоростей для всех рассматриваемых процессов.
5 1 января 2020 г.-31 декабря 2020 г. Наномеханика
Результаты этапа: 1. (Косьянчук В.В.)Исследована плоская задача о нестационарном свободномолекулярном течении газа в микроканале с рядом осциллирующих затворов. Показано, что изменение частоты колебаний, коичества завторов и расстояний между ними приводит к характерно различающимся профилям зависимости потока от времени, что может использовано в различных областях микро- и вакуумной техники. Также показано, что такое движение системы приводит к эффекту разделения при течении через канал газовой смеси, и более того, эффект разделению намного более выражен, чем при течении в стационарном режиме (изучавшемся в предыдущих работах), а также достигается при существенно меньших (на 1-2 порядка) частотах колебаний. Проведено моделирование течения разреженного газа в микроканале со стенками движущимися по закону стоячей волны, проведено сравнение с предыдущими результатами, где стенки двигались по закону бегущей волны. Получено, что режимы течения при движении стенок по двум разным законам характерно различаются. Однако движение стенок по закону стоячей волны не приводит к усилению эффекта разделения по сравнению с бегущей волной, как ожидалось. На основании данных квантово-механических и молекулярно-динамических расчетов, проведенных в прошлом году была получена новая двух-параметрическая зависимость для вероятности неупругого столкновения в модели обмена энергией Боргнакке-Ларсена. На примере задач о нольмерной релаксации было показано преимущество данной модели по сравнению с классическими (например моделью Бойда), а также недавней двух=параметрической моделью Жанга. Данная модель была встроена в методы прямо статистического моделирования Монте-Карло (DSMC) и событийного молекулярно-динамического моделирования (EDMD) и была решена задача об истечении разреженного газа (азота) в вакуум. Полученные результаты хорошо согласуются с экспериментом, что подтверждает корректность полученной модели. Также в ходе исследования было выведено интегральное соотношение, связывающее вероятность неупругого столкновения в модели Боргнакке-Ларсена и параметр вращательной релаксаций Zrot. На основании этого соотношения показано, что параметр Zrot в неравновесных процессах не может быть описан одной температурой, и требует двух температур (вращательной, поступательной для корректного описания). Помимо этого, была выявлена интересная особенность процесса вращательно-поступательного обмена энергией на молекулярном уровне – оказалось, что при больших скоростях вращения молекул интенсивность энергообмена увеличивается с ростом относительной скорости столкновения двух молекул, в то время как при малых скоростях вращения, энергообмен с ростом относительной скорости столкновения двух молекул ослабевает. 2. (Брюханов И.А.)С помощью молекулярно-динамического моделирования исследована подвижность краевой дислокации в сплавах твердого раствора Cu-Ni в широком диапазоне температур 100-1100 К и концентраций Ni 0-30%. Показано, что существует два разных режима влияния атомов Ni замещения на движение краевой дислокации. Когда движение дислокации обусловлено силами фононного трения, то атомы Ni действуют как барьеры и замедляют движение дислокации. С повышением температуры сила торможения, связанная с атомами Ni, уменьшается и становится незначительной, начиная с 700 К. С ростом напряжения движение дислокации переходит в режим «насыщения» и ее скорость приближается к скорости звука в материале. В этом режиме, атомы Ni увеличивают скорость движения дислокации, что связано увеличением модуля сдвига твердого раствора Cu-Ni с концентрацией Ni. Показано, что ширина дислокационной линии увеличивается с ростом температуры и уменьшается с увеличением концентрации никеля. На основе полученных результатов сделано предположение, что за фронтом ударной волны в сплавах меди с никелем может достигаться более высокая скорость пластической деформации, а, следовательно, и более быстрая релаксация напряжений, чем для чистой меди. 3. (Брюханов И.А.) Показано, что наличие подложки сжимает слои Pd пленки, тем самым увеличивая энергию диссоциации метана и диффузионный барьер атомов кислорода внутрь поверхности палладия. Показано, что самый минимальный барьер диссоциации CH4 достигается на дефектной подложке (100)γ-Al2O3 среди подложек из других оксидов (ZrO2, Al2O3, шпинели MeAl2O4). Установлено, что добавление атомов Co, Ni, Mn в подложку Al2O3 снижает барьер диссоциации CH4, причем степень влияния атомов соответствует порядку Co > Mn > Ni, в котором для атомов Ni рассчитано наибольшее снижение энергии диссоциации (увеличение скорости диссоциации, соответственно). Проведено сравнение результатов расчетов с экспериментом по окислению O2 на поверхности палладия. Найдено отличие в степени влияния атомов на скорость окисления Ni > Co > Mn. что связано с более высокой температурой окисления O2 для Pd/MnAl2O системы. 4. (совместно с лабораторией 109 Крупнов А.А., Погосбекян М.Ю.)Проведено исследование процессов адсорбции, десорбции и рекомбинации атомов N и O на поверхностях теплозащитных материалов SiO2 и SiC методами квантовой механики в рамках теории функционала плотности DFT (Х3LYP/cc-pVTZ). Для описания гетерогенных процессов использовался кластерный подход, в котором поверхность кристаллов SiO2 и SiC моделировалась различными кластерами, передающими стехиометрию кристалла и валентные состояния атомов, лежащих на верхнем слое. В ходе расчетов были определены ключевые структурные и энергетические характеристики молекулярных систем, описывающих взаимодействия атомов N и O с кремне содержащими поверхностями. В рамках метода переходного состояния с использованием квазистационарного приближения были рассчитаны константы скорости для процессов адсорбции и десорбции атомов N, O и молекул N2, O2, NO, а также ударной рекомбинации атомов N, O на рассматриваемых поверхностях в широком диапазоне температур 500-2200 K. Полученные константы скорости аппроксимированы обобщенной формулой Аррениуса. Таким образом, построена замкнутая кинетическая модель для исследования высокотемпературных физико-химических процессов, протекающих на поверхности тепловой защиты космических и гиперзвуковых летательных аппаратов на основе SiO2 и SiC.
6 1 января 2021 г.-31 декабря 2021 г. Наномеханики
Результаты этапа:
7 1 января 2022 г.-31 декабря 2022 г. Наномеханика
Результаты этапа:
8 1 января 2023 г.-31 декабря 2023 г. Наномеханика
Результаты этапа:
9 1 января 2024 г.-31 декабря 2024 г. Наномеханика
Результаты этапа:

Прикрепленные к НИР результаты

Для прикрепления результата сначала выберете тип результата (статьи, книги, ...). После чего введите несколько символов в поле поиска прикрепляемого результата, затем выберете один из предложенных и нажмите кнопку "Добавить".