“…Cell-centred and vertex-centred second-order finite volume (FV) methods are still the predominant techniques used in commercial and industrial computational fluid dynamics (CFD) solvers due to their robustness, easy implementation and relatively low cost [2,14,15,17,19,33]. Both, cellcentred and vertex-centred, FV techniques require a reconstruction of the gradient of the solution to ensure second-order convergence of the unknown and first-order convergence of the fluxes [1,3,11,12]. The accuracy of the scheme is therefore dependent on the accuracy of the reconstruction technique, which in turns depends on the quality of the mesh.…”