Метод Чепмена — Энскога (Bymk; Cyhbyug — |uvtkig)

Перейти к навигации Перейти к поиску

Метод Чепмена — Энскога (теория Чепмена — Энскога) — метод решения кинетического уравнения Больцмана[1]. На его основе могут быть получены уравнения газовой гидродинамики из уравнения Больцмана. Этот метод оправдывает феноменологические определяющие соотношения, возникающие в гидродинамических описаниях, таких как уравнения Навье — Стокса. При этом получаются выражения для различных коэффициентов переноса, таких как теплопроводность и вязкость, через молекулярные параметры. Таким образом, теория Чепмена — Энскога представляет собой важный шаг в переходе от микроскопического описания, основанного на частицах, к континуальному гидродинамическому описанию.

Теория названа в честь Сидни Чепмена и Дэвида Энскога, которые представили её независимо в 1916 и 1917 годах[2].

Отправной точкой метода Чепмена — Энскога является уравнение Больцмана для одночастичной функции распределения. :

где  — нелинейный интегральный оператор, моделирующий эволюцию при межчастичных столкновениях. Эта нелинейность затрудняет решение полного уравнения Больцмана и мотивирует разработку приближённых методов, таких как метод, предоставляемый теорией Чепмена — Энскога.

Учитывая эту отправную точку, различные предположения, лежащие в основе кинетического уравнения Больцмана, переносятся и на теорию Чепмена — Энскога. Самый простой из них требует разделения масштаба между продолжительностью столкновения и среднее свободное время между столкновениями  : . Это условие гарантирует, что столкновения являются чётко определёнными событиями в пространстве и времени, и выполняется, если безразмерный параметр мал, где  — диапазон межчастичных взаимодействий и  — плотность числа частиц[3]. В дополнение к этому предположению теория Чепмена — Энскога также требует, чтобы было намного меньше любых внешних временных масштабов . Это временные рамки, связанные с членами левой части уравнения Больцмана, которые описывают изменения состояния газа на макроскопических длинах. Обычно их значения определяются начальными/граничными условиями и/или внешними полями. Такое разделение масштабов означает, что столкновительный член в правой части уравнения Больцмана намного меньше, чем потоковые члены в левой части. Таким образом, приближённое решение можно найти из

Можно показать, что решение этого уравнения является гауссовым:

где  — масса молекулы и  — постоянная Больцмана[4]. Говорят, что газ находится в локальном равновесии, если он удовлетворяет этому уравнению.[5] Предположение о локальном равновесии приводит непосредственно к уравнениям Эйлера, которые описывают жидкости без диссипации, то есть с теплопроводностью и вязкостью, равными . Основная цель теории Чепмена — Энскога — систематическое получение обобщений уравнений Эйлера, учитывающих диссипацию. Это достигается путем выражения отклонений от локального равновесия в виде ряда теории возмущений по числу Кнудсена , что мало, если . Концептуально полученные гидродинамические уравнения описывают динамическое взаимодействие между свободным потоком и межчастичными столкновениями. Последние имеют тенденцию подталкивать газ к локальному равновесию, тогда как первые воздействуют на пространственные неоднородности, выводя газ из локального равновесия[6]. Когда число Кнудсена порядка 1 и более, газ в рассматриваемой системе нельзя назвать жидкостью.

Для первого порядка по получаются уравнения Навье — Стокса. Второй и третий порядки приводят соответственно к уравнениям Барнетта и суперуравнениям Барнетта.

Математическая формулировка

[править | править код]

Поскольку число Кнудсена появляется в уравнении Больцмана не явно, а неявно в терминах функции распределения и граничных условий, фиктивная переменная вводится для отслеживания соответствующих порядков разложения Чепмена — Энскога

Маленькое значение подразумевает, что столкновительный член доминирует над потоками , что аналогично утверждению о малости числа Кнудсена. Таким образом, подходящей формой разложения Чепмена — Энскога является

Решения, которые можно формально разложить таким образом, известны как нормальные решения уравнения Больцмана[7]. Этот класс решений исключает непертурбативные вклады (такие как ), которые возникают в пограничных слоях или вблизи внутренних ударных слоёв. Таким образом, метод Чепмена — Энскога ограничивается ситуациями, в которых такими решениями можно пренебречь.

