Из Википедии, бесплатной энциклопедии
Перейти к навигации Перейти к поиску
Место кинетического уравнения Больцмана на ступенях редукции моделей от микроскопической динамики к макроскопической континуальной динамике (иллюстрация к содержанию книги [1] )

Уравнение Больцмана или уравнение переноса Больцмана ( BTE ) описывает статистическое поведение термодинамической системы, не находящейся в состоянии равновесия , разработанное Людвигом Больцманом в 1872 году. [2] Классическим примером такой системы является жидкость с градиентами температуры в пространстве. заставляя тепло течь из более горячих областей в более холодные, за счет случайного, но смещенного переноса частицсоставляя эту жидкость. В современной литературе термин уравнение Больцмана часто используется в более общем смысле, имея в виду любое кинетическое уравнение, которое описывает изменение макроскопической величины в термодинамической системе, такой как энергия, заряд или число частиц.

Уравнение возникает не путем анализа отдельных положений и импульсов каждой частицы в жидкости, а, скорее, путем рассмотрения распределения вероятностей для положения и импульса типичной частицы, то есть вероятности того, что частица занимает заданную очень маленькую область пространства. (математически элемент объема ) с центром в этом положении и имеет импульс, почти равный заданному вектору импульса (таким образом, занимая очень небольшую область импульсного пространства ) в момент времени.

Уравнение Больцмана можно использовать для определения того, как изменяются физические величины, такие как тепловая энергия и импульс , когда жидкость перемещается. Можно также получить другие свойства, характерные для жидкостей, такие как вязкость , теплопроводность и электропроводность (рассматривая носители заряда в материале как газ). [2] См. Также уравнение конвекции – диффузии .

Уравнение представляет собой нелинейное интегро-дифференциальное уравнение , а неизвестная функция в уравнении представляет собой функцию плотности вероятности в шестимерном пространстве положения и количества движения частицы. Проблема существования и уникальности решений до сих пор полностью не решена, но некоторые недавние результаты весьма многообещающие. [3] [4]

Обзор [ править ]

Фазовое пространство и функция плотности [ править ]

