耗散粒子动力学的简单介绍和应用前景..docx

6

耗散粒子动力学的简洁介绍和应用前景

mg0424112 徐源

一．耗散粒子动力学的进展

耗散粒子动力学〔dissipativeparticledynamics〕是一种欣的计算机模拟和描述流体的方法，是对分子动力学〔MD〕和LGA模拟的继承和进展。分子动力学描述的精度较高，但是计算的代价较高，到目前为止，只能用来成功地处理一些简洁的流体。另外，分子动力学应用条件比较苛刻，只能处理两维问题。

LGA〔lattice-gasautomata〕是1986年Frisch,Hasslacher,Pomeau,Wolfram提出的描述流体行为的模拟方法，随后Rothman和Keller进展这个方法，使它能描绘不能互融的流体行为。但是LGA模拟有一个缺乏：模拟中，LGA引入一个重要的概念：格子(lattice),格子的存在导致伽利略不变性的消逝，因此，在描述压缩流体和多相流体时，误差较大。

耗散粒子动力学，集合了以上两种方法的优点：排解了虚拟格子的概念，从而避开了LGA方法的精度的麻烦，另一方面，保存了分别时间步骤的概念，简化了模型，加快了计算的过程。更重要的是，与前面两者相比较，耗散粒子动力学更简洁和准确的模拟了三维状态下流体的行为，因此，具有更重要的意义。

在耗散粒子动力学中，根本颗粒是“格子”，它表示流体材料的一个小区域，相当于MD模拟中我们所生疏的原子和分子。假设全部小于一个格子半径的自由

度被调整出去只保存格子间粗粒状的相互作用。在格子之间存在三种力，使得每个格子对保持格子数和线性动量都守恒：简谐守恒相互作用〔保守力〕，表示运动的格子之间的粘滞阻力〔耗散力〕和为保持不集中对系统的能量输入〔随机力〕。全部这些力都是短程力并具有一个固定的截止半径。通过选择适当这些力的大小，可得到一个相应于gibbs-cano系统的稳定态。对于格子运动方程积分可以产生一条通过系统相空间轨迹线，由它可以计算得到全部的热力学可观测量

〔如密度场，序参量，相关函数，拉伸张量等〕。

耗散粒子动力学通过记录粒子的位置和相应时刻来描绘粒子的行为。粒子在碰撞的过程中，动量和质量是守恒的，能量是不守恒的。和LGA模拟一样，模型在一个时间步的演化是分为两小步不断交替重复来进展的：〔1〕一个微小的碰撞时间〔impulsestep〕〔2〕一个传播时间△t(propagationstep).在碰撞步骤里，每个粒子〔Pi〕通过相互作用转变了方向和速度，然后进入传播时间，在这个步骤里，速度已经转变，每个粒子都沿着肯定的速度滑行，和其他粒子没有相互作用。

碰撞时间步骤：

?p???e?

i

j?i

ij ij

传播时间步骤：

?q??t(p??p)

i m i i

方程中：m是每个粒子的质量，q

i

代表粒子i的位置，e

?ij

?

是粒子j指向粒

子i的单位矢量，为了简化计算，以后我们假设每个粒子的质m量=1.?

if

转移的粒子i的动量，在HoogerbruggeandKoelman的原始模型中，?

代表从j

有如下

形式：

? 3 ? r??

if

if

if

p p ??

? ?1?

???

?w( ?

)e?

???r

2n??

r?? if

i j if?if?r?r

if ? c c?

ij c

其中：

?0???if_r?r

?if c

?

r?q?q 是两个粒子i和j之间的距离，n=N/V是系统的密度：系统

if i j

体积为V，由N个粒子。?

ij

是随机指数，使得整个数据代表碰撞的整体和随

机效果，影响了流体的压强，温度。(p?

i

p)e

?j ij

?

反映了流体的黏度。r是

c

临界半径，只有在这个范围内粒子间才会有相互作用，这是的粒子间的力定义为短

程力，而且使得极大的简化了计算步骤。

连续的时间步骤〔碰撞时间和传播时间〕使得系统满足微平衡的要求，并保

证了Gibbsian平衡的存在9。但是，Espanol和Warren数据分析说明特别当?t?1

时，Hoogerbrugge和Koelman的DPD模型不能保证正确的平衡分布。但是，

?t?1时，数据模拟和试验吻合较好。

为了使DPD模拟的平衡态满足经典平衡态的要求。Espanol和warren对Hoogerbrugge和Koelman的模型作了两点改进。

削减时间步骤长度〔timesteplength〕

2 在?

ij

的耗散项表达式中，引入一个修正因子2〔1?r r〕