Разработка эффективных комбинированных RANS/LES-методов для расчета сложных турбулентных струй тема диссертации и автореферата по ВАК РФ 01.02.05, кандидат технических наук Любимов, Дмитрий Александрович

• Любимов, Дмитрий Александрович
• кандидат технических наук
• 2008

Оглавление диссертации кандидат технических наук Любимов, Дмитрий Александрович

Введение

Глава 1 Выбор элементов для построения эффективных комбинированных RANS/LES-методов

1.1 Элементы численного метода, определяющие его точность

1.2 Влияние способа вычисления давления на точность метода

1.3 Требования, предъявляемые к разностным схемам для RANS и для LES

1.4 Модель турбулентности и SGS-модель

1.5 LES с неявной SGS-моделью

1.6 Выбор схем для построения эффективных комбинированных RANS/LES-методов

Глава 2. Численная реализация комбинированных RANS/LES-методов

2.1 Системы уравнений

2.2 Метод решения.

2.3 Модель турбулентности

2.4 Построение неявного оператора

2.5 Аппроксимация диффузионных членов уравнений

2.6 Аппроксимация конвективных членов уравнений Навье-Стокса для несжимаемых течений.

2.7 Аппроксимация конвективных членов уравнений Навье-Стокса для сжимаемых течений

2.8 Аппроксимация конвективных членов уравнений переноса

2.9 Метод DES для расчета несжимаемых течений

2.10 Комбинированный RANS/ILES-метод для расчета сжимаемых течений

2.11 Решение системы разностных уравнений

2.12 Повышение точности методов при расчетах на нерегулярных сетках

2.13 Граничные условия

2.13.1 Входная и выходная границы течения

2.13.2 Граничные условия вдали от струи для затопленных струй

2.13.3 Твердая стенка

2.13.4 Условие периодичности

2.13.5 Граница между блоками расчетной сетки

Глава 3. Расчет несжимаемых турбулентных струйных течений методом DES

3.1 Постановка задачи

3.2 Свободная прямоугольная струя

3.3 Пристеночная струя, истекающая из круглого сопла

Глава 4. Использование комбинированного RANS/ILES-метода для совместного расчета сжимаемых течений в круглых соплах и их струях.

4.1 Постановка задачи

4.2 Исследование влияния параметров потока на срезе круглого сопла на характеристики турбулентности в струе

4.3 Сверхзвуковые струи. Расчетное и нерасчетное истечение

Глава 5. Исследование с помощью RANS/ILES-метода влияния геометрии сопла на течение в струе.

5.1 Коническое шевронное сопло

5.2 Осесимметричное сопло двухконтурного ТРД

5.3 Сопло двухконтурного ТРД с шевронами на сопле газогенератора

5.4 Влияние малых изменений геометрии сопла двухконтурного ТРД на течение в струе

Введение диссертации (часть автореферата) на тему «Разработка эффективных комбинированных RANS/LES-методов для расчета сложных турбулентных струй»

В настоящее время наибольшее распространение для численного моделирования турбулентных течений жидкости и газа получил подход, при котором находится решение осредненных уравнений Навье-Стокса (Reynolds Averaged Navier-Stokes - RANS) с полуэмпирическими моделями турбулентности. Точность результатов, которую можно получить с помощью таких методов, в значительной степени зависит от типа течения. Так, при. расчете пограничных слоев отличие экспериментальных и расчетных результатов не превышает нескольких процентов. При расчете струйных течений с помощью RANS точность заметно ниже. Например, интенсивность расширения круглой затопленной струи в расчете с помощью лучших моделей турбулентности и в эксперименте различаются, на 30-300% [1]. Еще хуже модели турбулентности предсказывают многие «тонкие» эффекты, которые наблюдаются в экспериментах при анализе трехмерных струйных течений. В» частности, известно, что при развитии прямоугольной струи смешение: первоначально идет более интенсивно вдоль короткой стороны и поперечный« размер струи в этом направлении становится больше, чем вдоль длинной стороны. Это явление получило название «переворота осей», и его не удается описать с помощью традиционных моделей. Другой интересный эффект наблюдается при распространении круглой струи вдоль стенки. В этом случае смешение и растекание вдоль стенки идет значительно более интенсивно, чем в направлении перпендикулярном к стенке. Традиционные модели не описывают1 эти эффекты даже качественно, и только весьма сложные и специфические модификации моделей позволяют уловить эти эффекты.

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

