用N－S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二级时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型￣［１］。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

用N—S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二维时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

动态跨音速流动控制中自适应翼型的设计理念

通过对仿真几何模型的数值模拟，有针对性的消除了跨音速区的再压激波．改进了应用在非定常跨音速流动中的翼型，并介绍了一套设计直升机旋翼叶片可变形翼型的基本步骤，而且结合系统优化方法，还设计出了一种无激波翼型。同时还考虑了两种关于变形的分析模型，并应用二维、时间精确、隐式Navier-Stokes程序检验了这两种模型控制动态失速的能力。结果显示，对于直升机，受控变形翼型的发展潜力很大。应用这一翼型，能够大幅度的抑制大攻角时翼型的动态失速与小攻角、高马赫数时再压激波的强度，其结果将会帮助拓展直升机旋翼翼型的设计理念，使其能够进一步控制翼型的动态失速。     

低雷诺数下翼面局部振动增升机理研究

采用计算流体力学(CFD)方法研究低雷诺数下翼面局部振动对翼型气动特性及其流动特征的影响规律。建立局部振动激励的力学模型,并采用任意拉格朗日-欧拉坐标系下的特征线有限元(ALE-CBS)方法对局部振动激励下翼型绕流问题进行模拟,分析局部振动对非定常流动演化的影响规律,揭示其增升机理。研究结果表明:翼面局部变形的增加会有效降低翼型上表面的前缘压力;非定常流动分离中旋涡之间的距离及其演化频率与振动频率的关系是影响翼型翼面局部振动增升效果的重要因素。当流场主频率与振动频率相同,次要频率为主频率的2倍,即发生锁频时,翼面振动产生的移动分离泡能够使分离区从主流获取更多的能量,使翼型上表面保持较低的压力,有效提高翼型升力。     

翼型前缘变形对动态失速效应影响的数值计算

翼型或机翼的动态失速效应所引起的低头力矩和正气动阻尼限制了飞行器气动性能的提高,甚至可能诱导发生不稳定运动。应用于小尺寸机翼的前缘动态变形（DDLE）技术,通过实时改变前缘形状,能够改善翼型前缘区域的速度梯度,进而抑制动态失速效应。采用转捩剪切应力输运（SST）黏性模型结合分区混合动态网格技术,研究了这种前缘变形对机翼俯仰运动所引起的非定常流动的影响,得到通过小幅度前缘变形抑制和延迟动态失速的方法,从而提高翼型的气动性能。翼型NAC A0012的数值模拟结果与动态失速风洞试验结果比较表明：所使用的数值计算方法能够较为准确地模拟翼型在动态失速过程中升力系数与俯仰力矩系数的变化情况,可用于研究前缘变形对翼型俯仰运动所引起的非定常流动的影响。前缘动态变形翼型俯仰运动过程的非定常流场的数值模拟表明：在大迎角下不同幅度的前缘下垂运动能够抑制流动分离的发生,从而抑制动态失速,但在大迎角下小幅度高频率的前缘下垂变形能更高效地抑制动态失速;前缘变形幅度以及变形沿中弧线的分布对升力系数和俯仰力矩系数的影响并不明显。     

压气机二维叶栅涡脱落的数值模拟

对某压气机二维叶栅的非定常分离流场进行了数值模拟,通过对多种工况的计算,进行了频谱分析,对叶栅非定常流动的流场结构和流动机理做了初步的探讨。在来流均匀,定常边界条件下,叶栅内流动仍然表现出强烈的非定常性;在大攻角下有类卡门涡的脱落;旋涡脱落的频率,随着攻角和马赫数的变化而变化;分离点的位置随着攻角和马赫数的变化而变化。     

三维直叶栅非定常流动的并行计算研究

流动分离直接关系到压气机运行的安全性与效率,对分离流动的研究是叶轮机械真实流动研究中的一个重大课题。本文针对三维压气机单转子叶片中截面所构成的三维直叶栅跨音速分离流开发了通用数值计算程序,该程序基于B-L湍流模型及高精度差分方法。多种工况的数值计算显示本程序结果与实验值吻合比较理想,验证了程序的正确性。10°攻角下分离区脉动压力的频谱与实验结果的数量级吻合,说明本程序能够较好地模拟大攻角分离流这种非定常复杂流动。为了提高计算规模及计算速度,作者对程序进行了并行化并针对微机机群系统进行了并行优化。实际计算表明本程序具有较高的并行效率。      

模拟厚度翼型非定常运动的涡方法

本文研究二维厚度翼型的大攻角非定常运动。用沿翼型表面分布的源和面涡来模拟翼型的厚度效应和升力效应,应用离散涡方法模拟前缘分离涡层和尾涡层。计算了大攻角翼型作俯仰运动的情形,与实验结果符合较好。     

翼型和机翼阵风响应的数值模拟

