Кинетическое уравнение больцмана. Уравнение больцмана

МОСКОВСКИЙ ЭНЕРГЕТИЧЕСКИЙ ИНСТИТУТ

(технический университет)

ФАКУЛЬТЕТ ЭЛЕКТРОННОЙ ТЕХНИКИ

РЕФЕРАТ ПО ТЕМЕ

К ИНЕТИЧЕСКОЕ УРАВНЕНИЕ Б ОЛЬЦМАНА.

ВЫПОЛНИЛ:

Коркин С.В.

ПРЕПОДАВАТЕЛЬ

Шеркунов Ю.Б.

Вторая половина работы набита достаточно сложной математикой . Автор ([email protected], [email protected])не считает этот курсовой идеальным, он может служить лишь отправной точкой для написания более совершенной (и понятной) работы. Текст не является копией книги. Вспомогательную литературу см. в конце.

Курсовой принят с отметкой ОТЛ. (Окончательный вариант работы немножко затерялся. Предлагаю воспользоваться предпоследней «версией»).

Введение……………………………………………………………………………… 3

Условные обозначения………………………………………………………………. 4

§1 Функция распределения.

§2 Столкновение частиц.

§3 Определение вида интеграла столкновений

и уравнения Больцмана.

§4. Кинетическое уравнение для слабо неоднородного газа.

Теплопроводность газа.

Некоторые условные обозначения:

n - концентрация частиц;

d - среднее расстояние между частицами;

V - некоторый объём системы;

P - вероятность некоторого события;

f - функция распределения;

Введение.

Разделы физики термодинамика, статистическая физика и физическая кинетика занимаются изучением физических процессов, происходящих в макроскопических системах - телах, состоящих из большого числа микрочастиц. В зависимости от вида системы такими микрочастицами могут являться атомы, молекулы, ионы, электроны, фотоны или иные частицы. На сегодняшний день существуют два основных метода исследования состояний макроскопических систем - термодинамический, характеризующий состояние системы через макроскопические легко измеряемые параметры (например, давление, объём, температура, количество молей или концентрация вещества) и, по сути, не учитывающий атомно-молекулярную структуру вещества, и статистический метод, основанный на атомно-молекулярной модели рассматриваемой системы. Термодинамический метод не будет затрагиваться в данной работе. По известным законам поведения частиц системы статистический метод позволяет установить законы поведения всей макросистемы в целом. С целью упрощения решаемой задачи при статистическом подходе делается ряд предположений (допущений) о поведении микрочастиц и, следовательно, результаты, полученные статметодом, справедливы лишь в пределах сделанных допущений. Статистический метод использует вероятностный подход к решению задач, для использования этого метода система обязана содержать достаточно большое количество частиц. Одна из задач, решаемая статметодом, - вывод уравнения состояния макроскопической системы. Состояние системы может быть неизменным во времени (равновесная система) либо может изменяться с течением времени (неравновесная система). Изучением неравновесных состояний систем и процессов, происходящих в таких системах, занимается физическая кинетика.

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

§1 Функция распределения.

Для вывода кинетического уравнения Больцмана рассмотрим одноатомный идеальный газ, т.е. достаточно разряженный газ, состоящий из электрически нейтральных атомов или молекул. Единственным видом взаимодействия между частицами идеального газа являются столкновения между молекулами, происходящие, однако, настолько редко, что каждая молекула почти всё время движется как свободная. Рассматривая частицы газа как классические, можно утверждать, что на одну частицу приходиться объём. Число частиц в единице объёма есть концентрация. Значит среднее расстояние между частицами есть (предполагается достаточно большим по сравнению с радиусом действия межмолекулярных сил d). При получении уравнения Больцмана сделаем следующие предположения:

Частицы газа неразличимы (одинаковы);

Частицы сталкиваются только попарно (пренебрегаем столкновением одновременно трех и более частиц);

Непосредственно перед столкновением частицы движутся по одной прямой навстречу друг другу;

Столкновение молекул есть прямой центральный упругий удар;

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

Функция распределения молекул газа даётся в их фазовом:-пространстве.

есть совокупность обобщённых координат всех молекул; - совокупность обобщённых импульсов молекул. Соответственно

и. Обозначим через

элемент объёма фазового пространства молекулы. В заданном элементе фазового пространства находиться (в среднем) число частиц, равное (т.е. рассматриваются молекулы, значения q и p которых лежат в выделенных интервалах dq и dp). Функция распределения молекул газа выше была определена в фазовом пространстве, тем не менее, она может быть выражена через иные переменные, отличные от обобщённых координат и импульсов частицы. Произведём выбор аргументов функции f.

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

Поступательное движение классической частицы описывается координатами

центра тяжести частица и вектором скорости или вектором импульса (, где m – масса частицы). Для одноатомного газа поступательное движение – единственный вид движения частиц; число степеней свободы равно трём. Если частица представляет собой многоатомную молекулу, то возникают дополнительные степени свободы, связанные с вращением молекулы в пространстве и колебанием атомов в молекуле. Условиями применения квантовой механики являются малые массы и высокие концентрации частиц, а так же низкие температуры. Не рассматривая область низких температур, будем считать вращательное движение молекул газа классическим. Любое классическое вращательное движение описывается, прежде всего, вращательным моментом сил, действующих на тело. Под действием момента двухатомная молекула приходит во вращение в плоскости, перпендикулярной вектору момента. Кроме того, положение молекулы характеризуется углом поворота оси молекулы в плоскости вращения.

Рассмотрим молекулу водорода (или любую другую двухатомную молекулу) при Т=300 К. Согласно закону равнораспределения на каждую степень свободы (поступательную, вращательную или колебательную) в среднем приходится одинаковая кинетическая энергия, равная.

Пусть I - момент инерции молекулы, m - масса, d - среднее расстояние между атомами в молекуле.


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

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

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

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

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

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