Точность численного моделирования можно повысить, если вместо RANS и моделей турбулентности использовать методы прямого численного моделирования (ПЧМ). Использование методов ПЧМ и их различных вариантов (Direct Numerical simulation - DNS, Large Eddy Simulation — LES, Detached Eddy Simulation — DES) позволяет повысить точность и информативность описания турбулентных течений жидкости и газа. Благодаря развитию вычислительной техники и совершенствованию численных методов эти подходы находят все большее распространение. Работы, посвященные применению методов ПЧМ, появились более 20 лет тому назад. В нашей стране первые результаты по описанию с помощью ПЧМ перехода, от ламинарного течения к турбулентному в трубах были получены в группе профессора Б.Л. Рождественского в МИФИ. Заметный вклад в развитие ПЧМ внес коллектив под руководством академика О.Н. Белоцерковского в МФТИ. В последнее десятилетие важные результаты были получены в группах М.Х. Стрельца (ГИПХ, Санкт-Петербург) и Н.В. Никитина (НИИ Мех МГУ).

Анализ работ, посвященных ПЧМ, показывает, что наибольшее их число посвящено расчетам течений в пограничном слое и в каналах. Так, например, расчет пограничного слоя методом DNS был выполнен Ф. Спалартом еще в 1988 году. Позже его модификация метода LES — метод отсоединенных вихрей

DES) позволила существенно улучшить описание отрывных течений вблизи крылового профиля.

Моделирование струйных течений, в частности выхлопных струй ТРД, с помощью ПЧМ развивалось существенно медленнее, а точность таких расчетов весьма невысока. Между тем, как уже отмечалось выше, струйные течения плохо описываются и с использованием традиционных дифференциальных моделей турбулентности и RANS. По этой причине развитие методов ПЧМ для струй особенно актуально. Одни из первых попыток расчета затопленных турбулентных струй были описаны в работах [2, 3]. В последние годы большое внимание развитию ПЧМ для исследования в струйных течениях уделяется в Стенфордском центре (Lele, Bodony и др.), в университете Чалмерса (Davidson, Andersson и др.), в ГИПХе (Стрелец, Травин и др.). Оригинальный подход с использованием уравнений Эйлера был предложен А.Н. Крайко с сотрудниками.

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

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

На современном уровне развития вычислительной техники на все виды ПЧМ пригодны для практических приложений. Так метод DNS требует самой мелкой расчетной сетки. На современном уровне развития вычислительной техники он может быть использован для решения модельных задач только при сравнительно низких числах Рейнольдса. Для практических приложений возможно использование методов ПЧМ, использующих LES или комбинированные RANS/LES-методы. Особенности струйных течений, перечисленные выше, делают актуальным использование эффективного численного метода с высоким разрешением, который позволяет получить заданную точность результатов на сетке с относительно небольшим числом ячеек. Следует учитывать, что на точность результатов в методе LES также влияет используемая подсеточная (Sub-Grid Scale - SGS) модель турбулентности. Она отвечает за обмен энергией между крупными вихрями, которые разрешаются явно, и мелкими вихрями, размер которых меньше размера ячеек расчетной сетки. При выборе SGS-модели турбулентности нужно помнить, что те из них, что основаны на модели Смагоринского, в силу особенностей своего построения, удовлетворительно работают на сетках, которые близки к изотропным по всем координатным направлениям.

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

В настоящее время применяются два подхода к моделированию струйных течений с помощью LES. В первом случае выполняется расчет только струи, а параметры на срезе сопла задаются приближенно. Это позволяет избежать очень затратного расчета течения внутри сопла, но не дает возможности учесть реальное распределение параметров на срезе сопла. Струя становится нестационарной и турбулентной, как правило, только к концу начального участка, и течение в слое смешения не удается моделировать. Такую картину можно наблюдать в большинстве расчетов, выполненных с помощью LES, результаты которых приведены в обзорных публикациях [4, 5]. Кроме того, нарушается правильность описания физики течения на срезе сопла, поскольку в дозвуковом потоке возмущения могут передаваться вверх по потоку, а в такой постановке эта возможность отсутствует. Указанный упрощенный подход оправдан в случае, когда важно описать- основной участок струи, а точность описания течения около среза сопла не играет большой роли. К таким течениям, в частности, можно отнести упомянутые выше струи с деформацией поперечного сечения: пристеночную струю, вытекающую из круглого сопла, и свободную прямоугольную струю, в которых развитие течения определяется процессами, происходящими в самой струе. Течение в них не удается описать правильно даже качественно с помощью осредненных по Рейнольдсу уравнений Навье-Стокса и стандартных полуэмпирических моделей турбулентности [6, 7]. Исследование течения в первой из них особенно важно, поскольку ее результаты могут быть использовано при исследовании распространения выхлопной струи ТРД вдоль поверхности аэродрома.

