Численное моделирование обтекания моделей пассажирских самолетов в условиях ограниченного пространства и влияния элементов конструкции аэродинамической трубы тема диссертации и автореферата по ВАК РФ 05.07.01, кандидат технических наук Курсаков, Иннокентий Александрович

• Курсаков, Иннокентий Александрович
• кандидат технических наук
• 2011

Оглавление диссертации кандидат технических наук Курсаков, Иннокентий Александрович

Введение.

Глава 1. Описание численного метода. Тестовые расчеты.

1.1 Базовые уравнения.

1.1.1 Общая формулировка законов сохранения и численной схемы.

1.1.2. Система уравнений Рейнольдса и ее замыкание.

1.2 Аппроксимация системы уравнений Рейнольдса.

1.2.1 Конвективные члены схемы.

1.2.2 Диффузионные потоки схемы.

1.2.3 Источниковые члены схемы.

1.3 Реализация многосеточного подхода.

1.4 Особенности организации параллельных вычислений.

1.5 Валидация численного метода.

1.5.1 Тело вращения на хвостовой державке.

1.5.2 Обтекание компоновки крыло-фюзеляж-мотогондола.

1.5.3 Сверхзвуковое обтекание оживального тела под углом атаки.

Глава 2 Особенности моделирования нестандартных режимов работы силовой установки самолета в аэродинамической трубе с открытой рабочей частью

2.1 Особенности испытаний «реверсной» модели в АДТ Т-104 ЦАГИ.

2.2 Математическая модель аэродинамической трубы с открытой рабочей частью и стендом для исследования реверса тяги двигателя самолета.

2.3 Моделирование течения в рабочей части АДТ Т-104.

2.4 Влияние пилона на эффект попадания реверсных струй в двигатель самолета.

2.5 Влияние «бегущей дорожки» на эффект попадания реверсных струй в двигатель самолета.

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

3.1 Экспериментальная методика определения влияния поддерживающих устройств.

3.2 Особенности построения математической модели для задачи учета влияния поддерживающих устройств.'.

3.2.1 Выбор оптимальной топологии расчетной сетки.

3.2.2 Физические особенности течения в узле стыковки модели и хвостовой державки.

3.3 Определение влияния килевой державки.

3.4 Выбор формы килевой державки.

Выводы.

Введение диссертации (часть автореферата) на тему «Численное моделирование обтекания моделей пассажирских самолетов в условиях ограниченного пространства и влияния элементов конструкции аэродинамической трубы»

В настоящее время во всех областях деятельности человека активно применяются быстродействующие вычислительные машины. Как отмечает академик A.A. Самарский в одной из своих работ: «Широкое применение математических методов позволяет поднять общий уровень теоретических исследований, дает возможность проводить их в более тесной связи с экспериментальными исследованиями» [1] С 1970-х годов роль вычислительных технологий постоянно возрастала, что обусловлено как совершенствованием численных методов, так и уменьшением стоимости вычислительных мощностей. Это хорошо заметно в ЦАГИ, где взаимосвязь теоретических, экспериментальных и расчетных исследований положена в «фундамент института». Наиболее полный анализ успехов того времени проведен в монографиях, подготовленных под редакцией академика РАН Г.С. Бюшгенса [2], член-корреспондента РАН В.Я. Нейланда [3], член-корреспондента РАН В.Г. Дмитриева. Следует отметить, что именно в ЦАГИ в 1972 году разработан метод MUSCL [4], который определил «лицо» вычислительной аэродинамики в мире на много лет вперед. Это признал в своей работе [5] один из наиболее известных ученых США Van Leer. Другая фундаментальная работа сотрудника ЦАГИ А.Н. Минайлоса [6] остановила «бесполезный» в то время поиск эффективных методов высокого порядка и позволила сконцентрироваться ученым на доступных схемах второго порядка аппроксимации, уделив особое внимание свойству монотонности. Более тридцати лет развитие вычислительной аэродинамики шло по этому пути. И только в последние годы композиция двух методов известных российских ученых (академиков РАН Б.Г. Галеркина и С.К. Годунова) позволила сделать следующий шаг. Найден путь к созданию эффективного метода высокого порядка точности - метода Галеркина с разрывными базисными функциями [7]. Особый вклад в разработку этого метода внесли ученые ЦАГИ C.B. Ляпунов [8] и A.B. Волков [9].

