Авторефераты по всем темам  >>  Авторефераты по техническим специальностям

На правах рукописи

Волков Андрей Викторович

МЕТОД ЧИСЛЕННОГО ИССЛЕДОВАНИЯ ОБТЕКАНИЯ ПРОСТРАНСТВЕННЫХ КОНФИГУРАЦИЙ ПУТЁМ РЕШЕНИЯ УРАВНЕНИЙ НАВЬЕ-СТОКСА НА ОСНОВЕ СХЕМ ВЫСОКОГО ПОРЯДКА ТОЧНОСТИ

05.13.18. - Математическое моделирование, численные методы и комплексы программ

АВТОРЕФЕРАТ

Диссертации на соискание ученой степени доктора физико-математических наук

Москва - 2010

Работа выполнена в Федеральном государственном унитарном предприятии Центральном Аэрогидродинамическом институте им. проф. Н.Е. Жуковского (ФГУП ЦАГИ)

Официальные оппоненты: доктор физико-математических наук, профессор Толстых Андрей Игоревич, доктор физико-математических наук, в.н.с. Меньшов Игорь Станиславович, доктор технических наук, профессор Крицкий Борис Сергеевич Научный консультант - доктор физико-математических наук, Чернышев Сергей Леонидович Ведущая организация - ФГУП Центральный институт авиационного моторостроения им. П.И.Баранова г.Москва

Защита диссертации состоится 3 июня в 15.00 на заседании диссертационного совета Д 215.001.01 в Военно-воздушной инженерной академии имени профессора Н.Е.Жуковского по адресу 125190, г.Москва, ул.Планетная, д.3.

С диссертацией можно ознакомиться в библиотеке Военновоздушной инженерной академии имени профессора Н.Е.Жуковского

Автореферат разослан 2010 г.

Ученый секретарь Диссертационного совета Д 215.001.Кандидат физико-математических наук А.С. Ненашев

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность проблемы.

Решение проблем повышения качества и сокращения сроков проектирования новых ЛА всё в большей мере связывается с совершенствованием численных методов расчёта их аэродинамических характеристик. Модель течения, основанная на полных уравнениях НавьеЦСтокса или Рейнольдса, позволяет обеспечить новый уровень проектирования. В последние полтора десятилетия было разработано большое количество подходов к решению этих уравнений. Наиболее удачные численные схемы, реализованные в известных коммерческих вычислительных программах (FLUENT, CFX, STAR-CD, NUMECA), получили широкое распространение и относительно успешно используются при решении многих прикладных задач. В подавляющем большинстве современные программы базируются на методах конечных разностей (МКР) или конечного объёма (МКО). Для решения задач аэродинамического проектирования пространственных конфигураций ЛА современные расчётные схемы, базирующиеся на МКР или МКО, требуют использования значительных компьютерных и человеческих ресурсов, что приводит к заметному удорожанию результатов вычислений или, в некоторых случаях, к невозможности проведения численного эксперимента. Точность повсеместно используемых численных схем, как правило, не превышает второй порядок. Известно, что использование численных схем высокого порядка точности (выше второго) заметно снижает потребные вычислительные ресурсы и значительно расширяет круг решаемых задач. Однако, практическая реализация таких высокоточных схем встречает ряд принципиальных трудностей, решению которых посвящена данная диссертация.

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

Таким образом, актуальность работы определяется развитием авиационной техники и необходимостью разработки методов, направленных на повышение качества и сокращения сроков проектирования ЛА. Этой цели служит развитие численных схем высокого порядка точности, которые позволяют, с одной стороны, существенно уменьшить необходимые вычислительные ресурсы (память, время) по сравнению с существующими подходами и, с другой стороны, создают возможность построения алгоритмов, основанных на оптимальном автоматическом построении расчётных сеток с анизотропной адаптацией к особенностям течения.

Практическая значимость работы заключается в разработке эффективного сеточного метода решения уравнений Эйлера, НавьеЦ Стокса и Рейнольдса и в наборе практических рекомендаций по его применению. Созданные в процессе работы вычислительные программы используются при проведении расчётных исследований, как отдельных элементов ЛА, так и пространственной конфигурации в целом.

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

Научная новизна работы заключается в адаптации метода конечного элемента (МКЭ) к решению уравнений движения вязкого газа.

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

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

В качестве эффективного решателя системы сеточных уравнений, полученных при аппроксимации РМГ, впервые предложено использовать полиномиальный многосеточный подход с локально-неявным сглаживателем.

Метод РМГ высокого порядка точности для решения пространственных задач аэродинамики впервые практически реализован на гексаэдральных неструктурированных сетках.

Выполнено непосредственное сравнение требуемых компьютерных ресурсов при выполнении одних и тех же расчётов методами МКО и РМГ.

Автор защищает следующие результаты:

1. Метод численного решения пространственных уравнений Эйлера, Навье-Стокса и Рейнольдса при до и околозвуковых скоростях течения. Предложенный метод основан на конечно-элементном подходе и может быть использован для геометрически сложных областей течения с использованием концепции неструктурированных адаптивных сеток.

Метод использует меньшие компьютерные ресурсы по сравнению с традиционным подходом, основанным на методе конечного объёма второго порядка точности.

2. Универсальный способ устранения нефизичных осцилляций численного решения в областях течений с большими локальными неоднородностями. Способ основан на гибридном подходе, объединяющем немонотонную схему высокого порядка в основной области течения с монотонной схемой первого порядка точности в области локальной неоднородности.

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

4. Метод численного решения систем сеточных уравнений конечно-элементной аппроксимации течений, основанный на многоуровневом подходе с использованием различного набора полиномиальных базисных функций в реконструкции решения.

5. Результаты численных исследований течений вязкой жидкости и газа в модельных задачах аэродинамики и акустики, а также результаты расчётов обтекания компоновок ЛА.

Апробация работы. Разработанный метод прошёл тщательное тестирование путём сравнения результатов расчёта с имеющимися экспериментальными данными и/или аналитическими решениями.

Основные результаты диссертации содержатся в 12 статьях, опубликованных в ведущих рецензируемых научных журналах, представленных в перечне ВАК. Результаты докладывались на многочисленных международных и всероссийских научно-технических конференциях, в том числе: на Международном конгрессе по авиационным наукам ICAS (1996), 3-ей Европейской конференции по механике жидкости EUROMECH (1997), 6-ой международной конференции по генерации сеток для вычислительной аэродинамики (1998), 16-м Конгрессе Международной ассоциации математического и компьютерного моделирования IMACS (2000), конференциях Американского Института Авиации и Аэронавтики AIAA (2007, 2009), международной конференции по высокоскоростным течениям WEHSFF (2007), 5-ом Европейском конгрессе по численным методам в прикладной науке ECCOMAS (2008), франко-российском семинаре ONERA-ЦАГИ (2009), международной школе-семинаре Модели и методы аэродинамики (2004, 2006) и на 12-ти школах-семинарах ЦАГИ Аэродинамика летательных аппаратов (1998 - 2009).

В данную диссертацию включены результаты исследований, поддержанные РФФИ (Проекты №96-01-00209-л, № 96-01-00210-л, № 0001-00070-а, № 03-01-00236-а, №06-01-00283-а, №09-01-00243-а).

Объем работы. Диссертация состоит из введения, четырёх глав, заключения, списка литературы, включающего 156 наименований, и приложений. Общий объем - 189 страниц.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

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

Повышение порядка точности МКО приводит к необходимости расширения шаблона аппроксимации и к требованию повышения качества расчётных сеток. Таким образом, одно из решений проблемы эффективности численных схем связано с использованием схем высокого порядка точности, позволяющих выполнять расчёты на анизотропных адаптивных сетках. Схемы с высоким порядком дискретизации исследовались на компактном (Толстых 1990) и расширенном (Visbal & Gaitonde 2002) шаблонах с использованием стандартного метода конечных разностей. Для получения монотонного решения с высоким порядком аппроксимации в рамках МКР и МКО были разработаны также методы ENO/WENO (Harten et al. 1987; Liu et al. 1994;

Jiang & Shu 1996). Все эти подходы имеют солидный теоретический фундамент, но, к сожалению, применяются в основном на структурированных и неструктурированных сетках с почти равномерными распределением размеров и форм соседних ячеек.

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

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

