Local Mesh Refinement and Modeling of Faults and Pinchouts

Completely flexible grid block connections are necessary for local mesh refinement and modelling of faults and pinchouts. For practical purposes, it is necessary to combine such arbitrary connectivity with an adaptive implicit technique. In this paper, the methods used to discretize the flow equations for unusual geometries is discussed, and a rigorous analysis is given. A completely general sparse incomplete LU iterative solver is described which can solve Jacobians with arbitrary connectivity. Various test results are presented.