气体动理论硬核得分笔记

一、 气体动理论核心考点大纲

二、 高数运算指南与速度空间映射

气体动理论的本质是对速度空间（相空间）的概率密度函数进行多重定积分。考场上直接决定计算速度与准确率的核心在于熟练掌握雅可比行列式换元与高斯积分体系。

速度空间向标量速率空间的测度映射，依赖于微元体积的几何变换。在三维空间中，直角坐标系速度微元 d^3v = dv_x dv_y dv_z 必须映射到速率 v = |\vec{v}| 上。依据空间各向同性假设，物理系统的概率分布仅依赖于速率大小，通过对球坐标系立体角 d\Omega = \sin\theta d\theta d\phi 的积分，可将三维微元转化为球壳体积 4\pi v^2 dv 。在二维吸附气体模型中，微元映射为圆环面积 2\pi v dv ；在一维模型中，微元即为 dv_x。这种升维与降维映射是处理不同物理场景分布函数归一化的前置条件。

对于所有包含 e^{-\alpha v^2} 形式的积分，必须摒弃考场上的现场分部积分法，直接调用欧拉 - 泊松积分（高斯积分）的广义结论 。定义参数 \alpha = \frac{m}{2kT}，所有速度特征量均由母积分 I_n(\alpha) = \int_0^\infty v^n e^{-\alpha v^2} dv 决定。利用含参变量积分求导法与 Gamma 函数 \Gamma(x)，可得出普适解析解：I_n = \frac{1}{2} \alpha^{-\frac{n+1}{2}} \Gamma\left(\frac{n+1}{2}\right)。

常考幂次积分结果呈现严格的奇偶交替规律。偶数次幂积分与高斯曲线的面积、方差（能量）直接相关：

奇数次幂积分与系统的质心、平均绝对速率直接相关，其解析解不含 \pi：

[!example]- 高阶积分推导机制与考场防错校验

核心推导：基于 \Gamma(x+1) = x\Gamma(x) 及 \Gamma(1/2) = \sqrt{\pi}，推导 I_4。

I_4 = \frac{1}{2} \alpha^{-\frac{4+1}{2}} \Gamma\left(\frac{5}{2}\right) = \frac{1}{2} \alpha^{-5/2} \cdot \frac{3}{2} \cdot \frac{1}{2} \cdot \Gamma\left(\frac{1}{2}\right) = \frac{3}{8}\sqrt{\frac{\pi}{\alpha^5}}

降维推导法（Feynman Trick）：通过对参数 \alpha 连续求导降维。

I_2 = -\frac{\partial}{\partial \alpha} I_0 = -\frac{\partial}{\partial \alpha} \left( \frac{1}{2} \pi^{1/2} \alpha^{-1/2} \right) = \frac{1}{4} \sqrt{\frac{\pi}{\alpha^3}}

I_4 = -\frac{\partial}{\partial \alpha} I_2 = -\frac{\partial}{\partial \alpha} \left( \frac{1}{4} \pi^{1/2} \alpha^{-3/2} \right) = \frac{3}{8} \sqrt{\frac{\pi}{\alpha^5}}

偶数阶全空间积分校验：

注意题干积分限！若积分为 \int_{-\infty}^\infty，由于偶函数性质，结果为 I_{2n} 的两倍。

\int_{-\infty}^\infty e^{-\alpha v_x^2} dv_x = 2 I_0 = \sqrt{\frac{\pi}{\alpha}}

若积分限为 [-\infty, \infty] 的奇数阶积分，由于奇函数对称性，积分严格为 0。

三、 核心公式定理与物理升维映射

1. 微观压强公式与温度本质

理想气体的宏观压强是海量微观粒子与器壁发生弹性碰撞时，法向动量转移的时间与空间平均效应。微观动量通量通过牛顿第二定律与宏观压强建立映射关系 。

温度的微观本质是分子无规热运动剧烈程度的统计度量，且仅取决于分子的平均平动动能，与气体种类、质量无关。

综合可得理想气体状态方程微观表达形式 P = n k T，其中 n = \frac{N}{V} 为分子数密度，k 为玻尔兹曼常数 。

2. 能量均分定理与内能计算

处于热力学平衡态的系统，其总能量均匀分配在每一个独立平方项对应的自由度上，每一个自由度贡献 \frac{1}{2}kT 的平均能量 。理想气体的内能仅为其所有分子的动能（平动、转动、振动）之和，忽略分子间势能。

混合气体的内能满足线性叠加原理。系统的总内能计算公式为：

其中 \nu_j 为第 j 种组分的摩尔数。迈耶公式 C_{P,m} - C_{V,m} = R 揭示了定压加热过程中气体膨胀做功的额外能量消耗。在极高温度下，双原子分子的振动能级被激发，振动自由度 v=1 包含动能与势能两个平方项，将额外贡献 kT 的能量，此时等效自由度变为 t+r+2v = 7 。

3. 麦克斯韦速度与速率分布律

平衡态下，分子速度各个分量的取值互不影响，且空间没有择优方向。三维速度分布函数为三个一维高斯分布的乘积 。

经过向球坐标映射积分 4\pi v^2 dv，得到麦克斯韦速率分布函数，其表示单位速率区间内的分子数比例 f(v) = \frac{1}{N} \frac{dN}{dv}：

