可变形儒可夫斯基翼型非定常气动力的研究

对于翼面变形法向运动速度远小于来流速度的儒可夫斯基机翼,将解析解和离散涡方法相结合计算变形机翼的流场及非定常气动力,较详细地分析了变形机翼升力系数的准定常计算方法的误差来源,并给出修正方法.计算结果显示脱落涡尾迹对升力系数和机翼绕流环量的影响很小,变形机翼升力系数准定常计算方法的误差丰要来源于流体非定常运动引起的虚拟质量力,该非定常附加升力仅与当前时刻飞行姿态及翼犁形状和变形速率有关,与具体的变形历史过程无关,变形机翼的升力近似等于准定常计算结果叠加上相应的虚拟质量力.     

低Re数下变弯度机翼的非定常气动特性实验研究

利用变弯度机翼模型及相关的风洞实验平台,开展了以弯度变化速率影响为重点的机翼非定常特性研究。实验结果显示,在低Re数(~105)下,机翼弯度非定常变化得到的升阻力系数曲线与准定常条件下的结果存在显著差异。具体表现为:准定常状态下,曲线表现出明显的可逆性;而弯度非定常变化时,曲线在弯度递增区和递减区之间存在明显的迟滞效应,而且随着变形速率的增加,这种迟滞也越明显。流场显示结果表明,这种小St数下出现的流动迟滞是由于弯度变形导致的流动分离的分离点相对机翼运动迟滞所造成的。这说明弯度变化时,分离流场结构的响应时间尺度与弯度变化周期相当,也揭示了该条件下机翼弯度变化对流动的抑制作用主要是通过改变分离区的大小来实现的。     

亚音速气流中复合材料悬臂板的非线性振动响应研究

随着材料科学的发展,越来越多的新型材料应用到了工程实践中.在气流激励的作用下,对于以航空航天工程为背景、采用复合材料的板壳结构的非线性动力学问题仍是动力学领域的研究热点.本文研究了复合材料悬臂板在亚音速气流条件下的非线性振动和响应.根据理想不可压缩流体的流动条件和Kutta–Joukowski升力定理,基于升力面理论,利用涡格法计算了三维有限长平板机翼上的亚音速气动升力.将亚音速气动力施加到复合材料悬臂板上,利用Hamilton原理,考虑Reddy三阶剪切变形理论并引入冯·卡门非线性应变位移关系,建立了有限长平板的非线性动力学微分方程.利用有限元方法考察了不同几何参数下层合板悬臂板的固有特性,通过比较不同材料和几何参数的线性系统的固有频率,得到不同比例的内共振关系.利用Galerkin方法将偏微分方程截断为两自由度非线性常微分方程,在这里考虑了1:2的内部共振关系并利用多尺度法进行了摄动分析.对应多个选取参数,得到了频率响应曲线.结果展示了硬化弹簧型行为和跳跃现象.     

亚音速气流中复合材料悬臂板的非线性振动响应研究

随着材料科学的发展,越来越多的新型材料应用到了工程实践中.在气流激励的作用下,对于以航空航天工程为背景、采用复合材料的板壳结构的非线性动力学问题仍是动力学领域的研究热点.本文研究了复合材料悬臂板在亚音速气流条件下的非线性振动和响应.根据理想不可压缩流体的流动条件和 Kutta--Joukowski升力定理,基于升力面理论,利用涡格法计算了三维有限长平板机翼上的亚音速气动升力.将亚音速气动力施加到复合材料悬臂板上,利用Hamilton原理,考虑Reddy三阶剪切变形理论并引入冯$\cdot$卡门非线性应变位移关系,建立了有限长平板的非线性动力学微分方程.利用有限元方法考察了不同几何参数下层合板悬臂板的固有特性,通过比较不同材料和几何参数的线性系统的固有频率,得到不同比例的内共振关系.利用Galerkin方法将偏微分方程截断为两自由度非线性常微分方程,在这里考虑了1:2的内部共振关系并利用多尺度法进行了摄动分析.对应多个选取参数,得到了频率响应曲线.结果展示了硬化弹簧型行为和跳跃现象.      

UNSTEADY SEPARATION OF FLOW AROUND AIRFOIL WITH LOCAL ELASTIC STRUCTURE

对在低雷诺数下局部弹性翼型绕流中,局部弹性导致的自激振动所产生的复杂非定常流动分离现象和描述方法进行了分析．采用ALE—CBS方法数值模拟了具有可动边界的绕流流场问题,同时采用Galerkin方法求解局部弹性结构的控制方程．着重研究了翼型的局部弹性对流动分离和翼型性能的影响,并分别从Eulerian和Lagrangian的角度分析了局部弹性结构导致的不同非定常分离现象,其中Lagrangian角度可以方便地揭示出局部弹性翼型大幅度提高升力的机理和流动中的能量迁移．结果表明翼型的局部弹性对非定常分离和分离泡的演化过程有着明显的影响,可以使得流体质点由主流获取动量实现再附,并且在一定的攻角下可以将固定分离转变为移动分离,从而明显地提高了翼型的升力．     

