多相湍流流动基础和数值模拟西安讲学-3报告.ppt

西安理工大学讲学多相湍流流动基础和数值模拟(3）2005年7月7日-13日 多相湍流流动 模拟基础（1） 清华大学航天航空学院 工程力学系 周力行 教授 主要参考书 周力行，湍流气粒两相流动和燃烧的理论与数值模拟，科学出版社，北京，1994 周力行，多相湍流反应流体力学，国防工业出版社，北京，2002。 大 纲 基本守恒方程 液体/气体湍流模型— k-?模型 雷诺应力方程模型 代数应力模型 研究背景 流体湍流现象本身已经很复杂，多相湍流和湍流燃烧现象更复杂 连续相和离散相的脉动之间有强烈的相互作用 多相湍流流动和燃烧模型是CFD软件的理论基础 湍流模型，即周培源提出的封闭和求解雷诺方程的方法，目前仍然是解决工程问题的唯一有效方法 长期以来，湍流模型研究只针对单相流体本身 最近，双流体模型比颗粒拉氏轨道模型得到更快的发展 原来的双流体模型缺乏对颗粒脉动的确切认识 我们提出了“颗粒湍流”和“多相湍流模型” 的概念 基本守恒 方程组 层流/瞬态稀疏多相流动守恒方程组（1） 层流/瞬态稀疏多相流动守恒方程组（2） 针对湍流多相流动的Reynolds展开和取时平均 多相湍流反应流动的流体/气体方程组 雷诺平均后出现的二阶矩，即流体/气体雷诺应力、热流和物质流是未知的，要建立封闭模型 多相湍流反应流动的颗粒相方程组（只考虑阻力和重力） 雷诺平均后出现的二阶矩，即颗粒雷诺应力、热流和物质流是未知的，要建立封闭模型 液体/气体湍流模型 RANS模拟和湍流模型 湍流是流体力学中尚未完全解决的复杂理论问题 工程中早期常用的是取湍流粘性系数为常物性系数和射流积分法 周培源在50多年前建议求解Reynolds时均方程组及关联量输运方程的方法 30多年前Launder等人实现了湍流模型（模式）方法，即湍流的统观模拟，目前称为RANS 模拟 近代发展的DNS和LES，由于计算量很大，尚不能普遍用于工程问题 湍流模型的基本思想是根据物理概念用低阶关联量或平均流性质来模拟未知的高阶关联量，从而封闭平均量方程组或关联量方程组 雷诺应力方程的原始形式 模拟封闭的基本原则 分析每一项的物理意义 运用量纲分析法 在坐标变换中，模拟后的项与原来的项有相同的特性 对三阶关联量允许使用梯度模拟 各向同性耗散 各向异性向各向同性的回归 扩散项的梯度模拟-Daly-Harlow模型 压力应变项的IPCM模型（Isotropization of Production and Convection Model) k-ε双方程模型 混合长模型-零方程模型 湍流动能方程的简化-局部平衡：忽略对流和扩散：产生=耗散 其他的湍流模型 修正的 k-? 模型 -低Re 数k-? 模型-RNG- k-? 模型 -旋流和浮力流动的Richardson 修正的k-?模型 -可实现性 k-? 模型 和Reynolds 应力方程模型(DSM)相对应的代数应力模型（Algebraic Stress Model ，ASM) 多尺度k-? 和Reynolds 应力模型 用k-ε模型模拟无旋突扩回流流动 旋风分离器中强旋流动 突扩流动轴向速度（?实验，?DSM，---ASM，? k-?） 旋流突扩流动轴向速度（?实验，?DSM，---ASM，? k-?） 旋流突扩流动轴向速度(DSM) 旋流突扩流动轴向脉动速度(DSM) 旋流突扩流动轴向速度 环管中旋流流动 分离器中强旋流动（改进的 DSM ） 湍流模型应用的结论 不存在普遍而通用的湍流模型 对边界层、射流、管道、槽道和无分离的喷管流动等可以用混合长模型 对二维或三维无浮力，无旋或弱旋的突扩或钝体回流流动，可以使用各种k-ε模型 对强旋流动和浮力流动，用各种DSM或ASM模型 对应力方程的扩散项，压力应变项的模拟，耗散率方程的封闭等，针对流动的特点，还可以进行进一步研究 * * 瞬时量=平均量+脉动量. 产生项 Dij, ?ij, 和 ?ij 需要封闭模型 扩散项 压力-应变项 耗散项 各向同性的应力耗散模型 其他模型—GL 模型, SSG 和 FLT 非线性模型 Boussinesq 表达式 … k-ε ___ASM （a）k-ε（b）原 ASM（c）改进的 ASM *