УДК 62-137
Пути повышения энергоэффективности динамических насосов на основе современных компьютерных
технологий
© А.И. Петров, В.О. Ломакин, С.Е. Семенов МГТУ им. Н.Э. Баумана, Москва, 105005, Россия
Изложены разработанные на кафедре «Гидравлика, гидромашины и гидропневмоавтоматика» (Э-10) МГТУ им. Н.Э. Баумана методы оптимизации проточных частей динамических насосов. Приведены результаты оптимизации на примере выполненных работ. Кратко описаны основные методы, используемые для повышения энергоэффективности и надежности гидромашин. Рассмотрены результаты верификации расчетных зависимостей путем сравнения с экспериментальными данными.
Ключевые слова: динамические насосы, оптимизация, верификация, вычислительная гидродинамика.
Динамические и лопастные насосы — вид оборудования, совершенствование которого является одной из наиболее приоритетных задач стратегии энергосбережения как в России, так и за рубежом. По некоторым оценкам, до 25 % всей производимой человечеством энергии расходуется на привод лопастных насосов, из чего следует наличие значительного резерва экономии энергии в данной области техники. При этом в стоимости жизненного цикла насосного оборудования составляющая стоимости потребленной насосом за весь срок службы электроэнергии может достигать 80 % и более. Таким образом, необходим поиск возможностей экономии энергии как непосредственно в конструкции насоса, так и в сфере его эксплуатации.
В настоящее время в сфере эксплуатации насосного оборудования имеются значительные положительные сдвиги, связанные с внедрением новых методов регулирования насосов (в том числе частотного регулирования), а также систем дистанционного контроля режимов работы и состояния насосного оборудования. Это позволяет снизить эксплуатационные расходы и повысить долю насосов, работающих в режимах, близких к оптимальным.
Не меньшие резервы повышения энергоэффективности и надежности заложены в конструкции насосов. Большинство выпускаемых в настоящее время насосов были разработаны более 30 лет назад и имеют относительно низкие КПД и ресурс работы. Тем не менее такие насосы
востребованы на рынке, так как у предприятий-изготовителей часто нет средств для существенной модернизации конструкций насосов. Потребители насосного оборудования вынуждены выбирать между дешевым, но неэнергоэффективным отечественным оборудованием и дорогим импортным, которое не всегда удачно интегрируется в существующие системы контроля и ремонта оборудования.
Однако на российском рынке насосного оборудования существуют и позитивные тенденции развития. В связи с увеличением стоимости электроэнергии и ужесточением требований к энергоэффективности производства потребителям оборудования необходимы более совершенные насосы; те отечественные производители, которые вкладывают средства в модернизацию производства и конструкций насосов, более успешно конкурируют на рынке насосного оборудования.
Одно из основных направлений в сфере модернизации конструкций насосов — совершенствование геометрии его проточной части. Поскольку от характера течения жидкости в проточной части насоса зависят как КПД насоса (его энергоэффективность), так и гидродинамические нагрузки на роторе гидромашины (ресурс работы насоса), значительные резервы совершенствования насосного оборудования заложены именно в оптимизации геометрии проточной части насоса.
Метод оптимизации геометрии проточных частей насосов. На кафедре «Гидравлика, гидромашины и гидропневмоавтоматика» (Э-10) МГТУ им. Н.Э. Баумана в течение нескольких лет ведутся работы по модернизации и разработке геометрии проточных частей различных динамических насосов. В частности, были созданы линейки проточных частей нефтяных магистральных насосов типа НМ (шесть типоразмеров), грунтовых насосов большой мощности, струйных аппаратов различного назначения и др. Используемый метод оптимизации геометрии проточных частей насосов или разработки новых проточных частей включает в себя следующие этапы:
• параметризацию геометрии проточной части, позволяющую осуществлять автоматизированное построение как элементов проточной части (подводящее и отводящее устройства, рабочее колесо, предвклю-ченный шнек и т. д.), так и проточной части в целом;
• оптимизацию проточной части, проводимую на основе численного моделирования течений в проточной части насоса; при этом критерии оптимизации и метод оптимизации выбираются в зависимости от постановки задачи;
• комплексное гидродинамическое моделирование течений в проточной части насоса и ее элементах, причем задача может решаться как в стационарной, так и в нестационарной постановке;
• экспериментальную верификацию, основанную на изготовлении макетов насосов с помощью ЭБ-принтера, снятии интегральных пара-
метров макетов (нормальных, кавитационных, балансовых), а также на использовании дифференциального метода с замером давлений в 50—60 характерных точках проточной части насоса.
Выбор изменяемых геометрических параметров проточной части насоса и метод построения его 3D-модели зависит от конкретной задачи. Так, для струйного насоса (рис. 1) значения параметров можно указать и изменять непосредственно при создании 3D-модели в используемой CAD-системе (SolidWorks, CATIA и др.). Однако некоторые параметры (например, пропускная способность расчетного сечения отводящего устройства) не параметризуются в виде конкретных размеров, а сама форма проточной части автоматически не перестраивается в CAD-системе при изменении этого параметра. Поэтому приходится применять специальные системы проектирования проточных систем гидромашин. На кафедре Э-10 используются как специализированное программное обеспечение (программный пакет CFturbo), так и расчетные программы собственной разработки.
Рис. 1. Параметризация струйного насоса
Например, двухзавитковое спиральное отводящее устройство нефтяного магистрального насоса типа НМ имеет сложную форму и интегральные управляющие параметры, которые невозможно описать кон -кретным размером в 3D-модели. Отводящее устройство включает в себя следующие элементы с параметризуемой геометрией: спиральную часть; обводной канал; диффузор. Программный код, разработанный на кафедре Э-10, позволяет строить сечения каждого элемента на основе кривых Безье или 5-сплайнов (рис. 2) [1]. Кривые Безье — параметрические полиномы, координаты Х(и У(точек которых вычисляются по следующему алгоритму:
X (г) ^Рх,—^—г' (1 -1 Г;
1=0 '!(п -')!
У (*) = ¿^'^77——г. г' (1 - г)п ', '=0 '!(п -')!
где п — порядок кривой (порядок кривой на единицу меньше числа полюсов); РХ' и Ру' — координаты '-го полюса; г — параметр, изменяющийся в диапазоне значений 0.. .2.
Кривые Безье — частный случай В-сплайна
С (г) = 1 Щ (г)Р.
'=0
Ь2
а б
Рис. 2. Сечение спиральной части отводящего устройства насоса типа НМ, построенное с помощью кривых Безье (а), и сечение полуспирального подводящего устройства, построенное В-сплайном [2] (б)
Алгоритм поиска координат точек множества сечений отвода по значениям приблизительно 20 геометрических параметров реализован в виде кода программы на языке С++ (рис. 3, а). После построения сечения в виде набора точек оно записывается в файл и загружается в САО-системе как облако точек (рис. 3, б). Затем с помощью встроенных в САО-систему средств из облака точек формируется 3Б-модель проточной части насоса.
Такой алгоритм построения позволяет создавать 3Б-модели проточных частей гидромашин с различным набором значений геометри-
Рис. 3. Интерфейс программы расчета сечений отвода (а) и облако точек, загруженное в САБ-систему (б)
ческих параметров (впоследствии необходимых для оптимизации) в сжатые сроки (несколько десятков моделей в день), что существенно экономит затраты на квалифицированный труд.
Процесс оптимизации подразумевает выбор одного или нескольких критериев оптимизации и метода оптимизации. Критериями могут быть КПД насоса, гидродинамические нагрузки, кавитационные качества, скорость эрозии при перекачке абразивных сред и др. В качестве методов оптимизации можно использовать простой последовательный перебор различных вариантов проточных частей насоса без четко сфор-мулированных критериев качества (на основании интуитивных умозаключений) путем качественного анализа результатов численного мо-делирования гидродинамических процессов; сложный генетический алгоритм оптимизации, требующий существенных затрат вычисли-тельного времени; метод градиентного спуска; симплекс-метод; метод ЬР-т-поиска. На выбор конкретного метода также существенно влияют число критериев оптимизации и имеющиеся в распоряжении разработчика вычислительные мощности.
Следует отметить, что при однокритериальной оптимизации (например, в целях получения максимального КПД насоса) можно применять симплекс-метод или метод градиентного спуска.
В процессе проектирования отводящих устройств насосов типа НМ использовали два критерия: 1) увеличение КПД; 2) уменьшение радиальной гидродинамической нагрузки на ротор. В связи со сложностью составления целевой функции этих критериев применяли метод ЬР-т-поиска, суть которого заключается в следующем [3].
1. Выбор числа пробных точек Ы; в этой работе, исходя из времени расчета течения в проточной части отводящего устройства, было реше-
но осуществлять выборку по 20 точкам в пространстве варьируемых параметров.
2. Вычисление параметра
т = 1 +
1п' 1п2
' = 0, ..., N - 1.
3. Расчет значений точки для значений ц = 0, ..., М-1 (М — число варьируемых параметров)
Чи = Ё" к=1
0,5 Ё ((2'2-2)2г(2)2к-1-2)
2=к
2
1-к
где гцк — значение, взятое из соответствующей матрицы, приведенной в работе [4].
Полученные параметры нормированы и принимают значения 0.1. Чтобы рассчитать параметры уже в заданном диапазоне значений, необходимо вычислить их по формуле Хц = Хтт + (Хтлх + Хть)Чц.
После построения 3Б-моделей проточных частей, соответствующих пробным точкам со сгенерированным набором параметров, проводится численное моделирование в программном пакете БТАЯ-ССМ+. В полученном множестве решений отыскивается наиболее приемлемое решение.
В результате оптимизации параметров отводящих устройств насосов типа НМ удалось снизить потери энергии (указаны потери в спиральном отводе) и максимальную радиальную нагрузку на ротор [4].
Результаты оптимизации параметров отводящих устройств насосов типа НМ
т-т и Радиальная нагрузка
Потери напора Н, м
у у на роторе Рг, Н
НМ 10 000
прототип ..............
оптимизированный ..
НМ 7 000
прототип ..............
оптимизированный ..
НМ 3 600
прототип ..............
оптимизированный ..
13,1
8.74
7.75 5,6
15,34
10,24
8 643 310
926 526
2 345 59,5
При необходимости параметризации геометрии рабочих колес насосов можно использовать специализированный программный пакет СБШгЬо, позволяющий строить 3Б-модель рабочего колеса с большим числом изменяемых параметров (рис. 4).
Рис. 4. Параметризация геометрии рабочего колеса в программном пакете
СТШгЬо
Рассмотрим примеры построения моделей рабочего колеса грунтового насоса и результаты оптимизации геометрии колеса по углу охвата лопасти и углам установки лопастей в условиях жестких ограничений на габаритные размеры рабочего колеса (рис. 5).
В результате изменения угла охвата лопасти и углов установки лопасти гидравлический КПД грунтового насоса мощностью 1 250 кВт удалось увеличить с 74 до 79 %.
Оптимизация элементов проточных частей насосов происходит с использованием гидродинамического моделирования в программном пакете STAR-CCM+ с последующей верификацией результатов расчета на натурных или модельных испытаниях.
При численном моделировании используется метод контрольного объема, суть которого заключается в следующем [5]. Расчетная область разбивается на множество мелких ячеек, размеры которых колеблются в диапазоне значений 0,5.. .10 мм в зависимости от решаемой задачи и геометрического расположения в проточной части (например, внутри рабочего колеса расчетные ячейки измельчаются). Каждая ячейка представляет собой замкнутую область течения жидкости или газа (рис. 6), для которой осуществляется поиск полей макроскопических величин (например, скорости, давления), описывающих состояние среды во времени и удовлетворяющих определенным, сформулированным математически законам. Наиболее используемыми являются законы сохранения, записанные в эйлеровых переменных.
Для любой величины ф в каждой точке O(x, y, z, t) пространства, окруженной некоторым замкнутым конечным объемом, в момент вре-
Рис. 5. Исходное рабочее колесо (а), рабочее колесо с измененными углами установки лопасти на входе (б) и результат гидродинамического моделирования (в)
мени t существует следующая зависимость: общее количество величины ф в объеме может изменяться за счет
• транспорта количества этой величины через поверхность, ограничивающую контрольный объем — поток;
• генерации (уничтожение) некоторого количества величины ф внутри контрольного объема — источники (стоки).
Математически это описывается следующим образом:
Л
(¡X
д(Рф)
дг
У(ри ф) -У( Дфф) = 5с
ф'
о
х
где д (рф) / дг — изменение некоторой физической величины ф; У(ри ф) — конвективное слагаемое в абстрактном
Рис. 6. Ячейка, представляющая законе сохранения физической величи-собой замкнутую область тече- ны ф; У(£>фф) — диффузное слагаемое ния жидкости или газа в абстрактном законе сохранения фи-
зической величины ф; — источниковое слагаемое в абстрактном законе сохранения физической величины ф.
В случае расчета течения несжимаемой жидкости внутри насоса для каждой расчетной ячейки записываются следующие законы сохранения величин.
Уравнение неразрывности для несжимаемой жидкости
=
дУ дУу дУ?
дх ду дг
= 0,
где Ух ,Уу ,Уг — проекции вектора скорости на координатные оси [6]. Здесь и далее горизонтальная черта над величиной означает операцию усреднения величины по времени.
Уравнения сохранения количества движения, усредненные по времени (уравнение Навье — Стокса, усредненное по Рейнольдсу) [6],
дУх + У дУх + У дУх + у дУх . др + 2 д д? дх ду дг дх дх
дУх УУ V—--Р УУ
дх
Л
д_
ду
гдУх + дУу л
V
ду дх )
- рУхУу
д_
дг
V
ГдУх дУл
дх дх
- рУ'У::
дУ - дУ - дУ - дУ др д
+ рУх ^ + рУу ^ + рУ ^ = рЛ+ 2— д? дх ду дх ду ду
г
дх
V
'дУу + дУх "
^ дх ду )
ду
- рУУ
V
дУу
ду
- р УУ
дг
V
'дУу дУ ^ —— +—1
дг ду )
- рУуУХ
дУ + У дУг + У дУг + У дУ др + 2 д
р—+рУ—+рУ^Т+рУ^т = РЛ-т~+
д? дх ду дг дг дг
дУ УУ V—— рУУ
дг
д_
дх
V
гдУ2 дУх
Л
ч дх дг )
- рУУх
д_
ду
V
г дУу +дУ л
\
дг ду )
- р УУ
где р — плотность жидкости (р = 1000 кг/м3); ц — динамический коэф -фициент вязкости жидкости (ц = 0,001 Пах).
Члены вида рУ,У^)■ называются рейнольдсовыми напряжениями. Они появляются в уравнениях в результате формальной процедуры усред-
нения и делают систему уравнений незамкнутой. В результате возникает необходимость вводить дополнительные допущения и взаимосвязи в исходную систему уравнений.
Одним из основных допущений является принятие гипотезы Бус-синеска
(
-Р и и, = 2|
ди, + ди,
дх, дх.
1 дик
Л
3 дхк
-1 Рк§,-
Совокупность новых взаимосвязей и допущений образует модель турбулентности. В этой работе использовалась к-ю-ББТ-модель турбулентности, которая сочетает в себе преимущества к-ю- и к-8-модели. Для пристенной области применена к-ю-модель, а для ядра потока — к-8-модель.
Модель дополняется двумя уравнениями переноса кинетической энергии турбулентности и относительной скорости диссипации этой энергии [7]:
дк
дк
д
— + и,-= Рк-рк ш +-
дг
дх
дх,
, ч дк
(у+°к у т )
дш ТТ дш „2 _ 2 д — + и,-= а52 - Рш2 +-
дг
дх,
дх,
дш
-,л 1 дк дш
+2 (1 - Я )аШ2-——,
ш дх, дх,
где — коэффициент турбулентной вязкости, V т =
а1к
тах (а1ш, 5Р2)
Для численного моделирования используется структурированная призматическая расчетная сетка из пяти слоев вблизи твердых стенок и неструктурированная многогранная в ядре потока (рис. 7). В результате расчета можно получить поля распределения скоростей (рис. 8, а) и давления (рис. 8, б) в проточной части гидравлических машин, а также их рабочие характеристики.
Верификация данных, полученных при математическом моделировании проточных частей насосов, проводится на уменьшенных моделях гидромашин или на натурных образцах. На кафедре Э-10 выполняют следующие виды испытаний динамических насосов (рис. 9):
1) получение нормальных характеристик модельных насосов, а также значений интегральных параметров (напор и момент на валу насоса)
Рис. 7. Расчетная сетка вблизи твердой стенки (а) и сетка, полученная в программном пакете 8ТЛЯ-ССМ+ (б)
Рис. 8. Распределения скорости в сечении насоса типа НМ (а) и давления на границах проточной части насоса (б)
Рис. 9. Общий вид испытательного стенда
для сравнения их со значениями, найденными численным моделированием;
2) замеры давления в максимально большом числе точек проточной части насоса для сопоставления экспериментального распределения давления с расчетным распределением давления;
3) балансовые испытания моделей насосов в целях выделения основных составляющих КПД насоса для уточнения экспериментальных данных.
Для сокращения времени изготовления моделей и снижения массы конструкции вместо металлических литых деталей проточной части используются оболочечные модели
Рис. 10. Оболочечная ф°рма м°- (рис. 10), выполненные на ЗБ-принтере, дели насоса в процессе сборки
м ^ ^ с внешним их усилением для восприя-
тия сил от давления жидкости. При этом рабочее колесо изготовляется в виде полнообъемной модели или отливается из металла. Это позволяет создавать модели корпуса и рабочего колеса насоса непосредственно по ЗБ-моделям проточной части, полученным с помощью программ твердотельного моделирования, не тратя время на изготовление литейной оснастки и последующую обработку отливок, а также относительно легкие конструкции моделей, удобные при сборке.
Отверстия для замера давления выполняются в стенках корпуса в процессе его изготовления на ЗБ-принтере. Давление измеряется в разных точках проточной части насоса с помощью пластиковых трубок, одни концы которых подключены к имеющимся отверстиям для замера давления в различных точках корпуса макета насоса, другие — подсоединены к блоку датчиков.
Созданная система сбора информации позволяет замерять давления в 60 точках одновременно (интервал между замерами 2,5 мкс) с высокой частотой (до 20...100 кГц), что дает возможность регистрировать и пульсационные составляющие давления в различных точках.
Зависимости напора и момента от подачи приведены на рис. 11. Вблизи номинального режима работы насоса погрешность определения напора насоса и момента на валу насоса составляет 1.2 % [8], погрешность расчета давления в проточной части насоса составляет 2.4 %.
Выводы. В настоящее время на кафедре Э-10 создан метод ком -плексного расчета и оптимизации проточных частей динамических насосов, позволяющий рассчитывать параметры центробежных насосов, струйных аппаратов, предвключенных осевых колес и др. Метод разработан с использованием наиболее современных программных пакетов и дает возможность решать широкий спектр задач по проектированию и модернизации динамических насосов. Разработан метод изготовления макетов насосов с использованием ЗБ-принтера, позволяющий изготовлять и испытывать проточные части насосов, точно соответствующие по своей геометрии рассчитанным, причем сроки испытаний зна-
30
40
25
30 щ % 20
н" Р.
я о
20 | X 15
Я
10
10
5
80 120 Подача, м 3/ч а
Рис. 11. Зависимости напора насоса (1, 2) и момента на валу насоса (3, 4) от подачи, полученные экспериментально (1, 3) и численным моделированием (2, 4)
40 80 120 Подача, м3/ч б
чительно сокращаются по сравнению с традиционными методами изготовления опытных образцов и макетов.
Результаты расчетов с помощью такого метода неоднократно подтверждены экспериментально в процессе испытаний макетов и опытных образцов насосов.
ЛИТЕРАТУРА
[1] Ломакин В.О., Щербачев П.В., Тарасов О.И., Покровский П.А., Семенов C.E., Петров А.И. Создание параметризованных ЗБ-моделей проточной части центробежных насосов. Электронное научно-техническое издание «Наука и образование», 2012, № 4. URL: http://technomag.edu.ru/doc/354657. html (дата обращения 20.05.2013).
[2] Ломакин В.О., Петров А.И., Щербачев П.В. Разработка бокового полуспирального подвода с увеличенным моментом скорости жидкости на входе в рабочее колесо. Известия вузов. Сер. Машиностроение, спец. выпуск «Работы студентов и молодых ученых МГТУ им. Н.Э. Баумана», 2012, с. 3—5.
[3] Соболь И.М. Численные методы Монте-Карло. Москва, Наука, 1973, 311 с.
[4] Ломакин В.О., Артемов А.В., Петров А.И. Определение влияния основных геометрических параметров отвода насоса НМ10000-210 на его характеристики. Электронное научно-техническое издание «Наука и образование», 2012, №8. URL: http://technomag.edu.ru/doc/445666.html doi: 10.7463/0812.0445666 (дата обращения 25.04.2013).
[5] Патанкар С. Численные методы решения задач теплообмена и динамики жидкости. Москва, Энергоатомиздат, 1984, 124 с.
[6] Андерсон Д., Таннехил Дж., Плетчер Р. Вычислительная гидромеханика и теплообмен. Т. 1. Москва, Мир, 1990. 384 с.
[7] Белов И.А., Исаев С.А. Моделирование турбулентных течений. Санкт-Петербург: Балтийский государственный технический университет «Воен-мех», 2001, 109 с.
[8] Ломакин В.О., Петров А.И. Верификация результатов расчета в пакете гидродинамического моделирования STAR-CCM+ проточной части центробежного насоса АХ 50-32-200. Известия вузов. Сер. Машиностроение, спец.
выпуск «Работы студентов и молодых ученых МГТУ им. Н.Э. Баумана», 2012, с. 6—9.
Статья поступила в редакцию 08.06.2013
Ссылку на эту статью просим оформлять следующим образом:
Петров А.И., Ломакин В.О., Семенов С.Е. Пути повышения энергоэффективности динамических насосов на основе современных компьютерных технологий.
Инженерный журнал: наука и инновации, 2013, вып. 4. URL: http://engjournal.ru/ catalog/machin/hydro/689.html
Петров Алексей Игоревич родился в 1973 г., окончил МГТУ им. Н.Э. Баумана в 1996 г. Канд. техн. наук, доцент кафедры «Гидромеханика, гидромашины и гидропневмоавтоматика» МГТУ им. Н.Э. Баумана. Автор 37 научных работ, в том числе трех учебно-методических в области гидродинамики. e-mail: alex_i_petrov@mail.ru
Ломакин Владимир Олегович родился в 1985 г., окончил МГТУ им. Н.Э. Баумана в 2009 г. Канд. техн. наук, ассистент кафедры «Гидромеханика, гидромашины и гидропневмоавтоматика». Автор 10 научных трудов в области гидромашин, вычислительной гидродинамики. e-mail: lomakin_vladimir@list.ru
Семенов Станислав Евгеньевич родился в 1966 г., окончил МГТУ им. Н.Э. Баумана в 1989 г. Канд. техн. наук, доцент кафедры «Гидромеханика, гидромашины и гидропневмоавтоматика» МГТУ им. Н.Э. Баумана. Автор более 50 научных работ в области гидравлики и гидродинамики. e-mail: a1e10@rambler.ru