Для численного моделирования обтекания подвижных роторов вертикально-осевых ветроагрегатов и их элементов применяются осредненные по Рейнольдсу нестационарные уравнения Навье-Стокса несжимаемой жидкости. При замыкании исходных уравнений использованы три однопараметрические модели турбулентной вязкости SA, SARC, SALSA. Решение системы исходных уравнений получено с помощью неявного конечно-объемного численного алгоритма, который базируется на методе искусственной сжимаемости и многоблочных разномасштабных структурированных сетках. Верификация разработанных алгоритмов и комплекса программ выполнена на задачах о развитии течения в каверне, обтекании неподвижного и вращающегося цилиндра. Проведено сравнение однопараметрических моделей турбулентности SA, SARC, SALSA при обтекании профиля NACA 4412. Анализируются мгновенные линии тока, поля завихренности и гистерезисные кривые нестационарных аэродинамических характеристик при обтекании осциллирующего профиля NACA 0015. Показано преимущество модели турбулентности SALSA при закритическом обтекании профиля. Представлены результаты расчета роторов Дарье и Савониуса с различным количеством и геометрическими характеристиками лопастей. Обсуждаются физические особенности обтекания (отрыв пограничного слоя, взаимодействие лопастей, течение в следе и вблизи опорной башни) ветроколеса. Ключевые слова: ветроагрегат, ротор Дарье, Савониуса, отрыв потока, уравнения Навье-Стокса, модель турбулентности
