FlowVision — современный инженерный инструмент в исследовании газодинамических характеристик компрессоров

Сергей Барашков, Владимир Шмелев

Подготовка геометрии компрессора в CAD-системе. Импорт модели в FlowVision

Работа в препроцессоре

Модели расчета

База данных веществ

Интенсивно развивающаяся CAE-система нового поколения FlowVision находит широкое применение как в научно-исследовательских работах, посвященных изучению вопросов динамики жидкости и газа, так и в сугубо практических приложениях — проектных работах, связанных с созданием устройств, где рабочим веществом является жидкость или газ.

Программный комплекс FlowVision в настоящее время широко используется для решения прикладных задач газовой динамики. В данной статье рассматривается опыт, полученный в ходе совместных работ ОАО ПКО «Теплообменник» и ООО «ТЕСИС», по применению комплекса для расчета и анализа газодинамических характеристик компрессоров.

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

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

Подготовка геометрии компрессора в CAD-системе. Импорт модели в FlowVision

Подготовка геометрии в CAD-системе

FlowVision не имеет собственного геометрического процессора, поэтому геометрия компрессора должна импортироваться из CAD-системы. Пример геометрии компрессора, созданного в CAD-системе Unigraphics, приведен на рис. 1.

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

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

Перед сохранением модель необходимо конвертировать в форматы фасеточного представления геометрии VRML (Virtual Reality Modeling Language) или STL (STereo Lithography format). При конвертации пользователь должен задать ряд параметров, которые определяют точность представления поверхности совокупностью фасеток (треугольников). Так, например, в Unigraphics процесс конвертации реализуется через диалоговое окно, приведенное на рис. 3.

Импорт в FlowVision

Перед импортом геометрии проточной части компрессора в FlowVision в диалоговом окне «Предустановки» необходимо задать настройки препроцессора, главной из которых является значение допуска (рис. 4). Значение допуска должно соответствовать допуску (толерантности), с которым была создана в CAD-системе твердотельная модель.

Импорт геометрии в FlowVision задается командой «Создать» из меню «Файл». После этого геометрия модели из выбранного файла будет импортирована в расчетный вариант, проанализирована, разбита на поверхности и замкнутые подобласти. Структура нового варианта появится в окне препроцессора, а изображение геометрии — в графическом окне.

Если в процессе импортирования были обнаружены ошибки в геометрической модели, то они будут показаны в препроцессоре и графическом окне варианта красным цветом (рис. 5).

К ошибкам приводят ситуации, когда в модели присутствуют ребра, принадлежащие меньше или больше чем двум фасеткам, что приводит либо к открытым граням, либо к «многолистности». Вариант с ошибками в модели рассчитывать нельзя. Следует поправить геометрию, а затем вновь произвести ее импортирование.

Мощным инструментом по исправлению возникающих при экспорте из CAD-систем ошибок является входящий в состав программного комплекса модуль Flow3DVision. Не вдаваясь в подробное описание этого модуля, отметим, что Flow3DVision может исправлять ошибки в параметрических ( igs, prt, x_t ) и фасеточных ( wrl, stl ) форматах. Откорректированная модель при этом сохраняется в специальном формате, доступном для чтения из FlowVision.

Работа в препроцессоре

Подготовка задачи к решению

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

Первая область, в нашем случае имеющая название «Подобласть#1», соответствует объему колеса компрессора. Вторая, с названием «Ротор», — проточной части между корпусом и колесом, включая входной патрубок. В третью, под названием «Статор», входит проточная часть улитки с выходным патрубком.

Определив расчетные области, переходим к постановке задачи, которая состоит из нескольких этапов.

Выбор модели расчета

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

Применительно к нашей задаче область с названием «Подобласть#1» оставляем нерасчетной. Областям с названием «Ротор» и «Статор» ставим в соответствие модель расчета с названием «Полностью сжимаемая жидкость». Выбор этой модели в нашем случае обусловлен сверхзвуковым течением газа на периферии колеса из-за значительной скорости его вращения. Стандартный набор уравнений для рабочей области «Ротор» дополним «Моделью зазора». Эту модель мы применим к щели между торцевыми кромками лопаток колеса и корпусом компрессора.

Для выполнения описанных выше действий используем диалоговое окно «Выбор модели» (рис. 7), которое вызывается командой «Изменить модель…» из всплывающего меню каждой расчетной области.

Задание физических параметров

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

