超临界翼型低雷诺数流动分析及优化设计

以高空长航时无人机(UAV)翼型研究为背景,对超临界RAE2822翼型在高空高亚声速下的低雷诺数气动特性进行了数值模拟及优化设计研究。采用求解雷诺平均N-S方程的有限体积法,对典型低雷诺数下RAE2822翼型绕流进行数值模拟,验证了SST k-ω湍流模型的可靠性和准确性;基于不同高度不同雷诺数下RAE2822翼型的计算气动力对比分析,研究了高度增大所带来的低雷诺数效应;通过对低雷诺数下超临界翼型表面流场结构及流动机理的详细分析,提出了一种弱化激波的翼型设计思想,并通过优化算例验证了该思想的可行性。     

翼型低雷诺数层流分离现象随雷诺数的演化特征

层流分离现象是翼型低雷诺数条件下出现的典型流场特征。层流分离流动中包含流动分离、转捩、再附等非定常流动结构，层流分离流动的形成与演化会对翼型气动特性产生恶化作用。采用大涡模拟（LES）方法对低雷诺数范围内不同雷诺数下的翼型层流分离流动开展精细数值模拟，研究了雷诺数对翼型气动特性的影响规律及作用机理。LES方法采用隐式亚格子模型，基于结构化拼接网格，对流项离散和时间推进方法分别采用AUSM⁺格式以及双时间步方法。验证算例计算结果表明数值模拟方法的正确性及可靠性，雷诺数对翼型气动特性具有显著影响。随雷诺数降低，时均分离泡外形增大、位置后移，平均阻力系数增大，特别是在较低雷诺数下，翼型升阻力系数随时间出现振荡现象。进一步研究表明，造成不同时均分离泡形态和气动特性的原因在于翼型上表面分离剪切层的失稳与转捩特征。随雷诺数降低，流动黏性增大，导致分离剪切层速度梯度减小，流动发生转捩及再附位置后移，直至翼型表面不再发生转捩和再附。     

基于压力分布的开孔翼型气动特性数值模拟研究

以NACA0012为基准翼型,生成C型结构网格,采用SA湍流模型求解NS方程,获得基准翼型压力分布数值模拟结果。找出翼型上表面前缘低压区域和后缘高压区域,在两个区域进行开孔,让后缘高压区域气流通过导流管引流到前缘区域,加速翼型上表面气流流动,增大翼型速度环量。通过基准翼型与开孔翼型气动特性对比分析可知:小迎角下,基于压力分布的开孔翼型阻力系数略有增大,但能有效提高翼型的升力系数。     

超临界自然层流机翼设计及基于TSP技术的边界层转捩风洞试验

为了实现绿色航空节能减排的目标，层流设计技术成为飞行器设计者的研究热点。对于跨声速客机而言，超临界自然层流机翼设计技术将显著减小飞行阻力，提升气动性能，减少燃油消耗和污染物排放。首先，基于高精度边界层转捩预测技术耦合翼型优化设计系统，实现超临界自然层流翼型设计；经过合理的翼型配置，形成超临界自然层流机翼。转捩数值模拟分析结果表明，超临界自然层流机翼的层流流动特性良好。然后，以比例为1：10.4的试验模型在荷兰高速低湍流度风洞进行边界层转捩风洞试验，使用温度敏感材料涂层（TSP）技术拍照获得机翼表面在不同马赫数、雷诺数和迎角工况下的层流-湍流分布。最后，通过超临界自然层流机翼边界层转捩试验结果，探讨了该类型机翼的转捩特性随来流参数的变化规律，总结了超临界自然层流机翼设计的关键因素。此外，该模型也可用来验证边界层转捩预测技术在超临界、高雷诺数工况下的预测精度。     

NACA0012翼型低雷诺数绕流的实验研究

通过水槽氧气泡流动显示和PIV测速实验研究了NACA0012翼型在雷诺数为8200时的流动特性,重点炎注了翼型绕流结构随迎角的变化。研究发脱：分离点和分离翦切层形成旋涡的位置随迎角的增大而向上游移动,同时翼型上表面流动分离后形成的回流区尺寸随着翼利迎角的增加而增大。当流动再附于翼型上表面时,在再附点附近能够观测到展向涡的三维演化过程,并能观测到展向涡的局部配对现象。     

过失速薄翼增升流动控制方法

在低速风洞中对薄翼升力特性进行了试验研究。采用吸气装置在翼型上表面进行流动控制,利用外式天平测量翼型气动力,利用PIV测试设备获取翼型表面流场。试验来流速度为5m/s,雷诺数6.7×104。研究结果表明:过失速条件下,合适的吸气控制可以使翼型失速迎角延迟近7°,最大升力系数可增大近一倍;在翼型前缘进行吸气流动控制时,较小吸气流量即可延缓翼型失速,但当吸气流量达到一定值时后,随吸气流量增大翼型升力基本保持不变;流动控制参数存在优化空间,当吸气相对位置位于x/c=0.4附近时,吸气流量小于3%即可产生较大的升力增量。     

