The modeling of flow in fractured reservoirs leads to a problem which can rather naturally be divided into subdomains and the interface regions between subdomains. Greater flexibility in the numerical approximation is achieved if the subdomains can be independently meshed. This leads to meshes which do not conform at the interfaces between subdomains. The question of how to generate consistent approximation on nonmatching grids will be considered in this talk. In particular, mortar methods for coupling discretizations of different types in different subdomains will be considered. In addition, iterative algorithms for solving the composite system of equations will be discussed.