При втором подходе осуществляется совместный расчет течения в сопле и в струе. В пределе это позволяет получить реальное распределение параметров течения на срезе сопла, что существенно для сопел сложной формы, которые характерны для современных ТРД. Появляется возможность описать турбулентность в слое смешения сразу за кромкой сопла, учесть влияние геометрии сопла на течение в струе. Известен ряд публикаций, где был использован такой подход. Так, в работах [5,8-10] приведены результаты совместного расчета истечения изотермических и нагретых струй из осесимметричных сопел. В этих работах числа M на срезе сопла были в диапазоне 0.6-0.9, числа Re 5.0х104 - 1.2х106. Применялся как традиционный метод LES с явными SGS-моделями турбулентности [5, 9], так и вариант метода с неявной SGS-моделью турбулентности. Такие методы в англоязычной литературе называются Monotonically Integrated LES - MILES, более позднее название: Implicit LES - ILES. Метод ILES был использован в работах [8, 10]. В перечисленных публикациях расчеты выполнялись при помощи различных конечно-разностных схем: малодиффузионных схем 3-го порядка с разностями против потока [8, 9], схем с центральными разностями 4-го порядка [5] и схем с компактными разностями 6-го порядка [10]. Самая грубая сетка содержала

6 7

1.6x10 ячеек, а самая мелкая - 5.0x10 ячеек [10]. Однако и при использовании совместного расчета течения в сопле и в струе авторам не удалось повысить точность предсказания по сравнению с расчетами струй с помощью LES без расчета течения в сопле. Так в работе [10], где была использована самая мелкая расчетная тетка из всех перечисленных, представлено распределение максимальной интенсивности турбулентности в слое смешения по длине струи. Расчетные значения этого параметра сильно разнятся с известными экспериментальными данными.

Примеры расчетов течения в соплах двухконтурных ТРД с центральным телом с помощью LES описаны в работах [11 — 13]. Расчеты выполнялись с помощью различных вариантов LES, в которых применялись разные схемы конечно-разностной аппроксимации, на сетках с числом ячеек 6.0x106

2.0x10 . Сравнение с экспериментом отсутствует, поэтому явно судить о точности полученных результатов затруднительно, но представленный иллюстративный материал косвенно свидетельствует о том, что точность расчета струи около среза сопла невысока.

Анализ известных результатов совместных расчетов течений в соплах и их выхлопных струях с помощью LES показывает, что в настоящее время эта задача не решена с достаточной для практических приложений точностью. Главная причина этого состоит в том, что даже на самых мелких из использованных сеток не удается описать вихревые структуры в пограничном слое на стенках сопла и в слое смешения около среза сопла. Можно сделать приближенную оценку снизу потребных расчетных сеток для совместного расчета течения внутри сопла и в слое смешения при условии разрешения вихрей в пограничном слое. В авиационном сопле типичная толщина пограничного слоя 8 составляет около 3% от диаметра сопла D. Для правильного описания турбулентности в пограничном слое нужно разрешить вихри с поперечным размером менее 0.25. При использовании разностной схемы высокого порядка для описания вихря требуется не менее 5 ячеек. Простые вычисления показывают, что в азимутальном направлении сетка должна содержать, по крайней мере, 2600 ячеек. Расчеты с помощью LES выполняются на сетках близким к однородным, поэтому число ячеек расчетной сетки по другим направлениям будет того же порядка. Пусть по каждому направлению расчетная сетка содержит 1000 ячеек. Общее число ячеек в этом случае будет более 109 ячеек. Известно, что в турбулентном пограничном слое максимум пульсаций продольной компоненты скорости находится на границе ламинарного подслоя. При больших числах Рейнольдса, которые характерны для сопел ТРД, толщина ламинарного подслоя мала, и размеры вихрей, которые необходимо разрешить, будут еще меньше, чем в приведенной выше оценке, что потребует еще более мелкой сетки.

Альтернативой, позволяющей существенно снизить вычислительные затраты, может быть использование комбинированных RANS/LES-методов.

Один из первых вариантов таких методов - метод отсоединенных вихрей (DES) — был предложен в работе [14]. В этом случае течение около стенок сопла описывается с помощью нестационарного метода RANS, что позволяет использовать более грубые сетки для описания пограничных слоев. При этом на выходе из сопла получается реалистическое распределение осредненных параметров течения. Мгновенное поле течения удается моделировать только на некотором расстоянии от среза сопла. Точность и работоспособность RANS/LES-методов, размеры переходной зоны за срезом сопла, в значительной степени зависят от удачного выбора численного метода, схемы разностной аппроксимации, SGS-модели турбулентности.

