耗散粒子动力学中施加无滑移边界条件的新方法.pdfVIP

中国工程热物理学会 学科类别：传热传质学 学术会议论文 编号： 123455 耗散粒子动力学中施加无滑移边界条件的 新方法一 肖自锋，曹知红，易红亮，谈和平 （哈尔滨工业大学能源学院，哈尔滨150001） (Tel: 0451 Email: yihongliang@hit.edu.cn) 摘 要 耗散粒子动力学(DPD) 中固液边界条件常常会出现速度滑移现象，本文提出了一种新的处理方法 来实现无滑移边界条件。通过选取合适的壁面粒子与流体粒子之间的保守力系数可以控制近壁面 区粒子数，从而可以使密度波动降低至极小值，实现无滑移边界条件。通过该方法模拟泊肃叶流 动求得的流体性质与通过使用 Lees-Edwards 边界条件模拟剪切流动所求得的流体性质一致，证明 了该方法简单可靠。 关键词 耗散粒子动力学；无滑移边界条件；保守力系数 0前言 耗散粒子动力学(Dissipative Particle Dynamics, DPD)是由荷兰学者 Hoogerbrugge 和 Koelman 于 1992 年提出的一种介观的数值模拟方法[1]，可以看成是粗粒化的分子动力学方 法，每一个 DPD 粒子包含一团原子或者分子，可以使用更大的粒子尺寸和更长的时间步长， 因此比传统的分子动力学方法对计算机性能的要求更低，能应用到微观尺度和宏观尺度的问 题。另外，DPD 方法是一种基于粒子的方法，对于复杂流体的模拟有着得天独厚的优势。 近年来，耗散粒子动力学方法日益受到重视，逐渐应用到了各类复杂的流体流动，如相分离、 蛋白质等大分子悬浮、表面活性剂、胶体输运、聚合物溶液、生物薄膜、以及纳米悬浮液的 [2] 对流换热现象 。 DPD 方法的固液边界条件一般是通过冻结规则分布的 DPD 粒子代表固体壁面来实现。 由于DPD 方法中软势力的采用，并不能阻碍部分粒子逃逸模拟空间，因而须引入映像条件 把粒子重新引入模拟区域。为实现无滑移边界条件，一般采用反弹映像[3] ，即规定重新引入 的粒子速度方向和大小均与碰撞前相反。 DPD 方法中固液边界条件存在两个问题：第一，由于壁面粒子施加在耗散粒子上的保 守力随着粒子与壁面间距离而波动[4] ，这必然引起壁面附近流体密度的扭曲；第二，流体粒 子在壁面粒子所施加保守力的作用下远离壁面，那么具有降低流体粒子速度作用的耗散力也 随之降低，这将导致固体壁面处速度滑移。从原理上讲，只要粒子数分布均匀，无滑移边界 条件即可实现。 Duc Duong-Hong 等人[5]于 2004 年提出了调节两层固体壁面间的粒子距离的方法，用以 降低近壁面区流体粒子所受力的波动，从而实现无滑移边界条件。 Pivkin 等人[4]于 2005 年通过模拟得到了壁面粒子施加于流体粒子的保守力随两类粒子 之间保守力系数 abi 变化的函数，然后让其等于流体粒子本身所受其他流体粒子的保守力， 即可得到波动最小的保守力系数，并产生无滑移边界条件。但是此种方法所得的abi 过小， 一 国家自然科学基金 (No.51176040)，新世纪优秀人才支持计划 (NCET‐09‐0067) 壁面处粒子数密度过大，仍会导致较大的波动。 本文在 Pivkin 的基础之上直接模拟得到不同abi 下的粒子数密度波动，然后取其最小值 所对应 abi 即可得到无滑移边界条件。通过该方法模拟泊肃叶流动求得的流体性质与通过使 用 Lees-Edwards 边界条件模拟剪切流动所求得的流体性质一致，证明了该方法简单可靠。 1DPD基本思想 DPD 方法认为流体由大量相互作用的 DPD 粒子组成。每个 DPD 粒子又包含了大量的 微观流体分子或原子，分子或原子的无规则热运动通过 DPD 粒子的随机和耗散特性描述。 对于简单的流体，每个 粒子受到三种力作用：（）相互排斥的保守力；（）耗散力， D