Множество всех возможных положений r и импульсов p называется фазовым пространством системы; другими словами, набор из трех координат для каждой координаты положения x, y, z и еще трех для каждой компоненты импульса p x , p y , p z . Все пространство 6- мерное : точка в этом пространстве ( r , p ) = ( x, y, z, p x , p y , p z ), и каждая координата параметризована временем t . Небольшой объем («дифференциалэлемент тома ") написано

Поскольку вероятность N молекул , которые все имеют р и р в пределах находится под вопросом, в основе уравнения является величиной F , которая дает эту вероятность на единицу объема фазового пространства, или вероятность на единицу длины в кубе на единицу импульса в кубе, в момент времени t . Это функция плотности вероятности : f ( r , p , t ), определенная так, что, 

- количество молекул, каждая из которых имеет позиции, лежащие в элементе объема около r, и импульсы, лежащие в пределах элемента пространства импульсов около p , в момент времени t . [5] Интегрирование по области позиционного пространства и импульсного пространства дает общее количество частиц, которые имеют позиции и импульсы в этой области:

который представляет собой шестикратный интеграл . В то время как f связано с несколькими частицами, фазовое пространство предназначено для одной частицы (не для всех, что обычно имеет место в детерминированных системах многих тел ), поскольку речь идет только об одном r и p . Использование r 1 , p 1 для частицы 1, r 2 , p 2 для частицы 2 и т. Д. До r N , p N для частицы N не является частью анализа .

Предполагается, что частицы в системе идентичны (т.е. каждая имеет одинаковую массу m ). Для смеси более чем одного химического вещества требуется одно распределение для каждого, см. Ниже.

Основное заявление [ править ]

Тогда общее уравнение можно записать как [6]

где термин "сила" соответствует силам, действующим на частицы внешним воздействием (не самими частицами), термин "diff" представляет собой диффузию частиц, а термин "coll" представляет собой термин столкновения - с учетом сил действующие между частицами при столкновениях. Выражения для каждого члена с правой стороны приведены ниже. [6]

Обратите внимание, что некоторые авторы используют скорость частицы v вместо импульса p ; они связаны в определении импульса соотношением p = m v .

Условия силы и диффузии [ править ]

Рассмотрим частицы, описываемые f , каждая из которых испытывает внешнюю силу F, не связанную с другими частицами (см. Термин столкновения для последней трактовки).

Предположим, что в момент t некоторое количество частиц все имеют положение r внутри элемента и импульс p внутри . Если сила F мгновенно действует на каждую частицу, то в момент времени t + Δ t их положение будет r + Δ r = r + p Δ t / m, а импульс p + Δ p = p + F Δ t . Тогда в отсутствие столкновений f должна удовлетворять

Обратите внимание, что мы использовали тот факт, что элемент объема фазового пространства является постоянным, что может быть показано с помощью уравнений Гамильтона (см. Обсуждение теоремы Лиувилля ). Однако, так как происходят столкновения, плотность частиц в объеме фазового пространства  ' изменяется, так 

где Δ f - полное изменение f . Разделив ( 1 ) на  Δ t и взяв пределы Δ t → 0 и Δ f → 0, имеем 

Общие дифференциальный из F является:

где ∇ - оператор градиента , · - скалярное произведение ,

является сокращением для импульсного аналога ∇, а ê x , ê y , ê z - декартовы единичные векторы .

Заключительное заявление [ править ]

Разделив ( 3 ) на dt и подставив в ( 2 ), получим:

В этом контексте F ( r , t ) - это силовое поле, действующее на частицы в жидкости, а m - масса частиц. Термин справа добавлен для описания эффекта столкновения между частицами; если он равен нулю, то частицы не сталкиваются. Бесстолкновительное уравнение Больцмана, в котором индивидуальные столкновения заменяются агрегированными на больших расстояниях взаимодействиями, например кулоновскими взаимодействиями , часто называют уравнением Власова .

Это уравнение более полезно, чем основное уравнение, приведенное выше, но все еще неполное, поскольку f не может быть решено, если не известен член столкновения в f . Этот термин не может быть найден так же легко или обобщенно, как другие - это статистический термин, представляющий столкновения частиц, и требует знания статистики, которой подчиняются частицы, например, распределения Максвелла – Больцмана , Ферми – Дирака или Бозе – Эйнштейна .

Термин столкновения (Stosszahlansatz) и молекулярный хаос [ править ]

Термин столкновения двух тел [ править ]

Ключевой вывод, примененный Больцманом, состоял в том, чтобы определить член столкновения, возникающий исключительно в результате столкновения двух тел между частицами, которые, как предполагается, не коррелировали до столкновения. Это предположение было названо Больцманом « Stosszahlansatz » и также известно как « предположение молекулярного хаоса ». При таком предположении член столкновения может быть записан как интеграл импульсного пространства по произведению одночастичных функций распределения: [2]

где p A и p B - импульсы любых двух частиц ( для удобства обозначены как A и B ) до столкновения, p ′ A и p ′ B - импульсы после столкновения,

- величина относительных импульсов ( подробнее об этой концепции см. в разделе относительная скорость ), а I ( g , Ω) - дифференциальное сечение столкновения, в котором относительные импульсы сталкивающихся частиц превращаются на угол θ в элемент телесного угла d Ω из-за столкновения.

Упрощения термина столкновения [ править ]

Поскольку большая часть проблемы при решении уравнения Больцмана возникает из-за сложного члена столкновения, были предприняты попытки «смоделировать» и упростить член столкновения. Наиболее известное модельное уравнение принадлежит Бхатнагару, Гроссу и Круку. [7] Предположение в приближении BGK состоит в том, что эффект молекулярных столкновений должен заставить неравновесную функцию распределения в точке физического пространства вернуться к максвелловской равновесной функции распределения, и что скорость, с которой это происходит, пропорциональна частота столкновений молекул. Таким образом, уравнение Больцмана модифицируется к форме BGK:

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

Общее уравнение (для смеси) [ править ]

Для смеси химических веществ, обозначенных индексами i = 1, 2, 3, ..., n, уравнение для вида i имеет вид [2]

где f i = f i ( r , p i , t ), а член столкновения равен

где f ′ = f ′ ( p ′ i , t ), величина относительных импульсов равна

и я IJ это дифференциальное сечение, как и раньше, между частицами я и J . Интегрирование ведется по компонентам импульса в подынтегральном выражении (которые обозначены i и j ). Сумма интегралов описывает вход и выход частиц вида i в элемент фазового пространства или из него.

Приложения и расширения [ править ]

Уравнения сохранения [ править ]

Уравнение Больцмана может использоваться для вывода законов сохранения гидродинамики для массы, заряда, импульса и энергии. [8] : p 163 Для жидкости, состоящей только из частиц одного вида, числовая плотность n определяется как

Среднее значение любой функции A равно

Поскольку уравнения сохранения включают тензоры, будет использоваться соглашение Эйнштейна о суммировании, когда повторяющиеся индексы в произведении указывают на суммирование по этим индексам. Таким образом, и , где - вектор скорости частицы. Определим как некоторую функцию только от количества движения , которое сохраняется при столкновении. Предположим также, что сила является функцией только положения и что f равно нулю для . Умножение уравнения Больцмана на A и интегрирование по импульсу дает четыре члена, которые, используя интегрирование по частям, можно выразить как

где последний член равен нулю, поскольку A сохраняется при столкновении. Допуская массу частицы, интегральное уравнение Больцмана становится уравнением сохранения массы: [8] : pp 12,168

где - массовая плотность, - средняя скорость жидкости.

Допуская импульс частицы, интегральное уравнение Больцмана становится уравнением сохранения импульса: [8] : pp 15,169

где - тензор давления (тензор вязких напряжений плюс гидростатическое давление ).

Допустим , кинетическая энергия частицы, интегральное уравнение Больцмана становится уравнением сохранения энергии: [8] : pp 19,169

где - плотность кинетической тепловой энергии, - вектор теплового потока.

Гамильтонова механика [ править ]

В гамильтоновой механике уравнение Больцмана часто записывается в более общем виде как

где L - оператор Лиувилля (существует несовместимое определение между оператором Лиувилля, как определено здесь, и оператором, указанным в статье, ссылка на который есть), описывающим эволюцию объема фазового пространства, а C - оператор столкновения. Нерелятивистская форма L имеет вид

Квантовая теория и нарушение сохранения числа частиц [ править ]

Можно записать релятивистские квантовые уравнения Больцмана для релятивистских квантовых систем, в которых число частиц не сохраняется при столкновениях. Это имеет несколько приложений в физической космологии , [9] , включая формирование световых элементов в Большом Взрыве нуклеосинтеза , производство темной материи и бариогенезиса . Априори не ясно, может ли состояние квантовой системы быть охарактеризовано классической плотностью фазового пространства f . Однако для широкого класса приложений корректно определенное обобщение fсуществует, которое является решением эффективного уравнения Больцмана, которое может быть получено из первых принципов квантовой теории поля . [10]

Общая теория относительности и астрономия [ править ]

Уравнение Больцмана полезно в галактической динамике. Галактика при определенных предположениях может быть представлена ​​как сплошная жидкость; тогда его массовое распределение представлено буквой f ; в галактиках физические столкновения между звездами очень редки, и эффектом гравитационных столкновений можно пренебречь в течение времени, намного превышающего возраст Вселенной .

Его обобщение в общей теории относительности . [11] - это

где Γ α βγ - символ Кристоффеля второго рода (это предполагает отсутствие внешних сил, так что частицы движутся по геодезическим в отсутствие столкновений), с важной тонкостью, заключающейся в том, что плотность является функцией смешанной контравариантно-ковариантной ( x i , p i ) фазовое пространство в отличие от полностью контравариантного ( x i , p i ) фазового пространства. [12] [13]

В физической космологии для изучения космического микроволнового фонового излучения использовался полностью ковариантный подход. [14] В более общем плане при изучении процессов в ранней Вселенной часто пытаются принять во внимание эффекты квантовой механики и общей теории относительности . [9] В очень плотной среде, образованной первичной плазмой после Большого взрыва , частицы непрерывно создаются и аннигилируют. В такой среде квантовая когерентность и пространственное расширение волновой функции могут влиять на динамику, что делает сомнительным, что классическое распределение f в фазовом пространствекоторое появляется в уравнении Больцмана, подходит для описания системы. Однако во многих случаях можно вывести эффективное уравнение Больцмана для обобщенной функции распределения из первых принципов квантовой теории поля . [10] Это включает образование легких элементов в нуклеосинтезе Большого взрыва , производство темной материи и бариогенез .

Решение уравнения [ править ]

Доказано, что в некоторых случаях существуют точные решения уравнений Больцмана; [15] этот аналитический подход обеспечивает понимание, но обычно не используется в практических задачах.

Вместо этого численные методы (включая методы конечных элементов и решеточные методы Больцмана ) обычно используются для поиска приближенных решений различных форм уравнения Больцмана. Примеры приложений варьируются от гиперзвуковой аэродинамики в потоках разреженного газа [16] [17] до потоков плазмы. [18] Одним из приложений уравнения Больцмана в электродинамике является вычисление электропроводности - результат в первом порядке совпадает с полуклассическим результатом. [19]

Близкое к локальному равновесию решение уравнения Больцмана может быть представлено асимптотическим разложением по степеням числа Кнудсена ( разложение Чепмена-Энскога [20] ). Первые два члена этого разложения дают уравнения Эйлера и уравнения Навье-Стокса . У высших членов есть особенности. Проблема математического развития предельных процессов, которые ведут от атомистической точки зрения (представленной уравнением Больцмана) к законам движения сплошных сред, является важной частью шестой проблемы Гильберта . [21]

См. Также [ править ]

  • Уравнение Власова
  • H-теорема
  • Уравнение Фоккера – Планка
  • Уравнение Вильямса-Больцмана
  • Уравнения Навье – Стокса
  • Уравнение Власова – Пуассона.
  • Решеточные методы Больцмана

Примечания [ править ]

  1. ^ Горбань, Александр Н .; Карлин, Илья В. (2005). Инвариантные многообразия для физической и химической кинетики . Конспект лекций по физике (LNP, vol. 660). Берлин, Гейдельберг: Springer. DOI : 10.1007 / b98103 . ISBN 978-3-540-22684-0. Альтернативный URL
  2. ^ a b c d Энциклопедия физики (2-е издание), RG Lerner, GL Trigg, VHC publishers, 1991, ISBN (Verlagsgesellschaft) 3-527-26954-1, ISBN (VHC Inc.) 0-89573-752-3.
  3. ^ ДиПерна, RJ; Львов, П.-Л. (1989). «О задаче Коши для уравнений Больцмана: глобальное существование и слабая устойчивость». Анна. математики . 2. 130 (2): 321–366. DOI : 10.2307 / 1971423 . JSTOR 1971423 . 
  4. ^ Филип Т. Грессман и Роберт М. Стрейн (2010). «Глобальные классические решения уравнения Больцмана с дальнодействующими взаимодействиями» . Труды Национальной академии наук . 107 (13): 5744–5749. arXiv : 1002,3639 . Bibcode : 2010PNAS..107.5744G . DOI : 10.1073 / pnas.1001185107 . PMC 2851887 . PMID 20231489 .  
  5. ^ Хуанг, Керсон (1987). Статистическая механика (второе изд.). Нью-Йорк: Вили. п. 53 . ISBN 978-0-471-81518-1.
  6. ^ a b Энциклопедия физики Макгроу Хилла (2-е издание), CB Parker, 1994, ISBN 0-07-051400-3 . 
  7. ^ Bhatnagar, PL; Гросс, EP; Крук, М. (1954-05-01). «Модель столкновительных процессов в газах. I. Малоамплитудные процессы в заряженных и нейтральных однокомпонентных системах». Физический обзор . 94 (3): 511–525. Bibcode : 1954PhRv ... 94..511B . DOI : 10.1103 / PhysRev.94.511 .
  8. ^ а б в г де Гроот, SR; Мазур, П. (1984). Неравновесная термодинамика . Нью-Йорк: ISBN Dover Publications Inc. 978-0-486-64741-8.
  9. ^ а б Эдвард Колб и Майкл Тернер (1990). Ранняя Вселенная . Westview Press. ISBN 9780201626742.
  10. ^ а б М. Древес; К. Венигер; С. Мендизабал (8 января 2013 г.). «Уравнение Больцмана из квантовой теории поля». Phys. Lett. B . 718 (3): 1119–1124. arXiv : 1202.1301 . Bibcode : 2013PhLB..718.1119D . DOI : 10.1016 / j.physletb.2012.11.046 . S2CID 119253828 . 
  11. Ehlers J (1971) Общая теория относительности и космология (Варенна), RK Sachs (Academic Press NY); Thorne KS (1980) Rev. Mod. Phys., 52, 299; Эллис GFR, Treciokas R, Matravers DR, (1983) Ann. Phys., 150, 487}.
  12. ^ Деббаш, Фабрис; Виллем ван Леувен (2009). "Общее релятивистское уравнение Больцмана I: ковариантное рассмотрение". Physica . 388 (7): 1079–1104. Bibcode : 2009PhyA..388.1079D . DOI : 10.1016 / j.physa.2008.12.023 .
  13. ^ Деббаш, Фабрис; Виллем ван Леувен (2009). «Общее релятивистское уравнение Больцмана II: явно ковариантная трактовка». Physica . 388 (9): 1818–34. Bibcode : 2009PhyA..388.1818D . DOI : 10.1016 / j.physa.2009.01.009 .
  14. ^ Maartens Р, Т Gebbie, Эллис СКФ (1999). «Анизотропия космического микроволнового фона: нелинейная динамика». Phys. Ред. D. 59 (8): 083506
  15. ^ Филип Т. Грессман, Роберт М. Стрейн (2011). «Глобальные классические решения уравнения Больцмана без угловой обрезки». Журнал Американского математического общества . 24 (3): 771. arXiv : 1011.5441 . DOI : 10.1090 / S0894-0347-2011-00697-8 . S2CID 115167686 . CS1 maint: uses authors parameter (link)
  16. ^ Эванс, Бен; Морган, Кен; Хасан, Убей (01.03.2011). «Разрывное решение конечных элементов кинетического уравнения Больцмана в бесстолкновительной и BGK формах для макроскопических потоков газа» . Прикладное математическое моделирование . 35 (3): 996–1015. DOI : 10.1016 / j.apm.2010.07.027 .
  17. ^ Evans, B .; Уолтон, СП (декабрь 2017 г.). «Аэродинамическая оптимизация гиперзвукового космического корабля на основе решения уравнения Больцмана – БГК и эволюционной оптимизации» . Прикладное математическое моделирование . 52 : 215–240. DOI : 10.1016 / j.apm.2017.07.024 . ISSN 0307-904X . 
  18. ^ Pareschi, L .; Руссо, Г. (2000-01-01). «Численное решение уравнения Больцмана I: спектрально-точное приближение оператора столкновений». Журнал СИАМ по численному анализу . 37 (4): 1217–1245. CiteSeerX 10.1.1.46.2853 . DOI : 10.1137 / S0036142998343300 . ISSN 0036-1429 .  
  19. ^ HJW Мюллер-Кирстен, Основы статистической механики, Глава 13, 2-е изд., World Scientific (2013), ISBN 978-981-4449-53-3 . 
  20. Сидней Чепмен; Томас Джордж Коулинг Математическая теория неоднородных газов: описание кинетической теории вязкости, теплопроводности и диффузии в газах , Cambridge University Press, 1970. ISBN 0-521-40844-X 
  21. ^ «Тема выпуска« Шестая проблема Гильберта » » . Философские труды Королевского общества А . 376 (2118). 2018. DOI : 10,1098 / RSTA / 376/2118 .

