查看“︁玻尔兹曼方程”︁的源代码

{{NoteTA|G1=P}}
{{向量字體常規}}
{{about|描述热力学系统输运行为的方程|统计力学中熵和微观状态数量的关系|波茲曼熵公式|热力学中黑体辐射度与黑体本身的热力学温度T的关系|斯特藩－玻尔兹曼定律|统计力学中描述一定温度下微观粒子运动速度的概率分布|麦克斯韦-玻尔兹曼分布}}
[[File:StairsOfReduction.png|缩略图|300x300像素|波尔兹曼动力学方程在众多近似模型（从微观动力学到宏观连续介质动力学）中所处的位置<ref>{{cite journal |last1=N. Gorban |first1=Alexander |last2=V. Karlin |first2=Ilya |title=Introduction |journal=Invariant Manifolds for Physical and Chemical Kinetics |date=2005 |pages=1–19 |doi=10.1007/978-3-540-31531-5_1 |url=https://www.researchgate.net/publication/27246969_Invariant_manifolds_for_physical_and_chemical_kinetics |publisher=Springer Berlin Heidelberg |language=en |access-date=2016-10-07 |archive-date=2016-04-23 |archive-url=https://web.archive.org/web/20160423193634/https://www.researchgate.net/publication/27246969_Invariant_manifolds_for_physical_and_chemical_kinetics |dead-url=no }}</ref>。]]
'''玻尔兹曼方程'''或'''玻尔兹曼输运方程'''（'''Boltzmann transport equation'''，'''BTE'''）是由[[路德维希·玻尔兹曼|玻尔兹曼]]于1872年提出的一个方程，用于描述非[[热力学平衡|平衡]]状态[[热力学系统]]的统计行为<ref name="Encyclopaediaof">{{cite book |title=Encyclopedia of physics |publisher=VCH |isbn=3-527-26954-1 |edition=2nd}}</ref>。具有[[溫度梯度]]的[[流体]]即为这类系统的一个经典的例子：构成流体的微粒在系统中通过随机而具有偏向性的运动让[[热能|热量]]从较热的区域流向较冷的区域，而这一过程可用玻尔兹曼方程来描述。在现今的论文中，“玻尔兹曼方程”这个术语常被用于更一般的意义上，它可以是任何涉及描述[[热力学系统]]中宏观量（如能量，电荷或粒子数）的变化的动力学方程。

波尔兹曼方程并不去确定流体中每个粒子的[[位置向量|位置]]和[[动量]]，而是求出具有特定位置和动量的粒子的概率分布。具体而言，考虑某一瞬间，以位置矢量 <math>\mathbf{r}</math> 末端为中心的[[无穷小量|无穷小]]区域内，动量无限接近动量矢量 <math>\mathbf{p}</math>（即这些粒子在[[位置空间与动量空间|动量空间]]中也处于无穷小区域 <math>d^3\mathbf{p}</math>内）的粒子的概率分布。

波尔兹曼方程可用于确定[[物理量]]是如何变化的，例如流体在输运过程中的热能和动量；还可由此推导出其他的流体特征性质，例如[[黏度]]，[[熱導率]]，以及[[电阻率]]（将材料中的[[载流子]]视为气体）<ref name="Encyclopaediaof"/>，详见[[漂移–扩散方程|对流扩散方程式]]。

波尔兹曼方程是一个[[非線性]]的{{le|积微分方程|Integro-differential equation}}。方程中的未知函数是一个包含了粒子空间位置和动量的六维[[概率密度函数]]。方程解的[[存在性]]和[[唯一性]]问题仍然没有完全解决，但就最近发表的一些工作而言，对于解决这一问题还是有一定希望的。<ref>{{cite journal |last1=DiPerna |first1=R. J. |last2=Lions |first2=P. L. |title=Ordinary differential equations, transport theory and Sobolev spaces |journal=Inventiones Mathematicae |date=1989-10 |volume=98 |issue=3 |pages=511–547 |doi=10.1007/BF01393835}}</ref><ref>{{cite journal |last1=Gressman |first1=Philip T. |last2=Strain |first2=Robert M. |title=Global classical solutions of the Boltzmann equation without angular cut-off |journal=Journal of the American Mathematical Society |date=2011-09-01 |volume=24 |issue=3 |pages=771–771 |doi=10.1090/S0894-0347-2011-00697-8}}</ref>

== 概述 ==

=== 相空间与密度函数 ===
系统中所有可能的位置 <math>\mathbf{r}</math> 和动量 <math>\mathbf{p}</math> 组成的集合被称作此系统的[[相空间]]，其中位置坐标记为 <math>x,y,z</math>，动量坐标记为 <math>p_x,p_y,p_z</math>。整个空间是六[[維度|维]]的：空间中某一点的坐标可表示为 <math>(\mathbf{r}, \mathbf{p}) = (x,y,z,p_x,p_y,p_z)</math>，每个坐标均通过时间 <math>t</math> [[參數方程|参数化]]。微元（或微分体积元）可写作：
: <math> d^3\mathbf{r}\,d^3\mathbf{p} = dx\,dy\,dz\,dp_x\,dp_y\,dp_z\ </math>
波尔兹曼方程的核心是“<math>f</math>”函数，它表示的是在一段极短的时间内，每一相空间单位体积中的<math>N</math>个分子在微元 <math> d^3\mathbf{r}\,d^3\mathbf{p}</math> 中，位置都为 <math>\mathbf{r}</math> 且动量都为 <math>\mathbf{p}</math> 的概率。通过定义，我们可使概率密度函数 <math>f (\mathbf{r}, \mathbf{p}, t)</math> 满足以下条件：
:<math>dN = f (\mathbf{r}, \mathbf{p}, t)\,d^3\mathbf{r}\,d^3\mathbf{p}</math>

<math>dN</math> 被定义为在时间 <math>t</math>，位于 <math>(\mathbf{r}, \mathbf{p})</math> 的空间元 <math> d^3\mathbf{r}\,d^3\mathbf{p}</math> 中的粒子总数<ref name="Huang1987">{{cite book|author=Kerson Huang|title=Statistical mechanics|url=https://archive.org/details/statisticalmecha0000huan_02ed|year=1987|publisher=Wiley|isbn=978-0-471-81518-1}}</ref>{{rp|61-62}}。对坐标空间与动量空间的一个区域积分即可得该区域内所有具有对应位置和动量的粒子的总数：
:<math>N = \int\limits_\mathrm{positions} d^3\mathbf{r}\int\limits_\mathrm{momenta} d^3\mathbf{p}\,f (\mathbf{r},\mathbf{p},t) = \iiint\limits_\mathrm{positions} \quad \iiint\limits_\mathrm{momenta}  f (x,y,z,p_x,p_y,p_z,t)\,dx\,dy\,dz\,dp_x\,dp_y\,dp_z
</math>

虽然<math>f</math>是和一群粒子相关的，但此相空间是对于'''单一'''粒子的（而不是像多体系统中考虑全部粒子）。这里不使用 '''r'''<sub>1</sub>, '''p'''<sub>1</sub> 表示粒子1，'''r'''<sub>2</sub>, '''p'''<sub>2</sub> 表示粒子2，……，'''r'''<sub>N</sub>, '''p'''<sub>N</sub> 表示粒子N。

系统中的粒子被假定是[[全同粒子|相同]]的（因此他们均有相同的质量<math>m</math>）。对于具有超过一种化学组分的混合物，每一种成分都需要有一个分布函数，见下文。

=== 一般形式 ===
方程的一般形式可以写作：<ref name="McGrawHill">McGraw Hill Encyclopaedia of Physics (2nd Edition), C.B. Parker, 1994, ISBN 0-07-051400-3</ref>
:<math>\frac{d f}{d t} = \left(\frac{\partial f}{\partial t}\right)_\mathrm{force} + \left(\frac{\partial f}{\partial t}\right)_\mathrm{diff}+ \left(\frac{\partial f}{\partial t}\right)_\mathrm{coll}</math>
这里“force”一词指的是外部对粒子施加的力（而不是粒子间的作用），“diff”表示粒子的[[扩散作用|扩散]]，“coll”表示粒子的[[碰撞]]，指的是碰撞中粒子间相互的作用力。上述三项的具体形式将会在下文给出。<ref name="McGrawHill">McGraw Hill Encyclopaedia of Physics (2nd Edition), C.B. Parker, 1994, ISBN 0-07-051400-3</ref>

注意，一些作者会使用粒子的速度 <math>\mathbf{v}</math>，来代替上文的 <math>\mathbf{p}</math>；这两个物理量可以通过定义<math>\mathbf{p} = m\mathbf{v}</math>来联系。

== “force”项和“diff”项 ==
考虑一群以<math>f</math>分布的粒子。每个粒子均受到外力<math>\mathbf{F}</math>的作用（不包括粒子间作用力。粒子间的作用见后面对“coll”项的处理)。

假设在时间 <math>t</math>，一定数量的粒子都有位置 <math>\mathbf{r}</math>（于微元 <math> d^3\mathbf{r}</math> 内），和动量 <math>\mathbf{p}</math>（于微元 <math> d^3\mathbf{p}</math> 内）。如果此时有一个力<math>\mathbf{F}</math>在这一瞬作用在每个颗粒上，那么在时间 <math>t + \Delta\,t</math>，它们的位置将会是<math>\mathbf{r} + \Delta\,\mathbf{r} = \mathbf{r} + \mathbf{p} \Delta\,t/m</math>，动量将变成 <math>\mathbf{p} + \Delta\,\mathbf{p} = \mathbf{p} + \mathbf{F}\Delta\,t</math>。在没有碰撞的情况下，<math>f</math>必须满足


:<math>
f \left (\mathbf{r}+\frac{\mathbf{p}}{m} \Delta t,\mathbf{p}+\mathbf{F}\Delta t,t+\Delta t \right )\,d^3\mathbf{r}\,d^3\mathbf{p} =
f(\mathbf{r},\mathbf{p},t)\,d^3\mathbf{r}\,d^3\mathbf{p}
</math>


这里，注意到相空间元 <math> d^3\mathbf{r}\,d^3\mathbf{p}</math> 是恒定的这个事实可以从[[哈密顿方程]]（见[[刘维尔定理 (哈密顿力学)|刘维尔定理]]）得知。然而，由于存在碰撞，相空间元 <math> d^3\mathbf{r}\,d^3\mathbf{p}</math> 中的粒子密度是可变的，所以
{{NumBlk|:|<math>\begin{align}
dN_\mathrm{coll} & = \left(\frac{\partial f}{\partial t} \right)_\mathrm{coll}\Delta td^3\mathbf{r} d^3\mathbf{p} \\
& = f \left (\mathbf{r}+\frac{\mathbf{p}}{m}\Delta t,\mathbf{p} + \mathbf{F}\Delta t,t+\Delta t \right)d^3\mathbf{r}d^3\mathbf{p}
- 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>
|{{EquationRef|1}}}}
其中 <math>\Delta f</math> 指的是<math>f</math>的总变化量。（{{EquationNote|1}}）式除以 <math> d^3\mathbf{r}\,d^3\mathbf{p}\,\Delta t</math> 并取极限 <math> \Delta t\,\rightarrow 0</math> 和 <math> \Delta f\,\rightarrow 0</math> 可得
{{NumBlk|:|
<math>\frac{d f}{d t} = \left(\frac{\partial f}{\partial t} \right)_\mathrm{coll}</math>
|{{EquationRef|2}}}}
<math>f</math>的全微分：
{{NumBlk|:|
<math>\begin{align}
d f & = \frac{\partial f}{\partial t}dt 
+\left(\frac{\partial f}{\partial x}dx
+\frac{\partial f}{\partial y}dy
+\frac{\partial f}{\partial z}dz
\right)
+\left(\frac{\partial f}{\partial p_x}dp_x
+\frac{\partial f}{\partial p_y}dp_y
+\frac{\partial f}{\partial p_z}dp_z
\right)\\
& = \frac{\partial f}{\partial t}dt +\nabla f \cdot d\mathbf{r} + \frac{\partial f}{\partial \mathbf{p}}\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>
|{{EquationRef|3}}}}
其中''' ∇''' 为[[梯度]]算符，'''·''' 为[[点积]]，

: <math>
\frac{\partial f}{\partial \mathbf{p}} = \mathbf{\hat{e}}_x\frac{\partial f}{\partial p_x} + \mathbf{\hat{e}}_y\frac{\partial f}{\partial p_y}+\mathbf{\hat{e}}_z\frac{\partial f}{\partial p_z}= \nabla_\mathbf{p}f
</math>

是∇的动量类比的一个简写，'''ê'''<sub>''x''</sub>, '''ê'''<sub>''y''</sub>, '''ê'''<sub>''z''</sub> 为[[笛卡尔坐标系]]下的[[单位矢量]]。

=== 最终形式 ===
对（{{EquationNote|3}}）两边同除以''dt'' 并代入（{{EquationNote|2}}）可得：
: <math>\frac{\partial f}{\partial t} + \frac{\mathbf{p}}{m}\cdot\nabla f + \mathbf{F}\cdot\frac{\partial f}{\partial \mathbf{p}} = \left(\frac{\partial f}{\partial t} \right)_\mathrm{coll}</math>
这里，<math>F (\mathbf{r}, t)</math> 为流体中作用在粒子上的[[分子力场|力场]]，<math>m</math>为粒子[[质量]]。 右边的一项用于描述粒子间相互碰撞产生的影响；如果此项为零，则说明粒子之间没有碰撞。无碰撞情况下的玻尔兹曼方程常被称为{{le|弗拉索夫方程式|Vlasov equation}}。

这个方程比上一节“主要论述”中的一般形式更加有用。然而这个方程依旧是不完整的：除非已知<math>f</math>中的碰撞项，否则<math>f</math>是解不出来的。这一项并不像其他项一样可以简单地或一般地得到——这一项是表示粒子的碰撞的'''统计项'''，需要知道粒子遵守怎样的统计规律，例如[[麦克斯韦-玻尔兹曼分布]]，[[费米-狄拉克统计|费米-狄拉克分布]]或[[玻色–爱因斯坦统计|玻色–爱因斯坦分布]]。

== 碰撞项（Stosszahlansatz）和分子混沌 ==
[[路德维希·玻尔兹曼|玻尔兹曼]]的一个关键见解就是对碰撞项的确定。他假设的碰撞项完全是由假定在碰撞前不相关的两个粒子的相互碰撞得到的。这个假设被波尔兹曼称为“Stosszahlansatz”，也叫做“{{le|分子混沌假设|Molecular chaos}}”。根据这一假设，碰撞项可以被写作单粒子分布函数的乘积在动量空间上的积分：<ref name="Encyclopaediaof"/>
: <math>
 \left(\frac{\partial f}{\partial t} \right)_{\mathrm{coll}} = \iint gI(g, \Omega)[f(\mathbf{p'}_A,t) f(\mathbf{p'}_B,t) - f(\mathbf{p}_A,t) f(\mathbf{p}_B,t)] \,d\Omega\,d^3\mathbf{p}_A\,d^3\mathbf{p}_B.
</math>
其中 <math>\mathbf{p}_A</math> 和 <math>\mathbf{p}_B</math> 表示碰撞前任意两个粒子的动量（为了方便而标记为<math>A</math>和<math>B</math>）， <math>\mathbf{p}'_A</math> 和 <math>\mathbf{p}'_B</math> 表示碰撞后的动量
: <math>g = |\mathbf{p}_B - \mathbf{p}_A| = |\mathbf{p'}_B - \mathbf{p'}_A|</math>
指对应动量的大小（此概念参考[[相對速度|相对速度]]），<math>I(g, \Omega)</math> 是碰撞的[[截面 (物理)|微分散射截面]]。

== 对碰撞项的简化 ==
求解波尔兹曼方程时，许多挑战都来自于其复杂的碰撞项；因此我们会做一些对碰撞项“建模”和简化的尝试。现知最好的模型是由Bhatnagar，Gross和Krook作出的（BGK近似）<ref><cite class="citation journal">Bhatnagar, P. L.; Gross, E. P.; Krook, M. (1954-05-01). </cite></ref>。BGK近似中假设分子的碰撞会迫使一个物理空间中的某一点的非平衡分布函数回到麦克斯韦平衡分布函数，且其发生率正比于分子碰撞频率。于是，波尔兹曼方程可被写作以下的BGK形式：(也叫做“驰豫时间近似”，relaxation time approximation<ref>{{Harvnb|Grosso|2014|p=501}}.</ref>)
: <math>\frac{\partial f}{\partial t} + \frac{\mathbf{p}}{m}\cdot\nabla f + \mathbf{F}\cdot\frac{\partial f}{\partial \mathbf{p}} = \nu (f_0 - f)</math>
其中 <math>\nu</math> 是分子碰撞频率，和[[驰豫时间]] <math>\tau</math> 具有倒数关系：<math>\nu = 1/\tau</math>。<math>f_0</math>是此处局域的麦克斯韦分布函数，由空间中这一点的气体温度给定。

== 普适方程（对于混合物） ==
对于具有多种化学组分的[[混合物]]，我们以 i =1，2，3，……，n 标记各种成分。则对于组分i的方程是：<ref name="Encyclopaediaof"/>
: <math>\frac{\partial f_i}{\partial t} + \frac{\mathbf{p}_i}{m_i}\cdot\nabla f_i + \mathbf{F}\cdot\frac{\partial f_i}{\partial \mathbf{p}_i} = \left(\frac{\partial f_i}{\partial t} \right)_\mathrm{coll}</math>
其中 <math>f_i = f_i(\mathbf{r}, \mathbf{p_i}, t)</math>。碰撞项为
: <math>
 \left(\frac{\partial f_i}{\partial t} \right)_{\mathrm{coll}} = \sum_{j=1}^n \iint g_{ij} I_{ij}(g_{ij}, \Omega)[f'_i f'_j - f_if_j] \,d\Omega\,d^3\mathbf{p'}.
</math>
其中 <math>f' = f'(\mathbf{p_i'}, t)</math>，相对动量的大小是
: <math>g_{ij} = |\mathbf{p}_i - \mathbf{p}_j| = |\mathbf{p'}_i - \mathbf{p'}_j|</math>
I<sub>ij</sub> 是粒子i和粒子j之间的微分散射截面。此积分的和描述的是某一相空间元中，组分i粒子的进出。

== 应用和推广 ==
=== 守恒方程 ===

玻尔兹曼方程可用于推导流体动力学中的质量守恒，电量守恒，动量守恒，以及能量守恒定律<ref name="dG1984">{{cite book |last1=de Groot |first1=S.R. |last2=Mazur |first2=P. |title=Non-Equilibrium Thermodynamics |url=http://www.amazon.com/Non-Equilibrium-Thermodynamics-Dover-Books-Physics/dp/0486647412 |accessdate=2013-01-31 |year=1984 |publisher=Dover Publications Inc. |location=New York |isbn=0-486-64741-2 |archive-date=2013-03-28 |archive-url=https://web.archive.org/web/20130328114445/http://www.amazon.com/Non-Equilibrium-Thermodynamics-Dover-Books-Physics/dp/0486647412 |dead-url=no }}</ref>{{rp|p 163}}。对于只含有一种粒子的流体，粒子数密度 <math>n</math> 为：
:<math>n=\int f\,d^3p</math>

算符 A 的期望值由下式给出：
:<math>\langle A \rangle = \frac 1 n \int A f\,d^3p</math>

由于守恒方程中包含[[张量]]，以下使用[[爱因斯坦求和约定]]简化标记，即 <math>\mathbf{x}\rightarrow x_i</math> 且 <math>\mathbf{p}\rightarrow p_i = m w_i</math>，其中 <math>w_i</math> 为粒子速度矢量。定义某函数 <math>g(p_i)</math>，使得其唯一的自变量为动量 <math>p_i</math>（碰撞中动量守恒）。假设力 <math>F_i</math> 为位置的函数，且对于 <math>p_i\rightarrow\pm \infty</math>，<math>f</math> 为0。对玻尔兹曼方程两边同乘 <math>g</math> ，并对动量积分可得如下四项：

:<math>\int g \frac{\partial f}{\partial t}\,d^3p=\frac{\partial }{\partial t} (n \langle g \rangle) </math>

:<math>\int \frac{p_j g}{m}\frac{\partial f}{\partial x_j}\,d^3p=\frac{1}{m}\frac{\partial}{\partial x_j}(n\langle g p_j \rangle)</math>

:<math>\int g F_j \frac{\partial f}{\partial p_j}\,d^3p=-nF_j\left\langle \frac{\partial g}{\partial p_j}\right\rangle</math>

:<math>\int g \left(\frac{\partial f}{\partial t}\right)_\mathrm{coll}\,d^3p=0</math>

因为 <math>g</math> 在碰撞中守恒，所以最后一项为零。


令 <math>g=m</math>，即粒子质量，积分后的玻尔兹曼方程化为质量守恒方程<ref name="dG1984" />{{rp|pp 12,168}}：

:<math>\frac{\partial}{\partial t}\rho
+ \frac{\partial}{\partial x_j}(\rho V_j)
=0</math>

<math>\rho=mn</math> 为质量密度，<math>V_i=\langle w_i\rangle</math> 为平均流体速度。

令 <math>g=m w_i</math>，即粒子动量，积分后的玻尔兹曼方程化为动量守恒方程<ref name="dG1984" />{{rp|pp 15,169}}：

:<math>\frac{\partial}{\partial t}(\rho V_i) + \frac{\partial}{\partial x_j}(\rho V_i V_j+P_{ij}) - nF_i=0</math>

<math>P_{ij}=\rho\langle (w_i-V_i) (w_j-V_j) \rangle</math> 为压强张量（{{le|粘性应力张量|viscous stress tensor}}加上流体静力学[[压强]]）。

令 <math>g=\tfrac{1}{2}m w_i w_i</math>，即粒子动能，积分后的玻尔兹曼方程化为能量守恒方程<ref name="dG1984" />{{rp|pp 19,169}}：

:<math>\frac{\partial}{\partial t}(u+\tfrac{1}{2}\rho V_i V_i) + \frac{\partial}{\partial x_j} (uV_j+\tfrac{1}{2}\rho V_i V_i V_j + J_{qj}+P_{ij}V_i)-nF_iV_i = 0</math>

<math>u=\tfrac{1}{2}\rho\langle (w_i-V_i) (w_i-V_i) \rangle</math> 为动力热能密度（kinetic thermal energy density），<math>J_{qi}=\tfrac{1}{2}\rho\langle (w_i-V_i)(w_k-V_k)(w_k-V_k)\rangle</math> 热通量矢量。

=== 哈密顿力学 ===

在[[哈密顿力学]]中, 玻尔兹曼方程通常写作
:<math>\hat{\mathbf{L}}[f]=\mathbf{C}[f], \, </math>
其中 '''L''' 是[[刘维尔定理 (哈密顿力学)|刘维尔算子]](这里定义的刘维尔算子和链接文章中的定义不一致)，它描述了相空间体积的演化；'''C''' 是碰撞算子。非相对论下的'''L''' 写作

:<math>\hat{\mathbf{L}}_\mathrm{NR} = \frac{\partial}{\partial t} + \frac{\mathbf{p}}{m} \cdot \nabla + \mathbf{F}\cdot\frac{\partial}{\partial \mathbf{p}}\,.</math>

=== 量子理论和粒子数守恒的违背 ===
=== 广义相对论和天文学 ===
== 波尔兹曼方程的解 ==
直到2010年，波尔兹曼方程的准确解才在数学上被证明是{{le|良好|Pathological_(mathematics)#Well-behaved}}（well-behaved）的。这意味着，如果对服从波尔兹曼方程的系统施加一个微扰，此系统最终将回到平衡状态，而不是发散到无穷，或表现出其他的行为<ref><cite class="citation journal">Philip T. Gressman and Robert M. Strain (2010). </cite></ref><ref><cite class="citation web">[https://news.upenn.edu/news/university-pennsylvania-mathematicians-solve-140-year-old-boltzmann-equation-gaseous-behaviors?utm_source=feedburner&utm_medium=twitter&utm_campaign=Feed:%252520delicious/danmack/math%252520delicious/danmack/math "Mathematicians Solve 140-Year-Old Boltzmann Equation"] {{Wayback|url=https://news.upenn.edu/news/university-pennsylvania-mathematicians-solve-140-year-old-boltzmann-equation-gaseous-behaviors?utm_source=feedburner&utm_medium=twitter&utm_campaign=Feed%3A%252520delicious%2Fdanmack%2Fmath%252520delicious%2Fdanmack%2Fmath |date=20180111164938 }}. ''news.upenn.edu''<span class="reference-accessdate">. </span></cite></ref>。然而，这种[[构造性证明|存在性证明]]是无助于我们在现实问题中求解该等式的。 事实上，这个结论只告诉我们某种特定条件下的解是否存在，而不是如何找到他们。在实践中，数值计算方法被用于寻找各种形式的波尔兹曼方程的近似解，应用范围从稀薄气流中的[[高超音速|高超音速空气动力学]]<ref><cite class="citation journal">Evans, Ben; Morgan, Ken; Hassan, Oubay (2011-03-01). </cite></ref>，到[[等离子体]]的流动<ref><cite class="citation journal">Pareschi, L.; Russo, G. (2000-01-01). </cite></ref>中都可以见到。

== 参见 ==
* {{le|弗拉索夫方程|Vlasov equation}}
* [[H定理]]
* [[福克-普朗克方程]]
* [[纳维-斯托克斯方程]]
* {{le|弗拉索夫&#x2013;泊松方程|Vlasov–Poisson equation}}
* {{le|格子波尔兹曼方法|Lattice Boltzmann methods}}

== 注释 ==
<div class="reflist columns references-column-width" style="-moz-column-width: 40em; -webkit-column-width: 40em; column-width: 40em; list-style-type: decimal;">
<references /></div>

== 参考资料 ==
* {{en}}<cite class="citation book">Harris, Stewart (1971). [http://books.google.it/books/about/An_Introduction_to_the_Theory_of_the_Bol.html?id=KfYK1lyq3VYC&redir_esc=y ''An introduction to the theory of the Boltzmann equation''] {{Wayback|url=http://books.google.it/books/about/An_Introduction_to_the_Theory_of_the_Bol.html?id=KfYK1lyq3VYC&redir_esc=y |date=20150402170318 }}. Dover Books. p.&nbsp;221. </cite><cite class="citation book">[[国际标准书号|ISBN]]&nbsp;978-0-486-43831-3.</cite><span class="Z3988" title="ctx_ver=Z39.88-2004&rfr_id=info%3Asid%2Fen.wikipedia.org%3ABoltzmann+equation&rft.aufirst=Stewart&rft.aulast=Harris&rft.btitle=An+introduction+to+the+theory+of+the+Boltzmann+equation&rft.date=1971&rft.genre=book&rft_id=http%3A%2F%2Fbooks.google.it%2Fbooks%2Fabout%2FAn_Introduction_to_the_Theory_of_the_Bol.html%3Fid%3DKfYK1lyq3VYC%26redir_esc%3Dy&rft.isbn=978-0-486-43831-3&rft.pages=221&rft.pub=Dover+Books&rft_val_fmt=info%3Aofi%2Ffmt%3Akev%3Amtx%3Abook">&nbsp;</span>此书的推导始于[[刘维尔定理 (哈密顿力学)|刘维尔定理]]和{{le|BBGKY|BBGKY hierarchy}}，介绍了玻尔兹曼方程在现代物理框架下的地位。其他大部分统计力学的教科书，例如黄克孙的《Statistical Mechanics》，甚至原样照搬玻尔兹曼最初的演算过程。为了简化论述，这些教科书运用启发式的解释，避而不谈玻尔兹曼方程的应用范围以及方程中的某些假设，而这些假设正是玻尔兹曼方程与其他输运方程例如[[福克-普朗克方程]]或朗道方程组的不同之处。

* {{en}}<cite class="citation journal">Arkeryd, Leif (1972). "On the Boltzmann equation part I: Existence". ''Arch. Rational Mech. Anal''. '''45''': 1–16. </cite><cite class="citation journal">[[Bibcode]]:[[bibcode:1972ArRMA..45....1A|1972ArRMA..45....1A]]. [[DOI|doi]]:[[doi:10.1007/BF00253392|10.1007/BF00253392]].</cite><span class="Z3988" title="ctx_ver=Z39.88-2004&rfr_id=info%3Asid%2Fen.wikipedia.org%3ABoltzmann+equation&rft.atitle=On+the+Boltzmann+equation+part+I%3A+Existence&rft.aufirst=Leif&rft.aulast=Arkeryd&rft.date=1972&rft.genre=article&rft_id=info%3Abibcode%2F1972ArRMA..45....1A&rft_id=info%3Adoi%2F10.1007%2FBF00253392&rft.jtitle=Arch.+Rational+Mech.+Anal.&rft.pages=1-16&rft_val_fmt=info%3Aofi%2Ffmt%3Akev%3Amtx%3Ajournal&rft.volume=45">&nbsp;</span>
* {{en}}<cite class="citation journal">Arkeryd, Leif (1972). "On the Boltzmann equation part II: The full initial value problem". ''Arch. Rational Mech. Anal''. '''45''': 17–34. </cite><cite class="citation journal">[[Bibcode]]:[[bibcode:1972ArRMA..45...17A|1972ArRMA..45...17A]]. [[DOI|doi]]:[[doi:10.1007/BF00253393|10.1007/BF00253393]].</cite><span class="Z3988" title="ctx_ver=Z39.88-2004&rfr_id=info%3Asid%2Fen.wikipedia.org%3ABoltzmann+equation&rft.atitle=On+the+Boltzmann+equation+part+II%3A+The+full+initial+value+problem&rft.aufirst=Leif&rft.aulast=Arkeryd&rft.date=1972&rft.genre=article&rft_id=info%3Abibcode%2F1972ArRMA..45...17A&rft_id=info%3Adoi%2F10.1007%2FBF00253393&rft.jtitle=Arch.+Rational+Mech.+Anal.&rft.pages=17-34&rft_val_fmt=info%3Aofi%2Ffmt%3Akev%3Amtx%3Ajournal&rft.volume=45">&nbsp;</span>
* {{en}}<cite class="citation journal">Arkeryd, Leif (1972). "On the Boltzmann equation part I: Existence". ''Arch. Rational Mech. Anal''. '''45''': 1–16. </cite><cite class="citation journal">[[Bibcode]]:[[bibcode:1972ArRMA..45....1A|1972ArRMA..45....1A]]. [[DOI|doi]]:[[doi:10.1007/BF00253392|10.1007/BF00253392]].</cite><span class="Z3988" title="ctx_ver=Z39.88-2004&rfr_id=info%3Asid%2Fen.wikipedia.org%3ABoltzmann+equation&rft.atitle=On+the+Boltzmann+equation+part+I%3A+Existence&rft.aufirst=Leif&rft.aulast=Arkeryd&rft.date=1972&rft.genre=article&rft_id=info%3Abibcode%2F1972ArRMA..45....1A&rft_id=info%3Adoi%2F10.1007%2FBF00253392&rft.jtitle=Arch.+Rational+Mech.+Anal.&rft.pages=1-16&rft_val_fmt=info%3Aofi%2Ffmt%3Akev%3Amtx%3Ajournal&rft.volume=45">&nbsp;</span>
* {{en}}<cite class="citation journal">DiPerna, R. J.; Lions, P.-L. (1989). "On the Cauchy problem for Boltzmann equations: global existence and weak stability". </cite><cite class="citation journal">''Ann. of Math. (2)''. '''130''': 321–366. [[DOI|doi]]:[[doi:10.2307/1971423|10.2307/1971423]].</cite><span class="Z3988" title="ctx_ver=Z39.88-2004&rfr_id=info%3Asid%2Fen.wikipedia.org%3ABoltzmann+equation&rft.atitle=On+the+Cauchy+problem+for+Boltzmann+equations%3A+global+existence+and+weak+stability&rft.aufirst=R.+J.&rft.aulast=DiPerna&rft.au=Lions%2C+P.-L.&rft.date=1989&rft.genre=article&rft_id=info%3Adoi%2F10.2307%2F1971423&rft.jtitle=Ann.+of+Math.+%282%29&rft.pages=321-366&rft_val_fmt=info%3Aofi%2Ffmt%3Akev%3Amtx%3Ajournal&rft.volume=130">&nbsp;</span>

== 外部链接 ==
* {{en}}[http://homepage.univie.ac.at/franz.vesely/sp_english/sp/node7.html The Boltzmann Transport Equation by Franz Vesely] {{Wayback|url=http://homepage.univie.ac.at/franz.vesely/sp_english/sp/node7.html |date=20160924000636 }}
* {{en}}[https://web.archive.org/web/20151123214334/http://www.upenn.edu/pennnews/news/university-pennsylvania-mathematicians-solve-140-year-old-boltzmann-equation-gaseous-behaviors Boltzmann gaseous behaviors solved]
[[Category:偏微分方程]]
[[Category:统计力学]]
[[Category:傳輸現象]]
[[Category:热力学]]