Computational aeroelastic investigation of a transonic limit-cycle-oscillation experiment at a transport aircraft wing model

Aeroelastic measurements of a three-dimensional wing model, the so-called Aerostabil wing, were conducted in the Transonic Windtunnel Göttingen. This clean, backward-swept wing allowed the experimental investigation of limit cycle oscillations in a certain transonic parameter range. In this paper, a detailed insight into the observed physical phenomena, especially the measured limit cycle oscillations, is presented by means of CFD–CSM coupled simulations. These simulations on the basis of a detailed structural finite element model reveal the specific properties of the Aerostabil wing and furthermore allow investigating the unstable behavior of this windtunnel model for transonic flow settings. The aerodynamic characteristics include a two-shock system and large flow separation areas, further increasing the complexity of the aeroelastic problem. A structural single degree-of-freedom system is used for the prediction of the experimental stability range and the limit cycle oscillation investigations. Due to the good agreement of simulation and experiment the limit cycle oscillations can be explained by means of nonlinear aerodynamic effects.     

Limit cycle oscillation behavior of transonic control surface buzz considering free-play nonlinearity

The limit cycle oscillation (LCO) behaviors of control surface buzz in transonic flow are studied. Euler equations are employed to obtain the unsteady aerodynamic forces for Type B and Type C buzz analyses, and an all-movable control surface model, a wing/control surface model and a three-dimensional wing with a full-span control surface are adopted in the study. Aerodynamic and structural describing functions are used to deal with aerodynamic and structural nonlinearities, respectively. Then the buzz speed and buzz frequency are obtained by V-g method. The LCO behavior of the transonic control surface buzz system with linear structure exhibits subcritical or supercritical bifurcation at different Mach numbers. For nonlinear structural model with a free-play nonlinearity in the control surface deflection stiffness, the double LCO phenomenon is observed in certain range of flutter speed. The free-play nonlinearity changes the stability of LCOs at small amplitudes and turns the unstable LCO into a stable one. The LCO behavior is dominated by the aerodynamic nonlinearity for the case with large control surface oscillation amplitude but by the structural nonlinearity for the case with small amplitude. Good agreements between LCO behaviors obtained by the present method and available experimental data show that our study may help to explain the experimental observation in wind tunnel tests and to understand the physical mechanism of transonic control surface buzz.     

A theoretical investigation of the development of stationary crossflow vortices in the boundary layer on a swept wing

Crossflow instability plays very important role in the transition of the boundary layer on a swept wing, typical in the engineering applications. Experiments revealed that the linear stability theory well predicted the form of the crossflow vortices, but usually much overpredicted their growth rate. Using nonlinear theory of hydrodynamic stability, combined with some other considerations, we were able to obtain the growth rate in good agreement with experimental observations. The project supported by the National Natural Science Foundation of China, Grant No. 19572048     

Evolution of disturbances in a laminarized supersonic boundary layer on a swept wing

Experimental data on stability of a three-dimensional supersonic boundary layer on a swept wing are presented. The experiments are performed on a swept wing model with a lenticular profile with a 40° sweep angle of the leading edge at a zero angle of attack. The supersonic boundary layer on the swept wing was laminarized with the use of distributed roughness. A pioneering study of interaction of traveling and stationary disturbances is performed. Some specific features of this interaction are identified. The main reason for turbulence emergence in a supersonic boundary layer on a swept wing is demonstrated to be secondary crossflow instability. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 2, pp. 40–46, March–April, 2008.     

Flow-induced vibrations and the Landau equation

The initial and general transient temporal behaviours of flow induced vibrations were studied particularly with regard to the extent with which self-excited, flow-induced vibrations can be described by the Landau equation.Three different cases are studied experimentally, using bodies with generic shape. The first case represents a bluff body: resembling a simplified section of the Tacoma bridge, which has a single torsional degree of freedom. The second case is a two-dimensional airfoil in transonic flow having a heave and a torsional degree of freedom. The third case is an elastic half-wing model, also investigated in transonic flow.It is shown that in all three cases, beyond the critical point and at small initial amplitude the temporal development of the oscillations up to the limit cycle i.e. the envelope of the measured time functions, agrees with the corresponding curve, given by the solution of the Landau equation. For the Tacoma section and the airfoil the same agreement was demonstrated for initial amplitudes much larger than that of the limit cycle. In addition for both last cases the bifurcation behaviour was investigated and the Landau constants were determined. Finally an elementary physical explanation for the instability phenomenon was given.     