由此通过微积分推导出表征群体运动状态的三大统计特征速率：

三者大小顺序严格维持 v_p < \bar{v} < v_{rms}，其比值不随气体种类与温度变化，恒定为 1 : 1.128 : 1.224 。

4. 玻尔兹曼分布律

在保守外场中，粒子不仅具有动能分布，其在空间不同势能 \varepsilon_p(x,y,z) 处的数密度分布亦服从指数衰减律 。

当外场为均匀重力场且 \varepsilon_p = mgz 时，由态方程 P = nkT 即可推导出等温大气压强公式（Barometric formula）：

5. 碰撞频率与平均自由程

分子热运动过程中频繁碰撞，碰撞截面 \sigma = \pi d^2 决定了单分子在空间中扫过的等效碰撞管体积。考虑到两分子相互靠近时的相对运动，其相对速率均值为平动速率均值的 \sqrt{2} 倍 。平均碰撞频率 \bar{Z} 与平均自由程 \bar{\lambda} 公式为：

等温条件下，自由程与压强成严格反比关系；等容条件下，自由程与温度完全无关。

6. 输运现象的微观机制

当系统偏离平衡态但存在局部宏观梯度时，分子通过跨越平均自由程的热运动实现质量、动量和能量的转移 。三种输运现象的唯象定律与微观统计参数建立严格映射： 扩散系数 D （浓度梯度驱使质量输运，Fick 定律）：

黏滞系数 \eta （速度梯度驱使动量输运，Newton 黏滞定律）：

热导率 \kappa （温度梯度驱使能量输运，Fourier 导热定律）：

四、 极致高频母题库（高能计算预警）

母题 1：任意概率分布函数的截断与特征量求解

题目描述：设某容器内气体分子的速率分布函数在低速区域呈抛物线截断：当 0 \leq v \leq v_0 时，f(v) = C v^2 (v_0 - v)；当 v > v_0 时 f(v) = 0。（1）利用归一化条件求解常数 C；（2）计算该系统的平均速率 \bar{v} 与方均根速率 v_{rms}；（3）计算处于 \frac{1}{2}v_0 \sim v_0 速率区间的分子百分比。

破题思路：任何自定义分布函数的破题点均在于 \int f(v)dv = 1。截断函数将积分上限从 \infty 压制为截断点 v_0，随后分别计算 \int v f(v)dv 与 \int v^2 f(v)dv 即可完成高维物理量降维。

最终答案：（1）C = \frac{12}{v_0^4}；（2）\bar{v} = \frac{3}{5}v_0, v_{rms} = \sqrt{\frac{2}{5}}v_0；（3）31.25\%。

[!example]- 详细计算步骤与易错点分析

第一步：利用归一化条件定出常数 C。

\int_0^\infty f(v) dv = \int_0^{v_0} C (v_0 v^2 - v^3) dv = C \left[ v_0 \frac{v^3}{3} - \frac{v^4}{4} \right]_0^{v_0} = C \left( \frac{v_0^4}{3} - \frac{v_0^4}{4} \right) = C \frac{v_0^4}{12} = 1

解得：C = \frac{12}{v_0^4}。

第二步：求解平均速率 \bar{v}。

\bar{v} = \int_0^{v_0} v f(v) dv = \frac{12}{v_0^4} \int_0^{v_0} (v_0 v^3 - v^4) dv = \frac{12}{v_0^4} \left[ v_0 \frac{v^4}{4} - \frac{v^5}{5} \right]_0^{v_0} = \frac{12}{v_0^4} \left( \frac{v_0^5}{4} - \frac{v_0^5}{5} \right) = \frac{12}{v_0^4} \frac{v_0^5}{20} = \frac{3}{5} v_0

第三步：求解方均根速率 v_{rms}。

先求方均速率 \overline{v^2}：

\overline{v^2} = \int_0^{v_0} v^2 f(v) dv = \frac{12}{v_0^4} \int_0^{v_0} (v_0 v^4 - v^5) dv = \frac{12}{v_0^4} \left[ v_0 \frac{v^5}{5} - \frac{v^6}{6} \right]_0^{v_0} = \frac{12}{v_0^4} \left( \frac{v_0^6}{5} - \frac{v_0^6}{6} \right) = \frac{12}{v_0^4} \frac{v_0^6}{30} = \frac{2}{5} v_0^2

开方得到：v_{rms} = \sqrt{\frac{2}{5}} v_0。

第四步：计算特定速率区间的分子比例。

\frac{\Delta N}{N} = \int_{v_0/2}^{v_0} f(v) dv = \frac{12}{v_0^4} \left[ \frac{v_0 v^3}{3} - \frac{v^4}{4} \right]_{v_0/2}^{v_0}

上限代入：\frac{v_0^4}{3} - \frac{v_0^4}{4} = \frac{v_0^4}{12}。

下限代入：\frac{v_0 (v_0/2)^3}{3} - \frac{(v_0/2)^4}{4} = v_0^4 \left(\frac{1}{24} - \frac{1}{64} \right)= v_0^4 \frac{8-3}{192} = \frac{5 v_0^4}{192}。

