Физическая кинетика

Материал из свободной русской энциклопедии «Традиция»
Перейти к навигации Перейти к поиску

Физическая кинетика (греч. κίνησις — движение) — микроскопическая теория процессов в неравновесных средах. В кинетике методами квантовой или классической статистической физики изучают процессы переноса энергии, импульса, заряда и вещества в различных физических системах (газах, плазме, жидкостях, твёрдых телах) и влияние на них внешних полей. В отличие от термодинамики неравновесных процессов и электродинамики сплошных сред, кинетика исходит из представления о молекулярном строении рассматриваемых сред, что позволяет вычислить из первых принципов кинетические коэффициенты, диэлектрические и магнитные проницаемости и другие характеристики сплошных сред. Физическая кинетика включает в себя кинетическую теорию газов из нейтральных атомов или молекул, статистическую теорию неравновесных процессов в плазме, теорию явлений переноса в твёрдых телах (диэлектриках, металлах и полупроводниках) и жидкостях, кинетику магннитных процессов и теорию кинетических явлений, связанных с прохождением быстрых частиц через вещество. К ней же относятся теория процессов переноса в квантовых жидкостях и сверхпроводниках и кинетика фазовых переходов.

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

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

Кинетическое уравнение[править | править код]

Основной метод физической кинетики — решение кинетического уравнения Больцмана для одночастичной функциии распределения f ( x , p , t ) f(x, p, t) молекул в фазовом пространстве их координат x x и импульсов p p . Функция распределения удовлетворяет кинетическому, уравнению: f t + p m f x + F f p = S t f \frac{\partial f}{\partial t}+ \frac{\vec{p}}{m} \frac{\partial f}{\partial \vec{x}} + \vec{F} \frac{\partial f}{\partial \vec{p}}= St f

,где S t St — интеграл столкновений, определяющий разность числа частиц, приходящих в элемент объёма вследствие прямых столкновений и убывающих из него вследствие обратных столкновений. Для одноатомных молекул или для многоатомных, но без учёта их внутренних степеней свободы

S t f = ω ( f f 1 f f 1 ) d p 1 d p d p 1 St f = \int{ \omega \cdot (f' f'_1 - f f_1) dp_1 dp' dp'_1 }

,где ω \omega — вероятность столкновения, связанная с дифференциальным эффективным сечением рассеяния. ω d p d p 1 = | v v 1 | d σ \omega dp' dp'_1 = | v - v_1 | d\sigma

,где p p , p 1 p_1 — импульсы молекул до столкновения, v v , v 1 v_1 — соответственно скорости, p p' , p p — их импульсы после столкновения, f f f 1 f_1 — функции распределения молекул до столкновения, f f ' , f 1 f'_1 — их функции распределения после столкновения.

Для газа из сложных молекул, обладающих внутренними степенями свободы, их следует учитывать в функции распределения. Например, для двухатомных молекул с собственным моментом вращения M функции распределения будут зависеть также от M.

Из кинетического уравнени следует теорема Больцмана — убывание со временем H-функции Больцмана (среднего логарифма функции распределения) или возрастание энтропии, так как она равна H-функции Больцмана с обратным знаком.

Уравнения переноса[править | править код]

Физическая кинетика позволяет получить уравнения баланса средней плотностей вещества, импульса и энергии. Например, для простого газа плотность ρ \rho , гидродинамическая скорость V и средняя энергия E ¯ \bar{E} в удовлетворяют уравнениям баланса: ρ t + d i v ( ρ V ) = 0 \frac{\partial \rho}{\partial t} + div (\rho V) = 0 t ( ρ V α ) + β Π α β x β = 0 \frac{\partial }{\partial t} (\rho V_\alpha) + \sum_{\beta} {\frac{\partial \Pi_{\alpha \beta}}{\partial x_\beta}}= 0 t n E ¯ + d i v ( q ) = 0 \frac{\partial}{\partial t} n \bar{E} + div (q) = 0 Π α β = m V α V β f d p \Pi_{\alpha \beta} = \int {m V_\alpha V_\beta f dp}

