LIMIT-CYCLE OSCILLATIONS IN HIGH-ASPECT-RATIO WINGS

The paper presents results obtained for limit-cycle oscillations (LCOs) in high-aspect-ratio wings caused by structural and aerodynamic nonlinearities. The analysis is based on geometrically exact structural analysis and finite-state unsteady aerodynamics with stall. The results indicate that stall limits the amplitude of post-flutter unstable oscillations. At speeds below the linear flutter speed, LCOs can be observed if the stable steady state is disturbed by a finite-amplitude disturbance. A critical disturbance magnitude is required at a given speed and a critical speed is required at a given disturbance magnitude to initiate LCOs. The LCO initiation mechanism can be attributed to the change in structural characteristics of the wing with deformation. It is also observed that the LCO gets increasingly complex with increasing speed. Period doubling is observed at low speeds and as the speed increases the oscillations lose periodicity and become chaotic.     

The dynamics of two-dimensional cantilevered plates with an additional spring support in axial flow

In this paper, the dynamics of two-dimensional cantilevered flexible plates in axial flow is investigated using a fluid–structure interaction model. An additional spring support of either linear or cubic type is installed at various locations on the plate; its presence qualitatively affects the dynamics of the fluid–structure system. Without the spring, the cantilevered plate loses stability by flutter when the flow velocity exceeds a critical value; as the flow velocity increases further, the system dynamics is qualitatively the same: the plate undergoes symmetric limit cycle oscillations with increasing amplitude. With a linear spring, a state of static buckling is added to the dynamics. Rich nonlinear dynamics can be observed when a cubic spring is considered; the plate may be stable and buckled, and it may undergo either symmetric or asymmetric limit cycle oscillations. Moreover, when the flow velocity is sufficiently high, the plate may exhibit chaotic motions via a period-doubling route.     

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.     

Intermittency in a cantilever plate in a randomly fluctuating fluid flow

Fluid–elastic systems nearing dynamic instabilities are known to be sensitive to fluctuations in fluid flow. A cantilever plate in axial flow with random temporal fluctuations, is examined numerically for its dynamical behaviour. The numerical model comprises of a nonlinear structural model for the flexible plate, coupled with unsteady lumped vortex model for the fluid forces. As the mean flow velocity is increased, the system transitions to limit cycle oscillations from a state of rest, through a regime of intermittent oscillations. The conditions for onset and disappearance of intermittency are discussed and are interpreted using stochastic bifurcation theories. While the onset of intermittency is found to be unaffected by the time scales of the flow fluctuations, they are observed to affect the length of the intermittency regime. The effect of plate flexibility on intermittency is also discussed.     

Aeroelastic suppression of an airfoil with control surface using nonlinear energy sink

Aeroelasticity exists in airfoil with control surface freeplay, which may induce instability in an incompressible flow. In this paper, a nonlinear energy sink (NES) is used to suppress the aeroelasticity of an airfoil with a control surface. The freeplay and cubic nonlinearity in pitch are taken into account. The harmonic balance method is used to analytically determine the limit cycle oscillations (LCOs) amplitudes of the airfoil–NES system. Linear and nonlinear flutter speeds are detected from the airfoil with control surface freeplay. When NES is attached, both the linear flutter speed of airfoil without freeplay and the nonlinear flutter speed of airfoil with a freeplay are increased. Moreover, the LCO amplitude of airfoil is decreased due to NES. Then, the influences of NES parameters on the increase in flutter boundary of airfoil are carefully studied.     

Computational aeroelastic simulations of self-sustained pitch oscillations of a NACA0012 at transitional Reynolds numbers

