Numerical Investigation of a Dynamic Stall on a Single Rotating Blade

: Dynamic stall is a phenomenon on the retreating blade of a helicopter which can lead to excessive control loads. In order to understand dynamic stall and ﬁll the gap between the investigations on pitching wings and full helicopter rotor blades, a numerical investigation of a single rotating and pitching blade is carried out. The ﬂow phenomena thereupon including the Ω -shaped dynamic stall vortex, the interaction of the leading edge vortex with the tip vortex, and a newly noticed vortex structure originating inboard are examined; they show similarities to pitching wings, while also possessing their unique features of a rotating system. The leading edge/tip vortex interaction dominates the post-stall stage. A newly noticed swell structure is observed to have a great impact on the load in the post-stall stage. With such a high Reynolds number, the Coriolis force exerted on the leading edge vortex is negligible compared to the pressure force. The force history/vortex structure of the slice r / R = 0.898 is compared with a 2D pitching airfoil with the same harmonic pitch motion, and the current simulation shows the important role played by the swell structure in the recovery stage.


Introduction
Dynamic stall on the retreating blade of a helicopter is a factor that inhibits the progress of fast forward flight or manoeuvre flight; the subsequent aero-elastic response of the blade serves as a main source of vibrations and high structure loads. Since the discovery of dynamic stall, it has been a heated topic in the helicopter community. Understanding the phenomenon to implement effective flow control to ease the structural loads is the main purpose of the research.
Numerous experimental investigations [1][2][3][4][5] on oscillating airfoils; especially, the NACA 0012 airfoil was carried out in order to understand the mechanism of dynamic stall, as well as how the airfoil type, pitching amplitude, pitching rate, reduced frequency and the compressible effect impact aerodynamic coefficients. Ref. [6] used smoke as a flow indicator, showing the vortex shedding process over both small and large amplitude oscillations. Ref. [7] summarised early progress made in understanding dynamic stall, including the contribution of the dynamic stall vortex shedding to the surge of the normal force, the chord-wise force and the negative pitch moment; the role that reduced frequency or the pitching rate plays in the process; and the onset of vortex shedding as a result of the reverse flow upstream propagation from the trailing edge.
With the development of computational fluid dynamics in the late 1990s, many numerical simulations were carried out and compared with experimental data [8][9][10][11]. It was soon realized that 2D laminar simulations generally underestimate the unsteady load, while 2D simulations with turbulence models usually over-estimate the load. Experiments and numerical simulations on full rotor blades suggest a more complex flow structure, motivating research on vortex propagation [12] and the three-dimensional (3D) effect [13][14][15][16] on dynamic stall. Ref. [14,16] have done various experiments as well as numerical simulations as well as the validation of grid strategy and turbulence model on the hovering case of Caradonna-Tung blades are placed in Appendices A and B, respectively.

