In this study, a modelling methodology is proposed for RANS simulations of neutral Atmospheric Boundary Layer (ABL) flows on the basis of the standard k - ε model, which allows the adoption of an arbitrary shear stress model. This modelling methodology is first examined in the context of an open flat terrain in an empty domain to ascertain there are no substantial changes in the prescribed profiles. The results show that relatively good homogeneity can be achieved with this modelling methodology for various sets of inflow boundary profiles. In addition, to extend the solutions derived from the standard k - ε model to RNG k - ε model, the RNG k - ε model is in detail assembly and tuned. Finally, the topographic effects on surface wind speeds over a complex terrain are assessed with the combined use of the proposed methodology and the modified RNG model. The numerical results are in good agreement with wind tunnel testing results and long-term field observations. A discussion of the effects of horizontal homogeneity and turbulence models on the simulated wind flows over a complex terrain is also given.