两方程紊流模型在三维高速紊流流场数值模拟中的应用牛.PDFVIP

—— 篓：!：三旦竺璺生簦堕堡垄兰垒坚笙苎亘窒：!!坚兰!旦!!旦=!!旦 两方程紊流模型在三维高速紊流流场 数值模拟中的应用牛 马晓永李国君李亮李军 (西安交通大学能源与动力工程学院，陕西西安710049) 丫骶缩性引起的压缩性耗散和压力膨胀的影响可压缩k-CO模型和可压缩k-g模型，结合任意曲线坐标系下的三维雷 诺时均可压缩N-S方程组，对三维喷管内的跨一超音速流动和跨音速三维透平叶栅紊流流场进行了数值计算，并对不 Ⅲ紊流模型的计算结果进行了对比。结果表明在模拟三维高速紊流流场的一些复杂流动现象中，k-g模型要优于标准 的k—s模型和标准的k—go模型：可压缩修正后的两方程模型更加适台于跨音速紊流流场的数值模拟：本文所建立的 数值求解方法对三维高速紊流流动具有较高的模拟精度．能较好地模拟流场中流动参数的分布规律和分布特点。 关键词：紊流模型；可压缩性：三维高速紊流；数值模拟 0引言 紊流的数值模拟一直是计算流体动力学(CFD)的热点研究领域，目前面对工程应用的紊流数值 模拟研究仍必须采用引入紊流模型的方法。工程中应用紊流模型大多在以下两个方面存在问题：一是 逆压梯度下，有无分离流动的预测及低雷诺数情况下的模拟。另一个是可压缩流动，特别是高速流动 的预测偏差较大。两方程模型是目前预测精度较理想的模型，因为其具有以下特点：理论性强，使用 的经验公式少，而引入的物理概念多；Boussinesq假设有效避免使用6分量的雷诺应力模型，减少 了待求方程数目；在壁面上两方程有精确的边界条件，易于处理；在高速内流计算中表现出来的良好 性能。另外，对于目前透平叶栅内部紊流流场，大多是高速的跨音速流动，但计算所采用的紊流模型 仍延Hj在不可压缩流动基础上发展而成的紊流模型，而很少考虑到可压缩性的影响，这可能是各种紊 流模型对可压缩流动，特别是高速流动的预测偏差大的原因。 鉴于上述分析，本文选取了3种应用较多的两方程紊流模型和它们的可压缩修正模型：标准的 K-W K—G 型(comp model)。运用将可压缩紊流模型与任意曲线坐 model)和可压缩k-g模型(comp 标系下的三维雷诺时均可压缩N—S方程组相结合的计算方法，对两个算例：三维喷管内的跨一超音速 流动和跨音速三维透平叶栅紊流流场进行数值计算，比较了5种模型之间的优劣并分析其产生的差异 原因。 1 计算方法和紊流模型 采用通过Faver质量加权平均及无量纲化得到的任意曲线坐标系中矩阵形式的时均可压缩N-S方 法处理壁面。在计算过程中采用当地时间步长加速收敛。需要指出的是，采用两方程紊流模型时，需 求解的方程组有七个方程，若采用统一的高分辨率TVD格式，则特征矩阵为七阶矩阵，迭代求解时， ’基金项目：国家自然科学基金重点项目(No 640 两方程紊流模型在三维高速紊流流场数值模拟中的应用 这样大的矩阵的刚性较大，使得计算不易收敛。为此本文将基本控制方程组和紊流模型方程分开求解， 采用不同的差分格式和求解方法。这样避免了上述问题，也有利于数值计算的快速收敛。对紊流模型 方程求解时，先转化为任意曲线坐标系中无量纲化形式，差分离散左边的对流项采用二阶精度的QUICK 格式离散，右端的扩散项及生成 项采用中心差分格式，对于接理而得的线性方程组，采用交替方向隐式迭代法即ADI方法求解。 对于收敛准则，检查相邻两次迭代的流动参数相对变化量是否满足精度要求，即： △≯=霉野I眵嚣：一≯0，tJ／庐0，tf ，。 max陋≯ls ≯为流动参数：密度p和三个方向的速度“、v、W，占为给定的正小数。在实际的计算过程中，占 的选取需一定的经验。如果占取的过大，数值计算并未达到真实的收敛解，可能是假解；￡过小，会 使计算时间过长，而且最后阶段的收敛速度很慢，参数变化己不太明显。通常，以流动参数的最大相 对变化量下降2～3个量级确定的取值。 本文计算采用的标准的k．占模型和标准的k．∞模型方程这里不再给出，详情可查阅相关文献资 模型在物面附近数值