Возможность расчетного определения эрозионных зон осевого насоса, работающего на жидкометаллическом свинце

Свобода Дмитрий Геннадьевич; Жарковский Александр Аркадьевич; Иванов Евгений Александрович; Щуцкий Сергей Юрьевич; Дягилев Павел Юрьевич

УДК 621.671.001.024

ВОЗМОЖНОСТЬ РАСЧЕТНОГО ОПРЕДЕЛЕНИЯ ЭРОЗИОННЫХ ЗОН ОСЕВОГО НАСОСА, РАБОТАЮЩЕГО НА ЖИДКОМЕТАЛЛИЧЕСКОМ СВИНЦЕ

1 Санкт-Петербургский политехнический университет Петра Великого 2 АО «Центральное конструкторское бюро машиностроения»

Статья поступила в редакцию 25.03.2016

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

Ключевые слова: осевой насос, проточная часть, кавитация, характеристика, эрозия

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

• необходимость частых разборок проточной части и её отмывки от остатков жидкометаллического свинца;

• значительное время на извлечение модели насоса, связанные со сливом ТЖМТ и охлаждением стенда;

• трудности с сохранением чистоты ТЖМТ;

• значительное время на извлечение модели насоса, связанные со сливом ТЖМТ и охлаждением стенда;

• трудности с сохранением чистоты ТЖМТ;

• большие энергозатраты;

• необходимость обеспечения герметичности контура стенда.

Имеющийся опыт создания отечественных натриевых насосов для реакторных установок типа БН [1] показал, что экспериментальные характеристики насоса на жидкометаллическом натрии существенно не отличались от характеристик, полученных по результатам испытаний на воде. Наблюдался лишь некоторый рост КПД насоса при работе на натрии, что объяснялось меньшей вязкостью натрия. При этом натриевые насосы проектировались на рабочую жидкость

Свобода Дмитрий Геннадьевич, старший преподаватель кафедры «Турбины, гидромашины и авиационные двигатели». E-mail: [email protected]

Жарковский Александр Аркадьевич, доктор технических наук, профессор кафедры «Турбины, гидромашины и авиационные двигатели». E-mail: [email protected] Иванов Евгений Александрович, инженер кафедры «Турбины, гидромашины и авиационные двигатели». E-mail: [email protected]

Щуцкий Сергей Юрьевич, заместитель генерального директора -технический директор. E-mail: [email protected] Дягилев Павел Юрьевич, инженер-конструктор E-mail: [email protected]

«вода» методом гидродинамического моделирования.

Исследования, выполненные нами [2], также показали, что безразмерные экспериментальные эпюры для жидкометаллической эвтектики свинец-висмут при течении в цилиндрической трубе подобны расчетным эпюрам течения на воде, а расчетные напорные характеристики насоса, снятые на воде и свинце после пересчета по подобию совпадают. Отсюда был сделан вывод, что течение в каналах жидкометаллического свинца подчиняется таким же закономерностям, что и ньютоновская жидкость с плотностью и вязкостью жидкометаллического свинца. Разработанные нами проточные части осевых насосов испытывались на воде и жидкометаллическом теплоносителе. Гидравлические Н=0) и мощностные N=(0) характеристики насоса, полученные пересчетом по формулам подобия с воды на ТЖМТ, совпадали с характеристиками насоса, полученными на ТЖМТ.

Кавитационные качества осевых насосов, работающих на ТЖМТ, было решено попробовать отрабатывать на основании численных расчетов на воде, так как имеются сведения о том, что проточная часть с бескавитационной формой проточной части по результатам испытаний на воде, обладает высокими ка-витационными качествами и при работе на жидкоме-таллическом натрии.

Методика кавитационных исследований должна включает следующие этапы:

1) Численный расчет кавитационных срывных характеристик насоса для воды.

2) Получение экспериментальных кавитационных срывных характеристик на водяном стенде.

3) Сопоставление расчетных и экспериментальных характеристик.

4) Анализ возможности возникновения в модельном насосе кавитации и мест возможной кавитационной эрозии.

5) Исследование возможности возникновения кавитационной эрозии на водяном стенде с использованием экспресс-методов анализа на основе лакокрасочных покрытий [3] при обеспечении подобия располагаемого кавитационного запаса и величин относительных скоростей потока в рабочих колесах модельного и натурного насосов.

