Wind flow over forested areas is simulated using a 3-D Reynolds Averaged Navier–Stokes solver, employing the Wilcox k–ω model modified for atmospheric flows. Each tree is simulated individually with a two-region approach considering the leaf area and the trunk area. Each area affects the surrounding flow acting as a momentum sink term in the RANS equations. Leaf area sink terms are computed via Leaf Area Density which is described by an empirical relation. Predictions are compared against measurements and simulation results of other methods. A good agreement with measurements is observed as a result of a more detailed simulation of the canopy in comparison to simpler methods. A parametric analysis investigates the effect of the forest density as well as of the type and shape of the trees on the surrounding flow. A reduction of wind speed deficit up to 25% is predicted when the tree shape changes from cylinder to cone.