Естественно, что вычислительные технологии, прежде всего, востребованы в аэрокосмической отрасли. Это обусловлено чрезвычайной сложностью решаемых задач. Прежде всего, это многодисциплинарные задачи, объединяющие аэродинамику, прочность, динамику полета и другие отрасли знания. Подробный анализ современного состояния многодисциплинарных подходов проведен в работе Н.Г Бунькова [10]. Показано, что реальная потребность в вычислительных мощностях на два порядка превышает имеющиеся ресурсы. Другой класс решаемых задач произрастает из особой роли ЦАГИ, как крупнейшего экспериментального центра. В настоящее время появилась возможность осуществить численное моделирование течения в Аэродинамической Трубе (АДТ). При этом учитывается обтекание, как экспериментальной модели, так и элементов упомянутой АДТ. Можно указать, по крайней мере, два основных направления развития указанной технологии:

1. Исследование особенностей течения в АДТ и использование полученных результатов для совершенствования экспериментальной методологии, а, возможно, и для модернизации трубы, например [11, 12]

2. Учет влияния элементов конструкции АДТ, таких как перфорированные стенки, поддерживающие устройства различных типов, технологические полости и т.д. на результаты испытаний, проводимых в этой АДТ.

В ЦАГИ по инициативе и под научным руководством В .Я. Нейланда [13] разработана концепция «Электронной Аэродинамической Трубы (ЭАДТ)», которая объединяет оба указанных направления. Первые публикации в России [14] и за рубежом [15] имели положительный отклик. Особенно следует отметить, что ЭАДТ ни в коем случае не подразумевает отказа от экспериментальных исследований. Применение расчетных методов в качестве инструмента для восполнения экспериментальных данных оправдано с экономической и практической точек зрения [16]. Современные экспериментальные методики в сочетании с подходами ЭАДТ дополняют друг друга. Так, эксперимент позволяет получить результаты высокого уровня 8 точности, а расчеты восполнить их недостающими материалами. Например, ЭАДТ дает возможность осуществить корректный пересчет результатов эксперимента в АДТ на условия натурного полета. В настоящее время [17] можно указать четыре области успешного применения концепции ЭАДТ в технологическом цикле реального эксперимента. Это - 1) стадия подготовки модели; 2) стадия проведения испытаний; 3) стадия обработки экспериментальных данных; 4) стадия вторичной обработки экспериментальных данных (на этой стадии возможно восполнение экспериментальных данных в тех случаях, когда это необходимо). Так на этапе подготовки эксперимента обычно исследуется масштабный эффект и выбирается модель максимально большого размера для условий в данной АДТ. Кроме того, определяются критические зоны, в которых располагаются датчики повышенной точности, например, вакуумметры для измерения статического давления. В процессе проведения испытаний часто проводится адаптация проницаемости перфорации АДТ (A.B. Семенов в Т-128 [18]) с использованием предварительно рассчитанных полей течения в окрестности изолированной модели. При этом устраняется влияние стенок на результаты указанных испытаний (С.А. Глазков в Т-128 [19]). Наконец, проводится сопоставление расчетных результатов, полученных с учетом и без учета стенок АДТ, и вносятся необходимые коррекции в результаты эксперимента. Попутно, экспериментальные распределения давлений по поверхностям модели восполняются расчетными значениями давления в возмущенном поле. Очевидно, что реализация концепции ЭАДТ требует исключительно больших ресурсов ЭВМ. Кроме этого необходимо приложить усилия к созданию эффективной расчетной технологии. Эффективность технологии - понятие неоднозначное. Прежде всего, это удобство использования вычислительной техники и программного обеспечения. Но главным вопросом является скорость, другими словами, получение результатов в ограниченные сроки. Можно выделить основные группы факторов, ведущих к ускорению скорости счета: 9

Выбор оптимального алгоритма расчета;

Оптимизация кода;

Организация параллельных вычислений.

В данной работе основное внимание уделено первому пункту как наиболее трудоемкому и перспективному. Автор использует много сеточный подход, как наиболее подходящий для идеологии ЭАДТ. Кроме того, важно отметить, что существенный ресурс заложен в создании параллельного процесса. Параллельное программирование представляет собой отдельную задачу высокой степени сложности. Для того чтобы решить задачу на параллельном компьютере, необходимо распределить вычисления между процессорами системы, так чтобы каждый процессор был занят решением своей части. Кроме того, желательно оптимизировать межпроцессорные коммуникации, которые существенно тормозят передачу данных от одного процессора к другому. Часто возникают парадоксальные ситуации, когда с увеличением числа процессоров скорость расчета не увеличивается, а замедляется. Из-за сложности параллельных вычислительных систем и их существенного отличия от традиционных однопроцессорных компьютеров нельзя напрямую воспользоваться традиционными языками программирования. С точки зрения разработчика математического обеспечения, можно выделить несколько моделей параллельных вычислений: процесс/канал (Process/Channel), обмен сообщениями (Message Passing), параллелизм данных (Data Parallel), общей памяти (Shared Memory).

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

