空间发展湍流气粒两相平面混合层的大涡模拟与统计研究.pdfVIP

中国工程热物理学会 多相流 第十届年会 编号:016023 空间发展湍流气粒两相平面混合层 的大涡模拟与统计* 李勇 柳朝晖 郑楚光 华中科技大学煤嫌烧国家重点实验室 武汉430074 周向阳’ UniversityofCaliforniaatRiverside,USA Phone:027E-mail:liuzh@263.net 摘要:本文对空间发展的湍流气固两相平面混合层流动进行了大涡模拟研究。其中气相亚网格尺度 (SGS)使用结构函数模拟，气相控制方程组采用SHAPLE方法求解，固体颗粒运动用拉格朗日方法 计算。计算结果正确反映了流体祸结构的卷起、合并和破碎过程，以及小尺寸颗粒在涡边缘 (低涡度 区〕的局部富集现象.对直径分别为42m,72nm和135um分别进行了棋拟，并将统计结果和实脸测 量结果(Hishida,Ando和Maeda,1992)比较，表明两者的平均速度吻合很好，但颖粒数密度和脉动速度 存在较明显的差异，因此有必要对亚网格应力和颖粒之间的祸合作用展开深入研究. 关健词:大祸模拟，气粒两相流，平面混合层 1.引言: 气粒两相流的动力学在诸多工程应用和大气研究中十分重要。对于输运流中颗粒漂 移的实验与理论研究已经进行了很长时间[2.3.41。平面混合层是自由剪切层湍流研究的经 典算例，无论是数值还是实验方法都己经充分讨论。深入研究此类的湍流剪切流中颗粒 漂移机理，对于许多生产过程尤为重要，诸如锅炉燃烧、给粉器和喷雾燃烧室。 对于工程气粒两相流动的模式方法的研究，在近二十年有长足的进步，已分别从拉 格朗日和双流体两种途径建立了较系统的描述湍流两相流动的单向拐合、双向精合乃至 四向撰合作用的方法，但是由于两相湍流的内在复杂性和实验研究的困难，对描述两相 湍流相互作用的很多核心问题，如两相速度脉动的关联特性、颗粒对气体湍流的附加耗 散、颗粒的局部富集对两相湍流的影响等，都缺乏统一的或确切的认识。寻求更高级的 数学模型或实验手段对这些基本现象进行研究，是改进和发展两相湍流模式方法的迫切 需要。 在湍流的数值研究中，大涡模拟是一种+分有效的方法。可以以比直接模拟低得多 的计算代价得到大量的流动结构详细信息和工程上关心的平均参数.其基本思路是直接 ’本研究受国家3$点基础研究专项经费(G1999022207)f1国家自豁科学签金资助 第一作者:李勇，男，刀岁.硕士 联系人:柳栩晖，男，28岁.博士 156 计算大尺度结构流动，而对于数值分辨率以下的小尺度运动，需采用不同复杂程度的模 型模拟，如涡粘性膜性、结构函数模型和混合模型等，其中应用最广泛的是由 Smagorinsky于1963年提出的非线性涡粘性模型151。从Dearciddorflbl的第一篇大涡模拟 论文发表至今，己过30年，在这段期间，大涡模拟己相当成熟:基础理论进一步深入， 新的方法得到发展与验证，使用了更有效地数值方法。随着计算机运算速度和存储能力 的提高，大祸模拟可以应用于更多的流动问题，诸如包含化学反应的气相与气粒两相的 流动，压缩与不可压缩的流动等。 已有较多的文献对时间发展的平面混合层进行LES研究h-91，最近，Wang等还报道 了对时间发展的两相平面混合层数值模拟结果(101。本文将致力于空间发展的湍流气粒两 相平面混合层的大涡模拟。值得注意的是，空间发展模式的混合层同时间发展模式的不 同，后者假定在沿流线方向有周期性的边界条件191，但实际过程中混合层两侧卷吸的卷 吸是非对称的，必须考虑空间充分发展的情况。 本文的数值模拟将给出空间发展湍流平面混合层流场的瞬时空间结构信息和颗粒 的分布，并将大涡模拟的统计结果与实验数据[[1.10〕进行比较。 2.气粒两相流的大涡模型 在大祸模拟中大的能量输运尺度直接计算，只对较小的亚网格尺度量进行