Рис. 5. Распределение скорости
˜
W
z
= ˜
W
2
в поре для случая нерегулярной КЭ
сетки с применением методов Гаусса (
а, в
) и QMR (
б, г
) для случая равномерной
(
а
,
б
) и неравномерной КЭ сеток (
в, г
)
приведены результаты расчетов для продольной скорости
˜
W
z
в поре,
у которой имеется полость, отделенная от сквозного канала перешей-
ком. Вычисления показали (см. рис. 5,
г
и 6,
в
), что действительно в
таких порах формируется застойная зона, в которой скорость течения
практически нулевая, а основное движение газа происходит по сквоз-
ному каналу, причем значение продольной скорости
˜
W
z
слабо зависит
от геометрической формы застойной зоны и определяется только диа-
метром сквозного канала поры. На рис. 5 показано влияние типа КЭ
сетки (адаптивная или равномерная), а также общего числа КЭ на
точность расчетов. Как и ранее, использование адаптивной сетки в
сочетании с методом QMR приводило к лучшим результатам даже для
значительного числа КЭ (около 2000 элементов).
Выводы.
В настоящей работе продолжены исследования, начатые
в [5–10]. В результате дана вариационная формулировка локальной
задачи газовой динамики на ячейке периодичности и на ее основе
100
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2008. № 3
