Wind turbine layout optimization in wind farm is one of the most important technologies to increase the wind power utilization. This paper studies the wind turbine layout optimization with multiple hub heights wind turbines using greedy algorithm. The linear wake model and the particle wake model are used for wake flow calculation over flat terrain and complex terrain, respectively. Three-dimensional greedy algorithm is developed to optimize wind turbine layout with multiple hub heights for minimizing cost per unit power output. The numerical cases over flat terrain and complex terrain are used to validate the effectiveness of the proposed greedy algorithm for the optimization problem. The results reveal that it incurs lower computational costs to obtain better optimized results using the proposed greedy algorithm than the one using genetic algorithm. Compared to the layout with identical hub height wind turbines, the one with multiple hub height wind turbines can increase the total power output and decrease the cost per unit power output remarkably, especially for the wind farm over complex terrain. It is suggested that three-dimensional greedy algorithm is an effective method for more benefits of using wind turbines with multiple hub heights in wind farm design. (C) 2016 Elsevier Ltd. All rights reserved.