Известен ряд публикаций, посвященных расчетам струйных течений с использованием указанных методов. Так в работе [15] расчет течения в прямоугольном сопле и струе был выполнен с помощью метода DES [14]. Сетка содержала 4х106 ячеек. Число Рейнольдса Re было 2><105. Рассмотрены изотермическая и горячая струи с числом Мр0.9. Получено хорошее совпадение расчетных результатов с экспериментом для распределений осевых параметров: осредненной скорости и интенсивности турбулентности. В работе [16] с помощью комбинированного RANS/LES-метода были исследованы течения в коническом сопле и двухконтурном сопле с центральным телом. В первом случае сетка содержала 4X106 ячеек, во втором — 7.75х 10б ячеек. Для конического сопла совпадение осевой скорости с расчетными данными хорошее, но значения пульсаций скорости вдоль оси струи сильно расходятся с экспериментом. Для двухконтурного сопла не удалось получить удовлетворительного совпадения расчетных и экспериментальных данных. В перечисленных выше примерах за исключением [10], отсутствуют данные по уровню турбулентности в слое смешения струй при расчетах с помощью LES или комбинированных методов. Это косвенно свидетельствует о том, что указанными методами пока не удается описать характеристики турбулентности в слое смешения на начальном участке струй.

Анализ литературы показывает, что в настоящее время комбинированные методы находятся в стадии становления, и отсутствует ориентированный на практические приложения RANS/LES-метод для совместного расчета течений в соплах ТРД и их выхлопных струях на различных режимах.

Таким образом, актуальной задачей является повышение точности расчета струйных течений с помощью LES, совершенствование этих методов, создание эффективных комбинированных RANS/LES-методов для совместного расчета течения в соплах, включая сопла ТРД, и их выхлопных струях на различных режимах истечения, в том числе и сверхзвуковых. Для практики важно, чтобы эти методы были нетребовательны к качеству расчетных сеток, обладали высоким разрешением, позволяющим выполнять расчеты на сравнительно грубых сетках при сохранении приемлемой точности результатов.

Целью работы является разработка эффективных комбинированных RANS/LES-методов, позволяющих обеспечить высокую точность расчета турбулентных выхлопных струй различных сопел, включая сопла ТРД, при малых вычислительных затратах, и исследование с их помощью сложных турбулентных струй, в том числе сверхзвуковых струй со скачками уплотнения.

Основные задачи работы: анализ известных численных методов высокой точности для расчета сжимаемых и несжимаемых течений с целью выбора наиболее эффективных для создания комбинированных RANS/LES-методов расчета струйных течений; численная реализация метода DES высокого порядка разностной аппроксимации применительно к расчету низкоскоростных струйных течений; исследование с помощью разработанного варианта метода DES течения в свободной прямоугольной струе и пристеночной струе, вытекающей из круглого сопла; разработка эффективного комбинированного RANS/ILES-метода высокого порядка аппроксимации для расчета сжимаемых течений, включая течения со скачками уплотнении; оценка точности определения характеристик турбулентных выхлопных струй при совместном расчете течений в соплах разных типов и их струях с помощью разработанного RANS/ILES-метода; исследование влияния параметров течения на срезе сопла на характеристики турбулентности в струе; исследование влияния геометрии сложных сопел, включая сопла ТРД, на характеристики турбулентности выхлопных струй.

Диссертация состоит из введения, пяти глав, заключения и списка литературы.

Заключение диссертации по теме «Механика жидкости, газа и плазмы», Любимов, Дмитрий Александрович

Заключение

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

2. При расчете несжимаемых струйных течений предложены варианты методов ПЧМ с использованием метода, искусственной сжимаемости и высокого — 5-го порядка — разностной аппроксимации конвективных членов уравнений.

3. Для описания сжимаемых течений, включая течения со скачками1 уплотнения, разработан комбинированный 11АЫ8/Ш£8-метод, в котором при аппроксимации конвективных членов в уравнениях Навье-Стокса также удалось использовать монотонную схему 5-го порядка с дополнительно уменьшенной схемной вязкостью. Способность метода без дополнительных изменений описывать как

1 > ^ дозвуковые, так и сверхзвуковые течения повышает его практическую значимость.

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

5. Точность расчетов и надежная работа разработанных методов на нерегулярных криволинейных сетках обеспечивается использованием конечно-объемной схемы интегрирования по пространству, а также согласованным способом вычисления: параметров в схеме Роу; потоков в уравнении неразрывности и коэффициентов в конвективных членах уравнения переноса; производных в метрических коэффициентах и в диффузионных членах уравнений Навье-Стокса и переноса.

6. С помощью разработанных численных методов выполнены тестовые расчеты нескольких течений, проведено сопоставление с известными экспериментальными данными, анализ этих расчетов позволил получить некоторые новые результаты.

Так, с помощью метода ПЧМ для расчета низкоскоростных течений удалось правильно описать такие специфические особенности турбулентных струйных течений как переворот осей в прямоугольной струе и трансверсальное растекание пристеночной первоначально круглой струи. Последняя задача с помощью ПЧМ была решена впервые. Ее результаты могут быть использованы при исследовании распространения струи ТРД вдоль поверхности аэродрома.

7. Для пристеночной струи с помощью методов ПЧМ было получено значение константы Кармана близкое к 0,4. Это делает возможным использование «законов стенки» при-расчете турбулентных пристеночных течений с реальными числами Рейнольдса с помощью комбинированных RANS/LES-методов.

