Челябинский физико-математический журнал. 2020. Т. 5, вып. 3. С. 316-326.
УДК 517.958 БОТ: 10.47475/2500-0101-2020-15306
ДИФФЕРЕНЦИАЛЬНО ИНВАРИАНТНОЕ РЕШЕНИЕ УРАВНЕНИЙ ДИНАМИКИ ДВУХФАЗНОЙ СРЕДЫ
А. В. Панов
Челябинский государственный университет, Челябинск, Россия gjd@bk.ru
Рассматривается система уравнений, описывающая динамику газовзвеси в изотермическом случае. Выведена подмодель однородных движений среды, инвариантных относительно галилеевых сдвигов по первой координате. Данная подмодель является дифференциально инвариантной подмоделью относительно четырёхмерной подалгебры с базисом из операторов сдвига по всем пространственным координатам и оператора преобразования Галилея по первой координате. Подмодель приведена в инволюцию и проинтегрирована, найдено точное решение рассматриваемой системы. Проведено сравнение полученного решения с известным частично инвариантным решением. Относительно этой же подалгебры из оптимальной системы подалгебр алгебры Ли симметрий уравнений динамики газовзвеси, порождающей однородные движения среды, инвариантные относительно галилеевых сдвигов по первой координате, построено групповое расслоение уравнений динамики газовзвеси, разбивающее исходную систему на две части: автоморфную и разрешающую. Данное разбиение содержит все решения исходной системы и наоборот.
Ключевые слова: группа симметрий, дифференциально инвариантное 'решение, двухфазная среда.
Введение
Рассматривается система уравнений в частных производных, описывающая динамику двухфазной среды,
—Р1 + и = 0,
аЬ1
—Р2 , л- п
— + P2d.IV и2 = 0, —и 2 (1)
Р1—¡^ + т{7Р(Р1,Р2) = - — (и1 - и2), аъ1 т
Р2 —¡Т2 + Ш2^Р(Р1,Р2) = — (¡1 - ¡2), аъ2 т
здесь и1 = (м1,^1,^1), и2 = (п2,ь2,т2) — векторы скоростей первой и второй фаз, Ръ Р2 — парциальные плотности фаз, Р(р1, р2) — давление смеси (функциональный
параметр системы), т2 = -Р-2 — объёмная доля второй фазы, р22 — абсолютная
Р22
плотность второй фазы (постоянная величина), т1 = 1 — т2 — объёмная доля
. — д _ — д _
первой фазы, — = — + и1 • V, — = + и2 • V. аЬ1 дЬ аЬ2 дí
Работа выполнена при поддержке гранта РНФ (проект № 18-71-00041).
Система (1) была предложена в работе [1], её исследовали и развивали разные авторы (см., например, [2; 3]). Исследования системы (1) были направлены на решение различных проблем: корректность задачи Коши, задача о распаде произвольного разрыва, расчёт распространения ударных волн, расчёт разлёта сжатого объёма и другие задачи. Большинство этих задач решались численными методами. В современных исследованиях начали рассматривать полидисперсные смеси с существенно различными размерами дисперсных компонент среды [4], что приводит к дополнительным трудностям в численных расчётах. Как ранее некорректность задачи Коши, так и сейчас проблемы, связанные с жёсткостью уравнений полидисперсной среды, приводят к трудностям в решении задач многофазной газодинамики численными методами. Поэтому поиск точных решений, которые могут быть использованы для апробации численных схем, остаётся актуальной проблемой.
Одним из самых эффективных методов поиска точных решений нелинейных уравнений в частных производных является подход на основе свойств симметрии изучаемой системы. Среди основных монографий по данному направлению можно указать работы [5-8].
Ранее был найден ряд точных решений системы (1) в одномерном и трёхмерном случаях [9; 10]. Это были инвариантные либо частично инвариантные решения. Поиск таких решений основан на упрощении исходной системы за счёт свойств её симметрии. В инвариантном случае упрощения могут быть довольно различными: к алгебраической системе, к системе обыкновенных дифференциальных уравнений, к системе уравнений в частных производных с меньшим числом независимых переменных. В случае частично инвариантных решений исходная система упрощается, однако редуцированная система разбивается на две части — инвариантную и «активную» системы. Последняя включает как параметры функции из инвариантной системы и имеет нетривиальные условия совместности. Такие упрощения можно делать, основываясь на любой подалгебре из алгебры Ли симметрий исследуемой системы. Однако сопряжённые подалгебры будут давать подмодели (упрощения), переводимые друг в друга преобразованием из группы симметрий. Поэтому вывод и анализ подмоделей необходимо делать для несопряжённых подалгебр. Таким образом возникает задача о построении оптимальной системы подалгебр данной алгебры Ли симметрий. Можно сказать, что в этой оптимальной системе, т.е. в полном наборе несопряжённых подалгебр, содержится алгебраическая «структура» множества точных решений системы.
Некоторые подалгебры не порождают инвариантные подмодели в силу специфики действия соответствующей подгруппы в заданном пространстве зависимых и независимых переменных. В этом случае можно искать дифференциально инвариантные решения [11], что кажется более простой задачей по сравнению с поиском частично инвариантных решений. В данной работе предлагается пример поиска точного дифференциально инвариантного решения системы (1).
1. Дифференциально инвариантное решение
Алгебра Ли группы симметрий системы (1) с произвольной функцией давления Р имеет базис [12]
х — д х — д х — д Х1 — 7Т" , х2 — 7Г" , х3 — 7Т" ,
дх ду дг
д д д ддд д д д Х4 = --+ ---+ --, Х5 = ¿---+ —--+ --, Х6 = ¿—--+ —--+ --,
дх дп\ ди2 ду дь\ дь2 дг дгш\ дгш2
д д д д д д
X7 = УTT - Z— + Vi---Wi---+ --
oz ду дw1 ovi 0w2 0v2
д д д д д д X8 = z—--x——+ w^---щ —--+ w^---u2-
дх дг ди1 дт1 ди2 дт2
д д д д д д д Хд — Х---у— + П\—--VI---+ П^ —--, Хю — —.
ду дх дvl дщ дv2 ди2 дъ
Эта алгебра является алгеброй Ли группы Галилея, действующей в пространстве зависимых и независимых переменных системы (1). Оптимальная система подалгебр алгебры Ли группы Галилея была получена, например, в работе [13].
Рассмотрим подалгебру из оптимальной системы с базисными операторами Х1, Х2, Х3, Х4. Если рассматривать действие группы Ли, порождённое данной подалгеброй, в пространстве переменных — (Ъ, х, у, г, р1, р2, и1, и2), то полным набором функционально независимых инвариантов являются следующие функции:
1о — {Ъ, Р1, Р2, П1 — П2, VI, V2, ~Ш1, т2}.
Инвариантного решения не построить, так как недостаточно инвариантов — «зависимых переменных». Найдём дифференциально инвариантное решение системы (1). Для этого расширим пространство, на котором действует группа, до пространства 1-струй 31 — (Ъ,х,у,г,р1,р2, иь и2,Ри,Р1х,Р1у ,Ри ,Р21,..., и21, и2х, и2у, и2г ) — 44-мерное пространство с координатами: независимые переменные, зависимые переменные, все первые производные зависимых переменных по независимым. Продолжим действие группы на пространство 31, используя стандартные формулы продолжения (см., например, в [5]). Смысл данного продолжения заключается в том, что переменные, отвечающие за производные, должны преобразовываться так, чтобы производные преобразованной функции /(Ъ,х,у,г) совпадали с преобразованными производными функции /(Ъ,х,у, г). Итак, используя стандартную процедуру, получим продолжение данных операторов на первые производные
х 1 — — х 1 — д. х 1 —
1 дх 2 ду' 3 дг^
, д д д д д д д Х4 — + ^--+ ^--П1х"7,--П2х"7,--Vlx^--V2x^--
дх ди1 ди2 дии диц дvlt дv2t
д д д д -wix---w2x^--Pix7;--P2x
' 2x ^ f-* ix ^ f-* 2X r\ '
0wu 0w2i др-t др21
Теперь полный набор функционально независимых инвариантов состоит из ранее полученных функций и дифференциальных инвариантов:
Ii = Io U {fx, fy, fz, gitg2x g2tglx, hixpit рlxhlt, h2xp2t p2xh2t},
где f принимает значения из множества {pi,p2,u1,v1,w1,u2,v2,w2}, g принимает значения из множества {p,u,v,w}, а h из множества {v,w}. Стандартная техника поиска инвариантных решений предписывает построение следующих анзацев: необходимо приравнять инварианты с зависимыми переменными (pi, p2, ui,...) к новым неизвестным функциям, которые зависят от инвариантов, содержаших только независимые переменные. В случае дифференциально инвариантного решения получим равенства
pi = ri(t), vi = ^i(t), wi = Xi(t), u- -u2 = (p(t), u-x = a(t), u-y = ß(t), u-z = Y(t),...
(2)
Выполнение условий совместности приведёт к тому, что все неизвестные функции в уравнениях (2), оставшихся в точках, будут определяться из представленных уравнений. К примеру, неизвестная функция &(£) в уравнении г1у = &(£) определяется в силу уравнения г1 = 01^), а из уравнения и — и2 = р(£) и последних выписанных уравнений в (2) следуют равенства и2х = а(£), и2у = в(¿), и2-г = 7(¿). Подставим в систему (1) выписанные в (2) равенства. Получим уравнения
Г 2 1
Г- + па = 0, =--(01 — 02), 02 = -(01 — 02),
Г1Т т
х1 = — —(Х1 — Х2), х2 = 1(Х1 — Х2), (3)
г1т т
Г2 +1
р' + ар + (01 — 02)в + (Х1 — Х2 Ь = —-^
т
составляющие инвариантную часть системы (1), также остались уравнения на функцию и2
Р
^ = «24 + «2«2х + 02«2у + Х2«2.г--= 0,
Т (4)
О1 = «2х — а = 0, О2 = «2у — в = 0, Оэ = «2^ — 7 = 0,
составляющие «активную» часть системы (1). Функции р, 02, Х2, а, в, 7 входят в систему (4) как параметры. Поэтому уравнения ^ = 0, Ог = 0, г € {1, 2, 3}, должны удовлетворять условиям совместности, которые выражаются равенствами нулю скобок Якоби Сг] = 0, [Ог, О^-] = 0 на решениях системы (4).
Скобкой Якоби [14] двух уравнений ^ (х1,... , жп, и, иХ1, ..., иХп) = 0, О(х1,...,жп,и,иХ1, ..., иХп) = 0 на функцию и(х1, ..., жп) называется выражение
п
О] = £ ((^ + ихгО^ — (Охг + ихгОи) ^) .
г=1
Подсчёт скобок Ог] = 0 вместе с уравнениями Ог = 0 даёт уравнения на параметры а, в, 7, при которых система (4) имеет решение:
а' + а2 = 0, в' + ав = 0, 7' + а7 = 0.
Решения этих уравнений а = , в = , 7 = , либо а = 0, в = с2, 7 = сэ. В первом случае решение системы (3) имеет вид
Г1 = —-Г, Г2 = —-г, 01 = С5 + р,С4в т ,02 = С5 — С4в т ,
+ С11 + С11
_М±1 У
_М±1 у _ — (ц + 1)(С2С4 + СэСб)£ + С8 _м±11
Х1 = С7 + цсев Т , Х2 = С7 — Сев Т , р =-—-.-е т
К + С11
Подставляя найденные функции, в том числе а, в, 7, в уравнение ^ = 0 и решая полученное обыкновенное дифференциальное уравнение, найдём, что
и2 = —-г —(С5С2 + С7Сэ)£ + (С2С4 + СэСе)^е т--— е т + /(ж,у, г) .
К + С1| \ ц + 1 /
Наконец, функция /(ж, у, г) находится из оставшихся уравнений Ог = 0. Она имеет вид /(ж, у, г) = sgn(í + С1)ж + С2У + Сэг.
Во втором случае, когда а = 0, в = с2, 7 = с3, аналогично получим Г1 = Г°, Г2 = ^Г°, 01 = С5 + ^С4в т , 02 = С5 - С4в т ,
/ / -.ч/ ч. ч +
Х1 = С7 + ^Свв т , Х2 = С7 - Сев т = (-(^ + 1)(С2С4 + СзСе)£ + Сз) е т ,
I / \ / \ М+1 + С8 м+1^ «2 = -(С5С2 + С7Сз)^ + (С2С4 + С3С6 )£е т--— е т + С2у + Сз^
V ^ +1
Если вернуться к операторам Х1,Х2,Х3,Х4 и пространству переменных с набором инвариантов 1°, то здесь возможен поиск частично инвариантного решения в виде
Рг = Г^), V = 0Д£), ОД = Хг(^), «1 - «2 =
где функция и2, или и1, выражающая дефект инвариантности, является функцией переменных ¿, х, у, Здесь нет дополнительных предположений о производных, как ранее в дифференциально инвариантном случае. Это усложняет анализ и решение полученной системы, но найденное частично инвариантное решение обладает большим произволом — оно содержит произвольную функцию, такое решение найдено ранее в работе [12].
2. Групповое расслоение
В приведённом примере осуществлена редукция системы (1) к системе обыкновенных дифференциальных уравнений. Таким образом выделено одно из возможных решений системы (1), зависящее от некоторого набора констант интегрирования. Используя факт симметричности системы (1) относительно преобразований, порожденных операторами Х1, Х2, Х3, Х4, можно выполнить редукцию системы (1) так, что полученная система будет эквивалентна исходной системе (1), т. е. все решения исходной системы будут решениями новой системы и наоборот. Новая система называется групповым расслоением системы (1).
Несколько слов о конструкции, подробности алгоритма и основные теоремы есть в [5]. Данная группа Ли преобразований разбивает множество решений рассматриваемой системы на классы эквивалентных: два решения эквивалентны, если одно можно перевести в другое групповым преобразованием. Класс эквивалентных решений называется орбитой действия группы на множестве решений. Задача группового расслоение состоит в следующем: для заданной системы и группы Ли О её симметрий найти уравнения АО, определяющие орбиту произвольной функции, и найти уравнения Яб", решения которых задают всю совокупность орбит решений системы. Система АО является автоморфной, то есть любое её решение получается из одного при помощи преобразований допускаемой группы Ли О. Уравнения называют разрешающей системой.
Так как автоморфная система является инвариантной относительно действия группы О, то строить её надо так же, как инвариантное решение. Однако, чтобы учесть все решения системы (1), необходимо задавать многообразие в пространстве инвариантов размерности 4, размерность графиков решений системы (1) равна 4. Поэтому вместо (2) здесь будет анзац вида
V! = ^(¿,Р1,Р2,и), ^2 = ^(¿,Р1,Р2,и), = Ж1(г,Р1,Р2,«), ^2 = (¿, Р1, Р2, «),
Ргт = ЯГ(£,Р1,Р2,И), «гт = Ц™^, Р1, Р2, «), (5)
^т = Угт(£,рЬр2,и), ^¿т = Р1, Р2, «),
где за и обозначен инвариант и 1 — и2, I — 1, 2, т € {х, у, г}. Остались ещё уравнения р 1 гР2х — Р21Р1X — К(Ъ,р 1,Р2,и), и 1 Щх — и2ги 1 х — ^(Ъ, р 1,Р2,и),
V 1^2х — V2tV 1 х — V(Ъ,р 1,Р2,и), -и-2х — W2tW 1х — W(Ъ, р 1,Р2,и),
(6)
V 1хР 11 — р 1хV11 — Р (Ъ,р 1,Р2,и), т 1 хРи — р 1 х-и — Р2(Ъ,р 1,Р2,и),
V2xР2t — р2х^2Ь — Рз(Ъ,р 1,р2,и), -2x^21 — р2х—21 — Р4(Ъ, р 1,р2,и).
Уравнения (5), (6) образуют АО систему (аналог системы (4)), в которой функции, стоящие в правых частях равенств, считаются параметрами. Они находятся предварительно из разрешающей системы КБ — аналог системы (3) и уравнений на функции а, в, 7 — которая состоит из условий совместности АО системы и проекции системы (1) в пространство инвариантов (пространство переменных Д), учитывая равенства (5), (6).
Найдём разрешающую систему КБ. Выпишем проекцию системы (1) в пространство инвариантов с учётом равенств (5), (6). Для этого умножим первое уравнение в (1) на р2х и вычтем из него второе уравнение, умноженное на р 1 х, тогда вместо разности р 1 ф2х — рыр 1х в силу первого равенства в (6) запишем К, а вместо остальных слагаемых подставим их значения из равенств (5). Получим уравнение
К + ЩЩи + щку V — КУу ЩУ2 + ЩК\ Wl — К? К^2 + Щр 1 (их + V? + W1z) —
—кхр2 (их + V! + ^ух) — о,
поступая аналогично, выпишем первую часть разрешающей системы
и + тщи + иуи^ — иуиж + и?и^1 — и?и^ + —щ —2ип рх
т 1 ттх т2
—их--
р 2 р2
— — (щ — их) щ
V + VхVхи + КхКУV — КхцуV2 + VXV1 Wl — +(—Vх — — vA Ру —
\р l р2 )
— I 'рЛVх — Vх) ,
р1 ) т
W+WXWXu+WXW1УV — WXWУ+ WXWl1 Wl — W1XW21 №Г2 +1 —WX — — W1X) Р1
V рl р2 )
— — (* W2X — И?) Щ — Щ,
\Р1 ) Т
Р + VхКу — Кх^у) V + (VIК? — ВД1) Wl + Vхрl (их + Vу + Wl1) —
_—_Кхр — Р2 — ^2 р1 l у р l Т
¥2 + (WXКУ — ЩЩ) V + — КХГ ?) Wl + WXРl (их + Vу + WlZ) —
——КхР — р2 ^ — ^
р l р l Т
Рз + Vх ку — Кх^2у) V2 + (Vх к? — Кхух) + VI р2 (их + V! + W21) —
__У- У Р2 Д _ Т '
Л + ("Д - ) У2 + ("Д - ) " + ж2хр2 (и2Г + у2у + ") -
2 2 "> Р2 т
где Рт _ Рр1 Дт + РР2 Д™- Заметим, что полученные равенства являются алгебраическими уравнениями на функции из правых частей (5), (6). Вторая часть разрешающей системы состоит из условий совместности АС системы. Удобнее всего получать условия совместности на языке распределений Картана. В пространстве 1-струй 31 рассмотрим стандартные формы Картана в заданных координатах
Шг _ Фг - р^Ж - р^у¿у - р^^ - 2 _ 1, 2,
_ ¿м - м^ж - ¿у - ¿г - _ 3,4,
Шк _ ^ - ^¿ж - ¿у - ¿г - у^, к _ 5, 6,
шг _ ¿од - - ¿у - и^¿г - I _ 7, 8.
Для проверки совместности уравнений (5), (6) необходимо сузить формы Картана на поверхность, которую задают уравнения (5), (6) в пространстве 3и проверить является ли идеал, порождённый суженными формами Картана, дифференциальным. Условия, обеспечивающие последнее требование(дифференциальность идеала внешних форм), и будут условиями совместности АС системы.
Уравнения (5), (6) задают 8-мерную поверхность в 44-мерном пространстве 3 На этой поверхности, точнее в её касательном расслоении, после сужения форм Картана будет 8 уравнений, задающих распределение Картана. Выпишем их. Взяв за координаты на поверхности переменные (¿, ж, у, г, р1, р2, м1, м2) и подставив в формы Картана равенства (5), (6), получим
= - Д^ - (Я^Я? - Д^) ¿у - (Д^ - ЯД) ¿г - Д^ _ 0,
^2 = Ц^м - - (и2хиу - ЦЭД) ¿у - (ЦЩ - ) ¿г - _ 0, Пэ = УТ (У1р1 ¿р1 + У1р2¿р2 + У1«(^1 - ¿М2)) -
-Ух (У2Р1 ¿р1 + У2Р2¿р2 + - ¿М2)) -
- (УГУ* - УГУ2У) ¿у - (УГУ2 - УТУ") ¿г + (У2ГУи - УТУ,4 - У) ^ _ 0,
П4 = " (^1Р1 ¿р1 + ^1Р2¿р2 + ^«(¿М - ¿М2)) - "Г (^2Р1 ¿р1 + ^2Р1 ^1 +
(¿М1 - ¿М2)) - - ) ¿у - - "") ¿г+
+ ("а "и - "" - Ж) ^ _ 0,
П5 = УГ^р1 - ДГ (У1Р1 ¿р1 + У1Р2¿р2 + У„(^1 - ¿М2)) - (УГЯУ - ДГУУ) ¿у-
- (УГЯ1 - У") ¿г - (Д^ + л) ^ _ 0,
= ЖГ^р1 - ДГ (Ж1р1 ¿р1 + Ж1р2¿р2 + ^„(¿м - ¿М2)) - (ЖГДУ - ДТЖ) ¿у-- ("Д - Д?") ¿г - (Д^ + Л2) ^ _ 0, = УГ^р2 - ДГ (У2Р1 ¿р1 + У2Р2¿р2 + У2«(^М1 - ¿М2)) - (У'Д* - ДУ?)
- (У2гД2 - ДГУ2") ¿г - (Д^ + Лэ) ^ _ 0,
= - ДГ (Ж2Р1 ¿р1 + Ж2Р2¿р2 + - ¿М2)) - ("^ДУ - ДГЖ|) ¿у-
- (WxKz2 - R^WX) dz - RWzt + F4) dt = 0.
Интегральные многообразия распределения Картана на рассматриваемой поверхности должны иметь размерность 4, так как интегральное многообразие — это график произвольной функции переменных (t,x,y,z). Значит, среди уравнений Qi = 0,i = 1,... , 8, должны быть только 4 линейно независимых. Выписывая условия линейной зависимости, к примеру, последних четырёх форм от первых четырёх, получим первую часть условий совместности. Вторая часть уравнений совместности получится из требования замкнутости идеала, порождённого четырьмя независимыми формами Qi, относительно операции взятия внешнего дифференциала.
3. Заключение
Для системы уравнений двухфазной газодинамики найдено точное решение, являющееся дифференциально инвариантным относительно подалгебры, порождающей трансляции по координатам и галилеев сдвиг по первой координате. Сравнительно с поиском частично инвариантного решения на этой же подалгебре поиск дифференциально инвариантного решения оказался проще, что вызвано большим количеством накладываемых ограничений. Выведена часть уравнений группового расслоения системы (1) относительно этой же подалгебры. Оставшаяся часть уравнений может быть получена после выполнения стандартных выкладок, указанных в конце предыдущего пункта.
Автор выражает благодарность профессору Воронину Сергею Михайловичу, доктору физ.-мат. наук, профессору кафедры математического анализа ЧелГУ, за критику и обсуждения результатов работы.
Список литературы
1. Рахматулин, Х. А. Основы газодинамики взаимопроникающих движений сжимаемых сред / Х. А. Рахматулин // Прикл. математика и механика. — 1956. — Т. 20, № 2. — С. 184-195.
2. Сверхзвуковые двухфазные течения в условиях скоростной неравновесности частиц / Н.Н.Яненко, Р.И.Солоухин, А.Н.Папырин, В.М.Фомин. — Новосибирск : Наука, 1980. — 160 с.
3. Физико-математическое моделирование подавления детонации облаками мелких частиц / А.В.Федоров, П.А.Фомин, В.М.Фомин, Д. А.Тропин, Дж.-Р. Чен. — Новосибирск : НГАСУ, 2011. — 156 с.
4. Стояновская, О. П. Анализ численных алгоритмов расчёта быстрого обмена импульсом между газом и пылью при моделировании околозвездных дисков / О. П. Стояновская, Э.И.Воробьёв, В. Н. Снытников // Астроном. журн. — 2018. — Т. 95, № 7. — С. 487-500.
5. Овсянников, Л. В. Групповой анализ дифференциальных уравнений / Л. В. Овсянников. — М. : Наука, 1978. — 400 с.
6. Симметрии и законы сохранения уравнений математической физики / А. И. Бочаров, А. М. Вербовецкий, А. М. Виноградов и др. — М. : Факториал Пресс, 2005. — 380 с.
7. Олвер, П. Приложения групп Ли к дифференциальным уравнениям / П. Олвер. — М. : Мир, 1989. — 639 с.
8. Bluman, G. Symmetry and Intergration Methods for Differential Equations / G.Bluman, S.Anco. — New York : Springer-Verlag, 2002. — 419 p.
9. Федоров, В. Е. Инвариантные и частично инвариантные решения системы уравнений механики двухфазной среды / В.Е.Федоров, А.В.Панов // Вестн. Челяб. гос. ун-та. - 2011. - № 38 (253). - Физика. Вып. 11. - С. 65-69.
10. Panov, A. V. Invariant solutions and submodels in two-phase fluid mechanics generated by 3-dimensional subalgebras: barochronous flows / A. V. Panov // International Journal of Non-Linear Mechanics. - 2019. - Т. 116. - С. 140-146.
11. Хабиров, С. В. Классификация дифференциально инвариантных подмоделей / С. В. Хабиров // Сиб. мат. журн. - 2004. - Т. 45, № 3. - С. 682-701.
12. Панов, А. В. Точные решения уравнений динамики двухфазной среды. Коллапс газа и частиц в пространстве / А.В.Панов // Сиб. журн. индустр. математики. — 2017. - Т. 20, № 2. - С. 71-82.
13. Хабиров, С. В. Симметрийный анализ модели несжимаемой жидкости с вязкостью и теплопроводностью, зависящими от температуры / С. В. Хабиров. — Уфа : Гилем, 2004. - 37 c.
14. Камке, Э. Справочник по дифференциальным уравнениям в частных производных первого порядка / Э. Камке. - М. : Наука, 1966. - 260 c.
Поступила в 'редакцию 01.05.2020 После переработки 18.08.2020
Сведения об авторе
Панов Александр Васильевич, кандидат физико-математических наук, доцент кафедры математического анализа, Челябинский государственный университет, Челябинск, Россия; e-mail: gjd@bk.ru.
Chelyabinsk Physical and Mathematical Journal. 2020. Vol. 5, iss. 3. P. 316-326.
DOI: 10.47475/2500-0101-2020-15306
A DIFFERENTIAL INVARIANT SOLUTION OF TWO-PHASE FLUID DYNAMICS EQUATIONS
A.V. Panov
Chelyabinsk State University, Chelyabinsk, Russia
gjd@bk.ru
A system of differential equations, which describes the dynamics of a gas suspension in an isothermal case is considered. It is derived the submodel of homogeneous medium movements, which are invariant with respect to Galilean shifts at the first coordinate. This submodel is a differentially invariant submodel with respect to a four-dimensional subalgebra with a basis of shift operators over all spatial coordinates and a Galilean transformation operator over the first coordinate. This submodel has been led to an involution and integrated. An exact solution of the system is found. A comparison of the solution with a known partially invariant solution is done. A group bundle of the system is constructed with respect to the subalgebra from the optimal system of subalgebras of the Lie algebra of the symmetry group for the system of two-phase fluid dynamics. This group bundle divides the system into two parts: automorphic and resolution. All solutions of the system are contained in the set of solutions of the group bundle and vice versa.
Keywords: symmetry group, differential invariant solution, two-phase fluid.
References
1. Rakhmatulin Kh.A. Osnovy gazodinamiki vzaimopronikayushchikh dvizheniy szhimayemykh sred [Fundamentals of gasdynamics of interpenetrating motions of compressible media]. Prikladnaya matematika i mekhanika [Applied mathematics and mechanics], 1956, vol. 20, no. 2, pp. 184-195. (In Russ.).
2. YanenkoN.N., Soloukhin R.I., PapyrinA.N., FominV.M. Sverkhzvukovye dvukhfaznye techeniya v usloviyakh skorostnoy neravnovesnosti chastits [Supersonic two-phase flows under conditions of nonequilibrium of the velocities of the particles]. Novosibirsk, Nauka Publ., 1980. 160 p. (In Russ.).
3. FedorovA.V., FominP.A., FominV.M., TropinD.A., ChenJ.-R. Fiziko-matematicheskoe modelirovanie podavleniya detonatsii oblakami melkikh chastits [Physical and mathematical modeling of detonation quenching with clouds of fine particles]. Novosibirsk, Novosibirsk State University of Architecture and Civil Engineering, 2011. 156 p. (In Russ.).
4. Stoyanovskaya O.P., Snytnikov V.N., VorobyovE.I. Analysis of numerical algorithms for computing rapid momentum transfers between the gas and dust in simulations of circumstellar disks. Astronomy reports, 2018, vol. 62, no. 7, pp. 455-468.
5. Ovsyannikov L.V. Group Analysis of Differential Equations. New York, Academic Press, 1982. 432 p.
6. Bocharov A.I., Verbovetskii A.M., Vinogradov A.M. et al. Simmetrii i zakoni sokhraneniya uravneniy matematicheskoy fiziki [Symmetry and conservation laws of the equations of mathematical physics]. Moscow, Factorial Press Publ., 2005. 380 p. (In Russ.).
7. OlverP.J. Applications of Lie Groups to Differential Equations. New York, SpringerVerlag, 1986. 639 p.
The work is supported the Russian Science Foundation, project 18-71-00041.
8. BlumanG., Anco S. Symmetry and Intergration Methods for Differential Equations. New York, Springer-Verlag, 2002. 419 p.
9. FedorovV.E., Panov A.V. Invariantnye i chastichno invariantnye resheniya sistemy uravneniy mekhaniki dvukhphaznoy sredy [Invariant and partially invariant solutions of the system of equations of a two-phase medium mechanics]. Vestnik Chelyabinskogo gosudarstvennogo universiteta [Bulletin of Chelyabinsk State University], 2011, no. 38, pp. 65-69. (In Russ.).
10. Panov A.V. Invariant solutions and submodels in two-phase fluid mechanics generated by 3-dimensional subalgebras: barochronous flows. International Journal of Non-Linear Mechanics, 2019, vol. 116, pp. 140-146.
11. Khabirov S.V. Classification of differential invariant submodels. Siberian Mathematical Journal, 2004, vol. 45, pp. 562-579.
12. Panov A.V. Exact solutions of the equations of two-phase dynamics. Collapse of gas and particles in space. Journal of Applied and Industrial Mathematics, 2017, vol. 11, no. 2, pp. 263-273.
13. Khabirov S.V. Simmetriynyy analiz modeli neszhimayemoy zhidkosti s vyazkost'yu i teploprovodnost'yu, zavisyashchimi ot temperatury [Symmetry analysis of an incomressible fluid model with viscosity and heat conductivity that depends on temperature]. Ufa, Gilem Publ., 2004. (In Russ.).
14. Kamke E. Differentialgleichungen: Losungsmethoden und Losungen. I. Gewohnliche Differentialgleichungen. B. G. Teubner, Leipzig, 1977. xxvi+670 p.
Accepted article received 01.05.2020 Corrections received 18.08.2020