Прямое Монте-Карло моделирование

Прямое Монте-Карло моделирование (метод прямого статистического моделирования Монте-Карло) — метод вычислительной газодинамики, предназначенный для решения задач динамики разреженных газов. Метод может трактоваться как решение уравнения Больцмана.

Метод ПСМ основан на представлении газа множеством дискретных частиц (каждая из которых представляет собой большое количество реальных молекул), для которых задан стохастический процесс их столкновения друг с другом. Эволюция множества частиц описывается как равномерное прямолинейное движение, прерываемое в случайные моменты времени мгновенными актами парных столкновений, поэтому используются, как правило, модели столкновения с полным конечным сечением. Для упрощения алгоритма и существенного ускорения счёта, фазы перемещения и столкновения частиц разделены между собой и чередуются, а столкновительные партнеры выбираются только в пределах той же самой ячейки (без учёта взаимного расположения).

После достижения стационарного режима течения, макропараметры течения вычисляются осреднением параметров частиц в течение достаточно большого количества временных шагов.

Метод имеет три основных параметра дискретизации: временной шаг $\Delta t$, размер ячейки $\Delta h$ (столкновительные партнеры для каждой частицы выбираются только в пределах той же самой ячейки), число частиц в ячейке $N$. Временной шаг должен быть меньше времени между столкновениями $t_c$, размер ячейки должен быть меньше длины свободного пробега $\lambda$, число частиц в ячейке должно быть достаточно велико, чтобы вероятность повторных столкновений (когда две частицы сталкиваются друг с другом два раза подряд, не столкнувшись с другими частицами) была мала.

Имеет место сходимость второго порядка по $\Delta t / t_c$ (при условии, что частицы редко пролетают за временной шаг более одной ячейки за счёт теплового движения, иначе наблюдается первый порядок), второго порядка по $\Delta h / \lambda$, первого порядка по $\Delta h / (\lambda N)$.

Дисперсия накапливаемых макропараметров уменьшается обратно пропорционально количеству учтённых временных шагов (однако, слишком коротких временных шагов потребуется больше из-за временных автокоррелиряций параметров частиц в ячейке). Т.е. для уменьшения амплитуды погрешности вдвое, требуется рассчитать вчетверо больше временных шагов.

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

Трудоёмкость метода ПСМ непосредственно связана со степенью разреженности газа, которая определяется числом Кнудсена (отношение длины свободного пробега к характерному размеру рассчитываемой системы). Трудоёмкость быстро растёт при уменьшении числа Кнудсена, т.е. с повышением плотности газа, так как требуется мельчить сетку и наращивать число частиц. Ситуация осложняется тем, что установление стационарного режима в более плотном газе происходит дольше, в то время как временной шаг, наоборот, необходимо уменьшать. Как следствие, метод ПСМ применяется, в первую очередь, тогда, когда предположение о предельно малом локальном отклонении газа от равновесия не работает, соответственно, не применимы уравнения Навье‐Стокса, и требуется решение уравнений Больцмана.

История метода

В первые метод прямого статистического моделирования, использующий расщепление по процессам столкновения и переноса молекул был предложен Г. Бердом в 1963 году ^[1]. После этого была предложена схема Бёрда Time-Counter ^[2]. В начале 90-х годов практически все расчёты проводятся с применением схемы Бёрда Non-Time Counter ^[3], либо схемы мажорантной частоты.

Применимость метода

Поскольку разреженный газ -- это газ, в котором вероятность двойных столкновений много больше чем вероятность столкновений высокого порядка (тройных и т.д.), то метод применим для описания течений газа в свободно-молекулярном, переходном и континуальном режимах. Например, воздух удовлетворяет условию разреженности вплоть до давления в сотни атмосфер. Режим течения обычно определяется через число Кнудсена Kn.

Другое ограничение на применимость метода связано с нарушением условия о молекулярном хаосе, которое используется при выводе уравнения Больцмана. Возникновение статистической зависимости между моделирующими молекулами приводит необходимости увеличения числа моделирующих молекул. Для течений в около-континуальном режиме ($Kn<0.01$) этот фактор вынуждает использовать параллельные вычислительные системы

В настоящее время метод прямого статистического моделирования Монте-Карло применяется для исследования течений таких разных масштабов как обтекание космических аппаратов при входе в атмосферы планет, течения газа внутри микро- и наноустройств, течения газа при вакуумных технологических процессах.

Примечания

1. ↑ Bird G.A. Approach to translational equilibrium in a rigit sphere gas. Phys. Fluids. Vol.6, N 10, P 1518—1519 (1963)
2. ↑ Берд Г. Молекулярная газовая динамика. М.: Мир, 1981.
3. ↑ Bird G. A. Molecular gas dynamics and the direct simulation of gas flows. — Clarendon Press, Oxford. — 1994.

Ссылки

• Усовершенствованный метод прямого статистического моделирования для решения современных задач динамики разреженных газов
• Comments on DSMC by Graeme Bird
• DSMC Demo Applet

Рекомендуемая литература

• Берд Г. Молекулярная газовая динамика. М.: Мир, 1981.
• Леонтович М. А. Основные уравнения кинетической теории газов с точки зрения теории случайных процессов // Журн. эксперим. и теорет. физики. 1935. Т. 5. С. 75–79.
• Кац М. Вероятность и смежные вопросы в физике. М.: Мир, 1965.
• Иванов М. С., Рогазинский С. В. Экономичные схемы статистического моделирования течений разреженного газа // Мат. моделирование. 1989. Т. 1, № 7. С. 130–145.
• Галкин В.А., Осецкий Д.Ю. Случай больцмановского газа,приводящий к уравнению коагуляции Смолуховского. //ЖВМ и МФ. 2006.- т.46, N.3. с.536-549.

Библиогр.:9.

• В. А. Галкин, Д. Ю. Осецкий, “Математическое моделирование кинетики коагуляции”, Матем. моделирование, 18:1 (2006), 99–116