非结构网格在积冰数值模拟中的应用

基于非结构网格对翼型表面的积冰进行了数值模拟研究.通过求解气体欧拉方程得到翼型的绕流流场,使用穿透型壁面边界模拟水滴对翼型表面的撞击,采用迭代法对水滴控制方程进行数值求解,得到翼型表面的水滴收集特性;采用不同的积冰冻结模型模拟霜冰和明冰的冻结过程;利用冰层时间推进法模拟积冰过程,基于积冰法向生长假设生成积冰外形.对NA...     

充气式翼型绕流大涡模拟分析

针对多梁式充气机翼阿翼型,采用大涡模拟方法对其表面的复杂流动进行了数值仿真研究,得到了其迎角在0°-20°时的流动情况和气动力特性,并与其他仿真结果和风洞试验结果进行了比较分析。结果表明：采用大涡模拟方法得到的计算结果与风洞试验结果较为吻合,可以更好地模拟充气式翼型的表面流动情况。     

实体鼓包对超临界翼型的流动控制和减阻研究

采用数值模拟的手段,研究了实体鼓包对超临界翼型的流动控制和减阻作用。超临界机翼在非设计状态时由于波阻增加导致总阻力增大;在激波的波脚位置有效地使用实体鼓包,可以减小激波阻力;在中高升力系数情况下,使用实体鼓包可提高升阻比。     

基于OpenFOAM的超临界压力下低温射流大涡模拟

为了深入认识低温射流在超临界压力下的流动特性,将真实流体热物性模型写入开源CFD程序OpenFOAM平台中,开发出了用于模拟超临界压力下流动过程的均相求解器。针对液氮跨临界、超临界射流进行大涡模拟,研究了环境条件变化对于低温射流流动特性的影响。结果表明：采用的真实流体模型能够在广泛的温度与压力范围内准确计算流体的热物理性质,在此基础上开发的均相求解器能够准确描述超临界条件下的低温射流,其流动特性主要表现为稠密液体与环境气体之间的剪切层不稳定以及湍流混合,此类变密度射流在湍流充分发展区域同常密度射流一样具有自相似特性;超临界射流相比于跨临界射流具有更高的混合效率以及更短的液核穿透长度;环境温度升高以及环境压力降低均使得低温射流与周边气体之间的密度梯度增加,对剪切层中不稳定波动的发展起到抑制作用,导致了射流液核长度的增加以及扩张角度的减小。     

翼型结冰过程数值模拟验算与分析

应用FENSAP-ICE结冰计算软件,对NACA0012翼型进行了流动特性、水滴撞击特性以及冰型生成过程的计算;同时,对结冰前后的翼型进行气动力特性计算对比分析,其中包括升力特性对比、阻力特性对比、流场细节分析以及压力系数分布对比。计算结果表明：翼型前缘结冰后,导致翼型前缘气流提前分离,最大升力系数、失速攻角大幅减小,...     

基于遗传算法的压气机叶型优化设计

建立了一套将遗传算法和流场数值模拟技术相结合的压气机叶型优化设计系统,针对进口高亚音的压气机叶型进行优化。通过准三维叶栅通道计算程序进行流场数值模拟,评估叶型气动性能。选取了与叶型形状和气动性能密切相关的三个变量作为优化参数,按照可控扩散叶型的设计标准构造了优化目标函数。数值计算结果表明优化叶型的气动性能有了明显提高。对优化前后的叶型进行了叶栅吹风实验对比,进一步验证了优化效果。     

重型燃气轮机压气机高雷诺数前转捩叶型设计

针对重型燃气轮机压气机雷诺数高而导致的转捩位置前移,开发了一种比可控扩散叶型(CDA)损失更小、工作范围更宽的前转捩叶型.采用正问题优化设计方法,将叶型几何参数化、叶片到叶片流场分析与遗传算法相结合,实现了叶型的自动优化.优化目标综合权衡了叶型损失和攻角范围,为减少优化变量的数目,应用了一种特别的叶型几何模型,将厚度分布与中弧线之间进行了一定的关联.优化得到的前转捩叶型的主要特征是吸力面速度峰值的位置前移至距前缘约10%弦长处,叶型中后部的速度变化更为平缓.最后根据优化结果总结了前转捩叶型的设计规律.      

某型燃气涡轮起动机涡轮级与喷管一体化三维流场数值计算

通过三维N-S方程,对某型燃气涡轮起动机涡轮级和喷管内部的三维流场进行了数值模拟计算,分析了通道内的涡系产生和流动分离,揭示了该涡轮级和喷管内部流场的物理特征,为该起动机的改进设计提供了理论依据,并为整机的数值模拟奠定了基础。     

凝结相变对低温风洞雷诺数试验的影响