Numerical Method Flow Configuration and Numerical Setup
Currently, the chair of helicopter technology of the Department of Aerospace and Geodesy, TUM is preparing an experiment on measuring flow field and pressure on a single rotating blade. The simulation serves as a preliminary investigation of the phenomenon. According to the principle of simplicity, the blade is chosen to be, similar to the Caradonna-Tung blade, a rectangular blade with the NACA 0012 profile without a twist along the span. In order to mount pressure sensors thereon, the chord is selected to be c = 0.15 m in the first place; and in order to fit into the wind tunnel at Department of Mechanical Engineering, TUM, the span is selected as R = 0.8 m. The blade model used for the simulation is built by extruding the airfoil profile from r/R = 0.25 to r/R = 1; tips are rounded by revolving the half of the airfoil profile about its chord.
The advance ratio is set to be µ = 0.2 and rotation speed of the rotor be 275 rad/s or 2626.056 rpm; the rotor has 0 shaft angle and no flap, and the control is with a trim condition: thrust T = 900 N , C My = 0 and C Mx = 0, given by CAMRAD II [22] as: collective pitch θ c = 13.18 • , cyclic pitch θ s = −5.57 • and θ c = 7.74 • . We have to state that the pitch control is based on the lower-order aerodynamic model in CAMRAD II, which doesn't necessarily guarantee the desired operation condition in real case. This work is a preliminary investigation of the dynamic stall phenomena on the retreating blade prior to the experiment, hence the deviation is tolerant and the control yielded by CAMRAD II is accepted as the setting for the rotor. Mathematically speaking, the pitch follows: θ = 13.18 − 5.59 sin Ψ + 7.74 cos Ψ = 13.18 + 9.548 cos(Ψ + ϕ), ϕ = 35.84 • . (1) The flow condition for the rotor is summarised in Table 1. The computation is performed with the URANS code DLR-TAU [23]. Previous research has extensively validated that this code is capable of predicting dynamic stall on pitching wings [14,24], as well as for rotors [20]. For our case, a chimera technique is utilized, and the computation domain is split into three blocks: innermost blade block (Figure 1a) that resolves boundary layer on the blade and it does pure pitching; the outmost far-field block ( Figure 1c) with multi-resolution voxel grid refined with an inclined cylindrical ring volume that tracks the tip vortex; and the rotating block in-between (Figure 1b) with relatively uniform spaced cells, through which the tip vortex is convected to the far-field, and it does pure rotation. The flow information of each block is exchanged in the chimera region. The surface mesh on the blade is constructed to have a closer spatial resolution with simulations done by [21,24,25]. The spatial resolution of the blade block for the current simulation is summarized in Table 2. The blade block has a dimension of 3.2 c (3.2 times chord length) in radius and 6.76 c in length; it is constructed by extruding firstly the surface mesh normally 60 steps, and an O-type block is then assembled outside the boundary region; the overset mesh is constructed by extruding the outer surface mesh of the blade block algebraically 5 steps outward. The surface mesh has 11 million cells. The rotating block is a spherical volume with 5 c as radius, centering at the blade's rotation center; the voxel grid is generated between the outer surface of the previous overset mesh and the spherical surface domain. The average space in this block is around 5.9% c; the maximum space appears in the zipping region near the spherical surface 15% c. The chimera mesh is then constructed the same way by extruding equal spaced 6.67% c for five steps. The rotating block has 15.2 million cells. The farfield block extends upstream 15 R, downstream 35 R, upward 11.22 R, and downward 19.24 R, both sides horizontally 15 R. The voxel grid is utilized to generate finely spaced grid near the rotating blade and step by step coarsened till the farfield boundary. The possible tip vortex path is refined with an inclined cylindrical ring, whose center circle of upper surface coincides with rotor tip trajectory and has a width of 0.4 c, and a height of 1.5 R. The space within this tip vortex region is set to be 6.67% c. The farfield has altogether 11.5 million cells; the whole computation domain 37.7 million cells.  The unsteady problem is solved with a 2nd order central scheme inviscid flux discretization, 2nd order Roe upwind scheme for central convective turbulence flux, and implicit scheme for time stepping. The simulation is carried out first with 720 time steps per revolution ( ∆Ψ = 0.5 • ) for 2.5 revolutions, and then 1440 time steps/Rev (∆Ψ = 0.25 • ) for the subsequent revolutions. The inner iteration step is set to be at least 200 and maximum 400 to guarantee the density residual dropping two orders of magnitude. Furthermore, Cauchy convergence criteria are also implemented for saving computation time for each time step; if any of the force coefficients sequences C L (N), C D (N) and C my (N) in the last 30 inner iteration steps satisfies: and the inner iterations are considered converged for that time step. The residuals of the forces on blade, namely thrust F z and pitch moment M y , are plotted in Figure 2, where each circle represents the residual at an azimuth position after 400 inner iterations or reaches Cauchy convergence criteria. The ranges where no data points show up are where the residuals drop below 10 −6 . We perform simulations with both the Spalart-Allmaras model (SA model) and Menter k − ω SST turbulence models for the comparison of total forces, and the flow details with the SA model are researched thoroughly. The k − ω SST model was widely implemented on pitching wings [16,26] and blades in the rotor environment [21], and the SA model is reported to be inaccurate in predicting dynamic stall on pitching wings in some cases. However, the one-equation turbulence model (SA model) consumes less computation resources; the SA model shows a similar trend of the net force on the hub to the k − ω SST model; the averaged forces are nearly equivalent for the current case, as shown in Appendix A; the combination of such mesh strategy and SA turbulence model can capture the pressure distribution on the hovering rotor blade at the transonic Mach number and high Reynolds number, as shown in Appendix B. Considering all these aspects, the numerical result is still a worthwhile preliminary investigation into the dynamic stall on a pitching rotating blade.

Coordinate System
The Cartesian coordinate systems used in the analysis are illustrated in Figure 3, namely the inertial coordinate and body-fixed coordinate. The inertial coordinate is stationary, while the body-fixed coordinate rotates and pitches along with the blade. The azimuth angle Ψ is defined as 0 at the minimum x inertial , and it increases counterclockwise. The pitch angle is defined as positive when the blade's leading-edge is positioned to +z inertial . We only use azimuth angle Ψ and pitch angle θ to describe the position and orientation of the blade, respectively. Hence, If not specified, x, y, and z are referred to the body-fixed coordinate system.