§2 Столкновение частиц.

Рассмотрим столкновение молекул, одни из которых обладают значениями величин Г, лежащими в заданном интервале, а другие – в интервале. В результате столкновения молекулы приобретают значения величин Г в интервалах соответственно и. Далее для краткости будем говорить о столкновении молекул и с переходом

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

(штрихом обозначены конечные состояния, без штриха - начальные). Вероятность столкновения обладает важным свойством, которое следует из законов механики, относительно обращения знака времени. Если обозначить верхним индексом Т значения всех величин, получившихся при обращении знака времени, то будет иметь место равенство

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

Произведение дифференциалов представляет собой элемент фазового пространства, который не изменяется при обращении времени (дифференциалы в обеих сторонах равенства можно опустить) . Не изменяется так же потенциальная энергия молекул, и, следовательно, равновесная (больцмановская) функция распределения, которая зависит только от енергии:


(2)

V – макроскопическая скорость движения газа как целого. В силу закона сохранения энергии при столкновении двух молекул. Поэтому можно записать (3)

Отметим ещё тот факт, что сама функция вероятности в принципе может быть определена лишь путём решения механической задачи о столкновении частиц. Написанное выше равенства (1) , (2) и (3) дадут после сокращений в (1)

С учётом утверждения (*)

Интегрируя последнее равенство (для использования в дальнейшем) получаем соотношение:

§3 Вывод кинетического уравнения.

Рассмотрим производную от функции распределения по времени:

При движении молекул газа в отсутствии внешнего поля величины Г, как интегралы движения, не изменяются.


Выражение для производной примет вид: (6)

Пусть теперь газ находится во внешнем потенциальном поле, действующем на координаты центра тяжести молекул (например, в гравитационном поле). И пусть F – сила, действующая со стороны поля на частицу.


Правую часть равенства (6) обозначим через. Символ означает

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

есть отнесённое к единице времени изменение за счёт столкновений числа молекул в фазовом объёме. Полное изменение функции распределения в заданной точке фазового пространства запишется в виде:

(8)


Величина называется интегралом столкновений, а уравнение вида (8) – кинетическим уравнением. Реальный смысл кинетическое уравнение (8) примет только после определения вида интеграла столкновений.

§3 Определение вида интеграла столкновений и уравнения Больцмана.

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

Полное число столкновений при вышеуказанном переходе со всеми возможными значениями

При заданном, происходящих в единицу времени в объёме,определяется интегралом

В то же время происходят столкновения иного рода (называемые “приходом”), в результате которых молекулы, обладавшие до столкновения значениями величин, лежащими вне заданного интервала, попадают в этот интервал. Такие переходы могут быть обозначены следующим образом: (со всеми возможными значениями при заданном). Аналогично первому типу перехода полное число таких столкновений в единицу времени в объёме равно:


В результате всех столкновений изменение числа молекул в единицу времени в элементарном объёме определяется разностью между числом актов ухода и числом актов прихода:

(9) , где

Интеграл столкновений может быть определён как:

(изменение числа частиц в единицу времени в фазовом объёме dVdГ)

Из соотношений (8) и (9) получим вид интеграла столкновений


Заметим, что во втором члене подынтегрального выражения интегрирование по имеет

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

и кинетического уравнения

Полученное интегрально - дифференциальное уравнение носит название уравнения Больцмана.

Рассмотрим не зависящее от времени распределение в состоянии равновесия системы в отсутствии внешних воздействий. Такое распределение является стационарным (не зависит от времени) и однородным (не изменяется в области пространства, занимаемой системой). Наложенные условия обнуляют производную функции распределения по времени и трём координатам; левая часть кинетического уравнения обращается в нуль. Подынтегральное выражение обращается в нуль вследствие равенства (3). Следовательно, равновесное распределение в отсутствии внешних полей удовлетворяет кинетическому уравнению тождественным образом. Если газ находится в равновесном состоянии под действием внешнего потенциального (например, гравитационного) поля, то функция распределения и в этом случае удовлетворяет кинетическому уравнению. Действительно, равновесное распределение выражается через интеграл движения – полную энергию молекулы. Левая часть кинетического уравнения представляет собой полную производную, которая равна нулю как производная от функции, зависящей только от интегралов движения. Правая часть уравнения, как уже было указано, есть нуль. Таким образом, кинетическому уравнению удовлетворяет и функция распределения газа, находящегося в равновесии во внешнем потенциальном поле.

К указанным во “Введении” допущениям добавим ещё одно: столкновения молекул рассматриваются как мгновенные акты, происходящие в одной “точке” пространства. Кинетическое уравнение описывает процес, который протекает в интервале времени, много большем по сравнению с длительностью столкновений. В то же время, рассматриваемая область системы должна значительно превышать область столкновения частиц, которая имеет размеры порядка величины радиуса действия молекулярных сил d. Время столкновения по порядку величины может быть определено как (- средняя скорость движения молекул в газе). Полученные значения представляют собой нижний предел расстояния и времени, при рассмотрении которых допускается применение кинетического уравнения. Реальные физические задачи не требуют столь детального описания процесса; размеры системы и время наблюдения значительно превышают требуемый минимум.

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

- среднее расстояние между молекулами;


Величина - время свободного пробега. Для грубой оценки интеграла столкновений можно использовать:

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

§3 Переход к макроскопическим уравнениям. Гидродинамическое уравнение непрерывности.

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

В качестве примера рассмотрим рассмотрим способ получения гидродинамического уравнения.

Выражение определяет плотность распределения молекул газа в пространстве (концентрацию молекул газа). Произведение массы одной молекулы (предполагается, что газ состоит из одинаковых частиц) на плотность распределения молекул даёт массовую плотность газа: . Обозначим через макроскопическую скорость движения газа как целого, а через микроскопическую скорость молекул. Макроскопическая скорость (скорость движения центра масс) может быть определена как средняя величина от микроскопических скоростей молекул


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