Одним из наиболее перспективных подходов к высокоточной аппроксимации на базе МКЭ является метод Галеркина с разрывными базисными функциями (РМГ). В последние годы этот метод вызывает повышенный интерес многих исследователей вследствие его общности, гибкости и надежной теоретической обоснованности. Впервые метод был предложен Reed & Hill (1973) для решения уравнения, описывающего перенос нейтронов, а первый анализ был дан в работе Saint & Raviart (1974). Численное решение 2-D уравнений Эйлера и НавьеЦСтокса на треугольных неструктурированных сетках этим методом впервые представлено в работах Bassi & Rebay (1977). Наиболее полное теоретическое описание метода с решением 1-D и 2-D модельных задач приведено в работах Cockburn & Shu (1998, 2001). Подробное исследование РМГ можно также найти в работах Lyapunov & Wolkov (2000), Волкова и Ляпунова (2006, 2007, 2009) и в диссертации Ляпунова (2008).

Практическая реализация РМГ для решения нелинейных уравнений Эйлера, НавьеЦСтокса и Рейнольдса встречается с принципиальными трудностями, по причине которых этот метод ещё не получил широкого распространения. В настоящей диссертационной работе было выделено несколько основных проблем, решение которых и позволило успешно применить РМГ для расчёта пространственных конфигураций:

1) большое количество арифметических операций на одну степень свободы;

В отличие от МКО, в РМГ необходимо точное интегрирование не только потоков через грани контрольного объёма, но и вычисление объёмных интегралов внутри ячеек. Традиционный подход к интегрированию базируется на использовании квадратурных точек Гаусса. Суммирование значений подынтегральных функций в этих точках с определёнными весами обеспечивает получение точного значения в случае полиномиального представления подынтегральной функции. Чем выше порядок полинома, тем большее число гауссовых точек необходимо использовать. Присутствие объёмных и поверхностных интегралов в рассматриваемом методе приводит к необходимости использования чрезмерно большого количества этих точек на гранях контрольного объёма и внутри него.

2) немонотонность схемы высокого порядка в областях разрывов решений;

Решение задач обтекания может характеризоваться наличием в поле течения областей с большими локальными возмущениями, в частности, скачков уплотнения. Численное решение таких задач с использованием схем повышенного порядка точности приводит к нефизичным осцилляциям решения. В МКО для устранения нефизичных осцилляций используют ограничители градиентов решения, определяемые по значениям в ячейке и её соседях. Такой способ не применим в РМГ для решения задач неявным методом.

3) постановка граничных условий высокого порядка точности;

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

Порядок полинома, реконструирующего решение в ячейке, должен согласовываться с порядком кривой, её ограничивающей.

В традиционном МКО для получения второго порядка точности используют кусочно-линейное описание границы. Для решения уравнений НавьеЦСтокса при больших числах Рейнольдса требуется использование сеток с вытянутыми ячейками вдоль обтекаемого контура.

При этом реальная граница может пересекать ряд приграничных ячеек.

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

4) проблема решения больших нелинейных систем сеточных уравнений;

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

В настоящее время число 3-D реализаций РМГ для нелинейных уравнений законов сохранения весьма ограничено. Одна из первых успешных реализаций РМГ на неструктурированных гексаэдральных сетках представлена в публикациях автора диссертации.

В Главе 1 описаны известные сеточные методы аппроксимации уравнений Навье-Стокса и проанализированы их преимущества и недостатки. Обсуждаются некоторые технические подробности практического применения РМГ. Предложено два основных пути сокращения арифметических операций РМГ в расчёте на одну степень свободы: использование аналитического способа интегрирования потоков на низких многосеточных уровнях решателя и комплексная оптимизация квадратурных точек, используемых для расчёта поверхностных и объёмных интегралов.

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

В зз 1.3 и 1.4 первой Главы рассмотрена общая теория РМГ.

Предполагается, что вся расчётная область разбита на контрольные объёмы (элементы, ячейки). Решение в каждой ячейке хранится в примитивных переменных Q = , u, v, w, p. Система уравнений ( ) Рейнольдса представляется в примитивных переменных следующим образом:

Q + F U Q Fv U Q,U - S = 0, (1) ( ( ( ))- ( ( ) ) ) t здесь U = , u, v, w, E - вектор консервативных переменных; - ( ) плотность, u, v, w - компоненты вектора скорости течения, p - давление, E - удельная полная энергия; F U - невязкий, а F U,U вязкий ( ) ( )Ц v потоки; S - источниковый член, возникающий при использовании модели U турбулентности, матрица = есть якобиан преобразования от Q консервативных переменных к примитивным. В каждой ячейке сетки определяются локальные полиномиальные базисные функции j (x), линейная комбинация которых и определяет решение в ячейке:

K f Q(t, x) = (t) j (x) = u (t) (x), (2) u j j j j =где u (t) - вектор коэффициентов разложения, которые должны быть j определены в процессе решения, Kf - количество базисных функций в ячейке. Величина Kf связана с максимальной степенью базисного полинома K в пространственном случае через следующее соотношение:

K K = (K +1)(K + 2)(K + 3). При пространственной реализации РМГ ( ) f система базисных функций определена в физическом пространстве и имеет вид:

x ( - x0 y - y0 z - z) ( ) ( ) j x = ; j = 1,..., K ; 0 + + K, ( ) f h h h x y z где x0, y0, z0 - некая внутренняя точка элемента, а величины hx, hy, hz определяют размер ячейки вдоль соответствующих осей.

Система сеточных уравнений для коэффициентов u (t) из (2) j получается в соответствии со стандартной процедурой метода конечного элемента Галеркина, в котором требуется ортогональность невязки (левой части системы (1)) каждой базисной функции, используемой в реконструкции решения.

dui = -1M-1 i F - Fv n d + i F - Fv d + i S d. (3) ( ) ( ) - d t Здесь M - матрица интегралов от произведений различных комбинаций базисных функций, - контрольный объём, = - граница контрольного объёма. Уравнение (3) состоит из объёмных и поверхностных интегралов по границе ячейки. Значения зависимых переменных терпят разрыв на границе элементов, и ключевую роль здесь играют правила вычисления переменных и потоков на границах. Как и в МКО, в РМГ величина невязкого потока на грани между двумя ячейками определяется в результате решения задачи Римана о распаде произвольного разрыва. Проанализированы три приближённых подхода решения задачи Римана. Это линеаризованные методики Роу (Roe 1981), Ошера (Engquist & Osher 1981) и упрощённый вариант - схема ЛаксаФридрихcа (Lax 1954). Каждый из методов имеет свои преимущества и недостатки, а выбор конкретного подхода должен осуществляться в зависимости от решаемой задачи.

Правила аппроксимации градиентов примитивных переменных и вязких потоков в настоящей работе были реализованы в соответствии с рекомендациями работы Cockburn & Shu 1998.

В з1.5 обсуждаются особенности применения модели турбулентности Спаларта - Алмараса (Spalart & Allmaras 1992) для решения задач турбулентного обтекания. Выбор модели турбулентности был обусловлен кругом рассматриваемых задач - задач внешней аэродинамики при больших числах Рейнольдса. Отметим, что с точки зрения исследования эффективности использования схем высокого порядка точности, этот выбор не принципиален, и подобные исследования могли бы быть выполнены с любой другой моделью турбулентности. Добавление дополнительного уравнения, ответственного за модель турбулентности, значительно усложняет процесс получения общего решения, и эта сложность тем заметнее, чем выше предполагаемый порядок точности расчётной схемы. В диссертационной работе описаны некоторые технические детали применения этой модели, позволившие добиться надежной сходимости алгоритма поиска решения.

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

Следующий этап состоит в переводе этих фигур в канонические координаты, в которых положение их вершин строго детерминировано.

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

Интегрирование внутри элементарных фигур в канонических координатах осуществляется с использованием квадратурных формул Гаусса, которые сводят интегрирование к сумме значений подынтегральной функции в предписанных (гауссовых) точках с определёнными весами. Каждая квадратурная формула для определённой элементарной фигуры является уникальной и может быть найдена в специализированных справочниках (см., например, Strang & Fix, 1973), или в Интернете (например, www.cs.kuleuven.ac.be/~nines/research/ecf).

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