Подставив это разложение и приравняв коэффециенты при соотстветствующих степенях приводит к иерархии

где  — интегральный оператор, линейный по обоим аргументам, удовлетворяющий условию и . Решение первого уравнения является гауссовским:

для некоторых функций , , и . Выражение для предполагает связь между этими функциями и физическими гидродинамическими полями, определяемыми как моменты :

Однако эти два набора функций не обязательно одинаковы для (для они равны по определению). Действительно, систематически продвигаясь по иерархии, можно обнаружить, что аналогично , каждый также содержит произвольные функции и связь которых с физическими гидродинамическими полями априори неизвестна. Одно из ключевых упрощающих предположений теории Чепмена — Энскога состоит в том, чтобы предположить, что эти произвольные в остальном функции могут быть записаны через точные гидродинамические поля и их пространственные градиенты. Другими словами, пространственная и временная зависимость входит лишь неявно через гидродинамические поля. Это утверждение физически правдоподобно, поскольку малые числа Кнудсена соответствуют гидродинамическому режиму, при котором состояние газа определяется исключительно гидродинамическими полями. В случае , функции , , и предполагаются в точности равными физическим гидродинамическим полям.

Хотя эти предположения физически правдоподобны, остаётся вопрос, действительно ли существуют решения, удовлетворяющие этим свойствам. Точнее, нужно показать, что существуют решения, удовлетворяющие условиям

Более того, даже если такие решения существуют, остаётся дополнительный вопрос: охватывают ли они полный набор нормальных решений уравнения Больцмана, то есть не представляют ли они искусственное ограничение исходного разложения по . Одним из ключевых технических достижений метода Чепмена — Энскога является положительный ответ на оба этих вопроса[7]. Таким образом, по крайней мере на формальном уровне подход Чепмена — Энскога не теряет общности.

Установив эти формальные соображения, можно приступить к расчету . Результат[2]

где  — вектор и  — тензор, каждый из которых являются решением линейного неоднородного интегрального уравнения, которое можно найти явно с помощью полиномиального разложения. Здесь двоеточие обозначает двойное точечное произведение, для тензоров , .

Предсказания

[править | править код]

В первом порядке по числу Кнудсена тепловой поток установлено, что он подчиняется закону теплопроводности Фурье[8],

и тензор потока импульса  — ньютоновская жидкость[8],

где  — тождественный тензор, и  — теплопроводность и вязкость. Их можно рассчитать явно через молекулярные параметры путём решения линейного интегрального уравнения; в таблице ниже суммированы результаты для нескольких важных молекулярных моделей ( это масса молекулы и  — постоянная Больцмана)[9].

Таблица 1: Прогнозируемые выражения для теплопроводности и вязкости.
Model Notes
Жесткие упругие сферы диаметром С точностью до 3 десятичных знаков.
Молекулы с силой отталкивания обозначает Гамма-функция, а  — числовой коэффициент. Чепмен и Коулинг перечисляют несколько значений последнего, например. and .
Lennard-Jones potential: является функцией , которую можно вычислить численно. Варьируется от для до для .[10]

Из этих результатов можно получить уравнения Навье — Стокса. Учёт моментов скорости уравнения Больцмана приводит к точным уравнениям баланса гидродинамических полей , , и :

Как и в предыдущем разделе, двоеточие обозначает двойное точечное произведение, . Подставив выражения Чепмена — Энскога вместо и , приходим к уравнениям Навье — Стокса.

Сравнение с экспериментом

[править | править код]

Важным предсказанием теории Чепмена — Энскога является то, что вязкость, , не зависит от плотности (это можно увидеть для каждой молекулярной модели в таблице 1, но на самом деле это не зависит от модели). Этот противоречивый результат восходит к Джеймсу Клерку Максвеллу, который вывел его в 1860 году на основе более элементарных кинетических аргументов[11]. Это хорошо проверено экспериментально для газов обычных плотностей.

Таблица 2: Экспериментально измеренные значения для первых пяти благородных газов.[12]
Гелий 2,45
Неон 2,52
Аргон 2,48
Криптон 2,535
Ксенон 2,58