Убедимся в справедливости этого равенства следующим способом:

Интегрирование производится по каждой из переменых, а значит можно, не меняя интеграла, произвести переобозначение переменных, например, во втором интеграле:


Последнее выражение, очевидно, равно нулю и, следовательно, справедливым является равенство (14).

Запишем кинетическое уравнение и, предварительно умножив обе его части на массу частицы m , интегрируем его по:

Отсюда немедленно получаем гидродинамическое уравнение непрерывности:

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

§4. Слабо неоднородный газ. Теплопроводность газа.

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

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

Преобразуем интеграл столкновений в кинетическом уравнении (13): подстановка выражений функции распределения и поправки, обнуление интегралов столкновений,содержащих равновесную функцию распределения, сокращение членов, не содержащих малой поправки. Члены первого порядка дадут. Символ введен для обозначения линейного интегрального оператора


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

*************************************************

§4. Вычисление коэффициента теплопроводности одноатомного газа

Для вычисления коэффициента теплопроводности газа необходимо решать записанное выше уравнение с градиентом температуры.


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


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

И). Третий интеграл представляет из себя дополнительное условие на функцию g. Если кинетическое уравнение решено и функция

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

Если система находится в рановесии, то и этот интеграл равен нулю за счёт интегрирования по всем возможным направлениям в газе. При подстановке в () остаётся

В компонентах

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


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

Поток Q должен быть направлен в сторону, противоположную градиенту температуры, а величина соответственно должна быть положительна, что автоматически обеспечивается кинетическим уравнением (). В одноатомных газах скорость v- единственный вектор от которого зависит функция g (в многоатомных газах имеет место зависимость g не только от скорости v , но и отмомента M). Для одноатомных газов функция g имеет вид:

.

§5.Пример решения кинетического уравнения

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


а равновесная функция распределения примет вид.

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


В этой формуле r – произвольное, а s – целое положительное число либо нуль. В честноти

Свойство ортогональности этих полиномов при заданном индексе r и различных индексах s выаглядит следующим образом

Решение уравнения ищем в виде следующего разложения

Опустив в разложении член с s=0 , получим выражение адовлетворяющее () (нтеграл обнуляется в силу ортогональности полиномов с различными s). Выражение в скобках в левой стороне ()

есть. Уравнение () принимает вид


Для последнего выражения введены обозначения


Уравнение с l=0 отсутствует, поскольку в силу сохранения импульса

Коэффициент теплопроводности вычисляется подстановкой выражения () в интеграл (). С учётом условия () интеграл (с) может быть представлен в виде

В результате находим.

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

Заключение.

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

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

Список литературы.

1. Е.М.Лифшиц, Л.П.Питаевский. Физическая кинетика. Наука, М., 1979 г.

2. Ю.Б.Румер, М.Ш.Рывкин. Термодинамика, статистическая физика и кинетика.

КИНЕТИЧЕСКОЕ УРАВНЕНИЕ БОЛЬЦМАНА - интегродифференц. ур-ние, к-рому удовлетворяют неравновесные одночастичные функции распределения системы из большого числа частиц, напр, ф-ция распределения молекул газа по скоростям и координатам r , ф-ции распределения электронов в металле, в кристалле и т. п. К. у. Б.- осн. ур-ние мик-роскопич. теории неравновесных процессов (кинетики физической) , в частности кинетической теории газов . К. у. Б. в узком смысле наз. выведенное Л. Больцма-ном (L. Boltzmann) кинетич. ур-ние для газов малой , молекулы к-рых подчиняются классич. механике. К. у. Б. для квазичастиц в кристаллах, напр. для электронов в металле, наз. также кинетич. ур-ниями или ур-ниями переноса.

К. у. Б. представляет собой ур-ние баланса числа частиц (точнее, точек, изображающих состояние частиц) в элементе фазового объёма ; dr= =dxdydz )и выражает тот факт, что изменение ф-ции распределения частиц со временем t происходит вследствие движения частиц под действием внеш. сил и столкновений между ними. Для газа, состоящего из частиц одного сорта, К. у. Б. имеет вид

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


Скорости частиц до столкновения, - скорости тех же частиц после столкновения, - величина относит. скорости сталкивающихся частиц, - дифференц. эфф. сечение рассеяния частиц в телесный угол в лаб. системе координат, - угол между относит. скоростью и линией центров. Напр., для жёстких упругих сфер, имеющих радиус R , = , для частиц, взаимодействующих по закону центр. сил, (b - прицельный параметр, - азимутальный угол линии центров).

К. у. Б. учитывает только парные столкновения между молекулами; оно справедливо при условии, что длина свободного пробега молекул значительно больше линейных размеров области, в к-рой происходит столкновение (для газа из упругих частиц это область порядка диаметра частиц). Поэтому К. у. Б. применимо для не слишком плотных газов. Иначе будет несправедливо осн. предположение об отсутствии корреляции между состояниями сталкивающихся частиц (гипотеза молекулярного хаоса). Если система находится в статистич. равновесии, то интеграл столкновений (2) обращается в нуль и решением К. у. Б. является Максвелла распределение .

При более строгом подходе для построения К. у. Б. исходят из Лиувилля уравнения для плотности распределения всех молекул газа в фазовом пространстве, из к-рого получают систему ур-ний для ф-ций распределения одной, двух и т д. молекул (Боголюбова уравнения) . Эту цепочку ур-ний решают с помощью разложения по степеням плотности частиц с использованием граничного условия ослабления корреляций, заменяющего гипотезу молекулярного хаоса.

Решение К. у. Б. при разл. предположениях о силах взаимодействия между частицами - предмет кинетич. теории газов, к-рая позволяет вычислить кинетические коэффициенты и получить макроскопич. ур-ния для процессов переноса (вязкости, диффузии, теплопроводности) .

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

