УДК 532.517.4
DOI: 10.14529/ mmp230103
ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ВНУТРИКАМЕРНЫХ НЕСТАЦИОНАРНЫХ ТУРБУЛЕНТНЫХ ТЕЧЕНИЙ. ЧАСТЬ 2
А.М. Липанов1, А.А. Шумихин2, М.Р. Королева,2, А.И. Карпов2
Институт прикладной математики им. М.В. Келдыша РАН, г. Москва, Российская Федерация
2Институт механики УдмФИЦ УрО РАН, г. Ижевск, Российская Федерация
Течение потока газа в твердотопливном ракетном двигателе определяется особенностями физико-химических процессов, протекающих в камере сгорания, и процесса истечения газа из сопла. В работе предложена методика моделирования внутренних нестационарных турбулентных течений в ракетном двигателе с зарядом твердого топлива телескопического типа. Алгоритм разработан на основе системы уравнений сохранения гидромеханических параметров, описывающих поток сжимаемого вязкого газа. Система уравнений была записана с использованием цилиндрической системы координат. Представленная численная методика относится к классу алгоритмов, использующих подход Годунова. Основой разработанного алгоритма является схема расщепления векторов потоков, модифицированная для вязких течений. Предложенный алгоритм позволяет производить сквозной расчет течения продуктов сгорания твердого топлива по всему тракту ракетного двигателя. Результаты, полученные в ходе численного моделирования потока в ракетном двигателе, позволяют сделать анализ зависимости температуры газа на стенке двигателя от скорости горения низкотемпературной внешней шашки телескопического заряда.
Ключевые слова: внутрикамерные процессы; турбулентность; нестационарное течение; вычислительная гидрогазодинамика.
Введение
Современные твердотопливные ракетные двигатели (РДТТ) - это сложные объекты машиностроения. Этот аспект представляется важным и при проектировании двигателей, и при их отработке, и при их исследовании с использованием математических моделей. Ракетные двигатели, работающие на твердом топливе, характеризуются сложными физико-химическими и газодинамическими процессами, протекающими по всему тракту РДТТ [1-3]. Внутикамерные процессы в двигателе определяются особенностями движения газа в камере сгорания и при истечении газа из сопла. От точности методов, используемых для моделирования течения продуктов сгорания, в значительной степени зависит полнота изучения газодинамических процессов, протекающих в двигателе, и как следствие, энергомассовое совершенство и надежность проектируемых ракетных двигателей. Таким образом, для исследования течения продуктов сгорания в РДТТ важное значение имеет численное моделирование потоков вязкого сжимаемого газа. Предложенный алгоритм численного моделирования разработан на основе системы уравнений сохранения гидромеханических параметров потока газа. Система уравнений была записана с использованием цилиндрической системы координат [2,4-10]. Основой разработанной программы является схема Стегера - Уорминга. Алгоритм позволяет производить сквозной расчет течения продуктов сгорания твердого топлива по всему тракту ракетного двигателя. Конструкция заряда должна обеспечить максимальную защиту корпуса двигателя от воздействия продуктов сгорания топлива. Одним из решений, позволяющим защитить корпус от высокой температуры продуктов сгорания, является использование низкотемпературной газовой завесы. Для этого в зарядах телескопического типа внешнюю шашку изготавливают из низкотемпературного топлива. В работе приведены результаты исследований двигателя, использующего заряд данного типа.
1. Постановка задачи
Алгоритм моделирования разработан с использованием системы уравнений сохранения гидромеханических параметров потока вязкого сжимаемого газа в цилиндрических координатах, имеющей следующий вид [2,4-8]:
dQ dFx dFr 1 dFe 1 ( д т д д \ т
ir + 7Г1 + 7Г1 + —+ Feo = -[ тrrLx + —rLr + —Le )+Lee, 1 dt дх дг г дв r \дх дг дв )
где t - время; Q - вектор переменных; Fx, Fr, Fe, Feo - векторы конвективных потоков; Lx, Lr, Lo, Loo - векторы вязких напряжений; х, г, в - продольная, радиальная и угловая координаты. Векторы переменных и конвективных потоков описаны в работах [4,5,7]. Векторы вязких напряжений описаны в работах [6-8,10].
В формулах, описывающих компоненты векторов вязких напряжений, использовался коэффициент эффективной вязкости veff, который вычисляется как сумма коэффициента молекулярной вязкости v и коэффициента турбулентной подсеточной вязкости vsgs
Veff = V + Vsgs. (2)
Подсеточная вихревая вязкость vsgs, с использованием модели Смагоринского, обобщенной для сжимаемых течений [11,12], определялась как функция тензора скоростей деформаций и размера фильтра
Vsgs = Р (CsA)' IS|, (3)
í Dn ■ Зи' \
где \S\ = л/2SíjSíj, Su = - ( тг—- + т-^- ) - тензор скоростей деформаций, и - компо-v 2 \дxj дхг)
ненты вектора скорости, р - плотность, Cs - параметр Смагоринского, A - ширина фильтра. Ширина фильтра рассчитывается по формуле A = (AxiAx2 Ax3)1/3 = V1/3, где Ax1, Ax2, Ax3 - размеры ячеек сетки по соответствующим осям.
Из (3) видно, что величина подсеточной вязкости vsgs зависит от квадрата ширины фильтра (размеров расчетной ячейки). Таким образом, при уменьшении размеров ячейки вклад vsgs в суммарную вязкость veff также будет квадратично уменьшаться, т. е. при стремлении A ^ 0 метод крупных вихрей переходит в прямое численное моделирование.
2. Вычислительный алгоритм
Для выполнения численных исследований предложен вычислительный алгоритм 2-го порядка точности по пространству. Методика разработана на базе схемы Стегера - Уорминга (метода расщепления векторов потоков), модифицированной для вязких течений. Основной идеей данной схемы является расщепление вектора конвективных потоков, в частности компоненты Fx, на две части
Рх = F+ + F-, (4)
где вектор F+ связан с вектором, который формируется из положительных собственных чисел матрицы Якоби А+, и соответственно вектор Fx- связан с вектором, компоненты которого формируются из отрицательных собственных чисел А-. Схема Стегера - Уорминга описана в работах [13, 14], а также была использована в исследовании [5].
Система основных уравнений (1) по пространству дискретизировалась с использованием метода конечных объемов. В интегральной форме система уравнений записывается так:
5 [пМ-L. [ fdF* Л. dFr + ldF° + Г \ AV
— \ QdV + / — + — + —— + Fee dV = dt J J \ dx dr r dO J
Vi Vi
Vi
(5)
If д д д \ \ ( )
- -¿~rLx + — rLr + —-Le ) + Lee ) dV r \dx dr dO
где V - объем ¿-го конечного объема. С использованием теоремы Остроградского -Гаусса система (5) преобразуется к следующему виду:
J QdV + £ (^хпх + Ргпг + ^вп^ dS + J FeedV =
(6)
i I
(гЬхпх + гЬгпг + ЬеПе) dS+ LeedV, г г ] 3
З V
где пх, пг, Не - компоненты единичной внешней нормали п к поверхности ¿-го конечного объема Бг, гг - координата центра масс ¿-го конечного объема. Для течений с угловой скоростью = 0 (осесимметричное течение) выражение (6) приводится к виду
— J QdS+£ (Fxnx + Frnr)dC+^ FeedS = - ^ (гЬхпх + гЬгпг) dC+^ LeedS,
Зг сг Зг С г Зг
где Бг - площадь ¿-го конечного объема, Сг - контур ¿-го конечного объема. В дискретном виде уравнение (7) можно записать
х теЕг (8)
^ ^ {ХгтРХ:гтПх^т Тг,гт) ДСгт ¿б^г^г; г т&Ег
где т - номер соседнего элемента, Ег = {], к, 1} - множество соседних элементов ¿-го элемента, пх,гт, пг,гт - составляющие внешней нормали к ¿т-й стороне ¿-го элемента, ДСгт, ггт - длина и координата центра ¿т-й стороны ¿-го элемента, Fее,i, Ьее,г -значения компонент в центре ¿-го элемента. Элементы показаны на рис. 1.
Тогда приращение вектора переменных ДQi на текущей итерации по времени можно определить
^ ^ (^imLXiirnT^x,irn ТimLг,im) ДС\m Leeti J ■
i , , \ f'ini 1 J r in I I '
ri Si
m€Ei
Рис. 1. Параметры потока газа Q+•, Qij, Q+k, , Q¿, в центрах граней г-го конечного объема
Конвективные потоки на гранях объемов рассчитывались с применением алгоритма 2-го порядка точности. Для определения векторов параметров на грани г3, определялись параметры газа в центре грани. Значения р+ и вычислялись следующим образом:
др
др
Pi I I Tf ) • Ax+ + ( Tf
dx др
dr др
• Дг+,
Ax+j xij xi
Ari+ = rij - ri,
Pr:, = Pj + ( ^ ) • A^. + ( — ) • Ar
dr
Axjj — xij - xj, Ar,] = rij - rj,
где XI, г - координаты центра масс г-го конечного объема; х$, г$ - координаты центра масс з-го конечного объема; х^, г^ - координаты центра грани г3. Остальные параметры газа в центре грани рассчитывались таким же образом. Для определения параметров газа в центрах граней гк и И конечного объема производились аналогичные вычисления. Далее, учитывая (4), вычислялись собственные числа и потоки на гт-й стороне г-го элемента в (9)
^x,im(Qim )
Fx,jm — F+im(Q+m> )) + Fx,im(Qim, ^x,im(Qim)) ,
(Q+m)+|A x,im (Q+m )l ^ _ x,im (Qim)
x,im
2 ' r.iin \ ^im > 2
(10) (11)
Значения потоков Fx вычислялись с использованием выражений (10), (11) для всех граней контрольных объемов расчетной области. Вычисление значений потоков вектора Fr производились аналогичным образом. Расчет значения вектора Fqq производился с использованием значений вектора параметров в центрах масс контрольных объемов. Для определения вязких потоков в уравнении (9) рассчитывались компоненты тензора вязких напряжений [2,6-9]. В частности, компонента Txxjim на im-й стороне i-го элемента определялась с использованием выражений
(ди\ 2
Тхх,гт = 2/ie// (-J ~ -fleff (dlvQ)im ,
fdu\ 1 ffdu\ fdu\ \
~ 2 {{diJi + \d~x),J '
(12) (13)
+
j
j
^ Х г * теБ;
где их,гт - компонента скорости на гш-й стороне г-го элемента, их^ - компонента скорости в центре г-го элемента, ихт - компонента скорости в центре ш-го соседнего элемента. Остальные компоненты тензора напряжений вычислялись аналогично. Компонента гвв^ в центре г-го элемента рассчитывалась так:
V- 2
Ч = ^е//"1 - ^Ме// ((¿Ш^ , (16)
Г г 3
где ^ = (и; V; т) - вектор скорости.
Для дискретизации системы уравнений сохранения гидромеханических параметров потока газа (1) по времени использовалась явная схема 1-го порядка точности (схема Эйлера, одношаговая явная двухслойная схема). На текущем временном слое конвективные и вязкие потоки ГХП, ГП, ГП, ГПв и ЬП, Ц1, ЬП, ЬПв вычислялись с использованием значений вектора переменных также на текущем временном слое Qn. Далее на текущей итерации по времени рассчитывалось приращение вектора переменных AQn. В итоге определялись элементы вектора переменных Qn+1 и вычислялись все параметры газа на п +1 временном слое.
Qn+l = + AQn. (17)
При программировании применялась методика, описанная в работе [15]. Программа, разработанная с использованием данной методики, позволяет проводить параллельные вычисления на компьютерах, имеющих многоядерный центральный процессор.
3. Результаты расчетов
Расчеты проводились для двухрежимного модельного РДТТ. Схема двигателя показана на рис. 2. Исследуемый двигатель имеет составной заряд твердого топлива. Первая часть является зарядом торцевого горения (обеспечивает первый режим работы двигателя), вторая - зарядом телескопического типа (обеспечивает второй режим работы двигателя). Моделирование течения продуктов сгорания производилось для второго режима работы двигателя, т. е. моделировалось течение при горении телескопической части заряда, когда шашка торцевого горения уже выгорела. Торцевая поверхность шашек телескопического заряда была забронирована. Расчетная область показана на рис. 3.
Рис. 2. Схема твердотопливного двигателя
Рис. 3. Схема расчетной области
В расчетах шаг по времени принимался At = 2,0 • 10 8 с. Скорость горения твердого топлива вычислялась по формуле и = и0 (р/105 , где и0 - скорость горения при атмосферном давлении, $ - константа. Начальные условия в камере двигателя для температуры газа и давления задавалась соответственно Т0 = 2000 К и р0 = 106 Па. Площадь поверхности горения заряда и его геометрия в ходе расчетов не менялись. Расчеты проводились при трех значениях скорости горения низкотемпературной верхней топливной шашки и0 = 0, 002; 0, 005; 0, 008 м/с. Для остальных шашек скорость горения во всех расчетах составляла и0 = 0, 005 м/с. Константа в законе горения принималась $ = 0, 6. Температура продуктов сгорания низкотемпературной верхней шашки составляла Тс2 = 2000 К, остальных шашек Тс = 3300 К. Удельная газовая постоянная составляла Я = 275 Дж/(кг-К), показатель адиабаты 7 =1,18, параметр Смагоринского Оз = 0, 3. Также полагалось, что угловая скорость была равной т = 0 м/с во всех ячейках расчетной сетки.
Длина двигателя составляла 1,0 м. Радиус двигателя был равен 0,3 м. Радиус критического сечения сопла задавался 0,1 м. На входной границе расчетной области (поверхность горения шашек, с которой происходит вдув продуктов сгорания твердого топлива) начальных возмущений параметров течения не вводилось. Расчеты проводились на неструктурированной расчетной сетке (рис. 4), построенной из треугольных элементов. Расчет производился на сетках с различным количеством элементов.
(а)
(б)
Рис. 4. Расчетная сетка: (a) - вариант 6400 ячеек; (б) - вариант 267700 ячеек (фрагмент)
На рис. 4(а) показан грубый вариант сетки (6400 элементов), на рис. 4(б) показан фрагмент самой мелкой сетки (262700 элементов). Адаптации сетки к течению газа в процессе расчетов не производилось. Для создания расчетной сетки использовался бесплатный сеточный генератор GMesh (ver. 4.4.1), с открытым исходным кодом (http://geuz.org/gmsh/ ), имеющий GNU GPL лицензию (свободную copyleft лицензию для программного обеспечения).
На рис. 5 приведены распределения температуры для 3х скоростей горения низкотемпературной шашки и0 = 0, 002; 0, 005; 0, 008 м/c (сетка 262700 ячеек). На рис. 6 приведено распределение температуры для и0 = 0, 005 м/c (сетка 6400 ячеек).
(в)
Рис. 5. Распределение температуры при различных скоростях горения низкотемпературной шашки (£ = 0, 5 с, сетка 262700 ячеек): (а) и0 = 0, 002м/с; (б) и0 = 0, 005м/с; (в) ^о = 0, 008 м/с
Рис. 6. Распределение температуры при скорости горения низкотемпературной шашки и0 = 0, 005 м/с (£ = 0, 5 с, сетка 6400 ячеек)
Далее приводятся результаты расчетов при скорости горения низкотемпературной шашки и0 = 0, 005м/с (сетка 262700 ячеек). Линии тока течения газа для 3х моментов времени приведены на рис. 7.
(в)
Рис. 7. Линии тока течения газа для 3х моментов времени (сетка 262700 ячеек): (а) г = 0,4 с; (б) г = 0,45 с; (в) г = 0, 5 с
На рис. 8 приведено распределение завихренности вектора скорости. На рис. 9 приведены графики изменения температуры по времени в двух точках камеры сгорания двигателя. На рис. 10 приведена осредненная по времени температура потока в 2х сечениях при 3х скоростях горения низкотемпературной шашки. Сравнивая рис. 5 и рис. 6, можно отметить, что в случае проведения расчетов на грубой сетке (6400 ячеек) распределение температуры газа в потоке носит <слоистый> характер. В отличие от этого, при расчетах на мелкой сетке (262700 ячеек) в течении появляются вихревые структуры, приводящие к интенсивному перемешиванию потоков газа, имеющих различную температуру.
rotQ., с 1 1.000Е+04 7.143Е+03 4.286Е+03 1.429Е+03 -1.429Е-03 -4.286Е-03 -7.143Е-03 -1.000Е-04
Рис. 8. Распределение завихренности вектора скорости тог& (г = 0,5 с, сетка 262700 ячеек)
0.45 0.46 0.47 0.48 0.49 с 0.45 0.46 0.47 0.48 0.49 г, с
(а) (б)
Рис. 9. Изменение температуры по времени в 2х точках (при х = 0,4 м, сетка 262700 ячеек): (а) г = 0,15 м; (б) г = 0, 27м
2200 2400 2600 2800 3000 Т, К
(а)
Рис. 10. Осредненная по времени температура при 3х скоростях горения низкотемпературной шашки (1 - и0 = 0, 002 м/с; 2 - и0 = 0, 005 м/с; 3 - и0 = 0, 008 м/с) в двух сечениях: (а) х = 0, 2 м; (б) х = 0, 4 м
Заключение
В работе предложен вычислительный алгоритм для проведения численных исследований внутренних нестационарных турбулентных течений сжимаемого вязкого газа. С использованием данного алгоритма разработана программа и проведены тестовые расчеты. Моделирование течения в двигателе позволяет определить величину температуры газа у стенки двигателя в зависимости от скорости горения низкотемпературной шашки заряда. Из полученных результатов видно, что хотя горячий газ из ядра течения увеличивает температуру газа около поверхности корпуса двигателя за счет вихревого перемешивания, уменьшить температуру газа на поверхности корпуса в значительной степени можно, увеличивая скорость горения низкотемпературной шашки.
Литература
1. Соломонов, Ю.С. Твердотопливные регулируемые двигательные установки j Ю.С. Соломонов, А.М. Липанов, А^. Алиев, А.А. Дорофеев, B^. Черепов. - M.: Машиностроение, 2011.
2. Липанов, А.М. Численный метод расчета турбулентных течений и теплообмена в двигателях летательных аппаратов j А.М. Липанов, Ю.Ф. Кисаров, И.Г. Ключников jj Известия вузов. Авиационная техника. - 1988. - № 1. - С. 49-53.
3. Алиев, А^. Bнутренняя баллистика ракетных двигателей на твердом топливе j А^. Алиев, Г.Н. Амарантов. - M.: Машиностроение, 2007.
4. Годунов, С.К. Численное решение многомерных задач газовой динамики j С.К. Годунов, А^. Забродин, М.Я. Иванов, А.Н. Крайко, Г.П. Прокопов. - M.: Наука, 1976.
5. Липанов, А.М. Численное моделирование внутрикамерных нестационарных турбулентных течений. Ч. 1 j А.М. Липанов, С.Ю. Дадикина, А.А. Шумихин, М.Р. Королева, А.И. Карпов j j Bестник ЮУрГУ. Серия: Математическое моделирование и программирование. - 2019. - Т. 12, № 1. - С. 32-43.
6. Тадмор, З. Теоретические основы переработки полимеров j З. Тадмор, К. Гогос. - M.: Химия, 1984.
7. Липанов, А.М. Теоретическая гидромеханика ньютоновских сред j А.М. Липанов. - M.: Наука, 2011.
8. Коммисаров, Ю.А. Процессы и аппараты химической технологии. Ч. 1 j Ю.А. Комми-саров, Л.С. Гордеев, Д.П. Bент. - M.: Изд-во Юрайт, 2018.
9. Лыков, А^. Тепломассообмен j А^. Лыков. - M.: Энергия, 1971.
10. Шумихин, А.А. Численное моделирование течения вязкого сжимаемого газа в ракетных двигателей на твердом топливе с центральным телом j А.А. Шумихин, С.Ю. Дадикина jj Химическая физика и мезоскопия. - 2020. - Т. 22, № 2. - С. 184-196.
11. Pino, M.M. Subgrid-Scale Models for Compressible Large-Eddy Simulations j M.M. Pino, U. Piomelli, G.V. Candler j j Theoretical and Computational Fluid Dynamics. - 2000. -V. 13, № 5. - P. 361-376.
12. Bолков, К.Н. Моделирование крупных вихрей полностью развитого турбулентного течения в канале и сравнение моделей подсеточной вихревой вязкости j К.Н. Bолков j j Прикладная механика и техническая физика. - 2006. - Т. 47, № 3. - С. 31-42.
13. Steger, J.L. Flux Vector Splitting of the Inviscid Gasdynamic Equations with Application to Finite Difference Methods j J.L. Steger, R.F. Warming j j Journal of Computational Physics. - 1981. - V. 40, № 2. - P. 263-293.
14. Anderson, W.K. A Comparison of Finite Volume Flux Vector Splittings for the Euler Equations j W.K. Anderson, J.L. Thomas, B. Van Leer jj AIAA Journal. - 1986. - V. 24, № 9. - P. 1453-1460.
15. Шумихин, А.А. Моделирование турбулентных течений вязкого сжимаемого газа с использованием параллельных вычислений j А.А. Шумихин j j Химическая физика и мезоскопия. - 2021. - Т. 23, № 3. - С. 292-302.
Алексей Матвеевич Липанов, академик РАН, главный научный сотрудник, Институт прикладной математики им. М.В. Келдыша РАН (г. Москва, Российская Федерация), aml35@yandex.ru.
Андрей Александрович Шумихин, кандидат физико-математических наук, научный сотрудник, Институт механики УдмФИЦ УрО РАН (г. Ижевск, Российская Федерация), shumihin@udman.ru.
Мария Равилевна Королева, кандидат физико-математических наук, старший научный сотрудник, Институт механики УдмФИЦ УрО РАН (г. Ижевск, Российская Федерация), koroleva-ipm@yandex.ru.
Александр Иванович Карпов, доктор физико-математических наук, главный научный сотрудник, Институт механики УдмФИЦ УрО РАН (г. Ижевск, Российская Федерация), karpov@udman.ru.
Поступила в редакцию 25 апреля 2022 г.
MSC 76F65, 65N08 DOI: 10.14529/mmp230103
NUMERICAL SIMULATION INTRA-CHAMBER OF UNSTEADY TURBULENT FLOWS STIMULATE. PART 2
A.M. Lipanov1, A.A. Shumikhin2, M.R. Koroleva2, A.I. Karpov2
1Keldysh Institute of Applied Mathematics, Russian Academy of Science, Moscow,
Russian Federation
2Institute of Mechanics, Udmurt Federal Research Center, Ural Branch, Russian
Academy of Science, Izhevsk, Russian Federation
E-mail: aml35@yandex.ru, shumihin@udman.ru, koroleva-ipm@yandex.ru,
karpov@udman.ru
The flow of gas in a solid-fuel rocket engine is determined by the peculiarities of the physico-chemical processes occurring in the combustion chamber and the process of gas outflow from the nozzle. The paper proposes a method for modeling internal unsteady turbulent flows in a rocket engine with a solid fuel charge of a telescopic type. A system of defining equations written in a cylindrical coordinate system describing the flow of a compressible viscous gas is given. A computational algorithm is proposed that belongs to the class of methods using the Godunov approach, developed on the basis of a modified flow vector splitting scheme. The obtained results of a numerical study of the flow in a model rocket engine show the dependence of the gas temperature on the engine wall on the combustion rate of a low-temperature external charge of a telescopic charge.
Keywords: intra-chamber processes; turbulence; unsteady flow; computational fluid dynamics.
References
1. Solomonov Yu.S., Lipanov A.M., Aliyev A.V., Dorofeev A.A., Cherepanov V.I. Tverdotoplivnye reguliruemye dvigatelnye ustanovki [Regulable Solid-Propellant Rocket Engines]. Moscow, Mashinostroenie, 2011. (in Russian)
2. Lipanov A.M., Kisarov Yu.F., Klyuchnikov I.G. [Numerical Method for Calculating Turbulent Flows and Heat Transfer in Aircraft Engines]. Izvestiya vuzov. Aviacionnaya tekhnika, 1988, no. 1, pp. 49-53. (in Russian)
3. Aliyev A.V., AMarantov G.N. Vnutrennyaya ballistika raketnykh dvigateley na tverdom toplive [Internal Ballistics of Solid Propellant Rocket Engines]. Moscow, Mashinostroenie, 2007. (in Russian)
4. Godunov, S.K., Zabrodin A.V., Ivanov M.Ya., Kraiko A.N., Prokopov G.P. Chislennoe reshenie mnogomernyh zadach gazovoj dinamiki [Numerical Solution of Multidimensional Problems of Gas Dynamics], Moscow, Nauka, 1976. (in Russian)
5. Lipanov A.M., Dadikina S.Yu., Shumikhin A.A., Koroleva M.R., Karpov A.I. Numerical Simulation Intra-Chamber of Unsteady Turbulent Flows Stimulate. Part 1. Bulletin of the South Ural State University. Series: Mathematical Modelling, Programming and Computer Software, 2019, vol. 12, no. 1, pp. 32-43. (in Russian) DOI: 10.14529/mmp190103
6. Tadmor Z., Gogos C.G. Principles of Polymer Processing. New York, John Wiley and Sons, 1979.
7. Lipanov A.M. Teoreticheskaya gidromekhanika n'yutonovskih sred [Theoretical Hydromechanics of Newtonian Media]. Moscow, Nauka, 2011. (in Russian)
8. Kommisarov Yu.A., Gordeev L.S., Vent D.P. Processy i apparaty himicheskoj tekhnologii [Processes and Apparatuses of Chemical Technology. Part 1]. Moscow, URAIT, 2018. (in Russian)
9. Lykov A.V. Teplomassoobmen [Heat and Mass Transfer]. Moscow, Energiya, 1971. (in Russian)
10. Shumikhin A.A., Dadikina S.Yu. Numerical Simulation of a Compressible Viscous Gas Flow in Solid-Fuel Rocket Engine with a Central Body. Chemical Physics and Mesoscopy, 2020, vol. 22, no. 2, pp. 184-196. (in Russian) DOI: 10.15350/17270529.2020.2.18
11. Pino M.M., Piomelli U., Candler G.V. Subgrid-Scale Models for Compressible Large-Eddy Simulations. Theoretical and Computational Fluid Dynamics, 2000, vol. 13, no. 5, pp. 361-376.
12. Volkov K.N. [Simulation of Large Eddies of a Fully Developed Turbulent Flow in a Channel and Comparison of Subgrid Eddy Viscosity Models]. Journal of Applied Mechanics and Technical Physics, 2006, vol. 47, no. 3, pp. 31-42. (in Russian)
13. 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
14. 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
15. Shumikhin A.A. Numerical Simulation of a Compressible Viscous Gas Flow Using Parallel Calculations. Chemical Physics and Mesoscopy, 2021, vol. 23, no. 3, pp. 292-302. DOI: 10.15350/17270529.2021.3.26
Received April 25, 2022