Scale-resolving simulations are viewed as powerful means for predicting complex turbulent flows, as often encountered in aeronautical applications. However, since turbulent scales span over a considerable range from the smallest Kolmogorov scales to the largest of equivalence to configuration size, scale-resolving computations are often demanding on computational resources and, furthermore, on the underlying numerical methods used in the simulations. Nonetheless, hybrid RANS (Reynolds-Averaged Navier-Stokes)-LES (Large-Eddy Simulation) techniques are considered computationally accurate and affordable for aeronautical industry applications. This thesis explores and develops numerical methods suitable for hybrid RANS-LES. These methods are implemented in the Computational Fluid Dynamics (CFD) solver M-Edge
