基于涡粘性模型的翼型湍流流场熵产率计算

利用有限体积法实现了基于非正交同位网格的SIMPLE算法。基于熵分析方法,采用涡粘性模型求解湍流熵产方程,系统研究了湍流模型对二维翼型绕流流场熵产率的影响。通过计算NACA0012翼型在来流雷诺数为2.88×106时,0°攻角~16.5°攻角范围内的翼型表面压力系数分布和升阻力特性,验证了算法及程序的正确性。结果表明,选择不同湍流模型时,翼型流场熵产的计算结果存在差异,湍流耗散是引起流场熵产的主要原因;翼型流场的熵产主要发生在翼型前缘区、壁面边界层和翼型尾流区域,流场熵产率与翼型阻力系数线性相关;当产生分离涡时,粘性耗散引起的熵产下降。     

圆形坚冰影响翼型气动性能的数值分析

本文用FLUENT软件模拟了结冰后NACA 0012翼型周围流场的变化,并与结冰前NACA 0012翼型的气动性能进行了对比.工作中首先以未结冰的NACA 0012翼型(干净翼型)为标准模型进行了数值验证计算,再以经过检验的方法计算结冰模型,并与结冰风洞试验数据进行了对比.本文计算攻角为0°～20°,温度为250.37K,雷诺数为2,400,000,冰型为圆形坚冰.通过对比升力阻力性能,发现与干净翼型相比,结冰翼型的最大升力系数大约减少了50%,阻力系数增加了约65%,失速攻角降低了4°.结冰后翼型提前失速是造成气动性能恶化的主要原因.     

智能物面对非定常分离流的最优自适应控制

以低雷诺数二维大攻角翼型绕流为研究对象, 将非定常动边界计算流体力学方法与 最优控制方法有机结合, 研究二维不可压非定常流智能物面最优自适应流 动控制的理论与算法, 并将其用于固定攻角和俯仰振荡翼型绕流. 结果表明: 在给定合适的最优控制目标函数下, 智能物面可最优地实时改变形状, 得到能显著提高翼型性能的最优翼型. 最优翼型在非设计工况下的气动性能也比对照翼型有 所提高.     

基于涡粘性模型的翼型湍流流场熵产率计算

利用有限体积法实现了基于非正交同位网格的SIMPLE算法。基于熵分析方法,采用涡粘性模型求解湍流熵产方程,系统研究了湍流模型对二维翼型绕流流场熵产率的影响。通过计算NACA0012翼型在来流雷诺数为2.88×10⁶时,0°攻角~16.5°攻角范围内的翼型表面压力系数分布和升阻力特性,验证了算法及程序的正确性。结果表明,选择不同湍流模型时,翼型流场熵产的计算结果存在差异,湍流耗散是引起流场熵产的主要原因;翼型流场的熵产主要发生在翼型前缘区、壁面边界层和翼型尾流区域,流场熵产率与翼型阻力系数线性相关;当产生分离涡时,粘性耗散引起的熵产下降。     

新型等离子体激励器对翼型的减阻效果

论述了在西北工业大学低湍流度风洞中采用新型等离子激励器对NACA0015翼型的减阻实验.实验风速为35m/s,攻角范围取0°~20°.并参照压力分布和总压分布实验结果对减阻效果进行了对比分析.本文还进行了有关等离子体激励抑制翼型流动分离的数值模拟研究,基于等离子体激励器的简化模型将体积力以源项方式引入到N-S方程中求解,得到激励器工作时的流场分布.结果表明在新型等离子体激励器开启后:在小攻角范围内,尾耙的总压分布曲线与坐标轴的纵轴(尾耙高度轴)所围面积变化不大;当攻角≥12°时,尾耙的总压分布曲线与坐标轴的纵轴(尾耙高度轴)所围面积明显减小.从而说明该新型等离子体激励器能够有效地减少翼型的阻力.     

随机湍流工况低雷诺数风力机翼型优化研究

为改善小型风力机随机湍流工况适应性,以NACA0012翼型为研究对象,采用非嵌入式概率配置点法,获得随机湍流工况下小型风力机叶片翼型运行攻角分布规律;在气动优化中耦合层流分离预测,基于Transition SST模型、拉丁超立方试验设计、Kriging模型和带精英策略非支配排序遗传算法NSGA-II进行高湍流低雷诺数风力机翼型气动优化。结果表明,优化翼型叶片平均风能捕获效率分别提高3.01%和4.76%,标准差分别降低4.76%和14.93%,优化翼型湍流适应性增强。该方法将翼型设计与湍流风况相匹配,为湍流工况低雷诺数翼型及小型风力机设计提供参考。     