采用有限体积法来数值求解非定常Euler和NS方程,并在非定常计算中通过引入"网格速度"来计入阵风的影响,实现了对刚性翼型和机翼阵风响应的有效数值模拟。首先采用该方法来对翼型攻角突然增大的阵风、陡边型阵风这两种典型阵风情况下的翼型气动力响应进行了计算,结果表明:计算结果与现有的一些理论结果都符合较好,与文献中计算结果也吻合,NS方程的计算结果与Euler方程计算结果在变化趋势上是一致的,但具体的数值有差别。最后针对攻角突然增大的阵风,对不同马赫数下不同展弦比平直机翼的阵风响应进行了数值模拟,得出的计算结果与文献中的计算结果一致。     

基于零质量射流技术的旋翼翼型流场主动控制研究

旋翼翼型经常工作在大迎角来流条件下,流动分离普遍存在。通过在具有分离流动的翼型表面施加零质量射流控制,可以实现对流动分离的抑制,达到增升减阻、推迟失速的目的。本文通过对某旋翼翼型的流动控制的计算,采用数值模拟方法研究了零质量射流主动流动控制机理,并对影响零质量射流控制效果的射流喷角、动量系数和无量纲频率三个控制参数进行了分析比较,得出了一些有意义的结论。     

翼型大攻角绕流的数值模拟

以求解二维N-S方程数值模拟NACA0012翼型大攻角状态的可压绕流特性;N-S方程是在贴体坐标系中给出的,以代数方法生成C型网格系统。采用LU-ADI格式计算,为提高格式的稳定性在隐式和显式部分分别添加了2阶和4阶人工粘性项。应用BaldwinLomax湍流二层代数模型模拟了大攻角时前缘分离涡的形成,旋涡对流及其非定常现象。在某些Mach数和攻角下NACA0012翼型的湍流解具有周期性。通过比较,本文数值计算结果同实验及国外相应的数值计算结果基本吻合。     

Computations of unsteady transonic flow over airfoil at low Reynolds number

Transonic flow over a thin airfoil at low Reynolds number was studied numerically by directly solving two-dimensional full Navier-Stokes equations through 5th order weighted essentially non-oscillatory(WENO) scheme without using any turbulence model.A series of distinguished unsteady phenomena for a thin 2-D transonic airfoil flow were presented.Due to continuous adverse pressure gradient in the subsonic flow downstream of the sonic line, the unsteady separated boundary layer with main vortex and secondary vortex was developed at the rear of the airfoil.At the trailing edge,the vortex-shedding was characterized by periodical connection of the main vortex and secondary vortex on the other side of the airfoil.The unsteady separation and vortex-shedding occurred with the same period.On the airfoil surface,the average pulse pressure related to the unsteady supersonic region was obviously smaller than that related to the vortex-shedding at the trailing edge.With the attack angle increasing from 0° to 2°, the frequency of vortex-shedding decreases about 4.2%.At last, the turbulence intensity and many second-order statistics in the wake region were investigated.      

等离子体气动激励改善增升装置气动性能的试验

针对流动分离导致飞机增升装置气动性能下降的问题,进行了脉冲等离子体气动激励抑制增升装置流动分离的试验。研究了等离子体气动激励的频率、占空比及激励位置等参数对流动控制效果的影响。研究结果表明：等离子体气动激励通过加速近壁面附面层,增强附面层内的能量掺混,可有效抑制主翼和襟翼表面的流动分离,改善增升装置气动性能。在主翼前缘施加激励,可有效控制主翼表面大迎角下的失速分离,最大升力系数增大18.1%、临界失速攻角提高4°;在襟翼前缘施加激励,可有效抑制襟翼表面的流动分离,显著减小阻力,在4°迎角下,将试验模型阻力系数减小了28.7%,升力系数提高了7.1%。占空比对控制效果有较大影响,当占空比为10%～30%时,激励的非定常性更强,控制效果最好;占空比为50%的控制效果次之,占空比为100%时的控制效果最差。来流速度越高,逆压梯度越大,流动分离更难被抑制,控制效果也变差。该研究为在增升装置上应用等离子体流动控制技术提供了理论和方法的基础。     

低雷诺数下螺旋桨翼型非定常气动性能的比较

为得到适用于平流层螺旋桨桨叶叶素的翼型,采用数值方法,通过改变气流速度和攻角,比较研究了低雷诺数下8种典型翼型的非定常气动性能.结果显示,大攻角工况下翼型容易发生流动分离现象,翼型吸力面的分离越严重,气动力的振荡幅度越大.通过比较不同攻角和不同速度下翼型的升阻比大小和气动稳定性得出,SD 8000 - PT翼型具有较高的升阻比和较稳定的气动力,适合作为螺旋桨的翼型.     

低雷诺数下微型飞行器主动变形机翼非定常气动特性数值分析

以微型飞行器的气动外形设计为研究背景,通过数值求解N-S方程分析了主动变形柔性机翼的气动特性。空间离散格式采用中心有限体积法,时间推进采用双时间推进方法,其中子迭代过程由隐式LU-SSOR方法完成。在动网格技术的基础上,模拟分析了薄翼面上作行波运动的流场。计算结果发现,该行波翼型的气动特性相比于刚性翼型有明显改善,有效地抑制了大攻角下大尺度流动分离,同时升阻比提高了约35%,起到了显著的增升减阻的作用。本文的工作对于进一步的柔性机翼优化具有良好的理论研究与数值模拟的参考价值。     