8. Разработанный RANS/ILES-метод для расчета до- и сверхзвуковых турбулентных течений при больших числах Рейнольдса был применен к совместному расчету течений в соплах разных типов (коническом, коническом шевронном, сопле двухконтурного ТРД, сопле двухконтурного ТРД с шевронами на сопле газогенератора) и их выхлопных струях. Общее число ячеек расчетной сетки в этих расчетах было 0.8—1.3*106. При этом число ячеек расчетной сетки внутри сопла и вне его выше по потоку от среза сопла не превышало 16% от их общего числа. Точность результатов не уступала, а в большинстве случаев i превосходила достигнутую другими авторами на существенно более мелких к сетках.

9. Во всех рассмотренных случаях расчетов с помощью RANS/ILES-метода удалось получить хорошее совпадение с имеющимися экспериментальными данными, как по осредненным, так и по пульсационным характеристикам течения. Однако, на малых расстояниях (X/D<2) уровень пульсаций скорости и давления в слое смешения превышает на 10-20 % экспериментальные данные. Это расхождение связано со специфической особенностью комбинированных RANS/LES-методов. г

10.Для одноконтурных конических сопел было исследовано влияние акустического числа Маха Ма и температуры струи на срезе сопла на уровень пульсаций скорости в слое смешения струи. Установлено, что изменение температуры в диапазоне ^/Тоо=0.9-2.7 и числа Маха Ма=0.5-0.9 не приводит к заметному изменению уровня пульсаций. Полученный при этих расчетах уровень пульсаций давления на оси струи и в слое смешения близок к известным экспериментальным данным.

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

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

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

14.Для сопла двухконтурного ТРД с центральным телом с шевронами на сопле газогенератора даже на сетке с 1.3х106 ячеек удалось правильно описать.влияние шевронов на течение в струе, что подтверждается сравнением с экспериментальными данными.

15.В сопле двухконтурного ТРД с центральным телом было исследовано влияние несоосности сопел газогенератора и вентилятора. Установлено, что уровень энергии турбулентности в струе чувствителен даже к малой несоосности. Так, при несоосности 0.040с из-за более раннего слияния слоев смешения струй газогенератора и вентилятора, положение максимума энергии турбулентности смещается с ХЯ)С=10.0 с ХЛ)С=8.0. Пик энергии турбулентности увеличивается на 25-30% по сравнению с осесимметричным соплом. При этом азимутальная неравномерность в распределении энергии турбулентности достигает 70%.

108

Список литературы диссертационного исследования кандидат технических наук Любимов, Дмитрий Александрович, 2008 год

1. Bardina J.E., Huang P.G., Coakley T.J. Turbulence Modeling Validation, Testing and Development // NASA TM 110446. - 1997.

2. Freund J.B. Direct Numerical Simulation of the Noise from a Mach 0.9 jet // ASME P. 1999. -N. FEDSM99-7251.

3. Fureby C., Grinstein F.F. Motonically Integrated Large Eddy Simulations of Free Shear Flows // AIAA J. 1999. - V.37. - N. 5 - P.544-556.

4. Bodony D.J., Lele S.K. Review of the Current Status of Jet Noise Predictions Using Large Eddy Simulation (invited) // AIAA P. 2006. - N2006-468.

5. DeBonis J.R. Progress towards Large Eddy Simulation for Prediction of Realistic Nozzle Systems // AIAA P. 2006. - N2006-487.

6. Берч С.Ф., Лебедев А.Б., Любимов Д.А., Секундов А.Н. Моделирование трехмерных струйных и погранслойных течений // Изв. РАН. МЖГ. — 2001 -№5-С. 48-63.

7. Khritov К.М., Lyubimov D.A., Maslov V.P.,Mineev B.I.,Secundov A.N., Birch S.F. Three-dimensional wall Jets: Experiment, Theory and Application // AIAA H. 2002. - N2002-0732.

8. Lupoglazoff N., Biancherin A., Vuilot F., Rahier G. Comprehrensive 3D Unsteady Symulations of Subsonic and Supersonic Hot Jet Flow-Fields. Part 1: Aerodynamic Analysis // AIAA P. 2001. N2002-2599.

9. Anderson N., Eriksson L.-E., Davidson L. A Study of Mach 0.75 Jets and Their Radiated Sound Using Large Eddy Simulation // AIAA P. 2005. -N2004-3024.

10. Uzun A., Hussiani M.Y. High Frequency Noise Generation in the Near-Nozzle Region of a Jet // AIAA P. 2006. - N2006-2499.

11. Anderson N., Eriksson L.-E., Davidson L. LES Prediction of Flow and Acoustic Field of a Coaxial Jet // AIAA P. 2005. - N2005-2884.

