微液滴动力学特性的耗散粒子动力学模拟.pdfVIP

第57卷第7期2008年7月 物 理学 报 100l()．3290／2008／57(07)，3954．08 ACTAPHYSICASINICA @2008chin．Phys．soc． 微液滴动力学特性的耗散粒子动力学模拟* 常建忠1’ 刘谋斌2’ 刘汉涛1’+ I)(中J匕大学机电1：程学院，太原030051) 2)(新加坡南洋理工大!学化学与牛物医学工程系，新加坡649798) (2007年lo月25日收到；2007年12月12日收到修改稿) 对传统的耗散粒子动力学方法进行了改进．改进的耗散粒子动力学方法采用了包含远程吸引力和近距排斥力 的保守力势函数，从而使得用耗散粒子动力学方法模拟多相流动成为ar能．应用改进的耗散粒子动力学方法，对微 尺度F液滴的形成及液滴在微霞力下的大幅度振荡变形进行了数值模拟．计算结果表明，改进的耗散粒子动力学 (DPD)方法能够有效地描述微尺度下液滴的动力学特性，对研究复杂流体多相流动有着重要的意义． 关键词：多相流，微液滴，耗散粒子动力学(DPD)方法，保守力势函数 PACC：0340G，4710，4755K 容积等网格方法对模拟微液滴动力学特性也存在着 1．引 言 不少困难【7舟J．一方面，基于N．S方程的有限差分法 和有限容积法对介观和微观尺度的流体行为未必适 微尺度下液滴动力学，作为复杂流体多相流动 用；另一方面，有限差分法和有限容积法必须配备合 的一个重要分支，对喷墨打印、高分子及蛋白质悬浮 适算法以确定多相流动中的自由表面／运动界面．典 传输分析、微米／纳米粒子的合成以及药物运输等不 型的自由表面／运动界面追踪及捕捉算法有：1)粒子 同的工程领域有着非常重要的意义¨’2’．微尺度下液 滴动力学解析或半解析形式的理论研究，一般局限 在极少数非常特殊的情形中．而实验研究往往耗资 不菲，而且为了得到具有统计意义的数据导致实验 强．例如粒子追踪法需要大量的计算时间，复杂的界 周期漫长．近年来，越来越多的研究者把兴趣集中在 面重构使得VOF很难从二维推广到三维；LS往往伴 微尺度或者介观尺度下液滴动力学的数值研究上． 随质量的增加或减少，并对复杂几何区域适应性 微尺度及介观尺度下多相流体流动的计算模型 较差‘1引． 必须能够处理诸如自由表面和移动界面等运动边 介于纳米与微米之间的介观尺度模拟方法，已 界，必须能够处理不同流体间较大的密度比及黏度 经逐渐成为最热门的课题之一．这些方法包括：1)光 比，必须能处理由微尺度所带来的各种微效应．这些 滑粒子动力学(smoothed panicle 微效应包括低雷诺数(Reynolds肌mber，惯性一黏性 BoltzH聊n 比)，中高克卢森数(Knudsennumber，粒子尺寸．流动 区域比)以及较大的面积一体积比b—J．传统的分子 simulationMonte 动力学(molecular dyn咖ics，MD)方法b’6J，通过追踪(djssipativepanicle 各个分子／原子的运动(位置及动量)对整个流体系 点不同，应用各异．SPH源于天文物理学以模拟星系 统的特性进行模拟．其时间及空间尺度大致为纳秒 演变，最近被逐步推广到传统连续力学所适用的领 和纳米级．在现在的计算条件下，很难对介观尺度以 域，甚至进一步应用于微尺度下非连续介质的运动 上的流动区域进行模拟．而传统的有限差分和有限