Mathematical simulation of unsteady flow past a wing with jets

The various approximate approaches to the investigation of the unsteady aerodynamic characteristics of an airfoil with jet flap [1–3] are applicable only for an airfoil, low jet intensity, and low oscillation frequencies. In the present paper, the method of discrete vortices [4] is generalized to the case of unsteady flow past a wing with jets and arbitrary shape in plan. The problem is solved in the linear formulation; the conditions used are standard: no flow through the wing and jet, finite velocities at the trailing edges where there is no jet, and also a dynamical condition on the jet. The wing and jet are assumed to be thin and the medium inviscid and incompressible.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 139–144, May–June, 1982.     

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

对于翼面变形速度远小于来流速度情况下的儒可夫斯翼型亚音速绕流问题,通过仿射变换将可压缩流动转换成不可压缩流动,将解析解和离散涡方法相结合计算变形机翼的不可压缩流动速度场,再利用逆变换得到变形机翼的亚音速流动速度场,进而分析非定常气动力特性,建立变形机翼的准定常升力系数和非定常附加升力系数在可压缩和不可压缩两种状态下的简单近似对应关系。计算结果显示变形机翼的非定常气动升力近似等于准定常计算结果叠加上虚拟质量力导致的非定常附加升力,该非定常附加升力随翼型变形速率呈线性关系,由机翼当前时刻飞行姿态、翼型及其变形速率确定,与具体变形历史过程无关。低来流马赫数时虚拟质量力导致的非定常效应显著,高亚音速流动时准定常升力起主导作用。同时还分析了不同马赫数下机翼往复变形过程中升力的变化特性,指出尽管高亚音速变形机翼的气动升力近似等于准定常气动升力,但不能忽视非定常附加升力的影响,非定常附加升力将导致完成往复变形需要外界输入正比于Ma∞/[(1-Ma2∞)]的功。     

Rotary aerodynamic derivatives in the asymptotic theory of a high-aspect-ratio wing

The linear problem of inviscid incompressible flow around a high-aspect-ratio wing at an angle of attack and in the presence of steady pitching and rolling rotation is considered. The main integral equation of the problem is reduced to a sequence of one-dimensional integral equations without use of the matched asymptotic expansions method. The first few terms of the series for the circulation distribution over the wing surface are calculated. For an elliptic high-aspect-ratio wing the corresponding aerodynamic forces are calculated. The derivatives of the aerodynamic coefficients of the wing with respect to the angle of attack and the angular velocities are determined. The asymptotic expressions obtained are compared with the results of numerical calculations of the corresponding derivatives using the discrete vortex method.     

模型昆虫翼作非定常i运动时的气动力特性

基于Navier-Stokes方程的数值解,研究了一模型昆虫翼在小雷诺数(Re=100)下作非定常运动时的气动力特性.这些运动包括翼启动后的常速转动,快速加、减速转动,常速转动中快速上仰(模拟昆虫翼的上挥或下拍、翻转等运动).有如下结果在小雷诺数下,模型昆虫翼以大攻角(α=35°)作常速转动运动时,由于失速涡不脱落,可产生较大的升力系数.其机理是翼转动时,翼尖附近(该处线速度大)上翼面压强比翼根附近(该处线速度小)的小得多,因而存在展向压强梯度,同时存在着沿展向的离心力,此展向压强梯度和离心力导致的展向流动在失速涡的轴向方向,其可避免失速涡脱落.模型昆虫翼在快速加、减速转动和快速上仰运动中,虽然雷诺数小,但由于在短时间内产生了大涡量,也可产生十分大的气动力,例如在快速上仰运动中,升力系数可大于10.     

滑动蒙皮变后掠气动力非定常滞回与线性建模

针对低速不可压条件下滑动蒙皮方式变后掠过程中非定常动态气动特性开展了3方面的研究工作: (1)飞行器变形过程中非定常动态气动特性风洞试验技术研究;(2)变形过程中滞回效应研究和机理分析; (3)基于风洞试验结果开展变形过程中非定常动态气动力线性建模. 初步研究表明: (1)采用强迫振荡法可以有效地获取变形过程中非定常动态气动力滞回效应; (2)造成变形过程中气动滞回效应的机理有两个即``动边界效应'和``流场结构滞回效应', 造成滞回效应的机理可能主要在于后者; (3)引入升力系数和俯仰力矩系数随后掠角变化率的动导数概念, 可以建立变形过程中非定常动态气动力线性模型.      