为了研究凝结相变对低温风洞雷诺数试验能力的影响,基于Fluent软件建立了气-液两相凝结流动模型,相变模型采用考虑非等温效应修正的经典成核理论和Gyarmathy液滴生长理论;针对不同来流压力和不同试验雷诺数工况条件,对氮气绕流NACA 0012翼型进行了数值模拟。模拟结果表明:随着来流温度的降低,翼型附近区域气体膨胀越过气态饱和线并达到过冷状态,进一步降低来流温度则会在宏观层面上观测到凝结相变对当地流场的改变,与无凝结相变的流场相比,释放的潜热加热气流导致马赫数降低及压力系数的改变;在确保不破坏翼型气动性能试验的前提下,充分利用气体过冷区域来降低来流压力以此减少驱动功率和液氮喷入量是可行的,或者保持来流压力不变提升低温风洞的试验雷诺数。     

基于N-S方程的地面效应模拟

模拟了低速情况下二维翼型的地面效应问题，其中翼面流动利用一种改进的、基于不可压N－S方程的离散涡方法来模拟，地面模拟则利用一系列连续分布的涡板块来实现。利用上述模型研究了地面对翼型表面流动和压力分布的影响，并对有地面存在条件下的尾涡发展规律进行了讨论。计算结果表明，该方法能够有效地模拟低速情况下翼型的近地面流动，具有一定的工程应用价值。     

变来流速度下旋翼翼型非定常气动特性分析

为分析变来流速度状态下的旋翼翼型气动特性,提出了利用翼型平移来模拟来流速度变化的数值方法.在此方法基础上,采用基于隐式LU-SGS(lower upper symmetric Gauss-Seidal)方法的非定常雷诺平均N-S(Navier-Stokes)(RANS)方程,模拟了SC1095旋翼翼型在定迎角 变来流速度及变迎角 变来流速度状态下的非定常气动特性.通过对比分析发现:翼型在变速度-定迎角状态下会表现出明显的非定常现象,产生了前缘分离涡,气动特性会出现明显的迟滞效应及波动现象,脉动速度越大,非定常效果越明显.并且基准速度越大,翼型气动特性的峰值越大；翼型迎角越大,非定常涡出现的也越早.考虑直升机旋翼翼型实际工作环境,在变速度-动态失速状态下,翼型最大迎角处的气动力会得到一定程度的削弱,在小迎角下的气动力得到一定程度的增强,且脉动速度越大,翼型的非定常特性也越强.      

低速两段翼型的数值模拟

在多段翼型风洞实验中,很难观察翼型不同迎角下主翼、襟翼上的绕流及缝道流动。使用了O-H型混合的结构化网格,采用k-ε二方程湍流模型求解可压的N—S方程来模拟二元风洞中两段翼型的流场。计算结果与实验进行了比较,结果表明：计算与实验结果吻合良好,说明本方法可以较好地模拟两段翼型的绕流,并能很好地显示流场的变化情况。另外,从结果中可以看出：随着迎角的变化,主翼尾流和缝道间的流动对襟翼附面层分离有很大的影响作用。     

零质量射流推迟翼型失速的数值模拟

 通过在具有分离流动的翼型表面局部地引入零质量射流干扰，可以较小的能量代价实现对翼型绕流的全局性主动控制，从而达到抑制分离、推迟翼型失速的目的。本文发展了零质量射流与翼型绕流干扰的数值分析方法。在格心有限体积法的基础上，空间离散采用AUSM+up格式，时间推进运用含牛顿型LU-SGS子迭代的全隐式双时间法，且引入了预处理方法和多重网格法等技术。以TAU0015翼型失速主动控制为算例，通过在翼型上表面前缘附近引入零质量射流控制，使翼型失速特性得到改善。将计算结果与实验结果做了比较，验证了本文发展的数值模拟方法及程序的正确性。还以NACA0015翼型为例研究了射流动量系数、射流频率、射流偏角等关键参数对失速控制效果的影响，得出一些关于如何选择控制参数的结论。     

Application of Active Flow Control Technique for Gust Load Alleviation

A new gust load alleviation technique is presented in this paper based on active flow control. Numerical studies are conducted to investigate the beneficial effects on the aerodynamic characteristics of the quasi "Global Hawk" airfoil using arrays of jets during the gust process. Based on unsteady Navier-Stokes equations, the grid-velocity method is introduced to simulate the gust influence, and dynamic response in vertical gust flow perturbation is investigated for the airfoil as well. An unsteady surface transpiration boundary condition is enforced over a user specified portion of the airfoil’s surface to emulate the time dependent velocity boundary conditions. Firstly, after applying this method to simulate typical NACA0006 airfoil gust response to a step change in the angle of attack, it shows that the indicial responses of the airfoil make good agreement with the exact theoretical values and the calculated values in references. Furthermore, gust response characteristic for the quasi "Global Hawk" airfoil is analyzed. Five kinds of flow control techniques are introduced as steady blowing, steady suction, unsteady blowing, unsteady suction and synthetic jets. The physical analysis of the influence on the effects of gust load alleviation is proposed to provide some guidelines for practice. Numerical results have indicated that active flow control technique,as a new technology of gust load alleviation, can affect and suppress the fluid disturbances caused by gust so as to achieve the purpose of gust load alleviation.