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 parametric investigation of the propulsion of 2D chordwise-flexible flapping wings at low Reynolds number using numerical simulations

This paper presents a numerical investigation of the effects of chordwise flexibility on flapping wings at low Reynolds number. The numerical simulations are performed with a partitioned fluid–structure interaction algorithm using artificial compressibility stabilization. The choice of the structural dimensionless parameters is based on scaling arguments and is compared against parameters used by other authors. The different regimes, namely inertia-driven and pressure-driven wing deformations, are presented along with their effects on the topology of the flow and on the performance of a heaving and pitching flapping wing in propulsion regime. It is found that pressure-driven deformations can significantly increase the thrust efficiency if a suitable amount of flexibility is used. Significant thrust increases are also observed in zero pitching amplitude cases. The effects of the second and third deformation modes on the performances of pressure-driven deformation cases are discussed. On the other hand, inertia-driven deformations generally deteriorate aerodynamic performances of flapping wings unless the behavior of the wing deformation is modified by the presence of sustainable superharmonics in a way that produces slight improvements. It is also shown that wing flexibility can act as an efficient passive pitching mechanism that allows fair thrust and better efficiency to be achieved when compared to a rigid pitching–heaving wing.     

Velocity measurement of compressible air flows utilizing a high-speed video camera

PIV-measurements of the flow above a pitching airfoil were conducted in a transonic wind tunnel. An ultra high-speed video camera was used for separate recording of two exposures. The data was analysed using the cross-correlation method. The results show the applicability of the technique in high speed flows.     

On the role of thickness ratio and location of axis of rotation in the flat plate motions

Wind tunnel experiments were performed to characterize the flow-induced rotations and pitching of various flat plates as a function of the thickness ratio and the location of the axis of rotation. High-resolution telemetry, laser tachometer, and hotwire were used to get time series of the plates motions and the signature of the wake flow at a specific location. Results show that small axis offset can induce high-order modes in the plate rotation due to torque unbalance, and can trigger self-initiated pitching. The spectral decomposition of the flow velocity in the plate wake reveals the existence of a dominating high-frequency mode that corresponds to a static-like vortex shedding occurring at the maximum plate pitch. The associated characteristic length scale is the projected width at maximum pitching angle. The increase of the plate thickness ratio implies lower angular velocity in rotation cases. A simple model based on aerodynamic forces is used to explain the linear relation between pitching frequency and wind speed, the pitching frequency increase with axis offset, and the onset of pitching.     

Aerodynamic forces acting on an albatross flying above sea-waves

IntroductionIt was observed long time ago that the large sea bird such as albatross has excellent flightskills.The albatross can keep its soaring for a long time even in rough sea without flappingwings.This problemwas discussed byWood C.J.[1].He supposed …     

跨音速翼型上的激波／边界层干扰自适应控制计算

在激波区使用自适应壁对跨音速翼型的激波／边界层的相互作用（干扰）进行控制，可改变机翼的气动性能，这种被动控制可通过在翼型的激波区开一凹腔，其上覆盖一弹性橡胶膜柔壁来，本文给出用Ｎａｖｉｅｒ－Ｓｔｏｋｅｒ方程数值模拟这一自适应控制翼型的跨音速粘性绕流，提出了一个适应于本特殊情况（物面边界局部地区在求解过程中有变化）的处理办法。并探讨了自适应柔壁对当代跨音速翼绕流的影响。     

An efficient approach for calculation of pitching moment in nonlinear reduced frequency method at low Mach number transonic flows

In this research, an efficient methodology for calculation of pitching moment coefficient at low Mach number transonic flows by using the perturbed nonlinear reduced frequency approach is presented. The proposed approach uses the perturbation technique in the nonlinear frequency domain (NLFD) method to estimate the solution at high harmonics. In this approach, the density and velocity fields at high harmonics are perturbed about those at low harmonics. Perturbing the density and velocity fields, the semi‐linear form of the governing equations is obtained. The resulting solution vector and spatial operator are then approximated by discrete form of Fourier transformation and governing equations are solved by using the pseudo‐spectral approach. Numerical results show that the proposed approach predicts good pitching moment coefficient at low Mach number transonic flows with up to 50% savings in computational time. Copyright © 2011 John Wiley & Sons, Ltd.     

