Mortar coupling for heterogeneous partial differential equations

We are interested in the approximation of 2D elliptic equations with dominated advection and featuring boundary layers. In order to reduce the computational complexity, the domain is split into two subregions, the first one far from the layer, where we can neglect the viscosity effects, and the second one next to the layer. In the latter domain the original elliptic equation is solved, while in the former one, the pure convection equation obtained by the original one by dropping the diffusive term is approximated.