CC BY
479
УДК 621.671.001.024
РАСЧЕТНОЕ ОПРЕДЕЛЕНИЕ КАВИТАЦИОННЫХ ХАРАКТЕРИСТИК ЦЕНТРОБЕЖНЫХ НАСОСОВ
© 2011 В. А. Алексенский, А. А. Жарковский, П.В. Пугачев
Санкт-Петербургский государственный политехнический университет
Поступила в редакцию 12.03.2011
Представлены результаты расчётного исследования кавитационных характеристик центробежного насоса с использованием программного комплекса ANSYS CFX 12 на основе двухфазной модели течения и подхода кафедры гидромашиностроения СПбГПУ, основанный на анализе эпюр распределения давления на лопастях рабочих колес при бескавитационном обтекании. Оба подхода дали качественно близкие результаты и хорошее согласование с экспериментом. Второй подход требует значительно меньших затрат времени на проведение расчетов.
Ключевые слова: центробежный насос, вязкое течение, кавитация, математическая модель, прогнозирование
Вычислительная гидродинамика широко используется для расчёта течения в гидравлических машинах в условиях отсутствия кавитации, однако в связи со значительными трудностями численного и физического характера, возникающими при исследовании кавитации, вычислительная гидродинамика только недавно стала использоваться для расчёта течений с наличием кавитационных явлений в проточной части. Во многих вычислительных моделях, разработанных для исследования кавитации, используется баротропное уравнение состояния для выражения плотности смеси как функции локального давления. Этот подход к моделированию является достаточно привлекательным, т.к. он может быть без значительных трудностей включён в любую CFD-программу обладающую базовой функциональностью, однако подобные подходы предполагают, что термодинамическое состояние является равновесным, т. е. что при изменении условий течения двухфазная среда мгновенно переходит в равновесное термодинамическое состояние. На самом деле эффекты, связанные с конечной скоростью протекания кавитационных про-цесссов, являются важными. Эйсенбург в работе [3] приводит краткое описание динамики переходных процессов, происходящих при кавитации. Начало кавитации в используемых на практике жидкостях связывается с ростом пузырьков размером от 10-5 до 10-3 см.
Алексенский Вадим Александрович, аспирант. E-mail: vadim@aleksensky.com
Жарковский Александр Аркадьевич, доктор технических наук, заведующий кафедрой гидромашиностроения. E-mail: azharkovsky@pef.spbstu.ru Пугачёв Павел Владимирович, аспирант. E-mail: pugachxxx@mail.ru
Эти пузырьки содержат смесь пара и неконденсирующихся газов. Когда пузырек проходит через области, где давление падает ниже давления насыщенных паров, они, как правило, резко увеличиваются в размерах. Область развитой кавитации состоит из большого числа этих пузырьков. При дальнейшем прохождении пузырьков по потоку, когда они попадают в область более высокого давления, происходит их схлопывание. Процесс схлопывания является сложным и зависит от многих факторов, таких как поверхностное натяжение, эффектов связанных с вязкостью и наличия неконденсирующихся составляющих в среде. Именно со значительными силами и давлениями, возникающими на последней стадии схло-пывания, связана основная часть шума и повреждений, сопровождающих кавитацию.
Для получения расчетным способом срывной кавитационной характеристики H=f(Дh) разработчики программного комплекса Ашу8 СБХ 12 предлагают использовать двухфазную модель. Кавитация при этом рассматривается как фазовый переход [1], однако использование двухфазных моделей, как правило, связано со значительным ростом вычислительных затрат. Также, как показала практика, подобные методики не всегда дают хорошую сходимость с экспериментальными даны-ми при расчётах кавитационных характеристик центробежных насосов.
В данной работе предлагается альтернативный путь прогнозирования кавитационных характеристик центробежных насосов, основанный с некоторыми допущениями на по-доб-ной методике, разработанной и используемой для определения антикавитационных качеств
гидравлических турбин [2]. Суть методики состоит в следующем. Запишем выражение для гидравлического КПД РК насоса:
Р Ра
-<- и не сказывается на поле
Пг =
а -Мг
где Мг - крутящий момент относительно оси вращения рабочего колеса, ю - угловая частота вращения РК, Q - подача РК, Н - напор РК, р -плотность перекачиваемой жидкости.
Из представленной формулы следует, что изменение напора РК на фиксированном режиме работы ю=соп81:) может происходить за счет комплекса пГ ' Мг . Прогнозирование уменьшения напора более чем на 3% не имеет смысла, так как при развитой кавитации срыв работы центробежного РК происходит очень резко. Поэтому, если в представленном равенстве пренебречь изменением гидравлического КПД, то получим, прямую зависимость Н от Мг. Величина момента Мг определяется по формуле:
М2 = ^Ар-г-ёБ
кромки,
Р-Я
- удельная энергия.
Пусть, например, при некотором значении Д^! кавитация отсутствует и значение крутящего момента равно Мг1. Далее зададимся значением ДИ2< ДИ!. Эпюры давлений сместятся вниз на величину ДЛ=ДЛ!-ДЛ2. Будем считать, что на участке лопасти, где удельная энергия меньше удельной энергии соответствующей давлению парообразования, наличие кавитационной каверны как бы "срезает" бескавитационную эпюру давлений в зоне, где
Р'Я Р-Я
давлений на других участках поверхности лопасти и крутящий момент составил Мг2. В рамках указанных допущений можно по изложенной выше методике построить срывную характеристику Н=Г(ДЛ) и найти критическое значение кавитационного запаса ДЛкр, при котором произойдет снижение напора, например, на 3%.
где Ар - перепад давлений на элементарной меридиональной проекции площадки ёБ поверхности лопасти.
Из представленной формулы видно, что за счет перераспределения давлений, величина момента Мг будет меняться. Современные программные продукты, в частности Лшу8 СБХ 12 позволяют выполнить расчет пространственного течения вязкой жидкости в проточной части, по результатам которого есть возможность построить эпюры распределения давления вдоль контура лопасти в произвольном сечении. Пример такой эпюры представлен на рис. 1, где Б - криволинейная координата меридианной проекции линии сечения лопасти от входной к выходной
Рис. 1. Эпюра распределения давления вдоль контура лопасти в произвольном сечении
Для определения момента, действующего на лопасть РК целесообразно строить эпюру давлений в зависимости от меридианной длины линии сечения. При этом необходимо знать зависимость изменения радиуса от меридианной длины линии сечения. При изменении граничного условия - давления на входе в РК эпюра будет смещаться вдоль оси ординат на заданную величину. Если эпюра пересечет линию соответствующую давлению парообразования при данной температуре, то можно предполагать наличие кавитационных явлений в данном сечении лопасти. Для описанных методик определения кавитационных качеств насоса проведены верификационные расчетные исследования для центробежного колеса первой ступени питательного насоса с пБ =83, СКР «1200 (при Q«250 м3/ч), отличительной особенностью которого является перерасширенная горловина и, как следствие, наличие противотока на входе на всех режимах. На выходе задавался массовый расход, соответствующий режимам, по которым имелись экспериментальные данные. На рис. 2-7 представлены результаты проведенных исследований для трех режимов работы РК.
Из представленных результатов можно видеть, что срывные характеристики,
полученные на основе методики кафедры гидромашиностроения СПбГПУ [2], хорошо согласуются с экспериментальными данными до момента падения напора на ~ 4% (по результатам расчета течения с двухфазной моделью кавитации КПД РК в данном диапазоне практически не менялся). Что касается результатов, полученных по расчету с двухфазной моделью СБХ, то они хорошо показывают момент полного срыва работы РК. Однако время, затраченное на построение срывной характеристики по результатам расчета с двухфазной моделью СБХ, на порядок превышает время, необходимое для ее построения по методике [2].
Рис. 4. Срывные кавитационные характеристики РК питательного насоса с п^З (0=206 м3/ч)
Рис. 2. Срывные кавитационные характеристики РК питательного насоса с П5=83 (0=253 м3/ч)
Рис. 5. Поле меридианных скоростей в консольном РК питательного насоса с п^83 (0=206 м3/ч)
Рис. 3. Поле меридианных скоростей в консольном РК питательного насоса с nS =83 (0=253 м3/ч)
Рис. 6. Срывные кавитационные характеристики РК питательного насоса с ns=83 (0=94 м3/ч)
Рис. 7. Поле меридианных скоростей в консольном РК питательного насоса с п=83 ^=94 м3/ч)
СПИСОК ЛИТЕРАТУРЫ:
1. Zw art, P.J. A Two-Phase Flow Model for Predicting Cavitation Dynamics, ICMF 2004 International Conference on Multiphase Flow / P.J. Zw art, A.G. Gerber, T. Belamri // Yokohama, Japan. 2004. May 30-June 3.
2. Топаж, Г.И. Расчет интегральных гидравлических показателей гидромашин. -Л.: Издательство ЛГУ, 1989. 208 с.
3. Eisenberg, P. Cavitation. In Handbook of Fluid Mechanics. McGraw Hill, 1961.
CALCULATION DEFINITION OF CENTRIFUGAL PUMPS CAVITATIONAL CHARACTERISTICS
© 2011 V.A. Aleksensky, A.A. Zharkovsky, P.V. Pugachev
St.-Petersburg State Polytechnical University
Results of calculation research the centrifugal pump cavitational characteristics with use of program complex ANSYS CFX 12 on the basis of flow diphasic model and the approach of hydromechanical engineering department of StPSPU, based on the analysis of pressure distributions epures on blades of driving wheels at non-cavitational flow are presented. Both approaches have yielded qualitatively close results and the good coordination with experiment. The second approach demands considerably smaller expenses of time for carrying out the calculations.
Key words: centrifugal pump, viscous current, cavitation, mathematical model, forecasting
Vadim Aleksenskiy, Post-graduate Student. E-mail: vadim@aleksensky.com
Alexander Zharkovskiy, Doctor of Technical Sciences, Head of the Hydromachines Building Department. E-mail: azharkovsky@pef.spbstu.ru Pavel Pugachyov, Post-graduate Student. E-mail: pugachxxx@mail.ru
