# Numerical Investigation of the Effect of Airfoil Thickness on Onset of Dynamic Stall

## TL;DR

This study uses large eddy simulations to examine how airfoil thickness influences the onset of dynamic stall, revealing that stall mechanisms vary with thickness and involve laminar separation bubbles and boundary layer separation.

## Contribution

It provides new insights into the effect of airfoil thickness on dynamic stall onset mechanisms using high-fidelity simulations, highlighting a transition from laminar to turbulent boundary layer separation.

## Key findings

- Thinner airfoils exhibit two-stage transition not seen in XFOIL.
- Dynamic stall onset is marked by bursting of the laminar separation bubble.
- Thicker airfoils show stall triggered by separated turbulent boundary layer.

## 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.