湍流模型的发展及其研究现状.docx

湍流模型的发展及其研究现状 陈永辉 1 王 强 2 朴明波 3 （1 中国联合工程公司 浙江杭州 310022 2 哈尔滨工业大学 黑龙江哈尔滨 150001 3 辽宁工程技术大学 辽宁阜新 123000） 摘要 湍流是锅炉内燃烧中很常见的流动现象，其三维、随机特性很强。 通常采用时间平均的 N-S 方程来模拟湍流流动现象， 并引入湍流模型来解决 N-S 时均化过程中产生的方程封闭问题。 目前在工程计算中存在多种湍流模型，但大都有适用的范围，尚不 存在一个通用的湍流模型。 回顾湍流模型的发展，并对目前工程中常用的湍流模型做一个总结。 关键词 燃烧 湍流 湍流模型 雷诺平均 中图分类号：O357.5 文献标识码：A 文章编号：1672-9064(2009)02－0004－04 在锅炉内的气体流动情况十分复杂 ， 一 般 是 三 维 、 两 相 、 ρ 坠U +ρ 坠U （U U ＋u＇ u＇ ）＝ 坠P ＋ 坠 （2μS ） （2） 有化学反应的湍流流动 ，其内工质的各种物理 参 数 ， 如 速 度 、 坠t 坠xj j i j i 坠xi 坠xj ji压力、温度等都随着时间与空间发生随机变化。 由于湍流同时 具有扩散和耗散的特性， 前者意味着湍流可以使包括湍流本 身的所有流体信息在空间上传播， 充分利用湍流的扩散特性 具有很高的实用价值，在燃烧问题中，湍流扩散可以使各种物 体充分混合，提高燃烧效率；后者说明湍流会引起自身机械能 的损耗，因此从工程学的角度来看 ，湍流既有 优 点 ， 也 有 其 弊 端。 为了更加有效地研究燃烧问题，首先应该对湍流问题作深 入地研究 ji 目前在炉内的燃烧模拟中， 对湍流模型的选择考虑得不 多，大多是采用二方程 k-ε 模型或者更复杂的模型进行模拟， 而实际上目前在工程中存在多种湍流模型， 但尚不存在一个 通用的湍流模型。 本文的目的在于将其中常用的湍流模型做 一个总结，阐述其研究现状，并对其的发展提出自己的观点。 1 湍流的主要研究方法 原则上，N-S 方程组可以用于描述湍流运动的，但其强烈 的非线性使得用解析的办法精确描写湍流的三维时间相关的 全部细节变得尤为困难，甚至不可能实现；而且湍流流动的全 部细节对于工程实际来说，也是没有多少意义的，因为人们所 关心的只是其总效的、平均的影响。 这就决定了对湍流的研究 主要采用统计平均的方法。 目前在湍流的平均控制方程中，一 般 有 时 间 平 均 （ 雷 诺 平 均 ） 和质量加权平均 （Favre） 的 控 制 方 程，而后者主要运用于可压缩流动中 ，而在计 算 中 ， 大 多 数 仍 采用了雷诺平均的方法。 不可压的雷诺时均 N-S 方程[1]，其形式如下： 坠U =0 （1） 坠xi 其中 Ui，Uj，，P 为时均的速度与压力值，由于流体不可压，故密 度 ρ 认为是一个定值，u＇i，u＇j，p＇ 为相应的脉动值 ，Sji 为变形率 张量。 从中可以看出雷诺时均方程中引入了一个新的项 u＇ju＇i， 该项为脉动量的相关项，又称为雷诺应力项，其中包含了湍流 的细节信息。 一般时均项的值可以通过实验测得，但是脉动项 就难以得到，这就造成了方程不封闭。 引入湍流模型的目的就 是构造脉动项与时均项的关系式， 解决方程的封闭问题。 但 是，在模化过程中也必然会带来流场信息丢失，对计算结果产 生影响。 由于基于平均方程与湍流模型的研究方法只适用于模拟 小尺度涡的湍流运动，该方法不能够从根本上解决湍流问题。 为了使湍流求解更加准确，更能反映不同尺度的旋涡运动，研 究人员又发展了大涡模拟 LES 与直接数值模拟 DNS 方法。 但 是在工程实际中，由于 LES 与 DNS 对计算机硬件的要求相对 要高很多，对于炉内燃烧模拟来说 ，除了流动 问 题 ， 还 要 考 虑 到两相、化学反应过程，这样的计算量对现有的计算机来说仍 然难以接受， 兼之采用湍流模型加平均 N－S 方程的湍流模拟 能够给出工程实用的结果，故而目前在实际应用中，湍流模型 仍然扮演着一个很重要的角色。 2 湍流模型的类型 目前工程计算中常 用 的 是 2 类 湍 流 模 型〔2〕, 一 种 是 基 于 Boussinesq 的涡粘性假设 ，如零方程模型 ，一方程模 型 和 二 方 程模型； 另一类是引入二阶脉动项的控制方程而形成二阶矩 封闭模型，或称为雷诺应力模型。 目前在工程应用中，还是一 方程模型与二方程模型使用较多。 相对来说，它们对计算机 的要求较低，同时可以得到符合工程要求的计算结果。 随着研究的深入， 统计力学方法也愈来愈在湍流模型研 作者简介：陈永辉（1972~），男，浙江