Computation of shock-dominated turbulent flows using Reynolds-averaged turbulence models may encounter large numerical error. This is due to the fact that source terms in the turbulence model equations have nonconservative derivatives of mean flow quantities. The nonconservative error takes large values at flow discontinuities. In this paper, we study the numerical characteristics of the k − ω turbulence model in canonical interaction of a normal shock with homogeneous isotropic turbulence. Several cases with varying shock strength are computed, and the results are compared with exact solution. The effect of grid refinement is also reported. The governing equations are cast in an equivalent conservation form that gives physically consistent results at a shock wave. The predicted turbulence amplification match the exact solution obtained by direct integration of the simplified equations.