基于特征线分离算法的大涡模拟.pdfVIP

第27 卷第2 期 Vol.27 No.2 工 程 力 学 2010 年 2 月 Feb. 2010 ENGINEERING MECHANICS 41 文章编号：1000-4750(2010)02-0041-08 基于特征线分离算法的大涡模拟 *梁开洪，曹树良，陈 炎，祝宝山 (清华大学热能工程系水沙科学与水利水电工程国家重点实验室，北京 100084) 摘 要：将基于特征线的分离算法与大涡模拟相结合，推导了不可压流大涡模拟有限元离散方程组，并将该方法 应用于三维流场的层流及湍流非定常计算。将不同雷诺数下的三维顶盖驱动空腔流动计算结果与实验数据以及直 接数据模拟结果进行对比，吻合较好，验证了方法的可靠性和准确性。 关键词：大涡模拟；有限元；特征线法；分离算法；方腔流动 中图分类号：O241; O357.1 文献标识码：A LARGE EDDY SIMULATION COMBINED WITH CHARACTERISTIC-BASED SPLIT SCHEME * LIANG Kai-hong , CAO Shu-liang , CHEN Yan , ZHU Bao-shan (State Key Laboratory of Hydroscience and Engineering, Department of Thermal Engineering, Tsinghua University, Beijing 100084, China) Abstract: This paper introduces the characteristic-based split scheme into the large eddy simulation. The finite element discrete equations for the impressible flows are derived in the framework of large eddy simulation, and then these equations are applied to solve three dimensional unsteady flow problems. Numerical simulations of the flow in a lid-driven cubical cavity were conducted at different Reynolds numbers (Re=1000 and Re=10000). The results are in good agreement with experimental data and direct numerical simulation results. Key words: large eddy simulation; finite element; characteristic-based method; split scheme; lid-driven [4] 有限元法在数值模拟方面占有重要地位，但在 及特征线Galerkin 法 等。其中特征线Galerkin 法 流体动力学方面的应用没有有限体积法广泛。究其 沿着流体质点的运动轨迹(特征线)进行时域离散， 原因，主要存在两方面困难：首先动量方程中由于 以理性形式引入稳定项，具有比较明确的物理意 存在对流项导致方程不具备自伴随性，当对流项占 义，是一种处理对流占优问题的有效方法。 优时，应用标准Galerkin 离散过程求解将导致数值 为满足 LBB 条件，一般要求压力插值函数比 结果虚假振荡；其次对于不可压流体，由于压力变 速度插值函数低一阶，这一限制使研究者不能任意 量不出现在连续