Transonic limit cycle oscillation behavior of an aeroelastic airfoil with free-play

The limit cycle oscillation (LCO) behaviors of an aeroelastic airfoil with free-play for different Mach numbers are studied. Euler equations are adopted to obtain the unsteady aerodynamic forces. Aerodynamic and structural describing functions are employed to deal with aerodynamic and structural nonlinearities, respectively. Then the flutter speed and flutter frequency are obtained by V-g method. The LCO solutions for the aeroelastic airfoil obtained by using dynamically linear aerodynamics agree well with those obtained directly by using nonlinear aerodynamics. Subsequently, the dynamically linear aerodynamics is assumed, and results show that the LCOs behave variously in different Mach number ranges. A subcritical bifurcation, consisting of both stable and unstable branches, is firstly observed in subsonic and high subsonic regime. Then in a narrow Mach number range, the unstable LCOs with small amplitudes turn to be stable ones dominated by the single degree of freedom flutter. Meanwhile, these LCOs can persist down to very low flutter speeds. When the Mach number is increased further, the stable branch turns back to be unstable. To address the reason of the stability variation for different Mach numbers at small amplitude LCOs, we find that the Mach number freeze phenomenon provides a physics-based explanation and the phase reversal of the aerodynamic forces will trigger the single degree of freedom flutter in the narrow Mach number range between the low and high Mach numbers of the chimney region. The high Mach number can be predicted by the freeze Mach number, and the low one can be estimated by the Mach number at which the aerodynamic center of the airfoil lies near its elastic axis. Influence of angle of attack and viscous effects on the LCO behavior is also discussed.     

Wave forerunners of streamwise structures in a swept wing boundary layer

In a model experiment wave packets (forerunners) have been detected for the first time in the flow regions preceding the fronts of streamwise structures in a swept wing boundary layer. The characteristics of the wave packets and the generating streamwise structures and the dynamics of their downstream development are investigated. It is shown that the forerunners can transform into turbulent spots thus leading to laminar-turbulent transition. Certain components of the forerunners are compared with a periodic instability wave.     

Numerical investigation of the effects of structural geometric and material nonlinearities on limit-cycle oscillation of a cropped delta wing

This article presents numerical simulations of the limit-cycle oscillation (LCO) of a cropped delta wing in order to investigate the effects of structural geometric and material nonlinearities on aeroelastic behavior. In the computational model, the structural part included both the geometric nonlinearity that arises from large deflections, and the material nonlinearity that originates from plasticity. The Euler equations were employed in the fluid part to describe the transonic aerodynamics. Moreover, the load transfer was conducted using a 3-D interpolating procedure, and the interfaces between the structural and aerodynamic domains were constructed in the form of an exact match. The flutter and LCO behaviors of the cropped delta wing were simulated using the coupling model, and the results were compared with existing experimental measurements. For lower dynamic pressures, the geometric nonlinearity provided the proper mechanism for the development of the LCO, and the numerical results correlated with the experimental values. For higher dynamic pressures, the material nonlinearity led to a rapid rise in the LCO amplitude, and the simulated varying trend was consistent with the experimental observation. This study demonstrated that the LCO of the cropped delta wing was not only closely related to geometric nonlinearity, but was also remarkably affected by material nonlinearity.     

Nonlinear aeroelastic analysis of a two-dimensional wing with control surface in supersonic flow

