Pω增强型k-ω湍流模型在三角翼旋涡流动的应用

在三角翼旋涡绕流数值模拟中,标准 Wilcox k-ω湍流模型生成项未考虑旋度的影响而导致预测的旋涡强度较弱。通过引入探测因子区分剪切层和涡核,在旋涡流动的高旋度区域增加ω方程生成项的方法,基于结构化网格上的 RANS 求解器,加入了 Pω增强型 k-ω湍流模型,对绕尖前缘三角翼亚声速和跨声速旋涡流场进行了数值模拟。计算结果与 NASA 的 NTF 风洞和 DLR 的 DNW-TWG 风洞试验数据进行了对比分析,结果表明：不论在亚声速还是跨声速自由来流条件下,Pω增强型 k-ω湍流模型计算的压力分布、涡破裂位置均与试验数据吻合良好,准确地预测出了三角翼上翼面的主涡、二次涡结构,特别是跨声速条件下激波干扰导致的涡破裂的临界迎角及涡破裂位置,表明 Pω增强型 k-ω湍流模型在绕三角翼旋涡流动数值模拟中具有良好的适用性。     

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

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

小展弦比飞翼跨声速典型流动特性研究

小展弦比飞翼标模为国内自主设计的融合体飞翼通用研究模型,前缘后掠角为65°,展弦比为1.54。风洞试验结果表明小展弦比飞翼标模在跨声速迎角4°开始出现非线性升力,在迎角12°至16°范围内会出现升力突然下降、俯仰力矩突然上扬的现象。为了分析该现象的机理,通过数值模拟的方法研究了小展弦比飞翼标模在马赫0.9时的流动特性,分析了前缘涡的产生、发展直至破裂的整个过程,结果表明:小展弦比飞翼标模在迎角4°开始出现涡升力;随着迎角增加,前缘涡逐渐向内侧移动,涡强和背风面激波的强度也逐渐增加,前缘涡与激波发生交叉干扰并达到一个平衡流态;当前缘涡与激波无法维持既有平衡时则会发生涡破裂,流场急剧变化以达到新的平衡,从而导致升力突然下降并产生抬头力矩增量。     

高阶精度格式WCNS在三角翼大攻角模拟中的应用研究

采用5阶精度的加权紧致非线性格式(WCNS-E-5)数值模拟了65°三角翼的大攻角绕流流场,主要目的是考核高阶精度格式WCNS在大攻角旋涡流动方面以及跨声速流场的激波附面层干扰、涡破裂位置的模拟能力,重点研究不同网格规模和湍流模型对尖前缘三角翼涡系之间的相互作用的影响。通过求解任意坐标系下的雷诺平均N-S方程,采用5阶精度的加权紧致非线性格式(WCNS-E-5)和多块对接结构网格技术,两种湍流模型分别是一方程SA和两方程SST湍流模型,在与相应试验结果对比的基础上,详细研究了WCNS-E-5格式在跨声速大攻角旋涡流动中的表现,以及不同网格规模、两种湍流模型对主涡二次涡相互作用、涡破裂位置和表面压力分布的影响。本文的研究结果表明,高阶精度格式WCNS-E-5能成功应用于三角翼的跨声速大攻角流动,网格规模的增加进一步提高流场分辨率,SST湍流模型相对SA湍流模型在三角翼大攻角流动中具有更好的适用性。     

细长旋成体背风面非对称流场的DES模拟

利用脱体涡模拟DES方法对细长旋成体非对称绕流进行数值模拟,比较了基于S-A湍流模型的DES与RANS方法对背风面非对称分离涡的数值模拟能力,发现在大迎角非对称涡未破裂情况下,DES与RANS均能模拟出与实验相符的非对称绕流;当迎角增大,背风面非对称涡发生破裂时,RANS无法准确捕捉到背风面流场的非定常性,而DES能准确预测非对称涡的飘起与破裂,并与实验值接近。计算结果表明:与RANS方法相比,DES方法具有更好的模拟大攻角,大分离流动的能力,尤其在非对称流场的分离涡破裂模拟方面具有明显优势,能够更真实模拟出细长旋成体背风面分离涡破裂之后的非定常流动特征。     

马赫数和雷诺数对钝前缘三角翼流场的影响

