80°%2f65°双三角翼前缘涡破裂和自激摇滚的数值模拟研究.pdfVIP

第十四届全国计算流体力学会议论文集 800／650双三角翼前缘涡破裂及 自激摇滚的数值模拟研究 袁先旭，谢昱飞，陈坚强，邓小刚 (中国空气动力研究与发展中心，四川绵阳621000) 摘 要：本文构造了适当的网格拓扑，首先对70。后掠角单三角翼做了考核计算。以此为基础对 80。／65。双三角翼绕流作了初步计算研究，捕捉到了符合流动规律的涡破裂。还计算了双三角翼由初始 滚转角出发的自激摇滚，定性上与实验结果是符合的． 关键词：双三角翼；涡破裂；自激摇滚；数值模拟 0引言 双三角翼前缘涡破裂影响到其最大升力、升阻比等气动特性，因此一直是工程应用与流动机理 研究者都极为关注的课题u’。双三角翼一般由大后掠角的基本翼在前段加上更大后掠角的边条组成。 张涵信哺1曾给出旋涡沿轴向演化的流动拓扑分析，指出旋涡破裂一定出现在反压区，且泡型破裂对 应于鞍点起始，螺旋型破裂对应于正常点起始，另还存在一种可能的鞍结点型破裂对应于驻点起始。 这些结论得到了张树海、庞勇单三角翼绕流旋涡演化的数值验证。对于双三角翼，其边条前缘与基 本翼前缘各会形成一对集中涡，这些前缘集中涡的运动和相互干扰对涡自身的破裂有相当大的作用。 因此相比作为基本翼的单三角翼，双三角翼的绕流现象和规律更加复杂。 这种复杂性也反映在其动态行为上。唐敏中啦。的低速自激摇滚实验研究表明：800单三角翼只围 绕西F00平衡滚转角发生摇滚；而800／450双三角翼，例如在入射角矿=22。时，发生的是围绕办为约 验也表明：给一个初始滚转角妒600，随入射角仃增加，双三角翼可发生相图为点吸引子哼极限环 _混沌专极限环专点吸引子的变化；他们还指出，动态摇滚迟滞特性与静态滚转力矩曲线上的临界 点相关。张涵信H1以非线性动力学Hopf分叉和鞍结点分叉解释了上述现象。 运用数值模拟研究双三角翼绕流涡破裂，面临的最大挑战就是如何精确预测高廊数、大攻角情 况下大范围的分离流动。其主要技术难题有三：(1)精细网格。为了准确、细致地捕捉涡破裂及其 演化，计算分析和文献[6，7]都表明，网格的构造是关键之一，需要从网格拓扑上构造较为精细的 网格。(2)计算方法。计算方法包含三方面，都极为关键：一是格式精度。要分辨好涡破裂，需要 把握格式粘性的效应，需要应用和发展高精度格式。二是时间推进。涡破裂是非定常的过程，因此 需要应用时间准确的非定常计算，计算效率也需要提高。三是湍流处理，处理好湍流对于计算的准 确性十分重要。国际CFD技术的趋势是发展实用的RANS／LES混合算法，有取得成功的例子。(3)后 处理。计算得到的海量数据需要高效率的后处理，包括如何准确、直观地显示流场，捕捉到涡破裂 的细致结构：如何应用流动拓扑和非线性动力学理论描述涡破裂的规律，分析涡破裂的机理。 基于以上认识，我们在几方面作了探索研究，包括：参考文献建立了一种避免三角翼尖点网格 奇性的适当的网格拓扑结构；针对有详尽实验和大量数值模拟研究结果陋’7。的70。后掠角单三角翼 型(去中心体)进行了涡破裂绕流的初步计算研究；对于该外形还进行了一些状态的自激摇滚计算。 1 计算方法 757 第f蚀扁垒盈计算藏#力学台U论i集 对j一静态绕流．控制方程为完全气体的层漉Ns方程．在体轴系下计算。对于动态绕流．控制方 程为完全气体的屡流XS方程与剐体运动方程的耦台方程组．在风轴系下计算。N-9方程的对漉项 主要采用JST中心格式离散．桔性项采用中心格式离散。定常计算，时间方向采用LU--SGS方法推 进．以保证计算效率：非定常计算，时间方向采用双时间步方法推进，以保证时间精度。刚体运动 方程采用三点后差的差分离教。 边界条件物面上取为无滑移绝热壁．远场外边界都取为无反射开边界。 方法细节可参见文献[8]，整个计算在盎先旭编写的定常与非定常计算程序FLYOD上进行…。 2数值验证 66nz。使用DES方法对70。后掠角单三角翼进行了计算．并与实验结果进行了比较(圈1)． 计算条件是：来流马赫数Ma=0 法和同格。 削2是我们建立的70。后掠角单三角翼多块对接网格．共分6块．全流场网格数为6379016。 通过这样的拓扑结构布置，可以有效地避免头部尖点的网格奇性：由于H_c—H型的构造，可咀将背 风