When supersonic underexpanded jets are exhausted from the nozzle, complex shock cell configurations such as barrel shock, expansion fan, Mach disc, and exhaust-gas jet boundary are appeared repetitively. The shock cell is smeared by turbulence dissipation and disappeared in long distance from the nozzle. When underexpanded jet is suddenly impinged on a flat plate, it forms very complex flow structure. In this paper, we solve compressible Wavier-Stokes equation adapting finite volume method to obtain jet impingement flow structure and compare calculated data with experimental ones. It is shown that numerical simulation data are in good agreement with experimental one in a short distance between nozzle exit and flat plate and little influence of underexpanded ratio is appeared in jet impingement now distribution.