Method of calculating the aerodynamic efficiency of the axial fan.pdf В инженерных методиках, используемых для расчета осевых вентиляторов на стадии их проектирования, прибегают к определенным допущениям [1, 2], позволяющим, с одной стороны, упростить вычисления, а с другой - в целом правильно оценить основные характеристики вентилятора. Дальнейшее уточнение характеристик производится путем введения поправочных коэффициентов, получаемых из имеющейся экспериментальной базы данных [1, 3]. Для вновь создаваемых вентиляторов необходимо проведение серии экспериментов на специальных стендах по определению поправочных коэффициентов, в том числе и коэффициента полезного действия, а также по выявлению каких-либо недостатков, имеющихся в конструкции создаваемого устройства. С созданием современных вычислительных пакетов [4] появилась возможность проводить прямое численное моделирование (вычислительные эксперименты) уже имеющихся вентиляторов и адекватно предсказывать их характеристики, что ведет к существенному снижению материальных затрат [5]. Важной характеристикой осевого вентилятора является аэродинамический коэффициент полезного действия, представляющий собой отношение мощности потока воздуха, производимого вентилятором, к мощности момента сил действующего на лопатки рабочего колеса (ротора) вентилятора. Очевидно, что этот коэффициент будет зависеть от скорости вращения ротора, формы лопаток, их числа и угла расположения лопаток по отношению к направляющей вентилятора. В настоящей работе демонстрируется метод расчета аэродинамического коэффициента полезного действия на примере одноступенчатого осевого вентилятора со спрямляющим аппаратом с использованием прямого численного моделирования на основе вычислительного пакета Ansys-Fluent. Физическая модель осевого вентилятора Рассматривается осевой вентилятор, рис. 1, состоящий из ротора 1, спрямляющего аппарата (СА) 2, направляющей втулки ротора 3, направляющей втулки СА 4 и корпуса 5. Направляющие втулки ротора и СА имеют одинаковый радиус R1. Внешний корпус представляет собой трубу радиуса R2. На роторе с постоянным шагом расположено Nr рабочих лопаток, высота которых равна R2 - R1. Угол между окружным направлением и хордой профиля лопатки ротора равен а. Длина хорды лопатки ротора равна br, рис. 2, а) На СА с постоянным шагом расположено Ns направляющих лопаток, хорды профилей которых длиной bs параллельны образующей втулки СА, рис. 2, б). Рис. 2. Формы лопаток: а - рабочая лопатка ротора; б - направляющая лопатка СА Лопатки ротора, вращаясь с угловой скоростью ю, обеспечивают перемещение воздуха в направлении параллельном образующей втулки (продольное направление), при этом происходит повышение полного давления воздуха. В роторе поток воздуха закручивается и перемещается в СА, где, благодаря направляющим лопаткам, происходит выравнивание направления движения воздуха, при этом повышается его статическое давление [3, 6, 7]. Размеры осевого вентилятора и скорости вращения ротора принимаются равными соответствующему шахтному вентилятору местного проветривания, поэтому, в силу низких значений числа Маха, пренебрегаем эффектом сжимаемости и полагаем течение изотермичным. Математическая модель осевого вентилятора Поскольку области ротора и СА, свободные для прохода воздуха, имеют периодически повторяющиеся участки (угол периодичности для ротора 2 л/Nr и 2л/Ns для СА), то рассматривалась следующая геометрическая область, представленная на рис 3. Для моделирования течения воздуха через осевой вентилятор решается следующая система уравнений вязкого несжимаемого газа, записанная в общем случае, во вращающейся системе координат. (1) 2 Рис. 3. Геометрическая область моделирования: 1 - лопатка ротора; 2 - лопатка СА; 3 - направляющая ротора; 4 - направляющая СА; 5 - корпус; 6 - входная поверхность ротора; 7 - выходная поверхность ротора; 8 - входная поверхность СА; 9 - выходная поверхность СА Уравнение сохранения массы смеси V-(pVr ) = 0. Уравнение сохранения количества движения смеси V-(pVrV) + pco xV = -Vp + V-Teff, (2) где ur = ю x r - скорость вращения ротора; r - расстояние от оси вращения до рассматриваемой точки. Уравнение переноса турбулентной кинетической энергии V-(pkVr ) = Vц + ^л | Vk + Gk -ps. (3) Уравнение переноса скорости диссипации турбулентной кинетической энергии + cief Gk "c2£p7 k k ц+НошЬ |ve V-(pSVr ) = V (4) к где reff - эффективный тензор напряжений, т + Tturb; Tturb - турбулентный тензор напряжений, цШгЪ (VV + VVT)--pkU; т - тензор вязких напряжений, |(VV + VVT); ц - динамическая вязкость; |turb - турбулентная вязкость, к2 pC| -; Gk - производство турбулентной кинетической энергии вследствие сдвие гового течения (действие градиентов скорости), 2|turb S : S ; S - тензор деформаций, 2 (VV+VVT). Параметры в системе уравнений (1) - (4) задавались следующими: С1е = 1.44; С2е = 1.92; Сц = 0.09; стк = 1.0; сте = 1.3. Граничные условия для системы уравнений (1) - (4): на входе в область полагалось равенство полного давления газа атмосферному; на выходе из области полагалось равенство статическое давление газа атмосферному. На стенках корпуса, ротора и СА задавались условия прилипания: Vr = 0 . Турбулентные параметры потока на стенке и пристенной области определялись следующим образом: dk 0 ОГ - = 0 ; е р = --. dn кур Здесь к = 0.4187 - постоянная Кармана; индекс P - относится к центру пристенной ячейки разностной сетки. Метод решения системы уравнений (1) - (4) Для решения поставленной задачи привлекался вычислительный пакет Ansys Fluent. Система уравнений (1) - (4) решалась численно с использованием метода Патанкара. Конвективные члены уравнений аппроксимировались с помощью противопоточной схемы второго порядка точности. Согласование полей скорости и давления осуществлялось по алгоритму SIMPLE на сдвинутой сетке. Для СА система уравнений (1) - (4) решалась в неподвижной системе координат (ю = 0), а для ротора - во вращающейся с угловой скоростью ю относительно продольной оси системе координат. Пересчет параметров с выходной поверхности ротора на входную поверхность СА осуществлялось на основе концепции «поверхности смешения», реализованной в Ansys-Fluent. Согласно этой концепции, из области ротора в область СА передаются профили осредненных параметров полного давления, направляющих косинусов углов вектора скорости в радиальном, тангенциальном и осевом направлениях, турбулентной кинетической энергии и скорости турбулентной диссипации, которые затем используются как входные граничные условия для области СА. Расчетная область разбивалась на 340 тысяч шестигранных ячеек. Точность расчета контролировалась балансом потоков массы. Число ячеек по высоте лопаток равняется 31. Число ячеек вдоль направляющей втулки статора - 120, вдоль направляющей ротора - 130. Количество ячеек в окружном направлении ротора -46, в окружном направлении статора - 40. Сетка сгущалась вблизи лопаток статора и ротора. Результаты численного моделирования работы осевого вентилятора Расчеты параметров осевого вентилятора проводились при следующих значениях параметров: ц = 1.810-5 Пас; р = 1.225 кг/м3; R1 = 225 мм, R2 = 350 мм, br = 67 мм; bs = 80 мм; L1 = 400 мм; L2 = 100 мм; L3 = 600 мм; Nr = 15, Ns = 12; Скорость вращения ротора варьировалась от 500 об/мин до 5000 об/мин. Величина угла а, определяющего ориентацию рабочей лопатки ротора, варьировалась от 5 до 40° с шагом 5°. На рис. 4 и 5 приведены результаты расчета перепада статического давления в зависимости от частоты вращения ротора и угла положения рабочей лопатки вентилятора. Рис. 4. Зависимость перепада давления от частоты вращения ротора: 1 - а=15°; 2 - 30°; 3 - 40° Рис. 5. Зависимость перепада давления от угла поворота лопатки ротора: 1 - n = 1000 об/мин; 2 - n = 3000 об/мин; 3 - n = 5000 об/мин Видно, что увеличение частоты вращения ротора и увеличение угла поворота лопатки ведет к росту перепада статического давления. Аппроксимация результатов расчетов по методу наименьших квадратов показывает, что перепад давления пропорционален квадрату частоты вращения ротора, т.е. АР ~ п2. Рост угла поворота рабочей лопатки ротора до величины 35° ведет к росту перепада давления по закону АР ~ а1'46, после чего показатель степени становится меньше единицы, что говорит об уменьшении эффективности работы лопаток в качестве нагнетательных элементов вентилятора. На рис. 6 и 7 приведены результаты расчета объемного расхода воздуха в зависимости от частоты вращения ротора и угла положения рабочей лопатки вентилятора. Q, м3/с 1 ~ 20-- "И"2----- Рис. 6. Зависимость объемного расхода воздуха от частоты вращения ротора: 1 - а=15°; 2 - 30°; 3 - 40° 0 1000 2000 3000 4000 п, об/мин Q, м3/с 20 L 15 10 5 0 10 20 30 а, град Рис. 7. Зависимость объемного расхода воздуха от угла поворота лопатки ротора: 1 - п = 1000 об/мин; 2 - п = 3000 об/мин; 3 - п = 5000 об/мин Увеличение частоты вращения ротора ведет к линейному росту объемного расхода воздуха Q ~ а (см. рис. 6). Увеличение угла поворота рабочих лопаток вентилятора в исследуемом диапазоне также увеличивает объемный расход воздуха, создаваемый вентилятором. Следует отметить, что интенсивность прироста объемного расхода воздуха снижается для угла поворота лопаток в 40° (см. рис. 7). Одной из важных характеристик вентилятора является зависимость перепада давления от объемного расхода воздуха. Из рис. 8 видно, что на эту зависимость не влияет угол поворота рабочей лопатки вентилятора. Все множество расчетных данных аппроксимируются формулой АР = 12.117Q2 со средней относительной погрешностью 0.6 %, причем локальная относительная погрешность не превышает 2 %. Данная зависимость близка к той, которая получается из уравнения Бернулли: ар=^ q2- (5) где р - плотность воздуха, а S - площадь проходного сечения вентилятора. АР, Па 4000 Другой важной характеристикой вентилятора является мощность потока, создаваемая ротором. Ее величину можно определить как произведение объемного расхода воздуха на перепад давления: Np = Q-АР . Из рис. 4 - 7 следует, что мощность потока воздуха увеличивается с ростом частоты вращения ротора и с ростом угла поворота рабочей лопатки. При работе вентилятора на рабочие лопатки ротора действует сила сопротивления со стороны воздуха. При этом ротор относительно продольной оси испытывает момент сил. Расчеты показывают, что момент сил, действующий на рабочие лопатки ротора, увеличивается с увеличением частоты вращения и угла поворота рабочей лопатки. Причем момент сил, действующий на лопатки, пропорционален квадрату частоты вращения ротора: M\ ~ п2 и является функцией угла а, представляющей собой полином второй степени. 3000 2000 1000 А 1 О 2 -- +3 О 4 ■ 5 □ 6 А ▲ А А А 7 Д 8 } ■ А ? ifi ь jUt^ 0 5 10 15 Q, м3/с Рис. 8. Зависимость перепада давления, создаваемого в вентиляторе, от объемного расхода воздуха через вентилятор: 1 - а = 5°; 2 - а = 10°; 3 - а = 15°; 4 - а = 20°; 5 - а = 25°; 6 - а = 30°; 7 - а = 35°; 8 - а = 40° Зная полный момент сил, действующий на ротор, можно определить мощность сил сопротивления, которые испытывают лопатки ротора при обеспечении получаемых аэродинамических характеристик: Nv = MxNr 2nn/60 . Поскольку M\ ~ n , то мощность вентилятора пропорциональна кубу частоты вращения ротора. Поскольку мощность сил сопротивления пропорциональна моменту сил, действующего на лопатки ротора, то зависимость мощности вентилятора от угла поворота рабочей лопатки будет такой же, как и для момента сил. Отношение мощности потока воздуха, создаваемого вентилятором, к мощности вентилятора представляет собой эффективность (КПД) вентилятора, П = (Np / Nv) -100% или аэродинамический коэффициент полезного действия. На рис. 9 и 10 приведены результаты расчета КПД, выраженного в %. Видно, что 10 20 30 а, град П, % 0 Рис. 9. Зависимость КПД вентилятора от угла поворота рабочей лопатки: 1 - n = 500 об/мин; 2 - n = 1500 об/мин; 3 - n = 5000 об/мин 60 40 20 X r^-J \-п-b-О-Q э Г "А-1 .-Д-^ 1 [--3 -1 f 4-0-5 1-» 1 Об -1-ф- 4-" 1000 2000 3000 4000 n, об/мин Рис. 10. Зависимость КПД вентилятора от скорости вращения ротора: 1 - а = 5°; 2 - а = 10°; 3 - а = 15°; 4 - а = 20°; 5 - а = 25°; 6 - а = 30° с ростом угла поворота рабочей лопатки КПД вентилятора увеличивается и достигает своего максимального значения при значении угла поворота равного 30°. Дальнейшее увеличение угла поворота снижают КПД, рис. 9. Увеличение скорости вращения ротора также ведет к росту КПД вентилятора, рис. 10. Наиболее сильный его прирост происходит при возрастании скорости вращения ротора от нуля до 2000 об/мин. Дальнейшее увеличение скорости вращения ведет к незначительному изменению коэффициента полезного действия. Сравнение результатов численного моделирования с результатами инженерной методики На рис. 11 показана рабочая лопатка ротора, окружная скорость которой при вращении ротора вокруг оси равна и = aR. Воздух движется относительно лопатки со скоростью с, тогда абсолютная скорость движения воздуха есть w. Проекция этой скорости на продольную ось представляет собой среднерасходную скорость воздуха через вентилятор. В общем случае вектора скоростей w и u не совпадают. Различные варианты диаграммы скоростей для ротора осевого вентилятора при- Рис. 11. Диаграмма скоростей воздуха на рабочей лопатке Предположим, что скорость вращения ротора вентилятора и среднерасходная скорость воздуха такие, что u=w. В этом случае будет иметь место равенство u = и tg^). Тогда объемный расход воздуха через ротор можно вычислить следующим образом: Q = a-R - S - tg (а), (6) где a - угловая скорость вращения ротора; R - средний радиус лопатки, равный (R^ + R^/2; S - площадь, свободная для прохода газов через ротор, п (R,2 - R^) - N (Rc - Rm )5/cos (а); 5 - толщина лопатки равная 5 мм. Для проверки используемого допущения u = w на рис. 12 приведено сравнение зависимости расхода воздуха от угла поворота лопатки для скорости вращения ротора 3000 об/мин и 5000 об/мин, полученных с помощью численного моделирования и по формуле (6). Видно, что для углов поворота лопаток в диапазоне от 20 до 25°, формула (6) хорошо предсказывает значения объемного расхода. Для других значений углов поворота формула (6) дает как заниженные значения объемного расхода (при малых углах поворота лопатки), так и завышенные значения (при больших углах поворота лопатки). Причем расхождение в объемных расходах может достигать 40 %. Расчеты показывают, что отношение объемного расхода, полученное по формуле (6) Qh^, к объемному расходу, полученному при численном моделировании Q^, слабо зависит от скорости вращения ротора (отклонение не превышает 1%) и является только функцией угла поворота лопатки. 0 10 20 30 а,град Рис. 12. Зависимость объемного расхода воздуха от угла поворота рабочей лопатки: 1, 2 - численное моделирование; 3, 4 - расчет по формуле (6); 1, 3 - n = 3000 об/мин; 2, 4 - n = 5000 об/мин Используя средние значения для отношения объемных расходов q = Qh^/Q™, полученных при разных скоростях вращения ротора и аппроксимируя это отношение полиномом третьей степени в зависимости от угла поворота лопатки а в диапазоне от 5 до 40°, получим q(a) = 2.9784-10-05а3 - 2.0001-10-03а2 + 5.7287-10-02а + 3.9542-10-01. (7) Зависимость (7) является индивидуальной характеристикой лопаток ротора и для расчета реального объемного расхода воздуха через вентилятор можно рекомендовать формулу (6) с учетом поправочной функции (7): Q = ®-Д- S-tg(a)/q(а) . (8) Для определения перепада давления на роторе воспользуемся соотношением (5) и (6): .pjWja) . (8) 2 На рис. 13 показано изменение перепада давления в зависимости от угла поворота лопаток, полученное в результате численного моделирования (сплошные кривые) и по формуле (8) (штриховые кривые) для скорости вращения ротора 3000 об/мин и 5000 об/мин. Видно хорошее согласование результатов, полученных по формуле (8), и результатов численного моделирования для а = 15°^20°. Для меньших углов поворота лопаток формула (8) дает заниженные значения перепада давления в два раза, а для больших углов - завышенные значения. 0 10 20 30 а, град Рис. 13. Зависимость перепада давления от угла поворота рабочей лопатки: 1, 2 - численное моделирование; 3, 4 -расчет по формуле (8); 1, 3 - n = 3000 об/мин; 2,4 -n = 5000 об/мин Учитывая тот факт, что отношение перепада давления по формуле (8) АРи.ф, к перепаду давления, полученному в результате численного моделирования АРчм, слабо зависит от скорости вращения ротора (отклонение от среднего не превосходит 3 %) и является функцией угла поворота лопаток а, можно получить следующую аппроксимационную формулу для р=АРи.ф/АРчм, представляющую полином третьей степени от а: р(а) = 6.7832-10-05а3 - 4.0091-10-03а2 + 1.0557-10-01а - 1.4261 10-02. (9) Для расчета реального перепада давления можно рекомендовать следующую формулу: рю2R2tg2 (а) АР = --° v ' . (10) 2 р (а) Таким образом, на основе численного моделирования получены функции согласования д(а) и р(а) для формул зависимости перепада давления и объемного расхода воздуха для одноступенчатого осевого вентилятора от скорости вращения ротора и угла ориентации рабочих лопаток ротора. Предложено определение аэродинамического коэффициента полезного действия осевого вентилятора через отношение мощности создаваемого вентилятором потока воздуха к мощности силы сопротивления, действующей на рабочие лопатки ротора. Показано, что перепад давления в вентиляторе как функция объемного расхода воздуха не зависит от угла ориентации рабочей лопатки ротора.

Скачать электронную версию публикации