12. Mihaesku M., Gutmark E., Szasz R-Z., Fuchs L. Flow and Acoustics of a Nozzle: a Sensitivity Study to the Inlet Boundary Conditions // AIAA P. — 2006.-N2006-619.

13. Tristano I.H., Page G.J., McGuirk J.J. Large Eddy Simulation of Hot Coaxial Jets // AIAA P. 2006. - N2006-2497.

14. Spalart P.R., Jou W.-H., Strelets M., Allmaras S.R. Comments on the feasibility of LES for wings, and on a hybrid RANS/LES approach // First AFOSR International conference on DNS/LES. Rouston, Louisiana. — USA. — 1997. — P.669-687.

15. Paliath U., Morris P.J. Prediction of Jet Noise from Rectangular Nozzles // AIAA P. -2006. -N2006-618.

16. Abdol-Hamid Kh.S., Elmiligui A. Numerical Study of High Temperature Jet Flow Using RANS/LES and PANS Formulation // AIAA P. 2005. - N2005-5092.

17. Patankar S.V. and Spalding D.B. A Calculation Procedure for Heat, Mass and Momentum Transfer in Three Dimensional Parabolic Flows // Intern. J. of Heat and'Mass Transfer. 1972. - V.15. - P.1787-1806.

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

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

20. Крайко А.Н., Пьянков К.С. Течения идеального газа с отрывными зонами и нестационарными контактными разрывами сложной формы // Изв. РАН. МЖГ. 2006. №5. - С.41-54.

21. Годунов С.К. Разностный метод численного* расчета разрывных решений уравнений гидродинамики // Матем. сб. 1959. - Вып.47(89). — №3. -С.271-306.

22. Roe P.L. Approximate Riemann Solvers, Parameters Vectors, and Difference Schemes // J. of Сотр. Phys. 1981. - V.43. - P.357-272.

23. Van Leer B. Flux-Vector Splitting for the Euler Equations // Lecture Notes in Physics. 1982. - V.170. - P.507-512.

24. Liou M.-S., Steffen С JJr. A new Flux Splitting Scheme // J. of Сотр. Phys. -1993.-V.107-P.23-39.

25. Sukumar R., Chakravarthy S. R., Szema K.-Y. Euler Solver for Three-Dimensional Supersonic Flows with Subsonic Pockets // J. of Aircraft. — 1980. — V.24. — N2. P.73-83.

26. Васильев В.И., Волков Д.В., Любимов Д.А. Использование однопараметрической дифференциальной модели турбулентности в численных расчетах с помощью уравнений Навье-Стокса // ТВТ. — 1998. -Т.36.-Ш.-С.65-73.

27. Chorin A .J. A Numerical Method for Solving Incompressible Viscous Flow Problems // J. of Сотр. Phys. 1967 - V.2. - P.12-26.

28. Rogers S.R., Kwak D. Upwind Differencing Scheme for the Time-Accurate Incompressible Navier-Stokes Equations // AIAA J. 1990. - V.28. - N2. - P. 253-262.

29. Rogers S.R., Kwak D, Kiris C., Steady and Unsteady Solution of the Incompressible Navier-Stokes Equations // AIAA J. 1991 - V.29. - N4. - P. 603-610.

30. Козлов B.E., Любимов Д.А.,. Секундов, А.Н, Спаларт Ф.Р. Трансверсальное распространение турбулентности в пограничном слое. // Изв. РАН. МЖГ. -1998. №3. - С.77-84.

31. Lebedev А. В., Lyubimov D.A., Maslov., Mineev B.I., Secundov A. N. The Prediction of Three-Dimensional Jet Flows for Noise Applications // AIAA P. -2002. N2002-2422.

32. Козлов B.E., Лебедев А.Б., Любимов Д.А., Секундов А.Н. Некоторые особенности турбулентного течения в кромочном вихре // Изв. РАН. МЖГ. 2004. - №1. - С.78-85.

33. Birch S.F., Lyubimov D.A., Secundov A.N., Yakubovsky K.Ya. Numerical Modeling Requirements for Coaxial and Chevron Nozzle Flows // AIAA P. — 2003.-N2003-3287.

34. Birch S.F 2004, Lyubimov D.A., Secundov A.N., Yakubovsky K.Ya. Accuracy Requirements of Flow Inputs for Jet Noise Prediction Codes // AIAA P.-2004.-N2004-2934.

35. Birch S.F., Lyubimov D.A., Buchshtab P.A., Secundov A.N., Yakubovsky K.Ya. Jet-Pylon Interaction Effects // AIAA P. 2005. - N2005-3082.

36. Birch S.F., Lyubimov D.A., Maslov V.P., Secundov A.N. Noise Prediction for Chevron Nozzle Flows// AIAA P. 2006. - N2006-2600.