10 синхронизации - семафоры и блокировки процессов. Такая модель параллельных вычислений носит название PRAM (Parallel Random Access Machine, или абстрактная модель вычислений с параллельным случайным доступом). На практике время доступа к памяти на реальных машинах неоднородно из-за возможной иерархической организацией модулей памяти. Частота обращений к этой памяти может быть уменьшена за счет сохранения копий используемых данных в кэш-памяти, связанной с каждым процессором. Доступ к кэш-памяти намного быстрее, чем непосредственный доступ к общей памяти.

Используемая в работе модель параллельных вычислений накладывает ограничения на выбор архитектуры параллельных компьютеров. Это мультипроцессорные системы (Multiprocessing Systems) с общей или разделяемой памятью (Shared Memory). Естественно, наиболее оптимальным выбором являются системы, поддерживающие идеальную PRAM модель. Все процессоры в таких системах являются абсолютно одинаковыми, кроме того, имеется общая адресация оперативной памяти. Все процессоры также имеют одинаковое время доступа к оперативной памяти. Поэтому такие системы еще называют SMP (Symmetric Multiprocessing System) или UMA (Uniform Memory Access). SMP системы являются наиболее простыми, с точки зрения создания параллельных программ путем разбиения задачи на нити (Threads). При этом система автоматически распределяет нити по процессорам. Пользователь заботится лишь об оптимальном распределении узлов расчетной сетки по блокам, приписываемым соответствующим нитям. В рамках объектного подхода эта проблема решается достаточно просто путем описания блоков в виде классов. Узкое место таких систем - проблема согласования кэш-памяти. Действительно, в кэш-памяти двух процессоров может быть одновременно считано и модифицировано содержимое одной и той же ячейки общей памяти. Однако такая архитектура является очень сложной и дорогостоящей, с технологической точки зрения, и фактически ограничивается двух- и

11 четырехпроцессорными системами. Сходные возможности организации параллельных вычислений предоставляют компьютеры с виртуальной общей (разделяемой) памятью (Virtual sliared memory). У таких компьютеров общая память как таковая отсутствует. Каждый процессор имеет собственную локальную память и может обращаться к локальной памяти других процессоров, используя "глобальный адрес". Если "глобальный адрес" указывает на область за пределами локальной памяти, то доступ к памяти реализуется с помощью сообщений, пересылаемых по коммуникационной сети.

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

За рубежом существует успешный опыт реализации указанной концепции [20]-[26]. Первое упоминание об ЭАДТ появилось в литературе в 1981 году в статье [20], где появился термин Electronics Wind Tunnel. К сожалению, это статья не содержала серьезной проработки деталей алгоритмов и рекомендаций по практической реализации всей концепции. Настоящий прорыв произошел в 1997 году на конференции STAI (Supersonic Tunnels Association International). В процессе доклада ЦАГИ [27] было задано множество уточняющих вопросов, а на следующий день была организована специальная сессия для обсуждения

12 практических аспектов реализации концепции Е\¥Т.На этой сессии выступили такие ученые, как Э. АевсЬИшап (США) и Б. 81апш1апс1 (Великобритания). Это дало старт к развитию технологии моделирования испытаний в АДТ во многих странах мира. В настоящее время уже разработаны несколько компьютерных программ, при этом, особенно выделяются работы [21],[22], в которых описаны работающие прототипы систем, объединяющих экспериментальный и вычислительный подходы в условиях конкретных АДТ. Кроме этого существуют исследования, продвигающие отдельные направления и посвященные разработке сопутствующих алгоритмов и методов [18],[19],[28]. Например, важную роль имеет адаптация закона фильтрации Дарси [29], который используется для моделирования перфорации стенок АДТ.

В настоящее время существуют примеры успешного применения концепции ЭАДТ для решения практических задач. В работе [23] приведены результаты исследования модели перспективного высокоскоростного ЛА с выпущенной механизацией крыла. Расчеты проведены как для изолированной модели, так и для модели с поддерживающим устройством с учетом и без учета стенок трубы. Выполнена коррекция результатов эксперимента. Получено хорошее соответствие расчетных и экспериментальных характеристик модели во всем диапазоне чисел Мах и углов полета. Отмечено, что расчетные исследования дают хорошее представление о том, какое влияние на поток оказывают поддерживающие устройства. Работа [24] — это другой пример удачного синтеза экспериментального и численного подходов. Основная цель представленной работы заключается в создании методологии проектирования хвостового оперения модели. Экспериментально исследован широкий спектр интерференционных явлений, таких как отрыв на хвосте, падение эффективности руля направления и т.д. В рамках численного исследования обтекания модели ИЕМР1 [25] в конфигурациях «свободный поток» и модель на ножевой державке были определены зоны наибольшего влияния Поддерживающих Устройств (ПУ), что позволило дополнить, а впоследствии и

