Интегральные уравнения Больцмана

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

История развития

Предпосылки возникновения

В середине XIX века классическая термодинамика успешно описывала равновесные состояния, но не могла объяснить процессы перехода к равновесию. Джеймс Клерк Максвелл в 1867 году вывел распределение скоростей молекул в равновесном газе, однако динамика неравновесных систем оставалась неясной. Людвиг Больцман, развивая идеи Максвелла, в 1872 году опубликовал работу «Дальнейшие исследования теплового равновесия между молекулами газа», где впервые сформулировал кинетическое уравнение.

Формулировка интегрального уравнения

Больцман предложил уравнение, описывающее изменение функции распределения \( f(\mathbf{r}, \mathbf{v}, t) \) в фазовом пространстве координат и скоростей. В интегральной форме оно учитывает изменение числа частиц в элементе фазового объёма за счёт столкновений. Оператор столкновений \( \mathcal{C}[f] \) записывается как интеграл по скоростям и углам рассеяния, что делает уравнение интегро-дифференциальным. В 1916 году Сидни Чепмен и Дэвид Энског независимо разработали методы решения этого уравнения, что привело к созданию теории Чепмена — Энскога.

Развитие в XX веке

В 1930-х годах Николай Боголюбов предложил цепочку уравнений ББГКИ (Боголюбова — Борна — Грина — Кирквуда — Ивона), которая сводит многочастичные задачи к кинетическим уравнениям. В 1940-х годах Юджин Вигнер ввёл понятие функции распределения в квантовой механике, что привело к квантовым аналогам уравнений Больцмана. В 1950-х годах для плазмы были разработаны уравнения Власова — Больцмана, учитывающие коллективные взаимодействия. Современные вычислительные методы, такие как метод прямого моделирования Монте-Карло (DSMC), позволяют решать интегральные уравнения Больцмана для сложных задач аэродинамики и физики плазмы.

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

Основное уравнение

Интегральное уравнение Больцмана в безразмерной форме имеет вид: \[ \frac{\partial f}{\partial t} + \mathbf{v} \cdot \nabla_{\mathbf{r}} f + \frac{\mathbf{F}}{m} \cdot \nabla_{\mathbf{v}} f = \mathcal{C}[f], \] где \( f(\mathbf{r}, \mathbf{v}, t) \) — функция распределения, \( \mathbf{F} \) — внешняя сила, \( m \) — масса частицы. Оператор столкновений \( \mathcal{C}[f] \) в интегральной форме: \[ \mathcal{C}[f] = \int \int \left( f' f'_1 - f f_1 \right) \sigma(\Omega) \, |\mathbf{v} - \mathbf{v}_1| \, d\Omega \, d\mathbf{v}_1, \] где \( f = f(\mathbf{r}, \mathbf{v}, t) \), \( f_1 = f(\mathbf{r}, \mathbf{v}_1, t) \), \( f' = f(\mathbf{r}, \mathbf{v}', t) \), \( f'_1 = f(\mathbf{r}, \mathbf{v}'_1, t) \) — функции распределения до и после столкновения, \( \sigma(\Omega) \) — дифференциальное сечение рассеяния, \( d\Omega \) — элемент телесного угла.

Свойства оператора столкновений

Оператор столкновений обладает несколькими ключевыми свойствами:

• Сохранение числа частиц: \( \int \mathcal{C}[f] \, d\mathbf{v} = 0 \).
• Сохранение импульса: \( \int m \mathbf{v} \mathcal{C}[f] \, d\mathbf{v} = 0 \).
• Сохранение энергии: \( \int \frac{1}{2} m v^2 \mathcal{C}[f] \, d\mathbf{v} = 0 \).
• H-теорема: \( \frac{dH}{dt} \leq 0 \), где \( H = \int f \ln f \, d\mathbf{r} \, d\mathbf{v} \), что соответствует возрастанию энтропии.

Интегральные формы

Для стационарных задач (например, в теории переноса) уравнение Больцмана может быть сведено к интегральному уравнению Фредгольма второго рода: \[ f(\mathbf{r}, \mathbf{v}) = f_0(\mathbf{r}, \mathbf{v}) + \int K(\mathbf{r}, \mathbf{v}, \mathbf{r}', \mathbf{v}') f(\mathbf{r}', \mathbf{v}') \, d\mathbf{r}' \, d\mathbf{v}', \] где \( f_0 \) — свободный член, учитывающий начальные или граничные условия, а \( K \) — ядро, описывающее столкновения.

