新型湍流模型特性地研究.pdfVIP

中国工程热物理学会 热机气动热力学 2 004学术会议 编号：042085 新型湍流模型特性的研究 雷俊勇’，祁海鹰，由长福，徐旭常 (热能动力工程与热科学重点实验室，清华大学热能工程系，北京100084) 010Emaih 摘要：针对最近由高歌、熊焰发展的G．Y湍流模型，选择零压力梯度平板转捩作为对此jj!jj流模型 的考察算例，得到了与实验数据吻合很好的边界层速度分布．讨论了数值计算过程中G·Y模型所存 在的某些特征。通过简化分析，给出了G-Y模型中漂移流与湍流强度闻的定量联系，并以此计算， Bypass转捩，计算结果与实验吻合较好。最后，给出了发展G-Y模型所要进行的^丁作和可能的途径。 关键词：G．Y湍流模型，平板转捩，侧偏统计平均，数值计算 一、前言 众所周知，由于存在多尺度，宽频谱等特征，湍流是至今尚未解决的经典物理问 题之一。尽管可以使用DNS方法直接求解流动控制方程，但计算量难以令人承受。目 前一r程麻用可接受的仍然是以PANS为基础的各种湍流模型。由丁．对湍流机制了解的 欠缺以及数学描述的困难，以RANS为基础的各种湍流经验模型难以适用于各种场合。 因此，对现有湍流模型不断改进以及发展新模型仍然是当前揣流研究的重点之一。 动信息，并采用了符合物理判断的假设——正交各向异性，为湍流模型发展开辟了新 思路。针对一种新模型，从多方面详细考察。从而得到其数值计算特征和模型对各种 问题的适应性，无论是对模型的发展还是工程应用，无疑都具有重要意义。 在考察湍流模型的各种算例中，平板转捩具有特殊的意义。对于边界层转捩的数 值计算，由于在计算域中同时存在层流、过渡和湍流区，建立在高雷诺数假设上的各 种湍流模型在此处不能很好适应。尽管如此，出于计算成本等因素，目前求解转捩常 用方法仍然是使用修正的LRN(LowReynolds 修正的K-8模型计算了圆管边界层转捩；文献[4-6]也分别使用不同的修正项，对平板 转捩进行了数值预测。近期的研究者17,”对使用不同的LRN k-E模型计算平板转捩进行 了广泛考察，仔细研究了不同修正方法对计算结果的影响。得到结论：在定量方面， 各种模型在转捩计算上都不能给出精确的答案：转捩计算对k、E的初边值以及计算选 1作者简介：雷俊勇，男，博士研究生，从事气固两相湍流模型的研究 基金项目：离家重点基础研究专项经费资助(NoG1999022200) 184 择的开始点非常敏感：选用不同的经验参数极火地影响计算结果。文献[9】【10】分别使 用大涡、DNS计算了平板转捩问题。计算结果显示，除了DNS，其它计算或多或少的 需要添加各种经验成分。 综上所述，由于包含过多的经验系数，目前的湍流模型对转捩的模拟不能令人满 意，而不含经验系数的G．Y模型却有可能很好的计算转捩。因此，我们选择了转捩作 为考察G．Y湍流模型的算例。 二、G．Y模型： 不同于通常雷诺时间平均方法，G．Y湍流模型理论出发点在于“侧偏统计平均”。 即对湍流随机变量(脉动速度、压力等)进行统计的样本空间不是针对所有的可能值， 而是在某个特定范围内对随机变量进行统计平均，以此得到随机变量一阶统计量。经 过系列假设和模化，最终推导出一组无经验系数的描述湍流运动的方程组…J。G．Y模 型最终需要求解的方程组如下： 平均流连续方程： V·百=0 (1) 平均流动量方程 弦 跖 (2) ，一p ●一p 漂移流连续方程 和一 (3) 漂移流动量方程