A zonal RANS-LES method is presented and applied to the unsteady flow around an airfoil at high angle of attack. The attached boundary layer is simulated by RANS whereas the laminar separation bubble, the laminar-to-turbulent transition, and the trailing-edge separation are computed using LES. The general RANS domain is connected to the LES domain by overlapping regions. The RANS-to-LES transition for the turbulent boundary layer couples a turbulent inflow generation method with controlled forcing to ensure a fast and smooth transition. The transition from LES to RANS applies a reconstruction technique for the turbulent viscosity combined with a forcing layer to attain the correct boundary layer velocity profile in the RANS domain. This fully coupled zonal method is applied to predict the characteristics of an airfoil at high angle of attack, including the unsteady phenomena occurring near stall. The zonal method reduces the computational costs by a factor of 4 and its results show good agreement with the pure LES findings and experimental data.