三维气固两相混合层湍流拟序结构的直接数值模拟.pdfVIP

中国工程热物理学会 论文编号 燃烧学学术会议 樊建人·+ 郑友取 岑可法 764(0)，E—nail：Fanjr妇ail．hz 浙江大学能源系，杭州310027．TEL：(0571)7951 zj．cn ■厦 本文对三维气固两相混合层湍流拟序结构进行了直接数值模拟。气相流场应用 拟谱方法对卜-S方程组进行直接求解，通过模拟时间模式的混台层流动，分析流场失 稳后涡的卷起、配对、合并及撕裂过程，研究三维混合层湍流拟序结构的演变特征； 计算颓粒场时，采用LaF“gian方法，针对颗粒不同的stokes数，模拟了颗粒场的 瞬态分布．并分析流场三维大涡结构对颗粒分布的影响。 关■蜀『两相湍流拟序结构； t 骱守嬲剃堋粉布 1引言 研究气固两相湍流拟序结构以及颗粒在三维湍流混台层的分布机理对能源、化工、 环境等工业生产过程有着重要的指导作用。近年来，由于实验技术的改进及计算机的 发展，很多学者已通过实验、分析和直接数值模拟肄手段对纯气相混合层中涡系的详 细结构及其演化作了大量的研究，并获得了可喜的进展口卜HI。但是国内外对气固两相湍 流混合层的拟序结构的研究却很少有报道。 湍流的直接数值模拟是在不引进任何湍流模式的情况下，数值求解完整的N—S方程 组．解出不同尺度包含脉动运动在内的瞬时速度场。它的结果可用来检验和发展各种模 型，井为探索研究湍流的新方法提供依据。谱方法由于它具有精度高、收敛快、准确 的空间微分罅特点，是目前直接数值模拟采用最多的数值方法。但由于计算机资源的 局限性，目前它只能计算中等h扣olds数和简单几何边界的流动口l。尽管纯气相混合 层流动的数值模拟、实验研究和理论分析取得较一致的结论，但要对气固两相湍流混 合层的拟序结构的研究，特别是气固两相湍流混合层中湍流涡团的卷起、配对、合并 和撕裂过程，以及颗粒分布机理作出一个确切的动力学描述尚有困难，还有待于今后 国内外学者进一步研究，本文是作者研究三维气固两相湍流混合层拟序结构所取得的 初步研究结果。 。图采重大基础研究发展规划资助项目(编号：05) ”联系大 257 2数学模型 2．1基本方程 2．1．I藏量控■方程 圈l描绘了平面混台层流场示意图，混台层上下主流速度分别为仉，仉(们(纠。 ● 流体为不可压缩牛顿流体。无量纲形式的N．S方程组如下： ● =， 一 可 口+ 可 U 型打 上胎 = 9 0 U 式中：F=U×疗．口为总压，涡量甜=V×U，计算时约定U。=U，一U：和氏 为特征尺度(吼为2倍的混合层初始动量厚度)，鼬ynolds数胎=乩瓯／y。 2．1．2焉粒茜动方程 ． 本文在模拟颗粒场时，颗粒质量流率小于1，可以假定颗粒分散相为稀相，因此不 考虑颗粒之间的碰撞影响．无量纲的颗粒运动方程描述如下： 警-，／即∽ ㈩ 其中．y是颗粒的速度，U是颗粒所在点的流场速度。r是Stokes阻力系数的修正因 子． ，=1+o．15Re：”，Re：=Iu一矿ld，／V．stoks数定义为： ＆=旦孵，颗粒的速度和位置可以通过积分c。，式求得。 ● 2．2漉蟠韧始条件 初始条件由以下几部分组成，