ПОЛУЧЕНИЕ ИЗОТЕРМИЧЕСКИХ ХАРАКТЕРИСТИК, ПАРАМЕТРОВ УРАВНЕНИЯ СОСТОЯНИЯ ДЛЯ PETN МЕТОДАМИ РЕАКЦИОННОЙ МОЛЕКУЛЯРНОЙ ДИНАМИКИ И ТЕМОДИНАМИКИ

Abstract

Методом молекулярной динамики (МД) в программном пакете LAMMPS (Large-scale Atom-ic/Molecular Massively Parallel Simulator) с использованием реакционного силового поля ReaxFF-lg были рассчитаны изотермы нереагирующего вещества PETN в диапазоне давлений до 30 ГПа. Были получены значения коэффициента модуля всестороннего сжатия Ко = 9,6 ГПа и производной модуля сжатия Ко по давлению К'о = 8,0, которые можно использовать как параметры термического уравнения Берча-Мурнагана 3-го порядка. Были подобраны коэффициенты уравнения состояния (УрС) в форме Ми-Грюнайзена, часто используемого для моделирования теплофизических свойств вещества, в том числе при статическом и ударно-волновом сжатии. Для нахождения коэффициентов УрС применялся метод построения изохорно-изотермического потенциала твердых веществ в форме квазигармонического приближения Эйнштейна. Верификация полученных результатов показала хорошее согласие с экспериментальными данными в широком диапазоне изменения давления и температуры, в том числе вдоль ударной адиабаты.