Design of a non-linear hp-multigrid scheme with static near-wall p-adaptation for RANS simulations in a high-order flux reconstruction framework