比例值：\frac{12}{v_0^4} \left(\frac{v_0^4}{12} - \frac{5 v_0^4}{192} \right)= 1 - \frac{12 \times 5}{192} = 1 - \frac{60}{192} = 1 - \frac{5}{16} = \frac{11}{16} = 0.6875。

易错点：多项式积分下限不为零时，必须严格带入求差，极易在分式通分环节发生算术错误。

母题 2：二维麦克斯韦分布与表面吸附物理模型

题目描述：在金属表面物理实验中，单原子分子气体吸附在固体表面形成二维理想气体层（分子只能在 x-y 表面内自由平动）。系统温度为 T，原子质量为 m。（1）推导二维麦克斯韦速率分布函数 f_2(v)；（2）计算该二维系统的最概然速率 v_{p,2D}。 破题思路：二维速度空间的微元为圆环面积 2\pi v dv。采用等能面上的二维几率密度 C e^{-\frac{m(v_x^2+v_y^2)}{2kT}}，通过极坐标测度积分完成归一化。 最终答案：（1）f_2(v) = \frac{m}{kT} v e^{-\frac{mv^2}{2kT}}；（2）v_{p,2D} = \sqrt{\frac{kT}{m}}。

[!example]- 详细计算步骤与易错点分析

第一步：建立二维速度分布并归一化。

二维相空间分布概率正比于能量指数函数：F(v_x, v_y) dv_x dv_y = A e^{-\frac{m(v_x^2+v_y^2)}{2kT}} dv_x dv_y。

转换为极坐标 v^2 = v_x^2 + v_y^2，dv_x dv_y = v dv d\theta。积分消去角度项 \int_0^{2\pi} d\theta = 2\pi，速率微元变为 2\pi v dv。

速率分布函数为：f_2(v) dv = 2\pi A v e^{-\frac{mv^2}{2kT}} dv。

采用归一化条件：

\int_0^\infty f_2(v) dv = 2\pi A \int_0^\infty v e^{-\frac{mv^2}{2kT}} dv = 1

根据高斯积分公式 I_1(\alpha) = \frac{1}{2\alpha}，此处 \alpha = \frac{m}{2kT}，故积分值为 \frac{kT}{m}。

因此 2\pi A \cdot \frac{kT}{m} = 1 \implies 2\pi A = \frac{m}{kT}。

得到二维速率分布函数：f_2(v) = \frac{m}{kT} v e^{-\frac{mv^2}{2kT}}。

第二步：求解最概然速率。

对速率分布函数求导并令其等于零：

\frac{df_2(v)}{dv} = \frac{m}{kT} \left( 1 \cdot e^{-\frac{mv^2}{2kT}} + v \left(-\frac{mv}{kT}\right) e^{-\frac{mv^2}{2kT}} \right) = 0

提取公因式得到 1 - \frac{mv^2}{kT} = 0 \implies v_{p,2D} = \sqrt{\frac{kT}{m}}。

物理映射：三维的最概然速率为 \sqrt{\frac{2kT}{m}}，二维系统降低为 \sqrt{\frac{kT}{m}}，一维系统由于高斯分布峰值位于零点，其最概然速率为 0。空间维度的塌缩直接导致速度分布向低速区漂移。

母题 3：真空隙流（Effusion）的分子通量微观计算

题目描述：在温度为 T、压强为 P 的绝热容器壁上开有一极小孔洞，孔面积为 A（孔径远小于分子平均自由程，从而不破坏平衡分布）。推导单位时间内逸出小孔的分子数 \Phi 的理论表达式 。 破题思路：只考虑速度在 x > 0（向着孔洞）方向的分子。取速度相空间体积元进行全立体角与全速率域的半球面微积分。 最终答案：\Phi = \frac{1}{4} n \bar{v} A = \frac{P A}{\sqrt{2\pi m k T}}。

[!example]- 详细计算步骤与易错点分析

第一步：建立相空间分布函数 。 设孔面与 x 轴垂直，由于孔极小，内部气体仍服从麦克斯韦速度分布：

dN = n \left( \frac{m}{2\pi kT} \right)^{3/2} e^{-\frac{m(v_x^2+v_y^2+v_z^2)}{2kT}} dv_x dv_y dv_z

单位时间内，速度为 \vec{v} 的分子，能穿过小孔的有效空间是一个底面积为 A、斜高为 v dt、法向高为 v_x dt 的斜柱体（且 v_x > 0）。

单位时间逸出的小孔通量微元：d\Phi = A v_x dN。

第二步：速度空间半球面积分。

换成球坐标系 (v, \theta, \phi)，其中 \theta 为速度矢量与 x 轴正向（法向）的夹角。

v_x = v \cos\theta，空间体积元 dv_x dv_y dv_z = v^2 \sin\theta dv d\theta d\phi。

麦克斯韦速率分布函数 f(v) dv = 4\pi \left(\frac{m}{2\pi kT} \right)^{3/2} v^2 e^{-mv^2/2kT} dv。

考虑到空间各向同性，在 d\Omega = \sin\theta d\theta d\phi 立体角范围内的分子占比为 \frac{d\Omega}{4\pi}。

d\Phi = A (v \cos\theta) \left[ n f(v) dv \frac{\sin\theta d\theta d\phi}{4\pi} \right]