Ссылки [ править ]

  • Харрис, Стюарт (1971). Введение в теорию уравнения Больцмана . Dover Книги. п. 221. ISBN. 978-0-486-43831-3.. Очень недорогое введение в современный каркас (начиная с формального вывода из Лиувилля и иерархии Боголюбова – Борна – Грина – Кирквуда – Ивона (BBGKY), в которую помещено уравнение Больцмана). Большинство учебников статистической механики, таких как Хуанг, по-прежнему рассматривают эту тему, используя оригинальные аргументы Больцмана. Для вывода уравнения в этих книгах используется эвристическое объяснение, которое не раскрывает диапазон применимости и характерные допущения, которые отличают уравнения Больцмана от других уравнений переноса, таких как уравнения Фоккера – Планка или Ландау .
  • Аркерид, Лейф (1972). «Об уравнении Больцмана часть I: Существование». Arch. Rational Mech. Анальный . 45 (1): 1–16. Bibcode : 1972ArRMA..45 .... 1A . DOI : 10.1007 / BF00253392 . S2CID  117877311 .
  • Аркерид, Лейф (1972). «Об уравнении Больцмана, часть II: Полная начальная задача». Arch. Rational Mech. Анальный . 45 (1): 17–34. Bibcode : 1972ArRMA..45 ... 17А . DOI : 10.1007 / BF00253393 . S2CID  119481100 .
  • Аркерид, Лейф (1972). «Об уравнении Больцмана часть I: Существование». Arch. Rational Mech. Анальный . 45 (1): 1–16. Bibcode : 1972ArRMA..45 .... 1A . DOI : 10.1007 / BF00253392 . S2CID  117877311 .
  • ДиПерна, Р.Дж.; Львов, П.-Л. (1989). «О задаче Коши для уравнений Больцмана: глобальное существование и слабая устойчивость». Анна. математики . 2. 130 (2): 321–366. DOI : 10.2307 / 1971423 . JSTOR  1971423 .

Внешние ссылки [ править ]

  • Уравнение переноса Больцмана Франца Веселы
  • Больцмановское поведение газов решено