Kinematic Optimization of a Flexible Wing Undergoing Flapping and Pitching

In recent years, there has been widespread interest in the design of microair vehicles (MAVs) for flapping flight with high-aspect ratio wings due to their high efficiency and energy savings. However, the flexibility of a flapping wing causes the aeroelastic effect, which remains a subject of investigation. Generally, existing research simulates active bending and twisting of flexible wings under the assumption of neglecting flapping inertia. In this research, the kinematic optimization of a bionic wing with passive deformation in forward flight while undergoing flapping and pitching is considered. To this end, a computational aeroelasticity framework, which includes the three-dimensional unsteady vortex lattice method (UVLM) and the Newmark-β method, is constructed for flapping flight. Under the assumption of linear elastic deformation, this tool is capable of simulating attached flows over a thin wing and capturing unsteady effects of wakes. A bionic numerical wing with an aspect ratio of 6.5, chord Reynolds number of 1.9 × 105, and reduced frequency less than 0.1 is investigated in kinematic optimization. The computational aeroelasticity framework is combined with a global optimization algorithm to identify the optimal kinematics that maximize the propulsive efficiency under the minimum average lift constraint. Two types of numerical wings, rigid wing and flexible wing, are considered here to compare the influence of deformation on the aerodynamics of the flapping wing. The results show that the aeroelastic effect, which increases the flapping amplitude, yields a significant improvement in terms of propulsive efficiency. In addition, the optimization algorithm maximizes the thrust efficiency while satisfying the required lift. Moreover, the optimal kinematics of both the rigid wing and the flexible wing reach the maximum flapping angle, which indicates that a larger range of motions is needed for optimal kinetics when loosening the boundary conditions.

Efficient aerodynamic derivative calculation in three-dimensional transonic flow

ABSTRACTOne key task in computational aeroelasticity is to calculate frequency response functions of aerodynamic coefficients due to structural excitation or external disturbance. Computational fluid dynamics methods are applied for this task at edge-of-envelope flow conditions. Assuming a dynamically linear response around a non-linear steady state, two computationally efficient approaches in time and frequency domain are discussed. A non-periodic, time-domain function can be used, on the one hand, to excite a broad frequency range simultaneously giving the frequency response function in a single non-linear, time-marching simulation. The frequency-domain approach, on the other hand, solves a large but sparse linear system of equations, resulting from the linearisation about the non-linear steady state for each frequency of interest successively. Results are presented for a NACA 0010 aerofoil and a generic civil aircraft configuration in very challenging transonic flow conditions with strong shock-wave/boundary-layer interaction in the pre-buffet regime. Computational cost savings of up to 1 order of magnitude are observed in the time domain for the all-frequencies-at-once approach compared with single-frequency simulations, while an additional order of magnitude is obtained for the frequency-domain method. The paper demonstrates the readiness of computational aeroelasticity tools at edge-of-envelope flow conditions.