In this research, interaction of a shock wave with incoming turbulent flow is investigated. To this end, different turbulent flows with different intensities and integral length scales are generated and impact of these turbulent flows on a shock wave are examined. In this study, two-dimensional Navier-Stokes equation is numerically solved using a high order spatial-temporal discretization. For spatial discretization, two different methods are implemented. In stream wise direction, i.e. perpendicular to the shock wave, a sixth-order accurate essentially non-oscillatory method (ENO) has been used which is able to capture the shock wave. In spanwise direction, i.e. parallel to the shock wave, a sixth-order Pade scheme has been used which is able to accurately capture small scale flow field structures. Time integration is performed using a third-order Runge-Kutta method. Overall, it has been observed that the turbulent kinetic energy increases across the shock wave. Fluctuations with larger integral length scale show higher turbulent kinetic energy increase across the shock wave. Further, it has been observed that although the integral length scale of the upstream fluctuations does not influence the location of the shock, the intensity of the upstream fluctuations have a profound effect on the shock wave location.