方柱绕流大涡模拟.pdfVIP

方柱绕流大涡模拟 谢志刚1，许春晓2，崔桂香2 (I．航空六--Jk所。西安710065；2。清华大学，北京，100084) 摘要：采用有限体／有限元混合格式求解N．S方程，利用较稀疏的非结构网格对雷诺数为22，000的学方柱绕流进行 了人涡模拟，得到的结果与实验结果符合得很好，相对于以往的数值研究有较大改进。另外。还测试丫不同的弧格予模式和 边界条件以及低马赫数顶处理方法对数值模拟结果的影响，并模拟了绕两个并列方柱的流动，揭示了方柱间距不同时的不『r可 流态及流动特性。 关键词：大涡模拟，有限体／有限元混合格式，非结构网格，方柱绕流 O 引言 布，是目前为止最好的数值模拟结果，但是阻力系 数和回流区长度仍然与实验结果差别较人，尾迹内 湍流的大涡模拟一直是湍流研究的热点问题， 的恢复速度也明显大于实验值。以往的研究者都认 人们寄希望于在不久的将来火涡模拟能真正成为 为网格精度是决定计算结果的关键冈素，在Rodi 工程上预测复杂流动问题的有力工具。但是，目前 和Voke组织的两个研讨会中不同算例的网格数为 的大涡模拟研究主要针对一些典型的简单湍流流 动，大多采用规则网格和高精度的计算方法，而工 要130万个节点。如此规模的网格是难以应用于工 程问题的模拟经常需要采用非结构网格和低精度 程实际的，因此本文采用较稀疏的非结构网格模拟 的离散方法，对此种条件下大涡模拟的有效性和精 方柱绕流，并测试不同的亚格子模式和边界条件以 度的研究相对较少。因此，本文采用非结构网格和 及低马赫数预处理方法对结果的影响。 二阶精度的有限体／有限元混合格式，对方柱绕流进 行大涡模拟，以探索大涡模拟在j[程实际中的应用 1数值方法 前景。 本文采用有限体／有限元混合格式求解过滤后 方柱绕流包含多种复杂的流动现象，如剪切层 的可压缩的N．S方程，有限体厂有限元混合方法需要 的转捩，流动的分离和再附，周期性的涡脱落，尾 两套重合的网格，首先将三维计算域划分为四面体 迹的形成与发展等，这些流动现象与门然界和一r程 有限元网格，然后在每个顶点上构造一个有限体单 中的许多流动问题有关。Lyn【l】对Re=21．400的 元形成有限体网格。将过滤后的控制方程在-f『限体 方十}{绕流进行了细致的流场测量，给出了流动速度 单元内积分，利川ROE格式15J与MUSCL插值方 及其脉动量的分布，为数值模拟提供了事实依据和 法计算对流通量，利刚有限元方法计算粘性通餐， 考核标准。Rodi【21和Vokel31分别组织的两个研讨 然后采用四阶龙格库塔方法进行时问推进。本方法 会发现人涡模拟得到的平均阻力系数Cn和同流区 的精度为空间二阶、时间四阶。 K度船对计算条件很敏感、误差很大，尾迹中的 记方枉高度为d，米