A 2D multi-block high-speed compressible turbulent flow solver CFD2D based on the Jones and Launders two-equation k-ε turbulence model is developed. Method of solution employed is Finite Volume Method. Its basic algorithm is based on the approximate Riemann solver with the three-step Runge-Kutta time integration. Its additional feature includes Wilcox model for compressibility correction of k-ε turbulence model, Girmaji algebraic Reynolds stress (non-linear stress) model and linear stress model for evaluation of turbulent stresses. For validation purpose, code is applied to a 2D diamond aerofoil and a wedge ramp attached to a flat plate. CFD-predicted results are compared to the experimental results for shock wave and shock wave boundary layer interaction on the trailing edge of the fin. Contour plots are also compared to the Schlieren photographs. Flow simulation shows ability of the code to capture the physics of the flow both qualitatively and quantitatively.