对能飞出小孔的所有分子在半空间积分（v \in [0, \infty)，\theta \in [0, \pi/2]，\phi \in [0, 2\pi]）：

\Phi = \int_0^{2\pi} d\phi \int_0^{\pi/2} \sin\theta \cos\theta d\theta \int_0^\infty A n \frac{1}{4\pi} v f(v) dv

分离变量逐个击破：

1. \int_0^{2\pi} d\phi = 2\pi

2. \int_0^{\pi/2} \sin\theta \cos\theta d\theta = \left[\frac{1}{2} \sin^2\theta \right]_0^{\pi/2} = \frac{1}{2}

3. \int_0^\infty v f(v) dv = \bar{v}

   三者连乘：\Phi = 2\pi \times \frac{1}{2} \times \frac{1}{4\pi} A n \bar{v} = \frac{1}{4} n \bar{v} A。

第三步：转化为宏观可测参数表达。

代入态方程密度 n = \frac{P}{kT}，以及平均速率 \bar{v} = \sqrt{\frac{8kT}{\pi m}}：

\Phi = \frac{1}{4} \frac{P}{kT} \sqrt{\frac{8kT}{\pi m}} A = \frac{P A}{\sqrt{2\pi m k T}}

此公式是真空检漏器校准、航天器真空漏气速率计算的基础理论底座。

母题 4：同位素分离的级联网络与逸出气体的动能陷阱

题目描述：铀 -235 与铀 -238 的六氟化铀（UF_6）气体在混合容器中，通过多孔膜扩散（可视为多个独立隙流）。（1）证明两种同位素通过单层多孔膜的理论分离系数 \alpha = \sqrt{M_{238} / M_{235}}；（2）计算逸出小孔的气体分子的平均平动动能 \bar{\varepsilon}_{esc}，并与容器内部气体平均平动动能 \frac{3}{2}kT 进行严格对比 。 破题思路：基于母题 3 的流量公式 \Phi \propto 1/\sqrt{m} 给出分离系数。逸出分子的平均动能计算时，权重分布已突变为正比于 v \cdot f(v) 的通量分布。 最终答案：（1）流量比等于质量平方根的反比；（2）\bar{\varepsilon}_{esc} = 2kT（比内部气体高出 \frac{1}{2}kT）。

[!example]- 详细计算步骤与易错点分析

第一步：证明理论分离系数 。 两种同位素气体的逸出透过率分别为： \Phi_{235} = \frac{A P_{235}}{\sqrt{2\pi m_{235} k T}}，\Phi_{238} = \frac{A P_{238}}{\sqrt{2\pi m_{238} k T}}。 渗透后的浓度比（流量比）：\frac{\Phi_{235}}{\Phi_{238}} = \frac{P_{235}}{P_{238}} \sqrt{\frac{m_{238}}{m_{235}}}。 分离系数定义为前后浓度比值之比：\alpha = \frac{\Phi_{235} / \Phi_{238}}{P_{235} / P_{238}} = \sqrt{\frac{m_{238}}{m_{235}}}。这是二战曼哈顿计划中庞大的级联气体扩散分离工厂的核心方程。

第二步：推导逸出气体的平均动能（极致防错考点，难度极高）。

容器内部署麦克斯韦分布 f(v)，静止系平均动能 \int \frac{1}{2}mv^2 f(v)dv = \frac{3}{2}kT。

但在逸出小孔处，速度越快的分子，单位时间内击中孔面的几何概率越大。因此，逸出分子群中高速分子的比例被硬性抬高。

其分布函数重整化为通量加权分布：

f_{esc}(v) dv = \frac{v f(v) dv}{\int_0^\infty v f(v) dv} = \frac{v f(v)}{\bar{v}} dv

逸出分子的平均平动动能积分构建：

\bar{\varepsilon}_{esc} = \int_0^\infty \left(\frac{1}{2}mv^2\right) f_{esc}(v) dv = \frac{m}{2\bar{v}} \int_0^\infty v^3 f(v) dv

代入麦克斯韦速率分布 f(v) = 4\pi (\frac{m}{2\pi kT})^{3/2} v^2 e^{-\frac{mv^2}{2kT}}：

核心定积分项变为 \int_0^\infty v^5 e^{-\frac{mv^2}{2kT}} dv。

调用前述“高数运算指南”的奇数次高斯积分 I_5(\alpha)，其中 \alpha = \frac{m}{2kT}：

\int_0^\infty v^5 e^{-\alpha v^2} dv = \frac{1}{\alpha^3} = \left(\frac{2kT}{m}\right)^3

提取常数项处理全式：

\int_0^\infty v^3 f(v) dv = 4\pi \left( \frac{m}{2\pi kT} \right)^{3/2} \left(\frac{2kT}{m}\right)^3 = \frac{4\pi}{(2\pi)^{3/2}} \left( \frac{2kT}{m} \right)^{3/2} = \frac{4}{\sqrt{\pi}} \left(\frac{2kT}{m}\right)^{3/2}

除以平均速率 \bar{v} = \sqrt{\frac{8kT}{\pi m}} = \frac{2}{\sqrt{\pi}} \left(\frac{2kT}{m} \right)^{1/2}：

