Аннотация
Существует два метода получения информации о коэффициентах переноса газов и жидкостей: экспериментальное их измерение и молекулярное моделирование. Только для разреженного газа в простейших случаях построена кинетическая теория, которая завершает молекулярное моделирование явными формулами для коэффициентов переноса. Уже в плотном газе и в жидкости для расчета коэффициентов переноса приходится использовать метод молекулярной динамики. Метод этот требует, однако, огромных вычислительных ресурсов. Создание альтернативного менее затратного метода молекулярного моделирования поэтому весьма актуально. В данной статье предлагается стохастический метод молекулярного моделирования коэффициентов переноса в плотных газах. В основе метода лежит идея расщепления процесса движения молекул газа на свободный пролет и их столкновения. Свободный пролет молекул осуществляется в соответствие с их скоростями в данный момент времени, а столкновения разыгрываются случайным образом. Для этого сначала определяется может данная молекула столкнуться или нет. Затем, если столкновение реализуется, то также случайным образом из оставшихся молекул системы определяется молекула, с которой данная соударяется. В результате моделирования получаются данные о координатах и импульсах всех молекул рассматриваемой системы в последовательные моменты времени. Затем, используя методы неравновесной статистической механики, по этим данным вычисляются коэффициенты переноса. Они при этом получаются усреднением соответствующих динамических переменных по ансамблю независимых фазовых траекторий. Работоспособность алгоритма продемонстрирована на примере моделирования коэффициента вязкости аргона в широком диапазоне давлений (плотностей). Показано, что точность метода растет с увеличением числа молекул в моделируемой системе и числа используемых для усреднения фазовых траекторий.
Ключевые слова: процессы переноса, плотный газ, стохастическое моделирование, молекулярное моделирование, вязкость, статистическое моделирование