Численные исследования срывных кавитацион-ных характеристик проводились в программном комплексе ЛшуБ СБХ с использованием двухфазной модели течения: вода - воздух (пар). При СББ расчетах кавитация рассматривается как фазовый переход из жидкого в парообразное состояние при падении

давления в проточной части насоса до величины насыщенного пара. На входе в расчетную область задавалась полная энергия жидкости, на выходе - расход, соответствующий режиму, для которого имелись экспериментальные данные. Величина давления насыщенного пара для температуры составляла рп = 3 574 Па. Для напора и КПД модельной ступени были получены их графические зависимости от кави-тационного запаса на входе при постоянных значениях частоты вращения и подачи. В соответствии с ГОСТ 6134-2007 (ИСО 9906_1999) величина критического кавитационного запаса определялась по 3% падению величины напора.

Экспериментальные исследования моделей с к = 2 70 м м были выполнены на водяном стенде АО ЦКБМ. На рис. 1 представлены сравнительные расчетные и экспериментальные кавитационные характеристики модели осевого насоса, полученные на режиме . Напор ступени при построении кавитационных характеристик задавался в безразмерном виде, за 100% принимался напор, соответствующий экспериментальному кавитационному запасу Д к = 1 2 м. Из представленных на рис. 1 графиков следует, что полученные по результатам численного расчета течения с двухфазной моделью кавитации срывные кавитационные характеристики хорошо согласуются с экспериментальными. Таким образом, численный расчет методами CFD позволяет определять момент срыва параметров осевого насоса и находить значение критического кавитационного запаса Дк кр.

I

а.

Э 10

4 —

кр

образования паровых каверн. На стенде ФТ-4 НГТУ им. Р.Е. Алексеева были проведены экспериментальные кавитационные исследования модели проточной части с DpK = 200 мм, на натурной среде - жидкометалличе-ский свинец. На рис. 3 представлены частные кавита-ционные характеристики, полученные при разных значениях частоты вращения насоса.

Рис. 2. Расчетная срывная кавитационная характеристика модели насоса с Орк = 200 мм

Результаты эксперимента на жидкометалличе-ском свинце (рис. 3) показывают, что в отличие от расчетных кавитационных характеристик на воде (рис. 2), срыва параметров насоса на ТЖМТ не происходит, даже при снижении кавитационного запаса до величины близкой к нулю (абсолютное значение давления на входе в насос составляет рабс = 0,0 2 кгс/см2). Однако проведенная по завершении кавитационных испытаний ревизия проточной части, зафиксировала поверхностный эрозионный износ лопастей рабочего колеса (рис. 4).

Кавитационный запас Ah (м)

-•"Расчет —в—Эксперимент

Рис. 1. Срывные кавитационные характеристики модели насоса с Орк = 2 70 мм

Также были проведены численные исследования срывных кавитационных характеристик для модели с Орк = 200 мм на режиме @ = 1 7 0 м3/ч и частоте вращения п = 1000 об/мин (рис. 2), что соответствуют параметрам при испытаниях на свинцовом стенде. Напор ступени при построении кавитационных характеристик задавался в безразмерном виде, за 100% принимался напор, соответствующий кавитационному запасу м. Из рис. 2 следует, что критический

кавитационный запас модели насоса составил Д к

0 , 9 м. Согласно ГОСТ 6134-2007 (ИСО 9906-1999) величина коэффициента кавитационного запаса выбирается из диапазона Кдоп = 1 — 1, 3 . Тогда допустимый кавитационный запас модели с Орк = 200 мм:

Дкдоп = Кдоп ■ Д ккр = 1 , 2 ■ 0 , 9 = 1 , 0 8 м .

Свинцовый теплоноситель при рабочей температуре 400^600°С имеет давление насыщенного пара свинца в диапазоне Па, что на по-

рядки меньше, чем у воды и исключает возможность

Рис. 3. Экспериментальные частные кавитационные характеристики

Результаты ревизии рабочего колеса показывают, что со стороны входа потока жидкометалличе-ского свинца присутствуют следы «питтинга» (изъязвлений) лопастей рабочего колеса. Питтиннг также наблюдается в периферийной части всех лопастях и в области соединения бандажа РК. Питтинг ярко выражен в месте сопряжения лопастей и конусной втулки рабочего колеса со стороны всасывания. Максимальная глубина каверн питтинга достигает 1,5 мм. Поверхности лопаток рабочего колеса имеют волнисто-губчатый характер, что также характеризует наличие их эрозионного износа. Таким образом, экспериментальные кавитационные исследования насоса, перекачивающего ТЖМТ показывают, что, несмотря на отсутствие паровой кавитации в свинце и срыва его основных