\bar{\varepsilon}_{esc} = \frac{m}{2} \cdot \frac{ \frac{4}{\sqrt{\pi}} (\frac{2kT}{m})^{3/2} }{ \frac{2}{\sqrt{\pi}} (\frac{2kT}{m})^{1/2} } = \frac{m}{2} \cdot 2 \left( \frac{2kT}{m} \right) = 2kT

物理升维判据：逸出气体的平均平动动能为 2kT，明确指向一个非直觉结论——伴随隙流漏气发生，系统流失的每一分子平均携带的能量大于容器内部留守分子的平均内能。因此，在绝热条件下发生漏气，容器内部剩余气体的温度将发生剧烈下降！此机制即为前沿冷原子实验中磁阱蒸发冷却（Evaporative Cooling）的核心动力学基础 。

母题 5：平均自由程、碰撞频率与 ​\sqrt{2} 相对速度修正机制

题目描述：空气在 T=300 \text{K}，大气压 P = 1.01 \times 10^5 \text{Pa} 下，空气分子有效直径 d=0.35 \text{nm}。（1）计算分子数密度 n；（2）严格证明考量分子相对运动时，同种气体平均碰撞频率中 \sqrt{2} 因子的微观数学来源 。 破题思路：基于矢量合成法则与速度各向同性统计规律，严格推导均方相对速率并开方。 最终答案：（1）n = 2.44 \times 10^{25} \text{m}^{-3}；（2）证明见展开说明。

[!example]- 详细计算步骤与易错点分析

第一步：求解分子数密度。

根据态方程 P V = N k T，得出 n = \frac{N}{V} = \frac{P}{kT}。

极限代数代入：n = \frac{1.01 \times 10^5}{1.38 \times 10^{-23} \times 300} \approx 2.44 \times 10^{25} \text{m}^{-3}。

第二步：推导平均碰撞频率中的 \sqrt{2} 因子（全网顶尖高校防满分陷阱）。

设两个碰撞分子的速度向量分别为 \vec{v}_1, \vec{v}_2，其相对速度矢量定义为 \vec{v}_{rel} = \vec{v}_1 - \vec{v}_2。

相对速率平方的统计系综平均值为：

\overline{v_{rel}^2} = \overline{(\vec{v}_1 - \vec{v}_2) \cdot (\vec{v}_1 - \vec{v}_2)} = \overline{v_1^2} - 2\overline{\vec{v}_1 \cdot \vec{v}_2} + \overline{v_2^2}

由于大量分子在空间中进行各向同性的无规热运动，两分子的速度方向完全独立且不相关，因此速度点乘的统计平均必定为零：\overline{\vec{v}_1 \cdot \vec{v}_2} = 0。

对于同种气体，统计性质等价，存在 \overline{v_1^2} = \overline{v_2^2} = \overline{v^2}。

故得到均方关联：\overline{v_{rel}^2} = 2\overline{v^2}。

根据麦克斯韦分布函数的普遍性质，无论绝对速率还是相对速率的模，均满足特定的高斯类映射关系，可以严密过渡得到算术平均的线性映射（也可通过严格对二维相对速度分布积分得出）：

\overline{v_{rel}} = \sqrt{2}\bar{v}

第三步：推导平均自由程 \bar{\lambda}。

截面积为 \sigma = \pi d^2。在时间 \Delta t 内，单一分子相对于背景气体扫过的有效碰撞圆柱体体积为 V_{col} = \pi d^2 \overline{v_{rel}} \Delta t。

扫过体积内的目标分子均被视为发生碰撞，碰撞次数 Z \Delta t = n V_{col} = n \pi d^2 \sqrt{2} \bar{v} \Delta t。

从而平均碰撞频率：\bar{Z} = \sqrt{2} \pi d^2 n \bar{v}。

平均自由程的物理定义为连续两次碰撞间沿绝对轨道的平均路程：

\bar{\lambda} = \frac{\text{绝对参考系飞行距离}}{\text{碰撞总次数}} = \frac{\bar{v} \Delta t}{\bar{Z} \Delta t} = \frac{\bar{v}}{\sqrt{2} \pi d^2 n \bar{v}} = \frac{1}{\sqrt{2} \pi d^2 n}

易错点警示：极易在分子处将绝对飞行距离误写为 \overline{v_{rel}} \Delta t。平均自由程考查的是相对于绝对空间坐标系的质心位移，因此分子真实绝对飞行速率必定是 \bar{v}，由于分子被迎面而来的高速相对通量网“提前拦截”，碰撞频率增加，导致最终自由程分母出现 \sqrt{2}，路程被缩短。若气体采用范德瓦尔斯刚球模型进一步修正，公式常添加额外修正项 l = 1 - (d/2) 。

母题 6：混合气体的自由程网络与折合质量折减

题目描述：设某化学反应器中有组分 A 和组分 B 的混合气体。分子密度分别为 n_A, n_B；有效碰撞直径分别为 d_A, d_B；单分子质量分别为 m_A, m_B。建立 A 分子在混合气体中的总碰撞频率 \bar{Z}_A 与总平均自由程 \bar{\lambda}_A 理论方程 。 破题思路：A 分子的总碰撞频率是其与同类 A 分子碰撞的频率、加上与异类 B 分子碰撞频率的线性叠加。处理 A-B 异种碰撞时，必须重构有效碰撞直径，并基于热平衡下动能均分重算相对速率。 最终答案：\bar{\lambda}_A = \frac{\bar{v}_A}{\bar{Z}_{AA} + \bar{Z}_{AB}}，完整表达见展开项。

