CC BY
45
НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. И. Э. БАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 - 48211. Государственная регистрация №0421200025. ISSN 1994-0408
электронный научно-технический журнал
Применение гибридных алгоритмов глобальной оптимизации
к экстремальным задачам для гидромеханических систем
# 11, ноябрь 2013
Б01: 10.7463/1113.0604082
Сулимов В. Д., Шкапов П. М.
УДК 519.6
Россия, МГТУ им. Н.Э. Баумана spm@bmstu.ru fn3@bmstu.ru
Введение
В состав многих современных изделий и сооружений области высоких технологий, таких как летательные аппараты, реакторные установки АЭС и др., входят гидромеханические системы различного назначения. Разработка и последующая эксплуатация таких систем связаны с поиском решения экстремальных задач двух типов - оптимизации и диагностики. Задачи первого типа возникают, например, при выборе оптимальных параметров систем или реализации оптимального управления. Обеспечение безопасной и эффективной эксплуатации требует решения задач второго типа: коррекции математических моделей и диагностирования систем по результатам косвенных измерений. Входными данными для диагностирования являются результаты экспериментального определения следственных характеристик системы или процесса, включая регистрируемые параметры колебательных и ударных процессов [1]. Искомыми являются причинные характеристики, к которым относятся коэффициенты уравнений расчетной динамической модели, граничные условия, геометрические и ряд других характеристик. В задачах этого типа следует учитывать недифференцируемость и многоэкстремальность критериальных функций ввиду наличия кратных частот, а также неполноты и зашумленности информации, полученной при измерениях. Значительная трудоемкость численного решения обратных спектральных задач
обусловлена их некорректностью, которая чаще всего проявляется в неустойчивости решения относительно погрешностей входных данных [2].
Оптимизационное исследование сложных объектов основано на разработке и последующем уточнении их математических моделей [3-6]. Усложнение модели объекта, в свою очередь, вызывает необходимость создания новых, более эффективных методов оптимизации. Следует отметить, что в аспекте динамики спектры колебаний содержат существенную информацию об исследуемом объекте. Возникают задачи определения оптимальных собственных характеристик системы или процесса, а также использования собственных характеристик для коррекции моделей и диагностирования систем. Одним из активно развивающихся направлений, связанных с проблемой безопасности реакторных установок АЭС, является исследование двухфазных газожидкостных потоков. Некоторые результаты, как численного моделирования, так и экспериментальных исследований двухфазных потоков, представлены в работах [7-10]. Актуальной является задача идентификации аномалий фазового состава теплоносителя в контуре реакторной установки. Используемый далее общий подход основан на разработке и применении математических моделей систем, математических методов расчета основных динамических характеристик систем, методов теории обратных задач, методов глобальной оптимизации.
В первом разделе формулируется задача глобальной оптимизации. Краткое описание гибридных алгоритмов оптимизации дано во втором разделе; там же приведены сведения о реализации алгоритмов в виде прикладных программ. Численный пример применения гибридных алгоритмов рассматривается в третьем разделе, где представлены результаты решения модельной задачи диагностирования фазового состава теплоносителя в циркуляционном контуре реакторной установки.
1. Постановка задач
Некоторые стандартные постановки экстремальных задач для гидромеханических систем приведены в работе [11]. Так, задача коррекции расчетной динамической модели и диагностирования системы как обратная спектральная задача связана с поиском вектора переменных управления, при котором первые N собственных частот модели совпадают с составляющими некоторого заданного ограниченного спектра или близки к ним. Для оценки уровня рассогласования сравниваемых характеристик объекта используется векторный способ описания. Так
как информация о формах колебаний объекта зачастую отсутствует или является существенно неполной, ниже рассматривается только рассогласование между частотными составляющими нормального (при отсутствии аномалий) и заданного спектров. Возможные подходы основаны на минимизации квадратичной функции рассогласования или минимизации максимальной из функций рассогласования спектральных составляющих. Так, для попарно сравниваемых спектральных составляющих может быть построено следующее конечное множество критериев рассогласования
f (x) = Z(x)-Z*(x) x e X с Rn, i e /0 = {1, ... ,N}, где Zi (x), Z* (x) — собственные значения, относящиеся к исходному (текущему) и
заданному спектрам; через Rn обозначено n -мерное вещественное линейное пространство. Необходимо определить такой вектор переменных управления, который приводит к наименьшим отличиям между сравниваемыми спектрами, т.е. следует произвести настройку модели объекта на заданный спектр. Это эквивалентно одновременной минимизации всех N критериев рассогласования: требуется найти
min F (x).
xeX cR n
Здесь векторная целевая функция записывается в виде
F(x) = (fi(x),..., fN(x))T .
Задача коррекции расчетной динамической модели системы формулируется в следующем виде: определить вектор переменных управления x e X, который минимизирует максимальное значение критерия рассогласования, т.е. требуется найти
min max {f (x)}. (Pi)
xeX cR n ie/0
Решением сформулированной дискретной минимаксной задачи (Pi) является
* I * * у
такой вектор x = Ixi,...,xni , принадлежащий множеству допустимых значений, при котором скалярная критериальная функция f (x) = max {fi (x),..., fN (x)} принимает
минимальное значение. В случае, когда f (x ) = 0, спектр частот настраиваемой
модели полностью совпадает с заданным спектром по N низшим частотам. Последнее условие вследствие неполноты экспериментальных данных и погрешностей, полученных при измерениях, не выполняется. Далее рассматривается регуляризованная
задача (P1)S с многоэкстремальной не всюду дифференцируемой критериальной функцией f (х), где S> 0 — параметр регуляризации [11, 12]. Используется
предположение о том, что отношения приращений функций к приращениям аргументов не превышают некоторого порога. Последний определяется ограниченной энергией изменений в системе и может быть описан с помощью константы Липшица.
В обобщение постановок экстремальных задач выбора оптимальных параметров, коррекции расчетной динамической модели и диагностирования системы, принимается следующая формулировка задачи глобальной оптимизации: требуется найти
f (х* )= min f (х), (1)
хеХ с R п
где
X = {х е D : gi(х)< 0, i е I }, , (2)
D = {х е Rп : üj < х] < bj, j е J}. (3)
Использованы обозначения: f (х) — целевая функция; х — вектор переменных
управления; gi (х)— функции ограничений задачи, i е I; I = {1, ... ,m}— конечное
*
множество индексов; X — допустимая область; D — область поиска; х — глобальное решение; п — размерность задачи; J = {1, ... ,п}. Функции f (х), gi (х), i е I, задачи (1) - (3) предполагаются непрерывными липшицевыми. Предполагается также, что критериальная функция f : Rп ^ R является действительной, многоэкстремальной,
не всюду дифференцируемой и для нее задана вычислительная процедура, позволяющая определять значения функции в точках допустимой области. Необходимо также учесть возможную высокую трудоемкость вычисления критериальной функции, что может потребовать значительных вычислительных ресурсов.
2. Гибридные алгоритмы глобальной оптимизации
2.1 Алгоритм М-РСЛЬМБ
Детерминированные методы решения задач глобальной оптимизации многоэкстремальных функций к настоящему времени достаточно хорошо разработаны и находят широкое применение. Однако их эффективность существенно ограничена размерностью задачи. В случае большого числа переменных применяют алгоритмы
стохастической глобальной оптимизации. К ним относятся алгоритмы моделируемого отжига, генетические, управляемого случайного поиска и др. [13-15]. Вместе с тем, чувствительность к выбору параметров стохастических алгоритмов, устанавливаемых пользователем или определяемых содержанием задачи, сильно влияет на скорость сходимости итерационного процесса. В настоящее время значительное внимание уделяется проблеме гибридизации алгоритмов глобальной оптимизации [16, 17]. Каждый из рассматриваемых ниже гибридных алгоритмов объединяет какой-либо стохастический алгоритм сканирования пространства переменных, используемый при определении области, перспективной на глобальный экстремум, и детерминированный метод локального поиска. При этом вычислительная эффективность результирующего алгоритма по сравнению с исходными алгоритмами может значительно возрасти [18, 19].
Одним из наиболее мощных современных стохастических алгоритмов глобальной оптимизации является алгоритм PCA [20]. Существенным шагом алгоритма является сравнительная оценка качества решения, определяемого текущей и предшествующей конфигурациями системы. Пробное приближение принимается с определенной вероятностью, что исключает сходимость к локальному минимуму при поиске глобального решения. Позднее в работе [21] был предложен более мощный алгоритм M-PCA. Работа алгоритмов PCA и M-PCA основана на использовании аналогии с физическими процессами абсорбции и рассеяния частиц при ядерных реакциях. В алгоритме PCA для исследования области поиска используется одна частица. На начальном шаге выбирается пробное решение (Old_Config), которое затем модифицируется посредством стохастического возмущения (Perturbation()), что позволяет найти новое решение (New_Config). С помощью функции Fitness() определяются текущие значения минимизируемой функции и дается сравнительная оценка нового и предыдущего решений, на основании которой новое решение может быть принято или отвергнуто. Если новое решение отвергнуто, то происходит переход к функции Scattering(), реализующей схему Метрополиса. Для сканирования области, перспективной на минимум, применяются функции Perturbation() и Small_Perturbation(). Новое решение принимается, если оно лучше предыдущего (абсорбция); если найденное решение хуже предыдущего, то происходит переход в отдаленную область пространства поиска (рассеяние), что позволяет преодолевать локальные минимумы. Эффективность описанного поиска глобального решения алгоритмом PCA может быть значительно повышена за счет одновременного использования большого числа частиц.
Такой подход реализует алгоритм М-РСА, который непосредственно ориентирован на применение в среде параллельных вычислений. Наилучшее решение определяется с учетом данных о всех частицах, участвующих в процессе. Единственным задаваемым параметром для алгоритма М-РСА является число итераций.
Локальный поиск в гибридном алгоритме должен быть реализован с учетом предположения о недифференцируемости критериальной функции. Поэтому при выборе детерминированного метода локального поиска необходимо учитывать его способность работать с такими функциями. К числу современных методов, специально разработанных для решения задач недифференцируемой оптимизации, относятся, например, ЬипШе-метод с ограниченной памятью [22], субградиентный метод [23] и др. Альтернативный подход основан на построении сглаживающих аппроксимаций критериальных функций и последующем применении мощных классических алгоритмов оптимизации гладких функций [24]. В работе [25] представлен двухпараметрический метод построения сглаживающих аппроксимаций не всюду дифференцируемых функций и предложен вариант метода линеаризации ЬМБ со сглаживанием. Авторами разработан гибридный алгоритм, объединяющий стохастический алгоритм М-РСА сканирования пространства переменных и детерминированный метод ЬМБ локального поиска. Результирующий алгоритм М-РСАЬМБ реализован в виде прикладного программного обеспечения [26].
2.2 Алгоритм М-РСАМЫМ
Решение подзадачи локального поиска для не всюду дифференцируемой критериальной функции можно также получить с помощью методов, не использующих производные. В работе [27] представлен вариант симплекс-метода Нелдера-Мида, сходимость которого доказана теоретически. Кроме того, алгоритм, реализующий модифицированный метод Нелдера-Мида, является робастным для задач с разрывными или зашумленными критериальными функциями. Отмечена более высокая вычислительная эффективность модифицированной версии метода по сравнению с классической. Таким образом, может быть предложен гибридный алгоритм глобальной недифференцируемой оптимизации М-РСАМЫМ, объединяющий стохастический алгоритм М-РСА (общий поиск в пространстве переменных) и модифицированный симплекс-метод Нелдера-Мида (локальный поиск). Ниже представлен фрагмент псевдокода нового гибридного алгоритма М-РСАЫМЫ.
0 Generate an initial solution OldConfig BestFitness = Fitness (OldConfig) Update Blackboard For n = 0 to # of particles
For n = 0 to # of iterations Update Blackboard Perturbation ( )
If Fitness (New Config)) > Fitness (OldConfig) If Fitness (NewConfig) > Best Fitness
Best Fitness := Fitness (New Config)
End If
Old_Config := New_Config Exploration ( )
Else
Scattering ( )
End If End For End For
2. Exploration ( )
For n = 0 to # of iterations
Small_Perturbation ( ) Local search
using Modified Nelder-Mead Simplex Method Check stopping criterion: Find global solution Best Fitness Else continue
If Fitness (New Config) > Best Fitness
Best_Fitness := Fitness (New_Config)
End If
Old_Config := New_Config End For
Return
3. Scattering ( )
pscatt = 1 — ( Fitness (New_Config)) / (Best Fitness)
If psca,, > rand0m(0, 1)
Old Config := random solution
Else
Exploration ( )
End If
Return
В состав алгоритма M-PCAMNM входят также стандартные процедуры Perturbation( ) и Small_Perturbation( ) [21]. Другой гибридный алгоритм M-PCASFC, объединяющий стохастический алгоритм M-PCA и детерминированный метод кривой, заполняющей пространство, используемый при локальном поиске, представлен в работе [28].
3. Численные примеры
Рассматривается модельная задача диагностирования фазового состава теплоносителя в главном циркуляционном контуре серийного блока ВВЭР-1000 [8]. Переменными управления являются относительные значения скорости звука xi в теплоносителе на участках, соответствующих: зоне нагрева теплоносителя в напорном баке системы компенсации объема СКО (x1); выходному объему реактора (x2);
активной зоне реактора (x3); проточной части главного циркуляционного насоса циркуляционной петли с СКО (x4). При отсутствии в теплоносителе второй фазы нормальный спектр ю j определяется максимальными значениями скорости звука на выделенных участках.
*
В первом примере (таблица 1) аномальный спектр юj получен при наличии двухфазной смеси как в выходном объеме, так и в активной зоне реактора, при этом:
x* = 77,5%; x* = 88,0%; x* = 82,5%; x* = 100,0%.
Таблица 1 - Нормальный и аномальный спектры частот колебаний теплоносителя
j 1 2 3 4 5 6 7 8 9 10
ю j , Гц 0,89 6,77 9,82 15,44 15,96 18,94 24,57 26,69 27,07 30,52
* ю*, Гц 0,86 6,77 9,36 15,33 15,96 18,86 21,22 26,67 26,93 29,41
После определения области, содержащей глобальный минимум, завершающие итерации алгоритма M-PCALMS проводятся с использованием градиентной информации для сглаживающих аппроксимаций критериальной функции [25]. Сходимость решения иллюстрируют рисунки 1, 2 (Niter - число итераций).
1.00E+02
9.50E+01 J
9.00E+01
£
8.50E+01
8.00E+01
7.50E+01
— х(1) х(2) х(3) -х-х(4)
0 2 4 6 8 10 12 14
Niter
Рис. 1. Изменение значений переменных управления xt на завершающих итерациях алгоритма
1.00E+03
1.00E+02
S1
z
1.00E+01
1.00E+00
-f(x) ■Nr(p)
10
15
Niter
Рис. 2. Уточнение значений критериальной функции f ^)
и нормы вектора направления поиска Nr (p) на завершающих итерациях алгоритма
* *
Получено приближенное решение задачи: x1 «75,5%; X2 ~ 87,6%;
* *
Xз ~ 84,1%; X4 = 94,6%. Сравнение найденных значений переменных управления с соответствующими данными, определенными постановкой задачи (при наличии аномалии фазового состава теплоносителя), позволяет установить, что относительная погрешность вычисления значений переменных управления не превышает 2,0%, а точность настройки спектра частот имеет порядок 0,01 Гц.
*
Во втором примере (таблица 2) аномальный спектр Юj получен при наличии
двухфазной смеси как в выходном объеме, так и в активной зоне реактора, при этом:
^ ^
x1 = 76,0%; x2 = 87,5%; x3 = 84,0%; x4 = 93,5%. После определения области, содержащей глобальный минимум, завершающие итерации алгоритма М-РСАМЫМ проводятся с использованием модифицированного симплекс-метода Нелдера-Мида.
Таблица 2 - Нормальный и аномальный спектры частот колебаний теплоносителя
0
5
j 1 2 3 4 5 6 7 8 9 10
Ю j, Гц 0,89 6,77 9,82 15,44 15,96 18,94 24,57 26,69 27,07 30,52
* Ю*, Гц 0,86 6,77 9,36 15,33 15,96 18,86 21,22 26,67 26,93 29,41
Сходимость решения иллюстрируют рисунки 3, 4 (Niter - число итераций).
Получено приближенное решение задачи: ~ 75,5%; Х2 ~ 87,6%; Х3 ~ 84,1%;
*
Х4 = 94,6% . При этом относительная погрешность определения значений переменных управления не превышает 2,8%, а точность настройки спектра частот имеет порядок 0,01 Гц.
1.00E+02
9.50E+01
9.00E+01
Д 8.50E+01 и
8.00E+01
7.50E+01
7.00E+01
-о— х1 -»— х2 -А—х3 -X—х4
10
15 20 Niter
25
30 35
Рис. 3. Изменение значений переменных управления Х1 на завершающих итерациях алгоритма
0
5
Рис. 4. Уточнение значений критериальной функции f (х) на завершающих итерациях алгоритма
Выводы
Представлены два гибридных алгоритма глобальной оптимизации, объединяющие современный стохастический алгоритм M-PCA, используемый при сканировании пространства переменных, и методы локального поиска. В первом алгоритме M-PCALMS локальный поиск проводится методом линеаризации со сглаживающимися аппроксимациями критериальной функции. Во втором алгоритме M-PCAMNM при локальном поиске используется модифицированный симплекс-метод Нелдера-Мида. Разработано программное обеспечение, реализующее гибридные алгоритмы. Предложен подход к решению задачи вычислительной диагностики фазового состава теплоносителя в контуре реакторной установки с использованием гибридных алгоритмов. Модельные расчеты показали возможность идентификации аномалий фазового состава теплоносителя на выделенных участках контура с достаточной для приложений точностью.
Работа выполнена при финансовой поддержке Министерства образования и науки РФ (грант Президента РФ по поддержке научных исследований ведущих научных школ РФ, код НШ-4748.2012.8).
Список литературы
1. Ватульян А.О. Обратные задачи в механике деформируемого твердого тела. -М.: ФИЗМАТЛИТ, 2007. - 224 с.
2. Horvath M. Inverse spectral problems and closed exponential systems // Annals of Mathematics. - 2005. - V. 162, № 2. - P. 885-918.
3. Pulecchi T., Casella F., Lovera M. Object-oriented modelling for spacecraft dynamics: Tools and applications // Simulation Modelling and Theory. - 2010. - V. 18, № 1. - P. 63-86.
4. Lopez I., Sarigul-Klijn N. A review of uncertainty in flight vehicle structural damage monitoring, diagnosis and control: Challenges and opportunities // Progress in Aerospace Sciences. - 2010. - V. 46, № 7. - P. 247-273.
5. Bottcher M., KrnPmann R. Primary loop study of a VVER-1000 reactor with special focus on coolant mixing // Nuclear Engineering and Design. - 2010. - V. 240, № 9. -P. 2244-2253.
6. De Oliveira M.V., de Almeida J.C.S. Applications of artificial intelligence techniques in modeling and control of a nuclear power plant pressurizer system // Progress in Nuclear Energy. - 2013. - V. 63. - P. 71-85.
7. Poullikkas A. Effects of two-phase liquid-gas flow on the performance of nuclear reactor cooling pumps // Progress in Nuclear Energy. - 2003. - V. 42, № 1. - P. 3-10.
8. Kinelev V.G., Shkapov P.M., Sulimov V.D. Application of global optimization to VVER-1000 reactor diagnostics // Progress in Nuclear Energy. - 2003. - V. 43, № 1-4. -P. 51-56.
9. Pang S., Chen L., Zhang M., Yin Y., Chen T., Zhou J., Liao D.. Numerical simulation two phase flows of casting filling process using SOLA particle level set method // Applied Mathematical Modelling. - 2010. - V. 34, № 12. - P. 4106-4122.
10. Yang X., Schlegel J.P., Liu Y., Paranjape S., Hibiki T., Ishii M. Experimental study of interfacial area transport in air-water two-phase flow in a scaled 8x8 BWR rod bundle // International Journal of Multiphase Flow. - 2013. - V. 50. - P. 16-32.
11. Сулимов В.Д., Шкапов П.М. Методология решения экстремальных задач для механических и гидромеханических систем // Вестник МГТУ им. Н.Э. Баумана. Сер. «Естественные науки». Специальный выпуск. - 2012, № S8. - С. 17-34.
12. Li X.Y., Law S.S. Adaptive Tikhonov regularization for damage detection based on nonlinear model updating // Mechanical Systems and Signal Processing. - 2010. - V. 24, № 2. - P. 1646-1664.
13. Egea J. A., Vazquez E., Banga J.R., Marti R. Improoved scatter search for the global optimization of computationally expensive dynamic models // Journal of Global Optimization. - 2009. - V. 43, № 2. - P. 175-190.
14. Bertsimas D., Nohadami O. Robust optimization with simulated annealing // Journal of Global Optimization. - 2010. - V. 48, № 3. - P. 323-334.
15. Lu T.-C., Juang J.-C. Quantum-inspired space search algorithm (QSSA) for global numerical optimization // Applied Mathematics and Computation. - 2011. - V. 218, № 6. -P. 2516-2532.
16. La Cruz D., Noguera G. Hybrid spectral gradient method for the unconstrained minimization problem // Journal of Global Optimization. - 2009. - V. 44, № 2. - P. 193-212.
17. Thangaraj R., Pant M., Abraham A., Bouvry P. Particle swarm optimization: hybridization perspectives and experimental illustrations // Applied Mathematics and Computation. - 2011. - Vol. 217, № 12. - P. 5208-5226.
18. Gaviano M., Lera D., Steri A.M. A local search method for continuous global optimization // Journal of Global Optimization. - 2010. - V. 48, № 1. - P. 73-85.
19. Voglis C., Parsopoulos K.E., Papageorgiou D.G., Lagaris I.E., Vrahatis M.N. MEMSODE: A global optimization software based on hybridization of population-based algorithms and local searches // Computer Physics Communications. - 2012. - V. 183, № 2. -P. 1139-1154.
20. Sacco W.F., de Oliveira C.R.E. A new stochastic optimization algorithm based on particle collisions // Proceedings of the 2005 ANS Annual Meeting. Transactions of the American Nuclear Society, 2005. - V. 92. - P. 657-659.
21. Luz E.F.P., Becceneri J.C., de Campos Velho H.F.. A new multi-particle collision algorithm for optimization in a high performance environment // Journal of Computational Interdisciplinary Sciences, 2008. - V. 1. - P. 3-10.
22. Karmitsa N., Makela M.M. Limited memory bundle method for large bound constrained optimization: convergence analysis // Optimization Methods & Software. - 2010.
- V. 25, № 6. - P. 895-916.
23. Bagirov A.M., Jin L., Karmitsa N., Al Nuaimat A., Sultanova N. Subgradient method for nonconvex nonsmooth optimization // Journal of Optimization Theory and Applications. - 2013. - V. 157, № 2. - P. 416-435.
24. Beck A., Teboulle M. Smoothing and first order methods: a unified framework // SIAM Journal of Optimization. - 2012. - V. 22, № 2. - P. 557-580.
25. Сулимов В.Д. Локальная сглаживающая аппроксимация в гибридном алгоритме оптимизации гидромеханических систем // Вестник МГТУ им. Н.Э. Баумана. Сер. «Естественные науки». - 2010, № 3. - С. 3-14.
26. Сулимов В.Д., Шкапов П.М. Глобальная минимизация многомерной многоэкстремальной липшицевой целевой функции с использованием гибридного алгоритма V-PCALMS. Свидетельство о государственной регистрации программы для ЭВМ № 2012617546. Зарегистрировано в Реестре программ для ЭВМ 21 августа 2012 г.
- Федеральная служба по интеллектуальной собственности, патентам и товарным знакам, 2012.
27. Price C.J., Coope I.D., Byatt D. A convergent variant of the Nelder-Mead algorithm // Journal of Optimization Theory and Applications. - 2002. - V. 113, № 1. - P. 519.
28. Sulimov V.D., Shkapov P.M. Application of hybrid algorithms to computational diagnostic problems for hydromechanical systems // Journal of Mechanics Engineering and Automation. - 2012. - V. 2, № 12. - P. 734-741.
SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU
SCIENCE and EDUCATION
EL №FS77 - 48211. №0421200025. ISSN 1994-0408
electronic scientific and technical journal
Application of hybrid global optimization algorithms to extremum problems in hydro-mechanical systems # 11, November 2013 DOI: 10.7463/1113.0604082 Sulimov V.D., Shkapov P.M.
Bauman Moscow State Technical University, 105005, Moscow, Russian Federation
spm@bmstu.ru fn3@bmstu.ru
This article deals with problems of global optimization, model updating and diagnostics of hydro-mechanical systems. It was assumed that criterion functions were continuous, Lipschitzian, multiextremal and not always differentiable. Two novel hybrid algorithms were proposed; these algorithms use the modern stochastic Particle Collision Algorithm based on an analogy with absorption and scattering processes for nuclear particles, for scanning a search space. Local search was implemented with the use of a linearization method with smoothing approximations of criteria for the first algorithm, and a convergent variant of the Nelder-Mead simplex method was used for the second algorithm. Some results on solving model problems on computational diagnostics of the coolant phase constitution in the reactor primary circuit were presented.
Publications with keywords: global optimization, hydromechanical system, extremal problem, Lipschitzian function, criterion function, local search, hybrid algorithm Publications with words: global optimization, hydromechanical system, extremal problem, Lipschitzian function, criterion function, local search, hybrid algorithm
References
1. Vatul'yan A.O. Obratnye zadachi v mekhanike deformiruemogo tverdogo tela [Inverse problems in mechanics of deformable solids]. Moscow, Fizmatlit, 2007. 224 p.
2. Horvath M. Inverse spectral problems and closed exponential systems. Annals of Mathematics, 2005, vol. 162, no. 2, pp. 885-918.
3. Pulecchi T., Casella F., Lovera M. Object-oriented modelling for spacecraft dynamics: Tools and applications. Simulation Modelling and Theory, 2010, vol. 18, no. 1, pp. 63-86.
4. Lopez I., Sarigul-Klijn N. A review of uncertainty in flight vehicle structural damage monitoring, diagnosis and control: Challenges and opportunities. Progress in Aerospace Sciences, 2010, vol. 46, no.7, pp. 247-273.
5. Böttcher M., Krüßmann R. Primary loop study of a VVER-1000 reactor with special focus on coolant mixing. Nuclear Engineering and Design, 2010, vol. 240, no. 9, pp. 2244-2253.
6. De Oliveira M.V., de Almeida J.C.S. Applications of artificial intelligence techniques in modeling and control of a nuclear power plant pressurizer system. Progress in Nuclear Energy, 2013, vol. 63, pp. 71-85.
7. Poullikkas A. Effects of two-phase liquid-gas flow on the performance of nuclear reactor cooling pumps. Progress in Nuclear Energy, 2003, vol. 42, no. 1, pp. 3-10.
8. Kinelev V.G., Shkapov P.M., Sulimov V.D. Application of global optimization to VVER-1000 reactor diagnostics. Progress in Nuclear Energy, 2003, vol. 43, no. 1-4, pp. 51-56.
9. Pang S., Chen L., Zhang M., Yin Y., Chen T., Zhou J., Liao D. Numerical simulation two phase flows of casting filling process using SOLA particle level set method . Applied Mathematical Modelling, 2010, vol. 34, no. 12, pp. 4106-4122.
10. Yang X., Schlegel J.P., Liu Y., Paranjape S., Hibiki T., Ishii M. Experimental study of interfacial area transport in air-water two-phase flow in a scaled 8*8 BWR rod bundle. International Journal of Multiphase Flow, 2013, vol. 50, pp. 16-32.
11. Sulimov V.D., Shkapov P.M. Metodologiya resheniya ekstremal'nykh zadach dlya mekhanicheskikh i gidromekhanicheskikh sistem [Methodology of extremal problems solving for material and hydromechanical systems]. VestnikMGTUim. N.E. Baumana. Ser. Estestvennye nauki [Herald of the Bauman MSTU. Ser. Natural science], 2012, spec. iss. no.8, pp. 17-34.
12. Li X.Y., Law S.S. Adaptive Tikhonov regularization for damage detection based on nonlinear model updating. Mechanical Systems and Signal Processing, 2010, vol. 24, no. 2, pp. 1646-1664.
13. Egea J.A., Vazquez E., Banga J.R., Marti R. Improoved scatter search for the global optimization of computationally expensive dynamic models. Journal of Global Optimization, 2009, vol. 43, no. 2, pp. 175-190.
14. Bertsimas D., Nohadami O. Robust optimization with simulated annealing. Journal of Global Optimization, 2010, vol. 48, no. 3, pp. 323-334.
15. Lu T.-C., Juang J.-C. Quantum-inspired space search algorithm (QSSA) for global numerical optimization. Applied Mathematics and Computation, 2011, vol. 218, no. 6, pp. 2516-2532.
16. La Cruz D., Noguera G. Hybrid spectral gradient method for the unconstrained minimization problem. Journal of Global Optimization, 2009, vol. 44, no. 2, pp. 193-212.
17. Thangaraj R., Pant M., Abraham A., Bouvry P. Particle swarm optimization: hybridization perspectives and experimental illustrations. Applied Mathematics and Computation, 2011, vol. 217, no. 12, pp. 5208-5226.
18. Gaviano M., Lera D., Steri A.M. A local search method for continuous global optimization. Journal of Global Optimization, 2010, vol. 48, no. 1, pp. 73-85.
19. Voglis C., Parsopoulos K.E., Papageorgiou D.G., Lagaris I.E., Vrahatis M.N. MEMSODE: A global optimization software based on hybridization of population-based algorithms and local searches. Computer Physics Communications, 2012, vol. 183, no. 2, pp. 1139-1154.
20. Sacco W.F., de Oliveira C.R.E. A new stochastic optimization algorithm based on particle collisions. Transactions of the American Nuclear Society. Proceedings of the 2005 ANS Annual Meeting, 2005, vol. 92, pp. 657-659.
21. Luz E.F.P., Becceneri J.C., de Campos Velho H.F. A new multi-particle collision algorithm for optimization in a high performance environment. Journal of Computational Interdisciplinary Sciences, 2008, vol. 1, pp. 3-10.
22. Karmitsa N., Makela M.M. Limited memory bundle method for large bound constrained optimization: convergence analysis. Optimization Methods & Software, 2010, vol. 25, no. 6, pp. 895-916.
23. Bagirov A.M., Jin L., Karmitsa N., Al Nuaimat A., Sultanova N. Subgradient method for nonconvex nonsmooth optimization. Journal of Optimization Theory and Applications, 2013, vol. 157, no. 2, pp. 416-435.
24. Beck A., Teboulle M. Smoothing and first order methods: a unified framework. SIAM Journal of Optimization, 2012, vol. 22, no. 2, pp. 557-580.
25. Sulimov V.D. Lokal'naya sglazhivayushchaya approksimatsiya v gibridnom algoritme optimizatsii gidromekhanicheskikh sistem [Local smoothing approximation in hybrid
algorithm of optimization of hydromechanical systems]. VestnikMGTUim. N.E. Baumana. Ser. Estestvennye nauki [Herald of the Bauman MSTU. Ser. Natural science], 2010, no. 3, pp. 3-14.
26. Sulimov V.D., Shkapov P.M. Global'naya minimizatsiya mnogomernoy mnogoekstremal'noy lipshitsevoy tselevoy funktsii s ispol'zovaniem gibridnogo algoritma V-PCALMS [Global minimization of multidimensional multiextremal Lipschitz objective function using hybrid algorithm].The Certificate on state registration of computer program, no. 2012617546, 2012.
27. Price C.J., Coope I.D., Byatt D. A convergent variant of the Nelder-Mead algorithm. Journal of Optimization Theory and Applications, 2002, vol. 113, no. 1, pp. 5-19.
28. Sulimov V.D., Shkapov P.M. Application of hybrid algorithms to computational diagnostic problems for hydromechanical systems. Journal of Mechanics Engineering and Automation, 2012, vol. 2, no. 12, pp. 734-741.