Force and Moment
There is a marginal difference in both thrust and pitching moment coefficients in terms of the maximum values as shown in Figure 4. The stall event occurs later for the k − ω SST turbulence model, and the pitching-up moment in the post-stall phase is relatively smaller. This indicates differences of the detailed fluid structure in the stall regime; this will be studied in detail in the future but not within the scope of current study. Nevertheless, the overall trends of the curves are similar to each other, and we are using the numerical results of the SA turbulent model to illustrate the dynamic stall events on the rotating-pitching blade.   Figure 6 is the rotor map of both sectional normal force and moment coefficients. Force coefficients are non-dimensionalized by sound speed. The maximum CnM 2 appears in the region 0.6 < r/R < 0.85, 350 • < Ψ < 35 • , which is in the pitching-down phase. While the relative high value (orange area) begins at Ψ = 300 • , θ = 21.89 • ↑ (pitching up) and remains until Ψ = 60 • , θ = 12.21 • ↓ (pitching down). The negative normal force also exists near blade tip at azimuth angle 330 • < Ψ < 30 • , which is a result of the tip vortex; it appears at the blade root as well, which is the attributed reverse flow region and the strong root whirl flow, the latter a pure result of the abrupt cut of the blade model at the root. The minimum CmM 2 (pitch down moment) appears at r/R = 0.72, Ψ = 20 • , which is in the pitching-down phase well beyond the maximum pitch angle. The relatively small value (blue region) starts from Ψ = 300 • and ends near Ψ = 30 • , which is consistent with the change of normal force, and this indicates that the dynamic stall occurs in this region. The maximum CmM 2 (pitching-up moment) appears at blade tip, 90 • < Ψ < 110 • , where a normal shock exists on the suction surface, see Figure 7 . Cm = 0 is contoured with a bold line in the rotor map. The rugged line shows the recovery of Cm occurs inboard first while very late at radial locations 0.6 < r/R < 1.