где знак минус соответствует Ферми - Дирака статистике , а знак плюс - Бозе - Эйнштейна статистике, g - статистич. вес состояния (g = l для частиц со спином, равным нулю, и g=2 для частиц со спином),- импульс частицы. Ф-ции нормированы так, что представляют ср. число частиц в точке . Равновесные ф-ции распределения Ферми и Бозе обращают в нуль оператор столкновения (3).

Важным частным случаем К. у. Б. является кинетич. ур-ние для нейтронов, к-рые рассеиваются и замедляются ядрами среды. В этом случае внеш. сил нет и в ур-нии (1) надо положить F=0 . Плотность числа нейтронов обычно мала, так что можно пренебречь столкновениями между ними и учитывать лишь их столкновения с ядрами среды (см. Диффузия нейтронов, Замедление нейтронов) .

Процессы переноса, связанные с движением электронов в металле, также можно исследовать с помощью К. у. Б. В отсутствие решётки электроны свободно распространяются в металле н описываются , модулированными с периодом решётки и зависящими от k ; и номера энергетич. зоны l . Тепловое движение атомов решётки нарушает периодичность и приводит к рассеянию электронов (столкновениям между электронами и фононами). Ф-ция распределения электронов n(k, l, t )удовлетворяет К. у. Б. типа (1), в к-ром F = (E и Н - напряжённости электрич. и магн. полей, е - электрона), а интеграл столкновений имеет вид

где n=n(k ,l) , - волновые векторы и номера зон до и после столкновения, N= =N (f , s) - ф-ция распределения фононов, f и s - волновой вектор и поляризация фононов, - нач. и конечная энергии электрона при возбуждении фонона с энергией - дельта-ф-ция, - матричные элементы перехода электрона из состояния k , l в состояние , к-рые оценивают, исходя из определ. гипотез о механизме взаимодействия электронов с решёткой. Выражение (4) получено в предположении, что время свободного пробега электронов значительно больше неопределённости для времени столкновения. Теория электропроводности, термоэлектрич. и гальвано-магн. явлений в металлах и полупроводниках основана на решении К. у. Б.

В нек-рых случаях конденсиров. систем, когда известен характер теплового движения, можно построить К. у. Б. для элементарных возбуждений (квазичастиц). Напр., теория процессов переноса энергии в кристал-лич. решётке основана на ур-нии такого типа. Если в выражении для потенц. энергии решётки ограничиться квадратичными относительно смещений атомов членами, то тепловое движение атомов в кристалле описывается свободно распространяющимися фононами - квантами нормальных колебаний решётки. Учёт членов 3-й степени приводит к возможности столкновений между фононами. В результате ф-ция распределения фононов N (f , s ) будет изменяться во времени согласно кинетич. ур-нию

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

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

На решении этих ур-ний основана каскадная теория ливней.

Лит . см. при статьях Кинетическая теория газов . Кинетика физическая. Д. Я. Зубарев .

Материал из Википедии - свободной энциклопедии

Уравне́ние Бо́льцмана (кинети́ческое уравнение Больцмана ) - уравнение, названное по имени Людвига Больцмана , который его впервые рассмотрел, и описывающее статистическое распределение частиц в газе или жидкости. Является одним из самых важных уравнений физической кинетики (области статистической физики , которая описывает системы, далёкие от термодинамического равновесия, например, в присутствии градиентов температур и электрического поля). Уравнение Больцмана используется для изучения переноса тепла и электрического заряда в жидкостях и газах , и из него выводятся транспортные свойства, такие как электропроводность , эффект Холла , вязкость и теплопроводность . Уравнение применимо для разреженных систем, где время взаимодействия между частицами мало (гипотеза молекулярного хаоса).

Формулировка

Уравнение Больцмана описывает эволюцию во времени (t ) функции распределения плотности f (x , p , t ) в одночастичном фазовом пространстве , где x и p - координата и импульс соответственно. Распределение определяется так, что

f(\mathbf{x},\mathbf{p},t)\,d^3x\,d^3p

пропорционально числу частиц в фазовом объёме d³x d³p в момент времени t . Уравнение Больцмана

\frac{\partial f}{\partial t} + \frac{\partial f}{\partial \mathbf{x}} \cdot \frac{\mathbf{p}}{m} + \frac{\partial f}{\partial \mathbf{p}} \cdot \mathbf{F} = \left. \frac{d f}{d t} \right|_{\mathrm{coll}}. Здесь F (x , t ) - поле сил, действующее на частицы в жидкости или газе, а m - масса частиц. Слагаемое в правой части уравнения добавлено для учёта столкновений между частицами и называется интегралом столкновений . Если оно равно нулю, то частицы не сталкиваются вовсе. Этот случай часто называют одночастичным уравнением Лиувилля . Если поле сил F (x , t ) заменить подходящим самосогласованным полем, зависящим от функции распределения f, то получим уравнение Власова , описывающее динамику заряженных частиц плазмы в самосогласованном поле. Классическое же уравнение Больцмана используется в физике плазмы , а также в физике полупроводников и металлов (для описания кинетических явлений, то есть переноса заряда или тепла, в электронной жидкости).

\hat{\mathbf{L}}_\mathrm{GR}=\sum_\alpha p^\alpha\frac{\partial}{\partial x^\alpha}-\sum_{\alpha\beta\gamma}\Gamma^{\alpha}{}_{\beta\gamma}p^\beta p^\gamma\frac{\partial}{\partial p^\alpha},

Интеграл столкновений

Столкновения между частицами приводит к изменению их скоростей. Если W(\mathbf{v},\mathbf{v}^\prime)d^3v^\prime dt задает вероятность рассеяния частицы из состояния со скоростью \mathbf{v} в состояние со скоростью \mathbf{v}^\prime, то интеграл столкновений для классических частиц записывается в виде

