RANS simulations may not provide accurate results for all flow conditions. The interaction between a shock wave and a turbulent boundary layer is an example which may still be difficult to simulate accurately. Beside the inability to reproduce physical phenomena such as shock unsteadiness, the argument is put forward that the conventional numerical schemes, based on the Navier-Stokes equations, may be unable to generate a physically consistent turbulent stress tensor in the presence of large unresolved scales of motion. A large ratio between unresolved and resolved scales of motion, a sort of Knudsen number based on turbulent fluctuations, might introduce inaccuracies for which the turbulence model is not accountable. In order to improve the accuracy of RANS simulations, researchers have suggested various ad-hoc modifications to standard turbulence models which limit eddy viscosity or the turbulent stress tensor in the presence of strong gradients. Gas-kinetic schemes might be able to improve RANS predictions in shocklayers by removing or limiting the errors caused by the large scales ratio. These schemes are a class of their own; in the framework of a finite-volume or finite-elements discretizations, they model the numerical fluxes on the basis of the Boltzmann equation instead of the Navier-Stokes equations as is conventionally done. In practical terms, these schemes provide a higher accuracy and, more importantly, an in-built "multiscalar" mechanism, i.e. the ability to adjust to the size of unresolved scales of motion. This property makes them suitable for shock-capturing and rarefied flow. Gas-kinetic scheme may be coupled to a conventional RANS turbulence model; it is shown that the turbulent stress tensor is naturally adjusted as a function of the unresolved-toresolved scales ratios and achieves a higher physical consistency than conventional schemes. The simulations shown -well-known benchmark cases with strong shock-boundary layer interactions -have been obtained with a standard two-equation turbulence model (k-ω). It is shown that the gas-kinetic scheme provides good quality predictions, where conventional schemes with the same turbulence model are known to fail.