The operation of the wind turbines downstream is affected by the wake of the wind turbines upstream. Wind turbine wake flow is investigated by applying the actuator disc (AD) method. The modified k-ε turbulence model is proposed by using both the turbulent kinetic energy source term and the dissipation rate source term to improve the standard k-ε turbulence model for coordinating the generation and the dissipation of the turbulent kinetic energy. The dissipation rate parameter C 4ε that obeys a parabolic distribution is used, based on theoretical analysis. The force distributed on the AD is also used instead of a constant, as used in the classical AD method. The simulation results were consistent with the measurements that correspond to different kinds of wind turbines and conditions. The nacelle and the inflow turbulence intensity have great influences on accurately simulating the wake, so it is necessary to imitate the rotor along with the nacelle and accurately measure the inflow turbulence intensity. been continually developing [11]. The engineering "wake model" was established earlier and has a simple structure and short calculation time. The most famous Jensen model, which is also the basic model of the commonly used commercial wind resources software WAsP, was proposed to evaluate the wind turbine wakes on flat terrain [12]. The eddy viscosity model appeared as another analytical wind turbine wakes model and was used in the commercial wind resources software WindFarmer belonging to the Garrad Hassan company [13]. These models are based on assumptions of the flow behavior, such as axisymmetry and self-similarity. They work fast and are suitable for engineering applications, but they cannot predict wake accurately. The blade model can be constructed directly in the flow field and used to precisely calculate the load of the wind turbine and other aspects related to wind turbine performance [14]. Qian Y.R. used this method to predict the air-load with different wind speeds under yaw conditions [15]. However, this method needs more grids around the blade and a lot of computing resources; therefore, it does not apply to multiple wind turbines or wind farm scales. The actuator line model simplifies the blade into a line which is rotated as the blade, then constructs a functional relationship between the aerodynamic data of the blade's local section airfoil and body force, so that it can simulate the interaction between the blade and the inflow [16]. The actuator line model, combined with the large eddy simulation (LES) method, was used to investigate the near wake region and explain the reason why the velocity deficit and turbulence intensity increase nearby the rotor [17]. Storey R.C. coupled the model with multi-physics method in order to research the wind turbine performance under transient conditions [18]. The actuator line model also can demonstrate the wake of the wind turbine over complex terrain and reduce the cost of wind farm simulation to a certain extent [19], but it still takes a lot of computing res...