Основным критерием выбора явилось правило согласованности квадратурных формул по степени точности для контурных и объёмных интегралов. Это требование формализуется в следующем виде:

-rknrd = 0, r = x, y, z, k = 1, 2,3..., { } k rk d - где n, n,n - компоненты нормали к границе элемента. Данное x y z выражение является следствием теоремы о дивергенции.

Количество гауссовых точек в квадратурных формулах в основном и определяет общее количество арифметических операций всего метода.

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

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

Другой возможный подход к сокращению количества арифметических операций основан на аналитическом вычислении необходимых интегралов. В случае линейной зависимости потоков от зависимых переменных все интегралы сеточных уравнений (3) являются интегралами от полиномиальных функций, вычисления которых могут быть сделаны до начала основных итераций. Линейная связь потоков и зависимых переменных возникает, например, в задачах акустики, где решаются линеаризированные уравнения Эйлера: F U = f . Такое ( ) j j представление потоков позволяет свести их интегрирование к сумме интегралов от простых полиномиальных выражений:

i i j j j i j F(U)d = f d = f d, iF U n d = if nd = f ij n d.

( ) j j j Применение такого подхода возможно и в случае, когда поток является нелинейной функцией зависимой переменной. Для этих целей функция потока также может быть разложена в ряд по базисным M f функциям: F U = j. Здесь общее количество элементов ( ) f j j =разложения Mf может быть отлично от количества базисных функций в элементе. Коэффициенты разложения могут быть найдены по правилам нахождения коэффициентов в ряде Тейлора (с последующим переразложением по базисным функциям элемента). Как правило, для обеспечения необходимой точности число элементов разложения Mf должно превышать число базисных функций Kf в элементе. Однако число Mf можно уменьшить до Kf посредством решения следующей линейной системы уравнений относительно новых коэффициентов разложения f :

j K M f f f f f j j j j id = id, i = 1,..., K.

j=1 j= В реализованном РМГ на неструктурированных гексаэдральных сетках применение безквадратурного подхода оказалось эффективным лишь на грубых многосеточных уровнях, где основная цель решателя состоит в уменьшении низкочастотной составляющей ошибки решения на мелкой сетке. На верхнем многосеточном уровне, где используется наибольшее количество степеней свободы, интегрирование осуществлялось на базе квадратурного метода.

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

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

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

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

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

Нефизичные осцилляции приводят к невозможности получения решения на мелких и адаптированных сетках. Для их устранения используют либо добавление искусственной вязкости, либо применение монотонных схем первого порядка точности. Необходимость использования схем более высокого порядка точности в гладких областях течения обусловило создание так называемых гибридных схем или схем переменного порядка точности. Первая гибридная схема была предложена в работе Федоренко (1962), в которой также был описан алгоритм переключения между монотонной схемой и схемой высокого порядка. Другие гибридные схемы с плавным переходом от низкого порядка к высокому были предложены в работе Гольдин, Калиткин, Шишова (1965). Впоследствии, в рамках конечно-разностной идеологии построения численных схем, было разработано большое количество различных гибридных методов (обзор можно найти в книге Куликовский, Погорелов, Семенов, 2001).

В МКО для устранения нефизичных осцилляций в решении используют ограничители градиентов решения (лимитеры):

U t, x = U0 + grad U x - x0, предложенные Колганом (1972), а ( ) ( )( ) позднее Van Leer (1979). Для предотвращения осцилляций решения коэффициент уменьшается вплоть до нулевого значения. При этом восполнение решения в ячейке становится кусочно-постоянным, а точность схемы падает до первого порядка. Другой подход реализован в схемах ENO/WENO, в которых шаблон аппроксимации в процедуре реконструкции выбирается на основе локальной гладкости численного решения.

В з2.2 обсуждаются проблемы применения традиционных ограничителей решения в методе конечного элемента. Рассматривается u u одномерная схема РМГ для задачи переноса + = 0. В качестве t x набора базисных функций используются ортогональные полиномы Лежандра, определённые в промежутке от -1 до 1: 1 = 1, 2 = x, 3 = -1 Реконструкция решения в каждом элементе (3x ).

представляется выражением: U t, x = u11 x + u22 x + u33 x.

( ) ( ) ( ) ( ) Таким образом, коэффициент u1 является средним значением решения в ячейке, а u2 и u3 отражают вклад линейной и квадратичной составляющих соответственно. Традиционный ограничитель решения вводится в реконструкцию решения в виде коэффициента перед базисными функциями высокого порядка: U t, x = u11 x + u22 x + u33 x.

( ) ( ) ( ( ) ( )) При этом происходит ограничение вклада в решение составляющих высоких порядков, вплоть до использования кусочно-постоянного восполнения решения (при =0). Впервые в РМГ такой подход был применен в работе Cockburn, Hou, Shu (1990) при использовании явной схемы Рунге-Кутта интегрирования по времени. Главный недостаток такого подхода состоит в невозможности его применения в неявных схемах. Введение лимитеров, принимающих нулевые значения, приводит к вырождению матрицы Якоби для системы сеточных уравнений относительно неизвестных значений u1, u2 и u3. Действительно, система аппроксимирующих уравнений в ячейке представляется в виде системы нелинейных уравнений:

F1(u1, u2, u3) = + =.

F (u1, u2, u3) = i u u 0, i = 1,2,3 t x F (u1, u2, u3) = Здесь F1=F2=F3=0 - сеточные уравнения для определения u1, u2, u3 в рассматриваемой ячейке, выписанные в соответствии с алгоритмом РМГ.

Fi Очевидно, что при =0 Якобиан системы этих уравнений u j вырождается. Таким образом, устранение осцилляций решения возможно на основе использования других известных подходов, а именно, добавление искусственной вязкости и использование гибридной схемы.

Отметим, что метод добавления искусственной вязкости в РМГ опробован рядом исследователей Persson & Peraire (2006);

Barter & Darmofal (2007), в то время как разработка и исследование гибридной схемы впервые представлены в диссертации автора.

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

РМГ (K, )= РМГ (K) + (1- ) РМГ (0).

Здесь буквенная последовательность РМГ (K) обозначает классическую схему Галёркина с разрывными базисными функциями, описанную в з1.3. Особенность новой схемы заключается в возможности плавного регулирования порядка точности при плавном изменении коэффициента . В результате получается всюду дифференцируемая численная схема, с невырожденной матрицей Якоби, позволяющая непрерывно перейти от схемы высокого порядка к монотонной схеме низкого порядка точности.

Схему необходимо применять в областях разрывов решения и в местах неразрешенных особенностей течения.

Принципы построения конечно элементной схемы с плавным изменением порядка точности, описаны в з2.3 для случая использования сеток с треугольными элементами. Подчеркнём, однако, что форма элемента не принципиальна, и этот подход легко обобщается на элементы произвольной формы.

Для реализации РМГ (K, ) схемы каждый сеточный элемент разбивается на множество подэлементов: E = e1 e2 eNs, и с каждым подэлементом ассоциируется только одна узловая точка. Все базисные функции определяются в параметрическом пространстве для стандартного треугольника с вершинами P1(1, 0), P2(0, 1) и P3(0, 0).

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

Конкретные примеры разбиений треугольных элементов на подэлементы для случаев кусочно-линейной и кусочно-квадратичной реконструкций приведены на рис. 1. Здесь же приведены узловые точки, использованные при выборе базисных функций.

.

.

Y Y.

(0.5, 0.5) 6.

Y 6 Y YYY3 Y 4 Y...

3..1 б) а) Рис. 1. Разбиение треугольника в параметрическом пространстве для РМГ (K, ) схемы: а) - K=1; б) - K=В случае, когда максимальный порядок полинома в элементе ограничен первым порядком (кусочно-линейная реконструкция), треугольник разбивается на два треугольных и один четырехугольный элемент (Ns=3) (рис.1,a) и в данном элементе вводятся три линейные базисные функции:

1 = x, 2 = y, 3 = 1- x - y.

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

1 = -x + 2x2, 2 = - y + 2y2, 3 = 1- 3x - 3y + 2x2 + 2y2 + 4xy, 4 = 4x - 4x2 - 4xy, 5 = 4xy, 6 = 4y - 4y2 - 4xy.

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

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

