Numerical simulations of flow fields around the wind turbine rotor simplified as an actuator disk (AD) with zero thickness have been made to investigate the flow structure and wake development in different operation states. A N-S solver has been used and the energy extracted by the rotor is represented by a discontinuous pressure jump through the actuator disk. Axial pressure and velocity development from far upstream to far downstream is fully described by the simulations, which could never be obtained by the momentum theory. It is showed that there are significant differences in wake development between inviscid and viscous conditions. In inviscid simulations, the axial velocity keeps decreasing along the oncoming flow direction, which is consistent with the momentum theory. In viscous simulations, however, the axial velocity first decreases but then gradually recovers approaching to the undisturbed velocity, due to momentum transport from outer flow to wake flow by viscous shear effect. Based on the numerical analysis, the work of this paper is also focused on wake modeling. A new two-dimensional models based on nonlinear wake development has been developed, which is capable to describe the far wake more accurately.