约束大涡模拟方法在分离流中应用.pdfVIP

北京力学会第18届学术年会论文集：流体力学 约束大涡模拟方法在分离流中的应用 夏振华 史一蓬 肖左利陈十一 (北京大学工学院，100871) 摘要：在本文中，约束大涡模拟方法被用来模拟周期丘陵绕流和NACA002l翼型在60攻角 下的大规模分离绕流。在周期丘陵绕流的计算中，约束大涡模拟方法能相当准确地捕捉流 动的分离和再附位置，计算的结果不比分离涡模拟方法计算的结果差。在NACA0021翼型 的模拟中，我们能得到和实验符合得很好的压力分布，流动分离的主频等。从这些结果上 看，约束大涡模拟方法也可以被用来模拟分离流动． 关键词：约束大涡模拟方法，分离流，分离涡模拟方法 一、序言 高雷诺数壁湍流问题一直是科研和生活实践中常见的问题。对这类问题的研究，直 接数值模拟方法和传统的大涡模拟方法都冈为目前的计算能力的限制而尤能为力；雷诺 平均方法虽然在工程应用中备受青睐，但是它在非定常，分离等非平衡流动中的表现很 难得到满意。近些年，带譬面模型的大涡模拟方法和RANS／LES混合方法是研究这类问 题最多的方法，其中RANS，LES混合方法史被认为是今后数十年研究这类问题的主流方 法。虽然RANS／LES混合方法已经取得不少成果，但是它却还有一些问题，其中最突出 的就是在模拟附着壁湍流(如边界层湍流，槽道湍流等)时，出现的对数律不匹配问题 (LLM)。目前的一些解决方案都有明硅的非物理凶素，而且效果不够理想。 最近，史一蓬等人…针对这一问题，从约束的角度，提出了一种可模拟高雷诺数壁 湍流的计算方法：约束大涡模拟方法。约束大涡模拟方法的思想和传统的RANS／LES 混合方法的思想不一样，不再是基于雷诺平均方法和大涡模拟方法在方程及模犁上的高 度相似性，而是从大涡模拟方程的系综平均和雷诺平均方程的等价性出发，找到了大涡 模拟方法在模拟壁湍流时，亚网格模型应力应该满足的约束等式，从而从约束的角度， 准确计算了亚网格模型应力的平均部分。在约束大涡模拟方法里，外层使用传统的亚网 格应力模型，内层使用带雷诺应力约束的亚网格模型，虽然全场计算的是大涡模拟方程， 但是在内层，其隐含着求解了雷诺平均方程。在文掣u【21中，作者选用了槽道湍流，圆 柱绕流等基本算例，对该方法做了摹本验证。槽道湍流的结果显示，约束大涡模拟方法 能E,很z好的消除LLM现象，同时能够更准确的计算壁面摩擦系数，而且该方法对界面和 网格解析度有一定的适定性。在圆柱绕流的计算中，发现约束大涡模拟方法计算的结果 和分离涡模拟方法的结果很相近，在计算统计的时间间隔里，约束大涡模拟方法不比分 离涡模拟方法差。 二、算例 2．1周期丘陵湍流 SimulationofFlowAround 周期丘陵湍流是欧洲LESFOIL(Largecddy Airfoils)和 II．63 |匕京力学套第18月学术年套☆i集：m*^学 Research onFlowTurbuleneeand Community ERCOFTAC(EuropCan Combustion)等酬队朋 米检验存在分离现象时湍流模型车身及其模拟效果的一个杖计算例。在雷诺数达到定 程度之后．它存在严重的分离。就问题本身而言，它的几何外形简单(二维几何构掣)， 分离点的仃置小像后台阶流动样受儿何外形所托制，而且边界条件简单fH有周埘平u 嗣肇岍种条件1。这些特征·Ⅱ以JS量减少外界=jl‘确定…索封模掣和模拟问题的影响。奉 文中删期丘陵的训算嘲格如图1 2．2 60度攻角的NACA0021曩型绕滴 大攻角F的机翼绕流00题 血足RANS疗法的软肋，也是DES方法优势所在。60