Итоговые базисные функции схемы РМГ (K, ) определяются i как взвешенная сумма кусочно-линейных (или квадратичных) и кусочнопостоянных базисных функций: = i + 1- i. При этом решение в ( ) i каждом элементе E представляется как линейная комбинация базисных функций: U x = u x. Как было описано в з1.3, система сеточных ( ) ( ) j j уравнений для коэффициентов uj получается после умножения системы уравнений (3.1) на каждую из базисных функций и интегрирования по частям. Отметим, что базисные функции, а, следовательно, и потоки терпят разрыв на границах не только элементов, но и подэлементов. В связи с этим интегрирование должно выполняться по каждому подэлементу отдельно. Таким образом, для каждого элемента E получается следующая система уравнений:

K f d Ns u d + F - Fv n d ( ) i j j i ek ek j=1 = 0.

d t ek =- F - Fv d - Sd ( ) i i ek ek Свойства РМГ (K, ) схемы рассмотрены в з2.4 диссертации. Здесь на примерах решения задач конвекции и конвекцииЦдиффузии, имеющих гладкое решение, выполнены оценки порядка точности предложенной гибридной схемы. Показано, что сконструирована схема, в которой плавное изменение параметра приводит не к дискретному, а к плавному изменению порядка ошибки от монотонной схемы первого порядка точности при =0 до схемы с порядком K+1 при =1.

В з 2.5 описаны и проанализированы некоторые известные сенсоры осцилляций решения: Cockburn & Shu (1998), Krivodonova et al.

(2004), Persson & Peraire (2006), и предложен альтернативный сенсор областей сингулярности решения.

Новый сенсор основывается на разложении некой нелинейной комбинации примитивных переменных в ячейке по базисным функциям.

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

При реализации данного алгоритма в качестве нелинейной комбинации примитивных переменных использовалась величина полной энтальпии u2 + vтечения: H = E + -1. Коэффициенты ai первого разложения ( ) H x = aii x находятся в соответствии со стандартной процедурой ( ) ( ) - j МКЭ: ai = Mi,1 H d. Используя значения энтальпии H в Kf j j фиксированных точках элемента (его вершины, середины рёбер, центральная точка) можно получить другое разложение H x = bii x ( ) ( ) путём решения следующей системы линейных уравнений:

j bijij = H, j=1,Е, Kf.

Разница между этими коэффициентами разложения и составляет основу предложенного сенсора, а величина определяется из следующих соотношений:

1, если s2 < s0 - ( ) 1 s2 - s0 1+ 2 sin , = 2, s = g10 .

ai bi a1 если s0 - s2 s0 + 0, если s2 > s0 + Здесь s0 и - некоторые настраиваемые параметры, которые, как показала практика, достаточно консервативны по отношению к выбору решаемой задачи и типу расчетной сетки, и во всех выполненных расчётах использовались следующие параметры: s0=-1.5, =1.

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

Таким образом, настоящий алгоритм может быть применён в качестве индикатора зон скачкообразного изменения решения при расчёте РМГ (K, ) схемой на анизотропных адаптивных сетках.

В зз 2.6-2.8 приведены примеры использования РМГ (K, ) схемы на последовательности анизотропных адаптивных сеток при решении задач конвекции с разрывом, расчёте обтекания профиля транс- и сверхзвуковым потоком невязкого сжимаемого газа (решение уравнений Эйлера) и расчёте обтекания профиля трансзвуковым турбулентным потоком (решение уравнений Рейнольдса).

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

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

В з2.8.А представлена серия расчетов для известного тестового случая Case-6 Cook at. al. (1979) обтекания профиля RAE 2822 при числах M=0.725 и Re=6.5106. Экспериментальный угол атаки составляет 2.92, расчетный угол был выбран из условия равенства экспериментального и расчетного коэффициентов подъемной силы Cy=0.738. Расчеты производились на последовательности адаптивных сеток с использованием как кусочно-линейной РМГ (1, ), так и кусочноквадратичной РМГ (2, ) реконструкции решения в ячейке. Финальные расчетные сетки, полученные в процессе расчета с РМГ (1, ) и РМГ (2, ), представлены на рис. 2.

Со схемой второго порядка точности (K=1) было выполнено итераций адаптации, в то время как схема третьего порядка точности (K=2) приблизилась к такому же количеству степеней свободы за итераций (см. рис. 3). Результаты, полученные с использованием обеих схем, практически совпадают. Результат сравнения расчетного распределения давления с экспериментальными данными представлен на рис. 3.

Рис. 2. Финальная расчетная сетка, полученная в процессе расчета:

а) - РМГ (1, ); б) - РМГ (2, ).

Cy 0.8 -1.Cp -0.-0.0.РМГ К=РМГ К=0.0.0.0.103 104 10 0.2 0.4 0.6 0.8 X Число степеней свободы Рис.3. Рис.4.

Рис. 3. Зависимость коэффициента подъемной силы от количества степеней свободы расчетных схем РМГ (1, ), РМГ (2, ). Представлен расчет расчет турбулентного обтекания профиля RAE 2822 при M=0.725, Re=6.51Рис. 4. Распределение давления:

- экспериментальные результаты =2.92, Cy=0.738, Cx=0.0127;

ЧЧ - расчёт РМГ (1, ) =2.485, Cy=0.738, Cx=0.0151.

В з2.8.Б рассмотрено трансзвуковое отрывное обтекание профиля NACA64A010 при M=0.8, Re=2106, =6.2 (Johnson & Bachalo 1980).

В з2.9 приведён пример использования РМГ (K, ) схемы на неструктурированной гексаэдральной сетке для расчета крыла в трансзвуковом турбулентном потоке.

Представленные примеры решения тестовых задач демонстрируют возможности гибридной схемы РМГ высокого порядка точности подавлять нефизичные осцилляции на разрывах в решении при расчёте транс- и сверхзвуковых течений.

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

В з3.1 представлен явный алгоритм поиска решения, традиционно используемый в РМГ (см. Cockburn & Shu 1998). В частности, рассмотрено семейство так называемых оптимальных TVD схем РунгеКутта (Shu & Osher 1988), то есть схем с наибольшим допустимым числом Куранта, обеспечивающим при этом сохранение по времени полной вариации решения. В качестве примера использования явной TVD схемы интегрирования уравнений РМГ аппроксимации приводится решение задачи о распространении цилиндрической акустической волны малой амплитуды. Оценивается поле давления к моменту времени, когда акустическая волна подошла к границам расчётной области. Решение этой задачи требует одновременного использования схемы высокого порядка, как по времени, так и по пространству. Задача имеет аналитическое решение, что даёт возможность оценить реальный порядок точности разработанной численной схемы. Для реконструкции решения в ячейке использованы полиномиальные базисные функции четвёртого порядка (K=4). Интегрирование по времени выполнено шести шаговой TVD схемой Рунге-Кутта пятого порядка точности. Результаты, полученные на вложенной последовательности сеток, демонстрируют сходимость расчётного поля давления к аналитическому полю с порядком точности, превышающим пятый.

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

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

В диссертационной работе описываются некоторые подходы, позволяющие снизить использование требуемых вычислительных ресурсов. Неявные методы решения систем сеточных уравнений, полученных в результате аппроксимации методами МКР и МКО, получили широкое распространение (например, Venkatakrishnan 1988;

Егоров и Зайцев 1991). Эффективность неявного метода в рамках РМГ была апробирована на решении задачи расчёта многоэлементного профиля с элементами взлётно-посадочной механизации крыла на последовательности адаптивных сеток (Волков и Ляпунов 2006, 2007).

з3.3 посвящён описанию наиболее эффективного метода решения - многосеточного метода. В настоящее время в промышленных программах решения уравнений Навье-Стокса и Рейнольдса с целью ускорения итерационного процесса сходимости к стационарному решению используют многосеточный метод, предложенный Федоренко в 1961. Суть метода состоит в ускорении передачи информации между далеко расположенными друг от друга ячейками сетки. Это осуществляется через укрупнение сеточных ячеек путём их агломераций.

