Computational Fluid Dynamics (CFD) is a promising tool for the analysis and optimization of wind turbine positioning inside wind parks (also known as wind farms) in order to maximize power production. In this paper 3-D, time-averaged, steady-state, incompressible Navier-Stokes equations, in which wind turbines are represented by surficial forces, are solved using a Control-Volume Finite Element Method (CVFEM). The fundamentals of developing a practical 3-D method are discussed in this paper with an emphasis on some of the challenges that arose during their implementation. For isolated turbines, results have indicated that the proposed 3-D method attains the same level of accuracy, in terms of performance predictions, as the previously developed 2-D axisymmetric method and the well-known momentum-strip theory. Furthermore, the capability of the proposed method to predict wind turbine wake characteristics is also illustrated. Satisfactory agreement with experimental measurements has been achieved. The analysis of a two-row periodic wind farm in neutral atmospheric boundary layers demonstrate the existence of positive interference effects (venturi effects) as well as the dominant influence of mutual interference on the performance of dense wind turbine clusters.