气固两相流动中颗粒扩散的转捩现象.pdfVIP

讲学屯苏 第51卷第23期2006年12,El 论文 气固两相流动中颗粒扩散的转捩现象 罗坤 樊建人 岑可法 (浙江大学能源清洁利用国家重点实验室，杭州310027．E-mail：zjulk@zju．edu．cn) 摘要不引入任何湍流模型，采用有限容积方法对三维气固两相湍流射流进行了直接数值模拟，着重 考察湍流结构的转捩行为对不同Stokes数颗粒扩散的影响．为了得到高精度的结果，同时降低计算量 和存储量，对气相流体控制方程组的求解采用分步投影算法，对时间积分采用低存储、三阶精度的 Runge—Kutta积分格式；对颗粒控制方程的求解在拉格朗日框架中进行．流场的统计结果-9相关的实验 数据吻合良好，证实了数值算法的可靠性．对颗粒扩散的研究发现，展向涡结构对颗粒扩散的影响比较 突出．而在流场的演化过程中，观察到颗粒扩散的新行为，即中、小Stokes数的颗粒在流场中的分布出 现了“转捩”现象． 关键词 直接数值模拟气固两相流 拟序结构 颗粒扩散 转捩 气固两相流是一个只有30多年发展历史的新兴 现象，它表示流动状态由层流向湍流的过渡过程，以 学科，但它在许多工程技术领域中有着十分重要的 及由大尺度涡结构向小尺度涡结构的过渡过程．在 应用．仅以能源电力工业为例，从煤粉的制备、干燥、 转捩过程中，流场中不仅有大尺度涡结构，而且同时 输送、风粉混合和分离过程，到煤的流化、煤粉燃烧、 存在多种尺度的涡结构．这些不同尺度的涡结构与 煤粉再燃以及污染物的排放等一系列过程均是气固 颗粒间具有不同的耦合作用．为了深入理解多尺度 两相流的重要研究课题．近年来，我国北方一些城市 湍流流动中颗粒的扩散行为，本文采用直接数值模 和地区频繁爆发的“沙尘暴”现象，也是典型的气固 拟方法对气固两相湍流转捩射流进行了研究，着重 两相流问题． 考察湍流大尺度结构转捩为小尺度过程中不同 气固两相流动中离散固体颗粒在气相流场中的 Stokes数颗粒的扩散特性． 运动，被称为颗粒的扩散或弥散．颗粒的扩散往往成 1数学描述 为一个影响系统效率和稳定性的重要因素．预测并 1．1流动结构 控制气固两相流动中颗粒的扩散对于工程设备的优 化设计和高效应用具有重要的意义． 本文选取三维气固两相平面湍流射流进行研究， 其典型的结构示意图见文献[9】中的图1．为了观察流 Crowe教授领导的研究组经过一系列的数值模 动的转捩过程，基于流体人流速度砜和喷口宽度d的 拟和实验研究n巧】，率先应用颗粒的Stokes数，即颗 Iv=5000．整个计算 粒的空气动力学时间常数与流场大涡特征时间常数 初始流动雷诺数设置为Re=Und 的比值，来描述自由剪切流中不同颗粒的扩散特征． 整个人流喷口射人流场，而颗粒则仅从中心区边长为 Stokes数较小的颗粒容易被流体卷吸，在流场中的浓 d的正方形区域喷入计算区域．初始时刻，在射流人 度分布较均匀；Stokes数为1的量级的中等颗粒受大 涡结构旋转离心作用的影响最为突出，许多颗粒被甩 口高度为d的区域内给定如下速度剖面的剪切层： ，， ，， ／ 、 到大涡结构的外围，扩散率最高，浓度分布最不均匀； “：塑+垃tanhI上l， 而Stokes、数较大的颗粒由于自身惯性较大，扩散程