Журнал вычислительной математики и математической физики, 2020, T. 60, № 4, стр. 621-625
C.K. Годунов и кинетическая теория в ИПМ им. М.В. Келдыша РАН
С. З. Аджиев 2, Я. Г. Батищева 1, В. В. Веденяпин 1, 4, *, Ю. А. Волков 1, В. В. Казанцева 1, И. В. Мелихов 2, М. А. Негматов 3, Ю. Н. Орлов 1, Н. Н. Фимин 1, В. М. Чечеткин 1
1 Федеральный Исследовательский Центр Институт прикладной математики им. М.В. Келдыша РАН
125047 Москва, Миусская пл., 4, Россия
2 ГСП–1, МГУ им. М.В. Ломоносова, химический факультет
119991 Москва, Ленинские горы, 1/3, Россия
3 Акционерное общество
“Центральный научно-исследовательский институт машиностроения”
141070 Московская область, Королeв, ул. Пионерская, 4, Россия
4 Российский университет дружбы народов
117198 Москва, ул. Миклухо-Маклая, 6, Россия
* E-mail: vicveden@yahoo.com
Поступила в редакцию 14.11.2019
После доработки 14.11.2019
Принята к публикации 16.12.2019
Аннотация
Описана история развития сотрудничества коллектива ИПМ им. М.В. Келдыша РАН с С.К. Годуновым. В процессе этого сотрудничества получено много интересных результатов в теории кинетических уравнений и вычислительной математике. Библ. 20.
В 1971 г. вышел обзор (см. [1]). Обзор был первой статьeй в СССР по дискретным моделям уравнения Больцмана и содержал исследование модели с тремя скоростями. Как выяснилось позже, похожую модель выписал Больцман в 1872 г., но в настоящее время эта модель известна как модель Годунова–Султангазина [2], хотя за рубежом ее часто называют одномерной моделью Бродуэлла.
В обзоре было сформулировано много вопросов, задач и предположений. Одно из таких предположений было сформулировано следующим образом: “Выбранная нами модельная система на первый взгляд не намного сложнее, чем модель Карлемана; oднако детальный анализ этой системы показывает, что она не имеет квадратичных диссипирующихся интегралов, в связи с чем получение для нее теоремы существования, охватывающей все положительные времена $0 < t < \infty $, затруднительно”. Сотрудник ИПМ им. М.В. Келдыша B.B. Веденяпин на основе этого поставил ряд задач: 1) выявить разницу между двумя моделями, 2) понять существо диссипирующихся интегралов, 3) установить, как разница между моделями влияет на свойства уравнения Больцмана. В результате работы в этих направлениях была получена теорема о единственности $H$-функции Больцмана: для газа, описываемого уравнением Больцмана, энтропия является единственным экстенсивным (аддитивным) функционалом, возрастающим со временем. Также была выяснена причина единственности Н-функции как диссипирующего интеграла для моделей Бродуэлла и Годунова–Султангазина и неединственности для модели Карлемана. Это стало основным результатом кандидатской и докторской диссертаций В.В. Веденяпина, и опубликовано в УМН [3].
Сергей Константинович Годунов тесно сотрудничал с ИПМ им. М.В. Келдыша РАН. В одно из своих посещений ИПМ он обратил внимание отдела кинетических уравнений на дважды дивергентную форму гидродинамических следствий дискретных моделей (1) уравнения Больцмана:
В это же время вышла книга по кинетическим уравнениям [4]. Вышеприведенные формы Годунова были туда включены. Но С.К. Годунов не остановился на полученном, и предложил сделать то же самое для уравнения Власова. В.В. Веденяпин и М.А. Негматов эту проблему решили достаточно быстро, и решили опубликовать соответствующие результаты. Рецензент в отзыве на их работу написал, что основным результатом по новизне является как раз форма Годунова, в соответствии с чем предложил изменить название статьи. Авторы учли это пожелание, и даже написали вторую работу; так что были опубликованы практически одновременно две статьи с именем Годунова [5], [6]. Далее мы приведем полученную в этих работах форму Годунова для уравнений магнитной гидродинамики.
Двухжидкостная и одножидкостная МГД (или ЭМГД) с ненулевой температурой имеют много модификаций; для примера мы рассмотрим систему:
Обычно получают уравнения гидродинамического типа из системы кинетических уравнений, последовательно интегрируя и вводя моменты:
(1)
${{\sigma }_{{\alpha k\ell }}} = \int {({{p}_{k}} - {{P}_{{\alpha k}}})} ({{p}_{\ell }} - {{P}_{{\alpha \ell }}}){{f}_{\alpha }}({\mathbf{x}},{\mathbf{p}},t){{d}^{3}}{\mathbf{p}}\quad ({\text{тензор}}\;{\text{напряжений}}),$В этой системе первое уравнение представляет собой уравнение неразрывности, второе есть уравнение движения, a третье уравнение описывает изменение энергии хаотического движения. Эта точная система уравнений, но она не замкнута. Чтобы ее замкнуть, предполагают максвеллизацию либо из парных столкновений (добавляя интеграл столкновений), либо из взаимодействия со средой (добавляя линейный интеграл столкновений). Это означает, что моменты высшего порядка определяются через низшие с помощью максвелловского распределения
(2)
$f_{\alpha }^{0}({\mathbf{x}},{\mathbf{p}},t) = \frac{{{{n}_{\alpha }}({\mathbf{x}},t)}}{{\mathop {\left( {2{{k}_{{\text{B}}}}\pi {{T}_{\alpha }}{{m}_{\alpha }}} \right)}\nolimits^{3/2} }}exp\left( { - \frac{{\mathop {\left( {{\mathbf{p}} - {{{\mathbf{P}}}_{\alpha }}} \right)}\nolimits^2 }}{{2{{k}_{{\text{B}}}}{{T}_{\alpha }}{{m}_{\alpha }}}}} \right).$При этом оказывается, что:
(3)
$\begin{gathered} {{G}^{\alpha }}(\beta _{\mu }^{\alpha }) = \int {f_{\alpha }^{0}} (\beta _{\mu }^{\alpha }){{d}^{3}}{\mathbf{p}},\quad \mu = \overline {0,4} , \\ f_{\alpha }^{0}(\beta _{\mu }^{\alpha }) = \exp [\beta _{0}^{\alpha } + \beta _{1}^{\alpha }{{p}_{1}} + \beta _{2}^{\alpha }{{p}_{2}} + \beta _{3}^{\alpha }{{p}_{3}} + \beta _{4}^{\alpha }{{{\mathbf{p}}}^{2}}]. \\ \end{gathered} $(4)
$\frac{{\partial G_{{{{\beta }_{\mu }}}}^{\alpha }}}{{\partial t}} + \frac{{\partial G_{{{{\beta }_{\mu }}{{\beta }_{i}}}}^{\alpha }}}{{\partial {{x}_{i}}}} + K_{\mu }^{\alpha } = 0,\quad {\text{где}}\quad G_{{{{\beta }_{\mu }}}}^{\alpha } \equiv \frac{{\partial {{G}^{\alpha }}}}{{\partial {{\beta }_{\mu }}}}.$Такие уравнения являются источниками экономичных разностных схем “солверов” ($solvers$): система уравнений (4) coдержит только одну неизвестную функцию ${{G}^{\alpha }}$ (для каждого сорта частиц $\alpha $).
Вот еще одно замечание из вышеупомянутого обзора Годунова и Султангазина 1971 г., которое привело к новым результатам, полученным позднее в ИПМ им. М.В. Келдыша: “Следует заметить, что способ выбора дискретных скоростей, обеспечивающих описание течения газа с законами сохранения трех компонент импульса и с сохранением энергии, в общем случае не разработан. Эта разработка связана с трудностями комбинаторно-геометрического характера” [1]. Теперь можно уверенно сказать, что такой способ выбора дискретных скоростей разработан. В этой разработке помогла аналогия между гамильтоновой механикой и дискретными уравнениями Больцмана [7]–[11].
Задачи и перспективные направления исследований, поставленные в обзоре Годунова и Султангазина [1], стали развиваться и дальше, в том числе получены законы сохранения для уравнений химической кинетики и уравнения Лиувилля [12]–[18]. А результаты для уравнения Власова сейчас переносятся на уравнения Власова–Максвелла–Эйнштейна [19]–[21].
Список литературы
Годунов С.K., Султангазин У.M. О дискретных моделях кинетического уравнения Больцмана // Успехи матем. наук. 1971. Т. 26. № 3. С. 1–51.
Vasileva O.A., Dukhnovsky S.A., Radkevich E.V. On the nature of the local equilibrium of the Carleman and Godunov–Sultangazin equations // J. Math. Sci. 2018. V. 235. № 4. P. 392–454.
Веденяпин В.В. Дифференциальные формы в пространстве без нормы // Успехи матем. наук. 1988. Т. 43. № 1. С. 159–179.
Веденяпин В.В. Кинетические уравнения Больцмана и Власова. M.: Физматлит, 2001. Vedenyapin V.V. Boltzmann and Vlasov Kinetic Equations. Moscow: Fizmatlit, 2001.
Веденяпин В.В., Негматов M.A. О выводе и классификации уравнений типа уравнения Власова и магнитной гидродинамики. Тождество Лагранжа и форма Годунова // Теор. матем. физ. 2013. Т. 170. № 3. С. 468–480. Vedenyapin V.V., Negmatov M.A. Derivation and Magneto Hydrodynamics equations: Lagrange identity and Godunov’s form // Theoret. Math. Phys. 2013. V. 170. № 3. P. 468–480.
Vedenyapin V.V., Negmatov M.A. On the derivation and classification of equations of Vlasov-type and magnetic hydrodynamics. Lagrange’s identity, Godunov form and critical mass // J. Math Sci. (N.Y.) V. 202. № 5. P. 769–782.
Брюно А.Д. Ограниченная задача трех тел: плоские периодические движения. M.: Наука, 1990. Bruno A.D. The restricted 3-body problem: plane periodic orbits. Berlin: Walter de Gruyter & Co., 1994.
Арансон А.Б. Вычисление многогранника Ньютона, в сб. Материалы международн. конф. и Чебышевских чтений, посв. 175-летию со дня рождения П.Л. Чебышева, Т. 1. С. 32. Ред. Н.С. Бахвалов, С.Н. Кружков, К.Ю. Богачев и др. M.: Изд-во МГУ, 1996. Aranson A. Calculating the Newton polyhedron, in: Materials of the International Conference and the Chebyshev Readings Devoted to the 175th Birthday of P.L. Chebyshev. V. 1. P. 32. Ed. N.S. Bakhvalov et al. Moscow: Moscow Univ. Press, Moscow, 1996.
Веденяпин В.В., Орлов Ю.Н. О законах сохранения для полиномиальных гамильтонианов и для дискретных моделей уравнения Больцмана // Теор. матем. физ. 1999. Т. 121. № 2. С. 307–315. Vedenyapin V.V., Orlov Yu.N. On conservation laws for polynomial Hamiltonians and for discrete models of the Boltzmann equation // Theoret. Math. Phys. 1999. V. 121. № 2. P. 307–315.
Bobylev A.V., Vinerean M.C. Construction of discrete models with given invariants // J. Stat. Phys. 2008. V. 132. № 1. P. 153–170.
Bernhoff N., Vinerean M. Discrete velocity models for mixtures without nonphysical collision invariants // J. Stat. Phys. 2016. V. 165. № 2. P. 434–453.
Батищева Я.Г., Веденяпин В.В. Второй закон термодинамики для химической кинетики // Матем. моделирование. 2005. Т. 17. № 8. С. 106–110. Batishcheva Ya.G., Vedenyapin V.V. The second law of thermodynamics for chemical kinetics // Mat. Model. 2005. V. 17. № 8. P. 106–110.
Веденяпин В.В., Аджиев С.З. Энтропия по Больцману и Пуанкаре // Успехи матем. наук. 2014. Т. 69. № 6. С. 45–80. Vedenyapin V.V., Adzhiev S.Z. Entropy in the sense of Boltzmann and Poincare // Russian Math. Surv. 2014. V. 69. № 6. P. 45–80.
Аджиев С.З., Амосов С.А., Веденяпин В.В. Одномерные дискретные модели кинетических уравнений для смесей // Ж. вычисл. матем. и матем. физ. 2004. Т. 44. № 3. С. 553–558. Adzhiev S.Z., Amosov S.A., Vedenyapin V.V. One dimensional discrete models of kinetic equations for mixtures // Comput. Math. Math. Phys. 2004. V. 44. № 3. P. 523–528.
Аджиев С.З., Веденяпин В.В. Временные средние и экстремали Больцмана для марковских цепей, дискретного уравнения Лиувилля и круговой модели Каца // Ж. вычисл. матем. и матем. физ. 2011. Т. 51. № 11. С. 2063–2074. Adzhiev S.Z., Vedenyapin V.V. Time averages and Boltzmann extremals for Markov chains, discrete Liouville equations, and the Kac circular model // Comput. Math. Math. Phys. 2011. V. 51. № 11. P. 1942–1952.
Аджиев С.З., Веденяпин В.В., Волков В.В., Мелихов И.В. Обобщенные уравнения типа Больцмана для агрегации в газе // Ж. вычисл. матем. и матем. физ. 2017. Т. 57. № 12. С. 2065–2078. Adzhiev S.Z., Vedenyapin V.V., Volkov Yu.A., Melikhov I.V. Generalized Boltzmann-type equations for aggregation in gases // Comput. Math. Math. Phys. 2017. V. 57. № 12. P. 2065–2078.
Веденяпин B.B., Аджиев С.З., Казанцева B.B. Энтропия по Больцману и Пуанкаре, экстремали Больцмана и метод Гамильтона–Якоби в негамильтоновой ситуации // СМФН. 2018. Т. 64. № 1. С. 37–59. Vedenyapin V.V., Adzhiev S.Z., Kazantseva V.V. The entropy of the Boltzmann and the Poincare, the Boltzmann extremals, and the Hamilton–Jacobi method in the non-Hamiltonian context // Contemporary Mathematics. Fundamental Direct. 2018. V. 64. № 1. P. 37–59.
Веденяпин B.B., Негматов M.A., Фимин Н.Н. Уравнения типа Власова и Лиувилля, их микроскопические, энергетические и гидродинамические следствия // Изв. РАН. 2017. Т. 81. № 3. С. 45–82. Vedenyapin V.V., Negmatov M.A., Fimin N.N. Vlasov- and Liouville-type equations for their macroscopic, energetic and hydrodynamic consequences // Izv. Math. 2017. V. 81. № 3. P. 505–541.
Веденяпин B.B., Фимин Н.Н., Чечеткин В.М. Уравнение типа Власова–Максвелла–Эйнштейна и переход к слаборелятивистскому приближению // Ж. вычисл. матем. и матем. физ. 2019. Т. 59. № 11. С. 81–97. Vedenyapin V.V., Fimin N.N., Chechetkin V.M. Vlasov–Maxwell–Einstein-type equation and transition to a weakly relativistic approximation // Comp. Math. Math. Phys. 2019. V. 59. № 11. P. 81–97.
Веденяпин B.B., Фимин Н.Н., Чечеткин В.М. К вопросу о выводе уравнения Власова–Максвелла–Эйнштейна и его связь с космологическим лямбда-членом // Вестник МГОУ. Сер. Физика–Математика. 2019. № 2. С. 24–48.
Vedenyapin V.V., Fimin N.N., Chechetkin V.M. The system of Vlasov–Maxwell–Einstein-type equations and its nonrelativistic and weak relativistic limits // Internat. Journal of Modern Phys. D. 2020. № 1. P. 23. 2050006. Doi:10.1142/S0218271820500066
Дополнительные материалы отсутствуют.
Инструменты
Журнал вычислительной математики и математической физики