The phenomenon of low amplitude self-sustained pitch oscillations in the transitional Reynolds number regime is studied numerically through unsteady, two-dimensional aeroelastic simulations. Based on the experimental data, simulations have been limited in the Reynolds number range 5.0×10⁴<Re_c<1.5×10⁵. Both laminar and URANS calculations (using the SST k–ω model with a low-Reynolds-number correction) have been performed and found to produce reasonably accurate limit cycle pitching oscillations (LCO). This investigation confirms that the laminar separation of the boundary layer near the trailing edge plays a critical role in initiating and sustaining the pitching oscillations. For this reason, the phenomenon is being labelled as laminar separation flutter. As a corollary, it is also shown that turbulence tends to inhibit their existence. Furthermore, two regimes of LCO are observed, one where the flow is laminar and separated without re-attachment, and the second for which transition has occurred followed by turbulent re-attachment. Finally, it is established that the high-frequency, shear instabilities present in the flow which lead to von Kármán vortex shedding are not crucial, nor necessary, to the maintaining mechanism of the self-sustained oscillations.     

Limit cycle oscillations of rectangular cantilever wings containing cubic nonlinearity in an incompressible flow

Limit cycle oscillations (LCO) as well as nonlinear aeroelastic analysis of rectangular cantilever wings with a cubic nonlinearity are investigated. Aeroelastic equations of a rectangular cantilever wing with two degrees of freedom in an incompressible potential flow are presented in the time domain. The harmonic balance method is modified to calculate the LCO frequency and amplitude for rectangular wings. In order to verify the derived formulation, flutter boundaries are obtained via a linear analysis of the derived system of equations for five different cases and compared with experimental data. Satisfactory results are gained through this comparison. The problem of finding the LCO frequency and amplitude is solved via applying the two methods discussed for two different cases with hardening cubic nonlinearities. The results from first-, third- and fifth-order harmonic balance methods are compared with the results of an exact numerical solution. A close agreement is obtained between these harmonic balance methods and the exact numerical solution of the governing aeroelastic equations. Finally, the nonlinear aeroelastic analysis of a rectangular cantilever wing with a softening nonlinearity is studied.     

Nonlinear dynamics of an aeroelastic airfoil with free-play in transonic flow

Nonlinear dynamic behaviors of an aeroelastic airfoil with free-play in transonic air flow are studied. The aeroelastic response is obtained by using time-marching approach with computational fluid dynamics (CFD) and reduced order model (ROM) techniques. Several standardized tests of transonic flutter are presented to validate numerical approaches. It is found that in time-marching approach with CFD technique, the time-step size has a significant effect on the calculated aeroelastic response, especially for cases considering both structural and aerodynamic nonlinearities. The nonlinear dynamic behavior for the present model in transonic air flow is greatly different from that in subsonic regime where only simple harmonic oscillations are observed. Major features of the responses in transonic air flow at different flow speeds can be summarized as follows. The aeroelastic responses with the amplitude near the free-play are dominated by single degree of freedom flutter mechanism, and snap-though phenomenon can be observed when the air speed is low. The bifurcation diagram can be captured by using ROM technique, and it is observed that the route to chaos for the present model is via period-doubling, which is essentially caused by the free-play nonlinearity. When the flow speed approaches the linear flutter speed, the aeroelastic system vibrates with large amplitude, which is dominated by the aerodynamic nonlinearity. Effects of boundary layer and airfoil profile on the nonlinear responses of the aeroelastic system are also discussed.     

Extremely large-amplitude oscillation of soft pipes conveying fluid under gravity

In this work, the nonlinear behaviors of soft cantilevered pipes containing internal fluid flow are studied based on a geometrically exact model, with particular focus on the mechanism of large-amplitude oscillations of the pipe under gravity. Four key parameters, including the flow velocity, the mass ratio, the gravity parameter, and the inclination angle between the pipe length and the gravity direction, are considered to affect the static and dynamic behaviors of the soft pipe. The stability analyses show that, provided that the inclination angle is not equal to π, the soft pipe is stable at a low flow velocity and becomes unstable via flutter once the flow velocity is beyond a critical value. As the inclination angle is equal to π, the pipe experiences, in turn,buckling instability, regaining stability, and flutter instability with the increase in the flow velocity. Interestingly, the stability of the pipe can be either enhanced or weakened by varying the gravity parameter, mainly dependent on the value of the inclination angle.In the nonlinear dynamic analysis, it is demonstrated that the post-flutter amplitude of the soft pipe can be extremely large in the form of limit-cycle oscillations. Besides,the oscillating shapes for various inclination angles are provided to display interesting dynamical behaviors of the inclined soft pipe conveying fluid.     

