This study is aimed at modelling of the unsteady interacting boundary layer around wind turbine rotor blades. It is shown that nonconservative mechanisms are present in this system, affecting the definition and approximation of solutions to the problem. The objective is to overcome the inconsistency of classical (dissipative) schemes, and find a suitable approximation method. A path-consistent discontinuous Galerkin method is proposed to correctly model the dynamics of this system. In the future, the current implementation is to be joined with the existing part of the RotorFlow simulation code. This will enable the incorporation of unsteady aerodynamic phenomena in the early design stage of new and large wind turbine rotors, which is currently unfeasible due to the high computational cost.