基于Spalart-Allmaras(SA)湍流模型,求解定常RANS方程,数值模拟后掠角为65°的钝前缘三角翼的亚音速流场。分析不同雷诺数、马赫数下钝前缘三角翼上表面的复杂涡系结构及其发展过程,并与实验结果进行了对比。计算结果表明:在中度迎角下,所得的各弦向站位的物面压力系数与实验数据吻合良好;随着雷诺数的增加,三角翼前缘流动分离位置后移,其背风面涡系向内移动;而马赫数的增加则导致分离位置前移,背风面涡系向外移动。     

翼身相对厚度对小展弦比飞翼布局跨声速气动特性及流动机理的影响研究

为研究翼身相对厚度对小展弦比飞翼布局气动特性以及涡流特性的影响,基于已有试验结果的翼身相对厚度为0.16的65°后掠小展弦比飞翼布局,在保持前缘半径和外翼剖面形状相同情况下,通过降低飞翼布局的翼身厚度使其翼身相对厚度为0.08,在马赫数0.9条件下开展了翼身相对厚度影响的数值模拟研究。数值模拟结果表明,在相同迎角条件下,翼身相对厚度对飞翼布局前缘涡在翼面上形成的位置和涡强有较大的影响,翼身相对厚度较小时前缘涡形成的位置越靠近前缘;在前缘(约x/Cr=0.25之前)翼身相对厚度较小布局的涡核强度明显高于翼身相对厚度较大布局,且在前缘涡破裂之前,翼身相对厚度较小布局涡核强度沿弦向变化较为平缓,升力线斜率下降迎角较翼身相对厚度较大布局推迟约8°。研究结果还表明跨声速时,前缘涡的破裂主要与激波的干扰有关,当前缘涡穿过激波时,涡强和涡核轴向速度迅速降低,当涡核轴向速度降为0时,前缘涡破裂。     

后缘喷流对三角翼前缘涡的影响

本实验应用染色液流动显示技术和激光测速技术 (LDV)研究了 6 0°后掠三角翼在后缘差动喷流、对称喷流情况下前缘涡破裂位置、涡核的空间分布、涡核的速度分布以及三角翼背风面流动结构随迎角的演化等。实验结果表明 ,喷流增大了三角翼前缘涡涡核保持高速度的区域 ,推迟了涡核减速的位置 ,在大迎角情况下 ,对称喷流有助于消除由前缘涡振荡引起的“摇滚”现象。     

三角翼大迎角绕流的数值模拟

采用FINETM/HEXA非结构网格流场计算软件包对三角翼飞行器的低速大迎角绕流进行了数值模拟.结果表明,计算值与实验数据符合较好,采用Agglomeration 方法能大幅度提高收敛速度,而网格自适应方法用总体不大的网格单元数可较精细地模拟流场.本文还分析了前缘分离涡破裂前后的流动现象和旋涡横截面流线图谱的变化规律.     

考虑湍流影响的前缘分离涡运动与破碎特性的研究

本文以Ｎ－Ｓ方程和Ｈａｌｌ涡核模型假设为基础，导出了描述湍流的涡核运动方程。利用差分计算方法，对三角翼前缘分离涡运动及其破碎特性进行了数值求解，分析了涡核流场的结构和各物理量的变化特性，反映了三角翼前缘分离涡运动及其破碎特性的实质，说明了湍流涡核方程能更有地模拟三角翼前缘分离涡运动和准确地确定涡核破碎位置，并得出了攻钐雷诺数对破碎位置的影响，计算结果与实验数据十分吻合。     

三角翼跨声速动态失速与涡破裂特性研究

通过数值方法研究了高速气流中细长三角翼作匀速上仰机动飞行时的背风面分离流与涡破裂特性和气动力特性。结果表明,在三角翼匀速上仰非定常绕流中,涡破裂起始攻角与上仰速度关系比较复杂,在本计算速度范围内,涡破裂起始攻角变化很小;动态时涡破裂点发展速度随攻角变化规律与静态结果差别甚大;三角翼大攻角非定常绕流前缘涡涡轴附近的流动以不稳定涡和涡破裂为主要特征;在一定的上仰速度以后即使在中等攻角前非定常升力与定常升力仍有较大差别,与低速情况有所不同;上仰运动使涡破裂点发展速度被延缓,从而提高了失速攻角和最大升力,并使失速攻角与涡破裂起始攻角的间隔进一步拉大,同时不改变升阻比     