В РМГ такой подход малоэффективен, так как изначально предполагается использование грубых сеток с большим количеством переменных в каждой ячейке. В последнее время философия построения многосеточного метода была пересмотрена и расширена (Rnquist & Patera 1987). В работах Helenbrook et al. (2003), Жуков и др. (2004) метод был адаптирован к использованию в конечноэлементном методе аппроксимации течений. По аналогии с решением задачи на последовательности измельчённых сеток, в РМГ решается задача на фиксированной сетке, но с использованием различного набора базисных функций. В настоящей работе полиномиальный многосеточный метод впервые применён к решению пространственных уравнений НавьеСтокса в рамках РМГ. Предложены новые операторы перевода решения между различными уровнями многосеточного метода (операторы интерполяции и сборки), а также локально-неявный сглаживатель ошибки решения.

Предложен следующий алгоритм решателя, подробно описанный в з3.3.1. Определяется последовательность уровней многосеточного метода = 1,Е, L, такая, что уровень +1 в сравнении с уровнем имеет либо большее количество ячеек, либо большее количество базисных функций (степеней свободы). При традиционном сеточном измельчении (h-измельчение) используется последовательность вложенных гексаэдральных сеток, генерируемых алгоритмом агломерации соседних ячеек. При измельчении по степеням свободы (pизмельчение) используются наборы из базисных функций, описанных в таблице 1.

Таблица 1. Набор базисных функций в элементе K=0 K=1 K=2 K= Базис набора Базис наборов Базис наборов K=0 K=0, 1 K=0, 1, дополняется дополняется дополняется функциями: функциями: функциями:

x2, y2 z2, x3, y3, z3, 1 x, y, z xy, xz, yz, x2y, x2z, y2x, y2z, z2x, z2y, xyz Далее описание полиномиального алгоритма полностью совпадает с описанием традиционного многосеточного алгоритма по схеме полной аппроксимации (FAS), представленной в работе Brandt (1980).

На самом верхнем (мелком) многосеточном уровне L решается qL система уравнений Навье-Стокса или Рейнольдса: M + RL (qL ) = 0.

t Здесь q - искомый вектор решения, R(q) - вектор-функция невязки решаемых уравнений или результат пространственной аппроксимации законов сохранения без учёта производной решения по времени, M - блочно-диагональная матрица, состоящая из интегралов от произведений различных базисных функций элемента. Решение данной системы уравнений может быть выполнено явным или неявным способами, описанными выше. Эти способы - итерационные, и применение одной или нескольких итераций уменьшает лишь высокочастотную ошибку решения, но не позволяет решить задачу полностью. Поэтому данную операцию принято называть сглаживанием, а соответствующий алгоритм - сглаживателем. Применение сглаживателя позволяет получить коррекцию решения на верхнем слое L: qL = qL + qL.

На более грубом (расположенном ниже) многосеточном уровне = L -1 к решаемой системе уравнений добавляется правая часть:

q M + R (q ) = F, где F = R I +1 q +1 + I +1 F +1 - R +1(q +1).

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

q(new) = q + qCorr, qCorr = I -1 -1 - I -1 q.

[ ] { } q Здесь фигурные скобки означают осреднение значения в рассматриваемой ячейке и её соседях.

Как было указано выше, многосеточные уровни и -1 могут различаться либо разным количеством ячеек в сетке при одинаковом количестве базисных функций в реконструированном решении, либо разным порядком полиномиальной реконструкции, но на фиксированной сетке. В зависимости от этих обстоятельств варьируется и вид операторов переноса решения и невязки решения, подробно описанных в зз 3.3.1 и 3.3.2.

В случае перевода решения между уровнями с разным количеством сеточных ячеек используются традиционные операторы интерполяции и сборки. Для перевода решения между уровнями с разным количеством базисных функций в реконструируемом решении были разработаны новые операторы интерполяции I +1 и сборки I +1 :

-1 -I +1 q = Mpqq, I +1 q +1 = Mqpq +1, [ ] [ ] Mpp Mqq -I +1 R +1 = Mqp R +1.

[ ] Mpp Здесь M - массовые матрицы, каждый член которых определяется в соответствии со следующими правилами:

-1 -Mi, j = j-1d, Mi, j = jd, Mi, j = jd, Mi, j = j-1d.

qq i pp i qp i pq i Ключевым алгоритмом любого многосеточного метода является так называемый сглаживатель - алгоритм подавления низкочастотной составляющей ошибки решения.

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

В з 3.3.5 приведён пример оценки эффективности полиномиального многосеточного метода в задаче об обтекании сферы невязким потоком при числе M=0.15. Семейство неструктурированных гексаэдральных расчётных сеток было построено с использованием программы HEXPRESS (компании NUMECA) (см. рис. 5). Алгоритм агломерации позволил построить систему из трёх вложенных сеток, решения на которых получалось с использованием наименьшего, кусочно-постоянного базиса. Далее уровни многосеточного алгоритма измельчались путём обогащения базиса, как это было описано ранее.

Таким образом, расчёты РМГ c K=3 (кусочно-кубический базис) выполнялись с использованием шестиуровневого алгоритма. Типичная история сходимости (зависимость L2 нормы невязки решения от времени CPU) представлена на рис. 6 для случая кусочно-квадратичной аппроксимации решения. Рисунок 6 демонстрирует эффект ускорения решения при использовании многосеточного алгоритма. Отметим, что без применения многосеточного алгоритма при кусочноЦкубичном базисе получение решения за разумное время становится невозможным.

Рис. 5 Рис. Рис. 5. Фрагмент гексаэдральной сетки около сферы.

Рис. 6. История сходимости. Эффект ускорения решения при использовании многосеточного алгоритма.

В Главе 4 проведены численные исследования точности новой схемы аппроксимации на последовательности вложенных сеток. Дано сравнение результатов расчёта и потребных вычислительных ресурсов по предложенной схеме высокого порядка точности с теоретическими и (или) экспериментальными результатами и результатами, полученными с использованием стандартной схемы МКО второго порядка точности.

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

Известно, что наличие скачка уплотнения ограничивает порядок точности численной схемы до первого (см., например, Иванов & Крайко 1978). При этом возникает очевидный вопрос о целесообразности применения схем высокого порядка для решения задач с разрывом.

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

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

p u x un x = O hск + O hгл.

( )- ( ) ( ) ( ) h n n Здесь u x - аналитическое решение, un x - численное решение, ( ) ( ) h полученное на сетке n из вложенной последовательности сеток, hск - n размер ячеек в области скачка уплотнения, hгл - размер ячеек в гладкой n области течения, p - порядок точности численной схемы. Используя измельчённую сетку в окрестности скачка уплотнения, величину ошибки можно сделать сравнимой с ошибкой в гладкой области течения, где используется схема с высоким порядком точности, но на сетке с крупными ячейками. Такой подход приводит к существенной экономии количества расчётных узлов и, как следствие, к уменьшению компьютерных ресурсов.

В з 4.2 приводится оценка фактического порядка точности РМГ в задаче о ламинарном обтекании кругового цилиндра при M=0.15 и Re=40. Рассмотрены две вложенных последовательности из 4 сеток.

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

Другая последовательность сеток состояла из ячеек сложной формы, которые образовывались путем искусственной агломерации некоторых соседних треугольных элементов. Фрагменты таких сеток приведены на рис.7.

Предельное значение сопротивления, экстраполированное на бесконечно мелкую сетку, составило 1.5041881. Оценки порядков ошибки выполнены при использовании базисных функций третьего РМГ(K=3) и четвёртого РМГ(K=4) порядков. Величины сопротивления Cx и порядка ошибки p на последовательности из сеток треугольных ячеек и ячеек сложной формы отражены в таблице 2. Из представленных результатов видно, что, во-первых, несмотря на большие значения ошибок, полученных на сетках с агломерированными ячейками, порядок ошибки близок к порядку, полученному на треугольных сетках. Кроме того, для всех случаев полученные значения близки к теоретическому результату: K+1.

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