Vortex Structure
Since dynamic stall events are closely associated with the vortex generating and shedding on the upper surface of a blade, we show the vortex structures on the rotating blade at different azimuth angles in Figure 8. The location of vortex cores over the blade is plotted in Figure 9, which is extracted based on the eigenmodes of the velocity vectors [27]. This will be compared with previous research to illustrate the difference in rotating environment. An which is similar to the observations for a pitching finite wing by [14], and also similar to what [15] called the Λ-type arch vortex. A similar observation is also reported by [20]. Unlike on a pitching finite wing, as soon as the Ω-shape vortex is shed, the leading edge vortex (LEV) comes into interaction with the blade tip vortex to a great extend, and even an inclined arch vortex appears at the tip, see (e) Ψ = 345 • , θ = 22.1 • ↓.
In order to see the interaction of the leading-edge vortex and the tip vortex, we take slices at several chord-wise locations near the tip region (0.94 < r/R < 1.125), and integrate ω x ,±ω y in this region where ω x > 0, which is the main contribution to the tip vortex in cases that are free from strong interaction. The circulation in the y-direction indicates the mixture of LEV in the tip vortex region. This result is plotted in Figure 10 with selected azimuth angles.
The consequent strong interaction of the LEV and tip vortex dominates the flow behavior of the post-stall phase. There are mainly four stages in this phase, symbolized by the shedding of LEV at the tip area. The 1st stage begins at Ψ = 300 • , when the LEV accumulates at the blade-tip's leading-edge, which is shown as the increasing integration of negative y circulation −Γ y ; the LEV starts shedding at Ψ = 317 • as the peak of −Γ y moves rearward, and at the place where this peak locates, the circulation in the x-direction Γ x in the tip region is also relatively larger. As the separated DSV moves downstream and outboard, the tip vortex at x/c = 0.957 is elevated and distorted. Alongside the rearward movement of −Γ y , +Γ y grows in the vicinity of the tip surface. This explains why the negative C n exists at the tip region on the rotor map of C n M 2 . The second stage begins at Ψ = 0 • when the previous main DSV is transported off of the trailing edge; another LEV is accumulated at the leading edge. In this stage, the tip vortex at x/c = 0.952 is no longer obvious, but the net circulation in the x-direction increases. This means, in the tip region that there exists a more diffused vorticity field rather than the conventional concentrated structure. As the LEV grows over the chord, the tip vortex is not apparent; nevertheless, it increases continuously in this area. From this time on, the LEV sheds and entrains the tip vortex, forming a large separation region over the tip area. Finally, Ψ = 35.96 • marks the end of previous shedding of LEV and the beginning of the next shedding of LEV at the tip area. A reversely rotating tip vortex is generated beneath the mixture layer of the separated leading-edge vortex and tip vortex at x/c = 0.952. In the near wake, the tip vortex locates more outboard compared to the 1st stage, and a pair of counter rotating vortices is obvious. At azimuth Ψ = 53.96 • , the interaction of LEV and tip vortex begins to disentangle at x/c = 0.418. Another growing and shedding of LEV at the tip area is observed, but in this stage there is no strong interaction with the tip vortex. It enters the recovery stage. The interaction of the leading edge vortex with the tip vortex represents three important features: • The leading-edge vortex is not "pinned" by the tip vortex, which is contrary to the observations for a pitching wing; • There is a strong correlation of circulation in the x-direction and the chord-wise location of DSV. This can be explained as a result of lower pressure created by the DSV, which induced higher velocity around the tip. • The concentrated tip vortex is entrained into LEV during the pitching-down phase, and a pair of counter-rotating vortices is observed in the near wake of the blade.
Ref. [28] measured the tip wake in the near field of an oscillating wing; he discovered the hysteresis behaviour of the tip vortex, and a more diffused tip vortex during a down-stroke compared with the up-stroke. Current numerical simulation complies with the observation of pitching wing. Another significant feature shown in Figure 8 is the conical structured LEV on the rotating blade, which is very similar to the research on a rotating wing by [29]. This conical shaped structure is treated as a more stable LEV, "pinned" at the leading edge [30], if we look at slices taken along the chord-wise direction. Ref. [30] explains that the stability of LEV on rotating wing is due to the span-wise flow, which transports vorticity to the tip and inhibits the growth of LEV, hence the critical circulation beyond which the LEV will detach from the leading edge is not reached. In addition, this can explain the phenomenon observed by [17], who concluded that vortex evolution varies with radial positions on the blade, namely, at some locations, LEV sheds into dynamic stall vortex (DSV), while, in other places, only a secondary vortex is generated. When the phenomenon inspected from a three-dimensional perspective, this conclusion is just an incomplete description of the Ω-shaped vortex outboard and the conical vortex structure inboard. That is to say, the detached DSV is rather a part of the arch vortex, which attaches to the blade surface with two "legs"; and this is due to the stability of the conical LEV structure in which no DSV is shed and observed inboard.
The special phenomenon of the rotating and pitching blade is the vortex generated inboard which can be seen in Figure 8b 20.92 • ↓) as a swelling part within the conical LEV region inboard the blade, and it appears on different radial locations and grows in size while moving outboard. This swell structure creates a relative higher pressure on the blade, as is shown in Figure 11 where a shallow grey area at r/R = 0.4875 is encompassed by a darker area. This is a result of the alleviated LEV, when compared with the one at r/R = 0.625, where the LEV stays close to the upper surface of the blade. From the x plane slice, one can easily identify the detachment of the vortex structure. This structure corresponds to what [17] have observed on their retreating blade that inboard LEV didn't generate a DSV but a secondary vortex. As the swell structure moves outboard, it carries the vorticity that arched away from the surface and hence at that radial location no shedding of the LEV occurs. The span-wise flow also contributes to the outboard moving of the swell structure. This swell structure seems to be a result of the Coriolis force on the rotating blade at first thought, but we carefully examined the pressure force and the Coriolis force of the conical LEV slices at several spanwise locations and found out that the Coriolis force is too small to affect the vortex structure. The swell structure shown in Figure 8 exists at Ψ = 270 • near the root corner, prior to the formation of the Ω-type vortex. Indeed, the structure appears even earlier in this simulation at Ψ = 245 • , but it doesn't have any obvious effect on the pressure distribution on the blade's surface. The growth, however, does have an impact on the pressure distribution later on as shown in Figure 11. The reason why this structure appears first inboard is still unclear. However, it cannot be explained by a relative stronger Coriolis force towards blade root at low Reynolds number: a greater magnitude inertial Coriolis force that pulls the LEV core away from the leading edge exists when the a rotating wing is positioned closer to the rotation axis [19]. We evaluated forces in the x-direction on the slices of the vortex structure, as shown in Figure 12a. Following the divergence theorem, we evaluated the pressure force on the slices of the vortex core: The Coriolis force exerted on the vortex core can be evaluated as ∆F Corx = ρ(ΩV y )dA. The ratio of the Coriolis force to the pressure force is plotted in Figure 12b. At blade root, this value is as small as 10 −3 , and it decreases further as the non-dimensional spanwise location r/R increases. Although V y increases with r/R, this linear increase is not comparable to the quadratical increase of the pressure force, resulting in a continuously decreasing ratio as r/R increases. It seems that, at such a high Reynolds number, the Coriolis force has little effect on the vortex core structure. Possible explanations for the emergence of the swell structure can be the effect of the jet-like radial flow itself or the upward flow at the blade root. Further research with a rotor hub or a winglet at root can help to clarify the mechanism of the vortex structure's onset. Since the position of the swell structure influences the sectional pitching moment, we thereby extracted the vortex cores and plot its x position against azimuth angle Ψ in Figure 9. The lighter color represents an aft position of the vortex core, which results in a negative pitching moment component. In Figure 5 at a radial location r/R = 0.607, Ψ = 340 • , there is a small kink on both CnM2 and CmM2 curves; In Figure 9, we see that, at Ψ = −21 • (or 339 • ), the vortex core at r/R = 0.62 is located well beyond the line of the nearby LEV cores. The lighter color formed an edge in the contour map that represents the trajectory of the swell structure as it moves afterward as well as outward. We define this ridge as the position of the swell structure (r/R) swell when (1) ∂(x/c) VortexCores /∂(r/R) = 0, (2) ∂ 2 (x/c) VortexCores /∂ 2 (r/R) < 0, in radial locations (3) (r/R) < 0.8, where (x/c) VortexCores is the scatter plot of the normal projection of the vortex cores on the blade in Figure 9. This value (r/R) swell can be estimated as a quasi-uniformacceleration movement on the blade: (r/R) swell = 1 2 a swell ∆Ψ 2 + ηµ cos(Ψ)∆Ψ + r 0 .
The first term is the centrifugal effect, the second term is the contribution of the yawing effect, with η representing a lower flow speed in vicinity of the blade's upper surface, and the third term is the initial position of this vortex structure on the blade. In our case, the white dotted line in Figure 9 gives an approximation of the trajectory using the expression mentioned above for the movement, and a swell = 0.25, η = 0.8, r 0 = 0.3426.

