耗散粒子动力计算方法简介与应用.doc

耗散粒子动力计算方法简介及应用，附有计算机仿真出的相图 1楼 大 中 小 发表于 2006-11-6 02:10? 只看该作者 [转帖]耗散粒子动力计算方法简介及应用，附有计算机仿真出的相图 请注意 : 如果想要引用此帖的任何内容，请一定要在文章里加着原作者与发表的期刊刊名以及此帖后面所注明的reference，如有违者必须自行负上有关智慧财产权法的全部责任及法律责任。作者 : 邱佑宗出处 : 工业材料杂志　213 期　93 年9 月内容 : 前　言? ?? ???传统之分子动力计算(Molecular Dynamics; MD)(1)，长度尺度约在数奈米(Nanometer)至数十奈米左右，而每一时间步骤(Time Step)约在0.5~1 飞秒(Fentosecond)。以现今计算机之计算速度与记忆容量，合理的仿真步数，约在百万(106)至千万(108)，故其所计算问题之模拟时间，约在奈秒(Nanosecond)左右。然而，许多材料制程中个别事件发生所涉及之空间尺度及所需之时间，常不仅止于此。以胶体、悬浮微粒于流体中之运动等现象为例，单一尺寸颗粒即在10 奈米至1 微米之间，故以分子动力计算方法，搭配现今之计算机计算，尚无法仿真类似过程；在另一方面来说，连体力学(Continuum Mechanics)亦不适用于此类问题。故发展介于奈米尺度与微米尺度间之模拟方法，已成为近十年来的热门课题之一。? ?? ???介观尺度计算方法的相关研究者众(2)，比较常用的有直接模拟蒙地卡罗(Direct Simulation Monte Carlo; DSMC)、平滑粒子水动力学(Smoothed Particle Hydrodynamics;SPH)、网格波兹曼法(Lattice BoltzmanMethod; LBM)， 以及耗散粒子动力学(Dissipative Particle Dynamics; DPD)等。上述方法之源起与适用领域各有所不同，直接模拟蒙地卡罗适于描述稀薄气体(Rare Gas)；平滑粒子水动力学多用于天文学之星系模拟；网格波兹曼法则特别适于计算具复杂边界形状之流体行为。至于耗散粒子动力学，则多被用来模拟流体相分离(Phase Separation)、界面活性剂(Surfactant)等高分子于水中运动等现象。? ?? ???耗散粒子动力计算(3)为一介观尺度之水动力学计算方法，适合模拟奈米至微米尺度之复杂流体，例如：胶体、悬浮微粒于流体中之运动等现象。此外，耗散粒子动力计算属粗粒化方法(Coarse Grained Method)，其以单一粒子代表一个至数百甚至数千的原子，用以模拟不发生化学反应的复杂流体系统极为合适。本文即简介耗散粒子动力计算方法之原理，以及相分离、两性分子自聚集(Self Aggregation)等应用实例。原　理? ?? ? 耗散粒子动力计算(Dissipative Particle Dynamics; DPD)属粗粒化方法，其以单一粒子代表一个至数百甚至数千的原子，以牛顿力学加上统计力学的观念，计算系统中所有粒子的运动轨迹；再以此信息计算物质的巨观性质，如表面张力、压力，或观察分子排列的情形。? ?? ???耗散粒子动力计算之基本架构(4)与分子动力计算类似，均采用粒子作为计算单元。惟分子动力计算中之粒子为原子，其尺寸、重量及作用力场(Force Field)都有清楚的定义。然而，耗散粒子动力计算中，一个以豆子(Bead)为称呼的粒子，可代表数个至数千个原子。粒子与粒子间之作用力场，除与分子动力计算类似，具有保守力(Conservative Force)，另加了漫散力(Random Force)及与其对应之耗散力(Dissipative Force)。首先，每个粒子的运动状态由古典动力学之牛顿运动方程式决定： 上式中ri 表第i 个粒子之位置向量；vi 则表其速度向量。而粒子i 所受到之作用力向量fi 则为------------DPD公式 其中Fij C为两粒子间之保守力向量，多数文献采用容许粒子相互穿透的短程线性斥力 这样的作用力，除了本身数学型式简单、计算迅速外，尚有其它的优点，例如系统容易达到平衡状态、截断半径小、所需计算的粒子对少等。同时，当模拟粒子为一团(Lump)原子而非单一原子时，这样的假设也不致于失真。此外，改变aij 可调整两粒子间之斥力大小。除保守力外，Fij R与Fij D则为耗散粒子动力学独特之漫散力及与其对应之耗散力。漫散力为： 其中σ 控制漫散力的强度；θij 为一均值为1 的随机数函数、wR(rij)则依两粒子间距调整漫散力之大小。耗散力则为： ? ?? ?其中γ 控制耗散力的强度；wD(rij)则依两粒子间距调