Классификация и виды

По типу взаимодействия

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

По форме записи

• Интегро-дифференциальная форма: классическая запись с производными по времени и координатам.
• Интегральная форма: для стационарных задач, сводится к уравнению Фредгольма.
• Линеаризованная форма: для малых отклонений от равновесия, где оператор столкновений линеаризуется, что упрощает анализ.

По применимости

• Уравнение Больцмана — Власова: для плазмы с учётом самосогласованного электромагнитного поля, где столкновения пренебрежимо малы.
• Уравнение Больцмана — Ландау: для плазмы с кулоновскими столкновениями, где оператор столкновений заменяется на интеграл Ландау.
• Уравнение Больцмана — Фоккера — Планка: для систем с малыми случайными возмущениями, например, в теории броуновского движения.

Методы решения

Аналитические методы

• Метод Чепмена — Энскога: разложение функции распределения в ряд по малому параметру (числу Кнудсена), позволяет получить гидродинамические уравнения (Навье — Стокса, теплопроводности) из кинетического уравнения.
• Метод моментов: система уравнений для моментов функции распределения (плотность, скорость, температура), замыкаемая через приближения (например, метод Грэда).
• Метод преобразования Фурье: для линеаризованных задач, сводит интегральное уравнение к алгебраическому в пространстве волновых векторов.

Численные методы

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

Применение

Аэродинамика и газовая динамика

Интегральные уравнения Больцмана используются для расчёта течений разреженных газов, где число Кнудсена велико (например, в вакуумной технике, при моделировании верхних слоёв атмосферы Земли, в микрофлюидике). Метод DSMC применяется для проектирования гиперзвуковых летательных аппаратов (например, российского ракетного комплекса «Авангард»).

Физика плазмы

В термоядерном синтезе уравнения Больцмана — Власова и Больцмана — Ландау описывают поведение плазмы в токамаках (например, в российском проекте ИТЭР). Они используются для анализа неустойчивостей, переноса частиц и тепла.

Транспортные явления

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

Астрофизика

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

Критика и ограничения

Допущения модели

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

Вычислительная сложность

Интегральный оператор столкновений имеет размерность \( 6N \), где \( N \) — число частиц, что делает прямое численное решение крайне ресурсоёмким. Для реальных задач (например, с \( 10^{23} \) частицами) используются упрощённые модели, такие как приближение БГК (Бхатнагара — Гросса — Крука), где оператор столкновений заменяется на релаксационный член.

Альтернативные подходы

В современной физике для описания неравновесных систем применяются методы молекулярной динамики, решёточные уравнения Больцмана (LBM) и цепочки ББГКИ. LBM, в частности, широко используется в вычислительной гидродинамике для моделирования течений в пористых средах и многофазных систем.

Интересные факты

• В 1872 году Больцман вывел H-теорему, которая стала математическим обоснованием второго начала термодинамики. Однако в то время она подвергалась критике со стороны физиков, включая Эрнста Маха, который отрицал реальность атомов.
• Уравнение Больцмана является основой для расчёта теплозащиты космических аппаратов при входе в атмосферу Земли. Например, для российского корабля «Союз» используются модели, основанные на DSMC.
• В 2010 году российский математик Григорий Перельман, известный доказательством гипотезы Пуанкаре, также работал над вопросами, связанными с кинетическими уравнениями, в частности, с уравнением Больцмана.
• Решёточное уравнение Больцмана (LBM) было разработано в 1988 году и с тех пор стало популярным методом в вычислительной гидродинамике, особенно для моделирования течений в сложных геометриях.

Источники

• Больцман Л. «Лекции по теории газов». — М.: Гостехиздат, 1956.
• Чепмен С., Каулинг Т. «Математическая теория неоднородных газов». — М.: Иностранная литература, 1960.
• Боголюбов Н. Н. «Проблемы динамической теории в статистической физике». — М.: Гостехиздат, 1946.
• Коган М. Н. «Динамика разреженного газа». — М.: Наука, 1967.
• Ландау Л. Д., Лифшиц Е. М. «Физическая кинетика». — М.: Наука, 1979.
• Бёрд Г. «Молекулярная газовая динамика». — М.: Мир, 1981.
• Сухинов А. А. «Кинетические уравнения и их приложения». — Ростов-на-Дону: Изд-во ЮФУ, 2015.