A theoretical and experimental investigation of the effects of a steady angle of attack on the nonlinear flutter of a delta wing plate model

Limit cycle oscillations (LCO) of wings on certain modern high performance aircraft have been observed in flight and in wind tunnel experiments. Whether the physical mechanism that gives rise to this behavior is a fluid or structural nonlinearity or both is still uncertain. It has been shown that an aeroelastic theoretical model with only a structural nonlinearity can predict accurately the limit cycle behavior at low subsonic flow for a plate-like wing at zero angle of attack. Changes in the limit cycle and flutter behavior as the angle of attack is varied have also been observed in flight. It has been suggested that this sensitivity to angle of attack is due to a fluid nonlinearity. In this investigation, we study the flutter and limit cycle behavior of a wing in low subsonic flow at small steady angles of attack. Experimental results are compared to those predicted using an aeroelastic theoretical model with only a structural nonlinearity. Results from both experiment and theory show a change in flutter speed as the steady angle of attack is varied. Also the LCO magnitude increased at a given velocity as the angle of attack was increased for both the experiment and theory. While not proving that the observed sensitivity to angle of attack of LCO in aircraft is due to a structural nonlinearity, the results do show that a change in the aeroelastic behavior at angles of attack can be caused by a structural nonlinearity as well as a fluid nonlinearity. In this paper, only structural nonlinearities are considered, but an extension to include aerodynamic nonlinearities would be very worthwhile.     

Nonlinear Analysis of Store-Induced Limit Cycle Oscillations

Flight tests of modern high-performance fighter aircraft reveal the presence of limit cycle oscillation (LCO) responses for aircraft with certain external store configurations. Conventional linear aeroelastic analysis predicts flutter for conditions well beyond the operational envelope, yet these store-induced LCO responses occur at flight conditions within the flight envelope. Several nonlinear sources may be present, including aerodynamic effects such as flow separation and shock-boundary layer interaction and structural effects such as stiffening, damping, and system kinematics. No complete theory has been forwarded to accurately explain the mechanisms responsible. This research examines a two degree-of-freedom aeroelastic system which possesses kinematic nonlinearities and a strong nonlinearity in pitch stiffness. Nonlinear analysis techniques are used to gain insight into the characteristics of the behavior of the system. Numerical simulation is used to verify and validate the analysis. It is found that when system damping is low, the system clearly exhibits nonlinear interaction between aeroelastic modes. It is also shown that although certain applied forcing conditions may appear negligible, these same forces produce large amplitude LCOs under specific realizable circumstances.     

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.     

Flow field around the flapping flag

The flapping flag is a canonical fluid–structure interaction problem that describes a cantilever plate with flow along its elastic axis. When the flapping flag loses stability it enters a large amplitude Limit Cycle Oscillation (LCO). While theoretical models can accurately predict the flutter velocity and frequency, there are still discrepancies between the experimental observations and the theoretical predictions of the post-critical LCO response. This note provides recent flow field visualizations in a single longitudinal plane for a cantilevered aluminum plate in axial flow during its LCO. Particle Image Velocimetry (PIV) techniques are used to show that the flow over the midspan of the plate is attached even during the violent LCO motion. This observation suggests that potential flow aerodynamic models may be able to capture the essential features in the flow field.     

Amplification and amplitude limitation of heave/pitch limit-cycle oscillations close to the transonic dip

