EISSN 1726-3522
Язык: ru

ЯВНЫЕ СХЕМЫ ВЫСОКИХ ПОРЯДКОВ ТОЧНОСТИ ДЛЯ ЗАДАЧ МОЛЕКУЛЯРНОЙ ДИНАМИКИ (2021)

Рассмотрены явные симплектические разностные схемы Рунге-Кутты-Нистрема (RKN) с числом стадий от 1 до 5 для численного решения задач молекулярной динамики, описываемых системами с распадающимися гамильтонианами. Для числа стадий 2 и 3 параметры RKN-схем получены с помощью техники базисов Гребнера. Для числа стадий 4 и 5 новые схемы най дены с применением метода численной оптимизации Нелдера-Мида. В частности, для числа стадий 4 получены четыре новые схемы. Для числа стадий 5 получены три новые схемы в дополнение к четырем схемам, известным в литературе. Для каждого конкретного числа стадий найдена схема, являющаяся наилучшей с точки зрения минимума ведущего члена погрешности аппроксимации. Верификация схем осуществлена на задаче, имеющей точное решение. Показано, что симплектическая пятистадийная RKN-схема обеспечивает более точное сохранение баланса полной энергии системы частиц, чем схемы более низких порядков точности. Исследования устойчивости схем выполнены с помощью программного пакета Mathematica.

Тип: Статья
Автор (ы): Ворожцов Евгений Васильевич, Киселев Сергей Петрович
Ключевые фразы: молекулярная динамика, УРАВНЕНИЯ ГАМИЛЬТОНА, СИМПЛЕКТИЧЕСКИЕ РАЗНОСТНЫЕ СХЕМЫ, устойчивость

Идентификаторы и классификаторы

УДК
519.63. Численные методы решения дифференциальных уравнений с частными производными
532.51. Общие вопросы движения жидкостей. Гидродинамика в собственном смысле слова
eLIBRARY ID
46169015
Текстовый фрагмент статьи