[!example]- 详细计算步骤与易错点分析

第一步：计算同种碰撞频率 \bar{Z}_{AA}。

依据单组分结论直接写出：\bar{Z}_{AA} = \sqrt{2} \pi d_A^2 n_A \bar{v}_A。

第二步：推导异种碰撞频率 \bar{Z}_{AB}（高度综合考验）。

A 碰到 B 时，两球心间的极限距离为两者半径之和：r_{AB} = r_A + r_B = \frac{d_A + d_B}{2}。

有效碰撞截面积为 \sigma_{AB} = \pi r_{AB}^2 = \pi \left(\frac{d_A + d_B}{2}\right)^2。（绝对不能错误地将 d_A, d_B 平方后平均）。

相对速度均方值的重算依赖于独立组分热平衡：

两气体处于同温 T 混合，均方速率分别为 \overline{v_A^2} = \frac{3kT}{m_A}，\overline{v_B^2} = \frac{3kT}{m_B}。

异种分子的速度叉乘统计平均依然为零 \overline{\vec{v}_A \cdot \vec{v}_B} = 0。

相对速度均方：

\overline{v_{rel,AB}^2} = \overline{v_A^2} + \overline{v_B^2} = \frac{3kT}{m_A} + \frac{3kT}{m_B} = 3kT \left(\frac{m_A + m_B}{m_A m_B}\right)

引入经典力学的折合质量 \mu = \frac{m_A m_B}{m_A + m_B}，相对速度大小统计遵循折合麦克斯韦分布，由此得出相对平均速率：

\overline{v_{rel,AB}} = \sqrt{\frac{8kT}{\pi \mu}} = \sqrt{ \frac{8kT}{\pi} \left( \frac{1}{m_A} + \frac{1}{m_B} \right) }

巧妙利用单组分平均速率关系 \bar{v}_A = \sqrt{\frac{8kT}{\pi m_A}} 与 \bar{v}_B = \sqrt{\frac{8kT}{\pi m_B}}，将上式重组为极简形式：

\overline{v_{rel,AB}} = \sqrt{ \bar{v}_A^2 + \bar{v}_B^2 } = \bar{v}_A \sqrt{1 + \frac{m_A}{m_B}}

因此 A 撞击 B 的频率项为：\bar{Z}_{AB} = n_B \pi \left(\frac{d_A + d_B}{2}\right)^2 \bar{v}_A \sqrt{1 + \frac{m_A}{m_B}}。

第三步：组合求总平均自由程。

A 分子经历的整体碰撞频率 \bar{Z}_A = \bar{Z}_{AA} + \bar{Z}_{AB}。

\bar{Z}_A = \bar{v}_A \left

最终得到极为对称的异种混合物平均自由程：

\bar{\lambda}_A = \frac{\bar{v}_A}{\bar{Z}_A} = \frac{1}{ \sqrt{2} \pi d_A^2 n_A + \pi \left(\frac{d_A + d_B}{2}\right)^2 n_B \sqrt{ 1 + \frac{m_A}{m_B} } }

此推导彻底将动量、几何截面与统计几率融合，常作为研究生入学考试的压轴大题。

母题 7：非惯性系中的玻尔兹曼分布（高速离心机）

题目描述：一根长度为 L 的封闭圆柱形空心管，内部充有摩尔质量为 M_{mol} 的理想气体，温度恒定为 T。管子绕经过其一端且与管轴垂直的轴，以恒定角速度 \omega 旋转。若管中心转轴处的压强为 P_0，求管子末端处的压强 P_L。

破题思路：转动系统中存在随径向距离增加而下降的离心势能场。应用玻尔兹曼密度分布律替换势能项，结合理想气体状态方程即可求解。

最终答案：P_L = P_0 e^{\frac{M_{mol} \omega^2 L^2}{2RT}}。

[!example]- 详细计算步骤与易错点分析

第一步：建立旋转离心势能场。

在随管旋转的非惯性系中，质量为 m 的单个气体分子在距离轴心径向距离 r 处受到离心力 F_c = m \omega^2 r（方向向外）。

离心力属于保守力，可定义离心势能 \varepsilon_p(r)。取中心转轴处 r=0 势能为零：

\varepsilon_p(r) = - \int_0^r F_c dr = - \int_0^r m \omega^2 r dr = - \frac{1}{2} m \omega^2 r^2

易错点：离心势能是负值！因为质点受力方向指向势能减小的方向，即离心力向外排斥，势能向外递减。

第二步：应用玻尔兹曼分布律。

分子数密度分布函数为：

n(r) = n_0 e^{-\frac{\varepsilon_p(r)}{kT}} = n_0 e^{\frac{m \omega^2 r^2}{2kT}}

这表明，随着半径 r 增大，气体被“甩”向外侧，密度呈现指数型激增（而非重力场中的指数型衰减）。

