基于DES方法的三角翼激波-涡干扰流场数值模拟

采用基于Spalart-Allmaras湍流模型的脱体涡模拟（DES）方法,数值求解Navier-Stokes方程,模拟绕尖前缘三角翼的跨音速流动,并对三角翼上翼面的复杂激波-旋涡干扰流场进行了分析。与NASA兰利研究中心的NTF风洞实验结果对比分析表明,DES方法能很好地模拟跨音速三角翼上的旋涡流动。随着攻角由中度攻角增加到大攻角,支架附近的激波越来越强,对主分离涡的干扰作用越来越大,直至出现激波干扰导致的涡破裂。激波的形状、位置及涡破裂位置均与实验结果吻合良好。     

机翼／机身低亚声速最大升力系数预测方法

本文研究了大展弦比机翼／机身低亚声速最大升力系数的半经验预测方法。该方法采用便宜可靠的势流面元法／分离泡／边界层粘性耦合技术与最大升力状态下流动现象的观察经验（压力偏差律）相结合。其中势流面元法采用第二代低阶面元法,每个面元上常值源汇和常值仍极子奇点分布;机翼后缘分离流采用简单有效的虚拟固体边界包围分离区模拟方法。它包括二维和三维间的信息交换,以及粘性／势流间的耦合迭代。ＧＸ－Ⅱ机翼／机身算例表明该方法能够计算机翼／机身详细的亚声速气动载荷和表面压力分布。计算的Ｃｌ－α曲线结果与实验数据有较好的一致性。     

Numerical investigation of unsteady vortex breakdown past 80°/65° double-delta wing

An improved delayed detached eddy simulation （IDDES） method based on the k-x-SST （shear stress transport） turbulence model was applied to predict the unsteady vortex breakdown past an 80o/65o double-delta wing （DDW）, where the angles of attack （AOAs） range from 30° to 40°. Firstly, the IDDES model and the relative numerical methods were validated by simulating the massively separated flow around an NACA0021 straight wing at the AOA of 60°. The fluctuation properties of the lift and pressure coefficients were analyzed and compared with the available measurements. For the DDW case, the computations were compared with such mea-surements as the mean lift, drag, pitching moment, pressure coefficients and breakdown locations. Furthermore, the unsteady properties were investigated in detail, such as the frequencies of force and moments, pressure fluctuation on the upper surface, typical vortex breakdown patterns at three moments, and the distributions of kinetic turbulence energy at a stream wise section. Two dominated modes are observed, in which their Strouhal numbers are 1.0 at the AOAs of 30°, 32° and 34° and 0.7 at the AOAs of 36o, 38° and 40°. The breakdown vortex always moves upstream and downstream and its types change alternatively. Furthermore, the vortex can be identified as breakdown or not through the mean pressure, root mean square of pressure, or even through correlation analysis.     

绕跨声速三角翼的激波/涡干扰流场数值模拟

 在绕三角翼的跨声速流动中,随着迎角的增加,三角翼上的涡破裂位置会出现突然前移的现象。针对这一与亚声速下不同的流动现象,采用带曲率修正的Spalart-Allmaras(SAR)湍流模型,求解定常雷诺平均Navier-Stokes(RANS)方程,对不同迎角下绕65°后掠尖前缘三角翼的跨声速流动进行数值模拟,并在此基础上,采用基于SAR湍流模型的脱体涡模拟(DES)方法,对由激波干扰导致的前缘涡破裂位置的运动规律进行了初步探讨。模拟结果与试验结果对比表明:SAR湍流模型能准确地模拟出三角翼上的激波系统和旋涡结构,并能准确模拟出由于激波干扰导致的涡破裂位置突然前移的现象。此外,对涡破裂后流场的非定常数值研究发现,支架前端正激波的干扰作用使得涡破裂位置向下游移动比较突然,而向上游移动则相对缓慢。     

合成气旋流预混燃烧的大涡模拟

