结冰区流场特征一致的混合缩比翼型设计

提出了一种基于翼型前缘范围流场相似的多控制点混合缩比翼型优化设计方法:以前缘局部范围内的空间流场特征相同为目标,采用多目标遗传优化算法结合翼型参数化、多控制点、网格自适应技术及CFD,进行了不同工况下的混合缩比翼型设计,分别通过低湍流度风洞及结冰风洞对翼型的气动性能及积冰过程进行了实验.设计及实验结果表明:该混合缩比翼...     

不同雷诺数下翼型气动特性及层流分离现象演化

低雷诺数下空气黏性效应突出,翼型表面普遍存在层流分离现象,相比常规雷诺数情况气动特性显著恶化。采用带预处理的Roe方法求解非定常可压缩Navier-Stokes方程的数值模拟技术和低雷诺数低湍流度风洞油流显示试验技术,对FX63-137翼型不同雷诺数下气动特性和流动结构展开深入研究。通过风洞油流显示试验可以清晰获得低雷诺数层流分离流动的两道油流汇集线。数值模拟结果表明其分别为时均化主分离线和二次分离线,两种结果定性定量均吻合较好,证明了本文的研究方法有效可靠;雷诺数从500 000降至20 000,翼型气动特性和层流分离流动结构均发生显著的变化,伴随阻力系数剧增和升力系数剧降,时均化流动结构从附体至出现经典的长层流分离泡,并最终演化为后缘层流分离泡,相应的两种分离泡的非定常流动结构也存在显著差异;对于阻力系数和升力系数而言,存在不同的临界雷诺数,因为导致阻力系数剧增的机理在于经典长层流分离泡的产生使翼型压差阻力大增,而造成升力系数剧降的主要原因在于后缘层流分离泡使得等效翼型后部弯度减小;非定常结果显示正是由于翼型表面漩涡周期性的生成与脱落,才造成了低雷诺数下升力系数的周期性波动。翼型上表面主分离涡即将脱落时,流线在后缘附近再附,升力系数达到峰值;而当流体从下表面向上卷起二次分离涡时,尾部流线大尺度分离,升力系数降至谷值。     

涡襟翼在不同雷诺数下的控制分离特性研究

受鸟类抬起羽毛控制分离流的启发,涡襟翼成为翼型大迎角分离流的控制措施之一。采用数值模拟方法研究不同雷诺数下涡襟翼在控制翼型大迎角分离流动时的气动特性及其物理机制。结果表明：涡襟翼在低雷诺数下能够极大地改善翼型的大迎角升力特性,其物理机理是涡襟翼将翼型主分离涡的涡心位置控制在离翼型更近的区域,且涡心位置的涡量得到大幅提升,使得涡心附近的低压特性影响到翼型上表面,而且涡襟翼能够将翼型上方前区的低压与下游的高压隔开;但是在高雷诺数（对应常规飞机雷诺数）下涡襟翼改善翼型大迎角气动特性的效果远不如低雷诺数情况,由此解释了为什么鸟类能够通过羽毛抬起提高升力特性,而常规飞机的涡襟翼只能作为阻力板使用的原因。     

大型煤气柜风荷载的风洞试验及数值模拟

对某大型煤气柜进行了风洞测压试验及风压数值模拟。分析了试验模型表面风压分布及其脉动特性,并同数值计算结果、规范条文中类似断面结构的风压分布作对比。结果表明：风洞试验中由于结构表面分布的工字钢及表面粗糙度的处理,雷诺数效应对表面风压分布影响并不明显,但对表面绕流场分离区的风压值有一定影响。结构的均方根升力、阻力系数在频域表现为宽带谱;采用基于雷诺平均的RNGk-ε湍流模型能较准确地模拟表面平均风压分布,其计算结果同样可为结构抗风设计提供参考;在不同量级雷诺数下数值模拟得到的平均风压分布能反映出雷诺数效应的影响。     

Aerodynamic performance enhancement for flapping airfoils by co-flow jet

Introducing active flow control into the design of flapping wing is an effective way to enhance its aerodynamic performance. In this paper, a novel active flow control technology called Co-Flow Jet (CFJ) is applied to flapping airfoils. The effect of CFJ on aerodynamic performance of flapping airfoils at low Reynolds number is numerically investigated using Unsteady Reynolds Averaged Navier-Stokes (URANS) simulation with Spalart-Allmaras (SA) turbulence model. Numerical methods are validated by a NACA6415-based CFJ airfoil case and a S809 pitching airfoil case. Then NACA6415 baseline airfoil and NACA6415-based CFJ airfoil with jet-off and jet-on are simulated in flapping motion, with Reynolds number 70,000 and reduced frequency 0.2. As a result, CFJ airfoils with jet-on generally have better lift and thrust characteristics than baseline airfoils and jet-off airfoil when C_μ is greater than 0.04, which results from the CFJ effect of reducing flow separation by injecting high-energy fluid into boundary layer. Besides, typical kinematic and geometric parameters, including the reduced frequency and the positions of the suction and injection slot, are systematically studied to figure out their influence on aerodynamic performance of the CFJ airfoil. And a variable C_μ jet control strategy is proposed to further improve effective propulsive efficiency. Compared with using constant C_μ, an increase of effective propulsive efficiency by 22.6% has been achieved by using prescribed variable C_μ for NACA6415-based CFJ airfoil at frequency 0.2. This study may provide some guidance to performance enhancement for Flapping wing Micro Air Vehicles (FMAV).     

