Рис. 6. Линии тока вблизи цилиндра в момент времени
t
= 30
c
низкие скорости, при которых сжимаемостью воздуха можно прене-
бречь.
Как и в случае моделирования обтекания цилиндра, построение
сетки производилось средствами
snappyHexMesh
с предваритель-
ным построением грубой сетки из кубических ячеек в
blockMesh
.
STL-модель была также подготовлена в SALOME, причем для задания
геометрии использовался скрипт на языке Python в силу сложности
задания всех узловых точек профиля крыла вручную. Сетка, получив-
шаяся в результате преобразования ее в двумерную, содержит около
42 000
ячеек (рис. 7).
Поскольку в задаче число Рейнольдса Re
= 10
4
,
необходимо ис-
пользовать ту или иную модель турбулентности. Для проведения те-
стовых расчетов была выбрана стандартная
k
−
ω
SST-модель турбу-
лентности. Для моделирования процесса обтекания снова использо-
вался решатель
pisoFoam
.
Граничные условия задавались аналогич-
но предыдущей задаче: для скорости — условия I рода (
V
(
P
)
=
V
∞
)
на левой (и нижней при ненулевом угле атаки) границе, на остальной
части периметра — граничные условия
outletInlet
,
для давления —
естественные ГУ.
Рис. 7. Сетка в задаче об обтекании профиля крыла
160
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2012