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

молекулярная динамика, никель, потенциал погруженного атома, температура, термостат

Математическое моделирование в актуальных проблемах науки и техники