Based on the piston theory of supersonic flow and the energy method, a two dimensional wing with a control surface in supersonic flow is theoretically modeled, in which the cubic stiffness in the torsional direction of the control surface is considered. An approximate method of the cha- otic response analysis of the nonlinear aeroelastic system is studied, the main idea of which is that under the condi- tion of stable limit cycle flutter of the aeroelastic system, the vibrations in the plunging and pitching of the wing can approximately be considered to be simple harmonic excita- tion to the control surface. The motion of the control surface can approximately be modeled by a nonlinear oscillation of one-degree-of-freedom. The range of the chaotic response of the aeroelastic system is approximately determined by means of the chaotic response of the nonlinear oscillator. The rich dynamic behaviors of the control surface are represented as bifurcation diagrams, phase-plane portraits and PS diagrams. The theoretical analysis is verified by the numerical results.     

Effects of wing deformation on aerodynamic performance of a revolving insect wing

Flexible wings of insects and bio-inspired micro air vehicles generally deform remarkably during flapping flight owing to aerodynamic and inertial forces,which is of highly nonlinear fluid-structure interaction（FSI）problems.To elucidate the novel mechanisms associated with flexible wing aerodynamics in the low Reynolds number regime,we have built up a FSI model of a hawkmoth wing undergoing revolving and made an investigation on the effects of flexible wing deformation on aerodynamic performance of the revolving wing model.To take into account the characteristics of flapping wing kinematics we designed a kinematic model for the revolving wing in two-fold：acceleration and steady rotation,which are based on hovering wing kinematics of hawkmoth,Manduca sexta.Our results show that both aerodynamic and inertial forces demonstrate a pronounced increase during acceleration phase,which results in a significant wing deformation.While the aerodynamic force turns to reduce after the wing acceleration terminates due to the burst and detachment of leading-edge vortices（LEVs）,the dynamic wing deformation seem to delay the burst of LEVs and hence to augment the aerodynamic force during and even after the acceleration.During the phase of steady rotation,the flexible wing model generates more ver-tical force at higher angles of attack（40°–60°）but less horizontal force than those of a rigid wing model.This is because the wing twist in spanwise owing to aerodynamic forces results in a reduction in the effective angle of attack at wing tip,which leads to enhancing the aerodynamics performance by increasing the vertical force while reducing the horizontal force.Moreover,our results point out the importance of the fluid-structure interaction in evaluating flexible wing aerodynamics：the wing deformation does play a significant role in enhancing the aerodynamic performances but works differently during acceleration and steady rotation,which is mainly induced by inertial force in acceleration but by aerodynamic forces     

Experimental investigation on aerodynamic performance of a flapping wing vehicle in forward flight

The aerodynamic performance of a flexible membrane flapping wing has been investigated here. For this purpose, a flapping-wing system and an experimental set-up were designed to measure the unsteady aerodynamic forces of the flapping wing motion. A one-component force balance was set up to record the temporal variations of aerodynamic forces. The flapping wing was studied in a large low-speed wind tunnel. The lift and thrust of this mechanism were measured for different flapping frequencies, angles of attack and for various wind tunnel velocities. Results indicate that the thrust increases with the flapping frequency. An increase in the wind tunnel speed and flow angle of attack leads to reduction in the thrust value and increases the lift component. The aerodynamic and performance parameters were nondimensionalized. Appropriate models were introduced which show its aerodynamic performance and may be used in the design process and also optimization of the flapping wing.     

Reduced order nonlinear aeroelasticity of swept composite wings using compressible indicial unsteady aerodynamics

Nonlinear dynamic aeroelasticity of composite wings in compressible flows is investigated. To provide a reasonable model for the problem, the composite wing is modeled as a thin walled beam (TWB) with circumferentially asymmetric stiffness layup configuration. The structural model considers nonlinear strain displacement relations and a number of non-classical effects, such as transverse shear and warping inhibition. Geometrically nonlinear terms of up to third order are retained in the formulation. Unsteady aerodynamic loads are calculated according to a compressible model, described by indicial function approximations in the time domain. The aeroelastic system of equations is augmented by the differential equations governing the aerodynamics lag states to derive the final explicit form of the coupled fluid-structure equations of motion. The final nonlinear governing aeroelastic system of equations is solved using the eigenvectors of the linear structural equations of motion to approximate the spatial variation of the corresponding degrees of freedom in the Ritz solution method. Direct time integrations of the nonlinear equations of motion representing the full aeroelastic system are conducted using the well-known Runge–Kutta method. A comprehensive insight is provided over the effect of parameters such as the lamination fiber angle and the sweep angle on the stability margins and the limit cycle oscillation behavior of the system. Integration of the interpolation method employed for the evaluation of compressible indicial functions at any Mach number in the subsonic compressible range to the derivation process of the third order nonlinear aeroelastic system of equations based on TWB theory is done for the first time. Results show that flutter speeds obtained by the incompressible unsteady aerodynamics are not conservative and as the backward sweep angle of the wing is increased, post-flutter aeroelastic response of the wing becomes more well-behaved.     

