Многоцелевая оптимизация формы лопасти гидротурбины Текст научной статьи по специальности «Физика»

Вычислительные технологии

Том 11, № 5, 2006

МНОГОЦЕЛЕВАЯ ОПТИМИЗАЦИЯ ФОРМЫ ЛОПАСТИ ГИДРОТУРБИНЫ*

И. Ф. Лобарева, С. Г. Черный, Д. В. Чирков Институт вычислительных технологий СО РАН, Новосибирск, Россия

e-mail: cher@ict.nsc.ru

В. А. СкороспЕлов, П. А. Турук Институт математики СО РАН, Новосибирск, Россия e-mail: vskrsp@math.nsc.ru

Automatic runner blade optimization system proposed in our earlier work is based on the Genetic Algorithm for searching an extremum of the specific functional. This system is generalized by including an algorithm which enhance several blade performance criterions at once. Results of calculation of test optimization problems as well as an automatic optimization of the runner of Bratskaya hydro power plant are presented.

Введение

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

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

* Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (грант № 04-01-000246).

© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2006.

позволяющих оценить потери энергии, оказывается в связи с этим чрезвычайно дорогостоящим по затратам счетного времени. В то же время расчет потока в рабочем колесе в приближении идеальной жидкости с использованием соответствующих методов (см., например, [2, 3]) требует гораздо меньших затрат машинного времени. Расчет потока на сетке 40 х 20 х 20 занимает менее 3 мин на процессоре Pentium IV, 3 ГГц, поэтому в работе [1] моделирование течения в РК проводилось в рамках уравнений Эйлера идеальной жидкости. Явный расчет потерь энергии в этой постановке невозможен, поэтому были сформулированы несколько целевых функционалов, позволяющих косвенно учитывать различные виды потерь, оставаясь в рамках модели идеальной жидкости. Существенным недостатком описанного подхода является возможность оптимизировать лопасть лишь по одному из рассмотренных функционалов. Другие функционалы либо не принимаются во внимание, либо их величина контролируется в виде ограничений оптимизационной задачи. В связи с этим хотелось бы иметь алгоритм, позволяющий улучшать сразу несколько критериев качества (целевых функционалов). Необходимо отметить, что использование в процессе оптимизации трудоемких турбулентных расчетов позволяет рассчитать потери в явном виде и тем самым сократить количество функционалов. Однако и в этом случае желательно учитывать, например, кавитационные характеристики рабочего колеса, взаимодействие РК с отсасывающей трубой, а также эффективность работы турбины в других режимах.

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

В настоящее время большое распространение получили методы многоцелевой оптимизации, позволяющие одновременно минимизировать несколько целевых функционалов [4-9]. Вследствие того, что минимизируемые функционалы независимы, вообще говоря, не существует одного решения, дающего минимум сразу всех целевых функционалов. Решением такой задачи является целое семейство точек, называемое множеством (фронтом) Парето. Каждая точка фронта Парето — оптимальна в том смысле, что путем ее возмущения невозможно улучшить какой-либо целевой функционал без ухудшения остальных критериев качества. Таким образом, фронт Парето представляет собой множество оптимальных решений, каждое из которых лучше любого другого по крайней мере по одному целевому функционалу. Компромиссное решение (золотая середина) на фронте Парето выбирается человеком вручную исходя из вида поверхности фронта в пространстве функционалов и предпочтений, отдаваемых какому-либо функционалу. Во всех указанных выше работах по многоцелевой оптимизации для поиска фронта Парето используются варианты генетического алгоритма.

Генетические алгоритмы многоцелевой оптимизации начинают широко применять в задачах автоматического проектирования турбомашин [10-13]. Работа [10] посвящена оптимизации формы рабочего колеса поворотно-лопастной, [11, 12] — радиально-осевой гидротурбины. В [13] решается задача оптимизации сопловой лопатки паровой турбины. В перечисленных выше работах гидродинамические расчеты перебираемых конфигураций проточного тракта проводятся в рамках уравнений Навье — Стокса несжимаемой жидкости, а для моделирования турбулентности используется k — е-модель. Гидродинамические расчеты в [10, 12, 13] выполняются с использованием коммерческих пакетов программ

