气固两相三维湍流射流拟序结构及颗粒扩散的大涡模拟.pdfVIP

中国工程热物理学会 多相流 2002年学术会议 编号：026073 气固两相三维湍流射流拟序结构及颗粒扩 散的大涡模拟+ 罗坤金晗辉李文春樊建人”岑可法 (浙江大学热能工程研究所杭州310027) 摘要本文建立了气固两相矩形射流流场内气体、颗粒两相运动的数学模型并在此基础上完成了数值 模拟。其中对气流相运动采用大涡模拟．对颗粒相运动采用Lagrangian方法模拟．模拟结果显示了气 周两相矩形射流中气相流场在横向、展向和流向的大涡拟序结构。着重分析了流向涡拟序结构随对间 发展的演变过程；同时给出了气流相三维拟序结构作用下微细颗粒在整个流场空间中的运动扩散特性， 为工程应用提供了有益的参考． 关键词大涡模拟：矩形射流；拟序结构：颗粒扩散 1引言 在燃烧工程中，特别是利用切圆燃烧的电站锅炉，直流燃烧器的横截面经常做成矩 形。它的形状虽然简单，但其流动特性却属于复杂的三维两相矩形射流湍流流动。对矩 形射流进行系统的实验研究始于上个世纪60年代中叶．Sforza等人…最早用皮托管对 矩形射流进行测量，得到了流场的整体流动特征。Krothapalli等人“1采用热线仪对宽高 比AR=16．7的矩形射流进行了更为系统的测量研究，比较全面地考察7流场的三个发展 区域情况和各个截面上平均流速、脉动速度以及流场宽度的发展情况，并对流场大尺度 拟序结构进行了初步的分析．Tsuthiya等人“”则对宽高比AR5的矩形射流进行了研究。 Quinn等人”1在实验后也认为矩形射流相比平面射流、圆射流等二维射流来说，可以增 强与周围环境流体的混合与卷吸． 而在数值模拟方面，对矩形射流的研究也已经做了大量的工作．Mcguirk等人“1采 用双方程湍流模型对矩形射流流场进行模拟，认为射流的入流初始条件对流场影响很 大．Grinstein等人”1采用涡丝法对矩形射流喷口区的流场进行了模拟．Miller等人“1 的射流流向涡结构进行了分析． 上述研究大多都针对纯气相射流．对于气固鹿相矩形射流来说，由于颗粒相的存在， 致使两相射流的流动特性变得更为复杂，目前由于理论上和实验技术上的困难，在这方 面的研究尚不多见．而研究气固两相矩形滴流射流中大涡拟序结构的运动规律及其对颗 粒扩散方式的影响。具有重要的理论研究意义和工程应用价值． 基于此，本文运用大涡模拟的方法．研究较高雷诺数下矩形射流流场中不同截面上 展向涡量、横向涡量、流向涡量的分布特征，以及兰维拟序大涡结构作用下细微颗粒在 整个流场空间中的运动扩散特性，希望能为理论研究和工程应用提供有价值的信息。 2数学模型 2．1三维气相运动的亚网格应力(SGS)模型 }国家重点基础研究发展规划资助项目(编号： 本文米用Smagorinsky提出的亚网格尺厦(SGS，Sub—gridScale)模型“…。对各 流场瞬时变量在网格尺度上进行区域平均，可得到区域平均形式的三维N-S方程组： 连续性方程：丛趔：o 动量守恒方程： 鲁_=ux—j毒E可】_砉圈+毒ftr¨，[考+筹]] 在本文的计算中选取G--o．1。 2．2颗粒运动方程 本文在模拟颗粒运动时，颗粒质量流率小于1，可以假定颗粒分散相为稀相，忽略 颗粒间相互碰撞的影响。另外，还假定：流场中的所有的颗粒均为有规则的球体，并有 相同的密度pP和直径dP：颗粒的密度远大于流体的密度。根据牛顿第二运动定律，则 可将颗粒运动方程描述如下： m等=譬CoP：l(fi／-tip)l(uf-u．e)+m，g 进一步化简可得到： 鲁=如吲+g 式中：厅，代表颗粒所在位置的流场速度，i，代表颗粒的瞬时速度，颗粒的Stokes数 =pPd；／18,u，g为重力加速度，f为颗粒阻力系数的修iF_系数， 厂=l+o．15Reo：”7，飓为颗粒雷诺数，Re。=忙，一面。Id。／l，．4为颗粒的直径。 r o t o 。i ‘ 3数值方法和初值、边界条件 本文主要研究小宽高比矩形射流流场．矩形喷嘴的宽高比AR=2；喷口处射流的初始 来流速度Uo=lOm／s，并以台阶式分布