Recent results from flutter experiments of the supercritical airfoil NLR 7301 at flow conditions close to the transonic dip are presented. The airfoil was mounted with two degrees-of-freedom in an adaptive solid-wall wind tunnel, and boundary-layer transition was tripped. Flutter boundaries exhibiting a transonic dip were determined and limit-cycle oscillations (LCOs) were measured. The local energy exchange between the fluid and the structure during LCOs is examined and leads to the following findings: at supercritical Mach numbers below that of the transonic-dip minimum the presence of a shock-wave and its dynamics destabilizes the aeroelastic system such that the decreasing branch of the transonic dip develops. At higher Mach numbers the shock-wave motion has a stabilizing effect such that the flutter boundary increases to higher flutter-speed indices with increasing Mach number. Amplified oscillations near this branch of the flutter boundary obtain energy from the flow mainly due to the dynamics of a trailing-edge flow separation. A slight nonlinear amplitude dependency of the shock motion and a possibly occurring boundary-layer separation cause the amplitude limitation of the observed LCOs. The impact of the findings on the numerical simulation of these phenomena is discussed.     

Articulated Pipes Conveying Fluid Pulsating with High Frequency

Stability and nonlinear dynamics of two articulated pipes conveying fluid with a high-frequency pulsating component is investigated. The non-autonomous model equations are converted into autonomous equations by approximating the fast excitation terms with slowly varying terms. The downward hanging pipe position will lose stability if the mean flow speed exceeds a certain critical value. Adding a pulsating component to the fluid flow is shown to stabilize the hanging position for high values of the ratio between fluid and pipe-mass, and to marginally destabilize this position for low ratios. An approximate nonlinear solution for small-amplitude flutter oscillations is obtained using a fifth-order multiple scales perturbation method, and large-amplitude oscillations are examined by numerical integration of the autonomous model equations, using a path-following algorithm. The pulsating fluid component is shown to affect the nonlinear behavior of the system, e.g. bifurcation types can change from supercritical to subcritical, creating several coexisting stable solutions and also anti-symmetrical flutter may appear.     

Flow-induced oscillation of two circular cylinders in tandem arrangement at low Re

Results are presented for flow-induced vibrations of a pair of equal-sized circular cylinders of low nondimensional mass (m^*=10) in a tandem arrangement. The cylinders are free to oscillate both in streamwise and transverse directions. The Reynolds number, based on the free-stream speed and the diameter of the cylinders, D is 100 and the centre-to-centre distance between the cylinders is 5.5D. The computations are carried out for reduced velocities in the range 2≤U^*≤15. The structural damping is set to zero for enabling maximum amplitudes of oscillation. A stabilized finite element method is utilized to carry out the computations in two dimensions. Even though the response of the upstream cylinder is found to be qualitatively similar to that of an isolated cylinder, the presence of a downstream cylinder is found to have significant effect on the behaviour of the upstream cylinder. The downstream cylinder undergoes very large amplitude of oscillations in both transverse and streamwise directions. The maximum amplitude of transverse response of the downstream cylinder is quite similar to that of a single cylinder at higher Re beyond the laminar regime. Lock-in and hysteresis are observed for both upstream and downstream cylinders. The downstream cylinder undergoes large amplitude oscillations even beyond the lock-in state. The phase between transverse oscillations and lift force suffers a 180 jump for both the cylinders almost in the middle of the synchronization regime. The phase between the transverse response of the two cylinders is also studied. Complex flow patterns are observed in the wake of the freely vibrating cylinders. Based on the phase difference and the flow patterns, the entire flow range is divided into five sub-regions.     

Nonlinear dynamics of extensible fluid-conveying pipes,supported at both ends