13 скорректировать экспериментальные данные. В работе [26] выполнен вязкий расчет и проведено исследование интерференции ПУ и моделей JIA при испытаниях в АДТ. Поля возмущенного течения и поправки к интегральным характеристикам были получены при помощи сравнения расчетов с ПУ и без него. Рассмотрено два типа конфигураций ПУ В первой конфигурации моделируются низкоскоростные испытания модели транспортного самолета, установленной на цилиндрической опоре. Во второй — трансзвуковое обтекание модели административного самолета на Z-образной державке. Влияние державок исследуется для нескольких чисел М и углов атаки. Особое внимание уделено методике получения поправок к параметрам набегающего потока.

Для решения задачи моделирования потока в АДТ наряду с коммерческими программными продуктами широкого профиля (такими как ANSYS Fluent [30], ANSYS CFX [31], Numeca [32]), используются специализированные пакеты программ, разрабатываемые в крупных мировых центрах авиационной науки (in-house codes). Примерами таких продуктов являются elsA (Onera, Франция) [33], TAU [34], FLO Wer [35] (DLR, Германия). Использование собственных программных продуктов в этих центрах считается предпочтительным, поскольку в этом случае имеется полное представление об алгоритмах и способах получения числениого решения. Несмотря на существенные особенности указанных програмных продуктов, можно выделить общие подходы к построению численной методологии. В настоящее время вследствие ограниченности ресурсов решается математическая задача (задача Коши), сформулированная для осредненной по Рейнольдсу системы уравнений Навье-Стокса [36]. Для аппроксимации конвективных потоков используется метод Джеймсона [37] второго порядка точности с искусственной вязкостью, для диффузионных - центрально-разностная схема. Указанная система уравнений замкнута одной из известных дифференциальных моделей турбулентности: 1) #-со [38]; 2) SST [39]; 3) Спаларта-Алмараса [40],[41]. При

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

Ключевым вопросом при построении численной методологии является выбор способа интегрирования по времени. На сегодняшний день существует два принципиально разных подхода: явный и неявный [42]. Основным преимуществом неявных численных схем является их абсолютная устойчивость, это позволяет, в принципе, двигаться сколь угодно быстро по времени. Однако при моделировании нестационарных процессов большие шаги по времени должны приводить к большим погрешностям аппроксимации. Главным же недостатком явных численных схем является сильное ограничение на шаг по времени, что приводит к огромным временам счета при решении задач с отличающимися характерными линейными масштабами. Однако скорость получения решения может быть увеличена, например, за счет применения многошаговых процедур Рунге-Кутта [37], технологии «локального шага» по времени [43]. Значительное ускорение расчета обеспечивается применением многосеточных методов (МиШ§пс1). Впервые метод был описан в работе [44]. В задачах вычислительной аэродинамики многосеточный подход последовательно применялся в рамках уравнений полного потенциала [45],[46], Эйлера [37],[47-49] и в дальнейшем получил своё развитие применительно к системе уравнений Навье-Стокса [50],[51]. Этот подход базируется на свойстве волн (возникающих при несоответствии текущих полей течения заданным граничным условиям) распространяться со скоростью, зависящей от размеров ячеек. Длинные волны быстро продвигаются по крупной сетке, а мелкие ячейки отвечают за коротковолновую часть спектра. Построение серии вложенных сеток и продвижение возмущений по «крупной» сетке дает десятикратный выигрыш во времени при решении стационарной задачи. Особое внимание следует уделять «угловым» точкам. Желательно, чтобы все вложенные сетки выделяли эти точки своими узлами. При разработке

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

Численное моделирование обтекания модели самолета в условиях АДТ имеет ряд специфических особенностей. Для расчета течения газа в условиях АДТ разработаны специальные граничные условия [17]. Например, запуск «Электронной Аэродинамической Трубы» и выход на режим осуществляется итерационно. При этом контролируется величина скорости потока не на внешней границе расчетной области, а в специально обозначенной точке канала АДТ, называемой «Контрольной точкой». Именно в этой точке в процессе реального эксперимента измеряется скорость потока в АДТ. Аналогичным образом в этой точке задаются амплитуда и масштаб турбулентных пульсаций. При изменении угла атаки или скольжения модели невозможно применять простую процедуру, используемую для случая невозмущенного потока, а приходится перестраивать расчетную сетку, так как стенки АДТ при этом остаются неподвижными. Отметим, что на стенках нужно выполнять условие частичного протекания потока и обмена газа между рабочей частью АДТ и ее камерой давления. Отдельное внимание уделяется донным областям модели, в которых осуществляется стыковка с державкой. Как правило, здесь появляется нестационарное течение, которое требует реализации специальных расчетных процедур, таких как «Дробный шаг по времени» [52] или «Дуальный шаг по времени» [53]. Удобнее всего эти процедуры используются в рамках зонального подхода [54], который позволяет локализовать нестационарную часть задачи и существенно сократить время расчета.

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

