Оптимизация рабочих колес центробежных компрессоров методами вычислительной газодинамики: выбор параметров расчетной модели
В статье рассматриваются вопросы выбора параметров сеточной модели и расчетной области при решении задач оптимизации рабочих колес центробежного компрессора методами вычислительной газодинамики. При решении таких задач поиск и применение оптимальных параметров сеточной модели, расчетной области и настроек решателя позволяет обеспечить высокую точность моделирования при наиболее эффективном и производительном использовании машинного времени. Исследование проведено в комплексе Numeca Fine/Turbo для моделей турбулентности Spalart-Allmaras и Shear Stress Transport на примере двух рабочих колес: высоконапорное с ѱт =0,71 и Ф=0,064, и низконапорное ѱт =0,43 и Ф=0,06. Установлено, что выбор оптимальных параметров постановки задачи значительно сокращают время получения сошедшегося решения и облегчают дальнейшее решение оптимизационных задач.
рабочее колесо, центробежный компрессор, расчетная область, сеточная модель, оптимизация, Numeca Fine/Turbo
The article deals with the choice of the grid models and the computational domain parameters in the solution of centrifugal compressor’s impeller optimization problems with using computational fluid dynamics. In solving such problems searching and applying the optimal parameters of the grid model, computational domain, solver settings allows to provide high accuracy simulation with the most effective and productive use of the machine time. Research was conducted in the Numeca Fine/Turbo program complex for Spalart-Allmaras and Shear Stress Transport turbulence models on the example of two impellers: high-pressure with ѱт =0,71 and Ф=0,064, low-pressure with ѱт =0,43 and Ф=0,06. It is established that choose the optimum parameters significantly reduce the time of obtaining converging solutions and facilitate further the solution of optimization problems.
impeller, centrifugal compressor, computational domain, grid model, optimization, Numeca Fine/Turbo
Обозначения:
ѱт – коэффициент теоретического напора
Ф – условный коэффициент расхода
D2 – наружный диаметр рабочего колеса центробежного компрессора
b2 – ширина канала на выходе из рабочего колеса
b3 – ширина канала на входе в диффузор
y+ – сеточный параметр
Сокращения:
РК – рабочее колесо
БЛД – безлопаточный диффузор
КПД – коэффициент полезного действия, эффективность
CFD – Computational Fluid Dynamics (вычислительная гидрогазодинамика)
SA – модель турбулентности Spalart-Allmaras
SST – модель турбулентности Shear Stress Transport
Одним из основополагающих факторов, определяющих эффективность работы компрессорных установок, является качество проектирования проточной части. В современном машиностроении все большее внимание уделяется численной оптимизации проточных частей. Оптимизация по наиболее важным параметрам работы компрессора (коэффициент полезного действия, напор) подразумевает под собой проведение множества вариативных расчетов исследуемой проточной части. В отличие от оптимизационных циклов, построенных на алгоритмах решения обратной задачи, математическое моделирование прямой задачи позволяет оценивать вклад всех геометрических параметров в характеристики проточных частей. Опыт кафедры КВиХТ в оптимизации методами CFD и возможности данной процедуры отражены в работах [1, 2, 3, 4, 5].
Так как во многих программных пакетах численной гидрогазодинамики реализован полный цикл оптимизационных алгоритмов, то отпадает необходимость вручную перестраивать геометрию и сеточные модели. Таким образом, проблема остается только в больших объемах машинного времени, затрачиваемых на расчет. Так как данные расчеты эффективно производятся с помощью мощных вычислительных кластеров [6], то экономия машинного времени может стать значительным фактором. Поэтому увеличение эффективности решения подобных задач состоит в уменьшении времени на оптимизационный шаг. Этого можно добиться уменьшением сеточных моделей, выбором оптимальных расчетных областей и граничных условий. С другой стороны, необходимо выбрать такие параметры сеточной модели, при которых получаемые результаты будут корректно отражать общие результаты и изменения оптимизационных шагов. Таким образом, оптимизация геометрических параметров ступеней во многом зависит от точности применяемых расчетных моделей [7].
Цель работы и объект исследования
Целью работы является определение степени влияния выбора параметров расчетной области и сеточной модели на примере рабочих колес центробежных компрессоров на результаты и скорость получения решения задачи. Определяется наилучшая комбинация рассматриваемых параметров для использования при постановке задачи оптимизации рабочего колеса методами вычислительной газодинамики.
Объектом исследования являются два центробежных радиальных рабочих колеса с цилиндрическими лопатками:
РК1 – рабочее колесо со средней линией лопаток по дуге окружности. Условный расчетный коэффициент расхода Фр = 0,06; расчетный коэффициент теоретического напора ѱтр = 0,43; диаметр рабочего колеса D2 = 0,82м; относительный диаметр втулки = 0,23; угол выхода лопаток βЛ2 = 32°; количество лопаток z=13.
РК2 – рабочее колесо с средней линией лопаток, образованной по двум сопряженным параболам. Условный расчётный коэффициент расхода Фр = 0,064; расчётный коэффициент теоретического напора ѱтр = 0,715; диаметр рабочего колеса D2 = 0,402 м; относительный диаметр втулки = 0,3; угол выхода лопаток βЛ2 = 75,7°; количество лопаток z=24.
Построение сеточных моделей и проведение расчетов выполнены с помощью программного комплекса Numeca Fine/Turbo, преимущества которого описаны в работе [8], а основные аспекты работы приведены в [9]. Для всех расчетов задавались одинаковые граничные условия: стационарная постановка; рабочее тело – идеальный воздух; граничное условие входа – полное давление и температура; граничное условие выхода – массовый расход. Моделирование произведено для двух моделей турбулентности – Spalart-Allmaras и Shear Stress Transport. В расчете не учитывались области уплотнений по покрывающему и основному дискам. Таким образом, приведенные значения политропного КПД – гидравлические. Все расчеты проведены для 1 межлопаточного канала. Такой подход возможен, так как в других каналах течение периодически повторяется.
Результаты расчета отражены в значениях эффективности, которые рассчитаны по сечениям 0-0 (вход в РК) и 2’-2’ (выход из РК) согласно рисунку 1; в количестве итераций, необходимых для получения условно сошедшегося решения (критерии – падение среднеквадратичных невязок до 10-4, выход на неизменное значение параметров эффективности и напора в процессе решения); в обезразмеренном времени, приходящимся на 1 итерацию (значения отнесены к минимальному времени на 1 итерацию в каждой серии расчетов); в обезразмеренном машинном времени, необходимом для получения условно сошедшегося решения (значения отнесены к минимальному машинному времени в каждой серии расчетов). Время, приходящееся на 1 итерацию в общем случае зависит от размерности сетки и, в случае применения, от разбития решения на процессорные ядра. Машинное время в совокупности зависит от количества итераций для получения решения и от времени на 1 итерацию.
Исследование влияния положения выходной границы рабочего колеса
Отрывные течения и, так называемый, след после рабочего колеса может оказать существенное влияние на результаты расчета, если выходная область определена близко к выходу из рабочего колеса. Исследование проводилось при фиксированной входной границе путем удлинения участка безлопаточного диффузора с b3=b2 после рабочего колеса. Размер сеточных моделей увеличивался соответственно увеличению элементов на участке БЛД при его удлинении. Эскиз вариаций выходных границ расчетной области приведен на рисунке 1. Съем параметров для определения эффективности произведен в сечениях 0-0 и 2’-2’ ( =1,05D2). Результаты расчета для обоих рабочих колес представлены на рисунке 2 и 3.
Рисунок 1 – Схема контрольных сечений и вариации расчетной области
Рисунок 2 – Результаты оценки влияния выходной границы для низконапорного колеса РК1 по сечениям 0-0 – 2’-2’
Рисунок 3 – Результаты оценки влияния выходной границы для высоконапорного колеса РК2 по сечениям 0-0 – 2’-2’
Для низконапорного колеса разница полученных значений КПД при отдалении выходной границы расчетной области невелика, и при значении Dвых=1,15D2 уже выходит на постоянную величину. Количество необходимых итераций для сведения решения имеет степенную зависимость и значительно снижается при отдалении выходной границы на 1,3D2 и далее. Однако при отдалении границы возрастает и общее количество элементов, что ведет к увеличению времени, приходящемся на одну итерацию. В совокупности наиболее оптимальным вариантом с точки зрения машинного времени оказался вариант с выходной границей, расположенной на 1,3D2. Дальнейшее увеличение не ведет к значительному уменьшению необходимых для сходимости итераций при увеличении затрачиваемого времени на 1 итерацию.
В случае высоконапорного колеса положение выходной границы сильнее влияет на получаемую в результате решения эффективность рабочего колеса из-за большей неоднородности потока на выходе с лопаток. Поэтому даже по этой причине в высоконапорных рабочих колесах следует отодвигать выходную границу на 1,25D2 и более. С другой стороны, по результатам расчета сходимость решения резко ухудшается после 1,2D2. Количество необходимых итераций возрастает в 2,3 раза на модели турбулентности SA при удалении границы c 1,2D2 до 1,4D2. Это связано с развитыми отрывными течениями в диффузоре на удлиненных участках, и, следственно, с пульсациями и неоднородностями массового расхода на выходной границе, что замедляет схождение решения. На приведенных графиках ширина диффузора равна ширине канала на выходе из рабочего колеса, т.е. b3/b2=1. Большая ширина диффузора способствует возникновению отрывного течения в диффузоре, поэтому кроме исходного варианта был рассмотрен вариант с b3/b2=0,8. Сужение диффузора по результатам расчета не дало эффекта — сходимость после удаления выходной границы от 1,2D2 не только не улучшилась, но и общий уровень времени сходимости увеличился по сравнению с исходной постановкой задачи. В данном случае наиболее приемлемым вариантом видится выбор выходной границы 1,25D2, где в результате решения значение эффективности уже можно считать приемлемым, но и скорость сходимости еще не критично возросла.
Суммируя по влиянию выходной границы – для высоконапорных колес положение выходной границы более значительно влияет на значение КПД численного эксперимента; из-за большей неоднородности потока по сравнению с низконапорными на сведение задачи требуется большее количество итераций, и, следовательно, машинного времени; при этом на участке 1,05-1,2D2 количество необходимых итераций уменьшилось в 1,4 раза против 2,1 у низконапорного; при увеличении Dвых свыше 1,25D2 сходимость задачи резко ухудшается из-за отрывных течений и нестационарных явлений на участке БЛД. Для низконапорных колес влияние положения выходной границы на эффективность сравнительно мало; при удлинении безлопаточного участка сходимость значительно улучшается, что позволило добиться в рассматриваемом случае снижения машинного времени расчета более чем в 1,7 раза. Логически удлинение участка оказывается ограничено оптимальным соотношением необходимых числом итераций для сведения задачи и времени, приходящимся на 1 итерацию.
Для всех серий расчетов следует отметить следующее: модель турбулентности SA, являясь более простой, чем SST, показывает на 10-20% более быструю итерационную сходимость решений и на 5-15% меньшее время, приходящееся на итерацию. А модель SST, в свою очередь, показывает завышение КПД на 0,5-1% в абсолютных значениях по сравнению с моделью SA. Данное завышение было отмечено не только в данной работе, но и в других исследованиях, проводимых на кафедре КВиХТ СПбПУ в пакете Numeca. В целом кривые для обеих моделей турбулентности эквидистантны.
Исследование влияния положения входной границы рабочего колеса
При расчетах рабочих колес положение входной границы часто выбирается случайно или принимается с учетом геометрической формы участка, предшествующему рабочему колесу, но при этом не рассматривается входное устройство. Однако случайный или намеренный выбор длинного входного участка может быть не обоснован. В случае формального учета входной камеры необходимо задавать профили скоростей, давлений на входной границе, что в свою очередь проблемно в реализации и несет свои сложности. Иначе, задавая равномерные поля, не будет отражена реальная картина течения на входном участке, поэтому и рациональность увеличения сеточной модели за счёт удлинения входного участка не очевидна. Поэтому при решении задач, где фигурируют множественные расчеты, в целях экономии времени следует выбирать оптимальное положение входной границы.
В контексте оценки оптимальной длины входного участка перед рабочим колесом рассматривались значения длин 0…0,75D2. Результаты расчета приведены на рисунках 4 и 5.
Рисунок 4 – Результаты оценки влияния входной границы для низконапорного колеса по сечениям 0-0 – 2’-2’
Рисунок 5 – Результаты оценки влияния входной границы для высоконапорного колеса по сечениям 0-0 – 2’-2’
Для низконапорного колеса уменьшение длины входного участка привело к линейному увеличению эффективности на модели турбулентности SA (для рассмотренной задачи 0,3%), и к постоянному значению КПД для модели SST. При нулевой длине произошло резкое падение эффективности на 0,2% для обеих используемых моделей турбулентности. Данное падение вызвано тем, что профиль скорости на данной координате уже перестраивается для входа на лопатки и отличается от задаваемого однородным полем в случае расположения входной границы на этой координате. Это ведет к изменению обтекания лопаток и разнице в результатах. Поэтому, несмотря на продолжающееся ускорение сходимости, не следует выбирать слишком короткий входной участок. В общем, положение координаты, где профиль скорости значительно перестроится, зависит от комплекса факторов, например – ширины колеса, длины участка под уплотнения покрывающего диска, плавности меридионального поворота и разницы проходных площадей в сечениях 0-0 и 1-1. Это демонстрирует расчет высоконапорного колеса, где конструктивно предусмотрена более длинная система уплотнений покрывающего диска. Поэтому при нулевой длине входного участка профиль скорости еще не оказывает воздействия на обтекание лопатки, но отразился на резком ухудшении сходимости. Влияние на сходимость положения входной границы заключается в ускорении сходимости при ее приближении к рабочему колеса. Основной фактор ускорения сходимости — уменьшение общего количества элементов сетки, аналогично описанному далее при исследовании влияния общего количества сеточных элементов. Таким образом, если не моделируется входное устройство, то не следует выбирать длинный входной участок, увеличивая размерность сетки. Для обоих рассмотренных колес оптимальная длина с точки зрения результатов расчетов – 0,125D2.
Исследование влияния выбора сеточных параметров
Вопрос выбора размерности сетки при верификации численных экспериментов решается путем стандартного исследования на сеточную независимость. Такое исследование позволяет определить, начиная с какого количества узлов получаемое решение практически перестает зависеть от густоты сеточной модели. С точки зрения оптимизационных задач необходим более точный выбор сеточных параметров, обеспечивающих приемлемые результаты при минимальном затрачиваемом времени на решение. Сеточный генератор AutoGrid 5 программного комплекса Numeca позволяет автоматически генерировать блочно-структурированные параметризированные сетки. В автоматическом режиме для радиальных колес обычно можно использовать два вида топологии: HOH и H&I. Принципиальная разница между сетками с разными топологиями показана на рисунке 6 и заключается в том, что для сетки типа H&I, в отличие от типа HOH, характерно межсекторное разделение вблизи поверхности лопатки. В этой области разрешается пограничный слой, и течение здесь более сложное, чем в ядре потока. Такое разделение секторов может вести к менее интенсивной сходимости решения из-за осреднений, связанных с налагаемым граничным условием периодичности вблизи поверхности лопатки.
Рисунок 6 – сеточная топология HOH (слева) и H&I (справа)
В общем случае выбор варианта топологии обосновывается качеством получаемой сеточной модели: минимальной скошенностью элементов, максимальным коэффициентом расширения соседних элементов и т.д. Поэтому не для всех конфигураций рабочих колес можно применить оба варианта топологии. В случае, когда для рабочего колеса можно построить равноценные по качеству сетки для обеих топологий, следует производить выбор топологии с точки зрения скорости получаемого решения. Рассматривались сетки, построенные с помощью двух топологий для рабочих колес. Количество сеточных элементов сеток РК1 и РК2 – 1,15 млн., минимальный угол скошенности для РК1 – 45°, для РК2 – 36°. Результаты расчета для режима расчетного расхода представлены в таблице 1.
Таблица 1. Результаты расчета вариантов сеточной топологии при двух моделях турбулентности
РК1 |
РК2 | |||||||
SA |
SST |
SA |
SST |
|||||
HOH | H&I | HOH | H&I | HOH | H&I | HOH |
H&I |
|
Машинное время расчета, сек |
1538 |
1581 |
1630 | 1652 | 1611 | 1633 | 1890 | 1900 |
Итерация достижения сошедшегося решения | 696 | 835 | 830 | 1022 | 1070 | 1305 | 1390 |
1606 |
Таким образом, выбор оптимальной сеточной топологии позволяет по результатам исследования экономить до 20% машинного времени, что является ощутимым при проведении большого количества расчетов.
Стандартная процедура при проведении разнообразных численных экспериментов – исследование на сеточную независимость. Выводы по таким исследованиям приведены во многих работах [10, 11, 12].
С уменьшением размеров элементов сетки происходит увеличение их количества, в связи с чем при расчетах изменение диссипации энергии происходит более плавно, а расчетные потери КПД в элементах проточной части уменьшаются. Для рассматриваемых рабочих колес произведена оценка и выбор оптимального количества элементов для дальнейших оптимизационных расчетов. Результаты для низконапорного РК1 представлены на рисунке 7. Результаты для РК2 качественно повторяют результаты для РК1.
Рисунок 7 – Результаты оценки общего количества элементов сеточной модели для низконапорного колеса по сечениям 0-0 – 2’-2’
С ростом числа элементов сетки КПД низконапорного РК1 практически перестает меняться, начиная с 0,8 млн. элементов. Время, приходящееся на 1 итерацию растет линейно. Количество итераций для сведения решения имеет зависимость близкую к линейной. В итоге машинное время получения сошедшегося решения имеет нелинейный характер и значительно возрастает при размерах сеточных моделей, превышающих 1 млн. элементов. Поэтому для данного рабочего колеса будет рационален выбор сеточной модели размером 0,8-1 млн. элементов для одного лопаточного сектора – это практически не оказывает влияние на получаемые результаты при экономии вычислительных и временных ресурсов.
Сеточный параметр y+ — безразмерное расстояние от первого узла сетки до стенки. Данный критерий важен при выборе модели турбулентности – низко- или высокорейнольдсовых, которые определяют требование к сетке. Например, при использовании стандартных пристенных функций (высокоренольдсовые модели) необходимо избегать значений y+ ме́ньших, чем приблизительно 30, и бо́льших, чем 300 для задач расчёта потока газа. С масштабируемыми пристенными функциями (низкорейнольдосвые модели) значения этого параметра дают информацию о разрешающей способности сетки вблизи стенки. Использование масштабируемой пристенной функции обязывает иметь достаточно тонкий слой элементов, чтобы разложить пограничный слой. Рекомендуемый в литературе диапазон y+ для таких моделей находится в диапазоне от 0,001 до 1. Для определения данного параметра и его дальнейшей корректировки обычно проводится пробный расчет. Определяется значение критерия y+ во всей расчетной области, после чего корректируется глобально или локально размер первой пристенной ячейки для удовлетворения критерия. Влияние параметра рассмотрено в работе [13].
Считается, что низкорейнольдсовые модели точнее в предсказании точек отрывов, но, соответственно, такие модели требовательнее к размерности сетки. Рассматриваемые в работе модели – низкоренойльдсовые. Рассмотрено влияние параметра y+ в рекомендуемом диапазоне от 0,001 до 1 и при выходе за него (от 1 до 48) на решение и сходимость задачи.
На графике 8 приведены результаты для вариантов параметра y+. В таблице 2 для соответствующих вариантов приведено среднее значение y+ в межлопаточном канале рабочего колеса, максимальное значение во всей расчетной области и соответствующий размер первой пристенной ячейки.
Таблица 2. Вариации размера ячеек и значения параметра y+
Вариант | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
РК1 SA&SST | ||||||||
Размер ячейки, мм | 0,001 | 0,0015 | 0,005 | 0,01 | 0,03 | 0,08 | 0,15 | 0,3 |
Средний y+ | 0,17 | 0,25 | 0,82 | 1,75 | 5,5 | 12,6 | 19,7 | 30,4 |
Максимальный y+ | 0,59 | 0,86 | 3,03 | 6,23 | 14,6 | 25,2 | 34,5 | 47,7 |
РК2 SA&SST | ||||||||
Размер ячейки, мм | 0,0005 | 0,0007 | 0,0012 | 0,005 | 0,01 | 0,03 | 0,08 | 0,15 |
Средний y+ | 0,12 | 0,165 | 0,276 | 1,17 | 2,5 | 7,3 | 16,0 | 24,4 |
Максимальный y+ | 0,38 | 0,53 | 0,92 | 4,08 | 7,8 | 16,1 | 28,1 | 39,9 |
Рисунок 8 – Результаты оценки влияния параметра y+ сеточной модели для низконапорного колеса
Как видно из графиков, попадание y+ в рекомендуемый диапазон от 0,001 до 1 оказывает минимальное воздействие на сходимость решения. То есть скорость сходимости в данном диапазоне остается постоянной. Однако получаемые результаты эффективности достаточно сильно различаются – при увеличении среднего значения параметра y+ c 0,17 до 0,82 КПД низконапорного рабочего колеса увеличивается на 0,5%, а затем при y+=1,75 падает на 0,7%. Для высоконапорного колеса результаты качественно аналогичны за исключением некоторого ускорения сходимости для модели SST при стремлении y+ к 1,17.
Увеличение y+ за рекомендуемый интервал ведет к некорректному росту эффективности колеса вплоть до 99% КПД абс. Это связано с значительным уменьшением потерь, в том числе и в результате нивелирования отрывных течений из-за недостаточной размерности сетки для разрешения пограничного слоя. Вместе с некорректными результатами значительно возрастает и количество итераций, необходимых для получения сошедшегося решения. При примерно одинаковом времени, затрачиваемом на 1 итерацию, машинное время получения сведенного решения возрастает только пропорционально необходимому количеству итераций. В случае высокнапорного колеса РК2 результаты и время сходимости ведут себя аналогично низконапорному РК1. Однако решение для модели SST показало более жесткие требования к параметру y+: при y+>2,5 решение прерывается. Таким образом, для модели SST при расчетах высоконапорных рабочих колес прерывание процесса решения может быть связано не только с общим качеством сетки, но и с неприемлемо большим размером первого пристенного элемента. Модель SA в этом случае показала возможность расчета даже при значительном отклонении параметра y+ от рекомендуемого.
Моделируя ситуацию оптимизации исследуемых рабочих колес при выборе случайных и оптимальных параметров, получается следующее:
Случайные параметры: входная граница определена на удалении 0,5D2, выходная граница – 1,1D2; выбрана сеточная топология H&I; высота пристенной ячейки 0,01мм (соответствующее среднее значение y+=1,75); общее количество элементов 1,15млн.
Оптимальные параметры: входная граница определена на удалении 0,125D2, выходная граница – 1,3D2; выбрана сеточная топология HOH; высота пристенной ячейки 0,001мм (соответствующее среднее значение y+=0,17); общее количество элементов 0,92млн.
По результатам расчета для случайно заданных параметров время на 1 итерацию на одном процессорном ядре составило 3 секунды, а количество необходимых итераций для получения сошедшегося решения – 775. Следовательно, общее машинное время расчета для одного оптимизационного приближения составило бы 2322 секунды. Для оптимального набора параметров данные цифры соответственно равны 2,4 секунды на итерацию, 227 итераций для сошедшегося решения и общее время расчета – 535 секунд, что в 4,25 раза быстрее, чем при неоптимальных параметрах расчетной модели. Таким образом, 100 оптимизационных шагов заняли бы 64 часа для случайных и 15 часов для оптимальных параметров. Это существенно, учитывая тот факт, что визуально данные расчетные области и сеточные модели схожи. Также в комплексе Numeca реализованы специальные экспертные настройки CPU Booster, которые позволяют производить моделирование с числом Куранта, равным 1000. Число Куранта (называемое также число Куранта-Фридрихса-Леви) характеризует физический масштаб времени расчёта и определяется по формуле:
\[ CFL=\frac{u\Delta t}{\Delta x} \]
где \( {u} \) – скорость переноса, \( {\Delta t} \) – временной шаг, \( {\Delta x} \) – пространственный шаг.
Увеличение числа Куранта вручную позволяет ускорить процесс расчета задачи при большой серии вариантных оптимизационных расчётов. Данная опция показывает работоспособность не для всех видов задач. Однако, если удается произвести запуск с такими настройками, то получение сошедшегося решение значительно ускоряется. Для оптимальных модельных параметров с функцией решателя CPU Booster 100 оптимизационных шагов по результатам расчета заняли бы 11,5 часов.
Корректно подходя к задаче оптимизации проточной части — будь то рабочее колесо центробежного компрессора или ступень турбины, можно экономить вычислительные ресурсы и трудозатраты расчетчика. Это особенно актуально в настоящее время, так как перечень решаемых задач постоянно расширяется, а объемы вычислительных процессов неуклонно растут. Проведение не усложненных предварительных исследований позволяет выполнить в разы больше оптимизационных шагов или даже задач при тех же затрачиваемых вычислительных машинных ресурсах. А наработки по валидации расчетных моделей [14, 15, 16] позволяют более точно оценивать итоговые характеристики оптимизированной проточной части.
Заключение
Таким образом, в результате исследования выявлены наиболее значимые при проведении больших серий расчетов параметры расчетной модели рабочих колес центробежного компрессора – положение выходной границы и общая размерность сетки. Данные параметры вне зависимости от типа рабочего колеса при неоптимальном выборе увеличивают время получения сошедшегося решения в 2 и более раза, а также оказывают существенное влияние (до 1,7% КПД для выходной границы и 1,5% КПД для сеточной модели) на результаты решения. Это неприемлемо при оптимизации геометрической формы рабочих колес как с точки зрения получения аргументированного результата, так и с точки зрения рационального использования вычислительных ресурсов.
Кроме того, для адекватности получаемых результатов следует соблюдать сеточный параметр y+ для выбранной модели турбулентности и не допускать влияния на общие результаты расчётов и картину обтекания лопаток слишком близко расположенной к рабочему колесу входной границы. Слишком большие значения критерия расчетной сетки y+ оказывают наибольшее влияние на результаты решения. Из-за низкой разрешающей способности сетки в пограничном слое нивелируются отрывные течения и уменьшаются потери трения. Поэтому политропный КПД рабочего колеса может увеличиваться до 99% и более при неправильном выборе размера пристенной ячейки.
По результатам предварительной оценки определены оптимальные параметры расчетной модели при постановке задачи оптимизации. Для низконапорного колеса РК1 – выходная граница определена на 1,3D2; входная граница отодвинута от входа в рабочее колеса на расстояние 0,125D2; сеточная топология – HOH; размер первой пристенной ячейки – 0,002 мм; общий размер сеточной модели лопаточного сектора 0,9 млн. элементов. Для высоконапорного рабочего колеса РК2 – выходная граница определена на 1,25D2; входная граница отодвинута от входа в рабочее колеса на расстояние 0,125D2; сеточная топология – HOH; размер первой пристенной ячейки – 0,001мм; общий размер сеточной модели лопаточного сектора 0,7 млн. элементов. Данные наборы параметров модели обеспечивают наименьшее машинное время сходимости одного оптимизационного приближения.
СПИСОК ЛИТЕРАТУРЫ
- Kozhukhov Y.V., Danilishin A.M., Yun V.K. Multi-objective optimization for impeller shroud contour, width of vane diffuser and number of blades of the centrifugal compressor stage based on the CFD calculation. IOP Conference Series Materials Science and Engineering 08/2015; Volume 90(1):012046. DOI:10.1088/1757-899X/90/1/012046
- F.Frese, J. Einzinger, J.Will. Design optimization of an impeller with CFD and Meta-Model of optimal Prognosis (MoP). //10th International conference of turbocharges and turbocharging 2012, London. p. 121-135.
- R. Niemeier, S. Kunath, T. Most, J.Will, J. Einzinger. Robust Design Optimization and Operating Maps for Computational Fluid Dynamics. NAFEMS World Congress, San Diego 2015.
- A. Demeulenaere, J.-C. Bonaccorsi, D. Gutzwiller, L. Hu, H. Sun. Multi-Disciplinary Multi-Point Optimization of a Turbocharger Compressor Wheel. // ASME Turbo Expo 2015: Turbine Technical Conference and Exposition (Montreal, Canada)
- Dickmann H. P. 04/2013 Shroud contour optimization for a turbocharger centrifugal compressor trim family/ 10th European Turbomachinery Conference (Lappeenranta, Finland).
- Boldyrev, Y., Rubtsov, A., Kozhukhov, Y., Lebedev, A., Cheglakov, I., Danilishin, A. Simulation of unsteady processes in turbomachines based on nonlinear harmonic NLH-method with the use of supercomputers // CEUR Workshop Proceedings. Volume 1482, 2015, Pages 273-279
- В.Л. Блинов, Ю.М. Бродов, В.А. Седунин, О.В. Комаров. Выбор параметров расчетной модели при решении задач многокритериальной оптимизации плоских компрессорных решеток. // Компрессорная Техника и Пневматика, 2015. №1. С. 32-36
- Чеглаков И.В. Разработка методики проведения численного эксперимента и оптимизация неподвижных элементов малорасходных ступеней в среде моделирования Fine/Turbo на примере ступени СВД-22, дис. … магистра, СПбПУ, Санкт-Петербург, 2015.
- Батурин О.В. Методика цифрового моделирования осевых многоступенчатых турбин низкого давления с учётом неравномерного поля параметров на входе в турбину, трёхмерной структуры потока в лопаточных венцах и утечек через радиальные зазоры лабиринтных уплотнений. [Электронный ресурс]. Самара, 2012. 122 с.
- Свобода Д.А., Жарковский А.А. Экспериментальные и расчетные исследования осевого насоса с быстроходностью ns=560. // Известия Самарского научного центра Российской академии наук, 2013. т.15 №4(2). с. 573-578.
- Le Sausse, P., Fabrie, P., Arnou, D., & Clunet, F. (2013), CFD comparison with Centrifugal Compressor Measurements on a wide Operating Range, Johnson Controls Industries, In EPJ Web of Conferences, vol. 45, p. 01059.
- V. V. N. K. Satish Koyyalamudi and Quamber H. Nagpurwala, “Stall Margin Improvement in a Centrifugal Compressor through Inducer Casing Treatment,” International Journal of Rotating Machinery, vol. 2016, Article ID 2371524, 19 pages, 2016. doi:10.1155/2016/2371524
- AGHAEI TOG, R. – TOUSI, A. M. – TOURANI, A. Comparison of turbulence methods in CFD analysis of compressible flows in radial turbomachines. In: Aircraft Engineering and Aerospace Technology. Vol. 80, Iss. 6, 2008. pp. 657-665. ISSN 1748-8842.
- Яблоков А.М., Кожухов Ю.В., Лебедев А.А. Исследование течения в малорасходной ступени центробежного компрессора методами вычислительной газодинамики // Научно-технические ведомости СПбГПУ. Вып. 4(231)/2015. Изд. СПбГПУ. С. 59 – 69.
- Kozhukhov Y.V., Yun V.K., Reshetnikova L.V., Prokopovich M.V. Numerical Investigation of Different Radial Inlet Forms for Centrifugal Compressor and Influence of the Deflectors Number by Means of Computational Fluid Dynamics Methods with Computational Model Validation. IOP Conference Series Materials Science and Engineering 90(1):012047 • August 2015. DOI: 10.1088/1757-899X/90/1/012047
- KANG Shun, LIU Qiang, QI Ming-Xu .CFD validation of a high speed centrifugal compressor impeller[J] Journal of Engineering Thermophysics, 2005,V26(3): 400-404
REFERENCES
- Kozhukhov Y.V., Danilishin A.M., Yun V.K. Multi-objective optimization for impeller shroud contour, width of vane diffuser and number of blades of the centrifugal compressor stage based on the CFD calculation. IOP Conference Series Materials Science and Engineering 08/2015; Volume 90(1):012046. DOI:10.1088/1757-899X/90/1/012046
- F.Frese, J. Einzinger, J.Will. Design optimization of an impeller with CFD and Meta-Model of optimal Prognosis (MoP). //10th International conference of turbocharges and turbocharging 2012, London. p. 121-135.
- R. Niemeier, S. Kunath, T. Most, J.Will, J. Einzinger. Robust Design Optimization and Operating Maps for Computational Fluid Dynamics. NAFEMS World Congress, San Diego 2015.
- A. Demeulenaere, J.-C. Bonaccorsi, D. Gutzwiller, L. Hu, H. Sun. Multi-Disciplinary Multi-Point Optimization of a Turbocharger Compressor Wheel. // ASME Turbo Expo 2015: Turbine Technical Conference and Exposition (Montreal, Canada)
- Dickmann H. P. 04/2013 Shroud contour optimization for a turbocharger centrifugal compressor trim family/ 10th European Turbomachinery Conference (Lappeenranta, Finland).
- Boldyrev, Y., Rubtsov, A., Kozhukhov, Y., Lebedev, A., Cheglakov, I., Danilishin, A. Simulation of unsteady processes in turbomachines based on nonlinear harmonic NLH-method with the use of supercomputers // CEUR Workshop Proceedings. Volume 1482, 2015, Pages 273-279
- V.L. Blinov, Yu.M. Brodov, V.A. Sedunin, O.V. Komarov. Vybor parametrov raschetnoy modeli pri reshenii zadach mnogokriterialnoy optimizatsii ploskikh kompressornykh reshetok. // Kompressornaya Tekhnika i Pnevmatika, 2015. №1. S. 32-36
- Cheglakov I.V. Razrabotka metodiki provedeniya chislennogo eksperimenta i optimizatsiya nepodvizhnykh elementov maloraskhodnykh stupeney v srede modelirovaniya Fine/Turbo na primere stupeni SVD-22, dis. … magistra, SPbPU, Sankt-Peterburg, 2015.
- Baturin O.V. Metodika tsifrovogo modelirovaniya osevykh mnogostupenchatykh turbin nizkogo davleniya s uchetom neravnomernogo polya parametrov na vkhode v turbinu, trekhmernoy struktury potoka v lopatochnykh ventsakh i utechek cherez radialnyye zazory labirintnykh uplotneniy. [Elektronnyy resurs]. Samara, 2012. 122 s.
- Svoboda D.A., Zharkovskiy A.A. Eksperimentalnyye i raschetnyye issledovaniya osevogo nasosa s bystrokhodnostyu ns=560. // Izvestiya Samarskogo nauchnogo tsentra Rossiyskoy akademii nauk, 2013. t.15 №4(2). s. 573-578.
- Le Sausse, P., Fabrie, P., Arnou, D., & Clunet, F. (2013), CFD comparison with Centrifugal Compressor Measurements on a wide Operating Range, Johnson Controls Industries, In EPJ Web of Conferences, vol. 45, p. 01059.
- V. V. N. K. Satish Koyyalamudi and Quamber H. Nagpurwala, “Stall Margin Improvement in a Centrifugal Compressor through Inducer Casing Treatment,” International Journal of Rotating Machinery, vol. 2016, Article ID 2371524, 19 pages, 2016. doi:10.1155/2016/2371524
- AGHAEI TOG, R. – TOUSI, A. M. – TOURANI, A. Comparison of turbulence methods in CFD analysis of compressible flows in radial turbomachines. In: Aircraft Engineering and Aerospace Technology. Vol. 80, Iss. 6, 2008. pp. 657-665. ISSN 1748-8842.
- Yablokov A.M., Kozhukhov Yu.V., Lebedev A.A. Issledovaniye techeniya v maloraskhodnoy stupeni tsentrobezhnogo kompressora metodami vychislitelnoy gazodinamiki // Nauchno-tekhnicheskiye vedomosti SPbGPU. Vyp. 4(231)/2015. Izd. SPbGPU. S. 59 – 69.
- Kozhukhov Y.V., Yun V.K., Reshetnikova L.V., Prokopovich M.V. Numerical Investigation of Different Radial Inlet Forms for Centrifugal Compressor and Influence of the Deflectors Number by Means of Computational Fluid Dynamics Methods with Computational Model Validation. IOP Conference Series Materials Science and Engineering 90(1):012047 • August 2015. DOI: 10.1088/1757-899X/90/1/012047
- KANG Shun, LIU Qiang, QI Ming-Xu .CFD validation of a high speed centrifugal compressor impeller[J] Journal of Engineering Thermophysics, 2005,V26(3): 400-404
СВЕДЕНИЯ ОБ АВТОРАХ/AUTHORS
Неверов Владимир Валерьевич — Специалист, Санкт-Петербургский политехнический университет Петра Великого, без степени, без звания.
195251, Россия, г. Санкт-Петербург, Политехническая ул., 29.
E-mail: neverovpu@gmail.com
Neverov, V.V. Peter the Great St. Petersburg Polytechnic University.
29 Politechnicheskaya St., St. Petersburg, 195251, Russia.
E-mail: neverovpu@gmail.com
Кожухов Юрий Владимирович — Заведующий кафедрой, Санкт-Петербургский политехнический университет Петра Великого, канд. техн. наук, доц..
195251, Россия, г. Санкт-Петербург, Политехническая ул., 29.
E-mail: kozhukhov_yv@mail.ru
Kozhukhov, Y.V. Peter the Great St. Petersburg Polytechnic University.
29 Politechnicheskaya St., St. Petersburg, 195251, Russia.
E-mail: kozhukhov_yv@mail.ru
Яблоков Алексей Михайлович — ассистент, Санкт-Петербургский политехнический университет Петра Великого, без степени, без звания.
195251, Россия, г. Санкт-Петербург, Политехническая ул., 29.
E-mail: yablokovaleksey@mail.ru
Yablokov, Y.A.M. Peter the Great St. Petersburg Polytechnic University.
29 Politechnicheskaya St., St. Petersburg, 195251, Russia.
E-mail: yablokovaleksey@mail.ru
Лебедев Александр Анатольевич — доцент, Санкт-Петербургский политехнический университет Петра Великого, канд. техн. наук, без звания.
195251, Россия, г. Санкт-Петербург, Политехническая ул., 29.
E-mail: LAA7777@yandex.ru
Lebedev, A.L. Peter the Great St. Petersburg Polytechnic University.
29 Politechnicheskaya St., St. Petersburg, 195251, Russia.
E-mail: LAA7777@yandex.ru