Nonoverlapping domain decomposition method with preconditioner from asymptotic analysis of steady flow in high contrast media

Abstract

ABSTRACT We present a nonoverlapping domain decomposition method for steady flow in high contrast heterogeneous media modelled by an elliptic equation with coefficients that have very large amplitude variations on a small spatial scale. The linear system of equations resulting from matching the solution trace and the fluxes through the boundary of the subdomains is ill-conditioned, especially for fine meshes needed to capture the rapid variations of the solution. Our main contribution is to show with analysis and numerical simulations how to use an asymptotic approximation of the Dirichlet to Neumann map of the sub-domain problems to obtain a preconditioner for an efficient domain decomposition algorithm.