模拟设置
粒子数 N20
粒子半径 (px)10
恢复系数 e1.0
温度（初始速度）2.0
温度 (平均动能)
—
总能量
—
碰撞次数
0
能量守恒率
—
箱体内N个粒子相互碰撞的二维模拟。自由调节恢复系数、重力和粒子数，观察能量守恒和麦克斯韦-玻尔兹曼速度分布的形成过程。
模拟的核心是二维刚体粒子之间的碰撞处理。对于两个质量分别为 $m_1$, $m_2$ 的粒子，碰撞遵循动量守恒定律，但动能是否守恒取决于恢复系数 $e$。
$$ \begin{aligned}m_1 \vec{v}_1 + m_2 \vec{v}_2 &= m_1 \vec{v}_1' + m_2 \vec{v}_2' \quad \text{(动量守恒)}\\[5pt] e &= \frac{v_{2n}' - v_{1n}'}{v_{1n}- v_{2n}}\quad \text{(恢复系数定义)}\end{aligned}$$其中，$\vec{v}$ 和 $\vec{v}'$ 分别是碰撞前、后的速度矢量，下标 $n$ 表示沿两粒子中心连线方向（法向）的速度分量。当 $e=1$ 时，动能也守恒；$e<1$ 时，部分动能在碰撞中耗散。
系统达到热平衡后，粒子速度的统计分布服从麦克斯韦-玻尔兹曼分布。在二维情况下，速度大小 $v$ 的概率密度函数为：
$$ f(v) = \frac{m}{k_B T}v \, e^{-mv^2 / (2k_B T)} $$这里 $m$ 是粒子质量，$k_B$ 是玻尔兹曼常数，$T$ 是系统的热力学温度。这个分布描述了在给定温度下，具有某一速度的粒子出现的概率。模拟器右侧的分布图就是对这个理论分布的验证。
粉体工程与离散元法（DEM）：这是CAE中直接相关的应用。比如在模拟制药厂的药粉混合过程、农业中的谷物输送，或者水泥厂的原料搅拌时，工程师使用DEM软件。软件的核心就是计算数百万个微小颗粒之间像本模拟器这样的碰撞，从而预测粉体的流动、分离和混合效率。
气体动力学与航空航天：在高空稀薄气体中，空气分子的行为就类似于本模拟的粒子。研究卫星受到的微弱大气阻力，或设计火星探测器的进入舱时，需要用到基于此类分子碰撞理论的直接模拟蒙特卡洛（DSMC）方法。
材料加工与制造：在金属粉末冶金中，将金属粉末压制成型前，需要理解粉末的填充和流动特性。通过DEM模拟调整颗粒大小、形状和碰撞属性（恢复系数），可以优化模具设计，减少缺陷，提高零件密度和强度。
基础物理与化学教学：这个模拟器完美地可视化了从微观粒子碰撞到宏观热力学定律（如温度、压强）的桥梁。它帮助学生直观理解“温度本质是分子平均动能”以及平衡态统计分布是如何从动力学中自发涌现出来的。
首先，要注意“粒子数量越多越接近现实”的常见误解。虽然增加粒子数量确实能让模拟场景看起来更丰富，但计算负荷会呈爆炸式增长。例如，将粒子数从100增加到1000时，碰撞检测的计算量将增长约100倍！在实际工程应用中，基本原则是从能够捕捉现象本质的最小粒子数开始。其次，关于恢复系数e的设置：e=1代表“完全弹性”，e=0代表“完全非弹性”，但现实世界中几乎不存在“完全”状态。即使是金属球的碰撞，e值也在0.9左右；而像黏土这类材料，e值往往低于0.1。在理解教科书定义后，如何使该参数更贴近实际情况至关重要。最后，关于“启用重力后温度会下降？”的疑问：当粒子受重力下落并在底部聚集时，碰撞频率增加可能导致动能看似以热能形式“损失”。但这并非能量守恒定律的失效，其主要原因在于粒子势能与动能的转换，以及容器壁碰撞（这里恢复系数同样起作用！）导致的能量耗散。建议养成分析系统整体能量平衡的习惯。
本模拟器核心的粒子碰撞计算，构成了DEM（离散元法）数值分析方法的基石。DEM是粉体工程领域不可或缺的工具，例如用于预测制药厂中原料药与辅料的混合均匀性，或设计农业料斗中谷物顺畅流动的形态。更进一步，其原理与SPH（光滑粒子流体动力学）等粒子法模拟相通。SPH将流体行为建模为粒子集合而非连续介质，其应用范围涵盖飞溅水流、金属铸造过程，甚至天体物理学中的星系形成模拟。此外，它与分子动力学（MD）模拟在概念上具有亲缘性。MD虽计算原子/分子层面的相互作用，但同样基于牛顿运动方程和与此处碰撞计算类似的成对势函数力计算。可以说，通过本工具学习的“成对粒子碰撞”原理，已成为从微观到宏观所有基于粒子的数值模拟的通用语言。
建议首先从数学层面追踪“守恒定律”入手：在模拟器中开启“速度矢量显示”，记录特定两个粒子碰撞前后的动量（质量×速度矢量），验证矢量之和是否真正守恒。对于弹性碰撞（e=1），还可计算动能之和（$ \frac{1}{2}m_1v_1^2 + \frac{1}{2}m_2v_2^2 $）。这将帮助您直观理解数学公式与实际运动的关联。在此基础上，学习统计力学基础能深化认知：为何麦克斯韦-玻尔兹曼分布呈现特定形态（$ v^2 \exp(-v^2) $）？其关键在于“等概率原理”与“熵增原理”。建议通过在线资源检索这些关键词。最后，推荐探索“非球形粒子”与“摩擦/凝聚力的引入”等实践主题。现实中的粉体并非球体，且存在摩擦旋转和静电凝聚等现象。如何在DEM中建模这些效应是重要课题。掌握本工具基础后，若进一步研究“粒子旋转运动”“范德华力”等附加作用力，必将大幅拓展对模拟世界的理解。