Расчет разрушения вихревого следа за пассажирским

advertisement
ТРУДЫ МФТИ. — 2014. — Том 6, № 3
Аэромеханика
3
УДК 532.527
И. С. Босняков, Г. Г. Судаков
Центральный аэрогидродинамический институт им. профессора Н. Е. Жуковского
Расчет разрушения вихревого следа за пассажирским
самолетом с помощью метода моделирования больших
вихрей второго порядка аппроксимации
В работе обсуждаются физические эффекты, имеющие место при разрушении дальнего вихревого следа за самолетом. Представлены результаты расчета эволюции вихревого следа в турбулентной атмосфере с помощью метода моделирования больших
вихрей (LES). Описаны математические модели для задания начальных условий, обсуждаются проблемы задания граничных условий. Проведены методические расчеты
для выбора параметров сетки. Показано, что метод моделирования больших вихрей
LES описывает все наблюдаемые физические эффекты в вихревом следе без какихлибо дополнительных предположений, но с разной точностью.
Ключевые слова: вычислительная аэродинамика, моделирование больших вихрей, вихревой след, разрушение вихревого следа, численное моделирование, неустойчивость Кроу, турбулентность.
1.
Введение
При полете в атмосфере летательный аппарат создает вихревой след — когерентную
структуру, состоящую из пары вихрей с противоположным вращением потока. Вихревой
след достаточно долго сохраняется в турбулентной атмосфере и представляет опасность
для последующих самолетов. В воздушной зоне, близкой к аэропорту, количество самолётов может быть велико. Увеличение расстояния между самолетами из-за опасности попадания в вихревой след снижает пропускную способность аэропорта. Так как пассажиропоток с течением времени растет, возникает очень важная практическая и теоретическая
задача уменьшения предписанного безопасного расстояния между самолетами, что требует
детального знания процесса разрушения вихревого следа в турбулентной атмосфере.
Рис. 1. Фазы разрушения дальнего следа за самолетом В-757 при уровне турбулентности атмосферы
𝑞 = 0.5 м/с и масштабе турбулентности 𝐿 = 200 м. Расчет с помощью инженерной модели вихревого
следа [1], [9]
Процесс разрушения дальнего следа включает в себя следующие физические эффекты.
ˆ
Синусоидальная неустойчивость Кроу [2] (описывается линейной теорией устойчивости Кроу–Бейта). Этот вид неустойчивости приводит к касанию вихрей при
и образование вихревых колец.
𝑇 = 𝑇𝑙𝑖𝑛𝑘
Аэромеханика
4
ˆ
ТРУДЫ МФТИ. — 2014. — Том 6, № 3
«Потеря» циркуляции. Обусловлена взаимодействием вихрей следа и фоновой атмосферной турбулентности (перекачка продольных вихрей следа на периферию посредством их дробления и взаимодействия с фоновой турбулентностью).
В экспериментальных исследованиях выделяют две фазы разрушения следа (рис. 1).
Сначала идёт фаза турбулентной диффузии. Здесь темп диффузии определяется фоновой
турбулентностью атмосферы и, по-видимому, связан с дроблением вихрей от самолёта,
вихрями атмосферной турбулентности. Затем начинается фаза быстрого разрушения. В ней
темп диффузии определяется параметрами вихрей следа. Эта фаза характерна только для
малых уровней турбулентности атмосферы.
Так как существенным фактором, определяющим механизм разрушения дальнего следа, является турбулентность атмосферы, имеющая случайную природу, то прямое сравнение расчета с экспериментальными данными, полученными для единичного случая эволюции следа, не является окончательным доказательством правильности метода расчета.
Экспериментальные результаты для реального полета, полученные с помощью лидаров [6],
показывают значительные колебания по времени параметров вихря (циркуляции, максимальной скорости, положения вихрей). Для выхода из этой ситуации можно рекомендовать
использование данных, полученных с помощью инженерной модели следа, учитывающей
все необходимые физические факторы и опирающейся на усредненные данные по многим
реализациям, полученные как экспериментальным путем, так и с помощью расчетов методами вычислительной аэродинамики. Такая модель была создана в ЦАГИ [9] и признана
за рубежом [1]. На рис. 2–4 приведены основные параметры следа для самолета В-757,
полученные с помощью этой модели.
Рис. 2. Максимальная окружная скорость в ядре вихря для самолета В-757. Расчет с помощью
инженерной модели вихревого следа [9], [1]
В настоящей работе в качестве инструмента исследований вихревого следа использован метод моделирования больших вихрей (LES), который должен описывать все перечисленные выше физические эффекты без каких-либо дополнительных предположений. Этот
подход наравне с другими методами является сейчас одним из инструментов изучения вихревого следа.
При использовании LES для решения задачи о вихревом следе необходимо учитывать
влияние диссипативных свойств выбранного численного метода. Известно, что схемы второго порядка аппроксимации обладают существенной схемной вязкостью. Приемлемую
диссипацию можно достичь либо с использованием более подробных расчётных сеток, либо
с помощью использования расчётных схем высокого порядка аппроксимации. Схемы высокого порядка аппроксимации не полностью решают проблему схемной вязкости, так как
диссипация остается достаточно высокой для высокочастотных пульсаций [8].
За последнее время произошёл резкий рост производительности вычислений в связи с
широким распространением компактных суперЭВМ, поэтому разумным выходом из сло-
ТРУДЫ МФТИ. — 2014. — Том 6, № 3
И. С. Босняков, Г. Г. Судаков
5
жившейся ситуации может оказаться использование существенно более подробных расчётных сеток в сочетании со стандартными методами второго порядка аппроксимации.
Рис. 3. Время до касания вихрей для B-757. Расчет с помощью инженерной модели вихревого следа
[9], [1]
Рис. 4. Средняя циркуляция (5–15 м) вихрей следа B-757 в зависимости от времени. Расчет с
помощью инженерной модели вихревого следа [9], [1]
Имеется достаточно много работ в данном направлении, в которых использовались как
методы расчета высокого порядка аппроксимации (см., например, [8], [7], [5]), так и методы
второго порядка ([4]). Преимуществом настоящей работы является то, что она дополнена
анализом недостатков и ограничений метода моделирования больших вихрей применительно к задаче о вихревом следе, а также указывает на то, какие эффекты данный подход
описывает точно, а какие – недостаточно точно.
В настоящей работе предлагается использовать численный метод конечного объема второго порядка аппроксимации (промышленный пакет ANSYS CFX) на структурированных
сетках размером более 50 млн узлов. Первая часть работы посвящена исследованию диссипативных свойств численного метода на примере плоской нестационарной задачи о диффузии пары вихрей противоположного знака. Здесь дано описание структуры и положения пары вихрей за летательным аппаратом, для чего используется инженерная модель,
разработанная в НАСА [7]. На примере задачи об эволюции пары вихрей в плоскости произведены оценки для выбора расчётной сетки, так чтобы численная вязкость метода могла
считаться приемлемой. Во второй части работы описан генератор фоновой турбулентности
атмосферы, который обеспечивает задание начальных условий для задачи о разрушении
следа. Суперпозиция скоростей от пары вихрей и фоновой турбулентности позволяет задать начальные условия для трехмерной нестационарной задачи об эволюции пары вихрей
в турбулентной атмосфере, которая рассматривается в третьей части работы. Задача о
разрушении вихревого следа в турбулентной атмосфере решается с помощью метода мо-
Аэромеханика
6
ТРУДЫ МФТИ. — 2014. — Том 6, № 3
делирования больших вихрей LES с использованием стандартной модели турбулентности
Смагоринского и со специальными граничными условиями, справедливость которых обсуждается в настоящей статье. Показано, что данный подход описывает все перечисленные
выше механизмы разрушения следа без введения каких-либо дополнительных предположений, но с разной степенью точности для каждого из механизмов.
2.
Инженерная модель Проктора
Инженерная модель Проктора [7] для двухвихревой системы следа непосредственно за
самолетом (ближний след) на режиме посадки описывается следующими формулами:
4𝐺
,
Γ=
𝜋𝑙𝜌𝑉пол
(︁
)︁]︁ [︁
(︁
)︁]︁
Γ [︁
𝑉𝜏 = 1,4
1 − exp −10 (𝑟𝑐 /𝑙)0,75
1 − exp −1, 2527 (𝑟/𝑟𝑐 )2 ,
2𝜋𝑟
(︁
)︁]︁
Γ [︁
𝑉𝜏 =
1 − exp −10 (𝑟/𝑙)0,75 , 𝑟 > 𝑟𝑐 ,
2𝜋𝑟
𝑏=
𝜋𝑙
,
4
𝑟 < 𝑟𝑐 ,
𝑏 — размах вихрей, 𝑙 — размах крыла, 𝐺 — вес самолета, 𝜌 — плотность воздуха, 𝑉пол —
скорость полета, 𝑟𝑐 — радиус ядра вихря. В работе [7] величина 𝑟𝑐 не определена, и ее
где
рекомендуется брать из эксперимента.
Параметры вихревого следа, используемые в данной работе, взяты для самолёта B-757:
расстояние между центрами вихрей
вихря
3.
𝑟𝑐 = 1,8
𝑏 = 38,05
м, циркуляция
Γ = 335
2
м /с, радиус ядра
м.
Двумерная задача об эволюции двух вихрей в отсутствиe фоновой
турбулентности
Рассмотрим
0 < 𝑦 < 300
двумерную
нестационарную
задачу
об
эволюции
−100 < 𝑧 < 100 м. В качестве среды используется
𝑇 = 25 ∘ C. Уравнение энергии не используется.
м,
плотности при
следа
в
области
воздух постоянной
Размеры области определяются:
ˆ
необходимым расстоянием от течения с большими градиентами до внешних границ
расчётной области,
ˆ
размахом вихрей,
ˆ
расстоянием по вертикали, которое вихревая пара преодолевает за время своей жизни.
На всех границах использовано граничное условие симметрии. Начальное условие представлено полем скоростей двухвихревого следа за самолетом B-757 (см. п. 2). Решается
двумерная задача об эволюции двух вихрей при следующих предположениях:
ˆ
отсутствует фоновая турбулентность,
ˆ
используются уравнения Навье–Стокса для несжимаемой жидкости с ламинарной
вязкостью.
Шаг по времени выбирался из условия, чтобы максимальное число Куранта для каждой сетки не превышало величины
𝐶𝑢 < 1.
𝐶𝑢 < 2,
а среднее число Куранта по всем ячейкам
Результаты тестирования приведены на рис. 5. Результаты расчета максимальной
окружной скорости в ядре вихря должны лежать выше данных, полученных по инженерной модели следа для минимального уровня турбулентности.
Методы второго порядка аппроксимации требуют очень мелких сеток для корректного
описания данной задачи. Из рис. 5 видно, что даже при переходе от сетки
более густой
1800 × 1200
1200 × 800
к
сходимость по сеткам еще не достигается. Однако при работе с
ТРУДЫ МФТИ. — 2014. — Том 6, № 3
трехмерной задачей разбиение по оси
сеткам
И. С. Босняков, Г. Г. Судаков
7
𝑥 (не менее 100 ячеек) уже приводит к очень большим
1800 × 1200 × 100 = 96 · 106 .
Рис 5. Максимальная окружная скорость в вихре. Данные расчётов для разных сеток и численных
схем в сравнении с инженерной моделью ЦАГИ [1]
Таким образом, оптимальным будем считать выбор разбиения не более чем
1200 × 800.
Однако следует иметь в виду, что для этого разбиения схемная вязкость является доминирующей внутри ядра вихря. Переход от схемы второго порядка аппроксимации к схеме
третьего порядка в уравнении импульса (график Fluent_QUICK на рис. 5.) не дает существенного улучшения точности.
Переход от уравнений Навье–Стокса (ламинарная вязкость) к уравнениям Рейнольдса
(модель турбулентности SST с поправкой на кривизну линий тока) в данной задаче показал,
что рассчитанная максимальная скорости в ядре вихря практически совпадает для обеих
моделей. Это служит дополнительным доказательством того, что в ядре вихря схемная
вязкость для выбранных сеток доминирует над молекулярной и турбулентной вязкостью.
4.
Метод генерации турбулентного поля
Исходными данными для задания случайных скоростей турбулентных порывов ветра
являются форма спектральной плотности и два параметра – масштаб и дисперсия. В данной
работе используется модель спектральной плотности Кармана.
Энергетическая функция спектральной плотности предложена в [3] и имеет специальный вид:
55 2
(𝛼𝑘 𝐿Ω)4
𝜎 𝐿
,
27𝜋
[1 + (𝛼𝑘 𝐿Ω)2 ]17/6
⟨︀ ′2
⟩︀
2
′2
′2 — среднеквадратическое отклогде 𝐿 — масштаб турбулентности, 𝜎 = 𝑢 + 𝑣 + 𝑤
нение, Ω = (Ω𝑥 , Ω𝑦 , Ω𝑧 ) — вектор пространственной частоты, 𝛼𝑘 = 1,339 — постоянная
Кармана. Волновое число Ω𝑜 , соответствующее масштабу турбулентности 𝐿, определяется
𝐸(Ω) =
из вышеприведенной формулы следующим образом:
Ω0 𝐿 = 1.
Функция
𝐸(Ω)
задаёт величину компонент случайного спектрального вектора скоро-
сти порыва ветра. Корреляция и дисперсия этих компонент выражается корреляционным
тензором
𝑆𝑖𝑘 (Ω).
Он определяется по формуле [11]
(︀
)︀
𝐸(Ω) Ω2 𝛿𝑖𝑘 − Ω𝑖 Ω𝑘
𝑆𝑖𝑘 (Ω) =
·
.
4𝜋
Ω4
Аэромеханика
8
Следует подчеркнуть, что в
ТРУДЫ МФТИ. — 2014. — Том 6, № 3
Ω-пространстве
корреляция между частотами отсутствует
(следствие однородности и изотропности турбулентности), но имеются корреляции меж-
Ω-пространстве
ду фурье-компонентами скорости в отдельно взятой точке в
(следствие
уравнения неразрывности).
Наиболее общим способом задания реализаций является представление скоростей ветра
в виде рядов Фурье по дискретному набору частот:
𝑊𝑗 (𝑥, 𝑦, 𝑧) =
𝑀 ∑︁
𝐾 ∑︁
𝑁
∑︁
𝑗
𝐴𝑗𝑚𝑘𝑛 cos(2𝜋·(Ω𝑥𝑚 𝑥+Ω𝑦𝑘 𝑦+Ω𝑧𝑛 𝑧))+𝐵𝑚𝑘𝑛
sin(2𝜋·(Ω𝑥𝑚 𝑥+Ω𝑦𝑘 𝑦+Ω𝑧𝑛 𝑧)),
𝑚=1 𝑘=1 𝑛=1
где
𝑗
𝐴𝑗𝑚𝑘𝑛 , 𝐵𝑚𝑘𝑛
— независимые по индексам
𝑚, 𝑘 , 𝑛
гауссовские случайные величины,
описываемые корреляционной (ковариационной) матрицей
[︁
]︁
[︁
𝑆𝑗𝑙 (Ω) (𝑗, 𝑙 = 1, 2, 3):
]︁
𝑗
𝑙
𝑀 𝐴𝑗𝑚𝑘𝑛 𝐴𝑙𝑚𝑘𝑛 = 𝑀 𝐵𝑚𝑘𝑛
𝐵𝑚𝑘𝑛
= 𝑆𝑗𝑙 (Ω𝑥𝑚 , Ω𝑦𝑘 , Ω𝑧𝑛 )ΔΩ𝑥 ΔΩ𝑦 ΔΩ𝑧 ,
где
𝑀 [𝐴]
— математической ожидание случайной величины
𝐴.
При генерации поля описанным выше способом стоит учитывать, что число операций
и, как следствие, время их выполнения растёт по закону
𝑁 6,
где
𝑁
— характерное число
разбиений по каждой оси в физическом и в волновом пространствах.
В качестве примера на рис. 6 представлено поле
𝑤-компоненты
скорости для одной из
реализаций.
Рис. 6. 𝑤-компонента скорости. Параметры турбулентности: 𝑞 = 1 м/с, 𝐿 = 200 м
5.
Численный метод
В настоящей работе использован стандартный численный метод конечных объемов вто-
рого порядка аппроксимации как по пространству, так и по времени (промышленный
пакет ANSYS CFX-12.1). Расчет с помощью LES производился для случая несжимаемой жидкости. В качестве модели турбулентности использована стандартная модель LES
Смагоринского с константой
6.
𝐶𝑠 = 0.1.
Постановка задачи о разрушении вихревого следа в турбулентной
атмосфере
Предлагается рассмотреть эволюцию следа за самолетом B-757 в турбулентной атмо-
сфере в области
определяются:
0 < 𝑥 < 700
м,
0 < 𝑦 < 300
м,
−100 < 𝑧 < 100
м. Размеры области
ТРУДЫ МФТИ. — 2014. — Том 6, № 3
9
И. С. Босняков, Г. Г. Судаков
ˆ
по оси
𝑥 — длиной волны Кроу 𝐿𝐶𝑟𝑜𝑤 = 300 м [2]. Необходимо иметь хотя бы 2 · 𝐿𝐶𝑟𝑜𝑤 ;
ˆ
по оси
𝑦
— расстоянием, которое вихревая пара преодолевает по вертикали за время
своей жизни;
ˆ
по оси
𝑧
— размахом вихрей.
Параметры газа берутся для воздуха постоянной плотности при
𝑇 = 25 ∘
C, уравнение
энергии не используется. На основе данных предварительного тестирования п. 3 было вы-
560 × 320 × 480 = 86,016 · 106 ячеек по осям
Δ𝑦 = Δ𝑧 = 0.6 м, Δ𝑥 = 1,2 м. В этом случае
брано следующее разбиение расчетной области:
соответственно
𝑥, 𝑦 , 𝑧 .
Размер ячеек составил
удлинение ячейки имеет величину 2, что удовлетворительно для решения задачи с помощью LES. Такое разбиение приемлемо также для решения на имеющемся суперкомпьютере
(потребная оперативная память для этой задачи составляет 110 Гб при использовании пакета программ ANSYS CFX). Однако при таком разбиении густота сетки в поперечном
сечении достаточно далека от потребной.
В качестве начального условия были использованы суперпозиция фоновой турбулентности и поле от пары вихрей противоположного знака (модель Проктора, п. 2).
Фоновая турбулентность в соответствии с п. 4 была сгенерирована при следующих параметрах:
𝑞 = 1
м/с,
𝐿 = 200
м, что характерно для режима посадки для высот более
100 м. При генерировании фоновой турбулентности все пространство как физическое, так
и волновое разбивалось равномерно на
140 × 60 × 40
интервалов. Таким образом, мини-
мальные размеры турбулентных вихрей атмосферы в данной реализации имеют порядок
5 м. Меньшие размеры турбулентных вихрей генерировать не удается из-за ограничений
по времени счета (рост времени вычислений
𝑇 ∼ 𝜆−6 ,
где
𝜆
— размер вихрей).
𝑥
Опишем граничные условия. На торцах расчётной области по оси
использовалось
циклическое граничное условие (условие периодичности). На всех остальных границах задавалось постоянное (по времени) распределение скорости, полученное в начальный момент
с помощью генератора турбулентности, см. п. 4. Такое граничное условие обеспечивает отсутствие затухания фоновой турбулентности в течение всего времени расчета эволюции
следа. Правомерность использования «замороженного» граничного условия обосновывается тем фактом, что характерное время изменения граничных условий для параметров
𝐿𝐴 = 200
м,
𝑞 = 0,1÷1
м/с имеет порядок
𝜏∼
𝐿𝐴
= 200÷2000 c,
𝑞
что значительно больше времени существования вихревого следа (время жизни следа
𝑇𝐿 ∼ 20÷70
с в зависимости от уровня турбулентности
𝑞 = 0.1÷1
м/с).
В соответствии с результатами п. 3 переменный шаг по времени был выбран из условия,
что максимальное число Куранта Сumax
7.
=2
и составлял
Δ𝑡 ∼ 0,05
с.
Результаты расчета
Результаты расчета для положения вихрей следа в последовательные моменты времени
представлены на рис. 7. Касание вихрей происходит около
𝑇𝑙𝑖𝑛𝑘 = 40
c, что совпадает
с предсказанием линейной теории Кроу–Бейта [2] и инженерной модели следа [9]. При
𝑡 = 40
c начинается распад вихрей за счет коротковолновых возмущений, что указывает
на наличие фазы быстрого распада вихрей (второй). Эта фаза характеризуется резким
возрастанием в окрестности вихревых нитей амплитуды пульсаций скорости, что хорошо
видно на рис. 7.
В соответствии с линейной теорией Кроу–Бейта [2] вихри раскачиваются преимущественно в плоскости, параллельной оси
𝑥
и наклоненной под углом
𝜃 = 48∘
к оси
𝑧.
Инженерная модель следа позволяет оценить размер и расположение области опасности.
Аэромеханика
10
ТРУДЫ МФТИ. — 2014. — Том 6, № 3
Рис. 7. Синусоидальная неустойчивость Кроу. Вид сверху. Уровень турбулентности — 𝑞 = 1 м/с,
масштаб турбулентности — 𝐿 = 200 м
Рис. 8. Область вероятного положения вихрей в следе за самолетом B-757 при 𝑡 = 16 с в соответствии с инженерной моделью [9] (сплошные линии) и рассчитанные с помощью LES. Вид на
вихревой след вдоль оси 𝑥. Уровень турбулентности — 𝑞 = 1 м/с, масштаб турбулентности —
𝐿 = 200 м
Рис. 9. Средняя циркуляция Γ5−15
На рис. 8 показано положение вихрей в соответствии с инженерной моделью и рассчитанное
с помощью LES.
Результаты расчета средней циркуляции в кольце с внутренним и внешним радиусами соответственно 5 и 15 м
Γ5−15
приведены на рис. 9. Видно, что LES дает величины
ТРУДЫ МФТИ. — 2014. — Том 6, № 3
И. С. Босняков, Г. Г. Судаков
11
циркуляции вполне удовлетворительно коррелирующие с инженерной моделью.
Рис. 10. Максимальная окружная скорости в вихре следа
Рис. 11. Картины 𝑥-компоненты завихренности в сечении 𝑥 = 200 м в последовательные моменты
времени 𝑡 = 8, 32, 40, 48 c
На рис. 10 приведены результаты расчета максимальной окружной скорости вихрей
следа в сравнении с результатами, полученными с помощью инженерной модели. Видно,
что максимальная окружная скорость в ядре вихря следа также не противоречит оценке
инженерной модели следа. Кроме того, видно, что при
𝑡 ∼ 10
с имеет место локальное (по
времени) уменьшение темпов снижения максимальной окружной скорости, что иногда наблюдалось в летном эксперименте [7], но не получило разумного объяснения. Объяснение
этого эффекта состоит в следующем. Представляется, что из-за диффузии вихря радиус
ядра растет, а максимальная скорость в ядре вихря уменьшается. Этому процессу препятствует процесс растяжения вихревых нитей вследствие синусоидальной неустойчивости, в
результате чего темп увеличения диаметра ядра замедляется. В какой-то момент времени
при достаточно сильно выраженной деформации вихревых нитей возникает локальное (по
времени) уменьшение размеров ядра вихря и возрастание максимальной скорости в ядре
вихря. Этот эффект инженерной моделью не описывается.
Как было указано во введении, процесс «потери» циркуляции определяется механизмом
дробления вихря следа на более мелкие вихри и выносом этих вихрей за пределы области
влияния вихревой пары. Это положение проиллюстрировано на рис. 11.
Таким образом, расчет в рамках LES описывает все характерные механизмы разрушения вихрей, но с учетом ограничений по имеющимся вычислительным ресурсам с разной
степенью точности. Наиболее аккуратно описывается процесс синусоидальной неустойчивости Кроу, так как он инициируется большими длинами волн. Поэтому время образования
Аэромеханика
12
ТРУДЫ МФТИ. — 2014. — Том 6, № 3
вихревых колец и зоны опасности по вихревому следу получаются с хорошей точностью.
Менее точно описывается процесс «потери» циркуляции и связанные с ним характеристики
следа (например, максимальная окружная скорость вихря).
Работа была выполнена при поддержке гранта РФФИ 13-08-00346 и при финансовой
поддержке Министерства образования и науки РФ в рамках договора № 700013728 от
21.11.2012 «Разработка моделирующего комплекса реалистичного восприятия оператором
(летчиком) сложных режимов полета и оценки его психофизиологического состояния» по
218 постановлению правительства РФ.
Литература
1. Bobylev A.V., Vyshinsky V.V., Soudakov G.G. and Yaroshevsky V.A.
Aircraft Vortex Wake
and Flight Safety Problems // Journal of Aircraft. — March–April 2010. — V. 47, N 2.
2. Crow S.C
. Stability Theory for a Pair of Trailing Vortices in a Turbulent Atmosphere //
AIAA Journal. — 1970. — V. 8, N 12. — P. 2172–2179.
3. Etkin B.
4. Holzдpfel F., Misaka T., Hennemann I.
Dynamics of Atmospheric Flight. — New York : John Wiley & Sons Inc., 1972.
Wake-Vortex Topology, Circulation, and Turbulent
Exchange Processes // AIAA Paper 2010–7992. — 2010.
5. Proctor F.H., Hamilton D.W., Han J.
Wake Vortex Transport and Decay in Ground Effect:
Vortex Linking with the Ground // AIAA Paper 2000–0757. — 2000.
6. Sarpkaya T.
Decay of Wake Vortices of Large Aircraft // AIAA Journal. — 1998. — V. 36,
N 9.
7. Shen S., Ding F., Han J., Lin Y.-L., Arya S.P., Proctor F.H.
Numerical Modeling Studies
of Wake Vortices: Real Case Simulation // AIAA Paper 0755. — 1999.
8. Thornber
B., Mosedale A., Drikakis D.
On
the
implicit
large
eddy
simulation
of
homogeneous decaying turbulence // J. Comput. Phys. — 2007. — V. 226. — P. 1902–1929.
(doi:10.1016/j.jcp.2007.06.030)
9. Вышинский В.В., Судаков Г.Г.
Вихревой след самолёта в турбулентной атмосфере //
Труды ЦАГИ. — 2005. — № 2667. — М. : Издательский отдел ЦАГИ.
10. Гарбарук А.В., Стрелец М.Х., Шур М.Л.
Моделирование турбулентности в расчетах
сложных течений. — СПб. : Изд-во Политехнического университета, 2012.
11. Ландау Л.Д., Лифшиц Е.М.
Теоретическая физика. Т. 6. Гидродинамика. — 1988. —
M. : Наука.
Поступила в редакцию 14.11.2013.
Download