耗散粒子动力学对颗粒沉降问题的研究.pdfVIP

中国工程热物理学会 流体机械 学术学会论文 编q-．067007 耗散粒子动力学对颗粒沉降问题的研究 孔轶华，张楚华，席光 (西安交通大学能源与动力工程学院，710049，西安) 联系电话：029-maih xiguang@mail．xjtu．edu．cn 摘要：给出了耗散流体粒子运动的控制方程组、边界条件、数值方法等。使用耗散粒子动力学方法编 程计算了颗粒在重力作用下的沉降运动，观察到颗粒的质量和所受重力对颗粒运动轨迹的影响，且颗 粒质量越小．所受重力越小，颗粒运动所表现出的随机性越强烈。从而验证了所采用的数学模型、计 算方法在流动数值模拟中的可行性及潜在优势。 关键词：介观模型耗散粒子动力学颗粒沉降随机耗散理论 1 引言 Particle Dynamics)是一种介观模 流体流动的耗散粒子动力学【11(DPD，Dissipative 型方法。这种模型认为流体由大量相互作用的假想流体粒子(DPD粒子)组成，而每个 DPD粒子又包含大量的微观流体分子或原子。 颗粒在通道内的运动过程及其沉降规律是一个流固两相流的机理性问题，广泛应用 于生物医学、环境科学、能源动力等领域。本文使用DPD算法模拟颗粒在重力作用下 在由两垂直壁面限制的静止流体中的沉降现象。计算结果对于颗粒悬浮物在通道内的运 动过程、沉降规律及沉淀分布等实际问题具有指导意义。． 2 理论模型 DPD方法通过求解一组DPD粒子群的运动位移、速度、加速度得到流体运动规律。 对于任一个DPD粒子f，其位移及动量方程为： 要=谚，要=∑元+丘 u’(1) 云邓，i2毛厶¨c 亏和E分别为粒子f的位置和速度矢量，厶粒子／对粒子f作用力，C为外力，其中，兀 由三部分组成： f4=F；七F；+F； ‘2) 露为守恒力，是作用在中心线方向的软斥力： 霞F={口l『(1一勺)勺，勺名， (3) 9 【o，勺≥％ 其中嘞为粒子f和歹间最大斥力，弓=亏一乃，～=I‘l，0=乃／R1，当粒子间距离 大于k，粒子问作用力忽略不计。耗散力和随机力分别为： ：印=一7wD(白)(毛·％)弓 (4) (5) 劈=OWR(～)巳名=洲尺(勺)氲△t-1／2，,、驴 其中y和仃为耗散力和随机力强度，wD(％)和W月(勺)为权重函数，吃=谚一矿，， 只，为满足高斯分布的随机函数，幺；为随机数，均值为零，方差为l，对于每个时间步 长、每对粒子相互独立。 耗散力和随机力需满足平衡条件： 13 w％)=h)『’y=杀 (6) 权重函数与粒子间距离有关，’本文采用Groot和Warrenl21推荐形式： (7) w。(，．)=hR(，．)】2=【。