Comparative Analysis of the Characteristics of the Vortex Wake behind a Flapping Wing Performing Oscillations of Different Types

The wake characteristics of a custom-designed airfoil performing pitching oscillations, heaving oscillations, and a combination of pitch and heave oscillations are compared in this study. The influence of flapping parameters are investigated at a constant Reynolds number Re\(_{c} = 2640\) and is presented for the Strouhal numbers based on the oscillation amplitude, St_A, varying in the \(0.1 \leqslant {\text{S}}{{{\text{t}}}_{A}} \leqslant 0.4\) range. The generation of vorticity above and below the airfoil depends on the airfoil’s initial direction of motion and remains the same for all types of flapping oscillations investigated. The evolution of the leading-edge and trailing-edge vortices is presented. The heaving oscillations of the airfoil are found to have a greater influence on the characteristics of the leading edge vortex. The wake behind the combined pitch-heave oscillations appears to be governed by pitching oscillations below \({\text{S}}{{{\text{t}}}_{A}} = 0.24\), whereas it is driven by heaving oscillations above \({\text{S}}{{{\text{t}}}_{A}} = 0.24\). The force computations indicate that the mere existence of the reverse von Kármán street is not sufficient to develop the thrust on the airfoil. The periodic component of velocity fluctuations significantly influences the wake characteristics. The anisotropic stress field developed around the airfoil due to the periodic fluctuations of the velocity is presented. The coherent structures developed in the wake are identified using the proper orthogonal decomposition and a qualitative comparison of the structures for different flapping oscillations is presented. The energy transfer from the flapping airfoil to the fluid for different flapping oscillations is highest for heaving oscillations followed by combined pitch-heave oscillations and pitching oscillations.

     

Effect of flow regime change from subsonic to transonic on the air loads of an oscillating airfoil

In this research, the effect of flow regime change from subsonic to transonic on the air loads of a pitching NACA0012 airfoil is investigated. To do this, the effect of change in flow regime on the lift and pitching moment coefficients hysteresis cycles is studied. The harmonic balance approach is utilized for numerical calculation due to its low computational time. Verifications are also made with previous works and good agreements are observed. The assessment of flow regime change on the aforementioned hysteresis cycles is accomplished in the Mach number range of M=0.65–0.755. The reduced frequency and pitch amplitude also vary from k=0.03 to 0.1 and α₀=1–2.51°, respectively. Results show that the effect of increase in Mach number is to increase and decrease the lift coefficient during downstroke and upstroke, respectively, whereas at low reduced frequencies, the effect of increase in Mach number may lead to a reverse manner when airfoil moves toward its extremum angle of attack. Results also reveal that as the pitch amplitude varies, the shape of lift coefficient hysteresis cycle depends more on the pitch amplitude than on the appearance of shock. It is shown that as the Mach number increases, the incidence angles correspond to the extremum pitching moment, and depending on the reduced frequency, lie between zero and extremum angle of attack. These incidence angles shift toward the extremum angle of attack as the reduced frequency decreases. Results also show that the increase in pitch amplitude at low Mach number, in such a way that leads to the formation of shock around the extremum angle of attack, causes the extremum pitching moment to appear around these angles and at high Mach number, depending on the reduced frequency, the extremum pitching moment incidence angles would be between zero and extremum incidence angle.     

Energy harvesting of a heaving and forward pitching wing with a passively actuated trailing edge

