三维射流中颗粒碰撞的直接数值模拟.pdfVIP

中国工程热物理学会 多相流 学术会议论文 编号：076052 三维射流中颗粒碰撞的直接数值模拟 闫洁-，罗坤，。樊建人-，肖刚z l浙江大学能源清洁利用国家重点实验室，杭州310027 2东南大学洁净煤发电及燃烧技术教育部重点实验室。南京2l0096 联系电话：057l·879517“；E-mail：舀g逃@蓟』!四g：￡旦 摘要：为了进一步研究两相流动中颗粒的碰撞对其扩散特性的影响，对三维阿相射流进行了直接数值 模拟。其中对流场控制方程采用有限容积法和分步投影算法，对颗粒的罪踪采用拉格朗日方法，颗粒 问的碰撞采用硬球模型模拟。初步考察了sook韩数为10的中等颗粒的扩散情况，结果发现在流场中颗 粒局部浓度较高的区域碰撞频繁发生；由于射流初期太尺度涡结构的影响，颗粒的浓度分布最为不均， 因此碰撞次数在这一时期随时间呈线性增加，达到最大值后逐渐回落趋于平缓。此外，考虑颗粒间相 互碰撞作用以后，颗粒的分布更加均匀，颗粒沿横向和展向的扩散也都有所增加。 关键词：颗粒碰撞；硬球模型；直接数值模拟；气固两相三维射流 l引言 由于具有独特的流动特性，湍流射流一直受到广泛的关注。特别是气固两相混合射 流，更是具有广泛的工程应用背景，如能源领域内的气力输送过程、气固分离过程、燃 烧过程，化工领域内的气固混合反应过程以及环境工程领域内的污染物控制过程等。长 期以来，人们对单相和两相射流进行了大量的研究，取得了很多成果114J，其中数值模拟 研究大部分都是在忽略颗粒碰撞影响的基础上进行的。但是以往的研究表明，即使颗粒 的体积浓度只有4×104，颗粒碰撞的影响也不能被忽略睁7J；即使流动为稀疏流动，由于 中等颗粒的“选择聚集”效应，流场中某些局部区域的颗粒浓度也可以达到整个流场的25 倍pJ。在这些情况下，在颗粒高度聚集的地方，颗粒间的碰撞对颗粒运动和两相流动都 具有重要影响，需要进一步的研究。 目前，通过数值模拟来研究两相流动中颗粒间碰撞的基本模型有三种：1)硬球模 型p’1“，即认为颗粒之间的碰撞都是瞬时的，且在碰撞过程中颗粒本身不发生变彤，在 此基础上根据动量守恒原理得到颗粒碰撞后的运动状态。这种模型的突出特点是，只要 给定摩擦系数和恢复系数就能够求出碰撞后的颗粒速度。该模型比较适合低浓度颗粒流 的计算．当颗粒浓度过高时，由于碰撞频率太大，计算起来就比较困难。2)软球模型【l。“J， 也称离散单元法，即假设颗粒间的碰撞作用可以持续一段时间，将颗粒处理为有弹性的 球，而且允许同时有两个以上的颗粒发生碰撞。该模型对于两个以上颗粒同时碰撞机会 很少的低浓度两相流不太合适，但对碰撞频率较高的高浓度两相流，如稠密流化床的研 究却非常合适。3)直接模拟Montecano方法【l”。其特点是通过概率抽样的方法，研究 系统中的N个物理颗粒由n个样本颗粒所代替进行计算．这里每个样本颗粒代表性质相同 的一组物理颗粒，nN。这种处理有效节省了计算时间和计算机内存，但是也决定了它 ‘基垒顶目·田塞白蒜科学基盘(墒号，鲫50∞27)和东南丈学i胄净煤发电及崧烧技术教育部重点实验室开放研究萋盒 256 无法获得颗粒在碰撞过程的具体受力信息。 在．r述背景下，本文在前期研究的基础上对三维气固两相射流中的颗粒碰撞过程进 了射流区域内颗粒碰撞发生的空问分布、碰撞次数随时间的发展变化以及碰撞对统计特 性的影响。其中对流体运动的求解采用直接数值模拟方法，对离散颗粒的跟踪在拉格朗 目坐标系下进行，颗粒间的碰撞采用硬球模型。 2数值方法 2．1气相场的模拟 假定流体为不可压缩的牛顿流体，忽略体积力，则无量纲的控制方程可表示如卜-： 连续性方程：旦坠：o (1) az 动量方程：堕+业：一土望+上旦f盟+坐l (2) af