A numerical study of underexpanded jet and impingement on a wall mounted at various distances from the nozzle exit is presented. The 3-dimensional Wavier-Stokes equations and κ-ω turbulence equations are solved. The grids are constructed as overlapped grid systems to examine the distance effect. The DADI method is applied to obtain steady-state solutions. To avoid numerical instability such as the carbuncle phenomena that sometimes accompany approximate Riemann solver, the HLLE+ scheme is employed for the inviscid flux at the cell interfaces. A goal of this work is to apply a number of two-equation turbulence models based on the w equation to the impinging jet problem.