Π α β \Pi_{\alpha \beta} — тензор плотности потока импульса, m — масса частиц, n — плотность числа частиц, q = E V f d p q=\int{E V f dp} — плотность потока энергии.

Если состояние газа мало отличается от равновесного, то в малых элементах объёма устанавливается распределение, близкое к локально равновесному распределению Максвелла, с температурой, плотностью и гидродинамической скоростью, соответствующими рассматриваемой точке газа. В этом случае неравновесная функция распределения мало отличается от локально равновесной и решение кинетического уравнения даёт малую поправку к последней, пропорциональную градиентам температуры T \nabla T и гидродинамичой скорости V \nabla V , так как St f0=0.

С помощью неравновесной функции распределения можно найти поток энергии (в неподвижной жидкости) q = λ T q = - \lambda \nabla T , где λ \lambda — коэффициент теплопроводности, и тензор плотности потока импульса Π α β = ρ V α V β + δ α β P σ α β , \Pi_{\alpha \beta} = \rho V_{\alpha} V_{\beta} + \delta _ {\alpha \beta} P - \sigma'_{\alpha \beta}, где σ α β = μ [ ( V α V β + V β V α ) 2 3 δ α β d i v ( V ) ] \sigma'_{\alpha \beta} = \mu \left[ \left( \frac{\partial V_{\alpha}}{\partial V_{\beta}} + \frac{\partial V_{\beta}}{\partial V_{\alpha}} \right) - \frac{2}{3}\delta_{\alpha \beta} div (V)\right] — тензор вязких напряжении, η \eta — коэффициент сдвиговой вязкости, P — давление. Для газов с внутренними степенями свободы σ α β \sigma'_{\alpha \beta} содержит также член ζ δ α β \zeta \delta _{\alpha \beta} , где ζ \zeta — коэффициент «второй», объёмной вязкости, проявляющейся лишь при движениях, в которых d i v V 0 div V \ne 0 . Для кинетических коэффициентов λ \lambda , η \eta , ζ \zeta получаются выражения через эффективные сечения столкновений и, следовательно, через константы молекулярных взаимодействий. В бинарной смеси поток вещества состоит из диффузионного потока, пропорционального градиенту концентрации вещества в смеси с коэффициентом диффузии, и термодиффузионного потока, пропорционального градиенту температуры с коэффициентом термодиффузии, а поток тепла, кроме обычного члена теплопроводности, пропорционального градиенту температуры, содержит дополнительный, член, пропорциональный градиенту концентрации и описывающий эффект Дюфура. Кинетика даёт выражения для этих кинетических коэффициентов через эффективные сечения столкновений, кинетические коэффициенты для перекрёстных явлений, например термодиффузии и эффекта Дюфура, вследсвие теоремы Онсагера оказываются равными. Эти соотношения являются следствием микроскопической обратимости уравнений движения частиц системы, то есть инвариантности их относительно обращения времени.

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

Газы и плазма[править | править код]

Физическая кинетика позволяет исследовать явления переноса в разреженных газах, когда отношение длины свободного пробега λ к характерным размерам задачи L (то есть число Кнудсена λ/L) уже не очень малО и имеет смысл рассматривать поправки порядка 1/L (слабо разреженные газы). В этом случае кинетика объясняет явления температурного скачка и течения газов вблизи твёрдых поверхностей.

Для сильно разреженных газов, когда λ/L>1, гидродинамические уравнения и обычное уравнение теплопроводности уже не применимы и для исследования процессов переноса необходимо решать кинетическок уравнение с определёнными граничными условиями на поверхностях, ограничивающих газ. Эти условия выражаются через функцию распределения молекул, рассеянных из-за взаимодействия со стенкой. Рассеянный поток частиц может приходить в тепловое равновесие со стенкой, но в реальных случаях это не достигается. Для сильно разреженных газов роль коэффициента теплопроводности играют коэффициенты теплопередачи. Например, колчество тепла Q, отнесённое к единице площади параллельных пластинок, между которыми находится разреженный газ, равно Q = κ ( T 2 T 1 ) / L Q = \kappa (T_2 - T_1)/L , где T 1 T_1 и T 2 T_2 — температуры пластинок, L — расстояние между ними, κ \kappa — коэффициент теплопередачи.