16 структурированная сетка, адаптированная к особенностям геометрии. Применение структурированной сетки позволяет обеспечить максимальную точность результатов. Неструктурированная сетка обладает несомненным преимуществом удобства построения, но для достижения сопоставимой точности результатов расчета она требует существенно больше узлов. Следует отметить, что в настоящее время в ЦАГИ ведутся работы по развитию алгоритмов построения неструктурированных сеток на основе многогранников (polyhedral) с призматическим пограничным слоем. По утверждению авторов новой технологии [55], объем массивов памяти, необходимых для реализации указанных сеток, существенно меньше, чем в упомянутом выше случае. При этом время расчета приближается к соответствующему времени с использованием структурированных сеток (при сопоставимой точности). Для упрощения процесса построения структурированной сетки широко используется технология «Химера» [56]-[58], которая позволяет создавать независимые сетки для каждой части летательного аппарата и элементов экспериментальной установки.

Наиболее ответственный этап технологии - выбор расчетного метода. Существуют различные направления и школы. Как уже упоминалось ранее, в ЦАГИ ведутся работы [59] по развитию метода Галеркина применительно к неструктурированным адаптивным сеткам. В данной работе применяется классический подход, основанный на модификации метода Годунова-Колгана [4],[60],[61]. Он наиболее полно удовлетворяет требованиям, предъявляемым к решению сложных задач аэродинамики, включая задачи интерференции газовых струй с элементами планера и частями АДТ. Одна из задач посвящена исследованию распространения реверсных струй при торможении самолета на взлетно-посадочной полосе. Реверсирование тяги двигателя является одним из эффективных способов торможения самолета при послепосадочном пробеге. Однако при включении реверсивного устройства может возникнуть ряд проблем, обусловленных действием реверсивных потоков. Попадание

17 реверсных струй во входные устройства двигателей становится причиной искажения полей скоростей и температур на входе в двигатели, что является предпосылкой к возникновению помпажного режима работы компрессора и выключению двигателя [62]. Моделированию течения в реверсных устройствах посвящены работы [63-65]. При проведении экспериментальных исследований работы реверса тяги необходимо вносить ряд существенных ограничений. Например, при моделировании работы двигателя необходим подвод сжатого воздуха, который выдувается через сопло. Такой подвод осуществляется посредством специального пилона, который служит дополнительным ПУ и располагается в непосредственной близости от объекта исследования. Кроме того, моделирование в эксперименте движения самолета по взлётно-посадочной полосе чрезвычайно сложно с технической точки зрения.

Второй пример применения разрабатываемой методологии связан с вопросом корректного учета влияния державок на результаты эксперимента. Хвостовая державка рассмотрена более подробно. Это связано с ее большим распространением и интересными эффектами, возникающими в местах интерференции этой державки с фюзеляжем. Рассмотрена упрощенная модель цилиндрического тела с косым срезом, установленного на хвостовой державке. Особое внимание уделено структуре течения в зоне сочленения модели с державкой. В этой зоне из-за конструктивных особенностей крепления модели образуется полость, в которую попадает поток с внешней поверхности модели. Экспериментальное исследование течения в данной области затруднено [66]. Расчеты обтекания указанной конфигурации проведены в некоторых работах, например, [67,68]. При этом течение в полости узла стыковки ПУ и модели не учитывалось. Однако, как отмечено в [66],[68] наличие полости приводит к искажению, как интегральных характеристик, так и локальных полей течения в окрестности модели. В данной работе указанный недостаток устранен, расчеты выполнены с учетом каверны.

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

Заключение диссертации по теме «Аэродинамика и процессы теплообмена летательных аппаратов», Курсаков, Иннокентий Александрович

Выводы

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

1. Предложена модификация вычислительной программы EWT-LIAГИ путем внедрения в алгоритм расчета многосеточного алгоритма. Показано, что такая модификация позволяет ускорить время получения результата от двух до пяти раз.

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

3. Разработанная методология применена для моделирования работы экспериментальной установки «Реверсный стенд» в аэродинамической трубе Т-104 ЦАРИ. Полученный опыт позволяет заключить, что:

3.1.Наличие «бегущей дорожки», являющейся аналогом взлётно-посадочной полосы, приводит к сносу вихря, возникающего при взаимодействии реверсных струй с пограничным слоем, вниз по потоку;

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

3.3.Такие элементы экспериментальной установки как пилон для подвода воздуха к имитатору двигателя и неподвижный имитатор взлетно-посадочной полосы дают поправки разного знака к величине скорости «отсечки» реверсных струй, что позволяет применять результаты экспериментальных исследований с использованием указанного стенда без «бегущей дорожки» на практике.

4. Разработанная методология применена к расчету обтекания модели перспективного пассажирского самолета, установленной на килевой державке. Показано, что применение указанной методологии при обработке экспериментальных данных по влиянию килевой державки в АДТ Т-128 позволяет увеличить точность вычисления поправок к аэродинамическим коэффициентам до величин д \ АСха \ < 0,0001 ; 81А Су а \ < 0,001 на безотрывных режимах обтекания крыла модели по сравнению с методикой, основанной на линейных подходах.

5. Путём численного моделирования потока в области стыковки хвостовой державки с моделью показано, что течение в полости крепления державки создает нагрузки на модель, которые воспринимаются тензометрическими весами, но не имеют отношения к аэродинамическим нагрузкам при свободном обтекании. Дополнительные поправки к показаниям весов, обусловленные указанным влиянием, в рассматриваемом случае оцениваются величинами порядка АСха= 0,0007, АСуа = 0,01.

6. Разработанная методология внедрена в технологический цикл АДТ Т-104 и Т-128 ЦАГИ.

Список литературы диссертационного исследования кандидат технических наук Курсаков, Иннокентий Александрович, 2011 год

1. Самарский А.А., Вабищевич П.Н. Математическое моделирование и вычислительный эксперимент, «Математическое моделирование», 2000 г

2. Бюшгенс Г.С., Берджицкий E.JT. «ЦАГИ центр авиационной науки» // Москва, Наука 1993.

3. Под редакцией Нейланда В.Я. ЦАГИ основные этапы научной деятельности 1968-1993. М. Физматлит, 1996.

4. Колган В.П. Применение принципа минимальных значений производной к построению конечно-разностных схем для расчета разрывных решений газовой динамики // Ученые записки ЦАГИ. 1972. Т. 3, № 6, с. 68-77.

5. В. van Leer. Upwind and High-Resolution Methods for Compressible Flow: from Donor Cell to Residual Distrubution Schemes. Review Article. Commun. Comput. Phys., 1(2): 192-206, 2006.

6. Минайлос A.H. О значении монотонности конечно-разностных схем в методах сквозного счета. «ЖВМ и МФ», т. 17, №4, 1977, стр. 1058-1063.

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

8. Волков А.В. Особенности применения метода Галеркина к решению пространственных уравнений Навье-Стокса на неструктурированных гексаэдральных сетках // Ученые записки ЦАГИ Том XLI, №3, 2010

9. Ю.Дмитриев В.Г., Белоцерковский С.М., Буньков Н.Г. Системная роль математической компьютерной модели самолета в его жизненном цикле. «Техника воздушного флота», №4-5, 1998, стр. 7-17.127

10. Zhigang Yang, Qing Jia Assessment of Wind Tunnel Test Section Dimensions Using CFD I I 46th AIAA Aerospace Sciences Meeting and Exhibit 7-10 January 2008, Reno, Nevada

11. Joao B.P. Falcao Filho , Marcos A. Ortega Numerical study of the injection process in a transonic wind tunnel: The numerical details // Computers & Fluids 37 (2008) 1276-1308

12. Neyland V., Bosniakov S., Glazkov S., IvanovA., Matyah S., Mikhailov S., Vlasenko V. Conception of electronic wind tunnel and first results of its implementation // Progress in Aerospace Sciences. 2001. V. 37, Issue 2, p. 121-145.

13. Daniel Reckzeh and Heinz Hansen. High Reynolds-number windtunnel testing for the design of airbus high-lift wings. In New Results in Numerical and Experimental Fluid Mechanics V, pages 1-8. Springer, 2006.

14. Практические аспекты решения задач внешней аэродинамики двигателей летательных аппаратов в рамках осредненных по времени уравнений Навье-Стокса // Сборник статей, Труды ЦАГИ. 2007, вып. 2671

15. Neyland V.M., IvanovA.I., SemenovA.V., Semenova O.K., Amirjanz G.A. Adaptive-wall perforated test section for transonic wind tunnels // AGARD CP-585. June 1997. P. 16.1-16.16.