С другой стороны, теория предсказывает, что зависит от температуры. Для жёстких упругих сфер прогнозируемое масштабирование равно , в то время как другие модели обычно демонстрируют большее изменение в зависимости от температуры. Например, для молекул, отталкивающихся друг от друга с силой прогнозируемое масштабирование имеет вид , где . Принимая , соответствующий , показывает разумное согласие с экспериментально наблюдаемым масштабированием для гелия. Для более сложных газов согласие не столь хорошее, скорее всего, из-за пренебрежения силами притяжения[13]. Действительно, модель Леннарда — Джонса, включающую в себя притяжения, можно привести в большее соответствие с экспериментом (хотя и за счёт более непрозрачной модели). зависимость; см. запись Леннарда — Джонса в таблице 1)[14]. Для лучшего согласия с экспериментальными данными, чем то, которое было получено с помощью модели Леннарда — Джонса, использовался более гибкий потенциал Ми[15], дополнительная гибкость этого потенциала позволяет точно прогнозировать транспортные свойства смесей различных типов сферически-симметричных молекул.

Теория Чепмена — Энскога также предсказывает простую связь между теплопроводностью: и вязкостью, , в виде , где  — удельная теплоёмкость при постоянном объёме и  — численный фактор. Для сферически-симметричных молекул его значение, по прогнозам, будет очень близко к слегка зависящим от модели. Например, жёсткие упругие сферы имеют , а молекулы с силой отталкивания имеют (последнее отклонение в таблице 1 не учитывается). Частный случай молекул Максвелла (сила отталкивания ) имеет точно[16]. Поскольку , , и могут быть измерены непосредственно в экспериментах, простой экспериментальной проверкой теории Чепмена — Энскога является измерение для сферически-симметричных благородных газов. Таблица 2 показывает, что существует разумное согласие между теорией и экспериментом[12].

Расширения

[править | править код]

Основные принципы теории Чепмена — Энскога можно распространить на более разнообразные физические модели, включая газовые смеси и молекулы с внутренними степенями свободы. В режиме высокой плотности теория может быть адаптирована для учёта столкновительного переноса импульса и энергии, то есть переноса по диаметру молекулы во время столкновения, а не по средней длине свободного пробега (между столкновениями). Включение этого механизма предсказывает плотностную зависимость вязкости при достаточно высокой плотности, что также наблюдается экспериментально. Получение поправок, используемых для учёта переноса во время столкновения для мягких молекул (то есть молекул Леннарда — Джонса или Ми), в целом нетривиально, но успех был достигнут при применении теории возмущений Баркера — Хендерсона для точного описания этих эффектов с точностью вплоть до критической плотности различных смесей жидкостей[15].

Можно также провести теорию в более высоком порядке по числу Кнудсена. В частности, вклад второго порядка был рассчитан Барнеттом[17]. Однако в общих обстоятельствах эти поправки более высокого порядка могут не дать надёжных улучшений теории первого порядка из-за того, что разложение Чепмена — Энскога не всегда сходится[18]. С другой стороны, считается, что разложение, по крайней мере, асимптотично для решений уравнения Больцмана, и в этом случае усечение на низком уровне всё ещё даёт точные результаты[19]. Даже если поправки более высокого порядка действительно позволяют улучшить для данной системы, до сих пор ведутся споры по поводу интерпретации соответствующих гидродинамических уравнений[20].

Пересмотренная теория Энскога

[править | править код]

Распространение теории Чепмена — Энскога для многокомпонентных смесей на повышенные плотности, в частности, плотности, при которых кообъем смеси не пренебрежимо мал, было осуществлено в серии работ Э. Г. Д. Коэна и других[21][22][23][24][25] и была придумана пересмотренная теория Энскога (RET). Успешный вывод RET последовал за несколькими предыдущими попытками получить тот же результат, но которые дали несовместимые с неравновесной термодинамикой выводы. Отправной точкой для разработки RET является модифицированная форма уравнения Больцмана для -функции распределения частиц по скоростям,

где  — скорость частиц видов , на позиции и время ,  — масса частицы,  — внешняя сила, а столкновительные слагаемые

Отличие этого уравнения от классической теории Чепмена — Энскога заключается в операторе потока , в рамках которого оценивается распределение скоростей двух частиц в разных точках пространства, разделенных , где  — единичный вектор вдоль линии, соединяющей центры масс двух частиц. Ещё одно существенное отличие связано с введением факторов , которые представляют повышенную вероятность столкновений из-за исключённого объёма. Классические уравнения Чепмена — Энскога восстанавливаются, при заменах и .

