颗粒沉降运动的虚拟区域法直接数值模拟.pdf

中国工程热物理学会 多相流 第干—一届年会 编号：056089 颗粒沉降运动的虚拟区域法直接数值模拟 陈丽华邵雪明樊建人 浙江大学机械与能源工程学院，杭州，310027 摘要本文对二维情况下在方槽内775个颗粒的沉降运动采用基于四边形嗍格的分布式拉格朗日乘子 ／虚拟区域方法_i{f：行了直接数值模拟。展示了颗粒流沉降过程中流体的速度和涡量分布、流场雎力分布， 颖粒的运动轨迹和运动速度等参数。结果表明了颗粒流在沉降过程中在流场内形成人小不～的旋流区， 颖粒的屠涡对附近其他颗粒的运动有较大的影响。本文的计算说明拉格朗日乘子／虚拟区域方法对于模 拟存在有很多颗粒的流动是比较适合的。 关键词：拉格朗日乘予／虚拟区域方法，颗粒沉降，旋流区，尾涡 1引言 I劁体颗粒在重力作用F的沉降运动在自然界和工业应用中非常普遍。不同形状、不 同尺度的颗粒在不同雷诺数下具有不同的沉降运动特性，比如胶体颗粒微球依靠重力沉 降组装形成类似于原子或分子晶体的有序点阵，厌氧颗粒污泥依靠良好的沉降和自胤定 化技术从而有效地将微生物保留在厌氧反应器中以实现反应器高效稳定的运行等。这些 新技术的发展对两相颗粒悬浮流的沉降运动的机理研究提出更多更高的要求。 分布式拉格朗日乘子，虚拟区域方法是近年来出现的赢接数值模拟颗粒流的新方法， 特别是包含大量颗粒的流动。不同于以往的势流11．21和Stokes流[3．41等一些近似的方法， 拉格朗日乘子／虚拟区域方法的基本思想是口J：通过假设流体也充满颗粒占据的空间，这 样把一个儿何复杂区域(可能时随时间变化的)内的问题扩展到一个相对大的，但儿何 简单的区域(虚拟区域)中，控制流体流动的N--S方程和控制颗粒运动的刚体运动方 程通过无滑移边界条件和作．}fj在颗粒上的力和力矩耦合在一起进行求解，而刚体运动的 限制通过一个拉格朗日乘子来施加。这一方法的优点是：①作用在颗粒上的力和力矩不 再需要显式计算；②由于计算区域不再是随时间变化的，可以应用同定的结构化网格， 而不再象基于移动的非结构化网格的直接模拟方法那样，每一个时间步都需要重画网格 和投影，从而人大节约了计算时间：③由丁二扩展的计算区域一般都是几何简单的，这样 可以采用规则的网格和快速求解算法。 本文在二维情况下基于四边形网格对这一数值方法进行了实现，并对在二维方槽中 775颗圆形颗粒的沉降过程进行了模拟，对其沉降动力特性进行了讨论。 2数学模型 2，1控制方程 为简化起见，假设计算区域是矩形的，如图l所示。Q表示整个的区域．包括颗 为计算区域的边界。同时假设流体密度P，、 474 颗粒密度办和重力加速度g为常数。 颗粒和流体运动的控制方程分别如下 流体运动方程： 矶【争(u-V)u～尉V·盯 ‘ inQ＼P(t)， (1) L V·u=0inQ＼P(t)， (2) 图I围粒．流体两相流系统示意圈 盯=一pI-1-2qD[u](3) u=ur(，)Onr (4) U=U，+出，Xr』 oiloe,(f)，f-l，．，N， (5) ul，．o=uoinQ＼尸(o) (6) 颗粒运动方程： (7) M，盟：M，g+F』+一，d‘t 。。 ‘ ’’ ，，idfoj=z，