耗散粒子动力学方法

1、应用化工1252dpd发展历史发展历史dpd基本原理基本原理dpd基本操作基本操作dpd应用举例应用举例一二三四提出理论基础发展1992年，hoogerbrugge和koelman在结合分子动力学和气体格子法的优点的基础上， 找到了一个能够在介观的时间与空间尺度上模拟复杂流体的方法。dpd方法1995年espanol和warren提出dpd中耗散力和随机力中的权函数必须满足涨落耗散定理。1997年groot和warren通过自由能将高分子系统的flory-huggins理论与dpd方法相联系1999年，coveney和fabritiis将dpd粒子质量和大小定义成可变。2001年pagonab

2、arraga提出一种多体dpd方法。2003年groot将静电作用引入dpd模拟。2004年groot由界面张力获取相互作用参数。2005eric e. keaveny将md与dpd模拟相关联。一、基本概念：一、基本概念：珠子（珠子（beadbead）：）：截断半径（截断半径（r rc c, cut-off radius, cut-off radius）简化单位简化单位多个原子或分子的集合体。珠子的作用范围用单个珠子的质量、截断半径作为单位量纲。如：crrrimmm作用力：二、方法体系二、方法体系水珠子油珠子1、运动方程以及作用力运动方程：iidrvdtiiidvmfdtjirijdijcij

3、iffff作用力：jirijdijcijiffffpairwise dissipation forcespairwise random forcesthree forcethree force作用力：排斥力、耗散力、随机力说明：1995年espanol和warren提出dpd中耗散力和随机力中的权函数必须满足涨落耗散定理。1997年groot和warren通过自由能将高分子系统的flory-huggins理论与dpd方法相联系2004年groot由界面张力获取相互作用参数。(1)cijijijijfarr()dijdijijijfrvr 1/2rijrijijftr repulsive con

4、servative forces修正的velocity-verlet算法：二、方法体系二、方法体系2、数值积分方法 212,12iiiiiiiiiiiiiiir ttr ttvtf tv ttv ttf tf ttfr ttv ttv ttv tt f tf tt 说明：用粒子当前的位置、速度和力来计算下一个时刻的位置和速度, 然后再用新的位置和速度计算新的力, 进而修正速度, 每运行一步, 力就更新一次。温度稳定效果较好65. 0排斥力：二、方法体系二、方法体系3、排斥力、耗散力和随机力参数的选择jirijdijcijiffffcijf在dpd中它代表物理系统中原子和分子之间复杂作用力的全部

5、。可以通过dpd流体的状态方程来确定排斥参量，在materials studio中amorphous cell、blends、discover等模块配合使用可计算出排斥力参数；通过试验测量界面张力也可求得排斥力参数。随机力：耗散力：dijfrijf耗散力和随机力的值通过涨落耗散定理相关联, 两者中只能任取一个, 也就是说只设置耗散力参数。thwo参数参数wohtwoht排斥力、耗散力和随机力参数 排斥力参数：ijadpd流体flory-huggins理论自由能排斥力参数的计算模拟体系的相互作用与flory-huggins理论的映射相联系3253.50ijijax ： 1、最直接的方法是用flo

6、ry-huggins模型关联高分子系统的相平衡或其他热力学性质的实验结果而得到 参数； 2、利用分子力学方法得到粗粒化粒子系统的混合能而获得 参数。ijxijxijxdpd方法有助于解决在配方化学、药物设计、高分子科学、化学工程涉及到的复杂问题。它为复杂流体的研究提供了强大的理论工具。目前已经成功应用于复杂流体的研究，如：油/水/表面活性剂体系、聚合物、胶体溶液的化学形态、微观形貌、相分离和流体流变特性等方面的研究。1 1、对复杂流体到、对复杂流体到达平衡态的模拟达平衡态的模拟复杂流体中的聚合物或表面活性剂等的聚集行为研究多相流中聚合物的微相分离和形貌变化其它方面2 2、对非平衡态下、对非平衡

7、态下流动系统的模拟流动系统的模拟液滴的变形、破碎微通道内的流动其它方面非平衡态平衡态非平衡态平衡态1 1、对复杂流体到达平衡态的模拟、对复杂流体到达平衡态的模拟3%10%20%30%35%40%55%正十四烷基硫酸钠sts聚集形态表面活性剂体积百分数表面活性剂体积百分数1 1、对复杂流体到达平衡态的模拟、对复杂流体到达平衡态的模拟c=45%c=85%hydrogenated surfactant fluorinated surfactant氟表面活性剂聚集形态表面活性剂复配机理表面活性剂复配机理图一 界面张力随混合表面活性剂浓度变化图两种表面活性剂复配使界面孔隙减少，可以更加有效地降低界面张力。图二 表面活性剂在界面聚集形态图(a) 单一表面活性剂， (b) 两种表面活性剂复配(a)(b)1 1、对复杂流体到达平衡态的模拟、对复杂流体到达平衡态的模拟1 1、对复杂流体到达平衡态的模拟、对复杂流体到达平衡态的模拟2 2、对非平衡态下流动系统的模拟、对非平衡态下流动系统的模拟2 2、对非平衡态下流动系统的模拟、对非平衡态下流动系统的模拟利用dpd模拟囊泡破裂过程2 2、对非平衡态下流动、对非平衡态下流动系统的模拟系统的模拟利用dpd研究了毛细血管中红细胞的叠积机理。3 3、其它方面应用、