параметров, картина эрозионного разрушения схожа с картиной разрушений при паровой кавитации на воде. Наличие эрозионных разрушений при течении в осевом насосе жидкометаллического свинца можно объяснить тем, что в рабочей среде присутствует растворенный газ (кислород, водород, аргон), что вызывает появление в области низкого давления газовой кавитации. И если газовая кавитация в воде не обладает такой разрушительной силой, как паровая, то в случае течения свинца в область разрушившегося газового пузырька устремляется более тяжелый жидкометалличе-ский свинец и удар получается сильным. Поэтому в лопастном насосе, работающем на ТЖМТ интенсивность кавитационных эрозионных разрушений значительно выше, чем при работе на воде.

Рис. 4. Ревизия РК после кавитационных испытаний

Факт наличия интенсивных разрушений обтекаемых потоком поверхностей при значениях кавита-ционного запаса значительно превышающих его критические значение подтверждается многочисленными исследованиями. С.С. Руднев [4] ввел понятия первого, второго и третьего критических режимов кавитации и показал, что кавитационные каверны появляются уже

на начальной стадии возникновения кавитации. При этом на 1-м критическом режиме наличие каверн практически не сказываются на напоре и мощности насоса. Возникает вопрос о способах отработки геометрии проточной части лопастного насоса, работающего на ТЖМТ, которая предотвратит разрушения от кавитационной эрозии.

Известно, что для бескавитационной работы насоса необходимо, чтобы располагаемый кавитацион-ный запас превышал допускаемый кавитационный запас насоса:

Акр > Акдоп.

При работе лопастного насоса на воде, для его бескавитационной работы без изменения основных параметров, величина допускаемого кавитационного запаса находится через критический кавитационный запас и коэффициент запаса Кдоп = (1 ^ 1,3):

Акдоп = Кдоп ■ Аккр.

Исходя из полученных результатов кавитаци-онных исследований на ТЖМТ, значение коэффициента запаса должно быть таким, чтобы исключить появление кавитационных пузырьков, которые образуются при первом кавитационном режиме. Было решено исследовать области образования паровых каверн на обтекаемых поверхностях проточной части при численном моделировании на воде и связать их с местами возникновения питтинга при испытаниях насоса на ТЖМТ. При численных исследованиях в программном комплексе ЛшуБ СБХ были получены эпюры распределения давления вдоль контура лопасти при разных значениях кавитационного запаса. На рис. 5, 6 представлены эпюры для втулочного и периферийного сечений лопасти (Б - криволинейная координата в меридианной плоскости от входной (Б=0) до выходной (Б=1) кромки). Кавитационные каверны будут образовываться в тех областях лопасти, в которых местное статическое давление падает ниже давления насыщенного пара и где происходит фазовый переход рабочей жидкости из жидкого в газообразное состояние.

-йН=з,26 м (срыв параметров) ---давление парообразования {Рп=3574 па]

Рис. 5. Эпюры распределения давления во втулочном сечении лопасти

-¿11=3,26 м (срыв параметров) -—давление парообразования (Рп=3574 Па)

Рис. 6. Эпюры распределения давления в периферийном сечении лопасти

Анализ рис. 5, 6 показывает, что кавитационные зоны на лопастях рабочего колеса образуются при значениях кавитационного запаса, значительно превышающем критический кавитационный запас. В первую очередь каверны появляются на входной кромке и тыльной стороне лопасти. Причем во втулочном сечении каверна на тыльной стороне лопасти образуется при кавитационном запасе м, а в перифе-

рийном сечении при Д к = 5 , 3 3 м. С уменьшением ка-витационного запаса, каверны увеличиваются в размерах. Срыв параметров насоса происходит, когда кави-тационная каверна захватывает всю лопасть по высоте, каверны появляется на рабочей стороне лопасти. Расчетным путем установлено, что процессы кавитации в рабочем колесе насоса начинаются в периферийных сечениях на тыльной стороне лопасти при Д ккр. Для полного отсутствия кавитационных каверн необходимо обеспечить кавитационный запас Д к > 2 , 5 ■ Д к кр.