\left.\frac{\partial f}{\partial t}\right|_{coll}=\int_{\mathbf{v}^\prime} d^3v^\prime .

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

Приближение времени релаксации

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

f = f_0 + f_1,

где f_0(\mathbf{v}) - равновесная функция распределения, зависит только от скоростей частиц и известная из термодинамики, а f_1 - небольшое отклонение.

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

- \frac{f_1}{\tau} = - \frac{f-f_0}{\tau} ,

См. также

Напишите отзыв о статье "Кинетическое уравнение Больцмана"

Примечания

Ссылки

Литература

  • Черчиньяни К. Теория и приложения уравнения Больцмана. - М .: Мир, 1978. - 495 с.

Отрывок, характеризующий Кинетическое уравнение Больцмана

Движение человечества, вытекая из бесчисленного количества людских произволов, совершается непрерывно.
Постижение законов этого движения есть цель истории. Но для того, чтобы постигнуть законы непрерывного движения суммы всех произволов людей, ум человеческий допускает произвольные, прерывные единицы. Первый прием истории состоит в том, чтобы, взяв произвольный ряд непрерывных событий, рассматривать его отдельно от других, тогда как нет и не может быть начала никакого события, а всегда одно событие непрерывно вытекает из другого. Второй прием состоит в том, чтобы рассматривать действие одного человека, царя, полководца, как сумму произволов людей, тогда как сумма произволов людских никогда не выражается в деятельности одного исторического лица.
Историческая наука в движении своем постоянно принимает все меньшие и меньшие единицы для рассмотрения и этим путем стремится приблизиться к истине. Но как ни мелки единицы, которые принимает история, мы чувствуем, что допущение единицы, отделенной от другой, допущение начала какого нибудь явления и допущение того, что произволы всех людей выражаются в действиях одного исторического лица, ложны сами в себе.
Всякий вывод истории, без малейшего усилия со стороны критики, распадается, как прах, ничего не оставляя за собой, только вследствие того, что критика избирает за предмет наблюдения большую или меньшую прерывную единицу; на что она всегда имеет право, так как взятая историческая единица всегда произвольна.
Только допустив бесконечно малую единицу для наблюдения – дифференциал истории, то есть однородные влечения людей, и достигнув искусства интегрировать (брать суммы этих бесконечно малых), мы можем надеяться на постигновение законов истории.
Первые пятнадцать лет XIX столетия в Европе представляют необыкновенное движение миллионов людей. Люди оставляют свои обычные занятия, стремятся с одной стороны Европы в другую, грабят, убивают один другого, торжествуют и отчаиваются, и весь ход жизни на несколько лет изменяется и представляет усиленное движение, которое сначала идет возрастая, потом ослабевая. Какая причина этого движения или по каким законам происходило оно? – спрашивает ум человеческий.
Историки, отвечая на этот вопрос, излагают нам деяния и речи нескольких десятков людей в одном из зданий города Парижа, называя эти деяния и речи словом революция; потом дают подробную биографию Наполеона и некоторых сочувственных и враждебных ему лиц, рассказывают о влиянии одних из этих лиц на другие и говорят: вот отчего произошло это движение, и вот законы его.
Но ум человеческий не только отказывается верить в это объяснение, но прямо говорит, что прием объяснения не верен, потому что при этом объяснении слабейшее явление принимается за причину сильнейшего. Сумма людских произволов сделала и революцию и Наполеона, и только сумма этих произволов терпела их и уничтожила.
«Но всякий раз, когда были завоевания, были завоеватели; всякий раз, когда делались перевороты в государстве, были великие люди», – говорит история. Действительно, всякий раз, когда являлись завоеватели, были и войны, отвечает ум человеческий, но это не доказывает, чтобы завоеватели были причинами войн и чтобы возможно было найти законы войны в личной деятельности одного человека. Всякий раз, когда я, глядя на свои часы, вижу, что стрелка подошла к десяти, я слышу, что в соседней церкви начинается благовест, но из того, что всякий раз, что стрелка приходит на десять часов тогда, как начинается благовест, я не имею права заключить, что положение стрелки есть причина движения колоколов.
Всякий раз, как я вижу движение паровоза, я слышу звук свиста, вижу открытие клапана и движение колес; но из этого я не имею права заключить, что свист и движение колес суть причины движения паровоза.
Крестьяне говорят, что поздней весной дует холодный ветер, потому что почка дуба развертывается, и действительно, всякую весну дует холодный ветер, когда развертывается дуб. Но хотя причина дующего при развертыванье дуба холодного ветра мне неизвестна, я не могу согласиться с крестьянами в том, что причина холодного ветра есть раэвертыванье почки дуба, потому только, что сила ветра находится вне влияний почки. Я вижу только совпадение тех условий, которые бывают во всяком жизненном явлении, и вижу, что, сколько бы и как бы подробно я ни наблюдал стрелку часов, клапан и колеса паровоза и почку дуба, я не узнаю причину благовеста, движения паровоза и весеннего ветра. Для этого я должен изменить совершенно свою точку наблюдения и изучать законы движения пара, колокола и ветра. То же должна сделать история. И попытки этого уже были сделаны.
Для изучения законов истории мы должны изменить совершенно предмет наблюдения, оставить в покое царей, министров и генералов, а изучать однородные, бесконечно малые элементы, которые руководят массами. Никто не может сказать, насколько дано человеку достигнуть этим путем понимания законов истории; но очевидно, что на этом пути только лежит возможность уловления исторических законов и что на этом пути не положено еще умом человеческим одной миллионной доли тех усилий, которые положены историками на описание деяний различных царей, полководцев и министров и на изложение своих соображений по случаю этих деяний.

