CC BY
19
Вестник Челябинского государственного университета. 2012. № 30 (284). Физика. Вып. 14. С. 23-25.
МЕХАНИКА СПЛОШНЫХ СРЕД
Ю. М. Ковалев, В. Ф. Куропатенко
АНАЛИЗ ИНВАРИАНТНОСТИ ОТНОСИТЕЛЬНО ПРЕОБРАЗОВАНИЯ ГАЛИЛЕЯ МАТЕМАТИЧЕСКОЙ МОДЕЛИ «ЗАМОРОЖЕННОЙ» ГАЗОВЗВЕСИ
Проведен анализ инвариантности относительно преобразования Галилея математической модели «замороженной» газовзвеси. Показано, что математическая модель «замороженной» газовзвеси не является инвариантной относительно преобразования Галилея. Это приводит к появлению фиктивного источникового члена в уравнении энергии.
Ключевые слова: математическая модель, инвариантность, многокомпонентная смесь.
В связи с развитием современной вычислительной техники резко возросла роль математического моделирования физических процессов, используемого в науке и технике. Более того, есть такие проблемы, когда математическое моделирование является единственным средством предварительного изучения явлений. Поэтому с особой остротой встает вопрос адекватности математических моделей тем физическим процессам, которые они пытаются описывать. В природе практически нет чистых веществ, поэтому активно развиваются математические модели многокомпонентных сред [1-2]. Для верификации расчетов используют известные экспериментальные данные. Очень важно, чтобы условия проведения расчетов и экспериментов совпадали. В настоящей статье на примере анализа математической модели замороженной газовзвеси [3-4] мы покажем, к чему может привести ситуация, когда расчеты и эксперимент проведены в разных системах координат.
При решении поставленной задачи предполагалось, что частицы твердой фазы неподвижны и несжимаемы. Это означает, что вместо газовзвеси фактически рассматривается заполненная газом недеформируемая решетка. Твердые частицы имитируют ее узлы, а связи между узлами решетки не оказывают влияния на газодинамическое течение, т. е. используется модель «замороженной» газовзвеси, представленная в работах [3-4] при изучении ослабления ударных волн. Поскольку частицы неподвижны и несжимаемы, их объемная концентрация и, следовательно, объемная концентрация газа постоянны.
С учетом сказанного выше система уравнений из [3-4], описывающая в одномерном случае течение газа через решетку, имеет вид
Эр Э(ри)
дґ дх
д(ри) д(рии)
= 0,
дґ дх
д(рЕ) д(риЕ) д( Ри)
дґ
дх
дх
= ~в-
(1)
(2)
(3)
Здесь Р — давление, р — плотность, и — скорость, ґ — время, ^ — силы межфазного взаимодействия, Е—удельная полная энергия газа; Q — интенсивность теплообмена между газом и частицами. Функция ^ зависит от разности скоростей газа и частиц, функция Q — от разности температур газа и частиц. Функции ^ и Q не изменяются при переходе в новую систему координат.
Проведем анализ инвариантности системы уравнений (1)—(3) относительно преобразования Галилея. С этой целью перейдем в новую систему координат, которая движется с постоянной скоростью В относительно старой системы координат. Скорость в новой системе координат будет равна
и и = и + В,
координата определяется из уравнения хн = х + Вґ.
(4)
(5)
Производные по координате и времени определяются следующим образом:
В. (6)
дх дхг
Гд1 Гд1 Г д Ї
= +
дґ дґ дх„
V У V У \ н у
1 Работа выполнена при поддержке РФФИ, проект После перех°да в движущуюся систему к°°р-
№ 10-01-00032. динат значок Н будем опускать. Следовательно,
уравнение неразрывности газовой фазы (1) с учетом (4)-(6) принимает вид
dp + dp D + d(p(u - D) = 0
dt dx
dx
(7)
который после сокращения членов с противоположными знаками совпадает с (1).
Запишем теперь уравнение сохранения импульса газовой фазы (2) в новой системе координат:
др(и - Р) + др(и - Р) в + др(и - Р)2 + др + р = 0 дґ дх дх дх
После несложных преобразований оно принимает вид
дри дри2 дР „
+ —+ р = Ш1(Р), (8)
дґ дх дх
где
K)1(D) = D
ґ dp + dpu 2 dpu
dt dx
dx
dp dp
dx dx
(9)
Подставив (1) в (9) и сократив подобные члены, получим, что
ц(Р) = 0 (10)
и, следовательно, уравнение (8) совпадает с уравнением (2).
И, наконец, перейдем в новую систему координат в уравнении для удельной энергии газовой фазы (3). Учитывая, что
ІГ +и2
Е = £+----,
2
где є — удельная внутренняя энергия, запишем уравнение (3) в новой системе координат:
др^е + |(м -Р)2^ др^є + ^(и -Р)2
dt
dx
1
dp(u -D)l,£ + 2(u -D)2) dP(u -D) Q 0
+---------^------------ +—^------ + Q = 0.
dx
dx
Раскрыв скобки и сгруппировав члены, получим уравнение в новой системе координат для удельной полной энергии газовой фазы, распространяющейся с постоянной скоростью Р:
dp
.2 Л
є+
dpu
dt
где ю2 = — DF.
+
2
є+
dx
dPu
+ —— + Q = ®2,(11) dx
Как следует из уравнений (8), (10) и (11), для модели замороженной газовзвеси из [34] уравнение неразрывности газовой фазы и уравнение сохранения импульса газовой фазы являются инвариантными относительно преобразования Галилея, а уравнение энергии (3) не является инвариантным.
Оценим последствия неинвариантности уравнения энергии. В уравнении (11) исключим кинетическую энергию с помощью уравнения (2). Для этого умножим (2) на и и вычтем из (11), затем умножим (1) на £ и вычтем из полученного уравнения. Следствием этих преобразований является уравнение для внутренней энергии
дє дє P
-------+ u-------------—
dt dx p
dp dp
— + u — dt dx
u - D
F. (12)
Перейдем к субстанциональным производным, заменим плотность удельным объемом
V=
т
vpy
и сравним полученное уравнение
с уравнением для удельной внутренней энергии как функции энтропии и удельного объема
d є dV dS
— + P = T—.
dt dt dt
(13)
В результате из (12) и (13) получим уравнение производства энтропии газа:
Т^ = Р-(^(и - Р) - б). ш р
Если разделить энтропию на две части
^ = «рн + ^,
где SPH определяется «физикой» модели, а — Галилеевой неинвариантностью, то мы
получим уравнение производства энтропии
= F (u - D),
dt p
(14)
возникшее исключительно вследствие того, что авторы модели [3-4] пренебрегли фундаментальным принципом механики.
Рассмотрим возможность получения уравнения полной удельной энергии газовой фазы, инвариантного относительно преобразования Галилея. С этой целью из правой части уравнения полной удельной энергии газовой фазы (3)
вычитаем работу сил межфазного взаимодействия Ри, получаем
д(рЕ) д(риЕ) д( Ри)
дґ
дх
дх
= -Q - Ри. (15)
Переходя к новой системе координат, которая движется с постоянной скоростью Р, по формулам (4)-(6) получим уравнение полной удельной энергии газовой фазы в новой системе координат
д(р(є + и2 /2) + д(ри(є + и2/2)) +
дґ
дх + дРи
+-------= ^ - ри.
дх *
(16)
Из сравнения (15) и (16) видно, что уравнение полной удельной энергии газовой фазы (15) является инвариантным относительно преобразования Галилея.
Получим выражение для удельной внутренней энергии газовой фазы. С этой целью, как и ранее, исключим кинетическую энергию с помощью уравнения (2). Для этого умножим (2) на и и вычтем из уравнения (16). Далее, умножим (1) на е и вычтем из (16). В результате получаем уравнение для удельной внутренней энергии газовой фазы в виде
де де р .Эр Эр. _ .
+ и -^(^ + ид!) = -б / р. (17)
at дх р at дх
Как и ранее, переходя к субстанциональным производным и заменяя плотность удельным объемом V = (1/р), получим уравнение для удель-
ной внутренней энергии газовой фазы как функции энтропии и удельного объема (13). Вычитая из равенства (17) равенство (13), получим уравнение производства энтропии газовой фазы
Г- = 43 / р.
(18)
Уравнение (18) согласуется со вторым законом термодинамики. Видим, что производство энтропии газовой фазы определяется только меж-фазным теплообменом.
К сожалению, принцип инвариантности к преобразованию Галилея не выполняется в ряде моделей многокомпонентных сред, публикуемых в журналах. Такие модели не способны прогнозировать результаты тех физических процессов, для моделирования которых предназначены.
Список литературы
1. Нигматулин Р. И. Основы механики гетерогенных сред. М. : Наука, 1978. 336 с.
2. Куропатенко В. Ф. Новые модели механики сплошных сред // Инж.-физ. журн. 2011. Т. 84, № 1. С. 74-92.
3. Кругликов Б. С., Кутушев А. Г. Ослабление воздушных ударных волн экранирующими решетками // Физика горения и взрыва. 1988. № 1. С. 115-117.
4. Кругликов Б. С., Кутушев А. Г. Ослабление воздушных ударных волн слоями запыленного газа и решетками // Приклад. механика и техн. физика. 1988. № 1. С. 51-57.