Для определения мест возможной кавитацион-ной эрозии была выполнена визуализация объемного распределения пара. В местах образования кавитаци-онных каверн может происходить эрозионное разрушение поверхности лопасти, интенсивность которой зависит от степени развития кавитации, продолжительности ее воздействия, величины средней относительной скорости обтекания, агрессивности рабочей среды. На рис. 7-9 такая визуализация показана для тыльной стороны лопасти рабочего колеса при различных значениях кавитационного запаса.

Рис. 8. Паровые каверны при кавитационном запасе Д к = 3 , 7 7 м

На рис. 8 можно видеть, что при кавитационном запасе Д к = 3 , 7 7 каверна захватывает значительную часть тыльной поверхности лопасти рабочего колеса, что приводит к началу падения энергетических параметров насоса. При значении Д к = 3 , 2 6 м (рис. 9) паровая каверна нарушает кинематику обтекания лопасти потоком, что вызывает срыв параметров насоса.

Papour. Volume Fraction Contour 1

ß»

У

Рис. 7. Паровые каверны при кавитационном запасе Д к = 4, 0 8 м

Рис. 9. Паровые каверны при кавитационном запасе Д к = 3 , 2 6 м

Выводы: степень влияния кавитационных каверн на эрозию в условиях работы на ТЖМТ еще не достаточно изучены. Требуются дальнейшие исследования кавитационных качеств лопастного насоса по методике, предложенной выше. Использование численных методов вычислительной гидродинамики для среды «вода» позволяет значительно упростить и ускорить отработку проточной части осевого насоса.

СПИСОК ЛИТЕРАТУРЫ:

1. Митенков, Ф.М. Главные циркуляционные насосы АЭС / Ф.М. Митенков, Э.Г. Новинский, В.М. Будов. - М., Энерго-атомиздат, 1990. 376 с.

2. Свобода, Д.Г. Проектирование и исследование проточной части главного циркуляционного насоса низкой быстроходности / Д.Г. Свобода, A.A. Жарковский, П.Ю. Дягилев и др. // Тезисы МНТК «Инновационные проекты и техно-

логии ядерной энергетики», 7-10 октября 2014. - М.: ОАО «НИКИЭТ», 2014. С. 218-223.

3. Лунаци, Э.Д. Отечественный метод прогнозирования кавитационной эрозии лопастных насосов с применением легко разрушающихся от кавитации материалов и характеристик их относительной стойкости // Инженерный вестник. №8, август 2015.

4. Руднев, С.С. Кавитация в решетке профилей конечной толщины //Лопастные насосы. - Л.: Машиностроение, 1975. С. 198-209.

POSSIBILITY OF CALCULATIVE DEFINITION OF EROSIVE ZONES IN AXIAL PUMP, WORKING

AT METAL LIQUID LEAD

1 Peter the Great St. Petersburg Polytechnical University

2 JSC Central Design Bureau of Mechanical Engineering

Results of numerical and experimental studies of cavitational qualities of the axial pump, working at working liquid - metal liquid lead, are presented. It is shown that methods of computing hydrodynamics allow to define areas of possible cavitational erosion. Conditions of absence of cavitational erosion are established by calculative way.

Key words: axial pump, flowing part, cavitation, characteristic, erosion

Dmitriy Svoboda, Senior Teacher at the Department "Turbines, Hydromachines and Aircraft Engines". E-mail: [email protected]; Alexander Zharkovskiy, Doctor of Technical Sciences, Professor at the Department "Turbines, Hydromachines and Aircraft Engines". E-mail: [email protected];

Evgeniy Ivanov, Engineer at the Department "Turbines, Hydromachines and Aircraft Engines". E-mail: 20evgeniy08 @maiI.ru;

Shchutskiy, Deputy General Director - Technical Director. E-maiI: [email protected];

Pavel Dyagilev, Engineer-Designer. E-maiI: [email protected]

POSSIBILITY OF CALCULATIVE DEFINITION OF EROSIVE ZONES IN AXIAL PUMP, WORKING AT METAL LIQUID LEAD

Текст научной работы на тему «Возможность расчетного определения эрозионных зон осевого насоса, работающего на жидкометаллическом свинце»