Силы двунадесяти языков Европы ворвались в Россию. Русское войско и население отступают, избегая столкновения, до Смоленска и от Смоленска до Бородина. Французское войско с постоянно увеличивающеюся силой стремительности несется к Москве, к цели своего движения. Сила стремительности его, приближаясь к цели, увеличивается подобно увеличению быстроты падающего тела по мере приближения его к земле. Назади тысяча верст голодной, враждебной страны; впереди десятки верст, отделяющие от цели. Это чувствует всякий солдат наполеоновской армии, и нашествие надвигается само собой, по одной силе стремительности.
В русском войске по мере отступления все более и более разгорается дух озлобления против врага: отступая назад, оно сосредоточивается и нарастает. Под Бородиным происходит столкновение. Ни то, ни другое войско не распадаются, но русское войско непосредственно после столкновения отступает так же необходимо, как необходимо откатывается шар, столкнувшись с другим, с большей стремительностью несущимся на него шаром; и так же необходимо (хотя и потерявший всю свою силу в столкновении) стремительно разбежавшийся шар нашествия прокатывается еще некоторое пространство.

МОСКОВСКИЙ ЭНЕРГЕТИЧЕСКИЙ ИНСТИТУТ

(технический университет)

ФАКУЛЬТЕТ ЭЛЕКТРОННОЙ ТЕХНИКИ

РЕФЕРАТ ПО ТЕМЕ

КИНЕТИЧЕСКОЕ УРАВНЕНИЕ БОЛЬЦМАНА.

ВЫПОЛНИЛ:

Коркин С.В.

ПРЕПОДАВАТЕЛЬ

Шеркунов Ю.Б.

Вторая половина работы набита достаточно сложной математикой. Автор ([email protected], [email protected])не считает этот курсовой идеальным, он может служить лишь отправной точкой для написания более совершенной (и понятной) работы. Текст не является копией книги. Вспомогательную литературу см. в конце.

Курсовой принят с отметкой ОТЛ. (Окончательный вариант работы немножко затерялся. Предлагаю воспользоваться предпоследней «версией»).

Введение……………………………………………………………………………… 3

Условные обозначения………………………………………………………………. 4

§1 Функция распределения.

§2 Столкновение частиц.

§3 Определение вида интеграла столкновений

и уравнения Больцмана.

§4. Кинетическое уравнение для слабо неоднородного газа.

Теплопроводность газа.

Некоторые условные обозначения:

n - концентрация частиц;

d - среднее расстояние между частицами;

V - некоторый объём системы;

P - вероятность некоторого события;

f - функция распределения;

Введение.

Разделы физики термодинамика, статистическая физика и физическая кинетика занимаются изучением физических процессов, происходящих в макроскопических системах - телах, состоящих из большого числа микрочастиц. В зависимости от вида системы такими микрочастицами могут являться атомы, молекулы, ионы, электроны, фотоны или иные частицы. На сегодняшний день существуют два основных метода исследования состояний макроскопических систем - термодинамический, характеризующий состояние системы через макроскопические легко измеряемые параметры (например, давление, объём, температура, количество молей или концентрация вещества) и, по сути, не учитывающий атомно-молекулярную структуру вещества, и статистический метод, основанный на атомно-молекулярной модели рассматриваемой системы. Термодинамический метод не будет затрагиваться в данной работе. По известным законам поведения частиц системы статистический метод позволяет установить законы поведения всей макросистемы в целом. С целью упрощения решаемой задачи при статистическом подходе делается ряд предположений (допущений) о поведении микрочастиц и, следовательно, результаты, полученные статметодом, справедливы лишь в пределах сделанных допущений. Статистический метод использует вероятностный подход к решению задач, для использования этого метода система обязана содержать достаточно большое количество частиц. Одна из задач, решаемая статметодом, - вывод уравнения состояния макроскопической системы. Состояние системы может быть неизменным во времени (равновесная система) либо может изменяться с течением времени (неравновесная система). Изучением неравновесных состояний систем и процессов, происходящих в таких системах, занимается физическая кинетика.

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

§1 Функция распределения.

Для вывода кинетического уравнения Больцмана рассмотрим одноатомный идеальный газ, т.е. достаточно разряженный газ, состоящий из электрически нейтральных атомов или молекул. Единственным видом взаимодействия между частицами идеального газа являются столкновения между молекулами, происходящие, однако, настолько редко, что каждая молекула почти всё время движется как свободная. Рассматривая частицы газа как классические, можно утверждать, что на одну частицу приходиться объём. Число частиц в единице объёма есть концентрация. Значит среднее расстояние между частицами есть (предполагается достаточно большим по сравнению с радиусом действия межмолекулярных сил d). При получении уравнения Больцмана сделаем следующие предположения:

Частицы газа неразличимы (одинаковы);

Частицы сталкиваются только попарно (пренебрегаем столкновением одновременно трех и более частиц);

Непосредственно перед столкновением частицы движутся по одной прямой навстречу друг другу;

Столкновение молекул есть прямой центральный упругий удар;

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

Функция распределения молекул газа даётся в их фазовом:-пространстве.

есть совокупность обобщённых координат всех молекул; - совокупность обобщённых импульсов молекул. Соответственно

и. Обозначим через

элемент объёма фазового пространства молекулы. В заданном элементе фазового пространства находиться (в среднем) число частиц, равное (т.е. рассматриваются молекулы, значения q и p которых лежат в выделенных интервалах dq и dp). Функция распределения молекул газа выше была определена в фазовом пространстве, тем не менее, она может быть выражена через иные переменные, отличные от обобщённых координат и импульсов частицы. Произведём выбор аргументов функции f.

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

Поступательное движение классической частицы описывается координатами

