Mathematica模拟理想气体分子的热运动.docVIP

Mathematica模拟理想气体分子的热运动 ——论文 姓名：黄静 学院：物电学院 班级：1202班 学号：2012112010214 引言：自从计算机诞生后，科学家们努力尝试着从基本的力学过程来模拟分子和原子的热运动特性，这方面的成果是很多的，但是这些模拟计算太复杂了，因为它们往往要假设分子之间力的模型，还要追踪分子的位置，以及对是否碰撞给出判据，所推到导出的碰撞后分子的速度表示也十分复杂；对于分子与器壁的碰撞，要么没有涉及，要么假定分子撞击器壁之后以原能量返回，器壁的温度作用体现不出来。这样的工作，能揭示的信息不多，不太容易被大学生所接受。而这次我们将介绍一个新的模拟方法，，它不对作用力的细节进行假定，只关注分子碰撞的结果，同时可以考虑分子与器壁的碰撞，将器壁温度对气体能量的控制作用体现出来，能多方面考察分子热运动特性，是一种不错的模拟方案。 同种分子碰撞理论： 解得： 另一个分子的动量为： 器壁的作用：分子与器壁的碰撞元比分子之间的碰撞复杂，它不仅是不同种类的分子（原子）之间的碰撞，而且碰撞过程的能量交换还与温度有关。可以考虑一个简化的模型：分子与器壁的碰撞，使得分子丧失了原来的动量和能量，结果，分子以一个随机的动量和能量被反弹或放射。考虑到器壁的温度，可以设散射分子的动量有一个上限P。，这个上限与温度有关，根据固体热运动理论，动量的这个上限的平方该与温度成正比。 实验内容： 1.器壁的温度一定，分子初始动量较小 这时候，分子要从器壁活的能量，系统的总动量和平均能量都要随着碰撞次数的增多而增大，最后达到与器壁的热平衡，其标志就是分子的平均能量趋于稳定。分子数可以假设，实验取分子数n=1000，m为碰撞次数，Pt为一个分子的初始动量，P。为器壁给予分子的最大动量。Moment：分子的动量列表；position：预备抽取进行碰撞的分子序号；energy：分子能量列表；stat：体积一些状态下分子能量平均值，没隔50次碰撞统计一次这个次数可以依模拟分子的多少、碰撞次数的多少等情况而改变。 该图像告诉我们，在气体分子初始动量较小的情况下，分子要从器壁吸收能量，分子的平均能量要随时间增大，表现为能量曲线是上升的。经过与器壁的充分碰撞，气体分子最终与器壁达到能量平衡，能量曲线就达到水平状态，稍微有一些涨落。这时，气体的平均能量为33.78左右。 2. 器壁的温度一定，分子初始动量较大 模拟程序与上面的一样，要改变的是参数Pt，Pt=300，P。=10维持不变，得到图像： 如图，在开始的一段时间内，能量高的气体分子并没有雨器壁碰撞，其它能量小的分子反而从与器壁的碰撞中获得了能量，因而分子系统的平均能量石升高的；但是，当一个能量高的气体分子与器壁碰撞，丧失了自己的高能量，系统的平均能量就突然降低了。这样的过程，随后还会发生很多次，使分子系统的平均能量出现了几个跳跃式的下降，这样可以更快的趋于热平衡。 小结：我们首先导出了三维空间里理想气体分子碰撞的动量分配原理，在此基础上所设计的模拟方法，把器壁的温度成功地考虑了进来，用计算机模拟了气体分子热平衡的过程，证明了分子之间的碰撞时导致分子动量分布的根本原因，并验证了分子按动量的分布律很好地符合经典分布律。本模拟方法的理论很简单，图像清楚，抓住了分子运动和碰撞的主要特征，同时，也很好的运用了计算机，对mathematica的编程更加的了解了。所以，学好mathematica，对于提高物理研究水平有着实实在在的帮助。