工程师平台
首页仿真实验室分子热运动

分子热运动与理想气体仿真

数百个分子刚性弹性碰撞 · 拖动温度体积实时验证 PV = NkT

本站自研 · 秒开
硬球分子动力学
压强由壁面冲量实测
麦克斯韦速率分布
麦克斯韦速率分布(实测柱状 vs 理论曲线)
压强 P(壁面冲量实测)
0.000
P = Σ2m·v⊥ / (t·L)
PV / (NkT)
1.000
理想气体定律应 ≈ 1.000
平均动能 ⟨Ek⟩
0.000
二维 ⟨Ek⟩ = kT ∝ T
均方根速率 v_rms
0.000
v_rms = √(2kT/m)
平均速率 v_avg
0.000
v̄ = √(πkT/2m)
碰壁频率
0 次/s
实测 T = 300 K · V = 1.600
实验模式

四个参数自由调节,观察压强如何响应,副图显示实时速率分布

参数调节
1

拖动右侧活塞压缩气体,压强随之升高

300

温度升高 → 平均动能与速率同步增大(Ek ∝ T)

120

分子越多,压强越大,PV/(NkT) 的统计涨落越小

10

1≈H · 4≈He · 28≈N₂ · 40≈Ar;同温下 v ∝ 1/√m

💡 建议这样玩:
① 温度从 300 拖到 900 —— 平均动能读数正好涨 3 倍
② 切等温模式,体积从 1.0 压到 0.5 —— 压强翻倍,P·V 不变
③ 温度锁死,质量从 4 拖到 40 —— 速率降三分之一,动能不变
④ 分子数从 20 拖到 300 —— 看 PV/(NkT) 涨落越来越小
分子热运动与理想气体原理详解

分子热运动与布朗运动

一切物质都由大量分子组成,这些分子永不停息地做无规则运动,这就是分子热运动。它的直接实验证据是布朗运动:1827 年布朗在显微镜下发现悬浮花粉颗粒不停做折线式颠动,后来才明白这是周围液体分子从各个方向撞击颗粒、每一瞬间合力不为零造成的。颗粒越小、温度越高,布朗运动越剧烈。热运动的两个核心特征是永不停息完全无规则——单个分子的位置和速度无法预测,但巨大数量的分子在统计上呈现出极其稳定的规律,这正是统计物理的立足点。本仿真中每个分子都在做真实的自由飞行与碰撞,你会看到某个分子的轨迹毫无规律,但整体的压强、平均速率读数却相当稳定。

温度的微观本质:温度就是分子平均动能的标志

宏观的温度对应的微观量是分子的平均平动动能。三维气体中⟨Ek⟩ = (3/2)kT,本仿真是二维的,对应⟨Ek⟩ = kT(每个自由度贡献 ½kT,这叫能量均分定理)。由此可解出均方根速率,二维为v_rms = √(2kT/m),三维为 √(3kT/m)。这个式子说明两件事:温度只决定平均动能,与气体种类无关;同一温度下分子质量越大速率越小,按 1/√m 递减。绝对零度的意义也由此而来——那是分子平动动能为零的极限状态,因此温度不可能低于 0 K。

压强的微观起源:分子碰撞壁面的冲量

压强并不是气体「压着」墙的静态推力,而是极多分子频繁撞击器壁的冲量在时间上平均后的效果。一个分子垂直撞墙弹回,动量变化为Δp = 2m·v⊥,把单位时间内所有碰撞的冲量求和再除以面积,就得到压强。气体动理论由此推出P = (2/3)n·⟨Ek⟩(n 为分子数密度,三维),即压强由分子数密度和平均动能共同决定。本仿真严格按这个定义测量压强:累加每次碰壁的 2m|v⊥|,除以统计时间与容器周长,再做指数平滑抑制抖动。因此压强读数是独立测得的物理量,用它去验证 PV = NkT 才有意义。

理想气体状态方程 PV=nRT 与三大气体定律

把压强、体积、温度、物质的量联系起来的就是理想气体状态方程PV = nRT,其微观等价形式为PV = NkT,其中 R = 8.314 J/(mol·K),k = R/N_A = 1.380649×10⁻²³ J/K。三大实验定律都是它的特例:温度不变时P₁V₁ = P₂V₂(玻意耳定律);体积不变时 P/T 为常数(查理定律);压强不变时 V/T 为常数(盖-吕萨克定律)。这个方程建立在三条理想化假设上:分子自身体积可忽略、除碰撞外无相互作用、碰撞完全弹性。高温低压下真实气体与之吻合极好,低温高压下则需范德华方程修正。切到等温模式拖动活塞,副画布上的工作点会沿一条双曲线移动,这就是玻意耳定律的直观形态。

麦克斯韦速率分布与工程应用

同一温度下分子速率并不相同,而是服从麦克斯韦速率分布。二维情形下概率密度为f(v) ∝ v·exp(−mv²/2kT),曲线呈左陡右缓的偏态,存在一个最概然速率 v_p,且恒有 v_p < v̄ < v_rms。温度升高时曲线右移变平,说明高速分子比例增大。这一分布在工程上极其重要:真空技术要按分子平均自由程和碰壁频率来选择泵型与管径;气体分离与同位素浓缩利用不同质量分子速率差异(v ∝ 1/√m)实现扩散分离;内燃机与燃烧工程关心分布高速尾部的分子比例,因为只有动能超过活化能的分子才能引发反应,这也解释了温度对反应速率的指数级影响;气动元件、气弹簧和蓄能器的设计直接用状态方程计算压力行程关系;轮胎冷热压差、高原气压变化、压缩空气储能等日常与工业问题,都能用 PV = nRT 快速估算。

