Grid generation for reservoir simulation must honor classical key constraints and ensure boundary alignment such that control-volume boundaries are aligned with geological features including layers, shale barriers, fractures, faults, pinch-outs, and multilateral wells. Novel unstructured grid generation methods are proposed that automate control-volume and/or control point boundary alignment and yield perpendicular-bisector (PEBI) meshes both with respect to primal and dual (essentially PEBI) cells. In order to honor geological features in the primal configuration, we introduce the idea of protection circles that contain segments of key geological boundaries, while in order to generate a dual-cell feature aligned grid, we construct halos around key geological features. The grids generated are employed to study comparative performance of cell-centred versus cell-vertex flux-continuous control-volume distributed multi-point flux approximation (CVD-MPFA) finite-volume formulations using equivalent degrees of freedom and thus ensure application of the most efficient methods. The CVD-MPFA formulation (c.f. Edwards et al.) in cell-centred and cell-vertex modes is somewhat analogous and requires switching control-volume from primal to dual or vice versa, together with appropriate data structures and boundary conditions, however dual-cells are generated after primal grid generation. The relative benefits of both types of approximation, i.e., cell-centred versus vertex-centred, are contrasted in terms of flow resolution and degrees of freedom required
