Метод упругого отскока¶
Метод упругого отскока — численный алгоритм в вычислительной гидродинамике и физике частиц, используемый для моделирования взаимодействия дискретных элементов (частиц, молекул, узлов решётки) с твёрдыми границами или между собой. Метод основан на предположении, что при столкновении частица отражается от поверхности, сохраняя модуль своей скорости, а направление движения изменяется по законам зеркального отражения. Широко применяется в методах решёточных уравнений Больцмана (LBM), методе молекулярной динамики и при симуляции потоков в пористых средах.
¶Физическая основа
В классической механике упругий отскок описывается законом сохранения импульса и кинетической энергии. Для идеально гладкой и неподвижной стенки нормальная составляющая скорости частицы меняет знак, а тангенциальная составляющая остаётся неизменной. В численных методах это правило дискретизируется: частица, пересекающая границу ячейки или попадающая в узел решётки, «разворачивается» относительно нормали к поверхности.
В отличие от методов с поглощением или диффузным отражением, упругий отскок не вносит диссипации энергии в систему, что делает его пригодным для моделирования обратимых процессов и течений с малым числом Кнудсена.
¶Применение в методе решёточных уравнений Больцмана
В LBM метод упругого отскока реализуется на уровне функции распределения частиц. Когда пакет частиц достигает граничного узла, его значения перераспределяются в противоположном направлении. Формально это записывается как:
f_i(x, t+1) = f_j(x, t+1),
где i — направление, противоположное j. Такая схема обеспечивает выполнение условия прилипания (нулевой скорости) на стенке для изотермических течений. Погрешность метода составляет второй порядок точности по пространству, однако положение эффективной границы сдвигается на половину шага решётки, что требует коррекции при расчёте сил.
Разновидностью метода является «отскок с поправкой» (interpolated bounce-back), использующий интерполяцию для повышения точности на криволинейных границах. Также существуют модификации для моделирования частично проскальзывающих поверхностей, где вводится коэффициент аккомодации.
¶Метод отскока в молекулярной динамике
В молекулярной динамике упругий отскок применяется для упрощённого описания столкновений твёрдых сфер. При обнаружении пересечения двух частиц их скорости пересчитываются по формулам центрального удара. Для взаимодействия со стенкой используется зеркальное отражение с возможным учётом температуры стенки (термостат).
Данный подход позволяет существенно ускорить вычисления по сравнению с полным расчётом потенциалов взаимодействия, однако игнорирует деформацию частиц и обмен энергией с внутренними степенями свободы.
¶Преимущества и недостатки
К достоинствам метода относятся:
- простота реализации и высокая вычислительная эффективность;
- точное сохранение массы и импульса в системе;
- отсутствие эмпирических параметров, кроме геометрии поверхности.
Недостатки:
- невозможность моделирования адгезии, шероховатости и неупругих потерь;
- появление артефактов при моделировании течений в микроканалах с большим числом Кнудсена;
- необходимость коррекции положения границы в LBM.
¶Альтернативные подходы
Для учёта неупругих столкновений используются методы с коэффициентом восстановления, диффузное отражение (модель Максвелла) и комбинированные схемы. В LBM также применяется метод «отскок-назад» (backward bounce-back) для моделирования свободных поверхностей и метод «отскок-вперёд» для подвижных границ.
¶Литература
- Сухинов А. И., Чистяков А. Е. Математическое моделирование на основе методов решёточных уравнений Больцмана. — Таганрог, 2012.
- Krüger T. et al. The Lattice Boltzmann Method: Principles and Practice. — Springer, 2017.
- Allen M. P., Tildesley D. J. Computer Simulation of Liquids. — Oxford University Press, 2017.
- Хеер В. Методы молекулярной динамики. — М.: Мир, 1999.