A family of first-order accurate gradient schemes for finite volume methods

A new discretisation scheme for the gradient operator, suitable for use in second-order accurate Finite Volume Methods (FVMs), is proposed. The derivation of this scheme, which we call the Taylor-Gauss (TG) gradient, is similar to that of the least-squares (LS) gradients, whereby the values of the differentiated variable at neighbouring cell centres are expanded in truncated Taylor series about the centre of the current cell, and the resulting equations are summed after being weighted by chosen vectors. Unlike in the LS gradients, the TG gradients use vectors aligned with the face normals, resembling the Green-Gauss (GG) gradients in this respect. Thus, the TG and LS gradients belong in a general unified framework, within which other gradients can also be derived. The similarity with the LS gradients allows us to try different weighting schemes (magnitudes of the weighting vectors) such as weighting by inverse distance or face area. The TG gradients are tested on a variety of grids such as structured, locally refined, randomly perturbed, and with high aspect ratio. They are shown to be at least first-order accurate in all cases, and are thus suitable for use in second-order accurate FVMs. In many cases they compare favourably over existing schemes.

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

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

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

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

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

[6]  M. Darwish,et al.  The Finite Volume Method in Computational Fluid Dynamics: An Advanced Introduction with OpenFOAM® and Matlab , 2015 .

[7]  J. Tsamopoulos,et al.  A finite volume method for the simulation of elastoviscoplastic flows and its application to the lid-driven cavity case , 2018, Journal of Non-Newtonian Fluid Mechanics.

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

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

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

[11]  Ramakrishnan Thirumalaisamy,et al.  A new Green-Gauss reconstruction on unstructured meshes. Part I: Gradient reconstruction , 2020, J. Comput. Phys..

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

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

[14]  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..

[15]  Alexandros Syrakos,et al.  A critical analysis of some popular methods for the discretisation of the gradient operator in finite volume methods , 2016 .

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

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

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

[19]  M. A. Alves,et al.  Stabilization of an open-source finite-volume solver for viscoelastic fluid flows , 2017 .

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

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

[22]  R. Ma,et al.  Accuracy analysis of gradient reconstruction on isotropic unstructured meshes and its effects on inviscid flow simulation , 2019, Advances in Aerodynamics.

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

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