16. Glazkov S.A., Gorbushin A.R., IvanovA.I., SemenovA.V. Recent experience in improving the accuracy of wall interference corrections in TsAGI T-128 wind tunnel //Progress in Aerospace Sciences 37 (2001) 263-298

17. Baals D.D., Corliss W.R. "Wind Tunnels of NASA" Scientific and Technical Information Branch NASA, 1981.

18. Schwartz R.J., Fleming G.A. "Virtual Diagnostics Interface: Real Time Comparison of Experimental Data and CFD Predictions for a NASA Ares I-Like Vehicle," // Proceedings of ICIASF 07, R56, 2007

19. Shigeya Watanabe, Shigeru Kuchi-ishi, Takashi Aoyama "A Prototype System towards EFD/CFD Integration: Digital/Analog-Hybrid Wind Tunnel" // ICAS 2010

20. Wendy B. Lessard. Analysis of post-support and wind-tunnel wall interference on flow field about subsonic high-lift high-speed research configuration. Technical Report TP-2000-210555, NASA, November 2000.

21. Karl Pettersson and Arthur Rizzi. Estimating reynolds number scaling and windtunnel boom effects with the help of CFD methods. AIAA 2006-3162, 2006.

22. European Comission, editor. Aeronautics Research 2003-2006 projects, volume Project synopses volume 1 Research Projects from the first and second calls. Office for Official Publications of the European Communities, 2006

23. S. Mouton, Numerical investigations of model support interference in subsonic and transonic wind tunnels // 44ème Colloque d'aérodynamique appliquée AAAF,Nantes, March 23-25, 2009

24. Bosniakov S., Fonov S., Matyash S., Mikhailov S., Neyland V., Remeev N., Vlasenko V., Yatskevich N. CFD as a part of experimental technology in a Wind Tunnel // Proceedings of STAI meeting, Nashville, USA, 1997.

25. Collercandy R., Marques B, Lory J., Dbjay S., Espiau L. Application of CFD for wall and sting effects. HiReTT report HIRETTTNAFRCoWP2.2311002003, Airbus France, October 2003

26. Anderson D., TannehilJ., Fletcher R. Computational Fluid Mechanics and Heat transfer book. Washington DC: Hemisphere, 1984.

27. Jameson A,, Schmidt W., Turkel E. Numerical Solution of the Euler Equations by Finite Volume Methods using Runge-Kutta Time-Stepping Schemes // AIAA-81-1259 (1981).

28. Wilcox D.C, Turbulence modeling for CFD. 2nd edition. DCW Industries, 1998.

29. Menter F.R. Improved two-equation (k--^) turbulence models for aerodynamic flows // NASA TM-103975. 1992.

30. Spalart P.R., Allmaras S.R. A one-equation turbulence model for aerodynamic flows //AIAA Paper 92-439. Reno, NV. 1992.

31. Edwards J.R., Chandra S. Comparison of Eddy Viscosity-Transport Turbulence Models for Three-Dimensional, Shock-Separated Flowfields. // Journal of Aircraft, Vol. 34, No. 4, pp. 756-763, 1996.

32. BlazekJ. Computational Fluid Dynamics: Principles and Applications. Elsevier, 2001.

33. Jameson A. A perspective on computational algorithms for aerodynamic analysis and design // Progress in Aerospace Sciences. 2001. V. 37 N2

34. Федоренко Р.П. «Многосеточный метод для схем конечного элемента» // ЖВМ и МФ 1961 - Т1 N5

35. Jameson A. Acceleration of Transonic Potential Flow Calculations on Arbitrary Meshes by the MultiGrid Method // AIAA-79-1458-CP (1979).

36. Jameson A., CaugheyD. A Finite Volume Method for Transonic Potential Flow Calculations // pp. 35-54 in Proc. AIAA 3rd CFD Conf., Albuquerque, NM (1977).

37. Jameson A.,Yoon S. Multigrid Solution of the Euler Equations Using Implicit Schemes //AIAA-85-0293 (1985).

38. Mavriplis D.J. Three-Dimensional Unstructured Multigrid for the Euler Equations // AIAA-91-1549-CP (1991).

39. Brandt A. Multilevel adaptive computations in fluid dynamics// AIAA Journal, V.18,N 10, October 1980

40. Jameson. A., Martinelli L. Validation of a Multigrid Method for the Reynolds Averaged Equations // AIAA-88-0414 (1988).

41. Mavriplis D., Aftosmis M. Berger M. High Resolution Aerospace Applications Using the NASA Columbia Supercomputer // Proc. Supercomputing'05, Seattle, WA, November 12-18 (2005).

