玻尔兹曼方程

玻尔兹曼方程或玻尔兹曼输运方程（Boltzmann transport equation，BTE）是由玻尔兹曼于1872年提出的一个方程，用于描述非平衡状态热力学系统的统计行为。具有温度梯度的流体即为这类系统的一个经典的例子：构成流体的微粒在系统中通过随机而具有偏向性的运动让热量从较热的区域流向较冷的区域，而这一过程可用玻尔兹曼方程来描述。在现今的论文中，“玻尔兹曼方程”这个术语常被用于更一般的意义上，它可以是任何涉及描述热力学系统中宏观量（如能量，电荷或粒子数）的变化的动力学方程。
波尔兹曼方程并不去确定流体中每个粒子的位置和动量，而是求出具有特定位置和动量的粒子的概率分布。具体而言，考虑某一瞬间，以位置矢量 formula_1 末端为中心的无穷小区域内，动量无限接近动量矢量 formula_2（即这些粒子在动量空间中也处于无穷小区域 formula_3内）的粒子的概率分布。
波尔兹曼方程可用于确定物理量是如何变化的，例如流体在输运过程中的热能和动量；还可由此推导出其他的流体特征性质，例如黏度，热导率，以及电阻率（将材料中的载流子视为气体），详见对流扩散方程式。
波尔兹曼方程是一个非线性的。方程中的未知函数是一个包含了粒子空间位置和动量的六维概率密度函数。方程解的存在性和唯一性问题仍然没有完全解决，但就最近发表的一些工作而言，对于解决这一问题还是有一定希望的。
概述.
相空间与密度函数.
系统中所有可能的位置 formula_1 和动量 formula_2 组成的集合被称作此系统的相空间，其中位置坐标记为 formula_6，动量坐标记为 formula_7。整个空间是六维的：空间中某一点的坐标可表示为 formula_8，每个坐标均通过时间 formula_9 参数化。微元（或微分体积元）可写作：
波尔兹曼方程的核心是“formula_11”函数，它表示的是在一段极短的时间内，每一相空间单位体积中的formula_12个分子在微元 formula_13 中，位置都为 formula_1 且动量都为 formula_2 的概率。通过定义，我们可使概率密度函数 formula_16 满足以下条件：
formula_18 被定义为在时间 formula_9，位于 formula_20 的空间元 formula_13 中的粒子总数。对坐标空间与动量空间的一个区域积分即可得该区域内所有具有对应位置和动量的粒子的总数：
虽然formula_11是和一群粒子相关的，但此相空间是对于单一粒子的（而不是像多体系统中考虑全部粒子）。这里不使用 r1, p1 表示粒子1，r2, p2 表示粒子2，……，rN, pN 表示粒子N。
系统中的粒子被假定是相同的（因此他们均有相同的质量formula_24）。对于具有超过一种化学组分的混合物，每一种成分都需要有一个分布函数，见下文。
一般形式.
方程的一般形式可以写作：
这里“force”一词指的是外部对粒子施加的力（而不是粒子间的作用），“diff”表示粒子的扩散，“coll”表示粒子的碰撞，指的是碰撞中粒子间相互的作用力。上述三项的具体形式将会在下文给出。
注意，一些作者会使用粒子的速度 formula_26，来代替上文的 formula_2；这两个物理量可以通过定义formula_28来联系。
“force”项和“diff”项.
考虑一群以formula_11分布的粒子。每个粒子均受到外力formula_30的作用（不包括粒子间作用力。粒子间的作用见后面对“coll”项的处理)。
假设在时间 formula_9，一定数量的粒子都有位置 formula_1（于微元 formula_33 内），和动量 formula_2（于微元 formula_35 内）。如果此时有一个力formula_30在这一瞬作用在每个颗粒上，那么在时间 formula_37，它们的位置将会是formula_38，动量将变成 formula_39。在没有碰撞的情况下，formula_11必须满足
这里，注意到相空间元 formula_13 是恒定的这个事实可以从哈密顿方程（见刘维尔定理）得知。然而，由于存在碰撞，相空间元 formula_13 中的粒子密度是可变的，所以
- f(\mathbf{r},\mathbf{p},t)d^3\mathbf{r}d^3\mathbf{p} \\
& = \Delta f d^3\mathbf{r}d^3\mathbf{p} 
\end{align}</math>
其中 formula_44 指的是formula_11的总变化量。式除以 formula_46 并取极限 formula_47 和 formula_48 可得
formula_11的全微分：
\cdot d\mathbf{p} \\
& = \frac{\partial f}{\partial t}dt +\nabla f \cdot \frac{\mathbf{p}dt}{m} + \frac{\partial f}{\partial \mathbf{p}}\cdot \mathbf{F}dt
\end{align}</math>
其中 ∇ 为梯度算符，· 为点积，
是∇的动量类比的一个简写，ê"x", ê"y", ê"z" 为笛卡尔坐标系下的单位矢量。
最终形式.
对两边同除以"dt" 并代入可得：
这里，formula_52 为流体中作用在粒子上的力场，formula_24为粒子质量。 右边的一项用于描述粒子间相互碰撞产生的影响；如果此项为零，则说明粒子之间没有碰撞。无碰撞情况下的玻尔兹曼方程常被称为。
这个方程比上一节“主要论述”中的一般形式更加有用。然而这个方程依旧是不完整的：除非已知formula_11中的碰撞项，否则formula_11是解不出来的。这一项并不像其他项一样可以简单地或一般地得到——这一项是表示粒子的碰撞的统计项，需要知道粒子遵守怎样的统计规律，例如麦克斯韦-玻尔兹曼分布，费米-狄拉克分布或玻色–爱因斯坦分布。
碰撞项（Stosszahlansatz）和分子混沌.
玻尔兹曼的一个关键见解就是对碰撞项的确定。他假设的碰撞项完全是由假定在碰撞前不相关的两个粒子的相互碰撞得到的。这个假设被波尔兹曼称为“Stosszahlansatz”，也叫做“”。根据这一假设，碰撞项可以被写作单粒子分布函数的乘积在动量空间上的积分：
其中 formula_57 和 formula_58 表示碰撞前任意两个粒子的动量（为了方便而标记为formula_59和formula_60）， formula_61 和 formula_62 表示碰撞后的动量
指对应动量的大小（此概念参考相对速度），formula_64 是碰撞的微分散射截面。
对碰撞项的简化.
求解波尔兹曼方程时，许多挑战都来自于其复杂的碰撞项；因此我们会做一些对碰撞项“建模”和简化的尝试。现知最好的模型是由Bhatnagar，Gross和Krook作出的（BGK近似）。BGK近似中假设分子的碰撞会迫使一个物理空间中的某一点的非平衡分布函数回到麦克斯韦平衡分布函数，且其发生率正比于分子碰撞频率。于是，波尔兹曼方程可被写作以下的BGK形式：(也叫做“驰豫时间近似”，relaxation time approximation)
其中 formula_66 是分子碰撞频率，和驰豫时间 formula_67 具有倒数关系：formula_68。formula_69是此处局域的麦克斯韦分布函数，由空间中这一点的气体温度给定。
普适方程（对于混合物）.
对于具有多种化学组分的混合物，我们以 i =1，2，3，……，n 标记各种成分。则对于组分i的方程是：
其中 formula_71。碰撞项为
其中 formula_73，相对动量的大小是
Iij 是粒子i和粒子j之间的微分散射截面。此积分的和描述的是某一相空间元中，组分i粒子的进出。
应用和推广.
守恒方程.
玻尔兹曼方程可用于推导流体动力学中的质量守恒，电量守恒，动量守恒，以及能量守恒定律。对于只含有一种粒子的流体，粒子数密度 formula_75 为：
算符 A 的期望值由下式给出：
由于守恒方程中包含张量，以下使用爱因斯坦求和约定简化标记，即 formula_78 且 formula_79，其中 formula_80 为粒子速度矢量。定义某函数 formula_81，使得其唯一的自变量为动量 formula_82（碰撞中动量守恒）。假设力 formula_83 为位置的函数，且对于 formula_84，formula_11 为0。对玻尔兹曼方程两边同乘 formula_86 ，并对动量积分可得如下四项：
因为 formula_86 在碰撞中守恒，所以最后一项为零。
令 formula_92，即粒子质量，积分后的玻尔兹曼方程化为质量守恒方程：
formula_94 为质量密度，formula_95 为平均流体速度。
令 formula_96，即粒子动量，积分后的玻尔兹曼方程化为动量守恒方程：
formula_98 为压强张量（加上流体静力学压强）。
令 formula_99，即粒子动能，积分后的玻尔兹曼方程化为能量守恒方程：
formula_101 为动力热能密度（kinetic thermal energy density），formula_102 热通量矢量。
哈密顿力学.
在哈密顿力学中, 玻尔兹曼方程通常写作
其中 L 是刘维尔算子(这里定义的刘维尔算子和链接文章中的定义不一致)，它描述了相空间体积的演化；C 是碰撞算子。非相对论下的L 写作
波尔兹曼方程的解.
直到2010年，波尔兹曼方程的准确解才在数学上被证明是（well-behaved）的。这意味着，如果对服从波尔兹曼方程的系统施加一个微扰，此系统最终将回到平衡状态，而不是发散到无穷，或表现出其他的行为。然而，这种存在性证明是无助于我们在现实问题中求解该等式的。 事实上，这个结论只告诉我们某种特定条件下的解是否存在，而不是如何找到他们。在实践中，数值计算方法被用于寻找各种形式的波尔兹曼方程的近似解，应用范围从稀薄气流中的高超音速空气动力学，到等离子体的流动中都可以见到。