Важным моментом для успеха RET является выбор факторов , что интерпретируется как функция парного распределения, оцениваемая на расстоянии контакта . Здесь следует отметить важное замечание: для получения результатов, согласующихся с неравновесной термодинамикой, следует рассматривать как функционалы полей плотности, а не как функции локальной плотности.

Результаты пересмотренной теории Энскога

[править | править код]

Одним из первых результатов, полученных с помощью RET, который отличается от результатов классической теории Чепмена — Энскога, является уравнение состояния. В то время как из классической теории Чепмена — Энскога восстанавливается закон идеального газа, RET, разработанный для жёстких упругих сфер, даёт уравнение давления

которое согласуется с уравнением состояния Карнахана — Старлинга и сводится к закону идеального газа в пределе бесконечного разбавления (то есть когда )

Для коэффициентов переноса: вязкости, теплопроводности, диффузии и термодиффузии RET предоставляет выражения, которые точно сводятся к выражениям, полученным из классической теории Чепмена — Энскога в пределе бесконечного разбавления. Однако RET предсказывает зависимость теплопроводности от плотности, которую можно выразить в виде

где и  — относительно слабые функции состава, температуры и плотности, а  — теплопроводность, полученная из классической теории Чепмена — Энскога.

Аналогично полученное выражение для вязкости можно записать как

с параметрами и в виде слабо зависящих функций состава, температуры и плотности, а также значения, полученного из классической теории Чепмена — Энскога.

Для коэффициентов диффузии и коэффициентов термодиффузии картина несколько сложнее. Однако одним из основных преимуществ RET по сравнению с классической теорией Чепмена — Энскога является то, что предсказывается зависимость коэффициентов диффузии от термодинамических факторов, то есть производных химических потенциалов по составу. Кроме того, RET не предсказывает строгой пропорциональности

для всех плотностей, а скорее предсказывает, что коэффициенты будут уменьшаться медленнее с плотностью при высоких плотностях, что хорошо согласуется с экспериментами. Эти модифицированные зависимости плотности также позволяют RET предсказать зависимость коэффициента термодиффузии от плотности:

,

в то время как классическая теория Чепмена — Энскога предсказывает, что этот коэффициент, как вязкость и теплопроводность, не зависят от плотности.

Приложения

[править | править код]

Хотя пересмотренная теория Энскога даёт много преимуществ по сравнению с классической теорией Чепмена — Энскога, ценой этого является то, что её значительно сложнее применять на практике. В то время как классическая теория Чепмена — Энскога может быть применена к сколь угодно сложным сферическим потенциалам при наличии достаточно точных и быстрых процедур интегрирования для оценки требуемых интегралов столкновений, пересмотренная теория Энскога, в дополнение к этому, требует знания контактного значения парной функции распределения.

Для смесей газов в моделях твёрдых сфер эту величину можно вычислить без больших затруднений, но для более сложных межмолекулярных потенциалов получить её вообще нетривиально. Однако определённый успех был достигнут при оценке контактного значения парной функции распределения для жидкостей Ми (которая состоит из частиц, взаимодействующих посредством обобщённого потенциала Леннарда — Джонса) и использовании этих оценок для прогнозирования транспортных свойств плотных газовых смесей и сверхкритических жидкостей[15].

Применение RET к частицам, взаимодействующим посредством реалистичных потенциалов, также ставит перед проблемой определения разумного «диаметра контакта» для мягких частиц. Хотя они однозначно определены для твёрдых сфер, до сих пор не существует общепринятого значения, которое следует использовать для контактного диаметра мягких частиц.

Примечания