37. Forsythe J.R., Hoffmann K.A., Cummings R.M., Squires K.D. Detached-Eddy Simulation With Compressibility Corrections Applied to Supersonic Axisymmetric Base Flow // Trans. ASME, J. Fluid Engineering — 2002 — V. 124.-P.1-13.

38. Basu D., Hamed A., Das K., DES and^ Hybrid RANS/LES Models for Unsteady Separated Turbulent Flow Predictions // AIAA P. 2005. - N2005-503.

39. Strelets M. Detached Eddy Simulation of Massively Separated Flows // AIAA P.-2001.-N2001-879.

40. Tristano I.H., Li Q., Page G.J., McGuirk J.J. On the Effect of Convective Flux Formulation for LES of Compressible Flows Using Hybrid Unstructured Meshes // AIAA P. 2006. - N2006-3739.

41. Rider W.J. Margolin L. From Numerical Analysis to Implicit Subgrid Turbulence Modeling // AIAA P. 2003. - N2003-4101.

42. Mittal, R. and Moin, P., Suitability of Upwind-Biased Finite-Difference Schemes for Large-Eddy Simulation of Turbulent Flows // AIAA J. — 1997. -V. 35. — P.1415-1417.

43. Shur M.L., Spalart R., Strelets M.Kh. Noise Prediction for Increasingly Complex Jets. Part I: Methods and Tests // Int. J. of Aeroacousics. 2005. - V. 4. -N3&4. -P.213-246.

44. Menter F.R., Kuntz M., Bender R. A Scale-Adaptive Simulation Model for Turbulent Flow Prediction // AIAA P. 2003. - N2003-767.

45. Tucker P J. Davidson L. Zonal k-l based Large Eddy Simulations // AIAA P. -2003.-N2003-82.

46. Fan C.-C., Xiao X., Edwards J.R., Hassan H.A., Baurle R.A. Hybrid Large-Eddy/Reynolds-Averaged Navier-Stokes Simulation of Shock-Separated Flow // J. of Spacecraft and Rockets. 2004. - V.41. - N6. - P.897-906.

47. Chaouat B., Schiestel R. A New Partially integrated transport model for subgrid-scale stress and dissipation rate for turbulent developing flows // Phys. of Fluids. 2005 - V. 17. - article 065106.

48. Girimaji S.S. Partially-Averaged Navier-Stokes Model for Turbulence. A Reynolds-Averaged Navier-Stokes to Direct Numerical Simulation Bridging Method // Trans, of ASME, J. of Applied Mechanics. 2006. - V. 73. -R413-421.

49. Batten P., Goldberg U., Chakravarthy S. LNS An Approach towards Embedded LES // AIAA P. - 2002. - N2002-427.

50. Hedges L.S., Travin A.K. Spalart P.R. Detached-Eddy Simulations Over a Simplified Landing Gears // Trans, of ASME, J. of Fluid Engineering. — 2002. V. 124. -P.413-423.

51. Gerolymos G.A., Senechal D., Vallet I. Reynolds-StressModel-VLES Multiblock Implicit Solver using High-Order Upwind Schemes // AIAA P. -2006. N2006-3909.

52. Spalart P.R. ,Allmaras S.R.A One-Equation Turbulence Model for Aerodynamic Flows //La Recherche Aerospatiale. 1994. - № 1. - P. 5-21.

53. Nikitin N.V. Nicoud F., Wasisstho b., Squires K.D., Spalart P.R. An Approach to Wall Modeling in Large-Eddy Simulations // Physics of Fluids. — 2000. — V.12. N7. - P. 1629-1632.

54. Boris J.P., Grinstein F.F., Oran E.S., Kolbe R.J., New Insights into Large Eddy Simulations // Fluid Dynamics Research. 1992. - N10. - P. 199-228.

55. Grinstein F.F., Fureby C., DeVore C.R. On MILES based on Flux-Limiting Algorithms // Int. J. for Numerical Methods in Fluids. 2005. - V.47. -P.1043-1051.

56. Grinstein F.F., Fureby C. Recent Progress on Flux-Limiting Based Implicit Large Eddy Simulation. In European, Conference on Computational Fluid Dynamics ECCOMAS CFD 2006. 2006. - (Eds. P.wesseling, E. Onate, J. Peraux). -.Delft. - The Netherland.

57. Tucker P.G. Hybrid MLLES-RANS Method for more Dissipative Solvers and the use of Non-Linear LES // AIAA P. 2004. - N2004-71.

58. Smagorinsky J. General Circulation Experiments with the Primitive Equations // Monthly Weather Review. 1963. - V.91. - N3. - P.99-165.

59. Rider W.R. Effective Subgrid Modeling from the ILES Simulation of Compressible Turbulence. In European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006. 2006. - (Eds. P.wesseling, E. Onate, J. Peraux). - Delft. - The Netherland.

