УДК 621.336
О. А. СИДОРОВ, А. Н. СМЕРДИН, А. Е. ЧЕПУРКО (ОмГУПС)
Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Омский государственный университет путей сообщения», пр. Маркса, 35, Омск, Россия, 644046, тел 8-951-413-44-35, эл. почта: [email protected]
ПРИМЕНЕНИЯ СРБ-МЕТОДА ДЛЯ ОПРЕДЕЛЕНИЯ АЭРОДИНАМИЧЕСКИХ СИЛ ТОКОПРИЕМНИКА И ЕГО АЭРОДИНАМИЧЕСКОГО УСТРОЙСТВА
На сегодняшний день одним из самых перспективных направлений в моделировании воздействия внешней среды на объекты является вычислительная аэрогидродинамика (Computational Fluid Dynamics - CFD). Сегодня CFD-исследования могут быть эффективно применены на железнодорожном транспорте для решения задач моделирования обтекания подвижного состава воздушной и для расчета аэродинамических характеристик токоприемника и его аэродинамических устройств. Процесс моделирования CFD состоит из: генерации расчётной сетки; нахождения оптимального числового алгоритма;
моделирования турбулентных течений. Особенностью обтекания воздушным потоком токоприемника является то, что он представляет собой конструкцию, состоящую из большого количества элементов. При обдуве
спектры обтекания данных элементов могут пересекаться, т.е. их взаимное расположение обуславливает и действующие на них аэродинамические силы. Кроме того, расчетная сетка при СРВ-исследовании всего токоприемника крупнее, чем при СРВ-исследовании отдельных деталей или узлов, что сказывается на точности расчета (рис. 1).
Таким образом, анализируя конструкцию токоприемника (рис. 1), авторы статьи пришли к выводу, что для его СРВ-исследования целесообразно применять «поэлементный СРВ-метод». Для этого необходимо выявить элементы и узлы токоприемника с минимальной степенью влияния друг на друга при обтекании воздушным потоком, затем рассчитать аэродинамические силы, действующие на эти элементы и узлы. Далее, применяя принцип суперпозиции, определить аэродинамические силы, действующие на весь токоприемник.
Рис. 1. Расчетная сетка при С Р Э - исс лсдо ва н и и: а - при исследовании токоприемника, б - при исследовании верхнего узла токоприемника
Аэродинамическое воздействие на кон- дифференциальных уравнений Навье-Стокса в струкции и устройства описывается системой частных производных:
© Сидоров О. А. и др., 2013 ISSN2307-4221 Електрифтащя транспорту, № 6. - 2013.
Г//. dp „ 8 P—= Jx ■ P—— + 2--
dt x dx dx
м-
1JL
3 dx
м-
( du du du ^ —- +—- + —-
dx dy dz
du dp „ ô
P—- = iy ■ P— +2--
dt dy dy
M-
M-
_2_d_ 3 dy'
du
( du du du \
-JL +-+-—
dx dy dz
dp
dz
M-
2 d_
3 dz
M-
' du x duy +duz^ y dx dy dz y
dx
d + —
dx ( duуЛ
. дУ , +—
dy
du _2
dz d_
ç-
ç-
Ç-
d_
dy
K dx
d_
dz
'<K
dx
d_
dx
f(K
dx
M-
du„
duy\
dy dx
d_
dz
M-
dUX
dz
du
_Z_
dx
M-
du du ^ —- + —-
dy dz
'duy | du y dz dy
du du, Л
JL
dx
M-
f du du "ï
_У_ _|__x_
dx dy
M-
dy du
_Z
dx du„
dz _ dz
(1)
A
dy
M-
du
du ^ _y
dy dz
du ^ — + —2
dy dz
dp , d(pux) d(puy) d(puz)
dt
= 0,
dx dy dz
где и - вектор скорости потока с проекция-
ми il.
иу, II_ на соответствующие оси коор-
динат [м/с];
р - давление среды [Па]; р - плотность среды [кг/м3]; р - динамическая вязкость [Пас] при заданной температуре;
/х. Д.. /_ - проекции вектора объемной силы на соответствующие координатные оси [Н]; ^ - объемная вязкость [Пас]. Если принять, что // = const и С, = const, то система уравнений (1) значительно упрощается:
dur = Jx 1 dp + v d\ d\ d\
dt P dx dx2 dy2 dz2
duy = jy 1 dp + v ~d\ d\ 4- 7 d2ux 1 >
dt P dy dx1 dy2 dz2
du. 1 dp ~d2u. d2U, d2u.
~dt~ = Jz P ~z + V dx2 --f " dz2
+
+
+
fL.
KP
i£-\p
V b
V
3
\
V
3
du„ du, du.
dx
dx
du,,
dy
du h—-
dy
du,.
dz
1 ¿4 dz
du.
(2)
dx dy dz
dp | d(pux) | d(puy) ^ d(puz)
dt + + +
dx
dy
dz
2/
= 0,
где V - кинематическая вязкость [м /с], которая является отношением динамической вязкости к плотности воздуха.
Найти прямое решение данной системы уравнений на сегодняшний день не представляется возможным из-за того, что оно зависит от начальных и граничных условий. Поэтому при решении практических задач применяют численное моделирование, основной проблемой которого является моделирование турбулентных течений, для чего используют следующие СРБ-модели турбулентности [1]:
1) Прямое численное моделирование (DNS -Direct Numerical Simulation). Отсутствуют дополнительные уравнения. Нестационарные уравнения Навье-Стокса решаются с очень мелким шагом по времени, на мелкой пространственной сетке.
2) Осреднение по Рейнольдсу уравнений Навье-Стокса (RANS - Reynolds Averaged Navier-Stokes equations), замкнутых с помощью алгебраических или дифференциальных моделей турбулентности.
© Сидоров О. А. и др., 2013
3) Метод крупных вихрей (LES - Large Eddy Simulation). Влияние вихрей, которые по размерам меньше, чем размеры ячейки расчётной сетки, заменяется эмпирическими моделями.
4) Метод DES (Detached Eddy Simulation), который является комбинацией подходов RANS и LES. В зоне внешнего течения используется RANS, а в зоне отрыва потока с крупными вихрями - LES.
Наибольшая точность расчетов может быть получена при прямом численном моделировании (DNS). Однако при сегодняшнем уровне развития вычислительной техники отсутствует возможность решения задач с реальной геометрией с Колмогоровским масштабом длины и Колмогоровским масштабом времени при больших числах Рейнольдса [1]. Современные суперкомпьютеры способны решать задачи с
применением метода крупных вихрей (LES) и метода моделирования отдельных вихрей (DES), однако такая техника имеется только в крупных научно-исследовательских центрах. Поэтому для решения широкого спектра инженерных задач аэрогидродинамики широко используется метод осреднения по Рейнольдсу системы дифференциальных уравнений Навье-Стокса в частных производных (RANS).
Дополнительный модуль инженерного анализа Flow Simulation программного комплекса САПР SolidWorks 2010 применялся авторами для CFD-исследований. Для моделирования турбулентных течений в нем используется стандартная k-s модель турбулентности Лаун-дера-Сполдинга, основанная на решении системы дифференциальных уравнений в частных производных:
(
Р
ы
dx
+2-
Р
дх
ы
Mef
-ы-
д[их]
V & У
дЫ,т 15Ы1_ дМ
ду
д_
ду
К]
dz
дх
д[иу] д[иу]
о 5 +2--
ду
(
Р
ы
V
о д
+2—
dz
дх Я/
дх
Ы-
ф',]
ду
ду
д_
дх
ы
ду д\_ь
дх
у d f
dz Mef
У_ L V
д[их] d[uz]
dz
дх
dz
Я/
+
г и д\и 1 г ,
Ы-тг+Ы
N
S[uz]
V & У
ду
д_
дх
д[их] ду
dz
fö[ux\
dz
д[р]
dy /
£kP
dx d[p]
d_
dz
/V
дЫ ,
dz
dy
(3)
dz
S[uz]
dx
y + 8 f
y_ dy Vef V
дЫ, ды
dz
dy
где р - плотность среды, кг/м3;
\р\ - осредненное давление, Па;
[» ]. [»_ ] - осредненные проекции
скорости на координатные оси, м/с;
¡ле£- эффективная динамическая вязкость,
Пас.
Эффективная динамическая вязкость определяется по формуле [1]:
^=р\у + у}), (4)
где V - кинематическая вязкость, м2/с; V , - кинематическая вихревая вязкость, м2/с.
Кинематическая вихревая вязкость рассчитывается с помощью зависимости Прандтля-Колмогорова через величины кинетической энергии турбулентности к и диссипации этой энергии 8 :
к2
Vj =Cv-->
S
(5)
где Су= 0,09 - эмпирический коэффициент,
выбранный на основании экспериментов с близкими к равновесным турбулентными потоками [2].
Система дифференциальных уравнений, связывающих кинетическую энергию турбулентности к и энергию её диссипации 8 , записывается в виде:
© Сидоров О. А. и др., 2013
V
д([их\к) д{[иу}к) а([и,]-*)__а _
дх
(V + Vj) (дк_ дх
| | a(K]-g) а
дх
дх
JL
ду
ду
[ду
dz
дх
+Cj • — • — - С2 • — + С3 к р к
д_ dz
G2
(v + vjjöe &
к- р
где С - скорость генерации турбулентности; ак= 1,0; а = 1,3; Сг = 1,43; С2 = 1,92;
С3 = 0 - эмпирические коэффициенты, определенные на основе решения задач о плоской струе и слое смешения [2].
С помощью САПР 8оПс1\¥огк8 2010 авторами статьи было произведено СРБ-исследование токоприемника АИСТ и установлены компенсирующие аэродинамические силы, необходимые для получения рациональной аэродинами-
ческой характеристики асимметричного токоприемника [3, 4].
Аэродинамические силы, которые должно создавать аэродинамическое компенсирующее устройство для токоприемника, движущегося коленом вперед, во всем диапазоне рабочих высот и для скорости движения от 33,3 м/с до 69,4 м/с приведены в табл. 1.
Для токоприемника, движущегося коленом назад - в табл. 2.
Таблица 1
Аэродинамические компенсирующие силы необходимые при движении токоприемника коленом вперед
и, м/с 33,3 38,8 44,4 50 55,5 61,1 66,7 69,4
h, мм Рт у стр. ? H
400 1,8 2,8 3,8 5,5 7,1 9,6 11,7 13,5
1000 1,5 2,1 2,9 4,4 5,8 8,0 9,6 11,2
1200 1,4 1,9 2,5 3,9 5,3 7,5 9,0 10,7
1400 1Д 1,4 2,0 3,2 4,3 5,9 7,5 8,8
1600 0,6 1,2 1,5 2,6 3,4 4,8 6,3 7,6
1900 0,2 0,7 1,0 1,9 2,2 3,5 4,8 6,1
Таблица 2
Аэродинамические компенсирующие силы необходимые при движении токоприемника коленом назад
и, м/с 33,3 38,8 44,4 50 55,5 61,1 66,7 69,4
h, мм Рт устр., H
400 -4,6 -6,2 -8,1 -10,2 -13,6 -15,3 -16,5 -16,1
1000 -6,0 -7,7 -9,7 -12,0 -15,2 -17,2 -18,5 -18,3
1200 -6,4 -8,0 -10,1 -12,5 -15,8 -18,1 -19,5 -19,6
1400 -5,5 -7,3 -9,6 -12,0 -15,5 -18,0 -19,5 -19,8
1600 -7,0 -8,8 -11,4 -14,0 -17,6 -20,3 -21,9 -22,4
1900 -7,2 -9,1 -11,9 -14,6 -18,3 -21,0 -22,8 -23,6
© Сидоров О. А. и др., 2013 116 ISSN2307-4221 Електрифтащя транспорту, № 6. - 2013.
Полученные данные показывают, что модули компенсирующих аэродинамических сил различаются при движении в различных направлениях. Поэтому для компенсации аэродинамической подъемной силы асимметричного токоприемника целесообразно применять крылья определенного профиля.
Авторами статьи были проведены исследования аэродинамических характеристик крыльев, длина которых принималась равной 180 мм.
Аэродинамическое крыло устанавливалось в соответствии с рисунком 2. Такое положение было выбрано в соответствии с данными круговой обдувки [5]. На рисунке 2 приведена кинематическая схема, на которой изображен асимметричный токоприемник электроподвижного состава, состоящий из основания 1, пневмопривода 2, взаимодействующего с шарнирно сочлененной нижней штангой 3 и верхней штангой 4, на которой размещены каретки 7 с полозами 8, причем верхняя штанга 4 связана соединительной тягой 5 с основанием 1, направляющей штанги параллелограммного механизма 6 кареток 7, аэродинамического устройства 9. При этом аэродинамическое устройство 9 включает крыло 10, шарнирно соединенное с каретками 7 и способное изменять угол атаки за счет системы управляющих рычагов 11, приводимых в действие встречным потоком воздуха с помощью экранов 12. Упру-
гий элемент 13 обеспечивает заданный угол атаки в зависимости от скорости движения и осуществляет возврат механизма в первоначальное состояние при остановке поезда.
13 11
Рис. 2. Кинематическая схема токоприемника с установленным крылом
Спектры обтекания крыла профиля ЫА8А-2210, полученные в результате СРВ-исследования для скорости 69,4 м/с приведены на рис. 3.
Рис. 3. Спектры обтекания крыла профиля ЫАСА-2210 для угла атаки: а - 0°, б - 10°, в - 20°
© Сидоров О. А. и др., 2013
а
б
в
Спектры обтекания на рисунке 3, показывают, что для углов атаки 0° - 10° течение имеет ламинарный характер, без образования турбулентных вихрей. Профиль данного типа согласно [6] имеет критический угол атаки равный 16°, после чего коэффициент подъемной силы начинает резко снижаться, а лобовое сопротивление расти. При угле атаки в 20°, как видно из спектра обтекания, на задней кромке крыла образуется вихрь, течение приобретает турбулентный характер. Также с ростом угла атаки выше критического увеличивается лобовое сопротивление с одновременным снижением качества крыла.
Следовательно, регулирование аэродинамической подъемной силы токоприемника следует осуществлять при угле атаки крыла аэродинамического устройства ниже критического.
Авторами статьи было проведено СРБ-исследование крыльев с профилями заданного типа. Результаты этого исследования для скорости 69,4 м/с приведены на рисунках 4-6. Угол атаки изменяли от 0 до 45 градусов, продувка осуществлялась в обоих направлениях. Предполагается, что при движении асимметричного токоприемника коленом назад крыло аэродинамического устройства будет создавать опускающую силу, при движении в обратном направлении - подъемную силу.
20 H 10 5 0 -5 -10 -15 -20 -25 -30
„
» г с
! Ж
ж
f
X
с
X X
( х с
¡«движение коленом вперед ^движение коленом назад
а б
Рис. 4. Результаты СР Э -исслсдовани я крыла: а - профиль ЦАГИ-6-12%; б - зависимость аэродинамической подъемной силы крыла от угла атаки
15
н
5 0 -5 -10
-15
Рбт у-ф _2о -25 -30
Ж )
4 к ' Ж ) ' ж ; : *
ж
(
X
4 V
X
г X > С X 5 с х
Ж движение коленом вперед Хдвнженне коленом назад
а б
Рис. 5. Результаты СРБ-исследования крыла: а - профиль ЫАСА-2210: б - зависимость аэродинамической подъемной силы крыла от угла атаки
20 н 10 5 0 -5 -10 -15 -20 -25 -30
-Ж-it
ж—*—::
X-it—X"
Ждвижение коленом вперед X движение коленом назад
а б
Рис. 6. Результаты С Р Э - исс лсдо ва н и я крыла: а - профиль и8А-27; б - зависимость аэродинамической подъемной силы крыла от угла атаки
© Сидоров О. А. и др., 2013
Эффективное регулирование контактного нажатия может быть обеспечено при создании крылом аэродинамических подъемных сил, приведенных в таблице 1 и таблице 2 и осред-ненных по высоте. В наибольшей степени этому требованию соответствует крыло с профилем NACA-2210.
Выводы
1) Способ моделирования газодинамических процессов посредством ПЭВМ с использованием CFD-исследования с применением «поэлементного CFD-метода», является целесообраз-
СПИСОК ИСПОЛЬЗОВАННЫХ источников
1. Выбор подхода к определению турбулентной скорости распространения пламени в цилиндре газового двигателя [Текст] / Абрамчук Ф. И. и др. // Вестник национального технического университета «ХПИ». - 2009. - № 47. - С. 112 - 126.
2. SolidWorks 2007/2008. Компьютерное моделирование в инженерной практике [Текст] / А. А. Алямовский и др.// СПб.: БХВ-Петербург, 2008. -1040 с.
3. Чепурко А. Е. Регулирование аэродинамической подъемной силы токоприемника в процессе его взаимодействия с контактной подвеской [Текст] / А. Е Чепурко // Известия Транссиба. - 2012. - № 4. -С. 87 - 93
4. Сидоров О. А., Смердин А. Н., Емельянов М. В. Адаптация математической модели универсального измерительного токоприемника для исследования системы токосъема на линии Москва - Санкт-Петербург [Текст]/ О. А. Сидоров, А. Н. Смердин, М. В. Емельянов // Известия Транссиба. - 2012. -№4.-С. 51-58.
5. Радченко И. П. Круговая обдувка профиля NACA 23012 в аэродинамической трубе Т-103Н ЦАГИ. Технические отчеты ЦАГИ № 161, 1959.
6. Атлас аэродинамических профилей крыльев [Текст] / Ушаков Б.А. и др.//. Издание БИТ НКАП при ЦАГИ.-1940.
Поступила в печать 17.04.2013.
Ключевые слова: CFD-метод, модели турбулентности, аэродинамические силы, профиль крыла, система дифференциальных уравнений Навье-Стокса, спектры обтекания.
Внутренний рецензент Кузнецов В. Г.
ным при выполнении аэродинамических расчетов токоприемника.
2) Рассчитанные аэродинамические силы, создаваемые аэродинамическим устройством, позволили разработать устройство для автоматического регулирования нажатия.
3) Произведено CFD-исследование крыльев заданного профиля, на основании чего был сделан вывод, что для регулирования аэродинамической подъемной силы токоприемника необходимо использовать крыло с профилем NACA-2210.
REFERENCES
1. Abramchuk F.I. Vybor podkhoda k opredeleniyu turbulentnoy skorosti rasprostraneniya plameni v tsilin-dre gazovogo dvigatelya [The approach to the definition of the turbulent flame propagation velocity in the cylinder gas engine], Bulletin of the National Technical University "KPI", 2009, no. 47, pp. 112-126.
2. Alyamovskiy A. A. SolidWorks 2007/2008. Komp'yuternoe modelirovanie v inzhenernoy praktike [SolidWorks 2007/2008. Computer modeling in engineering practice], SPb., BKhV-Peterburg Publ., 2008, 1040 p.
3. Chepurko A. E. Regulirovanie aerodinamich-eskoy pod"emnoy sily tokopriemnika v protsesse ego vzaimodeystviya s kontaktnoy podveskoy [Regulation of aerodynamic lift susceptor during its interaction with catenary], Proceedings of the Trans-Siberian, 2012, no. 4, pp. 87-93.
4. Sidorov O. A., Smerdin A. N., Emel'janov M. V. Adaptacija matematicheskoj modeli universal'nogo iz-meritel'nogo tokopriemnika dlja issledovanija sistemy tokosema na linii Moskva - Sankt-Peterburg [Adapting mathematical model multimeter pantograph current collection system for studies on the Moscow - St. Petersburg], Proceedings of the Trans-Siberian, 2012, no. 4, pp. 51-58.
5. Radchenko I. P. Kmgovaja obduvka profilja NACA 23012 v ajerodinamicheskoj trabe [Circular profile NACA 23012 blowing in the wind tunnel], Technical Reports, 1959.
6. Ushakov B.A. Atlas ajerodinamicheskih profilej kryl'ev [Atlas of aerodynamic airfoils]. Izdanie BIT NKAP pri CAGI, 1940.
Внешний рецензент Саенко Ю. Л.
Контактное нажатие при взаимодействии токоприемника с контактной подвеской не остается постоянным. Одной из его дестабилизирующих составляющих является аэродинамическая подъемная сила токоприемника. Поэтому для обеспечения надежного и экономичного токосъема необходимо рассчитать рациональную аэродинамическую силу токоприемника, при которой износ контактирующих элементов будет минимальным с сохранением надежности контакта.
Проведение линейных испытаний токоприемника требует больших материальных затрат, как и его обдувка в аэродинамической трубе. Однако в настоящее время за счет широкой компьютеризации находят применение методы вычислительной гидрогазодинамики (СРР). На железнодорожном транспорте данные методы могут быть применены для расчета лобового сопротивления локомотива, для моделирования воздействия ветра на контактную подвеску и др. Авторами СРР-метод применялся для расчета аэродинамических характеристик токоприемника.
© Сидоров О. А. и др., 2013
Объектом исследования является метод применения CFD-анализа для расчета аэродинамических сил, действующих на токоприемник при высокоскоростном движении.
В статье авторами проведено исследование особенностей применения CFD-метода для моделирования аэродинамического воздействия на токоприемник и выбран профиль крыла для его последующего применения в аэродинамическом устройстве.
УДК 621.336
О. О. СИДОРОВ, А. М. СМЕРДШ, О. G. ЧЕПУРКО (ОмГУПС)
Федеральна державна бюджетна осв1тня установа вищоТ професшноТ осв1ти «Омський державний ушвер-ситет шлях1в сполучення», пр. Маркса, 35, Омськ, Роая, 644046, тел.: 8(951) 413-44-35, e-mail: [email protected]
ЗАСТОСУВАННЯ CFD-МЕТОДУ ДЛЯ ВИЗНАЧЕННЯ АЕРОДИНАМ1ЧНИХ СИЛ СТРУМОПРИЙМАЧА I ЙОГО АЕРОДИНАМ1ЧН1 ПРИСТР01
Контактне натискання при взаемоди струмоприймача з контактною пщвюкою не залишаеться постш-ним. Одшею з його дестабш1зуючих складових е аеродинамнна тдйомна сила струмоприймача. Тому для забезпечення надшного i економнного струмозшмання необхщно розрахувати рацюнальну аеродинамнну силу струмоприймача, при якш зношування контактуючих елементш буде мш1мальним ¡з збереженням на-дшносп контакту.
Проведения лшшних випробувань струмоприймача вимагае великих матер1альних витрат, як i його обдувания в аеродинамннш Tpy6i. Проте в даний час за рахунок широко! комп'ютеризаци знаходять застосу-вання методи обчислювальноТ пдрогазодинамки (CFD). На зал1зничному транспорт даш методи можуть бути застосоваш для розрахунку лобового опору локомотива, для моделювання впливу BiTpy на контактну пщвюку та ¡н. Авторами CFD-метод застосовувався для розрахунку аеродинамнних характеристик струмоприймача.
Об'ектом дослщження е метод застосування СРР-анал1зу для розрахунку аеродинамнних сил, що д1ють на струмоприймач при високошвидюсному pyci.
У CTarri авторами проведено дослщження особливостей застосування CFD-методу для моделювання ае-родинамтного впливу на пантограф i вибраний профшь крила для його подалыиого застосування в аеро-динамнному пристроТ.
Ключов1 слова: CFD-метод, модел1 турбулентносп, аеродинамтш сили, профшь крила, система дифере-
нц1альних ртнянь Нав'е-Стокса, спектри обт1кання.
Внутршнш рецензент Кузнецов В. Г. Зовншшш рецензент Саенко Ю. Л.
UDC 621.336
О. A. SIDOROV, А. N. SMERDIN, А. Е. CHEPURKO (OSTU)
Federal State Educational Institution of Higher Professional Education "Omsk State Transport University" 35 Marks Ave, Omsk, Russia, 644046, tel.: 8(951) 413-44-35, e-mail: [email protected]
USE CFD-METHOD TO DETERMINE THE AERODYNAMIC FORCES PANTOGRAPH AND AERODYNAMIC DEVICES
Contact pressure in the interaction with the pantograph catenary is not constant. One of its components is destabilizing aerodynamic lift trolley. Therefore, to ensure reliable and efficient current collection is necessary to calculate the rational aerodynamic force collector, in which the wear items will be minimal contact with reliable contact.
Conducting tests of linear pantograph requires high material costs, as well as its blowing in the wind tunnel. Currently, however, due to the wide computerization methods are used computational fluid dynamics (CFD). In rail transport, these methods can be used to calculate the drag of the locomotive, to simulate the effects of wind on the catenary, etc. The authors of CFD-method was used to calculate the aerodynamic characteristics of the pantograph.
The object of study is the method of application of CFD-analysis to calculate the aerodynamic forces acting on a current collector for high-speed motion.
In this paper the authors study the features of CFD-modeling method for aerodynamic effects on the current collector and selected airfoil for its use in the aerodynamic device.
Keywords: CFD-method, the turbulence model, the aerodynamic forces, the wing profile, the system of the Navier-Stokes differential equations, the spectra of the flow.
Internal reviewer Kuznetsov V. G. External reviewer Saenko U. L.
© Сидоров О. А. и др., 2013