Separation Points
The feature of dynamic stall is mainly the separation of the boundary layer, both at the trailing edge and the leading edge. Understanding, as well as modelling of dynamic stall need these separation points on the blade. Many methods, both Eulerian and Lagrangian, can be used to detect the flow separation in the unsteady flow. We present here pure Eulerian methods, the skin friction and shape factor criteria. The former theorem requiring wall shear stress is based on steady, laminar flows, but was recently extended by [31] to unsteady, turbulent and compressible flows; and the latter theorem requiring the near wall flow field is based on the boundary theory, and the separation criteria is proposed by [32] as H sep = 2.76 ± 0.23 for two-dimensional turbulent flows. We follow the rigorous skin friction criteria to extract the instantaneous separation points on the blade and then present the shape factor criteria resultant separation points to show how good the velocity field can be used to evaluate the separation points.
The skin friction lines are plotted in Figure 13 at selected azimuth angles. According to [33], the existence of a skin-friction line on a three-dimensional surface on which other lines converge is a necessary condition for the flow separation; and skin-friction line emerging from a saddle point indicates a global separation; otherwise, it is only a local separation. We don't distinguish here global or local separation, since we know that the leading-edge vortex that attached to the blade can also contribute to such a local separation, which is our interest as well. For the a blade section at a radial position, the aforementioned criteria can be expressed mathematically: Similarly, if c f x = 0 and ∂c f x /∂x > 0, this is the point of attachment. We extract the sectional separation points and attachment points for the blade in a rotor map in Figure 14. Note that the stagnation point is also an attachment point, hence we have excluded this point in our algorithm. The leading-edge separation of boundary layer starts outboard and inboard both at Ψ = 230 • , while, at mid-span r/R = 0.6, this starts at Ψ = 250 • . Together with the attachment points, we see that, between azimuth angle Ψ = 240 • and Ψ = 270 • , the flow re-attaches immediately beyond the separation points, which indicates the existence of attached leading edge vortex. When the blade enters the 3rd quadrant, the attachment points begin to move downstream, but the starting azimuth is highly dependent on the radial location. For example, at r/R = 0.89, the vortex begins to grow leeward very early while, at r/R = 0.5, the vortex remains at the leading edge for a long time. Note that the growth of LEV starts at both ends of the blade and propagates toward mid-span. A helicopter blade may not show the same character, for example simulation by [21], shows a rather inboard to outboard propagation of the growth of LEV, in which 7A rotor has a linear twist −8.3deg/R meaning larger pitch angle inboard. In our case, the pitch angle is the same for different radial stations, and we observed a double directional propagation of the growth of LEV. Ref. [24] mentioned that, in the vicinity of wing tip, the flow is attached due to tip vortex induced flow as the wing pitches. In our case, this holds for the azimuth angles before Ψ = 330 • . We have already seen in Figure 10 that there exists a strong correlation of tip vortex strength Γ x and the leading edge vortex strength Γ y− , and it is indeed this interaction of LEV with tip vortex makes the LEV no longer staying attached at tip. At Ψ = 90 • , we see a light-coloured region out board, where the separation occurs between x/c = 0.121 and x/c = 0.223, far away from the other leading edge separation positions. In Figure 15, we show the Mach contour, pressure coefficient c p and skin-friction coefficient in the x-direction c f x at radial location r/R = 0.898 at this azimuth angle. On the upper surface of the blade section, c f x drops from positive to negative and crosses 0 at x/c = 0.12, where the −c p shows a sharp decrease. This separation is obviously a product of shock wave. Based on the discussion above, we can categorise the rotor map into four regions, namely fully attached region (F.A.), leading-edge separation region (L.S.), fully separated region (F.S.) and a shock induced separation region (SI). These regions are plotted in Figure 16. Another separation criterion that can be used (see [21]) to detect separation in turbulence flows is based on the shape factor H i , namely the ratio of the momentum thickness over displacement thickness of the boundary layer for a two-dimensional flow, as illustrated in Figure 17: Ref. [32] stated that H sep = 2.76 ± 0.23 is characteristic of the boundary separation, yet they did not give a sufficient criteria. In order to imply the concept on our rotating blade, we transformed all the velocity to the blade-section-fixed coordinate and dealt with the velocity parallel to the airfoil at different radial sections, as shown in Figure 17, along n direction. The flow along radial direction is not considered when analysing the shape factor since observing from the blade, and the main component of the flow is perpendicular to the leading edge of the blade. The outer edge of the boundary layer is treated as the location where the velocity reached its maximum along the n ⊥ direction. To determine how good this criterion can predict separation using only velocity field, we plot in Figure 18 the shape factor and the skin-friction in the x-direction at azimuth position Ψ = 32.9 • . The lower bound H sep = 2.53 shows quite a good agreement with the x skin-friction contour on the blade surface. However, Ref. [32] shows only the correlation of separation and shape factor H in turbulent boundary layers, the utilisation of such criteria to detect separation should be considerably careful.