This study experimentally investigates the energy harvesting capabilities of an oscillating wing with a passively actuated trailing edge. The oscillation kinematics are composed of a combined heaving and forward pitching motions, where the pitching axis is well behind the wing center of mass. Passive actuation is attained by connecting the trailing edge with the wing body using a torsion rod. The degree of flexibility of the trailing edge is represented by the Strouhal number based on the trailing edge natural frequency. The trailing edge passive response is studied for oscillation Strouhal numbers of 0.017, 0.025 and 0.033. Instantaneous aerodynamic forces are measured in a closed loop wind tunnel at a Reynolds number of 40 000, based on the free stream velocity and the wing chord length. Measured results include the effective angle of attack induced by the trailing edge actuation as well as the lift and moment during the oscillation cycle. For the imposed kinematics in this study, the pitching motion has a positive contribution to the mean power output whereas the heaving motion has a relatively small but negative contribution. Additionally, by decreasing the natural frequency of the trailing edge closer to that of the imposed oscillation frequency, the magnitude of the lift and moment forces and hence the mean power output, increases. It is found that there exists a strong correlation between mean power output and the effective angle of attack, shown through the passive trailing edge response, resulting in an increase in energy harvesting potential.     

Performance of a flapping foil flow energy harvester in shear flows

Through numerical simulations, we investigate the energy harvesting performance of a heaving/pitching foil in shear flow. With two-dimensional Navier–Stokes simulations, we examined the energy harvesting efficiencies of such a system in linear shear flows and compared the results with those in uniform flows. It is found that in low shear rates, the performance of the system in linear shear flow is slightly higher than that in uniform flow, whereas the energy harvesting efficiency is greatly diminished if the shear rate is sufficiently high (this effect is more pronounced in higher frequencies). This is attributed to the effects of linear shear on the vorticity generation and the synchronization between fluid forcing and foil motion – when a strong shear flow is introduced the lift force induced by the leading edge vortex that is in phase with the heaving motion of the foil is diminished. Furthermore, by studying the instability of the wake behind the foil, we confirm that the optimal performance of the foil in linear shear flows is associated with the same physical mechanism that controls the performance of the foil in uniform flows, i.e. the excitation of the most unstable modes in the wake when the oscillation frequency of the foil is close to the frequencies of these modes.     

Profiling the Supersonic Part of a Plug Nozzle with a Nonuniform Transonic Flow

The effect of transonic flow nonuniformity on the profiling of optimal plug nozzles is studied in the inviscid gas approximation. Sonic and supersonic regions providing maximum thrust for given nozzle dimensions and a given outer pressure are designed for given subsonic contours and calculated nonuniform transonic flows. As in the case of uniform flow on a cylindrical sonic surface, the initial regions of the designed contours satisfy the condition that in these regions the flow Mach number is unity or near-unity. In all the examples calculated, the optimal plug nozzles produce a greater thrust than the optimal axisymmetric and annular nozzles with a near-axial flow for the same lengths and the same gas flow rates through the nozzle. It is established that contouring without regard for transonic flow nonuniformity can result in considerable thrust losses. However, these losses are due only to a decrease in the flow rate, while the specific thrust may even increase slightly.     

Biomimetic propulsion under random heaving conditions,using active pitch control

Marine mammals travel long distances by utilizing and transforming wave energy to thrust through proper control of their caudal fin. On the other hand, manmade ships traveling in a wavy sea store large amounts of wave energy in the form of kinetic energy for heaving, pitching, rolling and other ship motions. A natural way to extract this energy and transform it to useful propulsive thrust is by using a biomimetic wing. The aim of this paper is to show how an actively pitched biomimetic wing could achieve this goal when it performs a random heaving motion. More specifically, we consider a biomimetic wing traveling with a given translational velocity in an infinitely extended fluid and performing a random heaving motion with a given energy spectrum which corresponds to a given sea state. A formula is invented by which the instantaneous pitch angle of the wing is determined using the heaving data of the current and past time steps. Simulations are then performed for a biomimetic wing at different heave energy spectra, using an indirect Source-Doublet 3-D–BEM, together with a time stepping algorithm capable to track the random motion of the wing. A nonlinear pressure type Kutta condition is applied at the trailing edge of the wing. With a mollifier-based filtering technique, the 3-D unsteady rollup pattern created by the random motion of the wing is calculated without any simplifying assumptions regarding its geometry. Calculated unsteady forces, moments and useful power, show that the proposed active pitch control always results in thrust producing motions, with significant propulsive power production and considerable beneficial stabilizing action to ship motions. Calculation of the power required to set the pitch angle prove it to be a very small percentage of the useful power and thus making the practical application of the device very tractable.     