考虑间隙内流的二元机翼跨声速气动力分析

基于Navior-Stokes (N-S)方程和Spalart-Allmaras湍流模型,建立了含间隙二元机翼的跨声速绕流场计算模型和间隙内流计算模型.通过FLUENT数值仿真,分析了二元机翼的激波压力分布和马赫数分布,并与风洞试验 结果进行了对比验证.数值结果表明,连通波后高压区和波前低压区的间隙内流会提高机翼绕流场中激波前后的速度梯度和压力梯度.局部间隙内流会对机翼结构带来“向外吸”和“向前顶”的气动载荷作用.     

NPU-WA系列风力机翼型设计与风洞实验

针对兆瓦级大型风力机,研究发展了以具有更优良高雷诺数和高升力气动性能为特点的NPU-WA翼型族,风洞实验表明,该翼型族达到了在高雷诺数、高升力条件下实现高升阻比和外侧翼型对粗糙度不敏感的主要设计要求,为我国自主研发大型风力机提供了可以实际使用的翼型几何数据和雷诺数范围内1.0×106～5.0×106的风洞实验数据。     

低雷诺数翼型局部振动非定常气动特性

针对低雷诺数翼型特殊的气动特性,采用基于动网格的非定常数值模拟方法,研究翼型表面不同弦向位置的局部蒙皮以不同频率及振幅振动时对低雷诺数翼型气动特性及流场结构的影响,揭示蒙皮振动增升减阻的机理。研究表明,在低雷诺数条件下局部蒙皮振动可有效提高翼型气动特性,与刚性翼型相比蒙皮局部振动可使翼型升力系数提高,阻力系数降低,升阻比提高。振动位置对翼型气动特性及流场结构有显著的影响,振动表面位于翼型前缘附近或位于层流分离泡中心时可有效控制翼型层流分离,从而提高翼型气动特性。振动频率对翼型表面层流分离及转捩位置均有显著的影响,随着振动频率增加,翼型气动特性出现最优值。与刚性翼型相比,表面振动使翼型转捩位置略向上游移动,摩擦阻力增加,但振动使等效翼型相对厚度减小,压差阻力明显减小。在小幅振动范围内,随着振幅增加,流场非定常特性更加显著,翼型升阻比增加。     

高升力下二维构型阻力准确测量研究

翼型风洞试验阻力测量常使用尾迹流场测量积分求取阻力的方法,但各积分公式均建立在一定的假设基础上,有一定适用范围。在多段翼型流场N-S方程数值模拟和风洞试验的基础上,研究高升力情况下低速风洞阻力精确测量技术。通过N-S方程数值模拟求解多段翼型绕流场,分析尾迹流场的特点和常规风洞试验阻力计算公式推导时所作假设,提出新的更为准确的型阻计算公式;利用多段翼型绕流的数值模拟结果,积分表面压力和摩擦力求得翼型的气动特性,并利用计算得到的尾迹流场信息按照常规和新提出的风洞试验型阻计算公式计算阻力,将三者进行比较,检验提出的新型阻计算公式的准确性;通过风洞试验检验数值模拟得到的流场特点和新型阻计算公式。研究表明：在高升力条件下,传统型阻计算公式有很大的局限性,必须进行改进;提出的考虑尾迹区流动特点的新型阻计算公式能够得到更准确的阻力值。     

飞虫粘附翼型对翼型气动特性影响的实验研究

在西北工业大学NF-3低速风洞进行了飞虫粘附翼型对翼型气动性能影响的风洞实验研究。结果表明：模型表面粘附的飞虫数量累计到一定程度时,将会导致翼型表面50%区域以上的面积发生分离,引起翼型失速,并且实验前模型表面的飞虫数量会改变翼型的失速迎角,因此实验前必须将模型擦拭干净。鉴于飞虫粘着数量的不确定性,对于翼型在复杂环境下使用时其气动性能的变化需要加以关注。     

厚覆冰导线升力突变机理及数值模拟研究

针对新月形厚覆冰导线的升力系数在风攻角15°附近存在突变的问题,分别采用基于k-ωSST湍流模型的雷诺时均法和大涡模拟(LES)的数值方法对新月形厚覆冰导线在风攻角10°～20°范围进行了模拟。通过对比两种数值方法计算得到的覆冰导线气动力系数、流场结构和表面风压,发现LES方法能够更好地捕捉新月形覆冰导线表面的小尺度涡结构,得到的覆冰导线气动力参数计算结果与风洞试验数据高度吻合;而k-ωSST湍流模型难以模拟壁面上小尺度涡,捕捉不到升力系数的突变。根据覆冰导线不同壁面区域的压力分布,发现上侧壁面处的涡结构影响整体流场,并在下侧壁面曲率、来流夹角和壁面切线方向共同作用下导致升力系数突变。LES的气动力参数模拟结果可为覆冰导线防舞提供参考。     