42. Pervaiz M., Baron J.R. Spatiotemporal adaptation algorithm for two dimensional reacting flows // AIAA Journal. 1989. V. 27, №10.

43. Jameson A. Time dependent calculations using multigrid, with application to unsteady flows past airfoils and wings // AIAA 91-1596. 1991.

44. КажанЕ.В. О возможностях использования неявной схемы в рамках пакета EWT-ЦАГИ // Труды ЦАГИ. 2007, вып. 2671.

45. Peric М. Flow simulation using control volumes of arbitrary polyhedral shape // ERCOFTAC Bulletin. September 2004. No. 62

46. Benek, J.A., P.G. Buning and J.L. Steger. A 3-D Chimera Grid Embedding Technique; AIAA-85-1523 (1985)

47. Meakin, R.L. and N. Suhs. Unsteady Aerodynamic Simulations of Multiple Bodies in Relative Motion; AIAA-89-1996 (1989).

48. Dougherty, F.C. and J. Kuan. Transonic Store Separation Using a Three-Dimensional Chimera Grid Scheme; AIAA-89-0637 (1989).

49. Годунов C.K., Забродин A.B., Иванов М.Я., Крайко А.Н., Прокопов Г.П. Численное решение многомерных задач газовой динамики. М.: Наука, 1976.

50. Родионов А.В. Повышение порядка аппроксимации схемы С.К. Годунова. «ЖВМ и МФ», т.27, №12, 1987, стр. 1853-1860.

51. Акимов В.М., Бакулев В.И., Курзинер Р.И., Поляков В.В., Сосунов В.А., Шляхтенко С.М. Теория и расчет воздушно-реактивных двигателей. М.: Машиностроение, 1987.

52. J. Butterfield, Н. Yao, М. Price, С. Armstrong, S. Raghunathan,E. Benard, R. Cooper, D. Monaghan Enhancement of thrust reverser cascade performance using aerodynamic and structural integration // The Aeronautical Journal, December 2004

53. К.Вескег, J.Vassberg Numerical Aerodynamics in Transport Aircraft Design, Notes on Num. Fluid Mechanics, vol 100, pp 209-220

54. Keizo Takenaka, Kazuomi Yamamoto, Ryoji Takaki CFD validation study of NEXST-1 near Mach 1, ICAS 2004

55. Aurélia Cartieri, S. Mouton, G. Boyet Syudy of Support Iinterference Effects at S IMA Wind Tunnel within the "SAO" Project // ICAS 2010

56. Власенко В.В. О математическом подходе и принципах построения численных методологий для пакета прикладных программ EWT-ЦАГИ // Сборник статей. Труды ЦАГИ. 2007. Вып.2671. - С. 20-85.

57. Михайлов C.B. Объектно-ориентированный подход к созданию эффективных программ, реализующих параллельные алгоритмы расчета // Сборник статей. Труды ЦАГИ. 2007. Вып.2671. - С. 86-108.

58. С.М. Босняков, В.В. Власенко, И.А. Курсаков, C.B. Михайлов, Ю. Квест Задача интерференции оживального тела вращения с державкой аэродинамической трубы и особенности ее решения с использованием ЭВМ // Ученые записки ЦАГИ № 3 2011 г.

59. Dash S., Weilersteen G., Vaglio-Laurin R. Compressibility effects in free turbulent shear flows // TR-75-1436. AFOSR. 1975

60. Coakley T.J. Turbulence modeling methods for the compressible Navier-Stockes equations // AIAA-83-1693. 1983

61. Coakley T.J., Hsieh T. Comparison between implicit and hybrid methodsfor the calculation of steady and unsteady inlet flows // AIAA-85-1125. 1985

62. Шлихтинг Г. Теория пограничного слоя. М.: Наука, 1974

63. Коваленко В.В, Кравченко А.И, Старухин В.П, Щенников С.А. Применимость модели невязкого газа к описанию обтекания сверхзвуковым потоком тел большого удлинения. // Ученые записки ЦАГИ 1987, т. 17, №6 стр.23-29.

64. Под ред. Хемша М. и Нилсена Дж. Аэродинамика ракет. М.: Мир 1989.

65. Курсаков И.А. Опыт практического применения EWT к решению задачи расчета параметров потока на входе в силовую установку, установленную на подветренной стороне фюзеляжа JIA. «Труды ЦАГИ» Выпуск 2671 2007

66. Петров К.П. Аэродинамика тел простейших форм. Научное издание -М.: Факториал, 1998

67. Воеводин A.B., Прысев Б.Ф. Использование панельного метода расчета для исследования сходимости результатов испытаний в АДТ на ленточной подвеске и хвостовой державке. // Ученые записки ЦАГИ т. 39 №1-2 2008