Теория явлений переноса в плотных газах и жидкостях значительно сложнее, так как для описания неравновесного состояния уже недостаточно одночастичной функции распределения, а нужно учитывать функции распределения более высокого порядка. Частичные функции распределения удовлетворяют цепочке зацепляющихся уравнений (так называемых уравнений Боголюбова или цепочкой ББГКИ, т. е, уравнений Боголюбова—Борна—Грина—Кирквуда—Ивона). С помощью этих уравнений можно уточнить кинетическое уравнение для газов средней плотности и исследовать для них явления переноса.

Физическая кинетика двухкомпонентнои плазмы описывается двумя функциями распределения (для электронов f e f_e , для ионов f i f_i ), удовлетворяющими системе двух кинетических уравнений. На частицы плазмы действуют силы

F e = e ( E + v × B c ) F_e = -e (E + \frac {v \times B}{c}) , F i = Z F e F_i = -Z F_e


где Z e Z_e — заряд иона, E — напряжённость электрического поля, B — магнитная индукция, удовлетворяющие уравнениям Максвелла. Уравнения Максвелла содержат среднии плотности тока j и заряда ρ, определяемые с помощью функций распределения:

j = e v ( Z f i f e ) d p j = e \int {v (Z f_i - f_e) dp} , p = e ( Z f i f e ) d p p = e \int {(Z f_i - f_e) dp}

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

Конденсированные среды[править | править код]

Физическая кинетика неравновесных процессов в диэлектриках основана на решении кинетического уравнения Больцмана для фононов решётки. Взаимодействие между фононами вызвано ангармоническими членами гамильтониана решётки относительно смещения атомов из положения равновесия. При простейших столкновениях один фонон распадается на два или происходит слияние двух фононов в один, причём сумма их квазиимпульсов либо сохраняется (нормальные процессы столкновений), либо меняется на вектор обратной решётки (процессы переброса). Конечная теплопроводность возникает при учёте процессов переброса. При низких температурах, когда длина свободного пробега больше размеров образца L, роль длины свободного пробега играет L. Кинетическое уравнение для фононов позволяет исследовать теплопроводность и поглощение звука в диэлектриках. Если длина свободного пробега для нормальных процессов значительно меньше длины свободного пробега для процессов переброса, то система фононов в кристалле при низких температуpax подобна обычному газу. Нормальные столкновения устанавливают внутреннее равновесие в каждом элементе объёма газа, которьй может двигаться со скоростью V, мало меняющейся на длине свободного пробега для нормальных столкновений. Поэтому можно построить уравнения гидродинамики фононного газа в диэлектрике.

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

Физическая кинетика магнитных явлений основана на решении кинетического уравнения для магнонов. Она позволяет вычислить динамическии восприимчивости магнитных систем в переменных полях, изучить кинетику процессов намагничивания.

Физическая кинетика явлений при прохождении быстрых частиц через вещество основана на решении системы кинетических уравнений для быстрых частиц и вторичных частиц, возникающих при столкновениях, например для γ \gamma -лучей (фотонов) с учётом различных процессов в среде (фотоэффекта, комптоновского рассеяния, образования пар). В этом случае кинетика позволяет вычислить коэффициенты поглощения и рассеяния быстрых частиц.

Фазовые переходы[править | править код]

Физическая кинетика фазовых переходов первого рода, то есть со скачком энтропии, связана с образованием и ростом зародышей новой фазы. Функция распределения зародышей по их размерам (если зародыши считать макроскопическими образованиями, а процесс роста — медленным) удовлетворяет уравнению Фоккера — Планка: f t = t ( D f α A f ) \frac{\partial f}{\partial t} = \frac{\partial}{\partial t} \left (D \frac{\partial f}{\partial \alpha} - A f \right )

где а — радиус зародыша, D — «коэффициент диффузии зародышей по размерам», А пропорционально минимальной работе, которую нужно затратить на создание зародыша данного размера. Кинетика фазовых переходов 2-го рода в наиболее простом приближении основана на уравнении релаксации параметра порядка η \eta , характеризующего степень упорядоченности, возникающей при фазовом переходе: η t = γ Ω η , \frac{\partial \eta}{\partial t} = - \gamma \frac{\partial \Omega}{\partial \eta},