局部弹性翼型非定常分离的动力学特性

对在低雷诺数下局部弹性翼型绕流中, 局部弹性导致的自激振动所产生的复杂非定常流动分离现象和描述方法进行了分析. 采用ALE-CBS方法数值模拟了具有可动边界的绕流流场问题, 同时采用Galerkin方法求解局部弹性结构的控制方程. 着重研究了翼型的局部弹性对流动分离和翼型性能的影响, 并分别从Eulerian和Lagrangian的角度分析了局部弹性结构导致的不同非定常分离现象, 其中Lagrangian角度可以方便地揭示出局部弹性翼型大幅度提高升力的机理和流动中的能量迁移. 结果表明翼型的局部弹性对非定常分离和分离泡的演化过程有着明显的影响, 可以使得流体质点由主流获取动量实现再附, 并且在一定的攻角下可以将固定分离转变为移动分离, 从而明显地提高了翼型的升力.      

一种可用于机翼优化设计中的气动噪声预测模型

发展了一种可用于翼型/机翼外形设计中的气动噪声快速预测方法。相较于传统的半经验噪声预测方法,该方法以两方程非线性k-ε湍流模型模化雷诺应力的雷诺平均方程为背景,考虑了升力系数、三维流动效应以及机翼几何参数等因素对后缘噪声的影响。而相对于直接数值模拟或声类比拟方法,该方法虽不能准确预测噪声强度,但其计算量小,能给出不同翼型/机翼的相对总声压级,以及总声压级随升力系数的变化情况,易于应用于翼型/机翼气动外形优化设计中。通过计算分析二维NACA0012翼型几何参数或来流状态的改变所带来的气动噪声差异,与ANOPP软件及Brooks等计算结果进行对比,验证了该模型的可靠性。最后,计算分析二维、三维翼型/机翼气动噪声,凸显该方法在翼型/机翼气动外形优化设计中的应用价值。     

钝头翼型和机翼跨音速绕流的渐近展开分析

本文利用渐近展开匹配法分析钝头翼型的跨音速绕流,导出了描述前缘附近流动的一级近似、二级近似下的速位方程、边界条件及相应的近似解析解,并构成关于翼表面速度的一致有效合成解,消除了跨音速小扰动近似的前缘奇性,对于大展弦比后掠翼绕流,可利用翼型绕流分析结果,消除机翼前缘奇性。     

七孔探针可压缩流场测量研究

介绍了七孔探针用于亚音速可压缩流的标定方法。作为一种可以同时获得流动速度大小、流动偏角、总压和静压的气动测量装置,七孔探针被广泛应用于各种流动测量,包括可压缩流动。但是它的校准过程周期很长,代价昂贵,影响了探针的推广。本文以数值计算为手段,对七孔探针进行亚音速可压缩流校准与测量的研究。结果表明,其校准拟合精度流动角为2%,内外区的总静压相对标准偏差都没有超过3%,高于相同状态下的实验校准精度。在实际应用中,本方法用于指导传统实验标定方法,可以节约大量的标定时间和成本,使七孔探针在亚音速可压缩流的测量变得简单可行。     

Level Set方法在地面效应数值模拟中的应用

利用Level Set方法隐式捕捉界面并能处理复杂的物质界面及其拓扑结构变化（如合并、交叉、破碎等）情况的优势,基于可压缩N-S方程的预处理技术,采用基于界面两侧介质黎曼关系的修正虚拟流体方法对介质边界进行处理;将只能应用于可压缩流动界面模拟的修正虚拟流体方法推广应用于不可压缩流动问题的模拟,数值模拟了气-水运动界面与近水面飞行器的相互影响规律。通过与将水面视为对称面的简化模型对比发现：利用自由运动界面模型计算的NACA6409翼型升力较小,能更准确地反映飞行器流场对水面形状的影响,其结果更为符合实际;文中还对比了翼型距水面不同高度时的结果,得到了翼型受水面影响的强弱与翼型距水面的高度成反比的结论。     

翼型连续变形过程中非定常气动特性研究

可变形飞行器的优越性,使其成为近年来研究者关注的热点。当前对于可变形飞行器气动方面的研究工作主要集中于采用数值模拟或风洞实验,开展可变形飞行器气动布局概念和不同变形模态气动特性研究,而关于连续变形过程中非定常气动特性的研究则很少,几为空白。本文在无粘不可压(Ma=0.2)流场中,数值模拟了儒可夫斯基翼型连续改变厚度、弯度所引起的气动特性变化,分析了不同变形幅度、变形速度以及变形加速度对非定常升力系数的影响情况,这也为以后的三维数值模拟提供了基础。     

