颗粒碰撞的拉氏直接模拟的研究.pdfVIP

中国工程热物理学会 多相流 学术会议论文 编号：036081 颗粒碰撞的拉氏直接模拟研究 陈胤密，柳朝晖，郑楚光 (华中科技大学煤燃烧国家重点实验室，武汉430074) Tel：0278410 Eraail：chenyinmi@163．com 摘要：本文研究了三维空间内颗粒碰撞的拉各朗日模拟方法．成功的构建了算法并通过租重颗粒流的 碰撞算例对算法进行验证。网格排除法的优化措施，使计算量降低50％左右，而不增大相对误差．在 计算参数At，‰的敏感性考察中发现，模拟同一颗粒群的运动，当△f血眦，相对误差随At增大而 增失．而在满足不遗漏碰撞的情况下，相对误差与Ⅳ。。无关．颗粒碰撞的模拟对数据类型精度要求较 高．此算法为进一步研究细微颗粒的凝并过程莫定了基础． 关键词：气固两相流，颗粒碰撞，拉各朗日 0引言 颗粒间的碰撞和凝并过程普遍存在于自然现象和生产流程中，如大气中雾、雨的形 成和燃烧过程中微细颗粒物的形成和长大等。颗粒间的平均碰撞率是定量描述和模拟这 些过程的重要参数。常规的实验手段难以对颗粒碰撞过程进行测量以取得颗粒碰撞率， 而数值实验则是广为采用的替代手段。它对气相进行直接模拟，对颗粒相采用拉各朗日 模拟直接跟踪颗粒的轨迹并根据颗粒轨迹交叠与否来判断颗粒碰撞。 模拟颗粒的相互碰撞，计算量与％7成正比。限于计算能力，早期的气固两相流研 究，将颗粒视作质点来分析，不考虑颗粒碰撞的作用。随着计算性能的增强，借鉴分子 动力学的方法文f{1提出了各向均匀同性湍流中具有实际体积的颗粒的相互碰撞的拉各 朗日模拟方法。在此基础上，为研究细微颗粒间的平均碰撞率的数学模型，本文就算法 优化及计算参数的敏感性进行探讨。 1数值模拟方法 L1颗粒运动轨迹的模拟 在拉各朗日模拟中颗粒的运动轨迹由各个时间步长内的直线位移相连而成。颗粒的 位移和动量控制方程为： 坐：矿 dt≥’rr，y1。 u’(1) ％鲁砩掣+喜巧 局，％，怖，审和L诊别是颗粒的中心位置，速度，质量，迟豫时间和颗粒附近流 体速度的近似值。凡’是除了阻力以外作用在颗粒上各种外力，如布朗力，Saffrnan力等 371 等。 对颗粒轨迹的计算常_目{]阶Adams插值法，它的表达式为 X，(f+△r)=X，(，)+【c10(f+At)一cy,(t)]At (2) 文中a=0，c’=一I。 1．2颗粒碰撞的模拟 文[1】提及两种颗粒碰撞的模拟方法：预测法和追溯法。追溯法通过时间步长出终 止时颗粒的速度vp(t+AO反向计算过去一个时间步长△，内已经发生的碰撞。与预测法相 比，具有计算量小和易于实现的优势。虽然可能遗漏多元碰撞，但数密度较小的颗粒群 在较短的时间内多元碰撞的次数很少，追溯法是较好的选择。追溯法包括两个方面的内 容，即／(1)颗粒碰撞的判定和(2)相邻颗粒的范围的划分。 1．2．1颗粒碰撞的判定 假设三维空间内有颗粒PI和P2作相对运动。如图1(a)所示，假定没有发生碰撞， 坐标系的x轴与K(f)平行，Y轴位于由‘(f)和H(f)两个矢量确定的平面上。其中： ‘(f)=Xpl(f)一一z(f) 如果PI和P2在卜血到f这段时间发生了碰撞，则下式成立： ‘p一△‘)=(4+吐)／20s△‘sAt (5) 其中At是发生碰撞的时刻与f时刻之间的间隔。吐和吐分别是I和2的直径。实际过c P P 程中，因为各种外力的影响，Ⅵ(f)随时间变化。但为了计算上的简便，不妨假设在At内 颗粒速度为常值％l(t)和％2(t)，如此可得到颗粒中心相对距离的另一计算式： ‘(f—At,)=G(t)一w．(t)Ato