11 - 2004

Применение программного комплекса FlowVision для проектирования авиакосмических конструкций

А.А.Аксенов, В.Н.Коньшин

Геометрическая модель и импорт данных

Расчетная сетка

Подсеточное разрешение геометрии

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

Моделированиe течения вокруг крылатой ракеты

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

FlowVision применяется для решения следующих задач, возникающих при создании авиакосмических систем:

• обтекание ракетоносителя при старте и на этапе прохождения атмосферы;

• сверхзвуковое обтекание спускаемого аппарата при входе в верхние слои атмосферы;

• движение газов и жидкостей по магистралям;

• горение топлива в камере сгорания двигателя;

• течение в дозвуковой и сверхзвуковой частях сопла;

• теплообмен в системе охлаждения ракетного двигателя;

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

• вентиляция внутренних отсеков объектов;

• отделение головной части космического корабля.

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

Геометрическая модель и импорт данных

Геометрическая модель объекта, созданная в какой-либо системе САПР (например, в Unigraphics, Pro /ENGINEER, CATIA, SolidWorks и т.д.), импортируется во FlowVision с помощью нейтральных форматов передачи данных, таких как STL, VRML, IGES, Parasolid и т.д. (рис. 1).

В начало В начало

Расчетная сетка

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

В начало В начало

Подсеточное разрешение геометрии

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

В начало В начало

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

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

• полные уравнения Навье-Стокса (закон сохранения импульса);

• уравнение неразрывности (закон сохранения массы жидкости);

• закон сохранения энергии, записанный через полную энтальпию;

• уравнение состояния идеального газа;

• k-e-модель турбулентности.

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

Использование неявных конечно-разностных схем приводит к необходимости решения систем линейных алгебраических уравнений с большими разреженными матрицами нерегулярной структуры вследствие применения локально адаптивных схем аппроксимации. Данные матрицы имеют три различные структуры — симметричные положительные матрицы, несимметричные матрицы и блочные несимметричные матрицы. Матрицы первой структуры появляются при аппроксимации эллиптических уравнений, матрицы второй структуры — из аппроксимации скалярных параболических уравнений, а матрицы третьей структуры — из векторных параболических уравнений (например, векторных уравнений Навье-Стокса). Для решения систем алгебраических уравнений с матрицей первой структуры применяется метод сопряженных градиентов с предобуславливателем типа неполного разложения Холесского, а для решения задач с матрицами второй и третьей структуры используется метод Ланцоша с предобуславливателями типа неполного LU- и неполного блочного LU-разложения. Эти итерационные методы обладают высокой скоростью сходимости к решению и допускают параллельную реализацию для многопроцессорных вычислительных систем.

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

В начало В начало

Моделированиe течения вокруг крылатой ракеты

В качестве примера рассмотрим обтекание сверхзвуковым потоком с числом Маха 2,95 крылатой ракеты с прямоточным двигателем, а также течение внутри воздухозаборников с учетом отражения ударной волны от стенок.

Анализ сходимости результатов расчета в зависимости от количества расчетных ячеек приведен на рис. 4. Расчетная сетка адаптировалась внутри воздухозаборников и вокруг корпуса ракеты. Видно, что при использовании более 250 тыс. расчетных ячеек достигается удовлетворительная сходимость по коэффициенту подъемной силы.

На рис. 5 и 6 показано распределение числа Маха в непосредственной близости от ракеты, а также распределение давления, в том числе по поверхности ракеты, для угла атаки 4°.

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

На рис. 8 представлены зависимости аэродинамических коэффициентов подъемной силы и продольного момента ракеты от угла атаки. Из рисунка видно, что увеличение угла атаки приводит к увеличению коэффициентов подъемной силы (ось Y направлена вниз по отношению к ракете, поэтому сила отрицательна) и продольного момента (момент определялся относительно носовой части ракеты).

В настоящее время программный комплекс FlowVision используют многие ведущие предприятия авиакосмической отрасли России и Украины: РКЦ «Энергия» (г.Королев), ЦНИИМаш (г.Королев), КБ «Энергомаш» (г.Химки), ФГУП ММПП «Салют» (г.Москва), ОАО «Авиационный комплекс им. С.В.Ильюшина» (г.Москва), ЛИИ им. М.М.Громова (г.Жуковский), ГКБ «Южное» (г.Днепропетровск). Компания «ТЕСИС» является разработчиком программного комплекса FlowVision и осуществляет все работы по системной интеграции, техническому обслуживанию комплекса, а также по обучению и консультированию пользователей.

А.А.Аксенов

Канд. физ.-мат. наук, начальник отдела вычислительной гидродинамики компании «ТЕСИС».

В.Н.Коньшин

Канд. физ.-мат. наук, менеджер отдела вычислительной гидродинамики компании «ТЕСИС».

В начало В начало

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