Multiblock mixed finite element methods for singularly perturbed problems (Q1590694)

A generic multiblock mixed finite element method is constructed for the singularly perturbed problem \[ -\varepsilon^2\Delta p+ ap= f\text{ in the square }[0,1]^2. \] The problem domain is decomposed into non-overlapping subdomains which separate the boundary layer subregions from other subregions. The global convergence of \(O(N^{-1})\), independent of \(\varepsilon\) for \(p\) and \(\varepsilon\nabla p\) in \(L^2\)-norm are obtained for the lowest Raviart-Thomas spaces on the specially designed piecewise rectangular mesh, where \(N\) is the number of divisions in each direction. Numerical results using an efficient domain decomposition algorithm are presented, which show that the special domain decomposition performs much better than a general decomposition.