Separation bubble model for low Reynolds number airfoil applications

Midchord laminar separation bubbles, which act as a transition mechanism on low Reynolds number airfoils, make a contribution to wing section profile drag that becomes increasingly important at low Reynolds number. A model for the analysis of the boundary layer through the bubble is needed. The model developed here, which is based on Horton's method, provides a simple computationally efficient analysis that matches the integral boundary-layer analysis methods used on most existing boundary-layer codes. The bubble calculation is initiated by the detection of laminar separation. Transition location and boundary-layer growth in the laminar region are determined using Van Ingen's shortcut e" method and Schmidt's correlations, respectively. Following Horton, the turbulent region is calculated using an iterative scheme that also functions as a bursting criterion, but the original linear velocity distribution has been replaced by Wortmann's concave velocity distribution. Both computation efficiency and prediction accuracy were improved after this change. Testing against experimental data showed that the bubble model greatly improved the drag prediction accuracy of the analysis in the Reynolds number range from 0.2 to 1.5 x 10, especially in cases when the midchord bubble was dominant. The validity of the bubble model was further confirmed by accurate prediction of bubble size and reattachment velocity gradient.