а) сетка №-1: 10x5, б) сетка №-4: 80xТаблица 2. Порядки ошибок в задаче о ламинарном обтекании цилиндра Треугольные ячейки Ячейки сложной формы K=3 K=4 K=3 K=Cx p Cx p Cx p Cx p 1.5344122 1.4905387 2.1969095 1.9978410x1.5059944 4.06 1.5045982 5.06 1.5371651 3.75 1.5152094 4.20x1.5042885 4.17 1.5041962 5.66 1.5053013 4.69 1.5045512 4.40x1.5041953 3.80 1.5041883 5.40 1.5042349 4.4 1.5041995 4.80xАналогичные результаты получены и в задаче о турбулентном обтекании пластины при числах M=0.1 и Re=1106. Так же как и в предыдущем случае, порядок точности оценивался по величине сопротивления. С целью исключения сингулярной особенности течения в передней части пластины сопротивление вычислялось на серединном участке пластины. При использовании кусочно-линейной и квадратичной реконструкций решения в ячейках в данном численном эксперименте также удалось получить оптимальный порядок точности O[h(K+1)].

В з4.5 описана серия расчетов обтекания сферы дозвуковым потоком газа при M=0.15, которая была выполнена на последовательности вложенных сеток с использованием кусочнолинейного, квадратичного и кубичного полиномиального базисов. Все расчёты с высоким порядком точности (K > 2) выполнялись с учётом кривизны граней гексаэдров, непосредственно примыкающих к сфере. На рис. 8 представлена зависимость величины коэффициента сопротивления от числа переменных задачи. Известно, что в невязком течении величина сопротивления должна быть равна нулю, и, таким образом, величина отличия расчётного коэффициента сопротивления от нуля косвенно характеризует качество расчётной схемы. Анализ результатов показал, что скорость стремления к нулю коэффициента сопротивления по мере сгущения сеток выше, чем 2K+1, что, возможно, говорит о достижении суперсходимости (см. Barth 2004) конечно-элементного метода. На самой густой из использованных сеток, содержащей 8 148 ячеек, при применении кусочно-кубичного базиса для реконструкции решения коэффициент сопротивления четверти сферы составил величину 3.210-6.

Получить такое малое значение коэффициента сопротивления методом конечного объёма чрезвычайно затруднительно. Например, использование кода FLUENT, даже на сетке в 240 000 узлов, не позволило получить сопротивление меньше, чем 2.110-4.

|Cx| 10-10-10-РМГ (K=1) РМГ (K=2) 10-РМГ (K=3) 10-10-4.

.104.

5 10 10 5 2 N Рис. 8. Зависимость величины коэффициента сопротивления сферы от числа переменных в РМГ для кусочно-линейного, квадратичного и кубичного базисов В з4.4 представлено сравнение расчётов ламинарного обтекания пластины (M=0.35, Re=76000) методами РМГ и МКО в условиях эквивалентного количества степеней свободы. Расчёты выполнялись на гексаэдральных сетках разной густоты, имеющих только одну ячейку в zнаправлении. Некоторые параметры сеток, представленные в таблице 3, показывают: (1) - общее количество ячеек в сетке, (2) - расстояние от поверхности пластины до первого слоя ячеек, (3) - коэффициент увеличения расстояния между последовательными рядами сеточных точек и общее количество узлов (4) - в x и (5) - в y-направлениях.

Таблица 3. Параметры расчётных сеток Сетка (1) (2) (3) (4) (5) №1 6739 110-5 1.02 25 1№2 281 410-5 2.50 5 Отметим, что количество элементов в сетке №1 более чем в 20 раз превышает количество элементов в сетке №2. Такой выбор обеспечивает возможность проводить сравнения численных схем в условиях близкого количества степеней свободы. На рис. 9 представлены расчётные профили скорости и даны сравнения с аналитическим решением Блазиуса. Расчётные профили изображены комбинацией линий и маркеров. Каждая линия построена через множество точек в точном соответствии с формулой для реконструкции решения в ячейке. Маркеры же на каждой линии расположены приблизительно в серединах ячеек.

Таким образом, количество маркеров на линии отражает густоту сетки.

U U V Re Rex U U 1.а) б) 0.0.0.0.Решение Блазиуса 0.Сетка #1, МКО, =67N 0.Сетка #2, РМГ(K=3), N =281x20=560.0.0 2 4 6 0 2 4 6 Re Re Re Re Re y y y y y x x x x x Рис. 9. Сравнение U и V компонент профилей скорости с аналитическим решением Блазиуса (зелёная линия). Частота маркеров отражает густоту сетки Продольная U-компонента скорости, показанная на рис. 9.а), хорошо согласуется с аналитическим решением для обеих схем.

Поперечная V-компонента оказывается более чувствительной к качеству сетки (рис. 9.б). Результаты, полученные с РМГ K=3 на сетке №2, лучше согласуются с аналитическим решением, чем результаты, полученные традиционным МКО на самой мелкой сетке №1, несмотря на эквивалентное количество задействованных степеней свободы. Отметим, что осциллирующее поведение расчётной кривой, полученной МКО, является результатом использования недостаточно мелкой сетки.

Сравнение расчётов трёхмерного ламинарного течения в изогнутой трубке методами РМГ и МКО в условиях эквивалентного количества степеней свободы представлено в з4.6. Пространственное ламинарное течение воды внутри изогнутой под прямым углом трубки постоянного сечения было экспериментально исследовано Enayet et al.

(1982) с помощью лазеро-доплеровского измерителя скорости.

Результаты расчёта и эксперимента сравнивались при числе Рейнольдса Re=500, соответствующему ламинарному течению. Скорость потока на входе в трубу - 10 м/с. Установившееся течение в трубке определяется балансом вязких и невязких сил и характеризуется наличием пары вихрей, закрученных в противоположные стороны, образующихся вниз по потоку после прохождения трубочного изгиба. Имеются экспериментально измеренные профили скорости в пяти сечениях вдоль трубки: сечение 1 - 0.58 диаметра трубки вверх по потоку от изгиба, сечение 2 - 30 по изгибу, сечение 3 - 60 по изгибу, сечение 4 - 75 по изгибу, сечение 5 - один диаметр вниз по потоку от изгиба. Для численных исследований построено три неструктурированных гексаэдральных сетки. Грубая сетка №1 содержит 2 500 ячеек. При использовании РМГ K=3 на такой сетке задействованное количество степеней свободы равно 50 000. С целью сравнений результатов РМГ и МКО построена также сетка №2 с количеством ячеек 62 000 и сетка №3, содержащая 140 000 ячеек. Некоторые фрагменты сеток представлены на рис. 10.

Рис. 10. Фрагменты сеток Видно, что полигональная граница грубой сетки №1 существенно отклоняется от окружности. В расчётах вычисление интегралов граничных элементов выполнялось с учётом реальной границы, изображённой пунктиром.

На рис. 11 профили скорости вдоль сечений трубки, полученные с использованием схемы РМГ K=3 на сетке №1, сравниваются с аналогичными результатами, полученными МКО на сетке №3. Отметим, что результаты, полученные МКО на сетке №2 с числом элементов, обеспечивающим равное количество степеней свободы с РМГ, весьма далеки от экспериментальных данных. На рис. 11 экспериментальные результаты представлены кружками. Линия с квадратными маркерами демонстрирует результаты РМГ(K=3), в то время как пунктирная линия с треугольными маркерами показывает результаты МКО на сетке №3.

Рис. 11. Профиль скорости в 4-х сечениях трубки. Сопоставление эксперимента и численных результатов, полученных с РМГ K=3 на сетке №1 и МКО на сетке №Густота маркеров отражает густоту сетки, так как на одну ячейку приходится один маркер. Представленное сравнение демонстрирует эквивалентную точность расчётов по МКО и РМГ, а также их хорошее согласование с экспериментальными данными. Рассчитанный профиль скорости в сечении 1, как правило, во всех случаях находится в удовлетворительном согласовании с экспериментальными данными и поэтому здесь не представлен.

Таблица 4 представляет сравнение вычислительных затрат.

Использование схемы РМГ 4-го порядка точности на сетке №1 требует почти в 3 раза большего времени по сравнению со схемой МКО на сетке №2, обеспечивающей использование эквивалентного количества степеней свободы. Однако, принимая во внимание бльшую точность результатов РМГ схемы, целесообразно сравнить общие времена расчётных схем, требуемые для получения результатов с одинаковой точностью, которая, как видно из рис.11, обеспечивается МКО на сетке №3.