где γ \gamma — постоянный коэффициент, Ω \Omega термодинамический потенциал в переменных T и μ \mu , вблизи точки фазового перехода зависящий от μ \mu . Для этой зависимости используется разложение по степеням μ 2 \mu^2 и T T c T-T_c , где T c T_c — температура фазового перехода.

Явления переноса в жидкостях[править | править код]

Теорию явлений переноса в жидкостях также можно отнести к физической кинетике, хотя для жидкостей метод кинетических уравнений непригоден, но для них возможен более общий подход, основанный также на иерархии времён релаксации. Для жидкости время установления равновесия в макроскопически малых (но содержащих ещё большое число молекул) элементарных объёмах значительно больше, чем время релаксации во всей системе, вследствие чего в малых элементах объёма приближённо устанавливается статистическое равновесие. Поэтому в качестве исходного приближения при решении уравнения Лиувилля можно принять локально равновесное распределение Гиббса с температурой T (x, t), химическим потенциалом μ \mu (x, t) и гидродинамической скоростью V(x,t), соответствующими рассматриваемой точке жидкости. Например, для однокомпонентпой жидкости локально равновесная функция распределения (или матрица плотноси) имеет вид

f = 1 Z e x p ( β ( x , t ) [ H ( x ) μ ( x , t ) n ( x ) ] d x ) f = \frac {1}{Z} exp \left ( - \int{\beta (x, t) [H'(x) - \mu (x, t) n(x)] dx } \right )

, где β ( x , t ) = 1 / k T ( x , t ) \beta (x,t) = 1/kT (x, t) , H ( x ) = H ( x ) p ( x ) B ( x , t ) + 1 2 m n ( x ) V 2 ( x , t ) H' (x) = H(x) - p(x) B(x,t) +\frac{1}{2} m n(x) V^2(x,t) — плотность энергии в системе координат, движущейся вместе с элементом жидкости, H(х) — плотность энергии в неподвижной системе координат, p(х) — плотность импульса, n(x) — плотность числа частиц, рассматриваемые как фазовые функции, то есть функции от координат и импульсов всех частиц, например n ( x ) = j N δ ( x x j ) n(x) = \sum _{j} ^{N} {\delta(x-x_j)} . Приближённое решение уравнения Лиувилля для состояний, близких к статистически равновесному, позволяет вывести уравнения теплопроводности и Навье—Стокса для жидкости и получить микроскопические выражения для кинетических коэффициентов теплопроводности и вязкости через пространственно-временные корреляционные функции плотностей потоков энергии и импульсов всех частиц системы. Этот же подход возможен и для смеси жидкостей. Подобное решение уравнения Лиувилля есть его частное решение, зависящее от времени лишь через параметры β(x,t), ν(x,t), V x,t), соответствующие сокращённому гидродипамическому описанию неравновесного состояния системы, которое справедливо, когда все гидродинамические параметры мало меняются на расстояниях порядка длины свободного пробега (для газов) или длины корреляций потоков энергии или импульса (для жидкостей).

К задачам физической кинетики относится также вычисление обобщённой восприимчивости, выражающей линейную реакцию физической системы на включение внешнего ноля. Её можно выразить через функции Грина с усреднением по состоянию, которое может быть и неравновесным.

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

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

  • Лифшиц Е. М., Питаевский Л. П., Физическая кинетика.
  • Гуревич Л. Э., Основы физической кинетики.
  • Балеску Р., Равновесная и неравновесная статистическая механика.
  • К. Кейз, П. Цвайфель. Линейная теория переноса.*Б. Дэвисон, Теория переноса нейтронов.
  • А.М. Кольчужкин, В.В. Учайкин. Введение в теорию прохождения частиц через вещество.
  • Боголюбов Н. Н., Проблемы динамической теории в статистической физике.
  • Климонтович Ю. Л., Кинетическая теория неидеального газа и неидеальной плазмы.
  • Ферцигер Дж., Капер Г., Математическая теория процессов переноса в газах.

См. также[править | править код]