A numerical model for two-dimensional flows around a pitching foil in a viscous flow is presented. The model is numerically solved using the immersed boundary method and used to investigate the flow patterns of the foil pitching sinusoidally over a range of frequencies and amplitudes. A transition from the Karman vortex streets to the reverse Karman vortex streets are found, as the amplitudes of pitching motions increase. In the transition, the vortex streets undergo symmetry-breaking to the central lines of vortex streets. Those observations are in agreement with the previous experiment (Phys. Rev. E. 77 016308 2008). Furthermore, we examine the wake of the foils pitching with different frequencies. The transition from the Karman vortex streets to the reverse Karman vortex streets is also observed. An explanation is presented to the mechanism of the transition.
INTORDUCTIONThe pitching of tails is a common mode for animals to generate thrust. As a tail pitches vertically, thrust can be horizontally generated, depending on either amplitudes or frequencies of pitching motions. A simple model to understand the mechanism of pitching-generated-thrust is the Karman vortex streets (KVS) and the reverse Karman vortex streets (RKVS). At a certain range of Reynolds numbers, a KVS can be observed behind a steady foil without thrust generation. As either amplitudes or frequencies increase of pitching motions, the transition from the KVS to the RKVS can be observed and thus a thrust is generated. The present paper is devoted to numerically investigate the transition from the KVS to the RKVS. [2] were the first to observe that a flapping wing could generate thrust. Karman and Burgers [3] offered the first theoretical explanation of drag or thrust production based on the wake vortices, where the wake of the flow past bluff bodies is modeled by an infinite row of alternating vortices commonly known as von Karman Vortex Streets. In following years, many researchers investigated the RKVS in different situations, such as the wakes of fish [4] in experiment [5] and numerical simulations [6,7]. More recently, Godoy-Diana et al. [7] experimentally investigated the transition from the KVS to the RKVS in the wake of pitching foil. We will simulate the transition process using the immersed boundary method and compare the results obtained with Godoy-Diana et al' s experimental results.
Knoller [1] and Betz
NUMERICAL METHODTwo-dimensional incompressible Navier-Stokes equations are used209where v is the kinematic viscosity. The geometry of the pitching foil is described in Figure 1, where the chord C = 23 mm and the semicircle diameter D = 5 mm. The foil pitches with respective to the center of the semicircle. The control parameters are the velocity U , the oscillation frequency f and the peak to peak amplitude A. The main nondimension parameters are the Reynolds number Re, the pitching amplitude A D and the Strouhal number St, defined as:
The pitching motion of foils is a sinusoidal function of time: A(t) = −A D sin(2π f t).Th...