CC BY
29
Математика
УДК 532.5
АНАЛИЗ ИНВАРИАНТНОСТИ НЕКОТОРЫХ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ МНОГОКОМПОНЕНТНЫХ СРЕД1
Ю.М. Ковалев2, В.Ф. Куропатенко3
Проведен анализ инвариантности относительно преобразования Галилея математической модели «замороженной» газовзвеси. Показано, что математическая модель «замороженной» газовзвеси не является инвариантной относительно преобразования Галилея. Это приводит к появлению фиктивного источникового члена в уравнении энергии.
Ключевые слова: математическая модель, инвариантность, многокомпонентная смесь.
В связи с развитием современной вычислительной техники резко возросла роль математического моделирования физических процессов, используемых в науке и технике. Более того, есть такие проблемы, когда математическое моделирование является единственным средством предварительного изучения явлений. Поэтому с особой остротой встает проблема адекватности математических моделей тем физическим процессам, которые они пытаются описывать. В природе практически нет чистых веществ, поэтому активно развиваются математические модели многокомпонентных сред [1, 2]. Для верификации расчетов используют известные экспериментальные данные. Очень важно, чтобы условия проведения расчетов и экспериментов совпадали. В настоящей статье на примере анализа математической модели замороженной газовзвеси [3, 4] покажем, к чему может привести ситуация, когда расчеты и эксперимент проведены в разных системах координат.
При решении поставленной задачи предполагалось, что частицы твердой фазы неподвижны и несжимаемы. Это означает, что вместо газовзвеси фактически рассматривается заполненная газом недеформируемая решетка. Твердые частицы имитируют ее узлы, а связи между узлами решетки не оказывают влияния на газодинамическое течение, т.е. используется модель «замороженной» газовзвеси, представленная в работах [3, 4] при изучении ослабления ударных волн. Поскольку частицы неподвижны и несжимаемы, то их объёмная концентрация и, следовательно, объёмная концентрация газа постоянны.
С учётом сказанного выше система уравнений из [3, 4], описывающая в одномерном случае течение газа через решётку, имеет вид
Эр+^ри) = 0, (1)
д( дх
д(Ри) + д(рии) = _дР _ р (2)
дt дх дх
д(РЕ) + д(риЕ) + д(Ри) = (3)
дt дх дх Q. ( )
Здесь Р - давление, р - плотность, и - скорость, t - время, Е - силы межфазного взаимодействия, Е и £ - удельная полная и удельная внутренняя энергии газа; Q - интенсивность теплообмена между газом и частицами. Функция Е зависит от разности скоростей газа и частиц, функция Q - от разности температур газа и частиц. Функции Е и Q не изменяются при переходе в новую систему координат.
1 Работа выполнена при поддержке РФФИ, проект № 10-01-00032.
2 Ковалёв Юрий Михайлович - доктор физико-математических наук, профессор, заведующий кафедрой вычислительной механики сплошных сред, Южно-Уральский государственный университет.
E-mail: [email protected]
3 Куропатенко Валентин Фёдорович - доктор физико-математических наук, профессор, кафедра вычислительной механики сплошных сред, Южно-Уральский государственный университет.
4 Вестник ЮУрГУ, № 11, 2012
Ковалев Ю.М., Анализ инвариантности некоторых
Куропатенко В.Ф. математических моделей многокомпонентных сред
Проведем анализ инвариантности системы уравнений (1)-(3) относительно преобразования Галилея. С этой целью перейдем в новую систему координат, которая движется с постоянной скоростью О относительно старой системы координат. Скорость в новой системе координат будет равна
ин = и + О , (4)
координата определяется из уравнения
хн = х + Оґ. (5)
Производные по координате и времени определяются следующим образом:
—, (1)=Гі, +
Эх Эхн І Эґ у І Эґ
Г Э 1
Эхн у
О (6)
После перехода в движущуюся систему координат значок Н будем опускать. Следовательно, уравнение неразрывности газовой фазы (1) с учетом (4)-(6) принимает следующий вид:
Ъ+р + Э(р(и -О) = 0, (7)
ґ х х
который после сокращения членов с противоположными знаками совпадает с (1).
Запишем теперь уравнение сохранения импульса газовой фазы (2) в новой системе координат:
Эр(и - О) + Эр(и - О) о + др(и°)2 + др + р = о
дґ Эх Эх Эх
После несложных преобразований оно принимает вид
Эри Эри2 ЭР „ ...
-^- + ^— + ^- + Р = й!(О), (8)
Эґ Эх Эх
где
ЩО) = О(М+р -2р + О2Г^-р (9)
І Э ґ Эх Эх у І Эх Эх у
Подставив (1) в (9) и сократив подобные члены, получим
о\(О) = 0 (10)
и таким образом уравнение (8) совпадает с уравнением (2).
И, наконец, перейдём в новую систему координат в уравнении для удельной энергии газовой
и2
фазы (3). Учитывая, что Е = е + ~, запишем уравнение (3) в новой системе координат:
^ ~\2 | 1і п\2 1 г>\1~ , 1п\2
ре+-(и - о) ) + -(и- о) у о |Эр(и - О) [е+2(и - О) у + №(и - о)+в=о
Эґ Эх ёх Эх
Раскрыв скобки и сгруппировав члены, получим уравнение для удельной полной энергии газовой фазы в новой системе координат, распространяющейся с постоянной скоростью О,
Г и 2 1
Эр е + — Эри
V________У
Г и21
е + — 2
V у
ЭРи
, +-----^--------+1Г+е = (11)
at ах ах
где а>2 = - РЕ.
Как следует из уравнений (8), (10) и (11), для модели «замороженной» газовзвеси из [3, 4] уравнение неразрывности газовой фазы и уравнение сохранения импульса газовой фазы являются инвариантными относительно преобразования Галилея, а уравнение энергии (3) не является инвариантным.
Оценим последствия неинвариантности уравнения энергии. В уравнении (11) исключим кинетическую энергию с помощью уравнения (2). Для этого умножим (2) на и и вычтем из (11). Затем умножим (1) на £ и вычтем из (11). В результате получим уравнение для внутренней энергии:
Математика
д£ д£ Р (др др\ О (и -Р
— + и--------1 — + и — I + — =
дt дх р2 ^ дt дх) р
Е . (12)
р
Перейдём к субстанциональным производным, заменим плотность удельным объёмом V = V р и сравним полученное уравнение с уравнением для удельной внутренней энергии, как функции энтропии и удельного объёма:
£ = Т(13)
Л Л Л
В результате из (12) и (13) получим уравнение производства энтропии газа
Т ^ = -(Е(и - Р) - О .
Л р
Если разделить энтропию на две части
£ = $рН + ,
где «РН определяется «физикой» модели, а - Галилеевой неинвариантностью, то мы полу-
чим уравнение производства энтропии
Т^ = Е(и - Р), (14)
Л р
возникшее исключительно из-за того, что авторы модели [3, 4] пренебрегли фундаментальным принципом механики.
К сожалению, принцип инвариантности к преобразованию Галилея не выполняется в ряде моделей многокомпонентных сред, публикуемых в журналах. Такие модели не способны прогнозировать результаты тех физических процессов, для моделирования которых они предназначены.
Литература
1. Нигматулин, Р.И. Основы механики гетерогенных сред / Р.И. Нигматулин. - М.: Наука, 1978. - 336 с.
2. Куропатенко, В.Ф. Новые модели механики сплошных сред / В.Ф. Куропатенко // ИФЖ. -2011. - Т. 84, № 1. - С. 74-92.
3. Кругликов, Б.С. Ослабление воздушных ударных волн экранирующими решётками / Б.С. Кругликов, А.Г. Кутушев // ФГВ. - 1988. - № 1. - С. 115-117.
4. Кругликов, Б.С. Ослабление воздушных ударных волн слоями запыленного газа и решетками/ Б.С. Кругликов, А.Г. Кутушев // ПМТФ. - 1988. - № 1. - С. 51-57.
Поступила в редакцию 13 марта 2012 г.
6
Вестник ЮУрГУ, № 11, 2012
Ковалев Ю.М., Куропатенко В.Ф.
Анализ инвариантности некоторых математических моделей многокомпонентных сред
ANALYSIS OF THE INVARIANCE SOME MATHEMATICAL MODELS OF MULTICOMPONENT MEDIA
Yu.M. Kovalev1, V.F. Kuropatenko2
The analysis of the invariance under the Galilean transformation of the mathematical model of "frozen" gas suspension is done. It is shown that the mathematical model of the "frozen" gas suspension is not invariant under the Galilean transformations. This leads to appearance of a fictitious source term in the energy equation.
Keywords: mathematical model, invariance, multi-component mixture.
References
1. Nigmatulin R.I. Osnovy mekhaniki geterogennykh sred (Fundamentals of mechanics of heterogeneous media). Moscow, Nauka, 1978. 336 p. (in Russ.).
2. Kuropatenko V.F. Novye modeli mekhaniki sploshnykh sred (New models of continuum mechanics). Inzhenerno-Fizicheskii Zhurnal. 2011. Vol. 84, no 1. pp. 74-92. (in Russ.). [Kuropatenko V.F. New models of continuum mechanics. Journal of Engineering Physics and Thermophysics. 2011. Vol. 84, no. 1. pp. 77-99.]
3. Kruglikov B.S., Kutushev A.G. Fizika gorenija i vzryva. 1988. no. 1. pp. 115-117. (in Russ.).
4. Kruglikov B.S., Kutushev A.G. Oslablenie vozdushnykh udarnykh voln slojami zapylennogo gaza i reshetkami (Attenuation of air shock layers of dust and gas grills). Prikladnaja mekhanika i tekhnicheskaja fizika. 1988. no. 1. pp. 51-57. (in Russ.).
1 Kovalev Yury Mikhailovich is Dr. Sc. (Physics and Mathematics), Professor, Head of the Computational Continuum Mechanics Department, South Ural State University.
E-mail: [email protected]
2 Kuropatenko Valentin Fedorovich is Dr. Sc. (Physics and Mathematics), Professor, Computational Continuum Mechanics Department, South Ural State University.