Отметим, что в расчёте на 50 итераций многосеточного алгоритма обе схемы затрачивают примерно одинаковое время. Расчёт новым методом на грубой сетке позволяет использовать почти в два раза больший шаг по времени, а, следовательно, до сходимости одинакового уровня требуется в два раза меньшее количество итераций. Таким образом, использование РМГ схемы для решения этой задачи позволило сократить не только объём оперативной памяти, но и общее расчётное время.

Таблица 4. Сравнение вычислительных затрат различных схем аппроксимации Сет- Кол-во Время на Общее ка Кол-во степ. Отнош. 50 MG кол-во Отнош № ячеек своб. памяти итераций итераций врем.

(секунды) РМГ 1 2502 50040 0.36 1185 4375 0.K=МКО 2 62689 62689 0.45 496 7880 0.МКО 3 139902 139902 1.00 1115 9155 1.В з4.7 оцениваются возможности схемы РМГ при решении задачи о распространении сферической акустической волны. Известно, что решение задач аэроакустики требует применения численных схем высокого порядка точности, обладающих минимальными дисперсионными и диффузионными свойствами. Исследуемая схема, как оказалось, обладает такими свойствами и может быть использована в практических приложениях. Свойства схемы были изучены на примере классической задачи о распространении акустического импульса, имеющего сравнительно малое начальное возмущение. Задача рассматривается в кубической расчетной области. Решение сеточных уравнений осуществляется методом РунгеЦКутта 5-го порядка точности с глобальным шагом по времени. Представленные результаты соответствуют времени, когда передний фронт возмущения еще не дошел до границ расчетной области. На рис. 12 приводится сравнение результатов для РМГ K=1, 2, 3 схем и традиционной схемы МКО 2-го попядка точности (NUMECA), полученных на последовательности равномерных структурированных сеток. Количество ячеек в сетке определяется как произведение числа разбиений расчётной области по каждому из направлений NxNyNz. На представленных рисунках в подписи к соответствующей кривой отражается не только размерность использованной сетки, но и количество степеней свободы, задействованных при выполнении соответствующего расчёта. В 4-х символьном обозначении: NxNyNzm число m указывает количество использованных переменных в ячейке. Видно, что расчёты, представленные на рис. 12, выполнялись на индивидуальных сетках, подобранных таким образом, чтобы задействованное количество степеней свободы было эквивалентным. Линии снабжены маркерами, из расчёта один маркер на ячейку. Поэтому плотность маркеров также отражает густоту сеток.

Рис. 12. Сравнение точного и численного решений для акустического импульса.

Эквивалентное число степеней свободы Сравнение компьютерных ресурсов для этого случая приведено в таблице 5. Схема 4-го порядка точности обеспечила хороший результат в 5.6 раза быстрее, использовав при этом в 1.4 раза меньшую память.

Отметим, что получить такую же точность результатов со схемой РМГ K=1 и, тем более, с использованием традиционной конечно-объемной схемы не представилось возможным изЦза чрезмерно большой потребной компьютерной памяти.

Таблица 5. Сравнение компьютерных ресурсов в случае близкой точности результатов для схем РМГ K=2 и K= Кол-во Исп. Отн. Время Отн.

Сетка степеней память памяти CPU врем.

свободы (Kb) РМГ 60x60x60 2 160 000 743 1.4 13ч 48мин 5.K=РМГ 32x32x32 655 360 528 1.0 2ч 28мин 1.K=В з4.8 выполнено сравнение схем МКО второго порядка точности и РМГ (K=1) в задаче о турбулентном обтекании изолированного крыла дозвуковым потоком вязкого газа (M=0.4, =0.6, Re=7.17106). В качестве тестового случая был выбран известный эксперимент Mller at al. (1996) с крылом LANN. Расчёты выполнены на неструктурированных гексаэдральных сетках, обеспечивающих использование эквивалентного количества переменных. Грубая сетка с числом элементов 190 2использована для РМГ, в то время как расчёты МКО выполнялись на сетке с числом элементов 625 076. Сравнения результатов расчёта распределённых характеристик обтекания демонстрируют, что оба подхода обеспечивают примерно одинаковые результаты, несмотря на использование заметно более грубой сетки в РМГ.

Сравнение схем МКО и РМГ при решении задачи об обтекании конфигурации крыло + фюзеляж представлено в з4.9. Самолётная конфигурация крыло+фюзеляж DLR-F4, выбранная в качестве тестовой модели конференцией AIAA CFD Drag Prediction Workshop (DPW), была рассчитана РМГ(K=1) на сетке, содержащей 230 тысяч узлов. Объём использованной сетки был ограничен возможностями персонального компьютера с 2-мя гигабайтами оперативной памяти и 32- битной операционной системой. Однако новый метод позволил получить результаты, удовлетворительно согласующиеся с экспериментальными данными и с расчётами, выполненными по известным коммерческим кодам МКО на более густых расчётных сетках. Отметим, что по классификации DPW сетка с объёмом ячеек 1.5 2 млн. считается грубой, но даже для расчёта на таких сетках требуется вычислительная система с большей оперативной памятью и 64 битной операционной системой, позволяющей работать с памятью, превышающей ~ 3 Гб.

Расчётная сетка №1 для конфигурации DLR-F4, сгенерированная кодом HEXPRESSЩ (пакета NUMECA), содержит 229 739 ячеек и 261 752 узлов. Фрагмент поверхностной сетки и общий вид изображены на рис. 13.

в) б) Рис. 13. Фрагмент поверхностной сетки: а) конфигурация крыло+фюзеляж включая плоскость симметрии, б) внешняя граница расчётной области Средняя аэродинамическая хорда крыла рассчитываемой конфигурации составляет 0.1412 м., длина фюзеляжа - 1.92 м., а полуразмах - 0.5856 м. Внешняя граница расчётной области представляет собой полуцилиндр радиусом 5 м. и высотой 20 м. Ячейки, прилегающие к обтекаемой границе, имеют дополнительное разбиение, состоящее из 4х слоёв. Расстояние до первого слоя узлов составляет 210-5 м., а толщина последующих слоёв увеличивается с коэффициентом геометрической прогрессии 1.4.

Расчёты выполнены при числах М=0.75 и Re=3106. Начальный коэффициент турбулентной вязкости на бесконечности в модели турбулентности Спаларта-Алмараса был равен 5. Расчёты по МКО были выполнены кодом Fine/HexaЩ (пакета NUMECA) на сетке №2, содержащей более 1 млн. ячеек. Количество слоёв в приграничных ячейках использованной сетки было равным 12, а коэффициент геометрической прогрессии 1.22. На рис. 14 и в таблице 6 представлены сравнения результатов расчётов выполненных РМГ на сетке №1 и МКО на сетке №2. Видно, что результаты по схеме РМГ находятся в лучшем согласовании с экспериментом DLR.

Сопоставление расчётных (РМГ K=1) и экспериментальных интегральных характеристик обтекания по углам атаки представлено на рис. 15. Сплошная красная линия отражает результаты расчёта, а маркерами изображены экспериментальные значения (Redeker, 1994), полученные в различных АДТ (DLR, ONERA, DRA). Пунктирные линии ограничивают область результатов получаемых различными современными МКО кодами, представленными в материалах конференции DPW (Lee-Rausch et al. 2003).

Отклонение расчётной поляры от экспериментальной составляет величину 210-3 (20 каунтов). Такое различие может быть объяснено несовершенством используемой модели турбулентности и недостаточным количеством задействованных степеней свободы задачи.

Тем не менее, это отклонение минимально среди других результатов, полученных МКО на сетках с большим количеством ячеек и с использованием больших компьютерных ресурсов.

Таблица 6. Сравнение интегральных характеристик расчёта при = Cy Cx МКО 0.5728 0.03РМГ 0.4773 0.03Эксперимент DLR 0.4812 0.02Рис. 14. Сравнение эпюр распределений давления в сечениях крыла Анализ расчётов данной конфигурации при нулевом угле атаки, выполненный различными вычислителями, демонстрирует наличие вихря в области сочленения задней кромки крыла с фюзеляжем (рис.16).

Расчеты, проведённые автором, подтверждают этот эффект.

