The accomplishing of numerical flux conservation between initial grid and embedded grids highly enhances its reliability and efficiency in capturing discontinuous solution.