Comparison of Stall Events on a Pitching Airfoil and a Section of the Rotating Blade
In order to understand the differences of dynamic stall between a rotating blade section and a pitching airfoil with the same harmonic pitching, the force coefficients and the Q contours of both cases are plotted in Figure 19. This spanwise location is exactly the center of the Ω-type dynamic stall vortex. The numerical investigation of a pitching airfoil satisfying equation (1) is explained in detail by [34]. The main differences between the pitching airfoil (2D case) and the section of the three-dimensional rotating blade (3DR case) include: (1) the magnitude of the force coefficients (both C n and C m ) differ significantly in the pitching-up phase; (2) the onset of the moment stall of the pitching airfoil is, in terms of the pitch angle, earlier than the blade section. However, the C n curves are quite close in the pitching-down phase for both cases and the normal force stall occurs almost simultaneously.
The maximum normal force coefficient C n of the airfoil is 2.21 while the value of the blade section is only 1.16. In the process when the leading-edge vortex grows, there is an obvious increase of the slope C nα , while this is not observed on the blade section. The normal force stalls for both cases are near 20 • , and, in the post-stall stage, normal force coefficients C n are relative close to each other. The extreme value of the moment coefficient C m of the 2D case is −0.52, while the counterpart on the rotating blade is only −0.145. The moment stall of the blade section lags behind that of the pitching airfoil, and is relatively milder. Comparing the vortical structures at different dynamic stall stages, there seems to be a delayed LEV growth from (a) to (c). When the LEV is quite obvious on the pitching airfoil, the LEV is unable to be seen on the blade section at time (a). Time (b) for both cases are the stages before the shedding of LEV. However, both Q contours are comparable to each other at time (c), with both of them representing the aft-moving dynamic stall vortex. We see here that there is no phase of growing LEV on the blade section, as seen in time (b) of the 2D case. It seems that, at the moment when the LEV appears on the blade section, it pinches off immediately and moves aft-ward, which is significantly different from what happens on the 2D pitching airfoil. This strange behaviour of the LEV on the blade section is a consequence the following factors: (1) the effect of the induced velocity field in the rotating environment, which yields a smaller effective angle of attack (AoA) α e f f , and hence results in the delayed effect on the moment stall; (2) the outward transportation of y vorticity ω y by the radial flow, which reduces the strength of the leading-edge vortex, and results in the lack of growth phase of the LEV on the blade section; (3) the Ω-type vortex structure, which arches itself away from the surface, moves aft-ward and allows the vorticity that is transported from inboard to form a new LEV, and hence results in a strong increase in C n and a mild decrease of C m after the detachment of the LEV. Time (d) shows the pitching-down phase, where the vortex structures are similar to each other, apart from an obvious vortex in the vicinity of the blade section's trailing edge. This is a consequence of the aforementioned factor (2): this vortex is transported from inboard. Moreover, the trailing-edge vortex on the blade section is more diffusive than the pitching airfoil. The vorticity ω − y of the counter-clockwise rotating vortex over the upper surface of the airfoil in Figure 19 is integrated and plotted in Figure 20. This value interprets the strength of the LEV or DSV when the leading-edge separation occurs. The two curves have a similar trend, in the upward pitching phase, the curve of the circulation is quite flat. Then, after a light decrease they surge up and in the downward pitching phase, they attenuate in oscillation. There are still three major differences: (1) The maximum circulation of the DSV in 2D case is larger; (2) The oscillation of the 2D case in the downward pitching scenario is milder and the DSV strength of the 2D case is obviously smaller; (3) There is a delay of 5.2 • for the 3D rotating blade case when comparing the point of which the leading edge vortex begins to grow. The first point can be explained as the lack of a third dimension for the vortex to grow, resulting in a concentration of vorticity. The second point can be explained as the continuing transportation of the vorticity from inboard, where the swell structure grows with the radial flow on the blade, as shown in Figure 9. At azimuth angle Ψ = 90 • , the circulation of the 3D slice is clearly larger, which is also the effect of the continuous transportation of vorticity from inboard. Figure 16 shows that, at this azimuth angle, the leading-edge separation still exists inboard. The third point can be explained as the combined effect of the induced velocity field V i in the rotating environment and the outward transportation of vorticity by the radial flow. The induced velocity field means that the effective angle of attack on a rotor blade is different from that on a pitching airfoil: α e f f = θ(Ψ) − V i (r, Ψ)/U(r, Ψ); Additionally, the outward transportation of vorticity means that, unless an increase of Γ net is present, the circulation over the blade section will not increase. Here, "0.898−" is the side of the blade section towards root, and "0.898+" is the side of the blade section towards tip. Blade element momentum theory (BEMT) is the lower order model that is widely used to predict the loads on the rotor disk. The in-flow sub-model is used for determining the effective angle of attack (AoA) for the blade elements. If we take the Drees linear inflow model [35] and implement the momentum analysis on the rotor for our case, we can get the modelled induced AoA as shown in Figure 21a,b. At the radial location r = 0.898, the induced AoA α i has a maximum value of 2.6 • , which is smaller than 5.2 • . Hence, the evaluation of transported circulation is another modelling factor to relate the stall event on the blade section to that on the pitching airfoil.

