周文静
,
栾辉宝
,
孙杰
,
何雅玲
,
陶文铨
工程热物理学报
本文研究了广泛应用于耦合算法中的非周期性边界条件(NPBC),获得了耦合算法中粒子区边界附近原子所受边界力与边界处局部密度和温度的关系。基于此,本文拟合得到了边界力与局部密度和温度的关系式,与平衡态分子动力学模拟结果的对比验证了拟合公式的可靠性。最后本文将边界力的拟合公式应用到耦合算法中求解Poiseuille流动,结果表明粒子区边界附近的密度振荡大大减小,且速度和温度分布曲线能很好地与纯MD方法结果吻合。
关键词:
耦合
,
非周期性边界
,
边界力
,
有限容积法
,
分子动力学模拟方法