СЕХ-ТАБС/!от и СЕХ-5, оптимизация проводится по двум критериям. В [11] минимизируются потери в рабочем колесе для трех режимов: оптимального, максимальной загрузки, неполной загрузки, т. е. минимизируются три функционала. Оптимизационные расчеты в [11, 13] производились с использованием распараллеленного генетического алгоритма.

В настоящей работе предложен генетический алгоритм многоцелевой оптимизации, основанный на работах Фонсеки — Флеминга [4] и Хорна [5]. Работоспособность построенного алгоритма продемонстрирована на ряде ставших стандартными тестовых задач с числом варьируемых переменных от 2 до 30 и двумя или тремя целевыми функционалами [6-9]. Алгоритм применялся для многоцелевой оптимизации формы лопасти рабочего колеса Братской ГЭС. Как ив [1], расчеты течения в межлопастном канале РК проводились в рамках уравнений Эйлера. В качестве целевых функционалов рассматриваются кинетическая энергия на выходе из РК, размер области кавитации, отклонение потока в межлопастном канале от "осесимметричного" потока. Приводятся результаты оптимизационных расчетов, направленных на минимизацию двух и трех функционалов.

1. Описание алгоритма многоцелевой оптимизации

Задача многоцелевой оптимизации формулируется следующим образом.

Найти значения параметров х1;..., хN, обеспечивающих минимальные значения д функционалов:

шт (/1(х),... (х)), х =(х1,...,хм) € X С Еи, (1)

при наличии фазовых ограничений

X = {х : хь,1 < х* < хД;*} (2)

и целевых ограничений

^(х) < 0, э = 1,...,т. (3)

В данной работе для решения задачи (1)-(3) используется генетический алгоритм (ГА). В основу ГА положена теория эволюции, согласно которой популяция индивидуумов меняется в течение нескольких поколений путем рекомбинации (скрещивания) и мутации отдельных индивидуумов, подчиняясь при этом законам естественного или искусственного отбора относительно предписанного критерия. Применительно к проектированию лопасти рабочего колеса гидротурбины индивидуумом является рабочее колесо с лопастью какой-то конкретной формы. Индивидуумом мы также будем называть набор параметров х = (х1,... ,хм), определяющих эту лопасть.

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

Следующий этап работы генетического алгоритма — отбор. На этом этапе выбирается определенная доля Тг (0 < Тг < 1) наилучших индивидуумов для скрещивания. В случае одноцелевой оптимизации (д =1) лучшим будет индивидуум, которому соответствует наименьшее по популяции значение целевого функционала. То есть критерием качества

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

гапк(хг) = 1 + аг,

где аг — число индивидуумов текущей популяции, лучших по всем функционалам индивидуума хг [4]. Чем больше ранг, тем больше расстояние до фронта Парето. Отметим, что при решении задачи (1) индивидуумы должны удовлетворять ограничениям (3), однако при создании популяции не все индивидуумы им удовлетворяют. Ранг таким индивидуумам назначается в зависимости от того, насколько сильно не выполняются ограничения. Ранг лучшего из группы решений, для которых не выполняется хотя бы одно из ограничений, больше, чем самый высокий из группы индивидуумов, для которых ограничения выполнены.

Затем для каждого индивидуума популяции хг вычисляется критерий качества, основанный на ранге:

rank(x1)-1

Fr (хг) = 1+ £ n(k), k= 1

где n(k) — число индивидуумов ранга k. После этого проводится процедура, предложенная Голдбергом и Ричардсоном [14]. Цель этой процедуры — распределить популяцию около некоторого числа различных пиков в пространстве поиска, причем возле каждого пика сосредоточивается часть популяции, пропорциональная высоте пика [5]. Иначе говоря, процедура проводится во избежание преждевременной сходимости или сходимости к какому-либо подмножеству точек фронта Парето, а также для получения равномерного распределения точек вдоль фронта Парето. Для этого критерий качества индивидуума, основанный на ранге Fr, умножается на нишевое число т, вычисляемое для каждого индивидуума. Нишевое число тг — это оценка того, как много индивидуумов популяции располагается в окрестности (нише) индивидуума i. Это число рассчитывается по всей популяции, включая i-й индивидуум:

тг = Sh(P(i,j ))-

jePop

Здесь p(i,j) — расстояние между индивидуумами хг и xj, а Sh(p) — функция разделения такая, что Sh(0) = 1, Sh(p > ashare) = 0. В данной работе использовалась