Transonic wing stall of a blended flying wing common research model based on DDES method

Numerical simulation of wing stall of a blended flying wing configuration at transonic speed was conducted using both delayed detached eddy simulation(DDES) and unsteady Reynolds-averaged Navier-Stokes(URANS) equations methods based on the shear stress transport(SST) turbulence model for a free-stream Mach number 0.9 and a Reynolds number 9.6 × 10~6. A joint time step/grid density study is performed based on power spectrum density(PSD) analysis of the frequency content of forces or moments, and medium mesh and the normalized time scale0.010 were suggested for this simulation. The simulation results show that the DDES methods perform more precisely than the URANS method and the aerodynamic coefficient results from DDES method compare very well with the experiment data. The angle of attack of nonlinear vortex lift and abrupt wing stall of DDES results compare well with the experimental data. The flow structure of the DDES computation shows that the wing stall is caused mainly by the leeward vortex breakdown which occurred at x/x_(cr)= 0.6 at angle of attack of 14°. The DDES methods show advantage in the simulation problem with separation flow. The computed result shows that a shock/vortex interaction is responsible for the wing stall caused by the vortex breakdown. The balance of the vortex strength and axial flow, and the shock strength, is examined to provide an explanation of the sensitivity of the breakdown location. Wing body thickness has a great influence on shock and shock/vortex interactions, which can make a significant difference to the vortex breakdown behavior and stall characteristic of the blended flying wing configuration.     

带副翼偏转的三角翼自由滚转运动数值模拟

 通过耦合求解非定常Euler/Navier-Stokes方程和单自由度滚转运动方程,对带副翼偏转的65°后掠角尖前缘三角翼WI1-SLE自由滚转运动进行了研究,Navier-Stokes方程的求解采用基于Spalart-Allmaras湍流模型的脱体涡模拟(DES)。在多块结构网格上,应用基于弧长的无限插值理论(TFI)生成变形网格,实现副翼偏转,而三角翼的滚转运动则通过网格的整体旋转实现。结果表明:Euler方程和DES方法均准确地模拟出了三角翼在滚转运动过程中存在的3个平衡位置。出现平衡位置的原因分别是:①流动对称性;②机翼左侧发生涡破裂的分离涡与右侧分离涡相互平衡使得滚转力矩为0,并且平衡位置仅与三角翼两侧涡强的差有关;③副翼偏转和左右机翼不对称分离涡涡强差产生的滚转力矩相互平衡。此外,滚转运动对副翼偏角幅值很敏感,幅值的微小改变会影响最终的平衡位置和向平衡位置运动的路径。     

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

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

旋涡运动研究的某些进展

 本文简要介绍研究旋涡运动在以下问题上的某些结果:低速不同后掠角三角翼在各个迎角下的九种分离流类型及其边界;应用微分方程定性论与拓扑学对三维分离流与旋涡流的分析;旋涡破裂形态,对三角翼前缘涡破裂的实验研究与理论分析;受控分离与旋涡的干扰,二旋涡的位移、绕转与合并等。     

洞壁对过失速非定常三角翼涡破碎位置的影响

在南航低速风洞中用两组后掠角分别为６５°和７０°的三角翼模型进行了过失速非定常涡破碎位置测定实验。各组模型几何相似,展长与风洞宽度之比分别为０．１７５,０．３５和０．７。涡及其破碎点位置由ＴｉＣｌ４烟流显示并由相机记录。实验表明,同样攻角条件下,随着模型加大涡破碎点位置不断后移     

双三角翼前缘剖面形状对涡运动的影响

采用数值模拟和理论分析相结合的方法,研究了前缘剖面形状对双三角翼涡运动的影响,分析了前缘剖面形状对三角翼、双三角翼涡运动影响的不同机理 :对三角翼,尖前缘可以形成组织最好的涡结构,但对于双三角翼,圆前缘生成的旋涡结构较靠近翼面,涡结构紧密,诱导能力较强,可以形成有利的涡涡干扰,使内翼涡通过剪切层向外翼涡输入涡量更加容易,合并涡变得更加稳定,推迟了涡破裂,而且由于涡较靠近翼面,因而可以产生较高的非线性涡升力,这同传统的认识是不一致的。