О ПРИМЕНЕНИИ WENO-СХЕМ К МОДЕЛИРОВАНИЮ РЕАГИРУЮЩИХ ГАЗОВЫХ ПОТОКОВ

ПОТАПКИНА Ю. Ю., ПЕСКОВА Е. Е.1 О ПРИМЕНЕНИИ WENO-СХЕМ К МОДЕЛИРОВАНИЮ РЕАГИРУЮЩИХ ГАЗОВЫХ ПОТОКОВ

Аннотация. В работе исследованы WENO-схемы 5 и 7 порядка с оптимальными весовыми коэффициентами без учета индикаторов гладкости решения и с индикаторами гладкости. Показано, что при их применении к решению задач многокомпонентной газовой динамики с химическими реакциями, диффузией, вязкостью и теплопроводностью WENO-схема 5 порядка с оптимальными весовыми коэффициентами дает более точный результат при меньшей трудоемкости вычислений.

Ключевые слова: математическое моделирование, уравнения Навье-Стокса, WENO-схемы, дозвуковые реагирующие потоки.

POTAPKINA YU. YU., PESKOVA E. E.

ON THE APPLICATION OF WENO SCHEMES TO MODELING REACTING GAS FLOWS

Abstract. The 5th and 7th order WENO schemes with optimal weight coefficients without solution's smoothness indicators and with solution's smoothness indicators are investigated in the article. It is shown that the 5th order WENO scheme with optimal weight coefficients gives a more accurate result with less complexity of calculations to solve problems of multicomponent gas dynamics with chemical reactions, diffusion, viscosity and thermal conductivity.

Keywords: mathematical modeling, Navier-Stokes equations, WENO schemes, subsonic reactive flows.

Введение. В настоящее время широкое распространение получило исследование газодинамических течений в задачах химической промышленности, поскольку технологам для получения целевых продуктов необходимо знать большое количество параметров проведения реакции. Наиболее значимым методом исследования стало математическое моделирование, которое позволяет рассмотреть поведение различных реакций в разных условиях без проведения лабораторных экспериментов. Аппарат математического моделирования расширяется и улучшается в прямой зависимости от появления новых задач: создаются новые математические модели, эффективные вычислительные алгоритмы, повышается точность расчетов. Настоящая работа направлена на исследование применения схем высокого порядка точности (WENO-схема 5 и 7 порядка) к решению задач

1 Исследование Песковой Е. Е. выполнено за счет гранта Российского научного фонда № 23-21-00202.

многокомпонентной газовой динамики с учетом вязкости, диффузии, теплопроводности и химических реакций.

Математическая модель и численный алгоритм. Поскольку нас интересует исследование возможности повышения порядка точности расчетов с помощью WENO-схем, в настоящей работе мы рассматриваем одномерную модель уравнений вязкой дозвуковой динамики химически активной газовой смеси:

ди д

Векторы и, Р(и), Н(и), W(U) имеют следующий вид:

!рУш\ л (риуЛ АД /Яя

и = ( ри ),Р(и) = ( ри2 ),Н(и) = ( ТХХ )^(и) = ( 0 РЧ \phuj \Чх) ^ 0

Р¥™\ , . ¡Ри¥т\ , ч /1тх\ ч /К

I  т I . . I I Jтх . „ ч

и =

= \ 0 ).

р^ \ рЬи / \Чх/ ^ 0

В этой системе уравнений т = 1, ...,М, М - число компонент в газовой смеси, р - плотность смеси, Ут- массовая доля -ой компоненты смеси, и - скорость, И - энтальпия смеси, ]тх -диффузионный поток, Ит - скорость образования или расхода т-ой компоненты смеси, п = р — р0 - динамическая составляющая давления, р0 - термодинамическая составляющая давления, постоянная в области, тхх - вязкий поток, qх - потока тепла.

Система дополняется условием на дивергенцию вектора скорости:

V • и! = —

рсрТ \ I р М\тп

т

р^\Мыт СрТ '

Здесь Ср - теплоемкость смеси при постоянном давлении, Т - температура смеси, А -теплопроводность смеси, Отт1х - коэффициент диффузии, hт - энтальпия образования компоненты смеси, Мw - молекулярная масса смеси, Мwт - молекулярная масса компоненты смеси. Более подробно математическая модель представлена в работе [1].

Для построения вычислительного алгоритма используем равномерную сетку отрезков:

&Ах = = 1, .--,N^1= [х1-1,х1],1^1 = х{ —х{-1 = hх,hхNх = Ьх}

Численный алгоритм строится по схеме расщепления по физическим процессам [2]. На первом этапе в этом расщеплении решается система уравнений химической кинетики, на втором этапе интегрируются уравнения законов сохранения без учета давления, далее рассчитывается поле поправок к давлению из решения уравнения Пуассона, на последнем

этапе корректируются поле давления и поле скорости [2]. Поскольку мы исследуем применение WENO-схем, приведем подробно алгоритм второго этапа. Для интегрирования законов сохранения используется разностная схема вида:

ЦП+1 _ ЦП Р1+1/2 _ ¥1-1/2 Щ + 1/2 _ Щ-1/2 _

--\----— и.

Д t И.х кх

Здесь Н1+1/2, Н1-1/2 - диффузионные, вязкие и тепловые потоки на границе ¿и1 + 1, ¿_1и1 ячейками соответственно, которые рассчитываются по схеме с центральными разностями, ^1+1/2, ^\-1/2 - конвективные потоки на границе ¿и1 + 1, ¿_1и1 ячейками соответственно, которые рассчитываются с использованием потоков Русанова [3]:

¥+х/2 — и.5 (р(Ц[+1/2) + Р(и}+1/2) _ а(Щ+1/2 _ Ц}+1/2)), а — max(\u¡+l/2\,\uli+l/2\), где Щ+1/2 и Ц1+1/2 - значения вектора переменных и справа и слева от границы между г и г+1 ячейками, которые рассчитываем с использованием WENO схем 5 и 7 порядка точности [4; 5]. Их суть заключается в следующем. Для нахождения значения вектора переменных и справа и слева от границы между г и г+1 используется выражение:

