Interface solvers and preconditioners of domain decomposition type for multiphase flow in multiblock porous media

A multiblock mortar approach to modeling multiphase flow in porous media decomposes the simulation domain into a series of blocks with possibly different physical and numerical models employed in each block, Matching conditions along the interfaces are imposed through the use of mortar finite elements. A parallel domain decomposition algorithm reduces the algebraic nonlinear system to an interface problem which is solved via a nonlinear multigrid with Newton-GMRES smoothing.