In this paper we present a family of iterative methods to solve numerically second order elliptic problems with multi-scale data using multiple levels of grids. These methods are based upon the introduction of a Lagrange multiplier to enforce the continuity of the solution and its fluxes across interfaces. This family of methods can be interpreted as a mortar element method with complete overlapping domain decomposition for solving numerically multi-scale elliptic problems.