第三步：计算末端宏观压强。

结合局部态方程 P(r) = n(r) kT 与 P_0 = n_0 kT，可得到：

P(r) = P_0 e^{\frac{m \omega^2 r^2}{2kT}}

将分子质量 m 转换为摩尔质量 M_{mol}，利用 k = R/N_A 且 m N_A = M_{mol}：

将末端坐标 r=L 代入方程，获得最终压强：

P_L = P_0 e^{\frac{M_{mol} \omega^2 L^2}{2RT}}

此模型在医学生物化学与核工业中被广泛用于超高速离心机分离大分子蛋白或重同位素沉淀。

母题 8：混合气体系统等容吸热与自由度计算

题目描述：绝热刚性容器内装有 2 \text{mol} 的氮气（刚性双原子分子）与 1 \text{mol} 的氩气（单原子分子）混合。初始温度为 T_0。内部电热丝通电向气体释放焦耳热 Q。求（1）混合气体的等效定容摩尔热容 C_{V,mix} 及等效自由度 i_{mix}；（2）系统平衡后的末态温度 T_f 与压强 P_f。

破题思路：应用能量按自由度均分定理线性叠加计算系统总内能，反推等效自由度。结合热力学第一定律求温度。

最终答案：（1）i_{mix} = \frac{13}{3}，C_{V,mix} = \frac{13}{6}R；（2）T_f = T_0 + \frac{2Q}{13R}。

[!example]- 详细计算步骤与易错点分析

第一步：计算各组分独立内能与自由度。

氩气 Ar（单原子）：自由度 i_1 = 3，摩尔数 \nu_1 = 1，内能 U_1 = \nu_1 \frac{i_1}{2} RT = \frac{3}{2}RT。

氮气 N_2（双原子）：常温下为刚性，自由度 i_2 = 5，摩尔数 \nu_2 = 2，内能 U_2 = \nu_2 \frac{i_2}{2} RT = 2 \times \frac{5}{2}RT = 5RT。

混合气体总内能 U = U_1 + U_2 = \frac{13}{2}RT。

第二步：推导混合等效统计参数。 系统总摩尔数 \nu_{total} = \nu_1 + \nu_2 = 3 \text{mol}。 假设等效自由度为 i_{mix}，强制满足：U = \nu_{total} \frac{i_{mix}}{2} RT = 3 \frac{i_{mix}}{2} RT = \frac{13}{2}RT。 解出分数维自由度：i_{mix} = \frac{13}{3}。 定容摩尔热容为：C_{V,mix} = \frac{i_{mix}}{2}R = \frac{13}{6}R。 物理映射：混合气体的等效自由度常常出现分数，它不再具备单一分子空间平移与旋转的明确几何构型意义，仅代表热力学能量吸收权重的纯粹统计平均 。

第三步：计算吸热引发的态函数变化。

热力学第一定律微元形式：dQ = dU + dW。刚性容器体积恒定，对外做功 dW=0。

宏观热量积分：Q = \Delta U = \nu_{total} C_{V,mix} \Delta T = 3 \times \frac{13}{6} R (T_f - T_0) = \frac{13}{2} R (T_f - T_0)。

解得末态温度 T_f = T_0 + \frac{2Q}{13R}。

根据总态方程 P_f V = \nu_{total} R T_f = 3 R \left(T_0 + \frac{2Q}{13R} \right)，展开即得压强：P_f = \frac{3R T_0}{V} + \frac{6Q}{13V}。

母题 9：输运现象的动量通量（平板黏滞系数）

题目描述：在无限大相距为 L 的平行平板间充满压强为 P、温度为 T 的气体（分子质量为 m）。下极板静止，上极板以极小速度 u_0 沿 x 轴移动。气体层达到稳定的层流速度分布 u(z) = u_0 \frac{z}{L}。推导气体的内部黏滞系数 \eta 理论表达式 。 破题思路：计算穿过任一水平面 z 的分子上下穿梭所携带的宏观定向动量流（横向）。结合平均自由程推导通量，根据牛顿黏滞定律比对系数。 最终答案：\eta = \frac{1}{3} n m \bar{v} \bar{\lambda}。

[!example]- 详细计算步骤与易错点分析

第一步：建立跨层动量转移模型。

任取一高度为 z 的气层，其宏观流速为 u(z)。分子不仅做无规热运动（速度大小均值为 \bar{v}），还附加了宏观漂移速度 u。

考虑自上而下（从 z + \bar{\lambda} 处）穿过 z 面的分子。由于最后一次碰撞发生在 z + \bar{\lambda} 处，它们携带的是该层的宏观动量：m u(z + \bar{\lambda})。

单位时间内向下的分子通量为：\Phi_{\downarrow} = \frac{1}{6} n \bar{v}（利用三分之一自由度在 z 轴，正负向各半的各向同性简化近似模型，更严密积分给出 1/4 因子，但经典教科书此处均采用简化的 1/6 定性模型定出数量级）。

向下传输的总宏观动量流：J_{\downarrow} = \frac{1}{6} n \bar{v} \cdot m u(z + \bar{\lambda})。

第二步：计算净动量通量。

同理，自下而上（从 z - \bar{\lambda} 处）穿过 z 面的分子携带动量向上。向上的动量通量为：