翼型动态失速等离子体流动控制试验

针对动态失速引起的翼型气动性能恶化的问题,利用小型化的激励电源和介质阻挡放电等离子体激励器,借助动态压力测量和外触发式粒子图像测速（PIV）等手段开展了翼型动态失速等离子体流动控制试验研究。结果表明,等离子体气动激励能够有效控制翼型动态失速,改善平均气动力,提高翼型气动效率,减小气动力随迎角变化的迟滞区域。等离子体诱导出前缘附近的贴体翼面涡,促进分离流再附;增加了上翼面0.2~0.4弦长区域的吸力,减小了升力系数功率谱密度（PSD）分布的二、三、四阶能量幅值,在研究工况下实现了平均升力系数增加7.1%、失速迎角推迟1.3°和迟滞区域减小4.5%的明显控制效果;4°~9°迎角段,等离子体使得翼型平均阻力系数减小40%。此外,振荡频率增加使翼型绕流的非定常性增强,较高雷诺数下的翼型动态分离涡更加难以被抑制,均需要增加等离子体激励强度才能达到较好的控制效果。     

基于RANS/LES混合方法的分离流动模拟

飞行器在大迎角、快速俯仰机动时,流场中含有大尺度、非定常的涡结构,传统雷诺平均Navier-Stokes （RANS）模型不能准确模拟流场结构,根据国际上相关研究的发展趋势,需要采用混合RANS/大涡模拟（LES）模型来对复杂分离流动进行准确模拟。本文对基于分区混合与湍流尺度混合的双重RANS/LES混合计算模型进行发展与应用。通过典型简化模型的静、动态湍流大分离流动,测试和验证所采用的脱体涡模拟（DES）类方法,重点研究改进的延迟DES （IDDES）模型在动态问题应用中的正确性和有效性,并对所采用的数值模拟方法和相应的计算软件的可靠性、鲁棒性以及精度进行了考核验证。典型算例包括超声速圆柱底部流动、跨声速方腔流动、NACA0015机翼深失速分离涡模拟等。计算表明：发展的IDDES类混合计算模型可有效解决对数层不匹配的问题;对于定态非定常分离流动,DES、DDES、IDDES等模型计算结果差别不大,随着流动的非定常特性增强,IDDES模型的优势逐渐显现;对于动态非定常分离流动,则需要采用IDDES类模型。     

锯齿尾缘叶片气动特性和绕流流场的数值研究

以基于NACA 0018翼型的锯齿尾缘仿生叶片为研究对象，采用大涡模拟的方法研究锯齿相对齿宽与相对齿高对锯齿尾缘叶片的气动特性和非定常绕流流场的影响规律和机制.研究表明，尾缘锯齿参数对叶片气动性能的影响是复杂的非线性过程，在一定来流攻角范围内能提高升阻比，但失速提前.如在9.4°~14.8°来流攻角范围内，不同相对齿宽系列叶片的升阻比高于原始叶片，升阻比与锯齿相对齿宽之间没有线性关系.研究还表明，锯齿尾缘能延迟边界层分离，加速尾迹的流动掺混和能量扩散，改变非定常涡结构和涡脱落频率.相对齿高的变化对非定常流动特性的影响更为显著.尾缘锯齿诱导的二次湍流射流和吸力面侧反向涡对改变了原始叶片的绕翼环量，进而影响锯齿尾缘叶片的气动特性和绕流流场特性.      

介质阻挡放电等离子体对翼型流动分离控制的实验研究

在低速开口风洞中进行了等离子体激励器对NACA0015翼型流动分离控制的实验研究。采用PIV技术,对翼型绕流流场进行了测量,显示了施加等离子体激励后流场的变化。通过五分量天平对升力和阻力的测量,研究了激励电压和激励频率对翼型流动分离控制的规律。研究表明,低风速下在翼型前缘施加等离子体激励,能够有效地控制翼型流动分离,在来流为20m／s时,最大升力系数增加11％,失速迎角增加6°;在给定的流动状态下,激励电压和激励频率存在一个阈值,不同迎角下该阈值不同,迎角越大,分离越严重,对激励强度的要求也越高。     

基于RANS-LES混合方法的翼型大迎角非定常分离流动研究

使用雷诺平均Navier-Stokes方程-大涡模拟（RANS-LES）混合方法中的延迟分离涡模拟（DDES）方法,模拟了NACA 0015翼型在大迎角下的静态绕流和强迫振荡运动并与实验值进行了比较。在大迎角静态翼型大分离流动模拟中,DDES方法捕获了非定常RANS计算未能获得的机翼背风面的涡脱落现象。在所采用的RANS和DDES模型中,基于剪切应力输运（SST）湍流模型的SST-DDES混合方法给出的时均压力系数分布与实验吻合得最好。在大迎角强迫振荡翼型绕流模拟中,DDES方法得到的非定常气动载荷与实验值吻合得很好,正确地反映了最大迎角处阻力和俯仰力矩的阶跃性突变,而非定常RANS计算则给出了完全错误的趋势。