CC BY
77
DOI 10.5862/JEST.214.7 УДК 621.224.35
А.В. Семенова, Д.В. Чирков, В.А. Скороспелов
МНОГОЦЕЛЕВОЕ ОПТИМИЗАЦИОННОЕ ПРОЕКТИРОВАНИЕ ФОРМЫ ЛОПАСТИ РАБОЧЕГО КОЛЕСА ПОВОРОТНО-ЛОПАСТНОЙ ГИДРОТУРБИНЫ
A.V. Semenova, D.V. Chirkov, V.A. Skorospelov
MULTI-OBJECTIVE OPTIMIZATION OF KAPLAN RUNNER BLADE SHAPE
В мире широко развиваются методы оптимизационного проектирования формы лопасти рабочих колес гидротурбин, в основу которых положен генетический алгоритм (ГА). Авторы рассмотрели особенности проектирования поворотно-лопастных гидротурбин, которые затрудняют применение методов оптимизационного проектирования, разработанных для радиально-осевых турбин. Представлена методика двухрежимной многоцелевой оптимизации формы лопасти рабочего колеса поворотно-лопастной гидротурбины, направленная одновременно на повышение КПД в оптимуме универсальной характеристики и форсированном режиме и улучшение кавитационных качеств рабочего колеса в наиболее кавитационно-опасном режиме. На основе разработанной методики создан программный модуль проектирования поворотно-лопастной гидротурбины, который включен в программный комплекс оптимизационного проектирования «CADRUN-opt». При его использовании спроектированы рабочие колеса для поворотно-лопастных гидротурбин разной быстроходности.
ПОВОРОТНО-ЛОПАСТНАЯ ГИДРОТУРБИНА; ЛОПАСТЬ РАБОЧЕГО КОЛЕСА; МНОГОЦЕЛЕВАЯ ОПТИМИЗАЦИЯ; ГЕНЕТИЧЕСКИЙ АЛГОРИТМ.
Nowadays runner blade shape optimization methods based on Genetic Algorithm (GA) are being developed throughout the whole world. In the present paper several Kaplan turbine features that complicate the application of the optimization methods as compared to Francis turbines, are being analyzed.The method of multi-point multi-objective optimization of the Kaplan runner blade shape is presented. This method is aimed to increase the efficiency both in the best efficiency point (BEP), and in the full load point with improvement of cavitation characteristics of the runner in the most cavitation-dangerous point. Based on the developed method, the program module for Kaplan turbine design was created. This module was included in the fully integrated software system «CADRUN-opt». Kaplan runner blades of various specific speeds were designed by using this software system.
KAPLAN TURBINE, RUNNER BLADE, MULTI-OBJECTIVE OPTIMIZATION, GENETIC ALGORITHM.
Введение
Развитие методов вычислительной гидродинамики, методов решения оптимизационных задач, достижения в изучении рабочего процесса гидротурбины и повышение производительности компьютеров позволили автоматизировать процесс проектирования лопастной системы рабочего колеса (РК) гидротурбины. Для этого используют метод оптимизационного проектирования. Суть его заключается в пред-
варительной параметризации поверхности лопасти, последующем автоматизированном переборе различных комбинаций геометрических параметров (то есть, форм лопасти), гидродинамическом расчете течения для каждого варианта и выборе такой формы, которая обеспечит минимум одного или нескольких критериев качества. Для поиска, направляющего перебор, как правило, используется эволюционный (генетический) алгоритм оптимизации.
В данной работе представлена методика многоцелевого оптимизационного проектирования лопастной системы рабочего колеса поворотно-лопастной (ПЛ) гидротурбины, положенная в основу модуля программного комплекса «CADRUN-opt», разработанного в ОАО «Силовые машины» — «ЛМЗ» совместно с Институтом математики и Институтом вычислительных технологий СО РАН.
Особенности проектирования поворотно-лопастных гидротурбин
Несмотря на относительно простую форму лопасти рабочего колеса, поворотно-лопастные гидротурбины имеют ряд особенностей, которые затрудняют применение методов оптимизационного проектирования, разработанных для радиально-осевых (РО) турбин [1—3].
Во-первых, при проектировании лопастной системы ПЛ гидротурбины необходимо учитывать, что это машина двойного регулирования: регулирование расхода происходит не только путем изменения открытия направляющего аппарата, но и путем изменения угла установки лопасти. Благодаря этому обеспечивается эффективная работа гидротурбины в широком диапазоне расходов и напоров. Для повышения уровня средневзвешенного КПД лопастной системы проектирование необходимо проводить с учетом анализа работы гидротурбины на нескольких (по крайней мере, двух) режимах работы. Для ПЛ турбины наиболее важны режим оптимального КПД и форсированный режим [4]. Проектирование с учетом форсированного режима работы гидротурбины позволяет также получить требуемые кавитационные качества лопастной системы.
Далее, поворотно-лопастная турбина — быстроходная. Это накладывает при проектировании высокие требования на точность расчета энергетических характеристик, особенно потерь в отсасывающей трубе (ОТ). В настоящее время в программном комплексе «CADRUN-opt» используется комбинированная методика расчета потерь [5]. Суть этой методики в том, что потери в наиболее важных элементах проточной части определяются непосредственно при трехмерном моделировании турбулентного потока с использованием осредненных по Рейнольдсу уравнений Навье — Стокса и соответствующей модели тур-
булентности. Потери в оставшихся элементах проточной части рассчитываются по инженерным формулам расчетно-экспериментальной методики, не требующей проведения гидродинамического расчета в этих элементах. Такой подход позволяет при минимальных затратах вычислительных ресурсов с достаточной точностью прогнозировать потери энергии во всей турбине. Выбор элементов, потери в которых рассчитываются по формулам расчетно-экспериментальной методики, во многом определяется быстроходностью турбины. Так, потери в спирали и статоре хорошо оцениваются по инженерным формулам и для РО и для ПЛ гидротурбин. Кроме того, эти потери не зависят от формы рабочего колеса и могут быть оценены до оптимизационных расчетов. В некоторых работах по оптимизации поворотно-лопастных гидротурбин [6] при оптимизации формы РК течение в ОТ не рассчитывалось: предполагалось, что в точках комбинаторной зависимости потери в ОТ малы и не зависят от формы рабочего колеса. Однако, как показывают исследования, проведенные СКБ «Гидротурбомаш» совместно с институтами СО РАН [7], а также рядом зарубежных фирм [8], для быстроходных турбин расчет течения в отсасывающей трубе необходим для точного прогнозирования КПД турбины и положения оптимума на плоскости (01, я(). Поэтому в настоящей работе при оптимизации формы рабочего колеса ПЛ гидротурбины потери в направляющем аппарате, рабочем колесе и отсасывающей трубе рассчитываются непосредственно с использованием 3D моделирования турбулентного потока вязкой жидкости.
Следующая отличительная особенность ПЛ гидротурбин — наличие зазоров «лопасть — втулка», «лопасть — камера». Эти зазоры существенно влияют на профиль скорости на входе в отсасывающую трубу и, соответственно, потери в ней. С учетом сказанного необходимо рассчитывать энергетические характеристики при проектировании ПЛ гидротурбин с учетом зазоров [7]. Это требование существенно усложняет автоматическое построение расчетной сетки в области рабочего колеса, а также увеличивает затраты вычислительных ресурсов и время оптимизации.
Все перечисленные особенности учтены в излагаемой далее методике оптимизационного проектирования.
Параметризация поверхности лопасти РК ПЛ гидротурбины
Предварительный этап решения оптимизационной задачи — это параметризация геометрии. Параметризация лопастной системы РК должна, с одной стороны, обеспечить гибкое изменение геометрической формы лопастной системы, а с другой, — содержать по возможности меньшее число параметров. В настоящей работе варьируется только форма лопасти; форма втулки, положение оси поворота лопасти, а также число лопастей остаются неизменными. Форма лопасти рабочего колеса ПЛ гидротурбины представляется в виде RЫode(u,v) = r(u,v)+d(u,v)n(u,v),где м,уе[0,1]; (1)
г(м,у) = {Я(и,у), 1(и,у), Ф(и,у)} — серединная поверхность, представленная в цилиндрической системе координат; d(u,v) — функция распределения толщин; п(^ V) — орт нормали к серединной поверхности. Тогда модификация формы лопасти рабочего колеса может осуществляться путем вариации ее серединной поверхности и функции распределения толщин. В данной работе для формы лопасти рабочего колеса поворотно-лопастной гидротурбины реализована только вариация серединной поверхности лопасти, а именно:
путем вариации функции угловой координаты
за счет вариации формы входной и выходной кромок серединной поверхности (вариация В2 проекции).
Для Ф^^) используется «относительная параметризация», т. е. параметризируется и варьи-
руется в ходе оптимизации не сама функция Ф(^V), а отклонение Ф'^^) = —
где — функция угловой координаты ис-
ходной лопасти. Отклонение Ф'^^) задается в виде бикубического полинома от u и V, определяемого значениями 16 свободных параметров [1]. В силу того, что бикубический полином имеет бесконечную гладкость, у модифицированной поверхности гладкость такая же, как у исходной поверхности лопасти. Аналогичная относительная параметризация применена для В2 проекции: ее вариация осуществляется путем изменения значений 8 свободных параметров (см. рис. 1). Таким образом, модификация формы лопасти поворотно-лопастного рабочего колеса осуществляется путем вариации 24 геометрических параметров:
х = (xф, XВ2 = (x1, ■■■, x16, x17, ■■■, х24). (2)
Для каждого геометрического параметра х^ заданы фазовые ограничения — допустимый интервал его вариации [хг- т1п, х^ тах].
Численный метод и граничные условия
В работе [9] отмечено, что от выбора модели течения и постановки гидродинамического расчета при решении оптимизационной задачи зависит выбор целевых функционалов. Так, например, в случае, если отсасывающая труба не включена в расчетную область, в целевых функционалах необходимо косвенно учитывать влияние профиля скорости на выходе из рабочего колеса на потери в ОТ [3] либо подгонять профиль скорости за РК под заданный, как это сделано в [10].
18 А Входная кромка
Входная кромка
24^ Х22 + +
Рис. 1. Параметризация лопасти рабочего колеса
Выходная кромка
В настоящей работе для каждой модификации формы лопасти РК трехмерные расчеты проводятся в стационарной циклической постановке. Расчетная область включает в себя один межлопаточный канал направляющего аппарата (НА), один межлопастной канал РК, области зазоров и отсасывающую трубу. Течение в этой области рассчитывается путем численного решения трехмерных усредненных по Рей-нольдсу уравнений Навье — Стокса, замкнутых стандартной к— е-моделью турбулентности. Численный алгоритм представляет собой неявную конечно-объемную разностную схему, основанную на методе искусственной сжимаемости. Невязкие потоки через грани ячеек рассчитываются с использованием М^^-схемы 3-го порядка аппроксимации, вязкие потоки рассчитываются по центрально-разностной схеме 2-го порядка аппроксимации. Неявная дискретизация уравнений движения и линеаризация полученной схемы приводит к системе алгебраических уравнений большой размерности, которая разрешается методом приближенной LU-факторизации. Стационарное решение находится методом установления по псевдовремени. Численный алгоритм решения уравнений движения, основанный на неявном методе искусственной сжимаемости, подробно описан в [1]. Для передачи данных на границах расчетных областей «НА — РК» и «РК — ОТ» выполняется осреднение всех параметров потока (р, и, V, м>, к, е) по окружному направлению.
Расчет проводится для приведенной турбины с диаметром = 1 м, Н = 1 м без учета силы тяжести. На входе в НА держится угол входа по-
77 Р v2
тока и полная энергия потока Е =--1--, рав-
Р£ 2g
ная
Ещ = Н — ^
(3)
где Н8Р — потери в спирали и статоре, оцененные для каждого из режимов по инженерным формулам. Предполагается, что энергия постоянна по всему входному сечению. Предварительные расчеты в расширенной области, включающей расположенные выше по потоку спиральную камеру и статор, показывают, что это достаточно хорошее приближение. На выходе из ОТ держится энергия ЕоШ = 0 и постоянный градиент давления др / дг = pg. Таким образом, в ходе гидродинамического расчета остается фиксированным напор
на турбине (Н = 1 м), а расход 0 заранее не известен и определяется в ходе решения задачи. Эта постановка удобна, так как в большей степени соответствует реальному физическому процессу и не требует повторного приведения результатов при определении положения режимной точки на плоскости (0(, щ), в отличие от постановки с фиксированным расходом.
Метод оптимизационного поиска
Решение оптимизационной задачи сводится к нахождению минимума (максимума) одной или нескольких заданных целевых функций (или, иначе, функционалов). При проектировании лопастной системы рабочего колеса, как правило, целевые функционалы строят так, чтобы учесть энергетические и кавитационные требования. Вследствие того, что минимизируемые функции независимы, не существует единственного решения, дающего минимум сразу всех целевых функционалов. Решением такой задачи является целое семейство точек, называемое множеством (фронтом) Парето. Фронт Парето представляет собой множество оптимальных решений, каждое из которых лучше остальных возможных, по крайней мере по одному целевому функционалу. Выбор конкретного решения с фронта Парето осуществляется человеком вручную и является отдельной задачей.
Методы оптимизации можно разделить на градиентные методы (использующие понятие производной) и стохастические методы. Стохастические методы лучше подходят для решения многоцелевых задач и позволяют распараллелить процесс поиска, что особенно важно при оптимизации формы проточной части, где для определения целевых функций для каждой модификации формы необходимо проводить трехмерные расчеты течения.
В настоящей работе для нахождения фронта Парето используется многоцелевой генетический алгоритм, подробно описанный в [1]. Идея генетического алгоритма основана на теории эволюции живых организмов. В алгоритме используются аналоги механизмов генетического наследования и естественного отбора. Согласно теории популяция индивидуумов меняется в течение нескольких поколений путем рекомбинации (скрещивания) и мутации отдельных индивидуумов, подчиняясь при этом законам естественного или искусственного от-
бора относительно предписанного критерия. При проектировании лопастной системы рабочего колеса индивидуумом является форма лопасти, определяемая набором параметров к =
= (х1, ..., х24).
Целевые функционалы и ограничения
Турбулентная модель течения и выбранная постановка гидродинамического расчета позволяют использовать в качестве целевого функционала, подлежащего максимизации, КПД гидротурбины
ТЧГ М Ш
Р§ОН
(4)
где М — рассчитанный момент на валу рабочего колеса; ю — угловая скорость вращения рабочего колеса; О — рассчитанный расход через турбину; Н — заданный напор; пм — механический КПД; по — объемный КПД. Для ПЛ гидротурбин полагается, что пм = По = 1.
Требование минимизации потерь или максимизации КПД сводится к повышению КПД на двух режимах работы гидротурбины: оптимальном (01(1), п'(1)) ифорсированном (б1(2),щ'(2)).
Важная задача проектирования — обеспечить требуемые кавитационные качества лопастной системы. Для оценки кавитационных качеств лопастной системы используется подход, основанный на анализе распределения давления на лопасти в некавитационном потоке [11]. При таком подходе минимизируется вклад вращательного момента от области, где давление меньше давления ру парообразования. В настоящей работе используется идейно близкий подход, заключающийся в минимизации взвешенной относительной площади кавитации №гсш, на тыльной стороне лопасти, которая рассчитывается по формуле
xdSy - ydSx
Ж =- ,
<са | xdSy - ydSx '
Ssuc
(5)
Данный целевой функционал введен вместо применявшегося ранее для РО турбин целевого функционала Scav [12]. При его использовании нет необходимости вводить дополнительную весовую функцию, учитывающую местоположение области кавитации на тыльной поверхности лопасти.
Гидродинамические расчеты течения в проточном тракте при решении оптимизационной задачи проводятся при фиксированных значениях открытия а0 направляющего аппарата, угла Л установки лопасти рабочего колеса, частоты щ и полного напора турбины Н = 1 м. Расход находится в процессе решения. В такой постановке расход, пропускаемый через турбину модифицированной лопастной системой, может существенно отличаться от начального расхода, заданного в режимной точке (0[,щ). Корректное сравнение КПД и кавитационных качеств исходной и модифицированной лопастных систем возможно только при условии сохранения заданной режимной точки. Для соблюдения этого условия на каждом режиме вводится ограничение на величину расхода, полученного в результате расчета модифицированной лопастной си-
стемы:
01* -£ <О1 -£2,
(6)
где Scav — площадь области на тыльной стороне лопасти с давлением меньше ру; Ssuc — площадь всей тыльной стороны лопасти.
Наиболее кавитационноопасен форсированный режим, поэтому функционал №гсш, минимизируется только для форсированного режима.
/Т*
где О' — приведенный расход в заданной режимной точке; — приведенный расход, полученный по результатам расчета; е1, е2 — допустимые отклонения. Модифицированные лопастные системы, не удовлетворяющие этому ограничению, исключаются из оптимизации.
После проведения серии оптимизационных расчетов было принято решение о модификации целевых функционалов и введении дополнительных ограничений.
Модификация целевых функционалов
Возможное отклонение расхода модифицированных геометрий от заданного указывает на недостаток прямой максимизации КПД в форсированном режиме [9]. Действительно, в форсированном режиме комбинаторная зависимость ц(0[) имеет существенный наклон. И в ходе максимизации КПД может получиться геометрия, имеющая КПД выше, чем у исходного РК, но с меньшим расходом и лежащая ниже комбинаторной зависимости ц(0[) исходного РК. Эта
I
Рис. 2. Зависимость г|(0 для исходной Рис. 3. Зависимость 1¥сал>(()) для исходной
лопасти, используемая для оценки EffSlope лопасти и используемая для оценки WcavSlope
геометрия обозначена на рис. 2 серым квадратом. Проблема может быть решена уменьшением допустимого отклонения е1, е2, однако при жестких ограничениях велика вероятность вырождения поколения и, вследствие этого, замедления сходимости генетического алгоритма. Поэтому предложен и применен другой подход. Целевой функционал, отвечающий за КПД на форсированном режиме, модифицирован. Фактически максимизируется не КПД, а превышение точки (О', г) над прямой m, касательной к зависимости г(0{) исходной лопасти в точке О'2 (рис. 2). Тангенс угла наклона касательной т оценивается перед проведением оптимизационного расчета из экспериментальной или расчетной зависимости г(О') исходного варианта рабочего колеса. Далее будем обозначать этот целевой функционал EffSlope.
Аналогичная ситуация возникает и при оптимизации с использованием функционала Wcav [9]. В окрестности форсированного режима зависимость Wcav (0() имеет существенный наклон. Поэтому при прямой минимизации Wcav может получиться геометрия, у которой величина Wcav ниже, чем у исходного варианта, но на меньшем расходе (серый квадрат на рис. 3). Чтобы таких геометрий не получалось, используется модифицированный функционал WcavSlope, который максимизирует понижение точки (0, Wcav) от прямой I (см. рис. 3). Наклон прямой I, касательной к зависимости Wcav(0), определяется на основании предварительных расчетов исходной геометрии при различных расходах.
Ограничение на минимальный КПД гидротурбины
Необходимость одновременного улучшения энергетических и кавитационных качеств рабочего колеса предъявляет противоречивые требования к геометрическим параметрам, отвечающим за форму лопасти. Так, например, оптимальная с точки зрения энергетических показателей форма лопасти рабочего колеса может оказаться неприемлемой с точки зрения его ка-витационных качеств [11]. При оптимизации с тремя целевыми функционалами, описанными выше, фронт Парето, как правило, содержит большое число индивидуумов, имеющих хорошие кавитационные качества и низкий уровень КПД (рис. 4, 5). Эти индивидуумы заведомо неприемлемы для проектировщика. Тем не менее как равноправные члены множества Парето они дают потомков при создании нового поколения, большая часть из которых также будет лежать в области заведомо неприемлемых решений.
Для более эффективного использования вычислительных ресурсов и сосредоточения поиска в зоне приемлемых решений предложено ввести ограничение на минимально допустимый уровень КПД:
Г(г') >П ((= и). (7)
Лопасти, не удовлетворяющие этим ограничениям, исключаются из дальнейшей оптимизации. Введение данных ограничений существенно увеличивает количество анализируемых геометрий в области приемлемых решений и тем
П(2)/Пи
Зона приемлемых решений
.г
1
Заведомо неприемлемые решения
П(1)/Пи
Рис. 4. Трехцелевая оптимизация ПЛ 20 с целевыми функционалами ЕЛ1, ЕЛ81орв, WcavSlope
П(2)/Пи
Зона приемлемых решений
— 0.45
— 0.4
— 0.35
Зона неприемлемых решений
0.15
0.1
0.05
0.98 0.99 1 П(1)/Пис
Рис. 5. Трехцелевая оптимизация ПЛ 40 с целевыми функционалами ЕЛ1, ЕЛ81оре, WcavSlope
самым ускоряет процесс поиска формы лопасти рабочего колеса, отвечающей всем поставленным требованиям.
На рис. 6 и 7 представлены фрагменты фронтов Парето в зоне приемлемых решений оптимизационных расчетов для рабочего колеса с напором Ндо 20 м (ПЛ 20) с тремя целевыми функционалами (ЕЛ, EffSlope, WcavSlope), полученные без использования ограничения на минимально допустимый уровень КПД (рис. 6) и с использованием данного ограничения (рис. 7). Рассчитано одинаковое количество поколений. Из рисунков видно, что при решении оптимизационной задачи с использованием ограничения на минимально допустимый уровень КПД удалось получить индивидуумы с луч-
шими кавитационными качествами при сопоставимом уровне КПД.
На рис. 8 и 9 представлены аналогичные результаты для оптимизации рабочего колеса турбины с напором Н до 40 м (ПЛ 40). Из рисунков видно, что при решении оптимизационной задачи с использованием ограничения на минимально допустимый КПД удалось получить индивидуумы с лучшими кавитационными качествами и с более высоким уровнем КПД.
Оптимизация рабочего колеса поворотно-лопастной гидротурбины на напор 20 м
Описанная выше методика применена для оптимизации формы лопасти рабочего колеса ПЛ гидротурбины на напор Н до 20 м с числом
П(2)/Пи
\
0.45 0.4 0.35 0.3 0.25 0.2 0.15
1 1.005 1.01 П^УЛи
Рис. 6. Трехцелевая оптимизация ПЛ 20 без ограничения на минимальный КПД
П(2)/Пи
н
а
а
ЧР '1
Wcav
□ □ □ □ ¿3
0.6
0.45
0.4
0.35
0.3
0.25
0.2
0.15
0.1
0.05
■ н
1 1.005 1.01 П^УПи
Рис. 7. Трехцелевая оптимизация ПЛ 20 с ограничениями на минимально допустимый КПД
Ь А
Исх
■1
1
В
0.45
0.4
0.35
0.3
0.25
0.2
0.15
0.1
0.05
П(1)/Пи
П(2)/Пи
■ ■ Г и
Wcav
0.5
— 0.4 0.35
17_1/}?_И(
П(1)/Пи
Рис. 8. Трехцелевая оптимизация ПЛ 40 без ограничений на КПД
Рис. 9. Трехцелевая оптимизация ПЛ 40 с ограничениями на минимально допустимый КПД
исх
Рис. 10. Расчетная сетка с учетом зазоров «лопасть — втулка», «лопасть-камера»
лопастей рабочего колеса = 5. Расчетная сетка для CFD расчетов, включающая ОТ и зазоры, представлена на рис. 10. Общее число ячеек сетки — около 450 тыс. В качестве исходного взято РК, спроектированное вручную.
Перед решением оптимизационной задачи для выбора режимов проектирования были рассчитаны пропеллерные характеристики исходного РК на двух углах установки лопасти: на угле, близком к оптимальному, и угле, соответствующем форсированному режиму (рис. 11).
Для оптимизации выбраны режимы, соответствующие комбинаторным:
режим 1 — А = 0°, a0 = 26 мм, щ = 140 об/мин (О' ~ 1,11 м3/с);
режим 2 —А = 15°, a0 = 34 мм, щ = 150 об/мин (О ~2,08 м3/с).
Далее проведена двухрежимная трехцелевая оптимизация с помощью программного ком-
плекса «CADRUN-opt». В качестве целевых функционалов задавались:
И (ЕА) — максимизация КПД на режиме 1;
И2 (EffSlope) — максимизация КПД на режиме 2 с учетом зависимости г(О');
В3 (WcavSlope) — минимизация взвешенной относительной площади кавитации Wcav на тыльной стороне лопасти на режиме 2 с учетом зависимости Wcav(Q').
Ограничения на расход: на режиме 1 — О' е [1,06; 1,17], на режиме 2 — О'е [2,03; 2,13]. Ограничения на минимально допустимый уровень КПД: на режиме 1 — п > 0,99попт исх, на режиме 2 — п > 0,95попт исх. Варьировались все 24 геометрических параметра лопасти.
В результате решения оптимизационной задачи были рассчитаны 25 поколений генетического алгоритма по 120 индивидуумов в каждом. На рис. 12 показаны проекции геометрий из мно-
а)
п' , об/мин
160 155 150 145 140 135 130 125
а0=24 а0=26 а0=28
а0=18 а0=20 а0=22
0.8 0.0
12 13 0[, м2/с
б)
п1 , об/мин
160 155 150 145 140 135 130
а0=30 а0=32 а0=34
а0=24 а0=26 а0=28
111111111111111111111
1.4 1.5 1.6 1.7 1.8 1.!
I I I | I | I
2 2.1 2.2 2 3 м2/с
Рис. 11. Пропеллерные характеристики исходной лопасти на оптимальном угле Л = 0° (а)
и на форсированном угле Л = 15° (б)
П(2)/По
д25ч015
□ ©■ □ ■ □ -
рц,в ■ ■
^ ■ Шпя«, ■
п(1)/По,
0,95
Wcav
о;, м2/с
о;, м2/с
Рис. 12. Фронт Парето
жества Парето на плоскости (О', п2) и (О', 2), где «2» указывает на режим 2. Для дальнейшего анализа из фронта Парето проектировщик обычно выбирает несколько геометрий, представляющих собой компромиссные решения. Одна из таких геометрий, g25—i015, отмечена на рис. 12.
Для этой геометрии рассчитаны пропеллерные характеристики. На рис. 13 представлено сравнение разрезов пропеллерных характеристик, рассчитанных для исходной геометрии и для выбранной ^25—Ю15). Разрезы выполнены для п = 140 об/мин и п = 150 об/мин. На рис. 14 представлено сравнение кавитационных качеств исходной и оптимизированной лопастей рабочего колеса.
Из представленных результатов можно заключить, что повышение КПД на обоих режимах составило ~ 0,75 %. Кавитационные качества лопасти также были улучшены.
Важным требованием к проектированию поворотно-лопастного рабочего колеса является сохранение положения оптимума пропеллерной характеристики по частоте щ. Для проверки соблюдения этого требования была рассчитана пропеллерная характеристика выбранного рабочего колеса ^25—Ю15) на оптимальном угле установки исходного рабочего колеса (Л = 0°). Положение оптимума пропеллерной характеристики по оси п1 сохраняется с точностью до 3 об/мин.
П(2)/По:
0,99
0,9В
0,97
0,96
0,94
0,92
Я
пц=140 об/мин
1,6
-Исходная -д25-Ю15
Л=15"
пц=150 об/мин
■
Рис. 13. Пропеллерные характеристики исходной и оптимизированной лопастей на оптимальном и форсированном угле установки лопасти
Рис. 14. Распределение давления по тыльной поверхности исходной лопасти (слева) и оптимизированной лопасти (справа) в форсированном режиме работы. Серым показаны области пониженного давления
Выводы
Разработана методика двухрежимной (режим 1 — оптимальный, режим 2 — форсированный) многоцелевой оптимизации, направленной на повышение КПД в оптимуме и обеспечение пологости универсальной характеристики одновременно с улучшением кавитационных качеств рабочего колеса в наиболее кавитационноопасном режиме. Гидродинамические расчеты проводятся с учетом отсасывающей трубы и зазоров «лопасть — втулка», «лопасть — камера» для корректного вычисления баланса потерь в гидротурбине. При использовании предложенного целевого функционала №Ьау для оценки кавитационных качеств рабочего колеса нет необходимости вводить дополнительную весовую функцию, учитывающую местоположение области кавитации на тыльной поверхности лопасти. Модификация целевых
функционалов для форсированного режима дала возможность ослабить ограничения по расходу, а введение ограничений на минимально допустимый уровень КПД позволили повысить качество поиска в области приемлемых решений.
На основе методики, разработанной Институтом вычислительных технологий СО РАН и Институтом математики СО РАН совместно с ОАО «Силовые машины» — «ЛМЗ», создан программный модуль проектирования ПЛ гидротурбины. Модуль включен в оптимизационный комплекс «CADRUN-opt», внедренный в ОАО «Силовые машины» — «ЛМЗ». С его использованием спроектирован ряд рабочих колес для поворотно-лопастных гидротурбин разной быстроходности. Спроектированные рабочие колеса обладают повышенными энергетическими и кавитационными характеристиками.
СПИСОК ЛИТЕРАТУРЫ
1. Черный С.Г., Чирков Д.В., Лапин В.Н., Скоро-спелов В.А., Шаров С.В. Численное моделирование течений в турбомашинах. Новосибирск: Наука, 2006. 202 с.
2. Cherny S.G., Bannikov D.V., Chirkov D.V., Demi-anov V.A., Pylev I.M., Skorospelov V.A., Stepanov V.N.
Automatic optimal shape design of hydroturbine flow passage // Proc. of Hydro 2008 International Conference. Ljubljana, Slovenia. 2008. P 1-8.
3. Pilev I.M., Sotnikov A.A., Rigin V.E., Semenova A.V., Cherny S.G., Chirkov D.V., Bannikov D.V., Skorospelov V.A. Multiobjective optimal design of runner blade using efficiency and draft tube pulsation criteria // Proc. of the 26th IAHR Symp. 2012. Beij ing, China. P. 1-8.
4. Семенова А.В., Чирков Д.В. Методические вопросы проектирования проточных частей поворотно-лопастных рабочих колес // Сборник докладов 15-й Международной научно-технической конференции студентов и аспирантов. МГТУ им. Н.Э. Баумана. М., 2011. С. 163-166.
5. Sotnikov A., Cherny S., Chirkov D., Bannikov D. Two Approaches to Prediction of Hill Diagram for Francis Turbine // Proceedings of HYDRO VISION RUSSIA 2011 Int. Conf. Moscow, Russia, 27 30 Mach, 2011. CD-ROM. P. 1-15.
6. Lipej A., Poloni C. Design of Kaplan runner using multiobjective genetic algorithm optimization // Journal of Hydraulic Research. 2000. Vol. 38, №1. P.73-79.
7. Семенова А.В., Скороспелов В.А., Чирков Д.В.
Верификация численного прогнозирования энергетических характеристик поворотно -лопастных гидро -турбин // Тяжелое машиностроение. 2014. № 7. С. 29-33.
8. Skotak A., Obrovsky J. Analysis of the flow in the water turbine draft tube in Fluent and CFX // 25th CAD -FEM Users' Meeting. Dresden, Germany. 2007. P. 1-9.
9. Семенова А.В., Чирков Д.В., Лютов А.Е. Целевые функционалы при оптимизации рабочего колеса поворотно-лопастной гидротурбины // Научно-технические ведомости СПбГПУ. 2014. № 3(202). С. 97-106.
10. Skotak A., Obrovsky J. Shape Optimization of a Kaplan Turbine Blade // Proc. of the 23rd IAHR Symp. on Hydraulic Machinery and Systems. Yokohama, Japan. 2006. P. 1-8.
11. Топаж Г.И. Расчет интегральных гидравлических показателей гидромашин. Л.: Изд-во Ленинградского университета, 1989. 208 с.
12. Лобарева И.Ф., Черный С.Г., Чирков Д.В., Скороспелов В.А., Турук П.А. Многоцелевая оптимизация формы лопасти гидротурбины // Вычислительные технологии, 2006. Т. 11, № 5. С. 63-75.
REFERENCES
1. Chernyy S.G., Chirkov D.V., Lapin V.N., Skorospelov V.A., Sharov S.V. Chislennoye modelirovaniye tech-eniy v turbomashinakh. [Numerical simulation of fluid flows in turbomachines]. Novosibirsk: Nauka, 2006. 202 s. (rus.)
2. Cherny S.G., Bannikov D.V., Chirkov D.V., Demi-anov V.A., Pylev I.M., Skorospelov V.A., Stepanov V.N.
Automatic optimal shape design of hydroturbine flow passage. Proc. of Hydro 2008 International Conference. Ljubljana, Slovenia. 2008. P.1-8.
3. Pilev I.M., Sotnikov A.A., Rigin V.E., Semenova A.V., Cherny S.G., Chirkov D.V., Bannikov D.V., Sko-rospelov V.A. Multiobjective optimal design of runner blade using efficiency and draft tube pulsation criteria. Proc. of the 26th IAHR Symp. 2012. Beij ing, China. P. 1-8.
4. Semenova A.V., Chirkov D.V. Metodicheskiye vo-prosy proyektirovaniya protochnykh chastey povorotno-lopastnykh rabochikh koles. [Design problem of Kaplan turbine runner flow passage]. Sbornik dokladov 15-oy Mezhdunarodnoy nauchno-tekhnicheskoy konferentsii stu-dentov i aspirantov. MGTU im. N.E. Baumana. M., 2011.
S. 163-166. (rus.)
5. Sotnikov A., Cherny S., Chirkov D., Bannikov D. Two Approaches to Prediction of Hill Diagram for Fran-
cis Turbine. Proceedings of HYDRO VISION RUSSIA 2011 Int. Conf., Moscow, Russia, 27 30 Mach. 2011. CD-ROM. P. 1-15.
6. Lipej A., Poloni C. Design of Kaplan runner using multiobjective genetic algorithm optimizations. Journal of Hydraulic Research. 2000. Vol. 38, №1. P. 73-79.
7. Semenova A.V., Skorospelov V.A., Chirkov D.V. Verifikatsiya chislennogo prognozirovaniya energet-icheskikh kharakteristik povorotno-lopastnykh gidro-turbin. [Numerical prediction verification of Kaplan turbine efficiency]. Tyazheloye mashinostroyeniye. 2014. № 7. S. 29-33. (rus.)
8. Skotak A., Obrovsky J. Analysis of the flow in the water turbine draft tube in Fluent and CFX. 25th CADFEM Users' Meeting. Dresden, Germany. 2007. P. 1-9.
9. Semenova A.V., Chirkov D.V., Lyutov A.Ye. Tselevyye funktsionaly pri optimizatsii rabochego kolesa povorotno-lopastnoy gidroturbiny [Objective functions used for optimization of Kaplan turbine runner]. Nauch-no-tekhnicheskiye vedomosti SPbGPU. 2014. № 3(202). S. 97-106. (rus.)
10. Skotak A., Obrovsky J. Shape Optimization of a Kaplan Turbine Blade. Proc. of the 23rd IAHR Symp. on Hydraulic Machinery and Systems. Yokohama, Japan. 2006. P. 1-8.
11. Topazh G.I. Raschet integralnykh gidravlicheskikh pokazateley gidromashin [Computation of Integral Hydraulic Characteristics of Hydromachines]. L.: Izd-vo Leningradskogo universiteta, 1989. 208 s. (rus.)
12. Lobareva I.F., Chernyy S.G., Chirkov D.V., Sko-rospelov V.A., Turuk P.A. Mnogotselevaya optimizatsiya formy lopasti gidroturbiny [Multiobjective optimization of hydroturbine blade shape]. Vychislitelnyye tekhnologii. 2006. T. 11, № 5. S. 63-75. (rus.)
СВЕДЕНИЯ ОБ АВТОРАХ
СЕМЕНОВА Александра Владимировна — инженер-конструктор, ОАО «Силовые машины». 195009, Санкт-Петербург, ул. Ватутина, 3А. E-mail: Semenova_AV@lmz.power-m.ru
ЧИРКОВ Денис Владимирович — кандидат физико-математических наук старший научный сотрудник Института вычислительных технологий СО РАН. 630090, Новосибирск, пр. ак. Лаврентьева 6. E-mail: chirkov@ict.nsc.ru
СКОРОСПЕЛОВ Владимир Анатольевич — кандидат технических наук старший научный сотрудник, Института математики им. С.Л. Соболева СО РАН. 630090, г. Новосибирск, ул. Коптюга, 4. E-mail: vskrsp @math. nsc. ru
AUTHORS
SEMENOVA Aleksandra V. — OJSC «Power Machines». 3A, Vatutina str., St.Petersburg, Russia, 195009. E-mail: Semenova_AV@lmz.power-m.ru
CHIRKOV Denis V. — Institute of Computational Technologies SB RAS. Ak. Lavrentjev avenue 6, Novosibirsk, 630090. E-mail: chirkov@ict.nsc.ru
SKOROSPELOV Vladimir A. — Sobolev Institute of Mathematics. 4 Acad. Koptyug Ave., Novosibirsk, 630090, Russia. E-mail: vskrsp@math.nsc.ru
© Санкт-Петербургский политехнический университет Петра Великого, 2015
