In order to realize high order accuracy in an overset grid system, a high order interpolation method was developed based on Finite Volume Method (FVM). To improve stability in calculating a nonlinear discontinuity, it was combined with Multi-dimensional Limiting Process (MLP) to control oscillatory phenomena. Thus, it can be robustly applied to a discontinuous region as well as a continuous region and shows good results. The developed method was compared with bilinear interpolation method throughout some test cases and verified to be more accurate than conventional methods. Nomenclature c = speed of sound e = internal energy te = total energy M = Mach number p = pressure r = core radius u = velocity in the x direction θV = velocity in the θ direction v = velocity in the y direction α = MLP limiting variable ω = angular velocity γ = specific heat ratio