The Multi-point Optimization of Shock Control Bump with Constant-Lift Constraint Enhanced with Suction and Blowing for a Supercritical Airfoil

Both shock control bump (SCB) and suction and blowing are flow control methods used to control the shock wave/boundary layer interaction (SWBLI) in order to reduce the resulting wave drag in transonic flows. A SCB uses a small local surface deformation to reduce the shock-wave strength, while suction decreases the boundary-layer thickness and blowing delays the flow separation. Here a multi-point optimization method under a constant-lift-coefficient constraint is used to find the optimum design of SCB and suction and blowing. These flow control methods are used separately or together on a RAE-2822 supercritical airfoil for a wide range of off-design transonic Mach numbers. The RANS flow equations are solved using Roe’s averages scheme and a gradient-based adjoint algorithm is used to find the optimum location and shape of all devices. It is shown that the simultaneous application of blowing and SCB (hybrid blowing/SCB) improves the average aerodynamic efficiency at off-design conditions by 18.2 % in comparison with the clean airfoil, while this increase is only 16.9 % for the hybrid suction/SCB. We have also studied the SWBLI and how the optimization algorithm makes the flow wave structure and interactions of the shock wave with the boundary layer favorable.     

A local domain‐free discretization method for simulation of incompressible flows over moving bodies

This paper presents a local domain‐free discretization (DFD) method for the simulation of unsteady flows over moving bodies governed by the incompressible Navier–Stokes equations. The discretization strategy of DFD is that the discrete form of partial differential equations at an interior point may involve some points outside the solution domain. All the mesh points are classified as interior points, exterior dependent points and exterior independent points. The functional values at the exterior dependent points are updated at each time step by the approximate form of solution near the boundary. When the body is moving, only the status of points is changed and the mesh can stay fixed. The issue of ‘freshly cleared nodes/cells’ encountered in usual sharp interface methods does not pose any particular difficulty in the presented method. The Galerkin finite‐element approximation is used for spatial discretization, and the discrete equations are integrated in time via a dual‐time‐stepping scheme based on artificial compressibility. In order to validate the present method for moving‐boundary flow problems, two groups of flow phenomena have been simulated: (1) flows over a fixed circular cylinder, a harmonic in‐line oscillating cylinder in fluid at rest and a transversely oscillating cylinder in uniform flow; (2) flows over a pure pitching airfoil, a heaving–pitching airfoil and a deforming airfoil. The predictions show good agreement with the published numerical results or experimental data. Copyright © 2010 John Wiley & Sons, Ltd.     

风力机叶片翼型动态试验技术研究

风力机叶片动态振荡过程往往伴随着俯仰和横摆同时进行, 以前对许多动态问题不清楚的阶段, 工程上不惜以增加叶片重量为代价而采用偏安全的设计, 通常忽略横摆振荡的影响; 大型风力机设计对获取翼型更加全面、准确的动态载荷提出了更高要求, 研究横摆振荡对翼型动态气动特性的影响规律具有重要意义. 本文首次开展翼型横摆振荡动态风洞试验研究, 采用“电子凸轮”技术代替机械凸轮实现了振荡频率和振荡角度的无级变化, 基于设计的电子外触发装置实现了对动态流场的实时测量, 实现了风洞来流、模型角位移和动态压力数据的同步采集, 分别开展了翼型静态测压、俯仰/横摆动态测压、粒子图像测速和荧光丝线等试验研究, 试验结果准度较高、规律合理; 分析了动态试验洞壁干扰影响机制. 研究表明, 横摆振荡翼型的气动曲线也存在明显迟滞效应; 随着振荡频率升高, 翼型俯仰和横摆振荡下的气动迟滞性均增强; 翼型俯仰振荡正行程的动态失速涡破裂有所延迟; 洞壁与模型端部交界处的强三维效应对翼型压力分布影响较大; 建立的横摆振荡试验技术可为风力机动态掠效应的研究提供技术支撑.      

UNCONDITIONAL STABLE SOLUTIONS OF THE EULER EQUATIONS FOR TWO－AND THREE－D WINGS IN ARBITRARY MOTION

