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

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

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

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

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

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

绕三角翼纵向俯仰大迎角气动特性计算研究

 采用数值计算方法, 对三角翼从0°上仰至90°的动态流场结构进行了计算, 在此基础上, 对三角翼在上仰过程中受到横侧小扰动情况下的流场结构和气动力特性进行了计算研究。给出了三角翼纵向动态情况下的气动力系数变化, 特别是大迎角横侧力矩系数的变化特征, 并对受到横侧小扰动后横侧运动的稳定性进行了计算与分析。结果表明, 机翼的上仰运动延迟了机翼上翼面旋涡的破裂。同时, 随着机翼俯仰角速度的提高, 机翼抵抗旋涡非对称破裂的能力明显增强, 机翼运动的稳定性也明显提高。     

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

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

三角翼大攻角绕流数值模拟研究

采用Harten-Yee的二阶精度隐式TVD格式对N-S方程进行数值模拟,研究了后掠角为65°和70°的三角翼绕流的流场特性。计算结果表明：TVD格式中限制器函数的选取对于三角翼绕流结构的计算结果具有较重要的影响。结合张涵信的旋涡定性分析理论对旋涡的发展过程进行研究,此外,还对二次涡破裂现象和网格分布影响等问题进行了探讨。     

旋涡运动研究的某些进展

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

三角翼大迎角绕流特性数值模拟的网格处理技术研究

研究网格对三角翼大迎角绕流特性数值模拟的影响.结果表明,网格的生成需要结合流动现象,对网格拓扑结构和网格点分布进行选择与搭配.C-H型网格适宜模拟尖前缘分离涡流态,法向网格在一定范围内应等距增长,沿流向逆压梯度较大的区域内适当增加网格点,尾迹区网格则应做上翘处理.Euler方程具有模拟三角翼旋涡及预测涡破裂特性的能力,但对二次涡等粘性引起的流动细节把握能力不足.利用层流假设的N-S方程,通过合适的网格,也可得到满意的计算结果,但对涡破裂后的强烈非定常湍流流动模拟能力不足.采用旋涡螺旋度可准确反映主涡与二次涡流动,描述旋涡的破裂现象.用轴向速度迅速减小并小于来流速度的点作为涡破裂判据似应更合理.     

低背鳍对细长平板三角翼大迎角空气动力的影响

 对细长体平板三角翼和加上两个不同高度背鳍后的组合体在低速风洞进行了六分量天平测力实验，三角翼后掠角82.5°，背鳍当地高度与模型当地半展长比值分别为0.3和0.6，实验迎角范围12°～32°，包括1.66×10.6和2.33×10.6两个雷诺数。实验结果表明：0°侧滑角下，在翼面上发生旋涡破裂前，单独细长平板三角翼的横向力及横向力矩在实验迎角范围内始终为零；加上两个不同高度的低背鳍后，在一定的迎角下，三角翼的横向力及横向力矩开始不为零，流场定常；在更大的迎角下，流场变得非定常。实验结果初步验证了前人关于细长锥体分离涡的稳定性理论，并给出了旋涡失稳后，随着迎角的增大，流场进一步发展的状态。     

以外加辅助件推迟三角翼上涡破裂的初步水洞实验

对尖锐前缘、后掠角为60°的大攻角平板三角翼模型进行了水洞实验。利用附加的小辅助件,改变翼面上方的流场,使机翼前缘肌体涡推迟破裂。流动显示表明:在翼面上的适当位置安放圆弧形或三角形导流体,或在翼面上方另加一辅助小三角翼,能使涡破裂推迟的效果得到显著提高。此结果可供推迟涡破裂来改善飞机气动性能的研究工作参考。     

非定常扰动对涡破裂的影响

在旋涡发生器内进行了涡管内的涡破裂研究。试验表明,不论是泡状或是螺旋状破裂,趋近破裂点时流速均趋于零;非定常扰动使涡轴附近的轴向速度增大,从而延缓了涡破裂的发生;在定常流参数、扰动频率和振幅的适当组合下,扰动可促使破裂状态由泡状转变为螺旋状,大大推迟涡破裂的发生。初步的理论分析表明,上述现象可能来源于流动共振。文中最后还给出了水槽中前缘振动襟翼的实验结果。试验表明,襟翼振动可延缓破裂的发生,且存在一个最佳的频率。试验中还发现,大迎角时,若襟翼不振动,可交替地只在一侧产生前缘涡,襟翼振动可消除这一现象。     