Experimental study of pitching and plunging airfoils at low Reynolds numbers

Measurements of the unsteady flow structure and force time history of pitching and plunging SD7003 and flat plate airfoils at low Reynolds numbers are presented. The airfoils were pitched and plunged in the effective angle of attack range of 2.4°–13.6° (shallow-stall kinematics) and ?6° to 22° (deep-stall kinematics). The shallow-stall kinematics results for the SD7003 airfoil show attached flow and laminar-to-turbulent transition at low effective angle of attack during the down stroke motion, while the flat plate model exhibits leading edge separation. Strong Re-number effects were found for the SD7003 airfoil which produced approximately 25 % increase in the peak lift coefficient at Re = 10,000 compared to higher Re flows. The flat plate airfoil showed reduced Re effects due to leading edge separation at the sharper leading edge, and the measured peak lift coefficient was higher than that predicted by unsteady potential flow theory. The deep-stall kinematics resulted in leading edge separation that led to formation of a large leading edge vortex (LEV) and a small trailing edge vortex (TEV) for both airfoils. The measured peak lift coefficient was significantly higher (~50 %) than that for the shallow-stall kinematics. The effect of airfoil shape on lift force was greater than the Re effect. Turbulence statistics were measured as a function of phase using ensemble averages. The results show anisotropic turbulence for the LEV and isotropic turbulence for the TEV. Comparison of unsteady potential flow theory with the experimental data showed better agreement by using the quasi-steady approximation, or setting C(k) = 1 in Theodorsen theory, for leading edge–separated flows.     

A study on the mechanism of high-lift generation by an airfoil in unsteady motion at low reynolds number

The aerodynamic force and flow structure of NACA 0012 airfoil performing an unsteady motion at low Reynolds number (Re=100) are calculated by solving Navier-Stokes equations. The motion consists of three parts: the first translation, rotation and the second translation in the direction opposite to the first. The rotation and the second translation in this motion are expected to represent the rotation and translation of the wing-section of a hovering insect. The flow structure is used in combination with the theory of vorticity dynamics to explain the generation of unsteady aerodynamic force in the motion. During the rotation, due to the creation of strong vortices in short time, large aerodynamic force is produced and the force is almost normal to the airfoil chord. During the second translation, large lift coefficient can be maintained for certain time period and , the lift coefficient averaged over four chord lengths of travel, is larger than 2 (the corresponding steady-state lift coefficient is only 0.9). The large lift coefficient is due to two effects. The first is the delayed shedding of the stall vortex. The second is that the vortices created during the airfoil rotation and in the near wake left by previous translation form a short “vortex street” in front of the airfoil and the “vortex street” induces a “wind”; against this “wind” the airfoil translates, increasing its relative speed. The above results provide insights to the understanding of the mechanism of high-lift generation by a hovering insect. The project supported by the National Natural Science Foundation of China (19725210)     

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

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

风力机叶片的三维非定常气动特性估算

结合动量-叶素理论、非定常空气动力和动态失速模型来计算风力机叶片的二维非定常气动特性,并在此基础上经过适当修正后考虑三维旋转效应的非定常气动特性。分析比较二维和三维两种计算结果,给出更为合理的计算叶片非定常气动特性的方法。计算结果表明,风力机叶片的三维非定常气动特性计算结果与二维时的计算结果相比有较大改善。     

Low Reynolds unsteady flow simulation around NACA0012 airfoil with active flow control

This research numerically elucidates the effects of suction and blowing on the enhancement of unsteady aerodynamic characteristics of flows and their corresponding impact on stall delay over the well-known NACA0012 airfoil at various angles of attack (\( 12 \le {\text{AOA}} \le 20 \)) under low Reynolds numbers. For this purpose, an in-house solver written in C++ is developed. The numerical code utilizes the Jameson’s cell-centered finite volume numerical method accompanied by a progressive power-law preconditioning approach to suppress the stiffness of the governing equations. Many numerical simulations are performed over the suction-blowing control parameters, namely, the slot location (\( L_{j} \)), suction/blowing amplitudes (\( A_{j} \)), and suction/blowing angle (\( \theta_{j} \)). Most of the analyses are based on the measurements of the unsteady aerodynamic characteristics behaviors (such as lift, drag, moment coefficients, and stall phenomena) over the airfoil. The numerical results confirm that the unsteady behavior of the flow (vortex shedding) is weakened or approximately removed when suction is used, especially near the leading edge. In all of the test cases, the ratio of the average lift coefficient to the average drag coefficient increases with increasing suction and blowing amplitudes, except in the case of perpendicular blowing. Furthermore, the blowing is more sensitive to the blowing angle compared to the suction. From the suction and blowing results, it is concluded that the former has a more positive impact on the lift and drag characteristics, especially in the case of incompressible flow at Low-Reynolds regimes.     

