Abstract:Effect of airfoil thickness on onset of dynamic stall is investigated using large eddy simulations at chord-based Reynolds number of 200,000. Four symmetric NACA airfoils of thickness-to-chord ratios of 9%, 12%, 15%, and 18% are studied. The 3-D Navier Stokes solver, FDL3DI is used with a sixth-order compact finite difference scheme for spatial discretization, second-order implicit time integration, and discriminating filters to remove unresolved wavenumbers. A constant-rate pitch-up maneuver is studied with the pitching axis located at the airfoil quarter chord point. Simulations are performed in two steps. In the first step, the airfoil is kept static at a prescribed angle of attack ($=4^\circ$). In the second step, a ramp function is used to smoothly increase the pitch rate from zero to the selected value and then the pitch rate is held constant until the angle of attack goes past the lift stall point. Comparisons against XFOIL for the static simulations show good agreement in predicting the transition location. FDL3DI predicts two-stage transition for thin airfoils (9% and 12%), which is not observed in the XFOIL results. The dynamic simulations show that the onset of dynamic stall is marked by the bursting of the laminar separation bubble (LSB) in all cases. However, for the thickest airfoil tested, the reverse flow region spreads over most of the airfoil and reaches the LSB location immediately before the LSB bursts and dynamic stall begins, suggesting that stall could be triggered by the separated turbulent boundary layer. The results suggest that the boundary between different classifications of dynamic stall, particularly leading edge stall versus trailing edge stall are blurred. The dynamic stall onset mechanism changes gradually from one to the other with a gradual change in some parameters, in this case, airfoil thickness.
What problem does this paper attempt to address?
This paper attempts to study the effect of airfoil thickness on the onset of dynamic stall. Specifically, the authors used the Large Eddy Simulations (LES) method to study the behavior of NACA symmetric airfoils with four different thickness ratios (9%, 12%, 15%, and 18%) at the onset of dynamic stall under the condition of a chord - based Reynolds number of 200,000.
### Research Background
Dynamic stall is a non - linear fluid mechanics phenomenon, commonly seen in rapidly maneuvering airplanes, helicopter rotors, and wind turbines. This phenomenon is characterized by a significant increase in lift, drag, and pitching moment after the angle of attack exceeds the static stall value, far exceeding the values at static stall. Dynamic stall can be divided into two categories: mild stall and deep stall, mainly depending on the extent to which the angle of attack exceeds the static stall value.
### Research Methods
1. **Numerical Simulation Tools**: The widely - verified compressible Navier - Stokes solver FDL3DI was used for fluid flow simulation. This solver can handle multi - block overlapping grids and has a high - order interpolation method, which is suitable for high - precision simulation of complex geometries.
2. **Airfoil Selection**: Four NACA symmetric airfoils with different thickness ratios (9%, 12%, 15%, and 18%) were selected.
3. **Simulation Steps**:
- **Static Simulation**: First, the airfoil was fixed at an angle of attack of 4° to obtain a statistically stable solution.
- **Dynamic Simulation**: Then, the pitching rate of the airfoil was increased from zero to a selected value through a smooth ramp function and kept constant until the angle of attack exceeded the lift stall point.
### Main Findings
1. **Static Simulation Results**: Comparison with XFOIL showed that FDL3DI exhibited good consistency in predicting the transition position. For thinner airfoils (9% and 12%), FDL3DI predicted a two - stage transition, while this phenomenon was not observed in XFOIL.
2. **Dynamic Simulation Results**: The onset of dynamic stall is marked by the breakdown of the laminar separation bubble (LSB). However, for the thickest airfoil (18%), the reverse - flow region rapidly expanded to most of the airfoil and reached the LSB position immediately before the LSB breakdown, indicating that the stall may be triggered by the separated turbulent boundary layer.
3. **Change in Stall Mechanism**: As the airfoil thickness gradually changes, the onset mechanism of dynamic stall also gradually changes from one type to another, especially the boundary between leading - edge stall and trailing - edge stall becomes blurred.
### Conclusion
The airfoil thickness has a significant impact on the onset mechanism of dynamic stall. As the airfoil thickness increases, the stall mechanism gradually changes from the breakdown of the laminar separation bubble to the stall triggered by the separated turbulent boundary layer. These findings contribute to a better understanding of the physical mechanism of dynamic stall and provide theoretical support for the design of more efficient aircraft.