均匀各向同性湍流中颗粒扩散研究.pdfVIP

中国工程热物理学会 多相流 2002年学术会议 编号：026062 均匀各向同性湍流中的颗粒扩散研究 徐志锋 柳朝晖 郑楚光 (华中科技大学煤燃烧国家重点实验室，湖北武汉430074) TEL：027--8410 com E·mail：xuzf66@sohu 摘要： 均匀各向同性湍流中的颗粒扩散存在三个重要的效应：轨道穿越效应、惯性效应和连续性 效应，如何在颗粒所见流体速度的Langevin方程中考虑这些效应．是两相湍流PDF输运 方程模型封闭的关键问题。本文提出一种新的Langevin方程的漂移系数模型，以综台考虑 三种效应的影响，进行了相应的理论分析，并对文献中经典的颗粒湍流扩散实验进行了模 拟．结果正确地厦映了三个效应的影响． 关键词： 均匀再向同性湍流、两相流、扩散、PDF模型 ’ 1引言 大量的实验研究及数值计算结果表明，颗粒在湍流扩散过程中存在三种重要效应的 影响，即“轨道穿越效应”…颗粒扩散因为流体一颗粒相对速度滑移而减少，“惯性效应” …无体积力作用时颗粒的扩散随惯性的增大而增大．“连续性效应”…垂直和平行流体 一颗粒相对滑移速度方向上颗粒扩散系数不同。模型的主要任务是适当表示这些效应。 目前大量的计算颗粒扩散的模型均是在拉氏坐标系下进行的，它们可以归为两火 颗粒被假没同一系列的涡团相互作用，相互作用时间取涡生存时间和颗粒穿越涡团时间 的最小值，即tlnj=min(t“ay,tJ。。)：在TCM模型中，我们在跟踪某一颗粒时，还要跟踪和 该颗粒相对应的流体颗粒，该流体颗粒与颗粒从同一位置出发．通过采用拉氏时间关联 (沿着流体一颗粒轨道)和欧拉空间关联(填补流体和颗粒空间位置上的差距)获得颗 粒所见的流体速度。如果两者之间的距离大于流体涡团的长度尺度，则换用～个新的流 体颗粒。传统的模型一般仅考虑了轨道穿越效应，最近，研究者们分别尝试在ElM模 型【”和TCM模型l”中综合考虑三种效应的影响。 PDF方程模型是90年代以来新发展的两相湍流模型方法，它是基于湍流场随机性和 概率统计描述，该方法处理滴流问题有其独特的优势，非常适合于模拟那些必须考虑湍 流流动、复杂化学反应机理及其相互耦合的湍流反应流问题。两相湍流PDF输运方程模 型封闭的关键问题是，如何在颗粒所见流体速度的Langevin方程中考虑上述三个效应。 本文将在两相湍流PDF输运模型的基础上从封闭颗粒所见的气体瞬时速度的Langevin 方程山发。推导出～种新的漂移系数模型，综合考虑三种效应的影响．晟后和经典实验 结果进行对比。 2 模型描述 2．1 Simonin模型 颗粒所见的流体瞬时速度Langevin方程形式如下 奉项目受用家霞火基础研究专项经费fGl999022207)和国家自然科学基金酬J 等≮5蜀一万l瓦ap+砜a(u瓦a,；E,抄3。飞)磬 。—————r————一———百—一 +G。(％一i。)+嘭2q 。————百——一 其中u。，为颗粒所遇到的流体微团的瞬时速度，式(1)中，右边第一项反映流体平均速度 场对流体微团的加速度，第二项反映颗粒和流体微团轨道差异和流体平均速度梯度的附 加作用项．第三项则用于考虑粘性脉动、压力脉动和颗粒相对流体运动等的影响，其中 二阶张量G。。为相空间漂移系数，B。为相空间扩散系数，u．为一高斯白噪声。 我们将颗粒所见的流体瞬时速度Langevin方程中的漂移系数G。，写为： G㈣|=一6口}t” (21 需要注意的是，这里，t“为颗粒遇到的气体微团的拉氏积分时间尺度，是反映沿着颗 粒轨道的流体微团速度脉动相关性的衰减速度Simonin[10“j曾建议用如下的关系式来考 虑颗粒的“轨道交叉效应”和“连续性效应”对颗粒