ＵＮＣＯＮＤＩＴＩＯＮＡＬＳＴＡＢＬＥＳＯＬＵＴＩＯＮＳＯＦＴＨＥＥＵＬＥＲＥＱＵＡＴＩＯＮＳＦＯＲＴＷＯ－ＡＮＤＴＨＲＥＥ－ＤＷＩＮＧＳＩＮＡＲＢＩＴＲＡＲＹＭＯＴＩＯＮＧａｏＺｈｅｎｇｈｏｎｇ（高正红）（ＲｅｃｅｉｖｅｄＪａｎ．１２,１９９５,Ｃ...     

PrandtlPlane Joined Wing: Body freedom flutter,limit cycle oscillation and freeplay studies

Dynamic aeroelastic behavior of a joined-wing PrandtlPlane configuration is investigated herein. The baseline model is obtained from a configuration previously designed by partner universities through several multidisciplinary optimizations and ad hoc analyses, including detailed studies on the layout of control architecture. An equivalent structural model has then been adopted to qualitatively retain similar aeroelastic properties.Flutter and post-flutter regimes, including limit cycle oscillations (LCOs), are studied. A detailed analysis of the energy transfer between fluid and structure is carried out; the areas in which energy is extracted from the fluid are identified to gain insights on the mechanism leading to the aeroelastic instability. Starting from an existing design of control surfaces on the baseline configuration, freeplay is also considered and its effects on the aeroelastic stability properties of the joined-wing system are investigated for the first time.Both cantilever and free flying configurations are analyzed. Fuselage inertial effects are modeled and the aeroelastic properties are studied considering plunging and pitching rigid body modes. For this configuration a positive interaction between elastic and rigid body modes yields a flutter-free design (within the range of considered airspeeds).To understand the sensitivity of the system and gain insight, fuselage mass and moment of inertia are selectively varied. For a fixed pitching moment of inertia, larger fuselage mass favors body freedom flutter. When the moment of inertia is varied, a change of critical properties is observed. For smaller values the pitching mode becomes unstable, and coalescence is observed between pitching and the first elastic mode. Increasing pitching inertia, the above criticality is postponed; meanwhile, the second elastic mode becomes unstable at progressively lower speeds. For larger inertial values “cantilever” flutter properties, having coalescence of first and second elastic modes, are recovered.     

非定常俯抑振荡下的横向喷流数值模拟

采用高精度格式数值求解RANS方程,研究了定常状态下横向喷流流场,压力分布计算结果与实验结果基本吻合,并捕捉到喷流干扰流场中多种流动结构.在非定常计算过程中,飞行器的振动引起了法向力和俯仰力矩系数的相位滞后,推力放大因子随俯仰角周期变化.飞行器振动过程中,喷流流场的动态气动特性与稳态喷流有明显的区别,因此在利用横向喷流对飞行器进行姿态控制时,应该考虑由于飞行器姿态的变化对横向喷流所产生的非定常影响问题.     

Fourier time spectral method for subsonic and transonic flows

The time accuracy of the exponentially accu-rate Fourier time spectral method (TSM) is examined and compared with a conventional 2nd-order backward differ-ence formula (BDF) method for periodic unsteady flows. In particular, detailed error analysis based on numerical com-putations is performed on the accuracy of resolving the local pressure coefficient and global integrated force coefficients for smooth subsonic and non-smooth transonic flows with moving shock waves on a pitching airfoil. For smooth sub-sonic flows, the Fourier TSM method offers a significant accuracy advantage over the BDF method for the predic-tion of both the local pressure coefficient and integrated force coefficients. For transonic flows where the motion of the discontinuous shock wave contributes significant higher-order harmonic contents to the local pressure fluctuations, a sufficient number of modes must be included before the Fourier TSM provides an advantage over the BDF method. The Fourier TSM, however, still offers better accuracy than the BDF method for integrated force coefficients even for transonic flows. A problem of non-symmetric solutions for symmetric periodic flows due to the use of odd numbers of intervals is uncovered and analyzed. A frequency-searching method is proposed for problems where the frequency is not known a priori. The method is tested on the vortex shedding problem of the flow over a circular cylinder.