- •7. Особенности моделирования турбулентных течений (4 часа).
- •Уравнение энергии в дифференциальной и интегральной форме
- •Турбулентное движение и осреднение его параметров
- •Основные понятия
- •Правила осреднения
- •Вывод уравнения движения в форме Рейнольдса Уравнение движения вязкой жидкости в форме Навье-Стокса
- •Уравнение движения в форме Рейнольдса
- •Тензор рейнольдсовых напряжений
- •Турбулентность Физические основы турбулентности
- •Гипотеза Буссинеска
- •Длина пути смешения Прандтля
- •Турбулентный пограничный слой Строение турбулентного пограничного слоя
- •Профиль скорости в пограничном слое. Логарифмический профиль
- •Степенной профиль скорости в пограничном слое
- •Алгебраические мт
- •Свободно-турбулентные течения
- •Пристеночные течения
- •Однопараметрические мт
- •Уравнение кинетической энергии турбулентности
- •Диссипация энергии турбулентности
- •Двухпараметрические мт
- •Основные допущения двухпараметрических мт
- •Стандартная k- модель турбулентности
- •Стандартная k- модель турбулентности
- •Модель переноса касательных напряжений sst k-
- •Модели более высокого уровня мт с неизотропной турбулентной вязкостью
- •Модель Рейнольдсовых напряжений rsm
- •Мт с набором фильтрующих уравнений
- •Сравнение rans и les
- •Пример 1. Течение в u-образном колене
- •Пример 2. Теплообмен при обтекании тупоконечной пластины
- •Пример 3. Теплообмен при обтекании лопатки са
Пример 3. Теплообмен при обтекании лопатки са
Проведем моделирование теплообмена при обтекании турбулентным потоком сопловой лопатки с конвективной системой охлаждения. Геометрия расчетной области и заданные граничные условия представлены на рис. 7.15.
Характер обтекания профиля лопатки можно оценить, анализируя рис. 7.16 и 7.17. На рис. 7.16 представлено поле распределения числа Маха в расчетной области, а на рис. 7.17 - поле распределения статического давления в расчетной области
Распределение статического давления, полученного в ходе решения, и его сопоставление с экспериментальными данными представлено на рис. 7.18.
Рис. 7.15 Геометрия расчетной области и заданные граничные условия
Рис. 7.16 Поле распределения числа Маха в расчетной области
Рис. 7.17 Поле распределения статического давления в расчетной области
Рис. 7.18 Распределение статического давления по длине профиля и сопоставление результатов расчета с экспериментальными данными
Распределение температуры металла по длине профиля для вариантов расчета с различными МТ и сопоставление результатов расчетов с экспериментальными данными представлено на рис. 7.19.
Распределение температуры металла по длине профиля, полученное при расчете с применением v2-f МТ и сопоставление результатов расчета с экспериментальными данными представлено на рис. 7.20.
Распределения интенсивности турбулентности для разных МТ представлено на рис. 7.21.
Удовлетворительное совпадение расчетных и экспериментальных данных по распределению температуры стенки по профилю получено при использовании четырехпараметрической МТ v2-f.
Рис. 7.19 Распределение температуры металла по длине профиля: сопоставление результатов расчетов с применением различных МТ с экспериментальными данными
Рис. 7.20 Распределение температуры металла по длине профиля: сопоставление результатов расчета с применением v2-f МТ с экспериментальными данными
Это можно объяснить спецификой генерации кинетической энергии в районе входной кромки. Значительное завышение температуры стенки в районе спинки соответствует недостаточно физичному поведению модели турбулентности в трансзвуковой области и скачке уплотнения.
Рис. 7.21 Распределение температуры металла по длине профиля: сопоставление результатов расчета с применением v2-f МТ с экспериментальными данными
7-
