In the present work, verifying the applicability and potentiality of the IMERSPEC methodology for numerical and computational modeling of two-dimensional flows over airfoils and vertical axis wind turbines is proposed. It is a high-order convergence methodology with low computational cost when compared to other high-order methods, resulting from the coupling of the Fourier pseudo-spectral method and the immersed boundary method. To validate the proposed methodology, flow simulations are carried out over an airfoil NACA 0012 for a Reynolds number equal to 1000. From the spatial discretization procedure, there is convergence and good agreement of the lift and drag coefficients and the Strouhal number in relation to reference works. The behavior of the flows over the airfoil, as a function of the angle of attack, is evaluated by pressure and vorticity fields. From the analyzed flows, it is observed that the formation of different wake modes, constituted by swirling structures that vary their characteristic sizes, is influenced by the angle of attack. A case study is proposed based on the analysis of the main fluid dynamic aspects of flows over wind turbines with a vertical axis of three blades for a Reynolds number equal to 100. For this, a mathematical model responsible for the imposition of the rotational movement on the blades is presented in the turbine. Performance parameters, such as the coefficient of tangential force and normal force, and the analysis of velocity fields on the simulated turbine were presented and compared with other numerical methods. The good level of convergence and the accuracy of the obtained results show the promising capacity of the IMERSPEC methodology in solving problems of this nature.