耗散粒子动力学及其应用的新进展唐元晖.pdfVIP

·8 · 　　　　　　　高 　　分 　　子 　　通 　　报 20 12 年 1 月 　 综 　述 耗散粒子动力学及其应用的新进展 唐元晖 , 何彦东 , 王晓琳 (北京清华大学化学工程系化学工程国家重点实验室 ,北京 　100084) ( ) 　　摘要 :耗散粒子动力学 di ssip ative p article dynamics , D PD ,是 1992 年提出的一种针对于复杂流体的介观 模拟方法 ,凭着相对于其它模拟方法的高计算效率以及能够模拟流体水力学行为的优势 , 目前已经成为研究复 杂流体动力学行为的主要模拟方法之一 , 同时也是联系微观结构与宏观实验的重要桥梁 ,在辅助研究复杂流体 性质方面起到了重要的作用 ,本文对这一模拟方法在近几年的发展及其应用做了全面的总结 。首先介绍了近 年来方法本身的改进和完善 ,接着详细总结了目前的应用研究热点 。作为一种计算机实验工具 ,这种方法近年 来在对复杂流体的研究上进展迅速 ,既能用于预测平衡过程的结构和性质 ,又能用于研究非平衡过程的动力学 因素 。 　　关键词 : 耗散粒子动力学 ;介观模拟 ;复杂流体 引言 随着计算机科学技术的不断发展 ,越来越多的工程领域都采用计算机技术来帮助了解那些宏观实验 手段不能解决的问题 。流体是化学和化工等学科都非常关注的一类物质 ,按照空间尺度和运动的时间尺 度来分 ,它们可以在微观 、介观 、宏观多个尺度表现不同的性质如图 1 所示 ,从而催生了不同尺度下的研 ( ) ( ) 究方法 。对于微观尺度 0 1～10nm ,分子动力学 molecular dynamic ,MD 等分子模拟方法从上个世纪 ( ) [ 1～4 ] 50 年代起 ,就被应用于研究简单流体的平衡结构和 1～10p s 过程内的动力学 ;在宏观尺度 ,计算流 体力学辅助实验手段可以帮助人们很好地解决各种工程问题 ;但对于如胶体 、聚合物等这样的复杂流体 , ( ) 只有在介观尺度 10～1000nm ,10～1000n s 下才能更好地考察其结构和性质 。目前相关的模拟方法主 要有蒙特卡罗[ 5 ] (Mont e Carlo ,M C) 、自洽平均场理论[ 6 ,7 ] ( selfcon si st ent fiel d t heory , SCM F T) 、动态密 度泛函理论[ 8 ] (dynamic den sit y f unctional t heory ,DD F T) 、时间依赖的金兹伯格朗道方程[9 ] ( Ginzbur g L andau equation) ,但是这些方法在描述聚合物链运动时都存在限制于格子的问题 ,并且不适用于模拟枝 [ 10 ] ( 状聚合物 , 所 以出现 了在 1992 年 由荷兰 的 Hoo gerbru gge 和 Koelma 在格子气体法 L at ticega s auto ma , L GA) 和 MD 的基础上提出的 , 随后 由 War ren [ 11 ] 等发展的一种新的复杂流体的介观模拟方 ( ) 法 ———耗散粒子动力学 di ss