центра тяжести частица и вектором скорости или вектором импульса (, где m – масса частицы). Для одноатомного газа поступательное движение – единственный вид движения частиц; число степеней свободы равно трём. Если частица представляет собой многоатомную молекулу, то возникают дополнительные степени свободы, связанные с вращением молекулы в пространстве и колебанием атомов в молекуле. Условиями применения квантовой механики являются малые массы и высокие концентрации частиц, а так же низкие температуры. Не рассматривая область низких температур, будем считать вращательное движение молекул газа классическим. Любое классическое вращательное движение описывается, прежде всего, вращательным моментом сил, действующих на тело. Под действием момента двухатомная молекула приходит во вращение в плоскости, перпендикулярной вектору момента. Кроме того, положение молекулы характеризуется углом поворота оси молекулы в плоскости вращения.

Рассмотрим молекулу водорода (или любую другую двухатомную молекулу) при Т=300 К. Согласно закону равнораспределения на каждую степень свободы (поступательную, вращательную или колебательную) в среднем приходится одинаковая кинетическая энергия, равная.

Пусть I - момент инерции молекулы, m - масса, d - среднее расстояние между атомами в молекуле.


Средняя кинетическая энергия вращения молекулы;


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

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

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

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

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

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

§2 Столкновение частиц.

Рассмотрим столкновение молекул, одни из которых обладают значениями величин Г, лежащими в заданном интервале, а другие – в интервале. В результате столкновения молекулы приобретают значения величин Г в интервалах соответственно и. Далее для краткости будем говорить о столкновении молекул и с переходом

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

(штрихом обозначены конечные состояния, без штриха - начальные). Вероятность столкновения обладает важным свойством, которое следует из законов механики, относительно обращения знака времени. Если обозначить верхним индексом Т значения всех величин, получившихся при обращении знака времени, то будет иметь место равенство

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

Произведение дифференциалов представляет собой элемент фазового пространства, который не изменяется при обращении времени (дифференциалы в обеих сторонах равенства можно опустить) . Не изменяется так же потенциальная энергия молекул, и, следовательно, равновесная (больцмановская) функция распределения, которая зависит только от енергии:


(2)

V – макроскопическая скорость движения газа как целого. В силу закона сохранения энергии при столкновении двух молекул. Поэтому можно записать (3)

Отметим ещё тот факт, что сама функция вероятности в принципе может быть определена лишь путём решения механической задачи о столкновении частиц. Написанное выше равенства (1) , (2) и (3) дадут после сокращений в (1)

С учётом утверждения (*)

Интегрируя последнее равенство (для использования в дальнейшем) получаем соотношение:

§3 Вывод кинетического уравнения.

Рассмотрим производную от функции распределения по времени:

При движении молекул газа в отсутствии внешнего поля величины Г, как интегралы движения, не изменяются.


(последнее слагаемое в выражении производной обнуляется, т.к.)

(оператор набла)


Выражение для производной примет вид: (6)

Пусть теперь газ находится во внешнем потенциальном поле, действующем на координаты центра тяжести молекул (например, в гравитационном поле). И пусть F – сила, действующая со стороны поля на частицу.


Правую часть равенства (6) обозначим через. Символ означает

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

есть отнесённое к единице времени изменение за счёт столкновений числа молекул в фазовом объёме. Полное изменение функции распределения в заданной точке фазового пространства запишется в виде:

(8)


Величина называется интегралом столкновений, а уравнение вида (8) – кинетическим уравнением. Реальный смысл кинетическое уравнение (8) примет только после определения вида интеграла столкновений.

§3 Определение вида интеграла столкновений и уравнения Больцмана.

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

Полное число столкновений при вышеуказанном переходе со всеми возможными значениями

При заданном, происходящих в единицу времени в объёме,определяется интегралом

В то же время происходят столкновения иного рода (называемые “приходом”), в результате которых молекулы, обладавшие до столкновения значениями величин, лежащими вне заданного интервала, попадают в этот интервал. Такие переходы могут быть обозначены следующим образом: (со всеми возможными значениями при заданном). Аналогично первому типу перехода полное число таких столкновений в единицу времени в объёме равно:


В результате всех столкновений изменение числа молекул в единицу времени в элементарном объёме определяется разностью между числом актов ухода и числом актов прихода:

(9) , где

Интеграл столкновений может быть определён как:

(изменение числа частиц в единицу времени в фазовом объёме dVdГ)

Из соотношений (8) и (9) получим вид интеграла столкновений


Заметим, что во втором члене подынтегрального выражения интегрирование по имеет

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

и кинетического уравнения

Полученное интегрально - дифференциальное уравнение носит название уравнения Больцмана.

Рассмотрим не зависящее от времени распределение в состоянии равновесия системы в отсутствии внешних воздействий. Такое распределение является стационарным (не зависит от времени) и однородным (не изменяется в области пространства, занимаемой системой). Наложенные условия обнуляют производную функции распределения по времени и трём координатам; левая часть кинетического уравнения обращается в нуль. Подынтегральное выражение обращается в нуль вследствие равенства (3). Следовательно, равновесное распределение в отсутствии внешних полей удовлетворяет кинетическому уравнению тождественным образом. Если газ находится в равновесном состоянии под действием внешнего потенциального (например, гравитационного) поля, то функция распределения и в этом случае удовлетворяет кинетическому уравнению. Действительно, равновесное распределение выражается через интеграл движения – полную энергию молекулы. Левая часть кинетического уравнения представляет собой полную производную, которая равна нулю как производная от функции, зависящей только от интегралов движения. Правая часть уравнения, как уже было указано, есть нуль. Таким образом, кинетическому уравнению удовлетворяет и функция распределения газа, находящегося в равновесии во внешнем потенциальном поле.

К указанным во “Введении” допущениям добавим ещё одно: столкновения молекул рассматриваются как мгновенные акты, происходящие в одной “точке” пространства. Кинетическое уравнение описывает процес, который протекает в интервале времени, много большем по сравнению с длительностью столкновений. В то же время, рассматриваемая область системы должна значительно превышать область столкновения частиц, которая имеет размеры порядка величины радиуса действия молекулярных сил d. Время столкновения по порядку величины может быть определено как (- средняя скорость движения молекул в газе). Полученные значения представляют собой нижний предел расстояния и времени, при рассмотрении которых допускается применение кинетического уравнения. Реальные физические задачи не требуют столь детального описания процесса; размеры системы и время наблюдения значительно превышают требуемый минимум.

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