如何高效洁净地燃烧合成气是整体煤气化联合循环(IGcc)系统中面临的主要问题.利用旋流贫燃预混燃烧可以达到强化传热传质,提高燃烧反应效率,降低NOx等目的,但由于合成气中所含有的H2成分,同样面临燃烧不稳定问题.为了对合成气旋流预混燃烧室的燃烧不稳定性及污染物排放提供理论和数值上的指导,采用大涡模拟方法,对典型的旋流预混燃烧室中的合成气燃烧进行了模拟.结果显示,随着H2含量的增加,火焰缩短向上游发展,燃烧室下游涡破碎更加集中.较高的入口雷诺数,加剧了涡团破碎的程度和速度,同时影响中心回流区和几何回流区的大小及分布.随着当量比的减小,火焰褶皱程度加大,燃烧趋向不稳定.增大压力火焰长度会明显缩短,NOx排放量也会随之增大,但NOx的反应速率只有在较低压力下才会对压力变化比较敏感.     

使用RANS/LES混合方法对钝前缘三角翼进行数值模拟

本文阐述了建立在一方程湍流模型基础上的分离涡模型（DES）和延迟脱体涡模型（DDES）,它们很好的模拟壁面附近的小尺度流动。该方法在在大攻角分离流区域使用了亚格子应力的Smagorinsky大涡数值模拟,模拟亚音速旋涡流动。使用大攻角下的钝前缘65°大后掠三角翼模型研究了CFD对涡的发展、破裂以及复杂涡的发展和演变的模拟能力。本文的结果说明了该方法在模拟分离流的可行性。     

基于自适应混合网格的脱体涡模拟

基于混合网格和CGNS（CFD General Notation System）数据结构,建立了一种各向同性加密/稀疏的网格自适应方法。在悬空点的后处理中,让含有悬空点的单元转化为任意多面体,从而简化了自适应单元剖分模版,同时自适应网格单元之间可完全相容,自适应生成的网格能够直接用于可处理任意多面体的流场求解器。将该自适应方法与脱体涡模拟（DES）算法相结合,开展了65°后掠三角翼大迎角流动的数值模拟应用,并与初始网格的模拟结果进行了详细比较。对比表明：采用网格自适应方法适当增加局部网格量,能够以较小的成本迅速提高三角翼背风区的空间分辨率,增强数值模拟对小尺度涡系结构的解析能力,从而弥补了基于混合网格的脱体涡模拟中常用二阶格式计算的空间分辨率相对偏低、不利于湍流多尺度结构精细模拟的不足。     

超声速细长飞行器前置小翼展开过程的非定常特性

针对大长细比飞行器外形超声速大迎角条件下前置小翼展开过程引起的非定常问题,应用结构动网格技术和基于脱体涡模拟（DES）的非定常数值模拟技术进行了研究,获得了详细的小翼展开过程流场结构非定常变化特性,并分析了小翼展开引起的法向力、纵向压心系数等气动特性随展开角度的变化规律。研究结果表明：超声速大迎角条件下前置小翼展开过程对小翼附近区域以及尾舵区域产生了强烈的干扰影响。小翼完全展开后,压心前移4.1%,降低了飞行器静稳定性,法向力系数增加15.6%,提高了气动过载,对飞行器机动性能均产生有利影响。     

运输类飞机平尾失速敏感性及其试飞技术研究

从平尾失速理论研究出发,探讨了可能产生平尾失速的飞行状态,总结了可用于运输类飞机平尾失速敏感性研究的飞行试验技术,研究了平尾失速和机翼失速的相互区别,帮助飞行员在遇到机翼失速或平尾失速时进行必要的清晰判断,并采取适当的改出方法,从而保证飞机飞行安全.研究结果可为未来大型运输类飞机平尾失速敏感性试飞提供很好的技术支持.     

前缘缝翼开口分段对增升装置气动特性的影响

为了研究高升力系统中开口分段的前缘缝翼对增升装置气动特性的影响,对不同构型生成结构化网格,及其流场进行了大量的数值模拟。在同等状态下,结合不同构型的计算结果,分析了开口分段的前缘缝翼对增升装置气动特性的影响。研究表明,机翼失速不仅与开口宽度有关,而且与前缘缝翼上的开口位置有关。     

分区LES/DES混合方法及缝翼三维流场模拟

