为了计算更大尺度体系的导热系数,使用多粒子碰撞动力学(Multi-Particle Collision Dynamics, MPCD)方法模拟流体的运动。首先研究流体粗粒子质量、盒子尺寸、体系温度和旋转角度对导热系数计算结果的影响,然后采用MPCD方法计算Ar原子体系和水分子体系的导热系数,计算时间仅为3 min左右,相比分子动力学(Molecular Dynamics, MD)方法,大大减少了计算量。