A modified dual time integration technique for aerodynamic quasi-static and dynamic stall hysteresis

Simulation of the aerodynamic stall phenomenon in both quasi-static and dynamic conditions requires expensive computational resources. The computations become even more costly for static stall hysteresis using an unsteady solver with very slow variation of angle of attack at low reduced frequencies. In an explicit time-marching solver that satisfies the low Courant number condition, that is, CFL≤1
, the computational cost for such simulations becomes prohibitive, especially at higher Reynolds numbers due to the presence of thin-stretched cells with large aspect ratio in the boundary layer. In this paper, a segregated solver method such as the Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) is modified as a dual pseudo-time marching method so that the unsteady problem at each time step is reformulated as a steady state problem. The resulting system of equations in the discretized finite volume formulation is then reduced to zero or near-zero residuals using available convergence acceleration methods such as local time stepping, multi-grid acceleration and residual smoothing. The performance and accuracy of the implemented algorithm was tested for three different airfoils at low to moderate Reynolds numbers in both incompressible and compressible flow conditions covering both attached and separated flow regimes. The results obtained are in close agreement with the published experimental and computational results for both quasi-static and dynamic stall and have demonstrated significant savings in computational time.