- среднее расстояние между молекулами;


Величина - время свободного пробега. Для грубой оценки интеграла столкновений можно использовать:

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

§3 Переход к макроскопическим уравнениям. Гидродинамическое уравнение непрерывности.

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

В качестве примера рассмотрим рассмотрим способ получения гидродинамического уравнения.

Выражение определяет плотность распределения молекул газа в пространстве (концентрацию молекул газа). Произведение массы одной молекулы (предполагается, что газ состоит из одинаковых частиц) на плотность распределения молекул даёт массовую плотность газа: . Обозначим через макроскопическую скорость движения газа как целого, а через микроскопическую скорость молекул. Макроскопическая скорость (скорость движения центра масс) может быть определена как средняя величина от микроскопических скоростей молекул


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

Убедимся в справедливости этого равенства следующим способом:

Интегрирование производится по каждой из переменых, а значит можно, не меняя интеграла, произвести переобозначение переменных, например, во втором интеграле:


Последнее выражение, очевидно, равно нулю и, следовательно, справедливым является равенство (14).

Запишем кинетическое уравнение и, предварительно умножив обе его части на массу частицы m , интегрируем его по:

Отсюда немедленно получаем гидродинамическое уравнение непрерывности:

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

§4. Слабо неоднородный газ. Теплопроводность газа.

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

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

Преобразуем интеграл столкновений в кинетическом уравнении (13): подстановка выражений функции распределения и поправки, обнуление интегралов столкновений,содержащих равновесную функцию распределения, сокращение членов, не содержащих малой поправки. Члены первого порядка дадут. Символ введен для обозначения линейного интегрального оператора


Указанный интеграл обращается в нуль для функций вида


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

*************************************************

§4. Вычисление коэффициента теплопроводности одноатомного газа

Для вычисления коэффициента теплопроводности газа необходимо решать записанное выше уравнение с градиентом температуры.


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


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

И). Третий интеграл представляет из себя дополнительное условие на функцию g. Если кинетическое уравнение решено и функция

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

Если система находится в рановесии, то и этот интеграл равен нулю за счёт интегрирования по всем возможным направлениям в газе. При подстановке в () остаётся

В компонентах

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


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

Поток Q должен быть направлен в сторону, противоположную градиенту температуры, а величина соответственно должна быть положительна, что автоматически обеспечивается кинетическим уравнением (). В одноатомных газах скорость v- единственный вектор от которого зависит функция g (в многоатомных газах имеет место зависимость g не только от скорости v , но и отмомента M). Для одноатомных газов функция g имеет вид:

.

§5.Пример решения кинетического уравнения

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


Для одноатомного газа теплоёмкость. Положив уравнению () придадим вид

Линейный интегральный оператор, соответствующий интегралу столкновений () ,определяется формулой


а равновесная функция распределения примет вид.

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


В этой формуле r – произвольное, а s – целое положительное число либо нуль. В честноти

Свойство ортогональности этих полиномов при заданном индексе r и различных индексах s выаглядит следующим образом

Решение уравнения ищем в виде следующего разложения

Опустив в разложении член с s=0 , получим выражение адовлетворяющее () (нтеграл обнуляется в силу ортогональности полиномов с различными s). Выражение в скобках в левой стороне ()

есть. Уравнение () принимает вид


Умножим его с обеих сторон на и проинтегрируем по. Получим систему алгебраических уравнений, которая может быть решена на ЭВМ:


Для последнего выражения введены обозначения


Уравнение с l=0 отсутствует, поскольку в силу сохранения импульса

Коэффициент теплопроводности вычисляется подстановкой выражения () в интеграл (). С учётом условия () интеграл (с) может быть представлен в виде

В результате находим.

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

Заключение.

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

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

Список литературы.

1. Е.М.Лифшиц, Л.П.Питаевский. Физическая кинетика. Наука, М., 1979 г.

2. Ю.Б.Румер, М.Ш.Рывкин. Термодинамика, статистическая физика и кинетика.

Наука, М., 1972 г.


Газа а << 1 «хорошим» кинетическим уравнением является уравнение Больцмана, которое несовместимо с требованием факторизации. Мы видели, что вывод уравнения Больцмана по Боголюбову предполагает только факторизацию функции F2 в «бесконечном прошлом». Рассмотрим случай β = U0/T <<; 1, что соответствует горячему газу со слабым взаимодействием между частицами, который, однако, ...

Планка; NA – число Авогадро;  - приведенная масса молекулы электролита, г; сi – молярная концентрация ионов (сi= c0); с0 – исходная молярная концентрация электролита. 3. Изоморфизм уравнений При рассмотрении движения в поле вязких сил удобно ввести понятие подвижности b. Подвижность определяется как предельная скорость, приобретаемая телом под действием силы, равной единице, т. ...

Интенсивностью электрического поля и величиной градиента температуры, входящими в неоднородный член в левой части. Далее в этой главе мы будем отыскивать решения кинетического уравнения для различных случаев в порядке увеличения сложности. § 7. Электропроводность Пусть на систему наложено только электрическое поле E, и в «бесконечной» среде поддерживается постоянная температура. С учетом...

Серебро Медь Железо Олово Жидкости Ртуть Вода Ацетон Бензол 0 0 0 -3 0 0 0 0 0 20 16 22,5 0,1765 0,1411 0,0237 0,0226 403 86,5 68,2 35,6 0,190 0,167 0,158 6,9 3. Теплопроводность жидкостей и газов Теплопроводность, один из видов переноса теплоты (энергии теплового движения микрочастиц) от более нагретых частей тела к...