Далее в диалоговых окнах «Начальные значения» и «Параметры модели» определяются соответствующие величины. В окне «Параметры модели» на вкладках «k-e параметры» и «Сопротивление щели» выбираются модели турбулентности (рис. 9) и зазора и задаются их параметры. В этом же окне на вкладке «Общие» пользователь уточняет диапазоны изменения рассчитываемых переменных.

Окно «Начальные значения» позволяет определить в каждой расчетной области значения давления и температуры в приращении относительно опорных величин, а также характеристики турбулентности (пульсация и масштаб турбулентности) и значения компонентов скоростей по осям координат в начальный момент времени (рис. 10).

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

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

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

Задание характеристик вращения колеса

Расчет движения газа относительно вращающегося колеса в FlowVision проводится во вращающейся вместе с колесом расчетной сетке. Для задания характеристик вращения необходимо во всплывающем меню опции «Движение» расчетной области «Ротор» выбрать команду «Свойства» и в открывшемся диалоговом окне (рис. 13) задать:

• тип движения — «Вращение»;

• направление вектора угловой скорости («Ось»);

• координаты точки, находящейся на оси вращения («Точка»);

• величину угловой скорости («Скорость»);

• начальный угол поворота («Угол»).

Дополнительно установим флажок в окошке «Использовать абсолютную координатную систему», что дает возможность вести описание задачи в неподвижной системе координат.

Определение и задание граничных условий

На этом этапе пользователь определяет типы граничных условий и задает их на конкретных поверхностях. В нашей задаче в обеих расчетных областях определим четыре типа граничных условий:

• вход — входное сечение входного патрубка;

• выход — выходное сечение патрубка улитки;

• скользящая поверхность — граница между расчетными областями «Ротор» и «Статор»;

• стенка — колесо и корпус компрессора.

В опции «Граничные условия» препроцессора каждой расчетной области, используя контекстное меню опций, создаем указанные типы граничных условий. Далее, вызывая диалоговое окно каждого граничного условия через его контекстное меню, задаем его характеристики.

В качестве граничных условий на входном сечении патрубка компрессора (рис. 14) введем условие равенства абсолютных значений температуры и давления воздуха опорным величинам (равенство нулю приращений), а турбулентной энергии и масштаба турбулентности — их начальным значениям в расчетных областях.

В качестве граничных условий на выходе из улитки компрессора возьмем «Нормальную массовую скорость», определяемую как отношение заданного расхода воздуха к его площади (рис. 15). Отрицательный знак величины при этом означает, что воздух через данное граничное сечение выходит из расчетной области. Температура, турбулентная энергия и масштаб турбулентности при этом определим как «Нулевой поток», что означает обнуление градиента величины по нормали к поверхности границы.

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

Поскольку проточная часть компрессора была разделена на две расчетные области, которые к тому же были разнесены в пространстве, в каждой области появилась цилиндрическая поверхность. На этих поверхностях решения из обеих областей должны совпадать. Применим граничное условие «Скользящая поверхность», цель которого и заключается в связывании расчетных областей, к выходной скользящей поверхности области «Ротор» и входной скользящей поверхности «Статор». Установление связи между расчетными областями проведем, используя команду «Создать» контекстного меню опции «Связать граничные условия» препроцессора. Результат создания этой связи закрепляется появлением компонента «Связь ротора и статора» в вышеуказанной опции. Свойства этого компонента (рис. 16) показывают, какие расчетные области связаны, через какие поверхности и какова вычисленная площадь последних.

Граничное условие «Стенка», которое мы применим к корпусу компрессора, отражает особенности движения газа около твердой поверхности. Параметры этого условия в программе могут быть уточнены пользователем как за счет применения условия «Стенка, логарифмический закон», задающего свойства пограничного слоя, так и за счет введения шероховатости поверхности, что также влияет на характеристики пограничного слоя. В нашем случае — на всей поверхности компрессора в обеих расчетных областях применим условие «Стенка, логарифмический закон» с нулевой шероховатостью (рис. 17).

Остановимся на особенностях применения граничного условия «Стенка» к поверхности компрессора в расчетной области «Ротор», учитывая, что описание задачи ведется в абсолютной системе координат, а также то, что в этой области применяется модель зазора. В диалоговом окне граничного условия «Стенка» в этом случае появляется дополнительная закладка «Зазор», отражающая применение в этой области специальной модели. В этой закладке пользователь для корпуса компрессора и выделенных в отдельную группу торцевых кромок лопаток колеса выбирает вариант «Зазор» (рис. 18). Описание задачи ведется в абсолютной системе координат, что задается в закладке «Скорость» варианта «Вращающаяся стенка, логарифмический закон» для всех вращающихся поверхностей.

Продолжение следует

«САПР и графика» 12'2004