微型飞行器低雷诺数空气动力学

微型飞行器(MAVs)设计绝不是常规飞行器在尺度上的简单缩小,面临许多技术难题.其中微型飞行器低雷诺数空气动力学是其最为根本的技术瓶颈之一,也是当前受到广泛关注的热点之一.本文紧密结合微型飞行器技术,对这一领域中所面临的低雷诺数空气动力学问题和近两年来该方向国内一些新的进展进行了较为详细的介绍.按照MAVs飞行方式和结构特性进行分类,简单介绍微型飞行器研究中的低$Re$数空气动力学问题.首先介绍了二维和三维固定翼低雷诺数空气动力学问题:包括层流分离泡,翼型升力系数小攻角非线性效应,静态迟滞效应,以及低$Re$数小展弦比机翼气动特性.第2,介绍了拍动翼低雷诺数空气动力学方面的研究工作.包括前人提出的昆虫低$Re$数下获得高升力的多种非定常拍动翼飞行机制:Wagner效应、Weis-Fogh效应(clap-and-fling)、延迟失速效应(delayedstall)、Kramer效应(rotational forces)、尾迹捕获效应(wakecapture)、附加质量效应(addedmass)等.以及国内学者近几年在拍动翼方面取得的一些研究成果.第3,介绍了柔性翼低雷诺数气动问题.研究表明柔性翼对于固定翼微型飞行器提高抗阵风能力,拍动翼微型飞行器产生足够的升力和推力.最后简单介绍了可变形翼(morphingwing)微型飞行器方面的一些研究工作,指出微型飞行器技术可以通过采用可变形翼设计,突破众多的技术瓶颈.另一方面,可变形翼概念可以通过在低成本,低速的MAVs上进行飞行试验,获得非常好的验证平台.      

应用当地流活塞理论的大攻角升力面颤振气动力表达式

本文应用当地流活塞理论,给出了弹性的夺动翼面的非定常压分分布以及用模态坐标的广义气动力系数表达式。同时提供了配套的当地流参数计算公式和采用样条函数计算广义气动力系数数值积分的表达式。     

Aeroelastic instability of a composite wing with a powered-engine

In this paper the aeroelastic instability of a wing, modeled as an orthotropic composite beam with a concentrated mass subjected to the engine thrust, is investigated in an incompressible flow. The wing is modeled using classical beam theory. Wagner function is used to model the unsteady aerodynamic loads, while the engine thrust is modeled as a follower force and a concentrated mass is used to model the engine mass. The numerical results of the developed generic and simple model are compared with published results, and an excellent agreement is observed. The fiber orientation, engine thrust, mass magnitude and its location are reported to have had significant effects on the aeroelastic instability boundaries.     

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

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

两种湍流模型时域颤振计算方法研究

采用时域计算分析方法进行了机翼跨音速颤振特性研究。在结构运动网格的基础上,采用格点格式有限体积方法进行空间离散和双时间全隐式方法进行时间推进求解雷诺平均N-S方程。针对流动粘性分别应用了SST湍流模型和SSG雷诺应力模型,通过对跨音速标模算例AGARD445.6机翼的计算结果与实验值的对比分析,其中应用SST湍流模型得到的颤振速度与实验值最为接近,特别是在跨音速段平均相对误差在3%以内;并且计算结果整体上反映了跨音速颤振"凹坑"物理特性,验证了方法的有效性。     

旋转流场下模型大幅俯仰运动实验研究

在中航气动院FL-8低速风洞中,采用单自由度振荡机构进行了旋转流场下大幅俯仰运动的气动特性实验研究。模型在绕风速轴连续旋转的同时,进行给定频率和振幅绕体轴的俯仰振荡运动。测量了模型的动态气动特性,着重分析了不同运动参数对模型气动特性的影响。实验结果表明,旋转速度的存在使大振幅俯仰振荡实验中的滚转力矩和偏航力矩产生了明显的迟滞特性,但对俯仰力矩和法向力的迟滞特性影响不大。     

Efficient unsteady aerodynamic loads prediction based on nonlinear system identification and proper orthogonal decomposition

