A Modified Van Der Pol Oscillator Model for the Unsteady Lift Produced by a Flapping Flat Plate for Different Positions of the Rotation Axis

To understand the nonlinear interaction between unsteady aerodynamic forces and the kinematics of structures, we theoretically and numerically investigated the characteristics of lift coefficients produced by a flapping thin flat plate controlled by the rotation axis position. The flat plate was placed in a 2-D incompressible flow at a very low Reynolds number (Re = 300). We showed that the behavior of the unsteady aerodynamic forces suggests the existence of a limit cycle. In this context, we developed a Reduced Order Model (ROM) by resolving the modified van der Pol oscillator using the Taylor development method and computational fluid dynamics (CFD) solutions. A numerical solution was obtained by integrating the differential equation of the modified van der Pol oscillator using the fourth-order Runge–Kutta method (RK4). The model was validated by comparing this solution with the reformulated equation of the added mass lift coefficient. Using CFD and ROM solutions, we analyzed the dependency of the unsteady lift coefficient generation on the kinematics of the flapping flat plate. We showed that the evolution of the lift coefficient is influenced by the importance of the rotation motion of the Leading Edge (LE) or Trailing Edge (TE), according to the position of the rotation axis. Indeed, when the rotation axis is moved towards the LE, the maximum and the minimum values of the lift coefficient are proportional to the downward and upward motions respectively of the TE and the rotation axis. However, when the rotation axis is moved towards the TE, the maximum and the minimum values of the lift coefficient are proportional to the downward and upward motions respectively of the LE and the rotation axis.