热对流条件下双颗粒沉降的直接数值模拟.PDF

第卷第期 年月 ５９ ３ ２０１０ ３ 物　 理　 学　 报 Ｖｏｌ．５９，Ｎｏ．３，Ｍａｒｃｈ，２０１０ １０００３２９０／ ２０１０／ ５９ ０３ ／１　８７７０７ （） ＡＣＴＡ ＰＨＹＳＩＣＡ ＳＩＮＩＣＡ ２０１０ Ｃｈｉｎ．Ｐｈｙｓ．Ｓｏｃ．  热对流条件下双颗粒沉降的直接数值模拟 刘汉涛常建忠安康苏铁熊 　 　 　 　 （中北大学能源环境工程与计算流体力学实验室，太原 ） 　 ０３００５１ （２００８ 年７ 月２６ 日收到；２００９ 年７ 月１０ 日收到修改稿） 　 　 在任意拉格朗日欧拉（ＡＬＥ）算法模拟等温惰性颗粒两相流的基础上，增加对能量方程的联立求解，对热对流 条件下双颗粒的沉降进行了直接数值模拟．结果表明：双颗粒在等温流体、热流体和冷流体中运动形式是不同的， 造成不同的主要因素是热对流引起了颗粒不同的动态尾迹，在热流体中颗粒尾部形成了涡的脱落，在冷流体中颗 粒尾部形成了羽流． 关键词：直接数值模拟，双颗粒沉降，热对流，两相流 ＰＡＣＣ：４７１０，４７５５Ｋ，９５３０Ｌ ［］６ ． ＡＬＥ 和非线性的代数方程组 另一种是推广的 标 １ 引 言 准伽略金有限元法．把颗粒相受力方程和扭矩方程 Ｇａｌｅｒｋｉｎ 与连续相的动量方程一起进行有限元的 方 颗粒沉降在自然界和工业设备中广泛存在，对 法推导，这样颗粒相和连续相相互作用的力和扭矩 其流动现象、机理和过程的研究无论是在学术理论 就不必专门加以计算． 数值计算的网格通过求解 还是在工程应用上都有重要意义由于计算上的困． Ｌａｐｌａｃｅ方程得到网格移动速度，来满足由于颗粒相 难及研究方法的限制，目前对运动颗粒两相流动的 的运动而导致的整个求解流场的移动．分散相位置 研究往往只局限于无热对流的惰性颗粒，国内外对 的更新将由其速度决定时间步长由分散相的速度． 于有热对流的两相流动机理还未能很好的理解和 和加速度来自动调整．如局部网格发生畸变将重新 ［１—３］ 划分网格．方程的非线性部分由牛顿迭代求解，线 掌握 ． 实验方法中，激光粒子速度场仪ＰＩＶ （ｐａｒｔｉｃｌｅ ｉｍａｇｅ ｖｅｌｏｃｉｍｅｔｒｙ），ＰＴＶ （ｐａｒｔｉｃｌｅ ｔｒａｃｋｉｎｇ 性部分由ＧＭＲＥＳ算法来求解． ｖｅｌｏｃｉｍｅｔｒｙ ），粒子动态分析仪ＰＤＡ （ｐａｒｔｉｃｌｅ ｄｙｎａｍｉｃｓ ａｎａｌｙｚｅｒ）虽能获取流场结构，但无法获得 ２ 数值方法 运动颗粒间的受力，不能测量热对流情况下颗粒与 ［４，５］ 本文在应用任意拉格朗日欧拉（ ）方法模 流体间热量质量的传递 ．数值分析方面，直接数  ＡＬＥ ［，］ ７ ８ 值模拟由于将计算网格缩小到颗粒尺寸以下进行 拟