绕振荡水翼流动及其转捩特性的数值计算研究简

通过对比标准k-ω SST 湍流模型和基于标准k-ω SST 湍流模型修正的γ-Re_θ 转捩湍流模型对绕振荡NACA66 水翼流动的数值计算结果与实验结果,对水翼振荡过程的水动力特性和流场结构变化进行了分析研究. 结果表明：与标准k-ω SST 湍流模型的数值计算结果相比,基于标准k-ω SST 湍流模型修正的γ-Re_θ 转捩湍流模型能有效预测绕振荡翼型流场结构和水动力特性,捕捉流场边界层发生的流动分离和转捩现象;绕振荡水翼的流动过程可分为5 个特征阶段,当来流攻角较小时,在水翼前缘发生层流向湍流的转捩现象,水翼动力特征曲线出现变化拐点;随着来流攻角的增大,顺时针尾缘涡逐渐形成并向水翼前缘发展;当攻角较大时,前缘涡分离导致动力失速,水翼的动力特征曲线出现大幅波动;水翼处于顺时针向下旋转阶段,绕水翼的流动状态逐渐由湍流过渡为层流.     

低雷诺数下柔性翼型气动性能分析

基于流固耦合方法对吸力面5%至95%弦长处为三段柔性结构的NACA0012翼型绕流进行了数值模拟,研究了不同弹性模量下柔性翼型的气动性能和结构响应.结果表明:在大攻角下,翼面变形影响着翼型表面的非定常流场,起到延缓失速和提高升力的作用;失速后柔性翼的升力系数下降得较为缓慢,且柔性越大,升力系数下降得越平缓;适当减小弹性模量能够提高翼型的气动性能,然而弹性模量过小反而不利于翼型气动性能的提升,并且翼面会产生大幅度的振动.     

利用非结构化网格方法对翼型绕流的数值研究

利用同位非结构化网格上的压力加权修正算法 ,对翼型湍流绕流进行了数值分析。详细地给出了一孤立翼型在不同攻角下的分离流结构及翼型表面压力分布 ,为了显示非结构化网格方法在求解多连通流动区域的优越性 ,对双翼型绕流进行了数值计算。在数值分析中 ,对阵面推进法进行改进来生成三角形网格 ,采用有限控制体方法直接在物理空间中的非结构化网格单元上离散 Navier- Stokes方程及 k- ε方程 ,形成的代数方程组通过预条件矩阵共轭梯度平方法求解。计算结果表明 :当流动为附着流时 ,计算结果与实验值吻合程度令人相当满意 ;而在分离区内 ,计算结果与实验值存在一定的误差 ,需对分离区内的湍流模型做进一步的改进。     

对称翼型低雷诺数小攻角升力系数非线性现象研究

采用Rogers发展的三阶Roe格式,求解非定常不可压N-S方程,时间方向为二阶精度双时间步方法, 数值模拟了对称翼型SD8020低雷诺数(Re=40000,100000)条件下,流场层流分离涡结构和升力系数随攻角的变化．同试验比较证明了数值模拟的正确性．通过对数值模拟时均化流场结果的详细分析,发现对称翼型在小雷诺数0°攻角附近出现的层流分离泡,其内部结构和演化规律都不同于经典层流分离泡模型,从而提出了一种后缘层流分离泡模型．并应用该模型对对称翼型小攻角低雷诺数流场特性以及升力系数非线性效应的形成机理进行了研究和解释．     

Numerical simulation of flow over an airfoil in heavy rain via a two-way coupled Eulerian–Lagrangian approach

Airfoil performance degradation in heavy rain has attracted many aeronautical researchers’ eyes. In this work, a two-way momentum coupled Eulerian–Lagrangian approach is developed to study the aerodynamic performance of a NACA 0012 airfoil in heavy rain environment. Scaling laws are implemented for raindrop particles. A random walk dispersion approach is adopted to simulate raindrop dispersion due to turbulence in the airflow. Raindrop impacts, splashback and formed water film are modeled with the use of a thin liquid film model. The steady-state incompressible air flow field and the raindrop trajectory are calculated alternately through a curvilinear body-fitted grid surrounding the airfoil by incorporating an interphase momentum coupling term. Our simulation results of aerodynamic force coefficients agree well with the experimental results and show significant aerodynamic penalties at low angles of attack for the airfoil in heavy rain. An about 3° rain-induced increase in stall angle of attack is predicted. The loss of boundary momentum by raindrop splashback and the effective roughening of the airfoil surface due to an uneven water film are testified to account for the degradation of airfoil aerodynamic efficiency in heavy rain environment.     

