Уравнение Больцмана
Уравнение Больцмана или уравнение переноса Больцмана ( БТЕ ) описывает статистическое поведение термодинамической системы, не находящейся в состоянии равновесия ; он был разработан Людвигом Больцманом в 1872 году. [ 2 ] Классическим примером такой системы является жидкость с температурными градиентами в пространстве, заставляющими тепло перетекать из более горячих областей в более холодные за счет случайного, но смещенного переноса частиц , составляющих эту жидкость. В современной литературе термин «уравнение Больцмана» часто используется в более общем смысле, имея в виду любое кинетическое уравнение, описывающее изменение макроскопической величины в термодинамической системе, такой как энергия, заряд или число частиц.
Уравнение возникает не в результате анализа отдельных положений и импульсов каждой частицы в жидкости, а в результате рассмотрения распределения вероятностей положения и импульса типичной частицы, то есть вероятности того , что частица займет данную очень небольшую область пространства. (математически элемент объема ) по центру позиции , и имеет импульс, почти равный заданному вектору импульса (таким образом, занимая очень небольшую область импульсного пространства ), в момент времени.
Уравнение Больцмана можно использовать для определения того, как изменяются физические величины, такие как тепловая энергия и импульс , когда жидкость перемещается. Можно также получить другие свойства, характерные для жидкостей, такие как вязкость , теплопроводность и электропроводность (рассматривая носители заряда в материале как газ). [ 2 ] См. также уравнение конвекции-диффузии .
Уравнение представляет собой нелинейное интегро-дифференциальное уравнение , а неизвестная функция в уравнении представляет собой функцию плотности вероятности в шестимерном пространстве положения и импульса частицы. Проблема существования и единственности решений до сих пор до конца не решена, но некоторые недавние результаты весьма многообещающи. [ 3 ] [ 4 ]
Обзор
[ редактировать ]Фазовое пространство и функция плотности
[ редактировать ]Множество всех возможных положений r и импульсов p называется фазовым пространством системы; другими словами, набор из трех для каждой координаты положения x, y, z и еще трех для каждой компоненты импульса p x , py координат , p z . Все пространство шестимерно : точка в этом пространстве равна ( r , p ) = ( x, y, z, p x , p y , p z ) , и каждая координата параметризована временем t . Малый объем (« элемент дифференциального объема ») пишется
Поскольку вероятность существования N молекул, все из которых имеют r и p в пределах В основе уравнения лежит величина f, которая дает эту вероятность на единицу объема фазового пространства или вероятность на единицу длины в кубе на единицу импульса в кубе в момент времени t . Это функция плотности вероятности : f ( r , p , t ) , определенная так, что - это количество молекул, позиции которых лежат внутри элемента объема. о r и импульсах, лежащих внутри импульсного пространства элемента около р , в момент времени т . [ 5 ] Интегрирование по области пространства позиций и пространства импульсов дает общее количество частиц, которые имеют положения и импульсы в этой области:
что является 6-кратным интегралом . Хотя f связано с несколькими частицами, фазовое пространство предназначено для одной частицы (не для всех из них, что обычно имеет место в детерминированных системах многих тел ), поскольку только об одном r и p речь идет не является частью анализа . Использование r 1 , p 1 для частицы 1, r 2 , p 2 для частицы 2 и т. д. вплоть до r N , p N для частицы N .
Предполагается, что частицы в системе идентичны (поэтому каждая имеет одинаковую массу m ). Для смеси более чем одного химического вещества необходимо одно распределение для каждого, см. ниже.
Основное заявление
[ редактировать ]Тогда общее уравнение можно записать в виде [ 6 ]
где термин «сила» соответствует силам, действующим на частицы внешним воздействием (а не самими частицами), термин «разность» представляет собой диффузию частиц , а термин «колл» - это термин столкновения , учитывающий силы действие между частицами при столкновении. Ниже приведены выражения для каждого термина в правой части. [ 6 ]
Обратите внимание, что некоторые авторы используют скорость частицы v вместо импульса p ; они связаны в определении импульса соотношением p = m v .
Силовые и диффузионные условия
[ редактировать ]Рассмотрим частицы, описываемые f , каждая из которых испытывает внешнюю силу F, не связанную с другими частицами (см. термин столкновения для последней трактовки).
Предположим, что в момент времени t некоторое количество частиц находится в позиции r внутри элемента. и импульс p внутри . Если на каждую частицу мгновенно действует сила F , то в момент времени t + Δt их положение будет и импульс p + Δ p знак равно p + F Δ t . Тогда при отсутствии столкновений f должно удовлетворять
Обратите внимание, что мы использовали тот факт, что элемент объема фазового пространства является константой, что можно показать с помощью уравнений Гамильтона (см. обсуждение теоремы Лиувилля ). Однако поскольку столкновения все же происходят, плотность частиц в фазовом объеме меняется, поэтому
| ( 1 ) |
где Δf — общее изменение f . Разделив ( 1 ) на пределам ∆t → → 0 и ∆f и переходя к 0 , имеем
| ( 2 ) |
Полный дифференциал f : равен
| ( 3 ) |
где ∇ — оператор градиента , · — скалярное произведение , является сокращением аналога импульса ∇ , а ê 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 ] В приближении БГК предполагается, что эффект молекулярных столкновений заключается в том, чтобы заставить неравновесную функцию распределения в точке физического пространства вернуться к максвелловской равновесной функции распределения, и что скорость, с которой это происходит, пропорциональна частоте молекулярных столкновений. . Таким образом, уравнение Больцмана преобразуется в форму БГК:
где - частота столкновений молекул, а — локальная максвелловская функция распределения, заданная при температуре газа в этой точке пространства. Это также называется «аппроксимацией времени релаксации».
Общее уравнение (для смеси)
[ редактировать ]Для смеси химических веществ, обозначенных индексами i = 1, 2, 3, ..., n, уравнение для видов i имеет вид [ 2 ]
где ж я знак равно ж я ( р , п я , т ) , а член столкновения
где f' = f' ( p' i , t ) , величина относительных импульсов равна
I частицами ij — дифференциальное сечение, как и прежде, между i и j . Интегрирование производится по компонентам импульса в подынтегральном выражении (которые обозначены i и j ). Сумма интегралов описывает вход и выход частиц вида i в элемент фазового пространства или из него.
Приложения и расширения
[ редактировать ]Уравнения сохранения
[ редактировать ]Уравнение Больцмана можно использовать для вывода гидродинамических законов сохранения массы, заряда, импульса и энергии. [ 8 ] : 163 Для жидкости, состоящей только из одного вида частиц, плотность числа n определяется выражением
Среднее значение любой функции A равно
Поскольку в уравнениях сохранения используются тензоры, будет использоваться соглашение Эйнштейна о суммировании, где повторяющиеся индексы в произведении указывают на суммирование по этим индексам. Таким образом и , где – вектор скорости частицы. Определять как некоторая функция импульса единственное, которое сохраняется при столкновении. Предположим также, что сила является функцией только положения, и что f равно нулю для . Умножение уравнения Больцмана на A и интегрирование по импульсу дает четыре члена, которые с помощью интегрирования по частям можно выразить как
где последний член равен нулю, поскольку A сохраняется при столкновении. Значения A соответствуют моментам скорости (и импульс , поскольку они линейно зависимы).
Нулевой момент
[ редактировать ]Сдача в аренду , масса частицы, интегрированное уравнение Больцмана становится уравнением сохранения массы: [ 8 ] : 12, 168 где - массовая плотность, а – средняя скорость жидкости.
Первый момент
[ редактировать ]Сдача в аренду , импульс частицы, интегрированное уравнение Больцмана становится уравнением сохранения импульса: [ 8 ] : 15, 169
где – тензор давления ( тензор вязких напряжений плюс гидростатическое давление ).
Второй момент
[ редактировать ]Сдача в аренду , кинетическая энергия частицы, интегрированное уравнение Больцмана становится уравнением сохранения энергии: [ 8 ] : 19, 169
где - плотность кинетической тепловой энергии, а – вектор теплового потока.
гамильтонова механика
[ редактировать ]В гамильтоновой механике уравнение Больцмана часто записывается в более общем виде как где L — оператор Лиувилля (существует противоречивое определение между оператором Лиувилля, как он определен здесь, и определением в статье, на которую есть ссылка), описывающим эволюцию объема фазового пространства, а C — оператор столкновений. Нерелятивистская форма L есть
Квантовая теория и нарушение закона сохранения числа частиц
[ редактировать ]Можно записать релятивистские квантовые уравнения Больцмана для релятивистских квантовых систем, в которых число частиц не сохраняется при столкновениях. Это имеет несколько применений в физической космологии . [ 9 ] включая образование легких элементов в ходе нуклеосинтеза Большого взрыва , производство темной материи и бариогенеза . Априори не ясно, можно ли охарактеризовать состояние квантовой системы классической плотностью фазового пространства f . Однако для широкого класса приложений существует четко определенное обобщение f , которое является решением эффективного уравнения Больцмана, которое можно вывести из первых принципов квантовой теории поля . [ 10 ]
Общая теория относительности и астрономия
[ редактировать ]Уравнение Больцмана используется в галактической динамике. Галактику при определенных предположениях можно представить как непрерывную жидкость; его массовое распределение тогда представляется f ; в галактиках физические столкновения между звездами очень редки, а эффектом гравитационных столкновений можно пренебречь в течение времени, намного превышающего возраст Вселенной .
Его обобщение в общей теории относительности есть [ 11 ] где Γ а βγ — символ Кристоффеля второго рода (при этом предполагается отсутствие внешних сил, так что частицы движутся по геодезическим в отсутствие столкновений), с той важной тонкостью, что плотность является функцией смешанного контравариантно-ковариантного ( x я , pi ) контравариантного фазовое пространство в отличие от полностью ( x я , п я ) фазовое пространство. [ 12 ] [ 13 ]
В физической космологии полностью ковариантный подход использовался для изучения космического микроволнового фонового излучения. [ 14 ] В более общем плане изучение процессов в ранней Вселенной часто пытается принять во внимание эффекты квантовой механики и общей теории относительности . [ 9 ] В очень плотной среде, образованной первичной плазмой после Большого взрыва , частицы непрерывно создаются и уничтожаются. В такой среде квантовая когерентность и пространственное расширение волновой функции могут влиять на динамику, ставя под сомнение, подходит ли классическое распределение f в фазовом пространстве , которое появляется в уравнении Больцмана, для описания системы. Однако во многих случаях возможно вывести эффективное уравнение Больцмана для обобщенной функции распределения из первых принципов квантовой теории поля . [ 10 ] Это включает в себя образование легких элементов в ходе нуклеосинтеза Большого взрыва , производство темной материи и бариогенез .
Решение уравнения
[ редактировать ]Доказано, что в некоторых случаях существуют точные решения уравнений Больцмана; [ 15 ] этот аналитический подход дает понимание, но обычно непригоден для решения практических задач.
Вместо этого численные методы (включая методы конечных элементов и решеточные методы Больцмана для поиска приближенных решений различных форм уравнения Больцмана обычно используются ). Примеры применения варьируются от гиперзвуковой аэродинамики в потоках разреженного газа. [ 16 ] [ 17 ] плазменным потокам. [ 18 ] Приложением уравнения Больцмана в электродинамике является расчет электропроводности - результат в ведущем порядке идентичен квазиклассическому результату. [ 19 ]
Вблизи локального равновесия решение уравнения Больцмана может быть представлено асимптотическим разложением по степеням числа Кнудсена ( Чепмена–Энскога разложение [ 20 ] ). Первые два члена этого разложения дают уравнения Эйлера и уравнения Навье – Стокса . Высшие члены имеют особенности. Проблема математического развития предельных процессов, ведущих от атомистической точки зрения (представленной уравнением Больцмана) к законам движения континуумов, является важной частью шестой проблемы Гильберта . [ 21 ]
Ограничения и дальнейшее использование уравнения Больцмана
[ редактировать ]Уравнение Больцмана справедливо только при нескольких предположениях. Например, предполагается, что частицы точечны, т.е. не имеют конечного размера. Существует обобщение уравнения Больцмана, называемое уравнением Энскога . [ 22 ] Член столкновения изменен в уравнениях Энскога таким образом, что частицы имеют конечный размер, например, их можно моделировать как сферы с фиксированным радиусом.
Для частиц не предполагается никаких других степеней свободы, кроме поступательного движения. Если существуют внутренние степени свободы, уравнение Больцмана должно быть обобщенным и может иметь неупругие столкновения . [ 22 ]
Многие реальные жидкости, такие как жидкости или плотные газы, помимо упомянутых выше особенностей имеют более сложные формы столкновений, будут не только бинарные, но также тройные и более высокие порядки столкновений. [ 23 ] Они должны быть получены с использованием иерархии BBGKY .
Уравнения типа Больцмана также используются для движения клеток . [ 24 ] [ 25 ] Поскольку клетки представляют собой составные частицы , обладающие внутренними степенями свободы, соответствующие обобщенные уравнения Больцмана должны иметь интегралы неупругих столкновений. Такие уравнения могут описывать инвазию раковых клеток в ткани, морфогенез и хемотаксисом эффекты, связанные с .
См. также
[ редактировать ]Примечания
[ редактировать ]- ^ Горбань, Александр Н.; Карлин, Илья В. (2005). Инвариантные многообразия для физической и химической кинетики . Конспект лекций по физике (ЛНП, т. 660). Берлин, Гейдельберг: Springer. дои : 10.1007/b98103 . ISBN 978-3-540-22684-0 . Все URL
- ^ Перейти обратно: а б с д Энциклопедия физики (2-е издание), Р.Г. Лернер , Г.Л. Тригг, издатели VHC, 1991, ISBN (Verlagsgesellschaft) 3-527-26954-1, ISBN (VHC Inc.) 0-89573-752-3.
- ^ ДиПерна, Р.Дж.; Лайонс, П.-Л. (1989). «О задаче Коши для уравнений Больцмана: глобальное существование и слабая устойчивость». Энн. математики . 2. 130 (2): 321–366. дои : 10.2307/1971423 . JSTOR 1971423 .
- ^ Филип Т. Грессман и Роберт М. Стрейн (2010). «Глобальные классические решения уравнения Больцмана с дальнодействующими взаимодействиями» . Труды Национальной академии наук . 107 (13): 5744–5749. arXiv : 1002.3639 . Бибкод : 2010PNAS..107.5744G . дои : 10.1073/pnas.1001185107 . ПМЦ 2851887 . ПМИД 20231489 .
- ^ Хуанг, Керсон (1987). Статистическая механика (Второе изд.). Нью-Йорк: Уайли. п. 53 . ISBN 978-0-471-81518-1 .
- ^ Перейти обратно: а б Энциклопедия физики МакГроу Хилла (2-е издание), С. П. Паркер, 1993, ISBN 0-07-051400-3 .
- ^ Бхатнагар, Польша; Гросс, EP; Крук, М. (1 мая 1954 г.). «Модель столкновительных процессов в газах. I. Процессы малой амплитуды в заряженных и нейтральных однокомпонентных системах». Физический обзор . 94 (3): 511–525. Бибкод : 1954PhRv...94..511B . дои : 10.1103/PhysRev.94.511 .
- ^ Перейти обратно: а б с д де Гроот, СР; Мазур, П. (1984). Неравновесная термодинамика . Dover Publications Inc. Нью-Йорк: ISBN 978-0-486-64741-8 .
- ^ Перейти обратно: а б Эдвард Колб и Майкл Тернер (1990). Ранняя Вселенная . Вествью Пресс. ISBN 978-0-201-62674-2 .
- ^ Перейти обратно: а б М. Древес; К. Венигер; С. Мендисабаль (8 января 2013 г.). «Уравнение Больцмана из квантовой теории поля». Физ. Летт. Б. 718 (3): 1119–1124. arXiv : 1202.1301 . Бибкод : 2013PhLB..718.1119D . дои : 10.1016/j.physletb.2012.11.046 . S2CID 119253828 .
- ^ Элерс Дж. (1971) Общая теория относительности и космология (Варенна), Р.К. Сакс (Academic Press, Нью-Йорк); Торн К.С. (1980) Rev. Mod. Физ., 52, 299; Эллис СКФ, Трециокас Р., Матраверс ДР, (1983) Энн. Физ., 150, 487}
- ^ Дебаш, Фабрис; Виллем ван Леувен (2009). «Общее релятивистское уравнение Больцмана I: ковариантная трактовка». Физика А. 388 (7): 1079–1104. Бибкод : 2009PhyA..388.1079D . дои : 10.1016/j.physa.2008.12.023 .
- ^ Дебаш, Фабрис; Виллем ван Леувен (2009). «Общее релятивистское уравнение Больцмана II: явно ковариантная трактовка». Физика А. 388 (9): 1818–34. Бибкод : 2009PhyA..388.1818D . дои : 10.1016/j.physa.2009.01.009 .
- ^ Мартенс Р., Гебби Т., Эллис СКФ (1999). «Анизотропия космического микроволнового фона: нелинейная динамика». Физ. Преподобный Д. 59 (8): 083506
- ^ Филип Т. Грессман; Роберт М. Стрейн (2011). «Глобальные классические решения уравнения Больцмана без углового обрезания». Журнал Американского математического общества . 24 (3): 771. arXiv : 1011.5441 . дои : 10.1090/S0894-0347-2011-00697-8 . S2CID 115167686 .
- ^ Эванс, Бен; Морган, Кен; Хасан, Убе (01 марта 2011 г.). «Разрывное конечно-элементное решение кинетического уравнения Больцмана в бесстолкновительной форме и форме БГК для макроскопических потоков газа» . Прикладное математическое моделирование . 35 (3): 996–1015. дои : 10.1016/j.apm.2010.07.027 .
- ^ Эванс, Б.; Уолтон, СП (декабрь 2017 г.). «Аэродинамическая оптимизация гиперзвукового возвращаемого аппарата на основе решения уравнения Больцмана – БГК и эволюционной оптимизации» . Прикладное математическое моделирование . 52 : 215–240. дои : 10.1016/j.apm.2017.07.024 . ISSN 0307-904X .
- ^ Парески, Л.; Руссо, Г. (1 января 2000 г.). «Численное решение уравнения Больцмана I: спектрально точная аппроксимация оператора столкновения». SIAM Journal по численному анализу . 37 (4): 1217–1245. CiteSeerX 10.1.1.46.2853 . дои : 10.1137/S0036142998343300 . ISSN 0036-1429 .
- ^ HJW Мюллер-Кирстен, Основы статистической механики, глава 13, 2-е изд., World Scientific (2013), ISBN 978-981-4449-53-3 .
- ^ Сидней Чепмен; Томас Джордж Коулинг Математическая теория неоднородных газов: описание кинетической теории вязкости, теплопроводности и диффузии в газах , Cambridge University Press, 1970. ISBN 0-521-40844-X
- ^ «Тема выпуска 'Шестая проблема Гильберта' » . Философские труды Королевского общества А. 376 (2118). 2018. дои : 10.1098/rsta/376/2118 .
- ^ Перейти обратно: а б «Уравнение Энскога — обзор | Темы ScienceDirect» . www.sciencedirect.com . Проверено 10 мая 2022 г.
- ^ ван Нойе, TPC; Эрнст, МГ (3 июня 1997 г.). «Кольцевая кинетическая теория идеализированного гранулированного газа». arXiv : cond-mat/9706020 .
- ^ Шовьер, А.; Хиллен, Т.; Прециози, Л. (2007). «Моделирование движения клеток в анизотропных и гетерогенных сетевых тканях» . Американский институт математических наук . 2 (2): 333–357. дои : 10.3934/nhm.2007.2.333 .
- ^ Конте, Мартина; Лой, Надя (12 февраля 2022 г.). «Множественная кинетическая модель с нелокальным зондированием для миграции клеток в оптоволоконной сети с хемотаксисом» . Бюллетень математической биологии . 84 (3): 42. дои : 10.1007/s11538-021-00978-1 . ISSN 1522-9602 . ПМЦ 8840942 . ПМИД 35150333 .
Ссылки
[ редактировать ]- Харрис, Стюарт (1971). Введение в теорию уравнения Больцмана . Дуврские книги. п. 221. ИСБН 978-0-486-43831-3 . . Очень недорогое введение в современную структуру (начиная с формального вывода из Лиувилля и иерархии Боголюбова-Борна-Грина-Кирквуда-Ивона (BBGKY), в которую помещено уравнение Больцмана). Большинство учебников по статистической механике, таких как Хуанг, по-прежнему рассматривают эту тему, используя оригинальные аргументы Больцмана. Для вывода уравнения в этих книгах используется эвристическое объяснение, которое не раскрывает диапазон применимости и характерные допущения, которые отличают уравнение Больцмана от других уравнений переноса, таких как уравнения Фоккера-Планка или Ландау .
- Аркерид, Лейф (1972). «Об уравнении Больцмана, часть I: Существование». Арх. Рациональный механизм. Анал . 45 (1): 1–16. Бибкод : 1972АрРМА..45....1А . дои : 10.1007/BF00253392 . S2CID 117877311 .
- Аркерид, Лейф (1972). «Об уравнении Больцмана, часть II: проблема полного начального значения». Арх. Рациональный механизм. Анал . 45 (1): 17–34. Бибкод : 1972АрРМА..45...17А . дои : 10.1007/BF00253393 . S2CID 119481100 .
- Аркерид, Лейф (1972). «Об уравнении Больцмана, часть I: Существование». Арх. Рациональный механизм. Анал . 45 (1): 1–16. Бибкод : 1972АрРМА..45....1А . дои : 10.1007/BF00253392 . S2CID 117877311 .
- ДиПерна, Р.Дж.; Лайонс, П.-Л. (1989). «О задаче Коши для уравнений Больцмана: глобальное существование и слабая устойчивость». Энн. математики . 2. 130 (2): 321–366. дои : 10.2307/1971423 . JSTOR 1971423 .