Controlling the near vortex wake of a swept half-wing by means of mini-flaps

The effect of mini-flaps on the flow pattern in the near vortex wake behind a model swept half-wing is investigated. The distributions of the time-average flow velocity were measured in a subsonic wind tunnel, in a section normal to the freestream velocity vector located at a distance of 3.8 wing half-spans from its trailing edge. When mini-flaps are mounted on both upper and lower wing surfaces, two vortices (tip and auxiliary) of the same sign are observable in the above-mentioned flow section; they are separated by an extended region of vorticity of the opposite sign. The model angle-of-attack effect on the intensities of the tip and auxiliary vortices is estimated.     

扑翼柔性及其对气动特性的影响

以往对扑翼气动特性的研究基本上都是基于简单的匀速刚性模型,但是通过大量观察不同飞鸟的扑翼动作发现,该模型与鸟翼的实际扑动还有很大差别。鸟翼不但上扑段和下扑段所需时间不同,而且在扑动过程中,鸟翼的形状无论沿弦向或展向都存在着相当大的柔性变形。本文在原有匀速刚性模型的基础上,加入了扑动速率变化和形状变化的影响,得出新的变速柔性扑翼分析模型,使之更接近鸟翼柔性扑动的真实情况。通过对比计算发现,柔性变形对扑翼的升力与推力都有着显著影响,如果控制得当,柔性变形能大大改善扑翼的气动性能。     

Navier–Stokes computation of transonic vortices over a round leading edge delta wing

A 3D Navier–Stokes solver has been developed to simulate laminar compressible flow over quadrilateral wings. The finite volume technique is employed for spatial discretization with a novel variant for the viscous fluxes. An explicit three-stage Runge–Kutta scheme is used for time integration, taking local time steps according to the linear stability condition derived for application to the Navier–Stokes equations. The code is applied to compute primary and secondary separation vortices at transonic speeds over a 65° swept delta wing with round leading edges and cropped tips. The results are compared with experimental data and Euler solutions, and Reynolds number effects are investigated.     

外挂物对大展弦比直机翼颤振特性的影响

为提高带外挂物大展弦比直机翼的颤振速度,基于假设模态法提出一种带集中质量弯扭组合梁模态分析手段,结合片条理论考察外挂物不同质量及布置形式对机翼颤振特性的影响。首先,基于弯扭组合梁建立带外挂物大展弦比直机翼的结构动力学模型,并利用假设模态法得到其弯曲和扭转模态。其次,引入片条理论近似计算有限翼展升力面的气动力,调整外挂物的质量、数目及其在机翼展向和弦向的相对位置,得到外挂物对机翼颤振特性的影响规律。最后,利用在机翼前缘附近悬吊小质量外挂物可提高机翼颤振速度的优势,探究颤振速度恢复方法并提出颤振速度恢复的优化问题,使得携带外挂物的机翼与不携带外挂物的机翼颤振速度基本相同。研究结果表明,外挂物的不同悬挂方式可引起机翼颤振模态的跳转,在以俯仰为颤振主模态的机翼上可调整外挂物位置以恢复原机翼的颤振速度。     

On bimodal flutter behavior of a flexible airfoil

