CC BY
7
Механика
УДК 519.63 DOI: 10.14529/mmph230205
ПОЛУАНАЛИТИЧЕСКИЙ МЕТОД РЕШЕНИЯ УРАВНЕНИЙ ГАЗОВОЙ ДИНАМИКИ В ПЕРЕМЕННЫХ ЭЙЛЕРА
М.С. Жарылканова, Н.Л. Клиначева, А.П. Яловец
Южно-Уральский государственный университет, г. Челябинск, Российская Федерация E-mail: [email protected]
Аннотация. Представлен полуаналитический метод решения системы уравнений газовой динамики в эйлеровых координатах. Поскольку конечными разностями заменяются только пространственные производные, система уравнений газовой динамики сводится к системе обыкновенных дифференциальных уравнений на пространственной сетке. Приближенное аналитическое решение данной системы дифференциальных уравнений для малого временного интервала применяется для описания динамики газа во всем требуемом временном интервале. Проведена верификация на одномерных тестовых задачах о распаде произвольного разрыва и распространении стационарных ударных волн различной интенсивности. Для сравнения одномерных задач приведено решение тестовых задач простым в реализации базовым методом крупных частиц. Показано, что полуаналитический метод обладает высокой точностью вычислений, а также является наиболее универсальным методом для расчета прикладных задач.
Ключевые слова: полуаналитический метод; метод крупных частиц; ударная волна; распад произвольного разрыва.
Введение
В работе [1] описан метод решения системы уравнений механики сплошной среды, записанной в лагранжевых переменных, который отличается от традиционных методов тем, что конечными разностями заменяются только производные по пространственным переменным. В результате такого подхода уравнения в частных производных сводятся к системе обыкновенных дифференциальных уравнений на пространственной сетке. Для малого временного интервала было найдено приближенное аналитическое решение этой системы уравнений. Полученное аналитическое решение дает возможность описать динамику рассматриваемой системы во всем требуемом временном интервале. Применение данного метода для решения задач газовой динамики и расчета упругопластических течений в твердых телах показало, что данный метод обеспечивает высокую точность выполнения законов сохранения и характеризуется высокой устойчивостью решения, что обусловлено отказом от аппроксимации конечными разностями производных по времени.
Поскольку применение лагранжевых переменных для решения задач газовой динамики ограничено характером течений рассматриваемой системы, то возникает необходимость обобщения метода [1] на случай решения задач в эйлеровых переменных. В данной работе приводится описание применения метода [1] к решению задач газовой динамики в эйлеровых переменных, тестирование его на стандартных задачах путем сравнения результатов расчета предлагаемым методом с решениями, полученными базовым методом крупных частиц (МКЧ) [2], различные модификации которого широко используется для решения многих прикладных задач [3, 4].
Описание метода для эйлеровых координат
Рассмотрим систему уравнений газовой динамики в эйлеровых координатах в одномерном случае (законы сохранения массы, импульса и энергии):
dp d(pv)
dt dx
■ = 0, (1)
d(pv) | dp) dP^ (2)
dt dx dx'
Жарылканова М.С., Полуаналитический метод решения уравнений
Клиначева Н.Л., Яловец А.П. газовой динамики в переменных Эйлера
д(рЦ)]д(рЦу) (3)
Ы дх дх
Здесь и далее р - плотность, V - массовая скорость, и - внутренняя энергия единицы массы, Р - давление. Данную систему уравнений следует дополнить начальными и граничными условиями.
Чтобы данная система уравнений обеспечивала выполнение второго закона термодинамики, в ней применяется неравновесное давление, которое учитывает конечное время релаксации физической системы к равновесию [5]. Следствием конечного времени релаксации к равновесному состоянию является то, что при уменьшении объема газа с конечной скоростью он в каждый момент времени не находится в равновесии, и тогда неравновесное давление может быть представлено в виде
Р = Р0 + 0Р, (4)
где Р0 — равновесное давление, которое находится через уравнение состояния Р0 = Р0 (р,и),
8Р = с ор — неравновесная добавка, обусловленная локальным изменением массовой плотности среды при сжатии, с — скорость звука. Локальное изменение плотности определяется как 8р = ртге1, где тге1 с- время релаксации к равновесию, Е, - характерный линейный размер рассматриваемого объема среды. Поскольку изменение плотности среды связано с изменением объема соотношением р!р = —УIV, то из сказанного выше следует, что при сжатии среды неравновесная добавка к давлению имеет вид при V < О:
ЗР = -рс£Г/Г. (5)
Адиабатическое расширение среды происходит за счет убыли внутренней энергии, причем максимальная скорость убыли внутренней энергии определяется свойствами самой среды и реализуется при равновесном процессе. Примером может служить разлет газа в пустоту [5], который является равновесным процессом (энтропия сохраняется). Таким образом, в случае расширения среды, когда V > О, величина 5Р = 0 .
Для численного решения уравнений (1)-(3) на занятую область наносится неподвижная сетка х1, х2,... хг,...х1. Массовая скорость среды vi определена в узлах, а центрах ячеек определяются термодинамические характеристики среды р1+ц2, Щ+1/2, Р+1/2, и+1/2.
Для построения численной схемы проинтегрируем уравнения (1) и (3) в интервале от х до х/+1, а уравнение движения (2) - в интервале от х/—1/2 до х/+1/2. Для малых интервалов интегрирования можно записать соотношения:
х/+1 х/+1 х1+1/2 т+1/2 = \ р^=Р/+1 /2^/+1 /2, (ти)1+1,2 = \ Рис1х,(| рvdx, х/ х/ х/—1/2 где Ец+1/2 = х1+1 — х1 есть объем ячейки с единичной поперечной площадью.
Сумма диагональных элементов тензора скоростей деформации определяет скорость относительного изменения объема и, соответственно, массовой плотности при сжатии или разрежении. В плоской геометрии разностный аналог этого выражения будет иметь вид:
(V*).
_ 5И = у,-+1 - У,- = Ум/2
1+1/2 ^ХЛ+1/2 ^й-1/2 ^г+1/2
Таким образом, имеем У1+] 2 = '',+1 — Уг • Данное выражение характеризует изменение объема массы в ячейке только за счет процессов разрежения или сжатия без потоков массы через границы ячейки.
В результате интегрирования уравнения (1) получим:
^ = —{^+1 — ^), (6) где потоки массы через границы ячеек вычисляются как в [2]:
J, = v •
fPi—1/2, Vi > 0
LPi+1/2, Vi < 0,
Интегрирование уравнения (2) дает:
Ji+1 " Vi+1'
f Pi+1/2, Vi+1 > 0 lPi+3/2, Vi+1 < 0'
8и 8mi , ч , ч
Щ~д1 V = Р+Ш _р_1/2_/г"1/2 ^'
где потоки импульса через границы хг_1/2 и хг+1/2 вычисляются аналогично потокам массы:
Уг_1/2 > 0 \рм, Уг+1/2 > О
(7)
'i—1/2
= V
i—1/2
IPVi,
V,
i—1/2
< 0
, 'i+1/2
= V
i+1/2
Р/+1^+1> V+1/2 <0'
V±1/2 = 0,5(v±1 + vi ), Pi =
m,
i+1/2
m
i—1/2
Xi+1 X
, m = 0,5( m+1/2 + mi—1/2 )•
i—1
Результат интегрирования уравнения (3) может быть представлен в виде
d( mU )
i+1/2
dt
= —P+1/2 (V+1 — Vi ) —(3i+1 — 3i ),
(8)
где потоки внутренней энергии вычисляются аналогично потокам массы или импульса:
|Pi—1/2^—vi > 0
= V,
^i = vi 1 tt /л, "°i+1
[Pi+1/ 2 i+1 / 2, Vi < 0
Pi+1/2Ui+1/2, Vi +1 > 0
i+1
[Р+3/2^; +3/2' Vi+1 < 0
Уравнения (6)-(8) представляют собой систему обыкновенных дифференциальных уравнений на всей пространственной сетке. Будем интегрировать эту систему на некотором временном интервале от tn до .
Полагая, что в интервале tn < t < tn+1 потоки массы через границы ячеек постоянны, решение уравнения (6) и новая плотность запишется в виде
m.
i+1/2 (tn+1 ) = m+1/2 (tn ) + (J+1 — Ji ) Лt, P+1/2 (tn+1) = m+1/2 (tn+1 )/£+1/2-
(9)
Чтобы получить систему дифференциальных уравнений для скорости, продифференцируем уравнение (7) по времени, полагая, что потоки массы, импульса и неравновесная добавка (5) постоянны в рассматриваемом временном интервале, то есть
85Р^
d2m,
dt2
1 = 0 d'i+1/2
d t
= 0,
P+1/2
d t
= 0.
Выразим также частную производную от давления по времени через субстанциональную производную:
°ri+l/2 _ р0 i+1 / 2
dt
f дроЛ
dx
i+1/2
Полную производную по времени в (10) представим в виде
ро ri+1/2
(10)
(11)
где бг+1/2 = Р,+\ 2 с/+1 2 / ^1+] 2 ^ ^1+] 2 - эффективный объем ячейки, который учитывает, что при ее сжатии возмущается лишь слой. Таким образом, следуя [1], эффективный объем ячейки будет определяться выражением
£
i+1/2
[ с At, vi+1 - V; < 0
U-+1/2» Vi+1 - vi ^ 0
С учетом сделанных замечаний и выражений (10) и (11) получим
d 2v,
dv- ->
Я2" + 2Y "ST + = Li-dt dt
(12)
где
Y =
1 dmi m dt
= — (Qi+1/2 m
Qi—1/2 ), Li = —
m.
(0+1/2Vi+1 + Й—1/2^—1 ) + (V f)+1/2 " (V fx ^
i—1/2
Жарылканова М.С., Полуаналитический метод решения уравнений
Клиначева Н.Л., Яловец А.П. газовой динамики в переменных Эйлера
Для построения разностной схемы для вычисления входящих в последнее выражение производных удобно представить vдP / дх = д{уР') / дх — Pдv / дх. Тогда имеем
( дР^ 1
I V дР = -[ч+1 (Р+3/2 — Р+1/2 ) + V (Р+1/2 — Р—1/2 )]
V дх у/+1/2 2Ьг+1/2
дР 1
I V — = —-[V (Р+1/2 —Р—1/2) + V—1 (Р—1/2 — Р—3/2 )].
V дх у г—1/2 2Ьг—1/2
Здесь также уместно привести разностную аппроксимацию для вычисления неравновесной добавки к давлению:
яр - нср)/+1/2 (4+1 — V), 4+1 — V < 0 ОР+1/2 = 1 л • [ 0, vt+1 — V/ > 0
Таким образом, выражение (12) представляет собой систему дифференциальных уравнений, решение которой позволит найти поле скоростей во всех узлах эйлеровой сетки. Для решения
(12) необходимо задать начальные условия: vt (/п ) = V" и (д^- / дt)", которое находится из уравнения (7), где все величины в момент времени /и должны быть известны. Подставляя (6) в (7), можно записать второе начальное условие в виде
( ду Л" 1
Ь- У =— {(^+1/2 — ^—1/2 ) + (^+1/2 —1/2)}" ,
где
_ Г 0 Vt+1 /2 > 0 _ Гр—1СVt—1 — V )Vt—1/2, V—1/2 > 0
^+1/2 "1 / ч «,^—1/2 г,
+1 — V X+1/2, V+1/2 < 0 I 0 V—1/2 ^ 0
Для малого временного шага можно получить простое аналитическое решение системы (12). Полагая в (12) найденные для момента t" величины у., с2, Ц постоянными на шаге интегрирования, можно решение (12) записать в виде
£(Ш/) ( Ц
V ^ ) = V" + (£ У
"
а |
— V
1 — е
-уЫ.
С(АДГ) + у ?
а
(13)
где
гл2 2 2
а =с — у
ц—т , ч ео8(аAt), а2 > 0 , ч 8т(аAt), а2> 0 , а=ла2 ,с(а^Н /( t), (а^Н ( ', .
» (си(аAt), а2 < 0 [(ш/), а2< 0
Для решения уравнения (8) воспользуемся методом разделения по физическим процессам, для чего представим искомое решение в виде
и+1/2 ^) = и"+1/2 + 0иР+1/2 + ои—, (14)
где Ои^/2 — приращение внутренней энергии за счёт работы сил давления, 5и°0\12 — приращение внутренней энергии за счёт конвекции. Уравнение для Ои/+1/2 запишем в виде
д 1
■77^1/2 =--^-(15)
а <1/2
Отметим, что в уравнении (15) не учитываются потоки массы и энергии, что соответствует описанию в переменных Лагранжа. Подставляя выражение (14) в уравнение (8) и учитывая (15),
с»х тсот
запишем уравнение для Оиг+1/2 в виде
д(т5исо™) Гдт л -^^ = —(3/+1 —Д. ) — ^^] и" 1/2. (16)
Уравнение (16) записано в переменных Эйлера.
Начальные условия для уравнений (15), (16) имеют вид ОиР1/2 (/п ) = 8и°0\12 (/п ) = 0.
Интегрируя по времени выражение (11) от ^ до t < и учитывая (4), запишем выражение
Р+1/2 (t )=Р+1/2 (tn ) + ^Р+1/2_й+1/2 (V (t)_ V {гп )).+1/2. (17)
Подставляя (17) в (15), получим после интегрирования выражение
p 1
Щ+1/2
Ш,
—P (tn )ЛУ + Q
ЛУ 2
, (18)
i+1/2
/+1/2
где АГ1+1/2=Г1+1/2-М.
Интегрирование уравнения (16) с учетом (6) приводит к выражению
Ш™ = ^ДП {[р+3/2 (Ц+3/2 _ Ц+1/2 ) У+1 ]^<0 _[р_1/2 (и_1^2 _ и+1/2 ) V ]^>0}" . (19)
Таким образом, здесь приведен полный набор формул для описания динамики газа или жидкости в одномерной геометрии. Выражения (9) с приведенными в (6) потоками позволяют определить массы и плотности вещества в ячейках в новый момент времени. Выражение (13) позволяет вычислить новое поле скоростей, а выражения (14), (17) и (18) - внутренние энергии в ячейках. Новые равновесные давления и температура находятся с помощью уравнения состояния.
Основным приближением, которое было сделано при интегрировании системы (1)-(3), является постоянство бг+1/2 на шаге интегрирования Дt. Это приближение выполняется, если на шаге интегрирования относительное изменение объема мало, то есть
Л t
vi+1
+1/2
где £н _ некоторая малая величина. Данное условие позволяет найти величину шага интегрирования рассматриваемой системы. В случае состояний среды, близких к покою, или в случае гладких течений сформулированный выше критерий для выбора временного шага работать не будет. Поэтому для выбора временного шага при интегрировании уравнений применяется комплексный критерий, включая условие Куранта:
Л = min
V Vi +1 — VI c у
(20)
где £с ~ 1/3 _характерная константа в условии Куранта, £н ~0,01. Условие (20) проверяется для всех ячеек. Обычно оба условия дают близкие значения временных шагов, однако в области распространения ударных волн первое условие будет определяющим.
Точность интегрирования рассматриваемой системы можно повысить, решение (13) осуществить в два этапа. На первом этапе по формуле (13) для момента времени tn+1/2 = tn + Д /2 находится поле скоростей +1/2, Ц+12, а также потоки массы и энергии через границы ячеек. На
втором этапе по формуле (13), в которой вместо Ц подставлено Ц+1/2 , находим уП . По формулам (9) и (18) вычисляются новая масса, плотность и внутренняя энергия через потоки для момента ^+1/2. Тестовые расчеты
Для верификации полуаналитического метода в эйлеровых координатах были решены одномерные задачи о распаде произвольного разрыва, о распространении стационарных ударных волн различной интенсивности и проведено сравнение полученного решения с аналитическим и решением, полученным базовым методом крупных частиц.
Задача 1. Распад произвольного разрыва. В некоторой области, ограниченной отрезком [0;1] (см), находится газ, в начальный момент времени (t = 0) разделенный контактной границей на две подобласти. Показатель адиабаты у= 5/3. Геометрия сетки: N = 800 равномерно распределенных точек вдоль оси х.
Начальные распределения параметров газа вдоль пространственной координаты х равны
Жарылканова М.С., Клиначева Н.Л., Яловец А.П.
Полуаналитический метод решения уравнений газовой динамики в переменных Эйлера
р = 3,87 кг/м3; 21 = 611,5 К; ^ = 0 м/с; Р1 = 1,4-106 Па; х<0,5 см; р0 = 2,58 кг/м3; 20 = 91,7 К; ио = 0 м/с; Р0 = 1,4-105 Па; х >0,5 см. На рис. 1 представлены распределения параметров на момент времени £ = 4,5 мкс. Решение, полученное полуаналитическим методом (ПА метод), совпадает с аналитическим решением [6], что говорит о высокой точности метода. Исходя из графиков, видно, что в момент времени, отличный от нуля, влево начинает распространяется волна разрежения, вправо - ударная волна. Решение, полученное численными методами, имеет «размытие» в области контактного разрыва. Решение, полученное базовым МКЧ, дает небольшой скачок параметров на фронте ударной волны.
Л -МКЧ ---Анапит. решение • П/а метод
Рис. 1. Распад произвольного разрыва. Распределения скорости, давления, безразмерной плотности и температуры на момент времени £ = 4,5 мкс
Задача 2. Распространение сильной ударной волны (УВ). В некоторой области, ограниченной отрезком [0;1] (см), находится газ (показатель адиабаты у= 5/3), начальные значения параметров которого равны:
р0 = 1,29 кг/м3, 20 = 300 К, и0 = 0 м/с, Р0 = 2,29-105 Па; 0 <х< 1 см. Через границу (х = 0) в расчетную область втекает стационарный поток с параметрами:
Р = 10Р,, Ц, 21, р.
Параметры и, 2, Р рассчитываются из соотношений на разрыве (соотношения Ренкина-Гюгонио). Геометрия сетки: N = 800 равномерно распределенных точек вдоль оси х.
На рис. 2 представлены распределения параметров падающей сильной ударной волны на момент времени £ = 3 с. На рис. 3 представлены распределения параметров отраженной ударной волны на момент времени £ = 9 с. Полуаналитический метод дает хорошее совпадение с точным решением как для падающей, так и для отраженной ударных волн, в отличие от МКЧ. Решение, полученное базовым методом крупных частиц, имеет сильные осцилляции при отражении ударной волны от жесткой стенки. Это говорит о том, что МКЧ нуждается в введении дополнительных слагаемых, обеспечивающих устойчивость решения, не всегда имеющих физическую обоснованность. Одним из способов обеспечения устойчивости решения МКЧ является введение искусственной вязкости [7], в которую входят эмпирические константы. Использование таких слагаемых делает метод не универсальным. В то же время полуаналитический метод не требует введения никаких искусственных добавок.
Задача 3. Распространение слабой УВ. В некоторой области, ограниченной отрезком [0;1] (см), находится газ (показатель адиабаты у= 5/3), начальные значения параметров которого равны:
р0 = 1,29 кг/м3, 20 = 300 К, и0 = 0 м/с, Р0 = 2,29-105 Па; 0 <х< 1 см. Через границу (х = 0) в расчетную область втекает стационарный поток с параметрами
Р = 1,5Р„ Ц, 21, р.
Параметры и1, 21, Р1 рассчитываются из соотношений на разрыве. Геометрия сетки: N = 800 равномерно распределенных точек вдоль оси х .
-t2
—мкч
---Аналит. решение
• ПУа метод
4 6 8 1
X, мм
Рис. 2. Распространение сильной ударной волны. Распределения безразмерной плотности и давления на момент времени t = 3 мкс
-МКЧ
- Аналит. решение П/а метод
4 6
X, мм
Рис. 3. Распространение сильной ударной волны. Распределения безразмерной плотности и давления на момент времени t = 9 мкс
На рис. 4 представлены распределения параметров падающей ударной волны на момент времени t = 6 мкс. На рис. 5 представлены распределения параметров отраженной ударной волны на момент времени t = 24 мкс.
-мкч
---Аналит. решение
• П/а метод
1.8 с? 1.6 1.4
1.2
-МКЧ ---Аналит. решение • П/а метод
-
1 1
хЮ
4 6
X, мм
10
2
0 2 4 6
X, мм
Рис. 4. Распространение слабой УВ. Распределения безразмерной плотности и давления на момент времени t = 6 мкс
а 5
В
с£4
3
4 6
X, мм
10
Рис. 5. Распространение слабой УВ. Распределения безразмерной плотности и давления на момент времени t = 24 мкс
Исходя из полученных графиков для слабой ударной волны, видно, что решение, полученное базовым методом крупных частиц, не совпадает с аналитикой и показывает сильные осцилляции, то есть метод не устойчив при описании распространения слабых ударный волн. В свою очередь решение, полученное рассматриваемым методом, полностью совпадает с аналитическим, что говорит об универсальной применяемости данного метода для решения различных прикладных задач.
Заключение
В данной работе представлен полуаналитический метод в эйлеровых координатах для решения задач гидрогазодинамики. Представлены результаты верификации данного метода на одномерных задачах и сравнение с другим численным методом. Предложенный в данной работе метод показывает высокую точность для ударных волн различной интенсивности, что делает метод универсальным и применяемым для различных прикладных задач.
Данный метод применим для расчетов в многомерном случае, поскольку в этом случае описанный здесь алгоритм решения задачи сохраняется полностью.
Литература
1. Яловец, А.П. Расчет течений среды при воздействии интенсивных потоков заряженных частиц / А.П. Яловец // Прикладная механика и техническая физика. - 1997. - Т. 38, № 1. -С.151-166.
Жарылканова М.С., Полуаналитический метод решения уравнений
Клиначева Н.Л., Яловец А.П. газовой динамики в переменных Эйлера
2. Белоцерковский, О.М. Численное моделирование в механике сплошных сред / О.М. Бело-церковский. - М.: Физматлит, 1994. - С. 27-39.
3. Shestakovskaya, E.S. On one Method of Calculating Moving Boundaries in Euler Coordinates / E.S. Shestakovskaya, Ya.E. Starikov // Journal of Computational and Engineering Mathematics. - 2019. - T. 6, № 4. - P. 44-56.
4. Беляев, П.Е. Влияние экранирующего слоя газовзвеси на силовое воздействие ударной волны на жёсткую стенку / П.Е. Беляев, Н.Л. Клиначева // Вестник ЮУрГУ. Серия «Математика. Механика. Физика». - 2016. - Т. 8, № 4. - С. 49-55.
5. Ландау, Л.Д. Гидродинамика / Л.Д. Ландау, Е.М. Лифшиц. - М.: Наука, Главная редакция физико-математической литературы, 1988. - 733 с.
6. Куропатенко, В.Ф. Основы численных методов механики сплошной среды: монография / В.Ф. Куропатенко, Е.С. Шестаковская. - Челябинск: Издат. центр Южно-Уральского государственного университета, 2017. - 253 с.
7. Садин, Д.В. Модификация метода крупных частиц до схемы второго порядка точности по пространству и времени для ударно-волновых течений газовзвеси / Д.В. Садин // Вестник ЮУрГУ. Серия «Математическое моделирование и программирование». - 2019. - Т. 12, № 2. - С. 112122.
Поступила в редакцию 27 марта 2023 г.
Сведения об авторах
Жарылканова Мадина Салимжановна - студент, Южно-Уральский государственный университет, г. Челябинск, Российская Федерация, e-mail: [email protected]
Клиначева Наталия Леонидовна - кандидат физико-математических наук, доцент, кафедра «Вычислительная механика», Южно-Уральский государственный университет, г. Челябинск, Российская Федерация, e-mail: [email protected]
Яловец Александр Павлович - доктор физико-математических наук, профессор, кафедра «Вычислительная механика», Южно-Уральский государственный университет, г. Челябинск, Российская Федерация, e-mail: [email protected]
Bulletin of the South Ural State University Series "Mathematics. Mechanics. Physics" _2023, vol. 15, no. 2, pp. 32-40
DOI: 10.14529/mmph230205
SEMIANALYTIC METHOD FOR SOLVING GAS DYNAMICS EQUATIONS IN EULER VARIABLES
M.S. Zharylkanova, N.L. Klinacheva, A.P. Yalovets
South Ural State University, Chelyabinsk, Russian Federation E-mail: [email protected]
Abstract. This paper presents a semi-analytical method for solving a system of equations of gas dynamics in Eulerian coordinates. Since only spatial derivatives are replaced by finite differences, the system of gas dynamic equations is reduced to a system of ordinary differential equations on a spatial grid. An approximate analytical solution of this system of differential equations for a small time-interval is used to describe the dynamics of a gas in the entire required time interval. Verification was carried out on one-dimensional test problems on the decay of an arbitrary discontinuity and the propagation of stationary shock waves of various intensities. To compare one-dimensional problems, the solution of test problems is given by the simple-to-implement basic particle-in-cell method. It is shown that the semi-analytical method has high accuracy of calculations, and is also the most universal method for calculating applied problems.
Keywords: semi-analytical method; particle-in-cell method; shock wave; decay of an arbitrary discontinuity.
References
1. Yalovets A.P. Raschet techeniy sredy pri vozdeystvii intensivnykh potokov zaryazhennykh chastits (Calculation of Medium Flows under the Influence of Intense Flows of Charged Particles). Prikladnaya mekhanika i tekhnicheskayafizika, 1997, Vol. 38, no. 1, pp. 151-166. (in Russ.).
2. Belotserkovskiy O.M. Chislennoe modelirovanie v mekhanike sploshnykh sred (Numerical Modeling in Continuum Mechanics). Moscow, Fizmatlit Publ., 1994, pp. 27-39. (in Russ.).
3. Shestakovskaya E.S., Starikov Ya.E. On One Method of Calculating Moving Boundaries in Euler Coordinates. Journal of Computational and Engineering Mathematics, 2019, Vol. 6, no. 4, pp. 44-56. DOI: 10.14529/jcem190404
4. Belyaev P.E., Klinacheva N.L. Impact of Gas Suspension Shielding Layer on the Force Effect of Shock Waves on a Rigid Wall. Bulletin of the South Ural State University. Series of "Mathematics. Mechanics. Physics", 2016, Vol. 8, no. 4, pp. 49-55. (in Russ.). DOI:10.14529/mmph160406
5. Landau L.D., Lifshits E.M. Gidrodinamika (Hydrodynamics), Moscow: Nauka Publ., 1988, 733 p. (in Russ.).
6. Kuropatenko V.F., Shestakovskaya E.S. Osnovy chislennykh metodov mekhaniki sploshnoy sredy: monografiya (Fundamentals of Numerical Methods of Continuum Mechanics: Monograph), Chelyabinsk, Izdat. tsentr Yuzhno-Ural'skogo gosudarstvennogo universiteta Publ., 2017, 253 p. (in Russ.).
7. Sadin D.V. A Modification of the Large-Particle Method to a Scheme Having the Second Order of Accuracy in Space and Time for Shockwave Flows in a Gas Suspension. Bulletin of the South Ural State University. Series "Mathematical Modelling, Programming and Computer Software", 2019, Vol. 12, no. 2, pp. 112-122. DOI: 10.14529/mmp190209
Received March 27, 2023
Information about the authors
Zharylkanova Madina Salimzhanovna is Student, South Ural State University, Chelyabinsk, Russian Federation, e-mail: [email protected]
Klinacheva Nataliya Leonidovna is Cand. Sc. (Physics and Mathematics), Associate Professor, Department of "Computational Mechanics", South Ural State University, Chelyabinsk, Russian Federation, e-mail: [email protected]
Yalovets Aleksandr Pavlovich is Dr. Sc. (Physics and Mathematics), Professor, Department of "Computational Mechanics", South Ural State University, Chelyabinsk, Russian Federation, e-mail: [email protected]