🎓 学生自学
  • • 高中物理与大学「热学」章节预习
  • • 直观理解温度就是平均动能
  • • 亲手验证 PV=NkT 与玻意耳定律
👨‍🏫 教师授课
  • • 课堂现场压缩活塞,替代静态板书
  • • 速率直方图讲清麦克斯韦分布
  • • 演示统计涨落随分子数减小
🔧 工程师参考
  • • 建立气体压强与温度的物理直觉
  • • 理解真空、气动与气弹簧设计前提
  • • 精确校核请用本站工程计算器
常见问题 FAQ

温度的微观本质到底是什么?

温度不是"热量的多少",而是分子平均动能的标志。气体分子永不停息地做无规则热运动,每个分子的速度千差万别,但它们平均动能只由温度决定:二维时 ⟨Ek⟩ = kT,三维时 ⟨Ek⟩ = (3/2)kT,与气体种类、质量、体积都无关。所以温度是一个统计量,对单个分子谈温度没有意义。你可以在上面把温度从 300 K 拖到 900 K,会看到平均动能读数同步涨到 3 倍,分子颜色整体变红、运动明显变剧烈,这就是"温度升高=分子跑得更疯"的直接可视化。

气体压强是怎么产生的?

压强来自数量巨大的分子不断撞击器壁所产生的持续冲量。单个分子垂直撞墙再弹回,动量改变量是 2mv⊥,相当于给墙一次微小的敲击;每秒有天文数字量级的分子在敲,宏观上就表现为一个稳定而连续的压力。本仿真就是这样测压强的:把每次碰壁的 2m|v⊥| 累加起来,除以时间和周长得到 P,而不是套用公式反算。所以你看到的压强读数是"真测出来的",这也让 PV/(NkT) ≈ 1 成为一个真正的验证而非同一恒等式。

为什么同温度下轻分子跑得比重分子快?

因为温度决定的是平均动能而不是平均速率。既然 ½m⟨v²⟩ 对所有分子都相同,那么质量小的分子必须靠更大的速率才能凑够同样的动能,于是 v_rms = √(2kT/m)(二维)随质量按 1/√m 下降。氢气分子质量约为氧气的 1/16,同温下速率就快 4 倍。你可以保持温度不变,把分子质量从 4(氦)拖到 40(氩),会看到平均速率降到约三分之一,但平均动能读数几乎不动。这正是工业上用气体扩散法分离铀同位素、以及氢气容易泄漏的根本原因。

什么是理想气体,真实气体什么时候会偏离?

理想气体是一个模型:分子自身体积可忽略、除碰撞瞬间外分子间无相互作用力、碰撞完全弹性。满足这三条时严格成立 PV = nRT(微观形式 PV = NkT)。真实气体在高温低压下非常接近理想气体,误差常在百分之一以内;但在低温高压下,分子间距被压小,范德华引力和分子自身体积都不能再忽略,就会明显偏离,需要改用范德华方程等修正模型。在本仿真里把体积压到最小、分子数拉到 300,你会看到 PV/(NkT) 略微大于 1,这正是"分子占体积"带来的真实偏差。

麦克斯韦速率分布为什么不是左右对称的?

因为速率有下界没有上界。速率不能为负,v → 0 时分布必须趋于零;而高速端只受指数因子 exp(−mv²/2kT) 压制,理论上可以无限延伸。二维分布 f(v) ∝ v·exp(−mv²/2kT) 是低速的线性上升与高速的指数衰减相乘的结果,因此呈现左陡右缓的偏态长尾。这也导致三个特征速率不重合,恒有 v_p < v̄ < v_rms(副画布中的三条虚线标出了它们)。温度升高时曲线整体右移并变矮变宽,但曲线下面积始终为 1,因为它是归一化的概率密度。

这个仿真的物理计算准确吗,能用于教学吗?

仿真采用真实的硬球分子动力学:每对分子沿连心线交换法向速度分量,严格同时守恒动量与动能;压强由壁面冲量累加实测;速率分布是对当前分子实测分箱后与理论曲线对照,不是画一条预设曲线。因此 PV/(NkT) 通常稳定在 1.0 上下几个百分点,涨落主要来自有限分子数的统计波动——分子数越多越稳,这本身也是统计物理的重要结论。它在你的浏览器本地实时计算,不上传任何数据,可自由用于课堂演示、课前预习与概念验证;若需提交精确数值,仍请用状态方程手算核对。

📌 说明:本仿真由工程师平台自主开发,分子动力学在你的浏览器本地实时计算,不依赖任何第三方服务,不收集任何操作数据。 算法为二维硬球模型:分子与器壁及分子之间均做完全弹性碰撞,等质量两体碰撞沿连心线交换法向速度分量(动量与动能同时守恒),并对重叠做位置修正; 邻近判定采用均匀网格空间哈希(计数排序 + 半邻域遍历),把 N=300 时的 45000 次朴素配对降到每帧千次量级,因此可稳定 60fps。 压强不是套公式反算,而是累加分子碰壁的冲量 Σ2m·v⊥ 后除以统计时间与容器周长实测得到,再做指数平滑。 本页采用二维约定:⟨Ek⟩ = kT,v_rms = √(2kT/m),速率分布 f(v) ∝ v·exp(−mv²/2kT);三维教材中对应为 (3/2)kT 与 √(3kT/m)。 长度、质量、玻尔兹曼常数均使用缩放仿真单位(容器高度取 1,质量以 amu 计,k 取 0.02),故压强、体积为无量纲相对值,但 PV/(NkT) 这一无量纲比值与真实物理完全一致,应稳定在 1.0 附近,其余偏差来自有限分子数的统计涨落与分子自身体积。 仿真用于概念理解与定性验证,如需提交精确数值结果请用状态方程核对。

相关仿真与计算工具