低雷诺数仿生分离流翼型气动特性初探

为了探索适合低雷诺数微型飞行器的翼型形式,基于对自然界鸟类和昆虫滑翔飞行时翅膀形状的观察,设计出一种由前缘削尖平板和后缘圆弧翼型组合而成的仿生分离流翼型。数值研究结果表明,气流在削尖平板的前缘点强制分离,形成大范围低压分离流动,随后在后部圆弧翼上表面再附形成稳定低压涡流区,从而实现较高的气动效率和较强的抵抗大气湍流的能力。上削尖平板可以使流动分离点固定在削尖点。相对于单独平板,仿生分离流翼型的升力系数有大幅提高,迎角为4°时提高了112%。此外,仿生分离流翼型可以在较宽的迎角范围内（4°~20°）保持高升力,但是迎角增加,阻力也快速增大,因此小迎角情况下（小于4°）气动效率更优。      

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

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

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

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

低速两段翼型的数值模拟

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

Effects of wing flexibility on aerodynamic performance of an aircraft model

The low-speed wind tunnel experiment is carried out on a simplified aircraft model to explore the influence of wing flexibility on the aircraft aerodynamic performance. The investigation involves the measurements of force, membrane deformation and velocity field at Reynolds number of 5.4 × 10⁴–1.1 × 10⁵. In the lift curves, two peaks are observed. The first peak, corresponding to the stall, is sensitive to the wing flexibility much more than the second peak, which nearly keeps constant. For the optimal case, in comparison with the rigid wing model, the delayed stall of nearly 5° is achieved, and the relative lift increment is about 90%. It is revealed that the lift enhanced region corresponds to the larger deformation and stronger vibration, which leads to stronger flow mixing near the flexible wing surface. Thereby, the leading-edge separation is suppressed, and the aerodynamic performance is improved significantly. Furthermore, the effects of sweep angle and Reynolds number on the aerodynamic characteristics of flexible wing are also presented.     

Aerodynamic periodicity of transient aerodynamic forces of flexible plunging airfoils

For flapping micro-air vehicles, periodical aerodynamics in every plunging period is always desired in the design of aircraft controlling system. During the periodic plunges of rigid airfoil, either periodical or non-periodical aerodynamic forces can be generated. Real airfoils are usually flexible so that the effect of flexibility on unsteady characteristics of aerodynamic periodicity should be considered. In this study, a fluid structure interaction analysis is employed to systematically investigate the aerodynamic periodicity of a periodically plunging flexible airfoil. The influences of several dimensionless kinematic and structural parameters on aerodynamic periodicity are investigated. The results show that the aerodynamic periodicity of plunging flexible airfoil qualitatively resembles that of rigid airfoil as the mechanisms contributing to aerodynamic periodicity are kept for airfoil with various flexibilities. However, aerodynamic periodicity at small plunging amplitude and high Reynolds number can be improved as airfoil becomes medium flexible but worsened as airfoil is very flexible. The reason is that the deformation of the airfoil changes the strength and structure of the nearby vortices, thereby leading to weak vortex–vortex and vortex–airfoil interactions. The extent of improvement in aerodynamic periodicity is also found significantly influenced by incoming velocity.     

翼型加装格尼襟翼的低雷诺数气动特性实验研究

为了研究低雷诺数下格尼襟翼对翼型气动特性的影响,通过风洞试验研究了Eppler387翼型加装0.5%～5.0%弦长高度格尼襟翼后的气动特性变化,试验雷诺数1.49×105～2.31×105。试验结果表明：低雷诺数下Eppler387翼型加装格尼襟翼后,升力系数和力矩系数明显增大,襟翼高度大于2%弦长时阻力系数显著增大。格尼襟翼在高升力系数下能够起到增大升阻比的作用,适用于微小型飞行器工作在大载荷状态,而0.5%弦长高度的襟翼还能够兼顾中小升力系数下的气动效率,同样适合于微小型飞行器在巡航状态使用。与原翼型相比,加装襟翼后最大升阻比对应的迎角提前,随襟翼高度的增加,翼型升阻比曲线峰值变得不再突出。     

战斗机弹射救生系统数值模拟研究

对弹射救生系统这种几何外形复杂、不规则非流线型钝体,以Jameson有限体积法为基础发展了一套基于高质量混合网格的Euler方程解算器,在Euler方程解算器基础上,采用S-A湍流模型发展了一套同样基于混合网格的N-S方程解算器,初步建立了适合弹射救生系统这类几何外形复杂物体的数值模拟技术平台,并采用多种加速收敛措施和合适的耗散项模型使得所发展的解算器在软件实现上鲁棒性很强。通过选取合适的计算起始迎角,对弹射救生系统成功地进行了大迎角和大侧滑角绕流计算,获得了与风洞试验比较吻合的计算结果,并且可以清楚地观察流场中马赫数分布、弹射救生系统各部位压力分布等流场细节,对弹射救生系统空气动力特性分析具有重要的指导意义。该计算平台可作为风洞试验的合理补充、加强和完善。     

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.