CC BY
25
УДК 532.517.4 DOI: 10.14529/mmp190103
ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ВНУТРИКАМЕРНЫХ НЕСТАЦИОНАРНЫХ ТУРБУЛЕНТНЫХ ТЕЧЕНИЙ. ЧАСТЬ 1
А.М. Липанов1, С.Ю. Дадикина2, А.А. Шумихин2, М.Р. Королева2, А.И. Карпов2
1Институт прикладной математики им. М.В. Келдыша РАН, г. Москва, Российская Федерация
2Удмуртский федеральный исследовательский центр УрО РАН, г. Ижевск, Российская Федерация
В работе представлена методика моделирования трехмерных внутренних нестационарных турбулентных течений, в частности, течения продуктов сгорания в твердотопливном ракетном двигателе. Приведена система определяющих уравнений, описывающая поток сжимаемого вязкого газа, записанная в цилиндрической системе координат. Предложен вычислительный алгоритм, разработанный на основе модифицированной схемы расщепления векторов потоков, относящийся к классу методов, использующих подход Годунова. Данный алгоритм пригоден для сквозного расчета внутреннего течения по всему тракту ракетного двигателя, включающего как зоны дозвукового течения в камере, так и зону сверхзвукового течения в сопле. Полученные результаты численного моделирования внутреннего турбулентного течения газа в модельном ракетном двигателе показывают осциллирующий ударно-волновой характер процессов, протекающих в камере двигателя на начальном этапе работы. Определено время выхода двигателя на стационарный режим.
Ключевые слова: внутрикамерные процессы; турбулентность; нестационарное течение; вычислительная гидрогазодинамика.
Введение
Для исследования газодинамических процессов, протекающих в различных технических устройствах и установках, в том числе в твердотопливных ракетных двигателях (РДТТ), важное значение имеет численное моделирование трехмерных внутренних турбулентных потоков. Трудности, связанные с постановкой физического эксперимента, выдвигают на ведущий план развитие методов математического моделирования таких течений. Улучшение рабочих характеристик РДТТ, снижение материальных затрат при разработке новых двигателей, сокращение сроков их отработки в значительной степени зависят от корректности численных исследований процессов, протекающих в двигателе, от реализуемой точности проводимых расчетов. Моделирование процесса выхода РДТТ на стационарный режим работы представляет задачу высокой сложности. Сложность моделирования течения продуктов сгорания в РДТТ обусловлена турбулентным характером течения, нестационарностью протекающих в тракте двигателя процессов, влиянием градиентов давления [1-6]. Значительный интерес также представляет и изучение влияния вдува продуктов сгорания в камеру двигателя на параметры турбулентного потока в каналах с геометрически сложной формой. Разработка новых методик для численного решения задач гидрогазодинамики всегда является актуальной задачей. С одной стороны, это объясняется необходимостью совершенствования известных, хорошо зарекомендовавших себя методов, в плане более эффективного использования ими увеличивающихся вычислительных
возможностей компьютеров, с другой стороны, программной реализацией новых, обусловленных современными практическими потребностями более сложных постановок решаемых задач.
Конфигурации внутренних поверхностей реальных двигателей имеют достаточно сложную геометрию. Течение продуктов сгорания в РДТТ имеет существенно дозвуковые области, области высоких дозвуковых скоростей и область соплового сверхзвукового течения. Проведение численного моделирования потока газа во всех этих областях в рамках единого алгоритма является сложной задачей, но успехи современной вычислительной гидрогазодинамики и растущая вычислительная мощность компьютерной техники побуждают к разработке математических моделей и алгоритмов, обеспечивающих сквозной расчет внутренних течений в ракетных двигателях.
В данной работе для моделирования турбулентного течения в тракте модельного твердотопливного двигателя предложен вычислительный алгоритм, основанный на численном интегрировании уравнений сохранения, записанных в цилиндрической системе координат [5, 7-10]. Во всех существенно нестационарных течениях возникают пульсации параметров потока, что накладывает определенные требования к используемой методике. Предложенный алгоритм пригоден для сквозного расчета всего тракта двигателя, как дозвукового течения в камере двигателя, так и сверхзвукового в сопле, при наличии пульсаций давления без введения искусственной вязкости. Алгоритм основан на модифицированном методе расщепления векторов потоков, относящимся к классу методов, использующих подход Годунова [11-14]. Для обеспечения достоверности получаемых результатов при создании методик численного исследования, предназначенных для моделирования течений в зарядах реальных конфигураций, проводится проверка корректности математической модели и тестирование вычислительного алгоритма на модельных задачах. В работе приведены результаты исследований внутрикамерного нестационарного течения сжимаемого вязкого газа в модельном ракетном двигателе.
1. Постановка задачи
Для моделирования трехмерного течения продуктов сгорания в камере РДТТ использовалась система уравнений гидромеханики в цилиндрической системе координат. В векторном виде система записывается следующим образом
dQ т
dt dx
x dFr 1 dFe ^ 1
+ ^ + + Fee = -
dr r dO
д т д т д T T-rLx + — rLr + — Lq дх dr до
+ Le
(1)
где Q - вектор гидромеханических параметров (ГМП); Fx, Fr, Fe, Fee - векторы конвективных потоков; Lx, Lr, Le, Lee - векторы диффузионных потоков; х, r, O -продольная, радиальная и угловая координаты; t - время. Векторы ГМП, конвективных и диффузионных потоков имеют вид
Q=
р ри pv pw v
pu pu2 + p puv puw F -P , ree—- r uv
pv ) Fx puv , Fr= pv2 + p , Fe= pvw v2 — w2
pw puw pvw pw2 + p 2vw
pe puh pvh pwh vh
(2)
r
Lx
0
Txx Txr Гхв
UTxx +VTxr +WTx0
Lr
0
Trx Trr
Tre
UTrx +VTrr+WTre
(3)
0
Гвх
Le = Tgr
Tee
UTQx+VTQr +WTee_
Здесь u, v, w - компоненты вектора скорости, p - давление, p - плотность, e - удельная энергия, h - удельная энтальпия, Txx, Trx, Tex, Txr, Trr, Ter, Txe, Tre, Tee - компоненты тензора вязких напряжений.
Компоненты тензора вязких напряжений вычисляются так
0
1 0
, L<ee—- — Tee
r
Ter
0
du 2
тхх = 2ß----ßchvil,
dx 3
9 9V 2 rV О Trr = ¿11—---ß div IL ,
dr 3
/1 dw v\ 2 твв = ¿ß —H— — -ßchvil,
Tr
(dv du ^ \dx dr
r qq r
3
Txe — Tex — ß
Tre — Ter — ß
1 du dw
r dQ dx
1 dv dw
r dQ dr
(4)
где ^ - коэффициент динамической вязкости, П = (и; V; и) - вектор скорости. Дивергенция вектора скорости в цилиндрической системе координат определяется формулой
^ du dv 1 dw v dx dr r dQ r
(5)
Аналогичная система уравнений использовалась в работах [5, 7-10]. Температура и полная удельная энтальпия вычислялись по соотношениям
1
T
U2 + V2 + w2
e —
h
u2 + v2 + w2
Cv У 2 J 2
Уравнение состояния записывалось в следующем виде
Р- !м e -
+ Cp T.
p
u2 + v2 + w2
2
(6)
(7)
Здесь С - удельная теплоемкость газа при постоянном объеме, Ср - удельная теплоемкость газа при постоянном давлении, Т - температура, Я - удельная газовая постоянная, V - показатель адиабаты.
2. Вычислительный алгоритм
Для расчета векторов конвективных потоков Ях, Яг, Яе (2) использовалась модифицированная схема Стигера - Уорминга (схема расщепления векторов потоков),
предложенная в работах [11, 12]. Применяя уравнение состояния (7), из выражений для векторов (2) можно исключить давление [11]. С использованием взвешивания по Фавру / = р//р, вектор Рх будет записан
ри
(V - 3) и2 + (V - 1) V2 + (V - 1) и)2
р\(у - 1) е -
рщ)
рии рин
Вектор Рх можно представить следующим образом
Рх ЗхЯ-,
(8)
(9)
где Зх
сЖ
дЯ
матрица Якоби. Используя (8), матрица Зх определяется так [15]
Зх
V - 3
и2 +
0
и-1
2
-III)
-и) и) 0
1 0 0 0
- 3) и) -(V -1)V -(V-1) и) V -1
') и) 0 0
и) 0 и) 0
н 0 0 0
(10)
Собственные числа матрицы Зх равны [11]
Ах1 = Ах2 = Ах3 = и ) Ах4 = и + с
А.
х5 и с
(11)
где с = (^/р)1/2 - скорость звука. Вектор Рх может быть записан в следующем виде [11]
Рх ЗхЯ СхЛхС- Я ,
(12)
где Сх - матрица правых собственных векторов Якобиана Зх; С-1 - матрица левых собственных векторов. Собственные числа Якобиана (11) являются диагональными элементами матрицы Лх
Лх
и) 0 0 0 0
0 и) 0 0 0
0 0 и) 0 0
0 0 0 и + с 0
0 0 0 0 ис -
(13)
Тогда вектор Рх, с использованием собственных чисел может быть представлен [11]
2и
2^ - 1)А 1 + Ах4 + Ах5 2(v - 1)Ах1и + Ах4 (и + с) + Ах5 (и - с)
2(V - 1)Ах1) + Ах4'С + Ах5'С
2(v - 1)Ах1и) + Ах4и) + Ах5и) 2(v - 1)Ах1Н + Ах4с + Ах5с
(14)
2
2
2
2
Основой схемы Стигера - Уорминга, как и других алгоритмов, использующих подход Годунова, является расщепление вектора потоков (14) на две части
F = F + + F"
1 x 1 x + 1 x
(15)
где вектор Ях+ соотносится с вектором, сформированным из положительных собственных чисел матрицы ,]х. Соответственно вектор Ях соотносится с вектором, сформированным из отрицательных собственных чисел. Матрицу Лх можно расщепить аналогично
Лx = Л+ + Л- .
Тогда векторы F+, F", входящие в выражение (15), примут вид
F+ = (G+) Л+ (G+)"1 Q+, F" = (G") Л- (G")"1 Q
(16)
(17)
Все собственные числа Л |, Л- , являющиеся элементами матриц Л |, Лх , определяются следующим образом [11]
Лx — Л+ + Л-
Л+
Ax + 1 Ax 1
Л-
Ax — |Ax 1
(18)
Используя (14) - (18), векторы можно записать [11]
F+
х
JL
2и
2(v — 1)Л+1 + Л+4 + Л+5
- 5+)
р_
2и
2(v - 1)Л+1и+ + Л+4 (U+ + с+) + Л+5 (U+ 2(v - 1)Л+г+ + Л+4г+ + Л+5г+ 2(v — 1)Л+1г5+ + Л+4г5+ + Л+5{й+ 2(v - 1M+h+ + л+4h+ + л+5h+
2(v - 1)Л-1 + Л-4 + Л-5
2(v - 1)Л-1и" + Л-4 (U" + ) + Л-5 (U" - ) 2(v - 1)Л-1г" + Л-4г" + Л-5г"
2(v - 1)Л"1г5" + Л—W" + Л—W"
2(v - 1)Л"^" + Лх4^'~ + л"5h"
(19)
(20)
Для численного интегрирования системы основных уравнений (1) по пространству применялся метод конечных объемов. Для расчета векторов потоков Я—1/2 и Я+1/2 соответственно на левой и правой гранях /-го контрольного объема, необходимо определить векторы параметров Q^х 1/2 и Q^|1/2. Параметры на границах /-го контрольного объема Q~х1/2 и Qj+1/2 принимались равными в центре ячейки. На рис. 1 показан /¿-ый контрольный объем и параметры потока газа на его границах.
2
2
1/2 Л Qi+m
-•-•-
Рис. 1. Параметры потока газа Q¿X1/2 и Q|1/2 на гранях /¿-ого контрольного объема
Далее используя выражения (15), (19), (20), вычисляем значения потока Ях. Аналогичные процедуры производятся для векторов Яг и Яе. Векторы потоков Ьх, Ьг, Ье
определяются с использованием центрально-разностной схемы апроксимации уравнений (3) - (5). Затем рассчитываются элементы вектора переменных Q и, соответственно, определяются все параметры течения.
Турбулентность моделировалась методом крупных вихрей, являющимся одним из наиболее перспективных методов расчетов турбулентных течений. Данный метод широко используется в практике создания численных методик. Например, он применялся в работах [16-18]. Для моделирования пристеночных течений использовалась локальная модель вихревой вязкости (Wall-Adapting Local Eddy-viscosity, WALE) [17].
3. Результаты расчетов
Расчеты проводились для модельного РДТТ, представленного на рис. 2. Ракетный двигатель имеет заряд твердого топлива зонтичного типа. Длина двигателя составляла 2,5 м. В начальный момент времени температура газа в камере составляла T0 = 294 K и давление p0 = 106 Па. Шаг по времени задавался равным At = 5, 0• 10-7c. Скорость горения топлива определялась по формуле ш = ш0 (p/105) , где ш0 - скорость горения при атмосферном давлении, $ - константа. При апробации вычислительной методики полагалось w = 0.
На рис. 3 а), 3 б), 4 а) приведены графики изменения осевой скорости, радиальной скорости и давления в камере на начальном отрезке времени работы двигателя в четырех точках камеры двигателя. В начальной стадии работы двигателя на графиках хорошо заметен ударно-волновой характер процессов, протекающих в камере РДТТ, это связано с отражением волн от элементов камеры и поверхности топлива. Процессы, протекающие в двигателе на начальной стадии работы, существенно нестационарны. Величины параметров потока имеют сильно изменяющийся характер, особенно это заметно на графиках параметров в зазоре между передним днищем и зарядом. На рис. 4 б) приведено изменение давления по времени до выхода двигателя на стационарный режим работы. Время выхода на режим составило Ь ~ 1,1 с. Для этого момента времени на рис. 5, 7 представлены распределения газодинамических параметров течения продуктов сгорания по всему тракту двигателя, как в камере, так и в сопле. На рис. 5 приведены линии тока течения, из которых видно образование вихрей в зонтичной области камеры и областью за утопленной частью сопла. Распределения осевой скорости и давления приведены на рис. 6, 7, на которых наблюдаются большие градиенты параметров потока продуктов сгорания.
Рис. 2. Схема твердотопливного двигателя
а) б)
Рис. 3. Изменение осевой (а) и радиальной (б) скоростей по времени: 1 - в зазоре между передним днищем и зарядом; 2 - на оси симметрии переднего днища; 3 - за утопленной частью сопла; 4 - в зонтичной части
0.005
а) б)
Рис. 4. Изменение давления по времени: а) 1 - в зазоре между передним днищем и зарядом; 2 - на оси симметрии переднего днища; 3 - за утопленной частью сопла; 4 - в зонтичной части; б) изменение давления на оси симметрии переднего днища до выхода на стационарный режим работы
Таким образом, предложенная методика позволяет проводить сквозной расчет течения в РДТТ, т.е. моделирование течений в дозвуковых и сверхзвуковых областях проводится с использованием единого алгоритма.
Заключение
На основе модифицированного метода расщепления потоков предложен алгоритм и разработана программа, позволяющая проводить сквозное моделирование внутренних нестационарных турбулентных течений сжимаемого вязкого газа в тракте РДТТ. Проведены тестовые численные исследования выхода модельного твердотопливного ракетного двигателя на стационарный режим работы. Представлены результаты мо-
Рис. 5. Линии тока течения газа
и, м/с
I2.45E+03 2.42Е+03 2.30Е+03 2.00Е+03 1.60Е+03
— 1.20Е+03
— 9.00Е+02 6.00Е+02 3.00Е+02
=J 0.00Е+00
р, Па
7.3Е+06 7.2Е+06 7.1Е+06 6.7Е+06 6.1Е+06 З.ЗЕ+06 9.6Е+05 4.0Е+05 2.0Е+05 0.0Е+00
Рис. 7. Распределение давления
делирования внутрикамерных процессов, протекающих в тракте двигателя. В результате расчетов определено время выхода двигателя на стационарный режим.
Литература
1. Липанов, А.М. Численный эксперимент в теории РДТТ / А.М. Липанов, В.П. Бобрышев, А.В. Алиев, Ф.Ф. Спиридонов, В.Д. Лисица. - Екатеринбург: УИФ Наука, 1994.
2. Ciucci, A. Simulation of Rocket Motor Internal Flows with Turbulent Mass Injection / A. Ciucci, G. Iaccarino, R. Moser, F. Najjar, P. Durbin // Center for Turbulence Research. University of Stanford. - 1998. - P. 245-266.
3. Алиев, А.В. Нестационарные внутрикамерные процессы в твердотопливных регулируемых двигательных установках / А.В. Алиев, О.В. Мищенкова, И.В. Черепов // Вестник МГТУ им. Н.Э. Баумана. Серия: Машиностроение. - 2016. - № 4. - С. 24-39.
4. Apte, S.V. A Large-Eddy Simulations Study of Transition and Flow Instability in a Porous-Walled Chamber with Mass Injection / S.V. Apte, V. Yang // Journal of Fluid Mechanics. -2003. - V. 477. - P. 215-225.
5. Липанов, А.М. Численный метод расчета турбулентных течений и теплообмена в двигателях летательных аппаратов / А.М. Липанов, Ю.Ф. Кисаров, И.Г. Ключников // Прикладная математика и механика. - 1992. - № 6. - С. 49-53.
6. Vuillot, F. Combustion and Turbulent Flow Effects in 2D Unsteady Navier-Stokes Simulations of Oscillatory Solid Rocket Motors / F. Vuillot, N. Lupoglazoff // AIAA Paper. -1996. - № 96-0884. - 15 p.
7. Липанов, А.М. Адаптированные цилиндрические координаты для внутренних объемов элементов конструкции РДТТ / А.М. Липанов, М.Р. Королева, С.Ю. Дадикина // Труды института математики и механики УрО РАН. - 2012. - Т. 18, № 1. - С. 213-221.
8. Волков, К.Н. Теплообмен в каверне с вращающимся диском в турбулентном режиме / К.Н. Волков, П.В. Булат, И.А. Волобуев, В.А. Пронин // Научно-технический вестник информационных технологий, механики и оптики. - 2017. - Т. 17, № 3. - С. 514-524.
9. Липанов, А.М. Теоретическая гидромеханика ньютоновских сред / А.М. Липанов. - М.: Наука, 2011.
10. Fluent 6.3 User's Guide, 2006. - URL: www.sharcnet.ca/Software/Fluent6/html/ug/ main_pre.htm
11. Steger, J.L. Flux Vector Splitting of the Inviscid Gasdynamic Equations with Application to Finite Difference Methods / J.L. Steger, R.F. Warming // Journal of Computational Physics. - 1981. - V. 40, № 2. - P. 263-293.
12. Anderson, W.K. A Comparison of Finite Volume Flux Vector Splittings for the Euler Equations / W.K. Anderson, J.L. Thomas, B. Van Leer // AIAA Journal. - 1986. - V. 24, № 9. - P. 1453-1460.
13. Van Leer, B. Towards the Ultimate Conservative Difference Scheme V. A Second Order Sequel to Godunov's Method / B. Van Leer // Journal of Computational Physics. - 1979. - V. 32, № 1. - P. 101-136.
14. Суров, В.С. Метод С.К. Годунова для многоскоростной модели гетерогенной среды / В.С. Суров, И.В. Березанский // Вестник ЮУрГУ. Серия: Математическое моделирование и программирование. - 2014. - Т. 7, № 2. - С. 87-98.
15. Chung, T.J. Computational Fluid Dynamics / T.J. Chung. - Cambridge: Cambridge University Press, 2002.
16. Akselvoll, K. Large-Eddy Simulation of Turbulent Confined Coannular Jets / K. Akselvoll, P. Moin // Fluid Mech. - 1996. - V. 315. - P. 387-411.
17. Ducros, F. Wall-Adapting Local Eddy-Viscosity Models for Simulations in Complex Geometries / F. Ducros, F. Nicoud, T. Poinsot // Proceeding of the 6th ICFD Conference on Numerical Methods for Fluid Dynamic, Oxford, United Kingdom. - 1998. - P. 293-299.
18. Шумихин, А.А. Использование схемы WENO для моделирования турбулентного течения в канале с обратным уступом / А.А. Шумихин, М.Р. Королева, С.В. Дадикина, А.И. Карпов // Вестник Удмуртского университета. Математика. Механика. Компьютерные науки. - 2017. - Т. 27, № 3. - С. 460-469.
Алексей Матвеевич Липанов, академик РАН, доктор технических наук, профессор, главный научный сотрудник, Институт прикладной математики им. М.В. Келдыша РАН (г. Москва, Российская Федерация), [email protected].
Светлана Юрьевна Дадикина, младший научный сотрудник, Институт механики, Удмуртский федеральный исследовательский центр УрО РАН (г. Ижевск, Российская Федерация), [email protected].
Андрей Александрович Шумихин, кандидат физико-математических наук, научный сотрудник, Институт механики, Удмуртский федеральный исследовательский центр УрО РАН (г. Ижевск, Российская Федерация), [email protected].
Мария Равилевна Королева, кандидат физико-математических наук, старший научный сотрудник, Институт механики, Удмуртский федеральный исследовательский центр УрО РАН (г. Ижевск, Российская Федерация), [email protected].
Александр Иванович Карпов, доктор физико-математических наук, главный научный сотрудник, Институт механики, Удмуртский федеральный исследовательский центр УрО РАН (г. Ижевск, Российская Федерация), [email protected].
Поступила в редакцию 21 сентября 2018 г.
MSC 76F65, 65N08 DOI: 10.14529/mmp190103
NUMERICAL SIMULATION INTRA-CHAMBER OF UNSTEADY TURBULENT FLOWS STIMULATE. PART 1
A.M. Lipanov1, S.Yu. Dadikina2, A.A. Shumikhin2, M.R. Koroleva2, A.I. Karpov2
1Keldysh Institute of Applied Mathematics, Russian Academy of Science, Moscow, Russian Federation
2Udmurt Federal Research Center, Ural Branch, Russian Academy of Science, Izhevsk, Russian Federation
E-mails: [email protected], [email protected], [email protected], [email protected], [email protected]
The technique of 3D internal unsteady turbulent flows simulation is proposed in this work, in particular the flow of combustion products into the solid fuel rocket engine. The system of governing equations describing the flow of viscous compressible gas written in the cylindrical coordinate system is presented. The computational algorithm based on a modified scheme of splitting vectors belonging to the class methods use Godunov approach is proposed. This algorithm is suitable for end-to-end calculation of the internal flow around all the paths of the rocket engine, including both subsonic flow in the chamber zone, and the zone of supersonic flow in nozzle. The simulation results of gas flow in the model rocket engine show oscillating shock wave processes the occurring into the engine chamber at the early stage. The time-stationary working mode engine is defined.
Keywords: intra-chamber of the processes; turbulence; unsteady flow; computational fluid dynamics.
References
1. Lipanov A.M., Bobryshev V.P., Aliev A.V., Spiridonov F.F., Lisitsa V.D. Chislennyi eksperiment v teorii RDTT [Numerical Experiment in the Theory of Solid Propellant Rocket Motors]. Ekatirinburg, UIF Nauka, 1994. (in Russian)
2. Ciucci A., Iaccarino G., Moser R., Najjar F., Durbin P. Simulation of Rocket Motor Internal Flows with Turbulent Mass Injection. Center for Turbulence Research. University of Stanford, 1998, pp. 245-266.
3. Aliev A.V., Mishchenkova O.V., Cherepov I.V. Nonstationary Intra-Chamber Processes in Solid-Propellant Controlled Propulsion System. Herald of the Bauman Moscow State Technical University. Series: Mechanical Engineering, 2016, no. 4, pp. 24-39. DOI: 10.18698/0236-3941-2016-4-24-39
4. Apte S.V., Yang V. A Large-Eddy Simulations Study of Transition and Flow Instability in a Porous-Walled Chamber with Mass Injection. Journal of Fluid Mechanics, 2003, vol. 477, pp. 215-225. DOI: 10.1017/S0022112002002987
5. Lipanov A.M., Kisarov Yu.F., Klyuchnikov I.G. Numerical Method for Calculating Turbulent Flows and Heat Transfer in Aircraft Engines. Journal of Applied Mathematics and Mechanics, 1992, no. 6, pp. 49-53.
6. Vuillot F., Lupoglazoff N. Combustion and Turbulent Flow Effects in 2D Unsteady Navier-Stokes Simulations of Oscillatory Solid Rocket Motors. AIAA Paper, 1996, no. 96-0884, 15 p.
7. Lipanov A.M., Koroleva M.R., Dadikina S.Yu. Adapted Cylindrical Coordinates for Internal Volumes of Structural Elements of a Solid-Propellant Rocket Engine. Proceedings of the Steklov Institute of Mathematics (Supplementary Issues), 2012, vol. 18, no. 1, pp. 213-221.
8. Volkov K.N., Bulat P.V., Volobuev I.A., Pronin V.A. Heat Transfer in a Cavity with Rotating Disk in Turbulent Regime. Scientific and Technical Journal of Information Technologies, Mechanics and Optics, 2017, vol. 17, no. 3, pp. 514-524. (in Russian) DOI: 10.17586/22261494-2017-17-3-514-524
9. Lipanov A.M. Teoreticheskaya gidromekhanika n'yutonovskikh sred [Theoretical Fluid Mechanics of Newtonian Media]. Moscow, Nauka, 2011. (in Russian)
10. Fluent 6.3 User's Guide, 2006. Available at: www.sharcnet.ca/Software/Fluent6/html/ug/ main_pre.htm
11. Steger J.L., Warming R.F. Flux Vector Splitting of the Inviscid Gasdynamic Equations with Application to Finite Difference Methods. Journal of Computational Physics, 1981, vol. 40, no. 2, pp. 263-293. DOI: 10.1016/0021-9991(81)90210-2
12. Anderson W.K., Thomas J.L., Van Leer B. A Comparison of Finite Volume Flux Vector Splittings for the Euler Equations. AIAA Journal, 1986, vol. 24, no. 9, pp. 1453-1460. DOI: 10.2514/3.9465
13. Van Leer B. Towards the Ultimate Conservative Difference Scheme V. A Second Order Sequel to Godunov's Method. Journal of Computational Physics, 1979, vol. 32, no. 1, pp. 101-136. DOI: 10.1016/0021-9991(79)90145-1
14. Surov V.S., Berezancky I.V. Godunov's Method for a Multivelocity Model of Heterogeneous Medium. Bulletin of the South Ural State University. Series: Mathematical Modelling, Programming and Computer Software, 2014, vol. 7, no. 2, pp. 87-98. (in Russian) DOI: 10.14529/mmp140208
15. Chung T.J. Computational Fluid Dynamics. Cambridge, Cambridge University Press, 2002. DOI: 10.1017/CBO9780511606205
16. Akselvoll K., Moin P. Large-Eddy Simulation of Turbulent Confined Coannular Jets. Journal of Fluid Mechanics, 1996, vol. 315, pp. 387-411. DOI: 10.1017/S0022112096002479
17. Ducros F., Nicoud F., Poinsot T. Wall-Adapting Local Eddy-Viscosity Models for Simulations in Complex Geometries. Proceeding of the 6th ICFD Conference on Numerical Methods for Fluid Dynamic, Oxford, United Kingdom, 1998, pp. 293-299.
18. Shumikhin A.A., Koroleva M.R., Dadikina S.Yu., Karpov A.I. Application of WENO Scheme for Simulation of Turbulent Flow in a Channel with Backward-Facing Step. Vestnik Udmurtskogo universiteta. Matematika. Mekhanika. Komp'yuternye nauki [The Bulletin of Udmurt University. Mathematics. Mechanics. Computer Science], 2017, vol. 27, no. 3, pp. 460-469. (in Russian) DOI: 10.20537/vm170313
Received September 21, 2018
