Accurate and reliable computation of the aerodynamic characteristics of wind turbines is very important for the development of new efficient designs. The flow around a wind turbine is modeled by a permeable disc (PD), solved through the Unsteady Reynolds-Averaged Navier–Stokes equations (URANS), here named PD/URANS method. The finite volume method and a total variation diminishing (TVD) scheme solve numerically the flow governing equations. The turbulent flow in the wake of the wind turbine is simulated utilising a one-equation turbulence model. The Glauert correction calculation considers a uniform normal force distribution (CT) on the virtual permeable disc applied to the flow, while the axial induction factor is obtained directly from the numerical solution of the URANS equations. The numerical axial induction factor obtained agrees fairly well with Glauert correction, except if the flow behind the turbine is highly unsteady and Reynolds number dependent.