Published online by Cambridge University Press: 26 October 2017
We study the influence of shock and boundary layer interactions in transonic flutter of an aeroelastic system using a Reynolds-averaged Navier–Stokes (RANS) solver together with the Spalart–Allmaras turbulence model. We show that the transonic flutter boundary computed using a viscous flow solver can be divided into three distinct regimes: a low transonic Mach number range wherein viscosity mimics increasing airfoil thickness thereby mildly influencing the flutter boundary; an intermediate region of drastic change in the flutter boundary due to shock-induced separation; and a high transonic Mach number zone of no viscous effects when the shock moves close to the trailing edge. Inviscid transonic flutter simulations are a very good approximation of the aeroelastic system in predicting flutter in the first and third regions: that is when the shock is not strong enough to cause separation, and in regions where the shock-induced separated region is confined to a small region near the trailing edge of the airfoil. However, in the second interval of intermediate transonic Mach numbers, the power distribution on the airfoil surface is significantly influenced by shock-induced flow separation on the upper and lower surfaces leading to oscillations about a new equilibrium position. Though power contribution by viscous forces are three orders of magnitude less than the power due to pressure forces, these viscous effects manipulate the flow by influencing the strength and location of the shock such that the power contribution by pressure forces change significantly. Multiple flutter points that were part of the inviscid solution in this regime are now eliminated by viscous effects. Shock motion on the airfoil, shock reversal due to separation, and separation and reattachment of flow on the airfoil upper surface, also lead to multiple aerodynamic forcing frequencies. These flow features make the flutter boundary quantitatively sensitive to the turbulence model and numerical method adopted, but qualitatively they capture the essence of the physical phenomena.