Leading edge vortex formation and detachment on a flat plate undergoing simultaneous pitching and plunging motion: Experimental and computational study

This study focuses on the formation and detachment of a leading edge vortex (LEV) appearing on an airfoil when its effective angle of attack is dynamically changed, inducing additional forces and moments on the airfoil. Experimental measurements of the time-resolved velocity field using Particle Image Velocimetry (PIV) are complemented by a computational study using an URANS (Unsteady Reynolds-Averaged Navier–Stokes) framework. In this framework a transition-sensitive Reynolds-stress model of turbulence, proposed by Maduta et al. (2018), which combines the near-wall Reynolds-Stress model by Jakirlic and Maduta (2015) and a phenomenological transition model governing the pre-turbulent kinetic energy by Walters and Cokljat (2008), is employed. Combined pitching and plunging kinematics of the investigated flat plate airfoil enable the effective inflow angle to be arbitrarily prescribed. A qualitative assessment of flow fields and a quantitative comparison of LEV characteristics in terms of its center position and circulation as well as an investigation of the mechanism causing the vortex to stop accumulating circulation revealed close agreement between the experimental and simulation results. Further considerations of the lift contribution from the pressure and suction side of the airfoil to the overall lift indicates that the qualitative lift evolution is reproduced even if the pressure side contribution is neglected. This reveals important characteristics of such airfoil dynamics, which can be exploited in future experimental studies, where direct aerodynamic force and moment measurements are greatly inhibited by dominating inertial forces.     

排翼布局飞行器气动性能的实验研究

通过低速低湍流度风洞实验,研究了利用排翼布局改善充气飞机采用大厚度翼型机翼带来的气动效率偏低问题。首先比较了采用不同厚度翼型的单翼与排式双翼布局的气动特性。在此基础上,为了优化排翼布局的气动特性,研究了给后翼安装偏转角对排翼布局气动特性的影响。同时,基于NACA0030翼型,设计了波纹型外形的充气机翼,比较了此外形下单翼和排翼布局气动性能的差异。实验结果表明,采用排翼布局能够改善采用厚翼型单翼布局的气动性能,而给后翼安装一定偏转角可以进一步提高排翼布局的升力和升阻比。采用波纹外形和光滑外形机翼模型的对比结果表明,波纹外形能够在大迎角时改善充气机翼的失速性能。分析认为,造成这一现象的流动机理是由于波纹型机翼在实验条件下提前由层流转捩为湍流,使失速推迟,流动分离现象有所减弱。     

超临界翼型跨音速特性实验研究

针对新设计的超临界翼型,采用风洞实验方法验证和评估了其气动特性。在增压连续式跨音速风洞(NF-6风洞)开展了超临界翼型跨音速特性的实验研究,验证了该翼型设计的压力分布曲线特点。激波位置和波后压力平台区长度表明设计结果和实验结果基本一致,揭示了超临界翼型跨音速的气动特性;阻力发散马赫数达到期望的设计指标,探讨了雷诺数对该翼型气动特性的影响。最后采用升华法实现了翼型表面流动特性的显示。结果表明转捩点约在16%弦长位置。     

A study of long separation bubble on thick airfoils and its consequent effects

A parametric study has been performed to analyse the flow around the thick-symmetric NACA 0021 airfoil in order to better understand the characteristics and effects of long separation bubbles (LoSBs) that exist on such airfoils at low Reynolds numbers and turbulence intensities. In the article, the prediction capabilities of two recently-developed transition models, the correlation-based γ–Re_θ model and the laminar-kinetic-energy-based κ–κ_L–ω model are assessed. Two-dimensional steady-state simulations indicated that the κ–κ_L–ω model predicted the separation and reattachment process accurately when compared with published experimental work. The model was then used to study the attributes and the effects of LoSBs as a function of the angle of attack, freestream turbulence intensity and Reynolds number. It was observed that LoSBs considerably degrade the aerodynamic performance of airfoils and lead to abrupt stall behaviour. It is, furthermore, illustrated that the presence of the LoSB leads to an induced camber effect on the airfoil that increases as the airfoil angle of attack increases due to the upstream migration of the bubble. An increase in the Reynolds number or turbulence levels leads to a reduction in the bubble extent, considerably improving the airfoil performance and leading to a progressive trailing-edge stall.     

