三维方柱不可压缩绕流的大涡模拟计算 Large eddy simulation of 3-dimensional incompressible flow around a square cylinder.pdfVIP

第26卷第2期 空气动力学学报 V01．26．No．2 2008年06月 ACTAAERoDYNAMICASINICA Jun．，2008 文章编号：0258．1825{2008)02．0167．06 三维方柱不可压缩绕流的大涡模拟计算 邓小兵，张涵信，李 沁 (中国空气动力研究与发展中心，四川绵阳621000) 摘要：用非定常不可压缩流动求解方法对三维方柱绕流问题作了大涡模拟计算研究。求解程序采用虚拟压缩方 法和双时间步的时间推进方案，对流项的差分采用高阶精度的数值差分格式。通过与实验和文献结果比较分析， 验证了基于虚拟压缩方法的湍流大涡模拟计算方案对于复杂的三维钝体绕流计算是可行的。 关键词：湍流；大涡模拟；虚拟压缩方法；不可压缩流动 中图分类号：V211．3 文献标识码：A 为完整、可靠的实验结果。文献[7]汇总了一次以方 O 引 言 柱绕流和立方体绕流为研究对象的大涡模拟研讨会 非定常不可压缩流动的计算方法有投影法、SIM． 的结果，其中对方柱绕流问题给出了16个大涡模拟 PLE法、虚拟压缩方法等【1’2J。其中虚拟压缩方法虽计算结果的比较分析。 然早在20世纪60年代就已经提出【3J，但直到20世 湍流数值计算的主要方法有：雷诺平均方法 纪90年代双时间步方法提出后，才被用于非定常流 动的研究u，2I。虚拟压缩方法是类比于可压缩流动 (DNS)。雷诺平均方法需要用模型计算雷诺应力。 的求解方案，易于采用针对可压缩流动发展的各种高 这个方法只关心流场的时均量，对计算机的要求比较 精度差分格式，而湍流的大涡模拟研究通常需要采用 低，目前在工程实际中得到广泛应用。但由于对湍流 高阶精度的数值离散格式。由于传统的原因，用于湍 机理了解太少，目前各种各样的雷诺应力模型都只对 流大涡模拟研究的程序采用虚拟压缩方法的比较少， 一类特定的流动有效，对于在工程实际中经常面对的 这个方法用于湍流大涡模拟计算的可行性还没有被 复杂流动结果并不理想。湍流的直接模拟方法用足 很好地研究L4J。文献[4]用虚拟压缩方法作了湍流大 涡模拟计算的可行性研究，但没有涉及到像钝体绕流 方程。这个方法提供了流动的所有细致结构，是目前 这样的复杂的流动问题。文献[1]发展了一套基于虚 最准确的湍流计算方法，但由于计算量非常大(要求 拟压缩方法的三维非定常不可压缩流动计算程序。 的空间节点数和时间步数分别达到0(Re9／4)和 程序采用双时间步的时间推进方法，对流项采用五阶 O(Re∽)[8J，在可以预见的将来，直接模拟方法对几 迎风格式。本文将该程序应用到不可压缩湍流方柱 乎所有的工程问题都将是完全无能为力的。大涡模 绕流的大涡模拟计算研究。通过与文献结果比较分 析，表明基于虚拟压缩方法的不可压缩流动计算方 量分解为亚滤波尺度量(“小涡”)和可分辨尺度量 案，应用于复杂流动的湍流大涡模拟计算是完全可行 (“大涡”)。对可分辨尺度量的贡献直接用数值方法 的。 模拟，而亚滤波尺度量的贡献则通过模型来模拟。由 不可压缩三维方柱绕流问题几何外形相对简单， 于采用亚滤波尺度模型，大涡模拟的计算量比直接模 却涉及到物面分离、大尺度非定常涡脱落以及流动转 拟方法要少得多，目前已经在一些工程问题中得到应 捩等复杂的流动现象，是复杂湍流计算研究广泛采用 用。同时由于亚滤波尺度模型是针对具有较多普适 的标准算例。文献[5，6]用激光多普勒速度仪对雷性的小尺度运动，模型构造的困难比雷诺应力模型显 诺数为22000的方柱绕流作了细致的