Conclusions
A numerical investigation with URANS using SA and k − ω SST turbulence models was carried out on a single rotating flat rectangular blade with a NACA 0012 airfoil, and the flow feature of the SA result was researched in detail. The collective controls were chosen in order to create dynamic stall on the retreating blade, and the cyclic control was acquired with CAMRADII with the minimum-lateral-moment C My trim condition. Several features not yet reported in previous research were found.

1.
A strong interaction of the DSV and the tip vortex was observed and described in detail for this case. Contrary to observations of a pitching wing, the LEV on a rotating pitching blade is not pinned by the presence of the tip vortex, but rather the shedded DSV drastically interacts with the tip vortex, yielding a pair of counter rotating vortices in the near wake; the slices also indicate a diffused tip vortex in the near wake, whether this diffusion is caused by the pure pitch-down effect as observed for the pitching wing or if it is further diffused by the interaction of the tip vortex is not yet clear.

2.
Despite an Ω-shaped vortex outboard, a swell vortex structure is observed to generate mostly inboard and moves outboard while gaining vorticity and size, which plays an important role after the shedding of the first DSV.

3.
The presence of the swell structure is not an outcome of the Ro effect inboard, contrary to the respective low Reynolds number counterparts. The Coriolis force is three orders of magnitude smaller than the pressure force, and hence unable to be the main cause. The mechanism of the onset of the swell structure is not yet clear.

4.
By tracing the vortex cores over the suction side of the blade, we proposed to use a quasi-uniform-acceleration movement to describe the position of the swell structure, and the parameters were given for this special case.