Numerical investigation of transonic flow over deformable airfoil with plunging motion

In this article, the transonic inviscid flow over a deformable airfoil with plunging motion is studied numerically. A finite volume method based on the Roe scheme developed in a generalized coordinate is used along with an arbitrary Lagrangian-Eulerian method and a dynamic mesh algorithm to track the instantaneous position of the airfoil.The effects of different governing parameters such as the phase angle, the deformation amplitude, the initial angle of attack, the flapping frequency, and the Mach number on the unsteady flow field and aerodynamic coefficients are investigated in detail. The results show that maneuverability of the airfoil under various flow conditions is improved by the deformation. In addition, as the oscillation frequency of the airfoil increases, its aerodynamic performance is significantly improved.     

定常吹气对倾转旋翼机向下载荷的影响

为进一步提高倾转旋翼机悬停状态下的有效载重,开展了定常吹气流动控制对向下载荷的影响研究。首先应用延迟脱体涡模拟(DDES)方法对翼型-90°迎角下非定常大范围分离流动结构进行了数值分析;然后分别开展了前缘吹气、后缘吹气降载措施研究,揭示了吹气降载的机理,并对不同吹气口位置和吹气动量系数的影响进行了定量分析,最后开展了前、后缘同时吹气作用下降载数值模拟研究。计算结果表明:前缘最佳吹气位置在翼型的前缘点,而后缘吹气最佳位置位于襟翼弦长的15%处;前缘吹气的降载效果要优于后缘吹气,而且吹气动量系数对向下载荷的影响较小;相对于初始未施加流动控制构型,阻力系数减小量可达到32.72%。     

横摆振荡翼型动态气动掠效应试验研究

大型风力机设计对获取翼型更加全面、准确的动态载荷提出更高要求, 研究翼型横摆振荡动态气动特性具有重要意义. 借助"电子凸轮"技术和动态数据同步采集手段, 针对翼型动态“掠效应”首次开展了横摆振荡风洞试验研究, 研究表明: 横摆振荡翼型的气动曲线存在明显迟滞效应, 吸力面压力周期性波动是主要诱因, 且随着振荡频率、初始迎角和振幅的增大, 气动迟滞特性均增强; 升力和压差阻力随横摆角变化的迟滞回线呈"W"形, 俯仰力矩迟滞回线呈"M"形, 升力差量迟滞回线呈"$\infty$"形; 负行程下翼型气动力相对于正行程下的更高, 且负行程下翼型气动力随振荡频率的增大而略有增大, 正行程下则明显减小; 升力系数功率谱密度分布在振荡频率倍频处的能量集中的幅值随着振荡频率增大有增大趋势; 吸力面1.2%和40%弦长处压力的滞回特性较强, 是由于翼面剪切层涡和动态分离涡周期性发展、运动、破裂和重建; 振幅为$10^{\circ}$时, 升力迟滞曲线呈"$^{\wedge}$"形, 振幅为$30^{\circ}$ 时, 升力迟滞曲线呈"$^{\wedge\wedge\wedge}$"形.      

Effect of film cooling on the aerodynamic performance of an airfoil

The effect of film cooling on the aerodynamic performance of turbine blades is becoming increasingly important as the gas turbine operating temperature is being increased in order to increase the performance. The current paper investigates the effect of blowing ratio on the aerodynamic losses of a symmetric airfoil by pressure measurements and Particle Image Velocimetry (PIV). The test model features 4 rows of holes located on the suction side at 5%, 10%, 15% and 50% of the chord length. The Reynolds number based on the airfoil chord is 1.2 × 10⁵. Experiments are performed by varying the location of air injection, the angle of attack, and the mainstream velocity. The coolant air is injected at ambient temperature and the blowing ratio is varied from 0 to 1.91. It is observed that the losses due to film cooling increase with blowing ratio of 0 to 0.48, and the wake is shifted towards the suction side. Conversely, the aerodynamic losses decrease when the blowing ratio is increased further from 0.64 to 1.91. This trend has been observed for all the experimental configurations. The effect of blowing ratio on flow separation is investigated with the time-averaged velocity fields obtained from PIV measurements. It is observed that low blowing ratios, the separation point shifts upstream and at high blowing ratios the ejected coolant energizes the flow and delays separation. The pressure field around the airfoil is reconstructed from the integration of the Poisson equation based on the PIV velocity fields. The experimental results can be used for validation of numerical models for predicting losses due to film cooling.