In the present work, an efficient surrogate-based framework is developed for the prediction of motion-induced surface pressure fluctuations and integral force and moment coefficients. The model construction is realized by performing forced-motion computational fluid dynamics (CFD) simulations, while the result is processed via the proper orthogonal decomposition (POD) to obtain the predominant flow modes. Subsequently, a nonlinear system identification is carried out with respect to the applied excitation and the resulting POD coefficients. For the input/output model identification task, a recurrent local linear neuro-fuzzy approach is employed in order to capture the linear and nonlinear characteristics of the dynamic system. Once the reduced-order model (ROM) is trained, it can substitute the flow solver within unsteady aerodynamic or aeroelastic simulation frameworks for a given configuration at fixed freestream conditions. For demonstration purposes, the ROM approach is applied to the LANN wing in high subsonic and transonic flow. Due to the characteristic lambda-shock system, the unsteady aerodynamic surface pressure distribution is dominated by nonlinear effects. Numerical investigations show a good correlation between the results obtained by the ROM methodology in comparison to the full-order CFD solution. In addition, the surrogate approach yields a significant speed-up regarding unsteady aerodynamic calculations, which is beneficial for multidisciplinary computations.     

Nonlinear characteristics of a thin wing in a shock-free flow around the nose

An effective means of controlling wing leading-edge stall at high angles of attack is deflection of the nose in order to assure shock-free entrance of the stream. A numerical method of computing the angles of nose deflection and the aerodynamic characteristics of a thin wing of arbitrary planform for a shock-free entrance of the steady ideal incompressible fluid stream is elucidated in this paper on the basis of nonlinear wing theory [1]. The problem is solved by the method of discrete vortices. In the computations, the wing and its wake, replaced by a vortex sheet, are modeled by a system of discrete vortices which are nonlinear segments with constant circulation along the length. The angles of deflection of the nose and the aerodynamic characteristics of the wing, including shunting of the free vortices shed from the side and trailing edges, are determined during the computation. Examples of an electronic digital computer are presented.     

Computational study of a pitching bio-inspired corrugated airfoil

The phenomenon of insect flight has been of scientific interest for many years and is more recently inspiring modern engineering devices such as Micro Aerial Vehicles (MAVs). Insect flight is characterized by unsteady fluid dynamics at low Reynolds numbers. The importance of viscous effects to the successful flapping flight of insects has been identified and with the current state of computing and Computational Fluid Dynamics (CFD) these effects can now be studied in detail. The present work attempts to simplify this complex phenomenon by considering symmetric oscillating rotational motion of a wing (pitching). What is of interest in this study is how the shape of a corrugated idealized insect wing affects the performance and flow characteristics around the pitching wing. Two dimensional CFD on an oscillating wing has been performed and reported. Measurements were taken to ensure the accuracy of the computational solution and the results validated against experimental PIV results. A range of frequencies and rotational amplitudes have been investigated. Lift and drag coefficients have been analyzed for all cases to quantify the effects of unsteady flow features on the performance of the oscillating wing. It was found that the wing shape used in this study resulted in the viscous features formed on the top of the wing exhibiting high sensitivity to the oscillating conditions and these influenced the performance of the wing. The flow features formed on the bottom of the wing remained similar throughout the cases tested. In the pitching regime this wing profile did not perform as well as published results for smooth airfoils in terms of thrust and propulsion efficiency. However this may be due to reduced frequency effects becoming important at our high pitching amplitude which need to be investigated further. There may be other oscillatory regimes that more accurately represent flapping flight in which the corrugated foil outperforms a smooth counterpart but these are yet to be investigated. Further research in this area may help answer the question as to how evolutionarily significant other benefits of a corrugated wing, such as being light and strong, are compared to its aerodynamic properties, the present results seem to favor the former.     

Heaving motions of a ram wing translating above water

Some high-speed marine vehicles utilize wings operating in the air but in close proximity to the water. At small clearances under the wings, and with augmentation of the incident airflow by front air-based propulsors, water surface deformations can influence the aerodynamic characteristics of such lifting surfaces. A dynamic model is developed in this paper for unsteady forces on a two-dimensional ram wing. The current model is based on the extreme-ground-effect theory for airflow and a linearized potential flow theory for water surface deformations. Parametric results are presented for aerodynamic characteristics of a flat plate, including a static lift coefficient, sensitivity on height, and added mass and damping coefficients in heaving motions. Variable parameters include relative speed, airflow augmentation, and frequency and amplitude of heaving oscillations. The developed method can be applied for more general dynamic modeling of high-speed air-supported marine craft.     