к

(V) +1/2.

У=1

к

Ц+1/2 — ^ Ц(

Здесь К — 3 для случая схемы 5-го порядка, К — 4 для схемы 7-го порядка точности, П(у) -оптимальные весовые коэффициенты, Ц¿^^ - значение на границе, полученное на шаблоне 5 — {х1-у, ...,Х1,..., х£_у+к-1}:

к-1

ц

1 + 1/2 — / ,'Ури1_У+р.

р = 0

В случае использования весовых коэффициентов с индикаторами гладкости значения вектора переменных и справа и слева от границы между г и г+1 находятся из выражения:

К №

(V) —

Щ+1/2 — ^ ^ ЦИ1/2^ — уК^Гу*

УК=л°(уу [е + и^ЦР'

у=1 у 1

Здесь IЯ(у) - это индикаторы гладкости, £ - некоторое малое число, вводимое, чтоб предотвратить деление на ноль. Величины сУр, 1для К — 3, К — 4 приведены в

работах [4; 5].

Вычислительные эксперименты. Рассмотрим следующую одномерную постановку задачи. Размер области 0.2 м, шаг по пространству 2 • 10_3 м, шаг по времени 1 • 10_5 с, расчет ведем до 0.2 секунд. Принимаются следующие начальные данные: в области от 0 м до 0.1 м

температура газа 800оС; в области от 0.1 м до 0.2 м температура газа 1200оС; скорость потока 0.1 м/с; давление 101325Па; состав газовой смеси - метан 100% (СН4). Граничные условия: на границе слева со скоростью 0.1 м/с втекает метан с температурой 8000С, на границе справа задаются условия вытекания. Такая постановка задачи принята для того, чтобы выяснить какая из схем будет лучше считать в области перепада температур. В этой же области начнут происходить интенсивные реакции с расходом метана и образования продуктов реакции.

Рассматриваем 3 варианта расчета значений вектора и на границе ячеек: 1) схемой первого порядка точности, т.е. и'[+1/2 = и1+1, Щ+х/2 = ^; 2) WENO схемой 5-го порядка аппроксимации; 3) WENO-схемой 7-го порядка аппроксимации. На рисунках 1 и 2 представлены распределения температуры и метана. Из графиков можно сделать вывод, что схема первого порядка аппроксимации сглаживает решение в областях резкого изменения газодинамических параметров, разница температур составляет величину 100С. Такой результат является достаточно большим расхождением в случае проведения лабораторных экспериментов. Расчет по схемам WENO 5-го и 7-го порядка практически совпадает, что говорит о преимуществе схемы WENO5, поскольку она использует более компактный шаблон. На рисунках 9 и 10 представлены расчеты по WENO-схемам с оптимальными весами и с расчетом индикаторов гладкости. Из графиков можно сделать вывод, что для рассматриваемых задач использование WENO-схем с оптимальными весами является преимущественным, поскольку графики полностью совпадают, а они являются менее трудоемкими. Полученный результат полного совпадения графиков можно объяснить преобладанием диффузионных процессов над конвективным переносом.

1000 980 960

—Т. 1 order -Т. WEN05 -T.WEN07

940

920 900

880-

860 840

820

8ÜU

7800 5 То 15 20 25 30 35 40 45 50

х,т

Рис. 1. Распределение температуры. Расчет по схемам 1-го порядка (синяя линия), WENO5 (красная линия), WENO7 (зеленая линия).

Рис. 2. Распределение метана. Расчет по схемам 1-го порядка (синяя линия), WENO5 (красная линия), WENO7 (зеленая линия).

Рис. 3. Распределение температуры. Расчет по схеме WENO5 (красная линия) и схеме WENO5 с оптимальными весами (бирюзовый пунктир).

Рис. 4. Распределение температуры. Расчет по схеме WENO7 (зеленая линия) и схеме WENO7 с оптимальными весами (синий пунктир).

Выводы. В работе проведено исследование применения WENO-схем 5-го и 7-го порядка точности с оптимальными весовыми коэффициентами и с весовыми коэффициентами, рассчитанными с использованием индикаторов гладкости к задачам многокомпонентных

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

СПИСОК ЛИТЕРАТУРЫ

1. Жалнин Р. В., Пескова Е. Е., Стадниченко О. А., Тишкин В. Ф. Моделирование течения многокомпонентного химически активного газа на примере пиролиза углеводородов // Препринты ИПМ им. М.В. Келдыша. - 2017. - № 101. - 16 с.

2. Пескова Е. Е., Снытников В. Н. Численное исследование конверсии метановых смесей под воздействием лазерного излучения // Журнал Средневолжского математического общества. - 2023. - Т. 25, № 3. - С. 159-173.

3. Русанов В. В. Расчет взаимодействия нестационарных ударных волн с препятствиями // Журнал вычислительной математики и математической физики. -1961. - Т. 1, № 2. - С. 267-279.

4. Shu C.-W. Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws // Advanced Numerical Approximation of Nonlinear Hyperbolic Equations. - 2006. - Vol. 1697. - P. 325-432.

5. Евстигнеев Н. М. О построении и свойствах WENO-схем пятого, седьмого, девятого, одиннадцатого и тринадцатого порядков. Часть 1. Построение и устойчивость // Компьютерные исследования и моделирование. - 2016. - Т. 8, № 5. - С. 721-753.