Решение проблем расходимости RANS-моделей и оптимизация пограничных слоев в OpenFOAM (SnappyHexMesh)
Обсуждение распространенной проблемы расходимости RANS-расчетов в OpenFOAM, связанной с низким качеством расчетной сетки в пристеночных областях. Рассматриваются причины возникновения искаженных ячеек в призматических слоях SnappyHexMesh и предлагаются методы их устранения для достижения стабильности симуляций.
Решение проблем расходимости RANS-моделей и оптимизация пограничных слоев в OpenFOAM (SnappyHexMesh)
Одна из самых распространенных причин расходимости гидродинамических расчетов — низкое качество расчетной сетки вблизи твердых стенок, где градиенты скоростей максимальны [cite: 13]. Дискуссия из профессионального сообщества инженеров-расчетчиков подробно освещает эту проблему на примере OpenFOAM.
Задача
Достичь сходимости RANS (Reynolds-Averaged Navier-Stokes) симуляции для задач внешней аэродинамики. Утилита checkMesh показывает приемлемую неортогональность (< 65), но выявляет наличие сильно скошенных ячеек (skewness > 4). Проблема локализована в призматических слоях (inflation layers), генерируемых встроенным сеточным генератором SnappyHexMesh (SHM), качество которых выглядит визуально "ужасно" [cite: 13].
Анализ проблемы и Методы решения
Опытные инженеры указали на несколько аспектов:
- Визуальные артефакты постобработки: В ParaView необходимо использовать опцию Crinkle Cut при создании срезов и отключать "decompose polyhedra". Часто то, что выглядит как искаженная сетка, является лишь алгоритмом интерполяции графического движка [cite: 13].
- Физика пограничного слоя: Толстые и плотные слои без плавного перехода (refinement transition) к базовой сетке ведут к ошибкам [cite: 13]. Если ячейка слоя попадает на границу изменения уровня измельчения (refinement level), SHM генерирует искажения.
- Альтернативные методы: Для полного контроля над призматическими слоями предложен скрипт
refineWallLayer. Подход заключается в экструзии только одного толстого слоя через SHM, а затем его программном разделении на множество тонких слоев с заданным коэффициентом роста (bias) и последующем сглаживании узлов (smoothMesh) [cite: 13].
Результат
Корректный расчет $y^+$ (безразмерного расстояния от стенки) и обеспечение плавного перехода объема ячеек от призматического слоя к гексаэдральной сетке позволяет избежать скачков (overshoots) параметра skewness и стабилизировать решатели давления/скорости (например, SIMPLE или PISO) [cite: 13].
Источник
Теги
- #OpenFOAM
- #CFD
- #Meshing
- #RANS
- #SnappyHexMesh