The dynamic aeroelastic behavior of an elastically supported airfoil is studied in order to investigate the possibilities of increasing critical flutter speed by exploiting its chord-wise flexibility. The flexible airfoil concept is implemented using a rigid airfoil-shaped leading edge, and a flexible thin laminated composite plate conformally attached to its trailing edge. The flutter behavior is studied in terms of the number of laminate plies used in the composite plate for a given aeroelastic system configuration. The flutter behavior is predicted by using an eigenfunction expansion approach which is also used to design a laminated plate in order to attain superior flutter characteristics. Such an airfoil is characterized by two types of flutter responses, the classical airfoil flutter and the plate flutter. Analysis shows that a significant increase in the critical flutter speed can be achieved with high plunge and low pitch stiffness in the region where the aeroelastic system exhibits a bimodal flutter behavior, e.g., where the airfoil flutter and the plate flutter occur simultaneously. The predicted flutter behavior of a flexible airfoil is experimentally verified by conducting a series of systematic aeroelastic system configurations wind tunnel flutter campaigns. The experimental investigations provide, for each type of flutter, a measured flutter response, including the one with indicated bimodal behavior.     

Subcritical, nontypical and period-doubling bifurcations of a delta wing in a low speed wind tunnel

Limit Cycle Oscillations (LCOs) involving Delta wings are an important area of research in modern aeroelasticity. Such phenomena can be the result of geometric or aerodynamic nonlinearity. In this paper, a flexible half-span Delta wing is tested in a low speed wind tunnel in order to investigate its dynamic response. The wing is designed to be more flexible than the models used in previous research on the subject in order to expand the airspeed range in which LCOs occur. The experiments reveal that this wing features a very rich bifurcation behavior. Three types of bifurcation are observed for the first time for such an aeroelastic system: subcritical bifurcations, period-doubling/period-halving and nontypical bifurcations. They give rise to a great variety of LCOs, even at very low angles of attack. The LCOs resulting from the nontypical bifurcation display Hopf-type behavior, i.e. having fundamental frequencies equal to one of the linear modal frequencies. All of the other LCOs have fundamental frequencies that are unrelated to the underlying linear system modes.     

Numerical study of vortex-breakdown over a slender wing in transonic flow at high incidence

The transonic flowfields and vortex-breakdown over a slender wing with the angle of attack from 10° to 28° are studied numerically, and the emphasis is on the secondary separation and the charateristics of vortex-breakdown. The results indicated that: (a) TVD schemes have strong capability for capturing vortices in three-dimensional transonic separated flow at large angle of attack. (b) The development of secondary vortices is more complex than that of leading-edge ones, and is affected by wing's configuration, angle of attack and compressibility simultaneously, and the effect of compressibility is more severe at low angle of attack. (c) The starting angle of attack for vortex-breakdown (when vortex bursting point crosses trailing-edge) is about 18° forM∞=0.85, then the bursting point moves upstream quickly with increasing angle of attack. (d) At α=24°, breakdown occurs over most part of upper side, and the wing begins to stall. Therefore, there is a large lag of angle of attack between the beginning of vortex-breakdown and the stall of the wing. (e) This lag increase with the decreasing of Mach number.     

Viscous and inviscid interactions of an oblique shock with a flexible panel

The complex self-sustained oscillations arising from the interaction of an oblique shock with a flexible panel in both the inviscid and viscous regimes have been investigated numerically. The aeroelastic interactions are simulated using either the Euler or the full compressible Navier–Stokes equations coupled to the nonlinear von Karman plate equations. Results demonstrate that for a sufficiently strong shock limit-cycle oscillations emerge from either subcritical or supercritical bifurcations even in the absence of viscous separated flow effects. The critical dynamic pressure diminishes with increasing shock strength and can be much lower than that corresponding to standard panel flutter. Significant changes in panel dynamics were also found as a function of the shock impingement point and cavity pressure. For viscous laminar flow above the panel without a shock, high-frequency periodic oscillations appear due to the coupling of boundary-layer instabilities with high-mode flexural deflections. For a separated shock laminar boundary layer interaction, non-periodic self-excited oscillations arise which can result in a significant reduction in the extent of the time-averaged separation region. This finding suggests the potential use of an aeroelastically tailored flexible panel as a means of passive flow control. Forced panel oscillations, induced by a specified variable cavity pressure underneath the panel, were also found to be effective in reducing separation. For both inviscid and viscous interactions, the significant unsteadiness generated by the fluttering panel propagates along the complex reflected expansion/recompression wave system.