This chapter explains how the mathematical models from Chapter 4 are implemented and integrated to form a full simulator. To this end, we introduce data structures to represent fluid behavior, the reservoir state, boundary conditions, source terms, and wells. We then explain in detail how the two-point flux approximation (TPFA) scheme is implemented in MRST for general unstructured grids. We also outline the basic solver used to compute time-of-flight and tracer partitions. We end the chapter by presenting a few examples that demonstrate how to set up simulations in MRST and set appropriate boundary conditions, source terms, or well models. The examples include the famous quarter-five spot problem, a corner-point grid with four intersecting faults, and a model of a shallow-marine reservoir (SAIGUP).

The chapter explains how you can generate grid models to represent subsurface reservoirs. We outline a number of elementary grid types: structured/rectilinear grids, fictitious domains, Delaunay triangulations, and Voronoi grids. We then explain stratigraphic grids that are commonly used to model real subsurface formations, including in particular corner-point and perpendicular bisector (PEBI) grids. We explain how such grids are represented in MRST using a data structure for general unstructured grids, and we discuss how to compute geometric properties like volumes, face areas, face normals, etc. We end the chapter by presenting an overview of alternative gridding techniques, including composite grids, multiblock grids, and control-point and boundary conformal grids.