In this paper, the post-divergence behaviour of extensible fluid-conveying pipes supported at both ends is studied using the weakly nonlinear equations of motion of Semler, Li and Païdoussis. The two coupled nonlinear partial differential equations are discretized via Galerkin's method and the resulting set of ordinary differential equations is solved either by Houbolt's finite difference method or via AUTO. Typically, the pipe is stable at its original static equilibrium position up to the flow velocity where it loses stability by static divergence via a supercritical pitchfork bifurcation. The amplitude of the resultant buckling increases with increasing flow, but no secondary instability occurs beyond the pitchfork bifurcation. The effects of the system parameters on pipe behaviour as well as the possibility of a subcritical pitchfork bifurcation have also been studied.     

Aeroelastic behavior of a flag in ground effect

The aeroelastic behavior of a flexible plate subjected to a uniform axial flow is investigated in the presence of a rigid plane set parallel to the plate. It is shown that the ground effect reduces the flutter inflow velocity and strengthens the possibility of using the flag for extracting energy from winds and currents. The numerical analysis is carried out assuming that both the unsteady potential incompressible flow and the plate can be described with 2D models, i.e., a lumped vortex panel method and a nonlinear Euler–Bernoulli beam model, respectively, without losing the essential features of the fluid–structure interaction. Asymmetry of post-critical behavior (limit-cycle oscillations) and sensitivity of the results to the main flag parameters (distance from the ground, mass ratio and damping) are also considered, including also the energy distribution over the identified proper orthogonal modes. The investigated reduction of the flutter velocity in ground effect has been also confirmed with experimental tests relative to a polypropylene flag with and without the rigid panel mimicking the presence of the ground.     

Characterization of the gas pulse frequency,amplitude and velocity in non-steady dense phase pneumatic conveying of powders

Current modelling techniques for the prediction of conveying line pressure drop in low velocity dense phase pneumatic conveying are largely based on steady state analyses. Work in this area has been on-going for many years with only marginal improvements in the accuracy of prediction being achieved. Experimental and theoretical investigations undertaken by the authors suggest that the flow mechanisms involved in dense phase conveying are dominated by transient effects rather than those of steady state and are possibly the principal reasons for the limited improvement in accuracy. This paper reports on investigations on the pressure fluctuation behaviour in dense phase pneumatic conveying of powders. The pressure behaviour of the gas flow in the top section of the pipeline was found to exhibit pulsatile oscillations. In particular, the pulse velocity showed variation in magnitude while the frequency of the oscillations rarely exceeded 5 Hz. A wavelet analysis using the Daubechie 4 wavelet found that the amplitude of the oscillations increased along the pipeline. Furthermore, there was significant variation in gas pulse amplitude for different types of particulate material.     

超音速气流中受热曲壁板的非线性颤振特性

基于von Karman 大变形理论及带有曲率修正的一阶活塞理论, 用Galerkin方法建立了超音速气流中受热二维曲壁板的非线性气动弹性运动方程; 采用牛顿迭代法计算得到由静气动载荷和热载荷引起的静气动弹性变形; 根据李雅谱诺夫间接法分析了壁板初始曲率与温升对颤振边界的影响; 对二维曲壁板的非线性气动弹性方程组进行数值积分求解,分析了动压参数对受热二维曲壁板分岔特性的影响, 给出了典型状态下曲壁板非线性颤振响应的时程图与相图. 分析结果表明对小初始曲率的曲壁板, 温升对其静气动弹性变形影响较大, 且随着温升的增加其颤振临界动压急剧减小; 对具有较大初始曲率的曲壁板, 温升对其静气动弹性变形的影响较弱, 且随着温升的增加颤振临界动压基本保持不变. 初始几何曲率与气动热效应使得曲壁板具有复杂的动力学特性, 不再像平壁板一样, 经过倍周期分岔进入混沌, 而会出现由静变形状态直接进入混沌运动的现象, 且在混沌运动区域中还会出现静态稳定点或谐波运动, 在大曲率情况下, 曲壁板不会产生混沌运动, 而是幅值在一定范围内的极限带振荡.