破裂涡流中非定常现象与频率特性实验研究

通过流动显示、表面动态压力测量及热线测量等实验手段,对三角翼破裂涡流中的多种频率成分进行了分析。频谱分析确定了破裂点脉动和螺旋波的频率特征。实验结果表明,螺旋波主频随着弦向位置的增大先是迅速而后平缓减小。前缘涡破裂点振动具有准周期性,在不同的弦向位置上主频大小几乎没有改变,在靠近破裂点的位置有较大的振动能量。实验分析还表明,在破裂涡的流动状态下,虽然没有形成完全分离流,三角翼绕流流场中已经存在涡脱落的现象。     

高马赫数来流条件下斜激波与流向涡对相互作用

针对斜激波入射流向涡对形成的复杂流动干扰现象,在来流马赫数6条件下,利用后掠斜坡式涡流发生器生成流向涡对,采用数值模拟结合激波风洞实验,研究了斜激波与流向涡对相互作用引起的激波变形和旋涡演变。结果表明,斜激波与流向涡对相互作用的剧烈程度大致分为强、中、弱三种。在强相互作用下,激波面向上游凸起的程度最剧烈,呈现类“T”字型特征,其对称面处接近于一道正激波,波后出现回流区和局部流动滞止,使得流向涡对发生破碎;在中等相互作用下,对称面处存在局部正激波,但没有出现回流区;在弱相互作用下,不再出现局部正激波。进一步地分析发现,斜激波与流向涡对相互作用产生的流向涡对破碎现象,不符合经典的激波作用下单个流向涡破碎理论的预测。采用流向涡对对称面处的气流参数,对经典的流向涡破碎理论进行修正后,能够预测本文流向涡对在斜激波作用下是否发生破碎。本文为掌握激波作用下流向涡对破碎现象的规律,提供了重要参考。     

弹体头部形状与弹翼位置对弹翼涡破裂位置的影响

1.引言 近代战斗机和战术导弹,为获得高机动性能,一般均在大迎用下飞行.众所周知弹翼涡的破裂对大迎角导弹气动特性有着重要影响.因此影响涡破裂的因素是设计部门极为关     

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

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

扩压叶栅叶尖非定常流动的涡动力学机理

为了揭示压气机叶尖区旋涡结构与流动非定常性之间的关联，采用URANS对一亚声速平面扩压叶栅在不同攻角下的流场进行了求解，并借助Q判据提取了叶尖瞬态涡系结构。结果表明：泄漏涡的破碎现象能够通过诱导新的涡结构间接作用于相邻通道的叶尖流动，是导致叶尖流场失稳的关键因素。在-0.3°和+0.7°攻角下，叶尖泄漏涡发生了螺旋破碎，并伴随有非定常诱导涡的出现，诱导涡对相邻叶片载荷的影响使得叶尖泄漏涡发生周期性摆动；在+1.7°攻角下，泄漏涡破碎会导致反流涡的形成，反流涡的输运会给叶片载荷和来流攻角带来非定常扰动，反过来又会作用于泄漏涡的破碎和反流涡的生成，最终表现为一种自维持的非定常流动现象。     

有限长圆管内刚体旋转流的直接数值模拟

通过对三维不可压缩粘性流体的直接数值模拟,研究了有限长圆管内刚体旋转流在初始扰动下,扰动的增长和动力学演化过程.首先,通过改变流动的Reynolds数和旋转角速度ω,分析了轴对称流动中Reynolds数和ω对流动失稳的影响.对于高Reynolds数轴对称流动,利用数值模拟得到的不稳定结果与无粘线性稳定性理论预测的结果一致.最后,为进一步研究流动不稳定发展的三维效应,又进行了完全三维流动的模拟.数值结果表明,三维流动中不稳定的主导模态为螺旋模态,模态在线性段指数增长,导致流动产生螺旋型的漩涡破裂,之后经过非线性段的增长后达到饱和,流动最终发展为轴对称泡型旋涡破裂的稳定状态.