We perform the first systematic study of the nonlinear evolution and scale coupling in Richtmyer-Meshkov (RM) flows induced by strong shocks. The smoothed particle hydrodynamics code (SPHC) is employed to ensure accurate shock capturing, interface tracking and accounting for the dissipation processes. We find that in strong-shock-driven RMI the background motion is supersonic. The amplitude of the initial perturbation strongly influences the flow evolution and the interfacial mixing that can be sub-sonic or supersonic. At late times the flow remains laminar rather than turbulent, and RM bubbles flatten and decelerate. In the fluid bulk, reverse cumulative jets appear and “hot spots” are formed—local heterogeneous microstructures with temperature substantially higher than that in the ambient. Our numerical simulations agree with the zero-order, linear, weakly nonlinear, and highly nonlinear theoretical analyses as well as with the experiments and suggest that the evolution of RMI is a multi-scale and heterogeneous process with a complicated character of scale coupling.