( 1 — P/ashare, Р < ashare

\ 0, Р > ashare,

где ashare — параметр разделения, который подбирается эмпирически. Индивидуумы, находящиеся на расстоянии ashare, ухудшают (увеличивают) функции качества Fs = Frm друг друга.

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

Sh(p)

отрезках разной длины, функционалы предварительно нормируются таким образом, чтобы они отображали область поиска X в гиперкуб с единичным ребром [0,1]q. Расстояние между индивидуумами вычисляется в пространстве функционалов по формуле

Р(«>j) = ,max !/fc(x') - fk(x)|-

k=1,...,q

Отбор индивидуумов проводится исходя из значения функции качества Fs. Затем выполняются рекомбинация, мутация, клонирование. В результате получается новое поколение, заменяющее предыдущее, и этот процесс повторяется до тех пор, пока не будет рассчитано заданное число поколений NQen.

Процесс поиска решения включает, таким образом, следующие шаги.

1. Формирование начальной популяции, состоящей из p0 индивидуумов. Каждый индивидуум — это набор параметров x = (xi,... , Xn), определяющий объект оптимизации, при этом XL,i < Xi < . Начальная популяция создается случайным образом применительно ко всем параметрам.

2. Вычисление значений целевых функционалов /i,..., /q для для каждого индивидуума популяции.

3. Вычисление ранга rank каждого индивидуума популяции.

4. Вычисление критерия качества, основанного на ранге Fr для каждого индивидуума.

5. Вычисление расстояний р между индивидуумами.

6. Вычисление нишевых чисел m для каждого индивидуума.

7. Вычисление функции качества Fs.

8. Отбор. На этом этапе выбирается определенная доля Tr (0 < Tr < 1) наилучших индивидуумов для скрещивания. Это будут индивидуумы, которым соответствуют наименьшие значения Fs.

9. Построение новой популяции. На данном этапе строится новая популяция, состоящая из p индивидуумов, путем рекомбинации. После чего каждый из индивидуумов слегка меняется в процессе мутации [1].

10. Клонирование. Копирование без изменения всех индивидуумов, удовлетворяющих ограничениям, и ранга 1 из предыдущего поколения в новое поколение. После клонирования размер нового поколения будет равен p + pe, где pe — число клонированных индивидуумов.

11. Далее осуществляется переход на шаг 2, пока не будет рассчитано заданное число поколений Noen.

2. Тестовые расчеты

Работоспособность алгоритма проверялась на следующих тестовых задачах (здесь приводятся общепринятые названия): [6], [7], ТХК [8], 05У [9]. В задаче число минимизируемых функционалов д равнялось трем, в остальных тестовых задачах — двум; число варьируемых параметров п в равнялось 30, в — 10, в — 12, в ТХК — 2 и в 05У — 6. В задачах ТХК и 05У накладывались ограничения. В табл. 1 приведены функционалы, область изменения параметров и ограничения для упомянутых выше тестовых задач.

Критериями качества алгоритма, осуществляющего многоцелевую оптимизацию, являются:

— минимизация расстояния между полученным с помощью алгоритма множеством решений и действительным фронтом Парето;

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

Таблица 1. Тестовые задачи многоцелевой оптимизации

Задача

Целевые функционалы

Ограничения и число параметров

Д(х) = хь

/2(х) = д (1 - 77175 - (/1/5) в1п(10пД)

п

д(Ж2, . . . , Жп) = 1 + Хг/(п - 1)

0 < Хг < 1, п = 30

г=2

/1(х) = 1 — exp(—4х1) 8ш6(6пх1 ), /2(х)= д (1 - (/1/д)2),

(п

^Хг/(п - 1)

г=2

0 < Хг < 1, П = 10

0.25

/1(х) = (1 + д) есэ

/2(х) = (1 + д) есэ

/з(х) = (1 + д) эш

пх;

пх\

пхт

есэ

эт

пхо

пхп

0 < хг < 1,

п = 12

д(х3,..., хп) = ^(хг - 0.5)2, а = 100

г=з

/1 (х) = х1, /2(х) = х2

0 < хг < п,

х1

^1(х) = 1 + 0.1 есэ 16 аг^ап — -

V х2У

-х21 - х22 < 0, ^2(х) = (х1 - 0.5)2+

+(х2 - 0.5)2 - 0.5 < 0,

п=2

/1(х) = -(25(х1 - 2)2 + (х2 - 2)2+

+ (хз - 1)2 + (х4 - 4)2 + (х5 - 1)2),

/2(х) = ^г

г=1

0 < х1, х2, х6 < 10,

1 < х3, х5 < 5,

0 < х4 < 6,

^1(х) = 2 - х1 - х2 < 0, ^2(х) = х1 + х2 - 6 < 0, ^з(х) = -х1 + х2 - 2 < 0, ^>4(х) = х1 - 3х2 - 2 < 0, ^5(х) = (хз - 3)2 + х4 - 4 < 0, Рб(х) = 4 - (х5 - 3)2 - х« < 0, п = 6

2

2

2

2

2

2

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

Рис. 1. Задача ZDTз. Рис. 2. Задача ZDT6.

Рис. 3. Задача Т^К. Рис. 4. Задача ОБУ.

При решении всех тестовых задач использовались следующие значения параметров алгоритма: параметр разделения а8ьаге = 0.01, доля отбираемых для рекомбинации индивидуумов Тг = 0.3, параметр рекомбинации й = 0.7, мутационный параметр д = 0.1. Число индивидуумов в начальной популяции р0 ив последующих р бралось одинаковым. В задачах без ограничений p0 = р = 70, с ограничениями p0 = p =100. Число поколений в задачах равнялось 150, в Т^К ^сеп = 50, в остальных 100. На рис. 1-5

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

3. Оптимизация формы лопасти

В качестве объекта многоцелевой оптимизации выбрано радиально-осевое рабочее колесо гидротурбины Братской ГЭС. Для параметризации формы лопасти использовался

подход с представлением срединной поверхности лопасти в виде бикубического полинома с 16 свободными параметрами [1]. Варьировалась только угловая координата Ф(и,у), определяющая кривизну лопасти. Д^-проекции лопасти, проточного тракта и распределение толщины оставались неизменными. Для каждого набора геометрических параметров 3Б-расчеты поля течения в полученном рабочем колесе выполнялись в рамках модели невязкой несжимаемой жидкости с использованием алгоритма [2, 3].

Расчеты течения проводились в одном межлопастном канале рабочего колеса при предположении, что течения в остальных межлопастных каналах циклически повторяются (рис. 6). Расчеты всех конфигураций выполнялись для заданного оптимального режима работы гидротурбины (при фиксированных расходе Q и частоте вращения п). Для оценки качества лопасти использовались целевые функционалы, значения которых вычислялись на основе рассчитанного 3Б-поля течения. Их описание приводится ниже.

Рис. 5. Задача DTLZ4.

Рис. 6. Расчетная область и сетка в ней.

Рис. 7. Вид весовой функции для ограничения по кавитации.

3.1. Ограничения

В процессе оптимизации ставились ограничение на напор в рабочем колесе и кавитаци-онное ограничение [1]. Выполнение ограничения на напор гарантирует, что полученное в результате оптимизации рабочее колесо будет соответствовать заданному изначально режиму работы. Явление кавитации возникает в случае, если величина статического давления в точке потока становится меньше, чем величина давления насыщенного пара при данной температуре [1]:

р < р„. (4)

Известно, что наиболее подвержены кавитации тыльные поверхности лопастей. В случае, когда область кавитации занимает более 15 % площади всей тыльной поверхности лопасти и располагается сразу за входной кромкой, наблюдается резкое падение мощности и КПД турбины. Поэтому ставилось ограничение на размер и расположение области кавитации:

< 0.1,

51

где 5сау = ш£сау; 58ис — площадь тыльной поверхности лопасти; £сау = 1, если выпол-

няется условие возникновения кавитации (4), и £сау = 0, если кавитации нет; ш = ш(/) — весовая функция, вид которой представлен на рис. 7. Параметр I — это нормированное расстояние от входной кромки, измеренное вдоль поверхности лопасти в направлении потока.

5

3.2. Целевые функционалы

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

Рис. 8. Предельные линии тока на рабочей стороне лопасти.

1. Я — кинетическая энергия в выходном сечении рабочего колеса. Кинетическая энергия потока в выходном сечении 50 расчетной области вычисляется по формуле

Я = 1 У V2 (V ■ ¿3). (5)

йо

Здесь V — вектор абсолютной скорости; Ж — вектор нормали к элементарной площадке выходного сечения, равный по модулю площади этого сечения и направленный вовне расчетной области.

2. Я2 — отклонение линий тока от "осесимметричных". Интуиция опытных проектировщиков подсказывает, что более качественной является та лопасть, на поверхности которой предельные линии тока ближе к линиям тока "осесимметричного" потока. "Осесиммет-ричный" поток получается из рассчитанного пространственного потока в межлопастном канале его осреднением в окружном направлении. Соответствующий целевой функционал можно сформулировать следующим образом:

Я = 1^(1 - а(в)ес8в) (6)

где 5 — площадь поверхности лопасти; в — угол между предельной линией тока и линией тока осесимметричного потока (рис. 8), а весовая функция а имеет вид

(К / 1, в < п/2, а(в) = \ ао, в > п/2.

За счет выбора подходящего весового коэффициента а0 минимизация этого целевого функционала позволяет также выполнить одно из основных требований к лопастной системе — отсутствие угла атаки при обтекании лопасти (совпадение линии растекания жидкости с линией входной кромки лопасти).

3. — размер и расположение области кавитации. В минимизируемом функционале учитывается как размер зоны кавитации, так и ее расположение на тыльной стороне лопасти:

Я = . (7)

3.3. Результаты оптимизационных расчетов 3.3.1. Двухцелевая оптимизация

Проведена минимизация кинетической энергии на выходе из рабочего колеса F\ (5) и отклонения линий тока от "осесимметричных" F2 (6) при ограничениях на напор в рабочем колесе и на размер и расположение области кавитации. При поиске фронта Парето использовались следующие значения параметров оптимизационного алгоритма: ashare = 0.01, Tr = 0.3, d = 0.7, у = 0.1, p0 = 100, p = 70. Рассчитано 80 поколений, однако в данном случае сходимость к фронту Парето достигнута за 45 поколений. На рис. 9 в плоскости минимизируемых функционалов изображены лучшие индивидуумы 0-, 3-, 5-, 10-, 20- и 80-го поколений. Решение, составляющее лучший компромисс между функционалами, обозначен как лопасть Opt. На рис. 10 приведены меридиональные распределения скоростей на

Рис. 9. Лучшие индивидуумы в различных поколениях.

Рис. 10. Распределение скорости на выходе из рабочего колеса: слева — начальная, справа — лопасть Opt.

Рис. 11. Линии растекания и линии тока на рабочей стороне лопасти: слева — начальная, справа — лопасть Opt.

Рис. 12. Распределение давления и зона кавитации на тыльной стороне лопасти: слева — начальная (0%), справа — лопасть Opt (5%).

выходе из рабочего колеса, а на рис. 11 — линии растекания и линии тока на рабочей стороне лопасти для начальной лопасти и лопасти Opt. Видно, что для лопасти Opt поле меридиональной компоненты скорости более равномерное, а линия растекания полностью совпала с входной кромкой лопасти. Однако на оптимальной лопасти появилась зона кавитации, составляющая 5% от площади тыльной поверхности лопасти (рис. 12).

3.3.2. Трехцелевая оптимизация

Поскольку при оптимизации исходной лопасти относительно критериев F1 и F2 были ухудшены ее кавитационные качества, проведен оптимизационный расчет, направленный на минимизацию трех функционалов (Fi, F2 и F3) при ограничении на напор. При решении данной задачи число индивидуумов в популяции p было увеличено до 100. На рис. 13 изображена проекция 3Б-фронта Парето, полученного после расчета 150 поколений, на

Рис. 13. Проекция фронта Парето на плоскость функционалов F1 и F3.

Рис. 14. Лопасть 0р"Ъ1. Распределение скорости на выходе из рабочего колеса. Линия растекания и линии тока на рабочей стороне лопасти.

Рис. 15. Отклонение срединной поверхности лопасти от начальной: слева — лопасть Opt, справа — лопасть Opt1. Светлая — начальная срединная поверхность, темная — оптимизированная.

плоскость функционалов Fi и F3. Кружками обозначены решения задачи — фронт Паре-то, линиями — линии уровня функционала F2. Бескавитационная лопасть, которая имеет линию растекания, совпадающую с входной кромкой, и обеспечивает наиболее ровное поле скорости на выходе из рабочего колеса, обозначена Optl. На рис. 14 показаны соответствующие этой лопасти распределение скорости на выходе из РК и линии тока. Как видно из рис. 15, в процессе оптимизации лопасти Opt и Optl существенно отклонились от начальной.

Заключение

Статья посвящена применению генетического алгоритма многоцелевой оптимизации для проектирования лопасти рабочего колеса гидротурбины. На задачах с тестовыми целевыми функционалами показана эффективность предложенного алгоритма нахождения фронта Парето. При двух- и трехцелевой оптимизации лопасти с 16 варьируемыми параметрами геометрии требуется рассчитать несколько большее число поколений (80-150 поколений), чем для одноцелевой задачи (около 30 поколений, см. [1]). С помощью алгоритма удалось найти такую форму лопасти, которая обеспечивает достаточно равномерное поле скорости за рабочим колесом (плавное безударное обтекание рабочей и тыльной сторон) и не приводит к образованию кавитации в потоке. Планируется построить алгоритм многорежимной оптимизации, позволяющий удовлетворять критериям качества лопастной системы на нескольких режимах работы гидротурбины. В настоящее время для ускорения процесса оптимизации разрабатывается параллельная версия алгоритма.

Список литературы

[1] Лобарева И.Ф., Скороспелов В.А., Турук П.А. и др. Об одном подходе к оптимизации формы лопасти гидротурбины // Вычисл. технологии. 2005. Т. 10, № 6. C. 52-74.

[2] Грязин Ю.А., Черный С.Г., Шаров С.В., Шашкин П.А. Об одном методе численного решения трехмерных задач динамики несжимаемой жидкости // Докл. РАН. 1997. Т. 353, № 4. С. 478-483.

[3] Cherny S.G., Sharoy S.V., Skorospeloy V.A., Turuk P.A. Methods for three-dimensional flows computation in hydraulic turbines // Russ. J. Numer. Anal. Math. Modelling. 2003. Vol. 18, N 18. P. 87-104.

[4] Fonseca C.M., Fleming P.J. Genetic algorithms for multiobjective optimization: formulation, discussion and generalization // Proc. of the 5th Intern. Conf. on Genetic Algirithms. 1993. P. 416-423.

[5] Horn J., Nafpliotis N. Multiobjective Optimization Using The Niched Pareto Genetic Algorithm. IlliGAL Report 93005. Urbana: Univ. of Illinois, 1993.

[6] Zitzler E., Deb K., Thiele L. Comparison of multiobjective evolutionary algorithms: empirical results // Evolutionary Computation. 2000. Vol. 8, N 2. P. 173-195.

[7] Deb K., Thiele L., Laumanns M., Zitzler E. Evolutionary Computation Based Multi-Criteria Optimization: Theoretical Advances and Applications. Berlin: Springer-Verlag, 2005.

[8] Tanaka M., Watanabe H., Furukawa Y., Tanino T. GA-based decision support system for multi-criteria optimization // Proc. of the Intern. Conf. on Systems, Man and Cybernetics-2. 1995. P. 1556-1561.

[9] Osyczka A., Kundu S. A new method to solve generalized multicriteria optimization problems using the simple genetic algorithm // Structural Optimization. 1995. Vol. 10. P. 94-99.

[10] Lipej A., Poloni C. Design of Kaplan runner using multi-objective genetic algorithm optimization // J. of Hydraulic Research. 2000. Vol. 38. P. 73-77.

[11] Enomoto Y., Kurosawa S., Suzuki T. Design optimization of Francis turbine runner using multi-objective genetic algorithm // Proc. of 22nd IAHR Symp. on Hydraulic Machinery and Systems. 2004.

[12] Mazzouji F., Francois M., Tomas L. et al. Refinements in Francis turbine design // Hydropower & Dams. 2004. Issue one. P. 53-58.

[13] Fayaretto C., Funazaki K., Tanuma T. The development of a genetic algorithm code for secondary flow injection optimization in axial turbines // Proc. of the Intern. Gas Turbine Congress. Tokyo, Nov. 2003.

[14] Goldberg D.E., Richardson J.J. Genetic algorithms with sharing for multimodal function optimization // Proc. of the 2nd ICGA. NJ. 1987. P. 41-49.

Поступила в редакцию 11 мая 2006 г.