распространение введенных значений на каждый узел; задание граничных условий; вычисление скорости в каждом узле; задание оператора столкновения и вычисление функции Больцмана на всей решетке; 6) задание градиента внешней силы; 7) циклическое повторение пунктов 3-6, до тех пор, пока система находится в неравновесном состоянии. Все физические параметры обезразмерены в единицах вязкости, единичный шаг по времени означает, что частица проходит характерную толщину вязкого слоя = /, где – есть скорость звука в данной точке, ν – вязкость жидкости. Управляющими параметрами являются число Маха, число Рейнольдса (Re) и число Рэлея (Ra). Были смоделированы двумерные течения Куэтта и Пуазейля, течение в квадратной полости с подвижной верхней стенкой, конвективное течение в той же полости с горизонтальным градиентом температуры. Для иллюстрации работы метода, приведем результаты моделирования двух последних задач. Рассмотрим течение вязкой жидкости в прямоугольной полости с твердыми стенками, верхняя стенка которой движется с заданной постоянной скоростью: проекция компоненты скорости на ось Ox – = 0.2, на ось Oy – = 0, использовалась решетка Nx x Ny=100х100, Re=120. Результатом численного решения методом LBM является векторное поле скорости, приведенное на рис. 1. 2) 3) 4) 5) Рис. 1. Поле скорости в ограниченной полости с подвижной верхней границей Геометрия следующей задачи эквивалентна предыдущей. Рассматривается конвекция вязкой жидкости в замкнутой полости, с нагревом на левой стенке, верхние и нижние границы теплоизолированные. Re = 120, 5 Ra=10 . 98