Calculation and Selection of Airfoil for Flapping-Wing Aircraft Based on Integral Boundary Layer Equations

The flight of a migratory bird-like flapping-wing aircraft is characterized by a low Reynolds number and unsteadiness. The selection of airfoil profiles is critical to designing an efficient flapping-wing aircraft. To choose the suitable airfoil for various wing sections, it is necessary to calculate the aerodynamic forces of the unsteady two-dimensional airfoil with a Reynolds number in the range of 105. While accurate, calculating this by solving the Navier–Stokes equations is impractical for early design stages due to its high consumption of computing resources and time. The computational demands for extending it to 3D aerodynamic calculations are even more prohibitive. In this paper, a relatively simple method is proposed. The two-dimensional unsteady panel method is utilized to derive the inviscid flow field, the unsteady integral boundary layer method is utilized to solve the boundary layer viscous flow, and the eN transition model is adopted to predict the position of the transition. These models are coupled with the semi-inverse interaction method to solve the aerodynamics of the unsteady low-Reynolds-number two-dimensional airfoil. The unsteady aerodynamics of the symmetric and cambered airfoils at different wing sections are calculated respectively by the proposed method. Mechanism analysis of the calculation results is conducted, and a symmetrical airfoil or a slightly cambered airfoil is recommended for the wing tip, a moderately cambered airfoil is suggested for the outer-wing section, and a highly cambered airfoil is suggested for the inner-wing section.