[править | править код]
  1. Чепмена — Энскога метод // Физическая энциклопедия : [в 5 т.] / Гл. ред. А. М. Прохоров. — М.: Большая российская энциклопедия, 1999. — Т. 5: Стробоскопические приборы — Яркость. — 692 с. — 20 000 экз. — ISBN 5-85270-101-7.
  2. 1 2 Chapman, Sydney; Cowling, T.G. (1970), The Mathematical Theory of Non-Uniform Gases (3rd ed.), Cambridge University Press
  3. Balescu, Radu (1975), Equilibrium and Nonequilibrium Statistical Mechanics, John Wiley & Sons, ISBN 978-0-471-04600-4
  4. Cercignani, Carlo (1975), Theory and Application of the Boltzmann Equation, Elsevier, pp. 78—79, ISBN 978-0-444-19450-3
  5. Balescu, p. 450
  6. Balescu, p. 451
  7. 1 2 Grad, Harold (1958), "Principles of the Kinetic Theory of Gases", in Flügge, S. (ed.), Encyclopedia of Physics, vol. XII, Springer-Verlag, pp. 205—294
  8. 1 2 Bird, R. Bryon; Armstrong, Robert C.; Hassager, Ole (1987), Dynamics of Polymeric Liquids, Volume 1: Fluid Mechanics (2nd ed.), John Wiley & Sons, pp. 10—11
  9. Chapman & Cowling, chapter 10
  10. Chapman & Cowling, p. 185
  11. Maxwell, James (1860), "V. Illustrations of the dynamical theory of gases.—Part I. On the motions and collisions of perfectly elastic spheres", Philosophical Magazine, 19 (124): 19—32, doi:10.1080/14786446008642818
  12. 1 2 Chapman & Cowling p. 249
  13. Chapman & Cowling, pp. 230—232
  14. Chapman & Cowling, pp. 235—237
  15. 1 2 3 Jervell, Vegard G.; Wilhelmsen, Øivind (2023-06-08). "Revised Enskog theory for Mie fluids: Prediction of diffusion coefficients, thermal diffusion coefficients, viscosities, and thermal conductivities". The Journal of Chemical Physics. 158 (22). doi:10.1063/5.0149865. ISSN 0021-9606.
  16. Chapman & Cowling, pp. 247
  17. Burnett, D. (1936), "The Distribution of Molecular Velocities and the Mean Motion in a Non-Uniform Gas", Proceedings of the London Mathematical Society, 40: 382, doi:10.1112/plms/s2-40.1.382
  18. Santos, Andres; Brey, J. Javier; Dufty, James W. (1986), "Divergence of the Chapman–Enskog Expansion", Physical Review Letters, 56 (15): 1571—1574, Bibcode:1986PhRvL..56.1571S, doi:10.1103/PhysRevLett.56.1571, PMID 10032711
  19. Grad, Harold (1963), "Asymptotic Theory of the Boltzmann Equation", The Physics of Fluids, 6 (2): 147, Bibcode:1963PhFl....6..147G, doi:10.1063/1.1706716
  20. García-Cólin, L.S.; Velasco, R.M.; Uribe, F.J. (2008), "Beyond the Navier–Stokes equations: Burnett hydrodynamics", Physics Reports, 465 (4): 149—189, Bibcode:2008PhR...465..149G, doi:10.1016/j.physrep.2008.04.010
  21. López de Haro, M.; Cohen, E. G. D.; Kincaid, J. M. (1983-03-01). "The Enskog theory for multicomponent mixtures. I. Linear transport theory". The Journal of Chemical Physics. 78 (5): 2746—2759. doi:10.1063/1.444985. ISSN 0021-9606.
  22. Kincaid, J. M.; López de Haro, M.; Cohen, E. G. D. (1983-11-01). "The Enskog theory for multicomponent mixtures. II. Mutual diffusion". The Journal of Chemical Physics (англ.). 79 (9): 4509—4521. doi:10.1063/1.446388. ISSN 0021-9606.
  23. López de Haro, M.; Cohen, E. G. D. (1984-01-01). "The Enskog theory for multicomponent mixtures. III. Transport properties of dense binary mixtures with one tracer component". The Journal of Chemical Physics (англ.). 80 (1): 408—415. doi:10.1063/1.446463. ISSN 0021-9606.
  24. Kincaid, J. M.; Cohen, E. G. D.; López de Haro, M. (1987-01-15). "The Enskog theory for multicomponent mixtures. IV. Thermal diffusion". The Journal of Chemical Physics (англ.). 86 (2): 963—975. doi:10.1063/1.452243. ISSN 0021-9606.
  25. Van Beijeren, H.; Ernst, M.H. (March 1973). "The non-linear Enskog-Boltzmann equation". Physics Letters A (англ.). 43 (4): 367—368. doi:10.1016/0375-9601(73)90346-0.

Литература

[править | править код]
  • Хуанг К.[англ.]. Метод Чепмена — Энскога // Статистическая механика. — М.: Мир, 1966. — С. 142—153. — 520 с.