Рис. 15. Интегральные характеристики обтекания конфигурации DLR-FРис. 16. Течение в области стыка крыла с фюзеляжем ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ И ВЫВОДЫ 1. На базе РМГ разработана новая гибридная схема, позволяющая непрерывно варьировать порядок точности в зависимости от локальных особенностей течения. В совокупности с предложенным сенсором схема высокого порядка точности обеспечила возможность получения монотонных решений в задачах с локальными неоднородностями, такими как скачок уплотнения или пристеночная часть пограничного слоя при использовании достаточно редких сеток.

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

3. Предложено два основных пути сокращения арифметических операций РМГ в расчёте на одну степень свободы:

o использование аналитического способа интегрирования потоков на низких многосеточных уровнях решателя;

o комплексная оптимизация квадратурных точек, используемых для расчёта поверхностных и объёмных интегралов.

4. Впервые метод РМГ высокого порядка точности реализован для решения пространственных задач аэродинамики и акустики с использованием неструктурированных гексаэдральных сеток и приведено его широкое сравнение с МКО.

5. Проведены численные исследования точности новой схемы аппроксимации на последовательности вложенных сеток. В процессе решения уравнений Эйлера и НавьеЦСтокса на неструктурированных сетках продемонстрирована оптимальная сходимость результатов к точным значениям с порядком K+1, где K - максимальный порядок базисного полинома.

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

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

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

Продемонстрированное в диссертации успешное применение метода РМГ высокого порядка точности к решению различных прикладных задач даёт основание рассчитывать на то, что в ближайшем будущем он составит достойную конкуренцию существующему подходу, основанному на МКО второго порядка.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ ОПУБЛИКОВАНЫ В РАБОТАХ (в научных журналах, представленных в перечне ВАК) 1. А.В. Волков. Особенности применения метода Галеркина к решению пространственных уравнений Навье-Стокса на неструктурированных гексаэдральных сетках // Ученые записки ЦАГИ. т. XL №6, 2009.

2. А.В. Волков, С.В. Ляпунов. Монотонизация метода конечного элемента в задачах газовой динамики // Ученые записки ЦАГИ.

т. XL №4, 2009.

3. В.Д. Боксер, А.В. Волков, А.В. Петров. Применение тангенциального выдува струй для снижения сопротивления сверхкритических профилей при больших дозвуковых скоростях // Ученые записки ЦАГИ. т. XL №1, 2009.

4. N.B. Petrovskaya, A.V. Wolkov, S.V. Lyapunov. Modification of Basis Functions in High Order Discontinuous Galerkin Schemes for Advection Equation // Applied Mathematical Modelling. Volume 32, Issue pp.826Ц835, May 2008.

5. В.Д. Боксер, А.В. Волков, А.В. Петров, Г.Г. Судаков. Можно ли улучшить аэродинамику сверхкритического крыла? Общероссийский научно-технический журнал Полет, стр. 70-76, 2008.

6. N.B. Petrovskaya, A.V. Wolkov. The Issues of Solution Approximation in Higher Order Schemes on Distorted Grids // International Journal of Computational Methods. (IJCM) v. 4, No. 2 Page: 367 - 382, June 2007.

7. А.В. Волков, С.В. Ляпунов. Применение конечно-элементного метода Галеркина с разрывными базисными функциями к решению уравнений Рейнольдса на неструктурированных адаптивных сетках // Ученые записки ЦАГИ. т. XXXVIII, № 3-4, стр. 22-30, 2007.

8. А.В. Волков, С.В. Ляпунов. Исследование эффективности использования численных схем высокого порядка точности для решения уравнений Навье-Стокса и Рейнольдса на неструктурированных адаптивных сетках // ЖВМиМФ, Том 46, №10, стр. 1894-1907, 2006.

9. А.В.Волков, С.В.Ляпунов, А.Н.Храбров. Влияние установившегося вращения на аэродинамические характеристики профиля при наличии отрыва потока // Ученые записки ЦАГИ. Том XXXIV, № - 4, 2003.

10. А.В.Волков, С.В.Ляпунов. Метод расчета вязкого отрывного обтекания систем крыловых профилей // Ученые записки ЦАГИ, том XXIX №3-4, 1998.

11. A.V.Wolkov, S.V.Lyapunov. Numerical Prediction of Transonic Viscous Separated Flow Past an Airfoil // Theoretical and Computational Fluid Dynamics, vol.6, №1, 1994.

12. А.В.Волков, С.В.Ляпунов. Метод расчета трансзвукового вязкого обтекания профиля с учетом изменения энтропии на скачках уплотнения // Ученые записки ЦАГИ, т.XXIV, №1, 1993.

ПРЕДСТАВЛЕНЫ В ТРУДАХ КОНФЕРЕНЦИЙ:

13. Ch.Hirsch, A.Wolkov, B.Leonard. Discontinuous Galerkin Method on Unstructured Hexahedral Grids // AIAA paper 2009-01771.

14. А.В.Волков. Анализ компьютерных ресурсов необходимых при численном решении 3-D уравнений Навье-Стокса методом DG.

Материалы ХХ-ой школы-семинара УАэродинамика летательных аппаратовФ, ЦАГИ, 2009 г.

15. Charles Hirsch, Andrey Wolkov, Benoit Leonard. Discontinuous Galerkin Method on Unstructured Hexahedral Grids for 3D Euler and NavierStokes Equations. 5th. European Congress on Computational Methods in Applied Sciences and Engineering //ECCOMAS 200, held at the Lido Island in Venice (Italy) on 30 June - 4 July, 2008.

16. А.В. Волков. Учёт кривизны обтекаемой поверхности в методе высокого порядка точности. Материалы ХIХ-ой школы-семинара УАэродинамика летательных аппаратовФ, ЦАГИ, 2008 г.

17. A.Wolkov, Ch.Hirsch, B.Leonard. Discontinuous Galerkin Method on Unstructured Hexahedral Grids for 3D Euler and Navier-Stokes Equations // AIAA paper 2007-4078.

18. А.В. Волков. Применение p-многосеточного метода к ускорению решения конечно-элементного метода аппроксимации при решении пространственных уравнений Навье-Стокса. Материалы ХVIII-ой школы-семинара УАэродинамика летательных аппаратовФ, ЦАГИ, 2007 г.

19. V.Vlasenko, A.Wolkov, Ch. Hirsch, Computationally effective Discontinuous Galerkin scheme for Linearized Euler Equations // WEHSFFC 2007. Moscow, November 19-22, 2007.

20. А.В. Волков. Применение метода Галеркина с разрывными базисными функциями к решению пространственных уравнений Эйлера и Навье-Стокса. Труды VI международной школы семинар ММА Модели и Методы Аэродинамики, Евпатория 2006.

21. А.В.Волков, С.В.Ляпунов. Численный метод высокого порядка точности для решения системы уравнений Навье-Стокса и Рейнольдса. ММА (Модели и Методы Аэродинамики). Материалы Третьей Международной школы-семинара, Евпатория 2003.

22. S.V.Lyapunov, A.V.Wolkov. Application of Discontinuous Galerkin finite element method to the solution of partial differential equations. Part II.

System of nonlinear equations: Euler equations // 16th IMACS World Congress. Lausanne-Switzerland August 21-25, 2000.

23. A.V.Wolkov, S.V.Lyapunov. Application of Discontinuous Galerkin finite element method to the solution of partial differential equations. Part I. 2D scalar conservation laws // 16th IMACS World Congress. LausanneSwitzerland August 21-25, 2000.

24. V.S.Sakovich, A.M.Sorokin, A.V.Wolkov, S.V.Lyapunov. Anisotropic unstructured grid generation for 3D flow simulation problems. 6th Int.

Conf. On numerical grid generation in computational fluid simulation, 1998.

25. A.V.Wolkov, S.V.Lyapunov. A Separated Flow Calculation About Airfoils and High-Lift Systems on Basis of Viscous-Inviscid Interaction Methods. EUROMECH, 3rd European Fluid Mechanics Conference, 1997, Gottingen, Germany 26. S.V.Lyapunov, A.V.Wolkov. Application of Viscous-Inviscid Interaction Methods for a Separated Flow Calculation About Airfoils and High-Lift Systems. ICAS Proceedings, ICAS-96-1.10.2, 1996.

Авторефераты по всем темам  >>  Авторефераты по техническим специальностям