On the numerical approximation of the K-epsilon turbulence model for two dimensional compressible flows