提出一种分区大涡/脱体涡模拟混合方法,采用高阶空间有限差分方法和二阶隐式LU-SGS时间推进方法,对有限展长缝翼的三维流场进行了数值计算。该混合方法吸收了脱体涡模拟方法和大涡模拟方法各自的优点,在非核心区域使用脱体涡方法进行计算,相对于大涡模拟方法具有较高的计算效率;相对于分区脱体涡模拟中采用的雷诺平均方法,在计算量没有明显增加的条件下提高了对流场各尺度流动的模拟能力。计算结果与风洞实验数据吻合良好,同时缝翼流场主要观测点的压力脉动数据为后续缝翼噪声分析和低噪声优化提供了基础。     

涡升力乘波体发展研究综述

乘波飞行器在高超声速具有良好的气动性能,但偏离设计状态,气动性能则难以保持。为在宽速域范围能一直维持较好的气动性能,研究人员利用定平面形状乘波设计的优势,提出“涡升力”乘波设计。本文将涡升力乘波体的设计方法归纳为基于吻切理论的定前缘型线法、基于激波装配法的波导体法和基于给定激波面的投影法三大类,综述了涡升力乘波体在宽速域气动特性的相关研究进展,并对涡升力乘波体的后续研究提出建议。     

解析诸国争相研制“察打一体”飞翼隐身无人机

“察打一体”无人机以飞翼无尾气动布局、翼身融合体、复合材料结构、雷达吸波材料等特征,具有升阻比大、有效载重量大、隐身性能好等优点,能够对大纵深、高风险和高价值目标进行侦察与打击,诸国争相研制,后续发展空间巨大。     

鸭翼-前掠翼气动布局纵向气动特性实验研究

前掠翼布局由于其潜在的优势,在未来战斗机的研制中将占有日益重要的地位.本实验通过可变前掠翼和鸭式前翼布局的风洞测力实验,重点分析比较了平板机翼在不同掠角下的纵向气动性能以及鸭翼的影响.实验结果表明,前掠翼在大迎角时能有效提高模型的升力系数,小迎角时其升阻比也略优于后掠翼.前掠翼布局能有效推迟失速,具有良好的失速特性;前掠角较大时,升力系数曲线在失速迎角附近有一个升力系数的平台,该布局具有缓失速特性.距离主机翼较远的鸭式前翼(模型M2)在主机翼前掠和后掠情况下,均可改善整体布局的失速特性,增大失速迎角,增强前掠翼布局缓失速的特点.近距耦合鸭翼(模型M3)显著提高了模型在大迎角下的升力系数.另外,主翼前掠和鸭式前翼布局飞行器具有较好的机动性.     

解析诸国争相研制“察打一体”飞翼隐身无人机

“察打一体”无人机以飞翼无尾气动布局、翼身融合体、复合材料结构、雷达吸波材料等特征,具有升阻比大、有效载重量大、隐身性能好等优点,能够对大纵深、高风险和高价值目标进行侦察与打击,诸国争相研制,后续发展空间巨大。     

扰流片对三角翼脱体涡影响的实验研究

在染色法流动显示实验的基础上,进一步深入研究了扰流片在不同标高攻角和来流速度情况下对70°三角翼脱体涡的影响。PIV测试结果表明,在适当位置,扰流片对三角翼的脱体涡的破裂具有一定的延迟作用,在扰流片的诱导涡与三角翼脱体涡方向一致时,脱体涡破裂延迟越靠后。     

基于分离涡模拟的起落架气动噪声研究

用基于SA湍流模式的非定常RANS、分离涡模拟(DES)和延迟分离涡模拟(DDES)分别对四轮基本起落架模型进行了数值模拟,并根据所得的非定常流场计算了表面声压级分布和声压谱。三种方法所用的时间成本大致相同,URANS略低于其它两种,而DDES由于在附着流动区更好地保持RANS特性,故时间成本略低于DES。计算显示,非定常RANS在附着流动区能够得到合理的结果,但不能准确刻画分离流动的流动形态。DES和DDES都能较好地刻画起落架绕流的定常和非定常特性,DDES的结果略好于DES。因为起落架流动属于大分离流动,所以DES也能够得到相对正确的结果。研究结果验证了分离涡模拟在起落架大分离非定常流动预测与噪声预测中的可行性,对减小起落架噪声的方法研究具有一定的意义。