fene珠-簧链分子模型的分子动力学模拟.docxVIP

fene珠-簧链分子模型的分子动力学模拟 0 基于分子白砂糖动力学的流变学研究模型 高内部结构对变形运动的复杂动力学响应具有非线性滤波性质，导致奇怪的流动现象。由于加工流量诱导的液体结构在成品中被“分解”，以决定成品的质量。为了获得聚合物分子的微观运动信息，可以建立高聚物的微观结构和形状，建立流变学研究模型，并使用分子模型研究高聚物系统的流变性，阐明分子的微观结构、形状和宏观流变性之间的内在联系。这不仅可以获得速度场、压力场、磁体场和其他宏观量，还可以获得流程中分子的微观信息。 近年来,随着计算机水平的大幅度提高,分子Brown动力学的计算机模拟作为一种数值研究手段在流变学领域得以广泛应用.笔者根据聚合物加工流动过程的特点,提出基于分子Brown动力学原理的FENE(Finitely Extensible Nonlinear Elastic)珠-簧链分子模型,对聚合物在加工流动过程中的一些关键问题,如流动过程分子链的位形、流动应力进行模拟. 1 理论基础 1.1 弹簧长度的选择 笔者采用的分子模型是FENE珠-簧链结构如图1所示,珠-簧链包含M个珠子和Ns=M-1个具有一定长度的弹簧.珠子代表溶液中相互作用的质点,不计重量的弹簧代表链中的约束.其中:rn(n=1,2,…,M)表示珠子的位形向量;Qj(j=1,2,…,Ns)表示珠-珠之间的连接向量;rc表示珠-簧链的质量中心;RG表示珠-簧链的轨迹半径. 1.2 弹簧模型的建立 动态链的控制方程可以由作用在珠子上的外部力的总和获得 FΗiHi+FEiEi+FBiBi=0,i=1,2,3,…,M(1) 式中:下标i代表第i个珠子;FΗiHi,FEiEi,FBiBi分别代表液体拖曳力、弹簧约束力和布朗力. 液体拖曳力是作用在珠子上的Stokes拖曳,不考虑液体之间的相互作用和流体惯性,液体拖曳力由以下公式给出 FΗi=-ζ(˙ri-u∞i)FHi=?ζ(r˙i?u∞i)(2) 式中:˙γγ˙i是小珠的速度;u∞i是第i个珠子的速度;ζ是液体的阻力系数. 弹簧的有效约束力为 式中:Fsisi为第i根弹簧所受到的弹力,它依赖我们所采用的模型.笔者采用FENE珠-簧链分子模型,则弹簧的有效约束力表示为 FFEΝEi=ΗsQi1-Q2i/Q20(4) 式中:Qi=ri+1-ri是第i个弹簧的连接向量;ri是第i个小珠的位形矢量;Hs是弹簧常数;Q0是每根弹簧的最大长度. 布朗力用来表征珠子和溶液分子碰撞频度,满足以下条件 式中:kB是波尔兹曼常数;T是绝对温度;〈…〉代表系综平均;δ是单位张量;δij是delta函数. 将式(2)、(3)、(5)、(6)代入式(1)并整理后,可以得到以下控制珠子位置演化的方程 dri=[du∞i+FEiζ]dt+√2kBΤζdwi,(i=1,2,?Μ)(7) 若溶剂是牛顿流体,则溶液的偏应力张量可用Kramers型系综平均表示为 τ=τs+τp?τs=-ηs?˙γ(8) 式中:˙γ是变形率张量(κ+κT);ηs是溶剂黏度常数;τs,τp分别为溶剂分子和聚合物分子对聚合物溶液应力的贡献. 聚合物应力张量τp可用Kramers型系综平均表示为 τp=-n〈QFE〉+nMκBTδ, (9) 式中:n单位体积内聚合物分子数目;δ是单位张量. 2 弹簧所受力的fsit 对控制珠子位置演化的方程,我们采用半隐式预测误差法.对式(7)进行无量纲处理,可以得到 dri=[κ?ri+14FEi]dt+1√2dwi(10) 式中:距离、时间分别除√kΤ/Ηs、ζ/4Hs进行无量纲化;dwi=√Δtrni?κ是速度梯度张量的转置;在无量纲情况下简单剪切流κ为: κ(t)=˙γ(t) 对(10)式,如果分别取两个相邻方程相减,并经过无量纲处理可以得到: {dQ1=[(κ?Q1)+14(Fs2-2Fs1)]δt+√12(dw2-dw1)dQi=[(κ?Qi)+14(Fsi-1-2Fsi+Fsi+1)]δt+√12(dwi+1-dwi),1iΝsdQΝs=[(κ?QΝs)+14(FsΝs-1-2FsΝs)]δt+√12(dwΝs+1-dwΝs)(11) 无量纲情况下FENE珠-簧链中弹簧所受力Fsi可以表示为: Fsi=Qi1-Q2i/b(12) 式中:b=HsQ20/κBT是无量纲伸展参数,Q0是每根弹簧的最大长度. 在每一个时间步长内,计算位移连接向量Q分三步: (1)用显示方法计算初值 Q*i=Qni+[14(Fn,si-1-2Fn,si+Fn,si+1)]δt+(κ?Qni)δt+√12(dwi+1-dwi)(13) (2)用隐式方法离散中间变量 ˉQi+12(δt)ˉFsi=Qni+12(κ?Qni+κ?Q*i)δt+14(ˉFsi-1+Fs,ni+1)δt+√