Заседание семинара "Высокопроизводительные вычисления" 12.03.2020 в 11-00. "Уравнение состояния плотных реагирующих газов для описания ударных и детонационных волн", Прууэл Эдуард Рейнович

Cеминар ИВМиМГ СО РАН: 
Высокопроизводительные вычисления
Руководитель семинара: 
д.т.н. Глинский Б.М.
Дата / Время проведения: 
четверг, 12 марта, 2020 - 11:00
Место проведения: 
Конф.-зал ИВМиМГ СО РАН
Докладчик
Ф.И.О. докладчика: 
Прууэл Эдуард Рейнович
Ученая степень: 
К.ф.-м.н.
Место работы: 
ИГиЛ СО РАН
Название доклада: 
Уравнение состояния плотных реагирующих газов для описания ударных и детонационных волн
Аннотация доклада: 

Реализована модель уравнения состояния позволяющая проводить расчеты равновесных термодинамических параметров плотных газов и флюидов с учетом межмолекулярных взаимодействий. Протестированный диапазон параметров составляет по температуре -- от 100 до 10 000 K и до давлений 100 ГПа. Базовыми параметрами для задания условий являются плотность, температура и химический элементный состав исследуемой смеси. Для этих условий вычисляется равновесный химический состав, внутренняя энергия, теплоемкости, показатель адиабаты и скорость звука. В расчетах учтена возможность образования следующих химических компонент: Ar, Ne, He, Kr, N2, N, O2, O, H2, H, H2O, OH, NH3, CO, CO2, CH4, NH3, NO и конденсированной фазы углерода.

Для построения уравнения состояния используются методы молекулярной динамки и статистической физики. При этом, вещество рассматривается как набор точечных объектов (молекул) взаимодействующих с центральным парным потенциалом в форме exp6. Дополнительно молекулы обладают внутренними степенями свободы, энергия  которых зависит только от температуры. Для определения давления и полной энергии системы численно решается задача движения небольшого ансамбля частиц (NVT ансамбль из 50-100 молекул), при этом внутренняя энергия системы вычисляется как суммы кинетических энергий молекул и потенциальной энергии взаимодействия, давление вычисляется по Теореме о вириале. Подбор параметров парных потенциалов взаимодействия осуществлялся из условий наилучшего соответствия экспериментальным данным: таблицам термодинамическим величин Американского института стандартов, ударным адиабатам сжиженных газов и параметрам детонационным конденсированных взрывчатых материалов.

Разработанный на базе модели программный комплекс (http://ancient.hydro.nsc.ru/chem) позволяет определять равновесные термодинамические параметры смесей газов при заданной плотности и температуре, рассчитывать ударные адиабаты сжиженных газов и определять парамтры детонации конденсированных взрывчатых веществ. Для построения ударной адиабаты, численно решается нелинейное уравнение Гюгонио в переменных плотность и температура, при этом, все необходимые параметры (давление и удельная внутренняя энергия) вычисляются описанным выше методом. Параметры детонации определяются из условий Чепмена-Жуге - на ударной адиабате с энерговыделением находится точка с условиями D=u+c, где D - скорость фронта, u - массовая скорость, c - равновесная скорость звука.

Прикрепленные файлы: