The flow around thick airfoils at low Reynolds numbers and the possibility of using energy methods to control flowfield, increase lift and reduce drag are considered. For calculations, a series of propulsion airfoils is used, obtained by a combination of solving the inverse problem of aerodynamics and stochastic methods of global optimization. The calculations take into account the transition from the laminar flow regime to the turbulent one, which has a significant effect on the lift and airfoil drag. The applicability of various turbulence models for airfoil modelling at low Reynolds numbers is discussed. The results of numerical simulation make it possible to determine the dependence of the lift coefficient on the angle of attack and the optimal flow rate of air taken from the airfoil surface for each angle of attack. The results of numerical simulation are compared with the data of a physical experiment, and conclusions are drawn about the accuracy of various turbulence models.