Both experiments and numerical simulations pertinent to the study of self-similarity in shock-induced turbulent mixing often do not cover sufficiently long enough times for the mixing layer to become developed in a fully turbulent manner. When the Mach number of the flow is sufficiently low, numerical simulations based on the compressible flow equations tend to become less accurate due to inherent numerical cancellation errors. This paper concerns a numerical study of the late time behaviour of single-shocked Richtmyer-Meshkov Instability (RMI) and associated compressible turbulent mixing using a new technique that addresses the above limitation. The present approach exploits the fact that RMI is a compressible flow during the early stages of the simulation and incompressible at late times. Therefore, depending on the compressibility of the flow field the most suitable model, compressible or incompressible, can be employed. This motivates the development of a hybrid compressible-incompressible solver that removes the low-Mach number limitations of the compressible solvers, thus allowing numerical simulations of late time mixing. Simulations have been performed for a multi-mode perturbation at the interface between two fluids of densities corresponding to an Atwood number of 0.5, and results are presented for the development of the instability, mixing parameters and turbulent kinetic energy spectra. The results are discussed in comparison with previous compressible simulations, theory and experiments.