60. Vasiliev V.I., Volkov D.V., Zaitsev S.A. Lyubimov D.A. Numerical Simulation of Channel Flows by a One-Equation Turbulence Model // Trans. ASME J. Fluid Engineering. 1997. - V.119. -P.885-892.

61. Travin A., Shur M., Strelets M., Spalart P. Detached-Eddy Simulations past a Circular Cylinder // Flow, Turbulence ans Combustuin. 1999. — V.63. -P.293-313.

62. Gamier E., Mossi M., Sagaut P., Comte P., Devolle M. On the use of Shock-Capturing Schemes for Large-Eddy Simulation // J. Comp. Phys. 1999. - V. 153. — P.273-311.

63. Suresh A., Huynh H. T. Accurate Monotonicity—Preserving Schemes with Runge-Kutta Time Stepping // J. Сотр. Phys. 1997. - V.136. - P.83-99.

64. Корецкий B.B., Любимов Д.А. Модифицированный метод приближенной факторизации для расчета потенциальных пространственных течений в каналах // ЖВМ и МФ. 1990. Т.ЗО. - №10. - С.1553-1570.

65. Любимов Д.А. Эффективный метод расчета пространственного обтекания воздухозаборников на околозвуковых скоростях // ЖВМ и МФ. 1991. — Т.31. -№9. - С.1355-1368.

66. Любимов Д. А. Возможности использования прямых методов для: численного моделирования турбулентных струй- // Аэромеханика и газовая динамика. 2003; - №3. — С. 14-20.

67. Любимов ДА. Разработка и применение эффективного RANS/ILES-метода для расчета сложных турбулентных струй // — ТВТ. — 2008. — Т.46. — №2: — 12С.

68. Bui Т.Т. A Parallel, Finite-Volume Algorithm for Large-Eddy Simulation* of Turbulent Flows// NASA TM-206570. -1999.

69. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика: В 10 т. Т. VI, Гидродинамика М.: Наука, 1988. - 733с.

70. Bridges J., Brown С. A. Parametric Testing of Chevrons on Single Flow Hot Jets //AIAAP: 2004. - N2004-2824.

71. Birch S.F. A Review of Axisymmetric Jet Flow Data for Noise Applications // AIAAP. 2006. - N2006-2602.73..Bridges J., Wernet M. Measurements of the Aeroacoustic Sound Source in Hot Jets // AIAAP. 2003. -N2003-3130.

72. Lau J.C., Morris P.J., Fisher M.J. Measurements in Subsonic and Supersonic Free Jets Using a Laser Velocimeter // J: Fluid Mech. 1979. - V.93. -Nl. -P.'1.-27.

73. Arakeri V.Y., Krothopalli A., Siddavaram V., Alkislar M.B., Lourendo M. On the Use ofMicrojets to Suppress Turbulence in a Mach 0.9 Axisymmetric Jet // J. Fluid Mech. 2003. - V.490. -P.75-98.

74. Uzun A., Lyrintzis A.S., Blaisdell G.A. Coupling of Integral Acoustic Methods with LES for Jet Noise Prediction // Int. J. of Aeroacousics. 2005. - V.4. -N3-4. — P.297-346.

75. Lau J.C. Effects of Exit Mach Numbers and Temperature on Mean-Flow and Turbulence Characteristics in Round Jets // J. Fluid Mech. — 1981. — V.105. — P.193-218.

76. Абрамович Г.Н., Гиршович T.A., Крашенинников С.Ю., Секундов А.Н., Смирнова И.П. Теория турбулентных струй. М.: Наука, 1984. — 716с.

77. Jones B.G., Adrian R.J., Nithianandan С.К., Plachon H.P. Spectra of Turbulent Static Pressure Fluctuations in Jet Mixing Layers // ALAA J. 1979: - V.17. — N5. - P.449-457.

78. Engblom W.A., Khavaran, Bridges J. Numerical Prediction of Chevron Noise Reduction Using WIND-MGBK Methodology // ALAA P. 2004. N2004-2979.

79. Massey S.J., Thomas R.H., Abdol-Hamid K.S., Elmiligui A.A. Computational and Experimental Flow Field Analyses and Pylon Interaction // AIAA P. — 2003.-N2003-3212.

80. Doty M.J., Henderson B.S., Kinzie K.W. Turbulent Flow Field Measurements of Separate Flow Round and Chevron Nozzles with Pylon Interaction Using Particle Image Velocimetry // AIAA P. 2004. - N2004-2826.

81. Thomas R.H., Kinzie K.W., Pao S.P. Computational Analysis of a Pylon-Chevron Core Nozzle Interaction // AIAA P. 2001. - N2001-2185.

82. Birch S.F., Khritov К. M., Maslov V.P., Mironov A. K., Secundov A.N. An Experimental Study of Flow Asymmetry in Co-axial Jets // AIAA P. 2005. -N2005-2845.