Исследование структуры потока и прогнозирование характеристик секционного центробежного насоса низкой быстроходности

С использованием построенных геометрических моделей в ядре потока создавались неструктурированные расчётные сетки на основе тетраэдальных элементов с их сгущением у входных и выходных кромок РК и на поверхности интерфейса. Вблизи твердых стенок для описания течения в пограничном слое были созданы слои призматических ячеек. Общее количество элементов составило около 8 млн. Расчёты выполнялись в квазистационарной постановке. Параметры потока на поверхностях сопряжения элементов проточной части передавались через поверхность-интерфейс типа «stage» (т.е. проводилось осреднение параметров по окружному направлению).

Рис. 1. Геометрическая модель расчётной области

Были использованы следующие граничные условия: на входе в рабочее колесо задавалось полное давление (1 атм.); на выходе диффузора спирального отвода – массовый расход, соответствующий рассчитываемому режиму (вариант граничных условий, обеспечивающий наилучшую стабильность расчёта); задавалась также частота вращения рабочего колеса. Для замыкания уравнений На-вье-Стокса, осреднённых по Рейнольдсу, использовалась k-e модель турбулентности [1]. Моделирование кавитационных явлений не производилось. Во всей области задавалась постоянная плотность среды. Использовался программный комплекс ANSYS CFX 12.1. В процессе проведения расчёта контролировались не только среднеквадратичные и максимальные значения невязок, но также и целевые параметры – моменты на рабочих колёсах, общий напор и напор каждой из ступеней. Расчет проводился до прекращения существенных изменений каждого из параметров. Механический и объёмный КПД вычислялись по эмпирическим методикам [2]. Расчёт был выполнен на трех режимах, для которых имелись экспериментальные данные. Результаты сопоставления расчетных и экспериментальных характеристик представлены на рис. 2.

Рис. 2. Сопоставление расчётных и экспериментальных характеристик исходного варианта проточной части

Сопоставление показывает, что имеет место хорошая сходимость результатов расчёта с экспериментом на номинальном режиме (расхождение по напору – 0,25%, по КПД – 1%). На других режимах был получен несколько худший уровень сходимости. Визуализация результатов расчёта течения на номинальном режиме в рабочем колесе, направляющем аппарате и обратных каналах представлена на рис. 3.

Рис. 3. Визуализация течения в плоскости z = const: а – рабочее колесо и направляющий аппарат первой ступени; б – обратные каналы

На рис. 3,а можно наблюдать вихревое течение в направляющем аппарате при переходе к зоне поворота потока, также видно, что течение в межлопастном канале рабочего колеса имеет неравномерный характер с низкоэнергетической зоной потока у задней стороны лопасти. Течение в обратном направляющем аппарате (рис. 3,б) имеет более равномерный характер, вихревые структуры появляются только в районе выходных кромок. На рис. 4,а представлена визуализация течения в рабочем колесе в плоскости конформного отображения (вблизи покрывающего диска). Наблюдается заметная вихревая зона в межлопастном канале. На рис. 4,б можно наблюдать вихревую зоны в меридиональной проекции с осреднением параметров по окружному направлению (она выглядит как область низких скоростей у покрывающего диска). Всё это позволяет сделать вывод о том, что возможно увеличить эффективность данного насоса, изменив оптимальным образом параметры его проточной части.

Рис. 4. Визуализация течения в рабочем колесе: а – на поверхности тока в плоскости конформного отображения; б – в меридиональной проекции

С использованием методики проектирования кафедры гидромашиностроения СПбГПУ, основывающейся на квазитрёхмер-ном методе расчёта невязкого течения и пространственном пограничном слое и низкоэнергетическом следе для вязкой жидкости, был спроектирован новый вариант проточной части рассматриваемого насоса. Число лопастей рабочего колеса было увеличено до 9, уменьшена его ширина, вместо канального отвода был применен отвод лопаточного типа. Расчет течения вязкой жидкости был произведен по такой же схеме, как и для исходного варианта насоса. Геометрическая модель проточной части представлена на рис. 5, результаты сопоставления характеристик исходного и модернизированного вариантов – на рис. 6. Из расчетных характеристик можно видеть, что модернизированный вариант проточной части имеет более высокий КПД, чем исходный. Особенно заметное различие достигается на режимах малых подач. Форма напорной характеристики не имеет существенных отличий.

Визуализация результатов расчета течения в новом варианте проточной части: в рабочем колесе и направляющем аппарате первой ступени представлена на рис. 7. Можно видеть, что течение в каналах нового варианта рабочего колеса имеет более равномерный характер – без ярко выраженных низкоэнергетических зон, также отсутствует вихревая зона в лопаточном диффузоре при приближении к зоне поворота потока в обратный направляющий аппарат.

Рис. 5. Геометрическая модель расчётной области для исследования модернизированного варианта проточной части

Рис. 6. Сопоставление характеристик исходного и модернизированного варианта проточной части

Рис. 7. Визуализация течения в плоскости z=const в рабочем колесе и в направляющем аппарате первой ступени (модернизированный вариант)

что после проведения модификации потери в рабочем колесе снизились в большей степени, чем в отводе. Для промежуточных ступеней был получен аналогичный результат.

Вывод: использование в гидромашиностроении численных методов моделирования трёхмерного течения представляет значительный практический и научный интерес, т.к. позволяет не только предсказывать интегральные параметры проточных частей, но и даёт возможности исследовать влияние разных геометрических размеров и форм их элементов на структуру потока. При этом (по сравнению с использованием физического эксперимента) достигается значительное снижение расходов на разработку.

Визуализация показывает, что в результате проведённой модификации проточной части удалось добиться устранения вихря в межлопастном канале. Поэлементный анализ гидравлических потерь в первой ступени показал,