耗散粒子动力学.pptVIP

* 排斥力参数，全原子模拟和试验方法可获取F-H参数 * FH理论适用于研究液体-液体、液体-固体的高分子混合系统。依据FH理论,两成分的混合系统,其平均自由能可表示为： χ为A与B间的作用参数。即： 将FH理论与DPD方法比较,可推得参数α与χ间的关系为： 随机力与耗散力参数 耗散力和随机力的值通过涨落-耗散定理相关联, 两者中只能任取一个, 我们把这个参数取为耗散力。Groot和间发现高斯噪声与均匀噪声对体系的运行几乎没有影响, 因此选用简单的均匀噪声, 当噪声幅度大于σ=8时, 积分运行结果不稳定, 当σ=3时, 在温度区间kBT=1-10内, 弛豫过程快而合理。因此, 在选择耗散力参数时, 所选的温度和参数不能使噪声比σ = 3 大, 否则模拟结果不可靠。 由式（6）可得当kBT=1时，γ=4.5。 2.4 DPD方法改进 2.4.1耗散力的改进 为更好地改进 DPD 模拟系统的动力学行为，我们将耗散力权函数修正为更通用形式： 耗散黏度可以表示为关于s的函数如（下）式所示 2.4.1耗散力的改进 耗散力权函数 上图显示了权函数在不同 s 值时的变化趋势。当 s= 2时权函数是与传统 DPD 模拟中的权函数表达式是一致的.权函数 wD(r)以及它的斜率在r/rc =1时都是连续的，然而当s1时权函数在r/rc =1c时仍连续，但是它的斜率不是连续的.并且s越小跳动越大.但权函数斜率的不连续对系统的影响不明显.故我们只要求权函数自身的连续，并允许在较小范围内权函数斜率的不连续性. 2.4.1耗散力的改进 S对粘度的影响 2.4.1耗散力的改进 三种DPD的动力特性 从表中可以看出把耗散力权函数的指数从 2.0 变到 0.5 就能够明显增加系统 Schmidt 数目和耗散黏度。 表中也表明了耗散系数γ和截断半径rc对系统动力学性质的决定作用,Sc与γ2和rc8成比例。然而增加Sc对系统温度会造成大的波动。假使系统温度得到有效控制，可以通过增加Sc来改进系统性质，但是实施这种方法有很多困难。很显然最有效的增加 Sc的方法是增加rc，又考虑到对耗散力和随机力的计算要求是它们与rc3成比例。模拟复杂系统时计算耗时又是一个考虑的重要因素。在目前研究中，我们结合修正的权函数和适当地增大截断半径来达到合理的计算耗时。 2.4.2 保守力权函数的改进  传统的DPD方法采用一个简单线性形式的保守力权函数ωC(r)=1-r,该式所描述的是一个纯排斥性质的软作用力。 采用此种纯排斥性的保守力权函数的DPD模型理论上只能模拟类似气体的特性，因为纯排斥力导致DPD粒子相互排斥离散开来，占据整个计算区域（或容器体积）。 虽然传统的DPD模型被广泛应用于多组分多相系统，却不能够模拟单组分多相（气-液-固）系统，因为粒子间只有排斥力，而没有吸引力将粒子结合在一起。因此传统的DPD模型不能模拟带有自由表面的流体流动，中空带气泡的液体，液滴动力学特性等等含气-液甚至气-液-固共存的多相系统。 2.4.2 保守力权函数的改进 因此传统的DPD模型不能模拟带有自由表面的流体流动，中空带气泡的液体，液滴动力学特性等等含气-液甚至气-液-固共存的多相系统。考虑到气-液-固共存及相互转化的物理本质,在保守力势函数中引入吸引力是必需的。引入在DPD粒子非常趋近时表现为排斥形式，而当DPD粒子间距离趋远时表现为吸引形式的保守力势函数能够描述气-液共存的流体系统 2.4.2 保守力权函数的改进 引入在DPD粒子非常趋近时表现为排斥形式，而当DPD粒子间距离趋远时表现为吸引形式的保守力势函数能够描述气-液共存的流体系统。 构造这种近距排斥远程吸引的保守力权函数的方式较多。一种可能的方式是把不同作用强度和截距的光滑粒子动力学的光滑函数叠加，从而组合出满足特定要求的保守力权函数。光滑粒子动力学中最常用的是三次样条光滑函数 2.4.2 保守力权函数的改进 采用光滑粒子动力学中最常用的三次样条光滑函数 该式对应的系数在一维、二维、三维系统中分别为3/2h、10/7πh2、1/πh3 * 2.4.2 保守力权函数的改进 三次样条光滑函数所定义的是非负的单调递减函数，在原点及截距光滑过渡。采用不同的作用强度系数及截距，能够组合构造出如下部分区域为正，部分区域为负的函数 U（r）以作为DPD模型中粒子间相互作用的保守力势函数。 式中W1（r）和W2（r）是非正则化的三次样条光滑函数，rc1和rc2，以及A和B分别是对应于W1（r）和W2（r）的截距及作用强度系数。 2.4.2 保守力权函数的改进 不同的rc1和rc2以及A和