J_{\uparrow} = \frac{1}{6} n \bar{v} \cdot m u(z - \bar{\lambda})。

z 面受到的向下的净动量流（即每层施加给下方流体层的摩擦切应力 \tau）：

\tau = J_{\downarrow} - J_{\uparrow} = \frac{1}{6} n \bar{v} m \left[ u(z + \bar{\lambda}) - u(z - \bar{\lambda}) \right]

第三步：泰勒展开定系数。

运用一阶泰勒展开，因为 \bar{\lambda} 极小：

u(z + \bar{\lambda}) - u(z - \bar{\lambda}) \approx \left(u(z) + \frac{du}{dz}\bar{\lambda} \right)- \left( u(z) - \frac{du}{dz}\bar{\lambda} \right)= 2 \bar{\lambda} \frac{du}{dz}。

代回原式：

\tau = \frac{1}{6} n m \bar{v} \left( 2 \bar{\lambda} \frac{du}{dz} \right) = \frac{1}{3} n m \bar{v} \bar{\lambda} \frac{du}{dz}

与牛顿黏性定律宏观形式 \tau = \eta \frac{du}{dz} 比对，直接剥离出微观参数：

\eta = \frac{1}{3} n m \bar{v} \bar{\lambda}

防错陷阱：气体黏滞系数随温度升高而增加（正比于 \sqrt{T}，因 \bar{v} \propto \sqrt{T} 且 \bar{\lambda} 恒定），这与液体（随温度升高黏度下降）的直观经验完全相反！因为气体黏性的本质是分子热扩散导致的动量交换，热运动越剧烈，动量交换越强。

母题 10：能量分布函数与最概然能量的非线性漂移

题目描述：已知麦克斯韦速率分布函数，推导分子平动动能 \varepsilon 的分布函数 f(\varepsilon)。并计算最概然动能 \varepsilon_p，证明它不等于最概然速率对应的动能 \frac{1}{2} m v_p^2 。 破题思路：利用微元概率相等的换元原则 f(\varepsilon) d\varepsilon = f(v) dv，将自变量 v 替换为 \varepsilon。求极值寻找最概然能量。 最终答案：f(\varepsilon) = \frac{2}{\sqrt{\pi}} (kT)^{-3/2} \varepsilon^{1/2} e^{-\frac{\varepsilon}{kT}}，最概然动能 \varepsilon_p = \frac{1}{2} kT。

[!example]- 详细计算步骤与易错点分析

第一步：变量替换生成动能分布函数。

已知动能公式 \varepsilon = \frac{1}{2} m v^2 \implies v = \sqrt{\frac{2\varepsilon}{m}}。

微元对应关系：d\varepsilon = m v dv \implies dv = \frac{d\varepsilon}{mv} = \frac{d\varepsilon}{\sqrt{2m\varepsilon}}。

速率分布函数为：f(v) = 4\pi \left(\frac{m}{2\pi kT} \right)^{3/2} v^2 e^{-\frac{mv^2}{2kT}}。

动能分布函数为：f(\varepsilon) = f(v) \frac{dv}{d\varepsilon}。

代入映射：

f(\varepsilon) = 4\pi \left( \frac{m}{2\pi kT} \right)^{3/2} \left( \frac{2\varepsilon}{m} \right) e^{-\frac{\varepsilon}{kT}} \cdot \frac{1}{\sqrt{2m\varepsilon}}

对系数进行约分合并：

f(\varepsilon) = \frac{4\pi}{(2\pi)^{3/2}} \frac{m^{3/2}}{(kT)^{3/2}} \frac{2}{m} \frac{1}{\sqrt{2m}} \varepsilon^{1/2} e^{-\frac{\varepsilon}{kT}} = \frac{2}{\sqrt{\pi}} \frac{1}{(kT)^{3/2}} \varepsilon^{1/2} e^{-\frac{\varepsilon}{kT}}

第二步：对能量求导寻找极值。

要求最概然能量 \varepsilon_p，令 \frac{df(\varepsilon)}{d\varepsilon} = 0。

\frac{d}{d\varepsilon} \left = \frac{1}{2} \varepsilon^{-1/2} e^{-\frac{\varepsilon}{kT}} - \frac{1}{kT} \varepsilon^{1/2} e^{-\frac{\varepsilon}{kT}} = 0

提取非零项 e^{-\frac{\varepsilon}{kT}} \varepsilon^{-1/2}，得：

\frac{1}{2} - \frac{\varepsilon}{kT} = 0 \implies \varepsilon_p = \frac{1}{2} kT

物理谬误粉碎：

考场上最容易想当然地认为最概然动能就等于最概然速率对应的动能。

然而，最概然速率对应的动能是：\varepsilon(v_p) = \frac{1}{2} m (\sqrt{\frac{2kT}{m}})^2 = kT。

计算结果 \varepsilon_p = \frac{1}{2} kT \neq \varepsilon(v_p)！

这是因为概率密度函数在非线性换元（v \to v^2）时，坐标轴被非均匀拉伸，dv 宽度的微元在映射到 d\varepsilon 时引发了概率“峰值”的左移扭曲。

五、 终极考前 Cheat Sheet