Parametric studies on relationships between flutter derivatives of slender bridge （Ⅱ）

Based on curve fitting of coefficients of three component forces of the Messina Straits Bridge, and the previously proposed semi-analytical expressions of flutter derivatives of flexible structure, the change of flutter derivatives of slender bridge cross-section with respect to its aerodynamic center, rotational speed and angle variation is studied using a parametric method. The calculated results are compared with the measured ones, and expressions of flutter derivatives of the Messina Straits Bridge are derived. The intrinsic relationships existing in flutter derivatives are validated again. It is shown that the influence of the rotational speed on flutter derivatives is not negligible. Therefore, it provides an additional semi-analytical approach for analyzing flutter derivatives of the bridge with streamlined cross-section to get its aerodynamic information.     

Analytical formulation of 2-D aeroelastic model in weak ground effect

This paper deals with the aeroelastic modeling and analysis of a 2-D oscillating airfoil in ground effect, elastically constrained by linear and torsional springs and immersed in an incompressible potential flow (typical section) at a finite distance from the ground. This work aims to extend Theodorsen theory, valid in an unbounded flow domain, to the case of weak ground effect, i.e., for clearances above half the airfoil chord. The key point is the determination of the aerodynamic loads, first in the frequency domain and then in the time domain, accounting for their dependence on the ground distance. The method of images is exploited in order to comply with the impermeability condition on the ground. The new integral equation in the unknown vortex distribution along the chord and the wake is solved using asymptotic expansions in the perturbation parameter defined as the inverse of the non-dimensional ground clearance of the airfoil. The mathematical model describing the aeroelastic system is transformed from the frequency domain into the time domain and then in a pure differential form using a finite-state aerodynamic approximation (augmented states). The typical section, which the developed theory is applied to, is obtained as a reduced model of a wing box finite element representation, thus allowing comparison with the corresponding aeroelastic analysis carried out by a commercial solver based on a 3-D lifting surface aerodynamic model. Stability (flutter margins) and response of the airfoil both in frequency and time domains are then investigated. In particular, within the developed theory, the solution of the Wagner problem can be directly achieved confirming an asymptotic trend of the aerodynamic coefficients toward the steady-state conditions different from that relative to the unbounded domain case. The dependence of flutter speed and the frequency response functions on ground clearance is highlighted, showing the usefulness of this approach in efficiently and robustly accounting for the presence of the ground when unsteady analysis of elastic lifting surfaces in weak ground effect is required.     

Parametric studies on relationships between flutter derivatives of slender bridge (Ⅱ)

Based on curve fitting of coefficients of three component forces of the Messina Straits Bridge, and the previously proposed semi-analytical expressions of flutter deriva-tives of flexible structure, the change of flutter derivatives of slender bridge cross-section with respect to its aerodynamic center, rotational speed and angle variation is studied using a parametric method. The calculated results are compared with the measured ones, and expressions of flutter derivatives of the Messina Straits Bridge are derived. The in-trinsic relationships existing in flutter derivatives are validated again. It is shown that the influence of the rotational speed on flutter derivatives is not negligible. Therefore, it provides an additional semi-analytical approach for analyzing flutter derivatives of the bridge with streamlined cross-section to get its aerodynamic information.     

Calculation of the aerodynamic characteristics of a helicopter lifting rotor subjected to air gusts

A numerical method of calculating the unsteady aerodynamic characteristics of a helicopter lifting rotor smoothly enveloped by an air gust arriving from an arbitrary direction is proposed. The problem is solved in the non-linear formulation for an incompressible medium. The flow over each blade is simulated by a system of discrete vortices. The deformation of the vortex trail from the blades under the influence of the gust is taken into account. The results of the calculations are reported.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 153–158, March–April, 1987.     

Aerodynamic characteristics of a deltawing in the case of harmonic oscillations with respect to the roll and yaw angles at high angles of attack

The aerodynamic characteristics of a delta wing in the case of harmonic oscillations with respect to the roll and yaw angles are obtained in a subsonic low-speed wind tunnel and analyzed. It is shown that at near-critical angles of attack the aerodynamic derivatives of the roll moment considerably depend on the reduced oscillation frequency. It is established that this dependence is due to a variation in the slip angle. A mathematical model that involves an ordinary linear differential first-order equation is used to describe the aerodynamic characteristics of the wing for the problems of aircraft flight dynamics at high angles of attack.