A critical analysis of some popular methods for the discretisation of the gradient operator in finite volume methods

Finite volume methods (FVMs) constitute a popular class of methods for the numerical simulation of fluid flows. Among the various components of these methods, the discretisation of the gradient operator has received less attention despite its fundamental importance with regards to the accuracy of the FVM. The most popular gradient schemes are the divergence theorem (DT) (or Green-Gauss) scheme, and the least-squares (LS) scheme. Both are widely believed to be second-order accurate, but the present study shows that in fact the common variant of the DT gradient is second-order accurate only on structured meshes whereas it is zeroth-order accurate on general unstructured meshes, and the LS gradient is second-order and first-order accurate, respectively. This is explained through a theoretical analysis and is confirmed by numerical tests. The schemes are then used within a FVM to solve a simple diffusion equation on unstructured grids generated by several methods; the results reveal that the zeroth-order accuracy of the DT gradient is inherited by the FVM as a whole, and the discretisation error does not decrease with grid refinement. On the other hand, use of the LS gradient leads to second-order accurate results, as does the use of alternative, consistent, DT gradient schemes, including a new iterative scheme that makes the common DT gradient consistent at almost no extra cost. The numerical tests are performed using both an in-house code and the popular public domain PDE solver OpenFOAM.

[1]  Monica Oliveira,et al.  The Finite Volume Method in Computational Rheology , 2012 .

[2]  Jian Wu,et al.  Similarity and Comparison of Three Finite-Volume Methods for Diffusive Fluxes Computation on Nonorthogonal Meshes , 2013 .

[3]  Avram Sidi,et al.  Review of two vector extrapolation methods of polynomial type with applications to large-scale problems , 2012, J. Comput. Sci..

[4]  Dimitri J. Mavriplis,et al.  Revisiting the Least-squares Procedure for Gradient Reconstruction on Unstructured Meshes , 2003 .

[5]  C. Ollivier-Gooch,et al.  A high-order-accurate unstructured mesh finite-volume scheme for the advection-diffusion equation , 2002 .

[6]  Boris Diskin,et al.  Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Inviscid Fluxes , 2011 .

[7]  D. K. Walters,et al.  Validation of OpenFOAM numerical methods and turbulence models for incompressible bluff body flows , 2015 .

[8]  Robert Schneiders,et al.  Quadrilateral and Hexahedral Element Meshes , 2002 .

[9]  Carl Ollivier-Gooch,et al.  Accuracy analysis of unstructured finite volume discretization schemes for diffusive fluxes , 2014 .

[10]  Joachim Sch.  NETGEN An advancing front 2D/3D-mesh generator based on abstract rules , 1997 .

[11]  S. A. Mohsen Karimian,et al.  Discretization and parallel performance of an unstructured finite volume Navier–Stokes solver , 2006 .

[12]  I. Demirdžić.  On the Discretization of the Diffusion Term in Finite-Volume Continuum Mechanics , 2015 .

[13]  John G. Bartzis,et al.  Numerical experiments on the efficiency of local grid refinement based on truncation error estimates , 2012, J. Comput. Phys..

[14]  S. Muzaferija,et al.  Adaptive finite volume method for flow prediction using unstructured meshes and multigrid approach , 1994 .

[15]  C. Kiris,et al.  Gradient Calculation Methods on Arbitrary Polyhedral Unstructured Meshes for Cell-Centered CFD Solvers , 2014 .

[16]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[17]  Qian Wang,et al.  Compact high order finite volume method on unstructured grids III: Variational reconstruction , 2017, J. Comput. Phys..

[18]  Timothy J. Barth,et al.  The design and application of upwind schemes on unstructured meshes , 1989 .

[19]  J. Tsamopoulos,et al.  Transient flow of gravity-driven viscous films over substrates with rectangular topographical features , 2016 .

[20]  Yannis Kallinderis,et al.  A priori mesh quality estimation via direct relation between truncation error and mesh distortion , 2009, J. Comput. Phys..

[21]  Alexandros Syrakos,et al.  Estimate of the truncation error of a finite volume discretisation of the Navier-Stokes equations on colocated grids , 2015, ArXiv.

[22]  S. Muzaferija,et al.  Finite-Volume CFD Procedure and Adaptive Error Control Strategy for Grids of Arbitrary Topology , 1997 .

[23]  Volker Seidl,et al.  AN IMPLICIT FINITE-VOLUME METHOD USING NONMATCHING BLOCKS OF STRUCTURED GRID , 1997 .

[24]  Boris Diskin,et al.  Notes on accuracy of finite-volume discretization schemes on irregular grids , 2010 .

[25]  T. Barth.  A 3-D upwind Euler solver for unstructured meshes , 1991 .

[26]  J. Tsamopoulos,et al.  A quasi-elliptic transformation for moving boundary problems with large anisotropic deformations , 2003 .

[27]  C. Ollivier-Gooch,et al.  An efficient implicit unstructured finite volume solver for generalised Newtonian fluids , 2016 .

[28]  F. Bramkamp,et al.  An adaptive multiscale finite volume solver for unsteady and steady state flow computations , 2004 .

[29]  Eiji Shima,et al.  New Gradient Calculation Method for MUSCL Type CFD Schemes in Arbitrary Polyhedra , 2010 .

[30]  Hiroaki Nishikawa,et al.  Effects of high-frequency damping on iterative convergence of implicit viscous solver , 2017, J. Comput. Phys..

[31]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

[32]  Anthony G. Straatman,et al.  An accurate gradient and Hessian reconstruction method for cell‐centered finite volume discretizations on general unstructured grids , 2009 .

[33]  Yannis Dimakopoulos,et al.  On the elliptic mesh generation in domains containing multiple inclusions and undergoing large deformations , 2009, J. Comput. Phys..

[34]  Georgios C. Georgiou,et al.  Solution of the square lid-driven cavity flow of a Bingham plastic using the finite volume method , 2013, ArXiv.

[35]  Weiwei Zhang,et al.  Accuracy preserving limiter for the high-order finite volume method on unstructured grids , 2017 .

[36]  Joe F. Thompson,et al.  Numerical grid generation: Foundations and applications , 1985 .

[37]  James L. Thomas,et al.  Accuracy of Gradient Reconstruction on Grids with High Aspect Ratio , 2013 .

[38]  J. Tsamopoulos,et al.  Transient flow of gravity-driven viscous films over 3D patterned substrates: conditions leading to Wenzel, Cassie and intermediate states , 2017 .

[39]  Kwan-Liu Ma,et al.  A Comparison of Gradient Estimation Methods for Volume Rendering on Unstructured Meshes , 2011, IEEE Transactions on Visualization and Computer Graphics.

[40]  Didier Vigneron,et al.  An implicit high order finite volume scheme for the solution of 3D Navier-Stokes equations with new discretization of diffusive terms , 2008 .

[41]  Hiroaki Nishikawa.  Robust and accurate viscous discretization via upwind scheme – I: Basic principle , 2011 .

[42]  Yves Marcel Ahipo,et al.  A robust and efficient finite volume scheme for the discretization of diffusive flux on extremely skewed meshes in complex geometries , 2009, J. Comput. Phys..

[43]  Audra E. Kosh,et al.  Linear Algebra and its Applications , 1992 .

[44]  F. Pinho,et al.  Numerical simulation of non-linear elastic flows with a general collocated finite-volume method , 1998 .

[45]  Hrvoje Jasak,et al.  Error analysis and estimation for the finite volume method with applications to fluid flows , 1996 .