顺桨风力机塔叶干涉翼型非稳气动力时频特性

由于风力机叶片与塔筒流场相互干涉,实际气动力与理想情况存在较大差异,这种干涉作用造成的气动力差异给叶片与塔筒结构可靠性带来不可忽视的影响.以翼型DU91-W2-250为研究对象,采用瞬态数值分析与本征正交分解方法,考虑叶片和塔筒流场相互干涉作用,分析顺桨工况翼型非稳气动力时频特性及其影响规律,量化不同雷诺数下塔叶相对位置及几何参数对气动力均值、波动幅度和频率的影响程度,通过流场模态能量分布形态分析,揭示流场干涉对气动力的影响机制.结果表明,翼型气动中心至塔筒几何中心的垂直距离、水平距离以及塔筒直径相对于翼型弦长的无量纲参数y*,x*和D*对气动力均有不同程度影响,其中y*对升阻力系数均值影响最大,对频率无明显影响,y*绝对值越大,C_l均值越接近单翼型C_l值,y*绝对值越小升阻力系数波动幅度越大,y*从-12增大到12,升力系数均值最小值为-0.48,最大值为1.16;x*减小和D*增大,反向阻力均值增大,波动幅度增大,波动频率略有下降,当x*小于临界值5时,带塔翼型阻力均值反向;在计算范围内,带塔翼型升力系数均值相对于单翼型升力系数最大偏差为...     

Aerodynamics and mechanisms of elementary morphing models for flapping wing in forward flight of bat

Large active wing deformation is a significant way to generate high aerodynamic forces required in bat's flapping flight. Besides the twisting, elementary morphing models of a bat wing are proposed, including wing-bending in the spanwise direction, wing-cambering in the chordwise direction, and wing area-changing. A plate of aspect ratio 3 is used to model a bat wing, and a three-dimensional unsteady panel method is used to predict the aerodynamic forces. It is found that the cambering model has great positive influence on the lift, followed by the area-changing model and then the bending model. Further study indicates that the vortex control is a main mechanism to produce high aerodynamic forces. The mechanisms of aerodynamic force enhancement are asymmetry of the cambered wing and amplification effects of wing area-changing and wing bending. Lift and thrust are generated mainly during downstroke, and they are almost negligible during upstroke by the integrated morphing model-wing.     

Effect of vortex shedding in unsteady aerodynamic forces for a low Reynolds number stationary wing at low angle of attack

This study elucidates the relation between wake vortex shedding and aerodynamic force fluctuations for a low Reynolds number wing from time resolved particle image velocimetry (TR-PIV) experimental measurements. The results reveal a periodic lift and drag variation within the shedding cycle and resolve the frequencies of those fluctuations from a proper orthogonal decomposition (POD) and power spectral density (PSD) analysis. To show the effect of vortex shedding on the body force fluctuations, the evolution of instantaneous aerodynamic forces is compared to the pressure field of the fluid flow and to the vortical structures in the wake of the airfoil. A six step model describing the vortex-force relation is proposed. It shows that changes in lift such as maximum lift and minimum lift are associated with the detachment of a vortex. It also shows that the minimum or local minimum drag value is obtained at the onset formation of a vortex on the airfoil wake. Similarly, the maximum or local maximum drag is obtained at the onset formation of the saddle on the airfoil wake. The model further explains the asymmetry observed in the unsteady drag force evolution. The model can be used to optimize flow control and fluid-structure interaction applications.     

低雷诺数翼型蒙皮主动振动气动特性及流场结构数值研究

针对低雷诺数(Re)翼型气动性能差的特点,文章通过对翼型柔性蒙皮施加主动振动的方法,提高翼型低Re下的气动特性,改善其流场结构.采用带预处理技术的Roe方法求解非定常可压缩Navier-Stokes方程,对NACA4415翼型低Re流动展开数值模拟.通过时均化和非定常方法对比柔性蒙皮固定和振动两种状态下的升阻力气动特性和层流分离流动结构.初步研究工作表明在低Re下柔性蒙皮采用合适的振幅和频率,时均化升阻力特性显著提高,分离泡结构由后缘层流分离泡转变为近似的经典长层流分离泡,分离点后移,分离区缩小.在此基础上,文章更加细致研究了柔性蒙皮两种状态下单周期内的层流分离结构及壁面压力系数分布非定常特性和演化规律.蒙皮固定状态下分离区前部流场结构和压力分布基本保持稳定,表现为近似定常分离,仅在后缘位置出现类似于卡门涡街的非定常流动现象.柔性蒙皮振动时从分离点附近开始便产生分离涡,并不断向下游移动、脱落,表现为非定常分离并出现大范围的压力脉动.蒙皮振动使流体更加靠近壁面运动,大尺度的层流分离现象得到有效抑制.