5.
By comparing the force coefficients on an outboard radial location with pitching airfoil, we show that a large difference exists between the two cases. The induced velocity field of the rotating environment, the radial flow due to the centrifugal effect, and the 3D effect created Ω-shape vortex altogether caused a simultaneous normal force stall and a delayed and mild moment stall on the blade section. The 2D simulation showed a stronger dynamic stall vortex that magnified both the normal force and the moment. 6.
By integrating the counter-clockwise rotating vorticity above the upper surface of the pitching airfoil and blade section, we showed the strong circulation in the recovery stage of the latter case, which is a consequence of the outward transportation of vorticity by the radial flow. The aforementioned swell structure also plays an important role when it crosses the blade section. The delayed increase of the circulation is assumed to be a result of the induced velocity field in the rotating environment and the transportation of vorticity due to the radial flow.
Based on the conclusions above, a more complex and unique vortex system exists on a pitching rotating blade. We thus suggest that, in order to model dynamic stall more precisely force modelling over a rotating blade, a better 3D vortex model is needed in addition to the blade-element-momentum theorem. Adjusting modelling parameters according to 2D pitching airfoil numerical simulations or experiments may not improve the accuracy of the model in predicting loads of a rotating blade.

Appendix A. Grid Convergence Study
In order to show the convergence of current grid, we created a coarse grid and medium grid with the same strategy as described in Section 2, and we list the spacing of each block for these grids in Table A1. h is the average of the edge length among the cells, and V is the average of the cell volumes. During coarsening the original fine grid, we keep the height of the first layer off the blade surface having the same value which satisfies y+ = 1. In addition, the number of connectors in three directions are all reduced according the level of the coarse grid. We present here only simulations with the Spalart-Allmaras turbulence model.

Coarse Grid Medium Grid Fine Grid
Pitch Block h = 5.74%c; Far-field Block h = 1.32c ; We choose the average of cell edge lengths in pitch block as the spacing indicator, and following [36], we can obtain the exact value of a partial differential equation as: where c is a constant and h is some measure of grid spacing, p is the order of convergence. We plot the grid convergence curve in Figure A1. The fitted curve of Equation (A1) is f = f h=0 − 1.251e − 5 · h 1.566 , which means the grid has an order of convergence of 1.566. We have used a 2nd order scheme for both mean-flow flux and turbulence flux with a Spalart-Allmaras model; theoretically, this value should be 2. Based on Equation (A1), the exact value of C T is 0.007942, and the current grid has an error of 0.4%. The superposition of C T and C My in different revolutions is plotted in Figure A2, and we can see that, after the 3rd revolution, the forces have already converged to a good extent. We have compared the average C T and C My of the SA model and the Menter k − ω SST model, and found that C T is exactly the same, while the averaged C My for SA model is −0.001328, and, for k − ω model, it is −0.001676. Figure A1. Grid convergence study with coarse, medium, and fine grids.

Appendix B. Validation of grid strategy and turbulence model on the hovering case
The Caradonna-Tung rotor model is adopted, and the surface mesh follows the same space criteria as described in Section 2. The only difference of the blade block is the root part, where unstructured pyramids' prisms and tetrahedra are used, as shown in Figure A3a. The farfield block consists of a region that is uniformly spaced as the rotating block in current simulation and the voxel grid growing from small to large scale on a farfield boundary. The region right below the rotor is refined with a point source in order to catch tip vortex trajectories, as shown in Figure A3b. The flow condition is summarised in Table A2, and the comparison of grid for the simulation and validation case are summarised in Table A3. The numerical method for the hovering case is the same as the current simulation, except that the rigid motion herein includes only rotation of the whole block. The Spalart-Allmaras turbulence model is adopted. Figure A4 shows the thrust coefficient along the simulation revolutions, and Figure A5 shows the comparison of the pressure coefficients on r/R = 0.5, r/R = 0.8, r/R = 0.96 span-wise locations. Figure A6 shows the tip vortex trajectory using Q-criterion shaded with vorticity magnitude.     The thrust coefficient of the experiment is C T = 0.00473, and the validation case gives a value of 0.00566, which is 19.6% larger than experiment data, which is consistent with the pressure coefficient, where, at r/R = 0.8, the maximum Cp has an error of 11.55%. Comparing the simulation result of the hover case with the Experiment data from [37], we conclude that the SA turbulence model is capable of qualitatively catching the main characteristic of the pressure distribution on the blade in rotation environment, and the voxel grid with point source as refinement is able to keep track of the tip vortex outside the fine grid region. The same mesh strategy and turbulence utilised for the current simulation is thus considered to give a close estimation of the flow phenomena on the pitching rotating blade.