均匀各向同性湍流中气粒两相流湍流变动附加项的模型与模拟.pdf

中国工程热物理学会 多相流 学术会议论文 编号：066105 均匀各向同性湍流中气粒两相流湍流变动 附加项的模型与模拟 何森棋柳朝晖林钊 华中科技大学煤燃烧国家重点实验室 武汉430074 Tel：027-8410Email：heqi4321@163．com 摘要： 气粒两相流中湍流变动的原理较为复杂，颗粒的出现会使气相湍流的双方程模型中出现反映流 体一颗粒湍流相互作用的的附加项，从而使流体湍动能增加或减少。本文分别从视密度加权平均和概率 密度输运方程两种途径出发，基于量纲分析，分别构造了两相关联耗散率与拟湍动能或两相关联湍动 能的显式代数关系式，井将其同时用于两相关联湍动能和气体耗散率方程的封闭。使用新构造的模型， 对文献中直接模拟的衰减均匀各向同性湍流(E19hobashi2003)进行了模拟，检验和分析了新构造的模 型对两相关联拟动能和气体耗散率的模拟效果。 关键词：气粒两相流湍流变动两相关联耗散率湍流模型 1．引言 湍流气粒两相湍流中包含很多复杂的过程：湍流对颗粒的扩散，颗粒对湍流的反向 耦合作用，颗粒一颗粒碰撞，颗粒一壁面碰撞等。要构建全面合理的气粒两相湍流模型， 对上面各方面的机理和子模型的研究都是极其重要的，其中最基本的又是两个问题：颗 粒的扩散与湍流变动。颗粒扩散是采用单向耦合的方法，在不考虑颗粒对气体作用的情 况下研究气体湍流对颗粒扩散的作用；而湍流变动是采用双向耦合的方法，同时考虑了 颗粒对流体的反向作用。 在湍流变动中，颗粒的加入会使K—e方程中出现由颗粒所产生的附加项。K—e方 程中由颗粒所产生的附加项与湍流变动中两相间的相互作用有关。对于附加项的处理， r 最简单的方法是认为耗散率的附加项与湍能的附加项成正比，即s口=孚k，通常五 ‘￡ 采用流体的特征时间尺度。为了能更全面的考虑两相流湍流变动中颗粒与流体之间的相 互作用对时间尺度的影响，本文采用颗粒所见流体微团的拉氏积分时间尺度气，来代替 556 流体特征时间尺度五，对模型进行改进，将颗粒扩散的“三个效应”的影响引入到模型 中来。 本文采用Zaichik的模型封闭湍流动能方程中的附加项，同时对于各向同性湍流中 计算颗粒所见流体的拉氏积分时间尺度时的两相相对平均滑移的速度I¨I，根据DNS的 关联耗散率50的改进模型进行数值验证· 2．K--e方程中附加项的封闭 2．1K--￡方程中的附加项 湍流变动的结果会给流体的K和e方程带来关于颗粒作用的附加项。 气粒两相流中K和e的通用输运方程为“： 婴，尝：户一D+T (1) 一 Dt‘Dt 其中P为产生项。D为耗散项，T为输运项。在均匀各向同性湍流中，不存在各种 平均梯度，因此产生项和输运项为零，只剩耗散项，该方程可简化为： 车：一s+战 (2) 窘一e：考¨缸 ㈤ 其中，址，A6为颗粒作用的附加项： 从一警(k—k)，船一型rr,(％) (4) 其中相关项匕=三霭为