A multi-moment finite volume method for incompressible Navier-Stokes equations on unstructured grids: Volume-average/point-value formulation

A robust and accurate finite volume method (FVM) is proposed for incompressible viscous fluid dynamics on triangular and tetrahedral unstructured grids. Differently from conventional FVM where the volume integrated average (VIA) value is the only computational variable, the present formulation treats both VIA and the point value (PV) as the computational variables which are updated separately at each time step. The VIA is computed from a finite volume scheme of flux form, and is thus numerically conservative. The PV is updated from the differential form of the governing equation that does not have to be conservative but can be solved in a very efficient way. Including PV as the additional variable enables us to make higher-order reconstructions over compact mesh stencil to improve the accuracy, and moreover, the resulting numerical model is more robust for unstructured grids.We present the numerical formulations in both two and three dimensions on triangular and tetrahedral mesh elements. Numerical results of several benchmark tests are also presented to verify the proposed numerical method as an accurate and robust solver for incompressible flows on unstructured grids.

[1]  Hrvoje Jasak,et al.  OpenFOAM: Open source CFD in research and industry , 2009 .

[2]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .

[3]  A. Chorin.  Numerical solution of the Navier-Stokes equations , 1968 .

[4]  Hrvoje Jasak,et al.  A tensorial approach to computational continuum mechanics using object-oriented techniques , 1998 .

[5]  Esteban Ferrer,et al.  A high order Discontinuous Galerkin Finite Element solver for the incompressible Navier-Stokes equations , 2011 .

[6]  R. Akoh,et al.  A CIP/multi-moment finite volume method for shallow water equations with source terms , 2008 .

[7]  Feng Xiao,et al.  CIP/multi-moment finite volume method for Euler equations: A semi-Lagrangian characteristic formulation , 2007, J. Comput. Phys..

[8]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[9]  P. Frederickson,et al.  Higher order solution of the Euler equations on unstructured grids using quadratic reconstruction , 1990 .

[10]  H. Kuhlmann,et al.  Accurate three-dimensional lid-driven cavity flow , 2005 .

[11]  B. V. Leer,et al.  Towards the Ultimate Conservative Difference Scheme , 1997 .

[12]  B. P. Leonard,et al.  A stable and accurate convective modelling procedure based on quadratic upstream interpolation , 1990 .

[13]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[14]  Chaowei Hu,et al.  No . 98-32 Weighted Essentially Non-Oscillatory Schemes on Triangular Meshes , 1998 .

[15]  Riccardo Rossi,et al.  Direct numerical simulation of scalar transport using unstructured finite-volume schemes , 2009, J. Comput. Phys..

[16]  A. Gosman,et al.  Solution of the implicitly discretised reacting flow equations by operator-splitting , 1986 .

[17]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[18]  Shahrouz Aliabadi,et al.  Development of a hybrid finite volume/element solver for incompressible flows , 2007 .

[19]  Zhi J. Wang,et al.  Spectral (Finite) Volume Method for Conservation Laws on Unstructured Grids. Basic Formulation , 2002 .

[20]  Marek Stastna,et al.  A short note on the discontinuous Galerkin discretization of the pressure projection operator in incompressible flow , 2013, J. Comput. Phys..

[21]  Roland Glowinski,et al.  Wall-driven incompressible viscous flow in a two-dimensional semi-circular cavity , 2006, J. Comput. Phys..

[22]  J. P. V. Doormaal,et al.  ENHANCEMENTS OF THE SIMPLE METHOD FOR PREDICTING INCOMPRESSIBLE FLUID FLOWS , 1984 .

[23]  Feng Xiao,et al.  Shallow water model on cubed-sphere by multi-moment finite volume method , 2008, J. Comput. Phys..

[24]  I. Demirdzic,et al.  Fluid flow and heat transfer test problems for non‐orthogonal grids: Bench‐mark solutions , 1992 .

[25]  O. Friedrich,et al.  Weighted Essentially Non-Oscillatory Schemes for the Interpolation of Mean Values on Unstructured Grids , 1998 .

[26]  B. Vanleer,et al.  Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection , 1977 .

[27]  D. Spalding,et al.  A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows , 1972 .

[28]  Chi-Wang Shu.  Total-variation-diminishing time discretizations , 1988 .

[29]  Feng Xiao,et al.  A 4th-order and single-cell-based advection scheme on unstructured grids using multi-moments , 2005, Comput. Phys. Commun..

[30]  Hong Li,et al.  A hybrid vertex-centered finite volume/element method for viscous incompressible flows on non-staggered unstructured meshes , 2012 .

[31]  C. Ross Ethier,et al.  A high-order discontinuous Galerkin method for the unsteady incompressible Navier-Stokes equations , 2007, J. Comput. Phys..

[32]  A. D. Gosman,et al.  Error Analysis of the Finite-Volume Method with Respect to Mesh Type , 2010 .

[33]  Takashi Yabe,et al.  A universal solver for hyperbolic equations by cubic-polynomial interpolation I. One-dimensional solver , 1991 .

[34]  F. Xiao,et al.  Numerical simulations of free-interface fluids by a multi-integrated moment method , 2005 .

[35]  R. Rannacher,et al.  Benchmark Computations of Laminar Flow Around a Cylinder , 1996 .

[36]  Chi-Wang Shu,et al.  The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. IV. The multidimensional case , 1990 .

[37]  T. Yabe,et al.  The constrained interpolation profile method for multiphase analysis , 2001 .

[38]  T. Yabe,et al.  Completely conservative and oscillationless semi-Lagrangian schemes for advection transportation , 2001 .

[39]  Feng Xiao,et al.  Unified formulation for compressible and incompressible flows by using multi-integrated moments II: Multi-dimensional version for compressible and incompressible flows , 2006, J. Comput. Phys..

[40]  P. Wesseling,et al.  Benchmark solutions for the incompressible Navier–Stokes equations in general co‐ordinates on staggered grids , 1993 .

[41]  J. Hesthaven,et al.  Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications , 2007 .

[42]  Z. Wang,et al.  Spectral (Finite) Volume Method for Conservation Laws on Unstructured Grids , 2002 .

[43]  Jungwoo Kim,et al.  An immersed-boundary finite-volume method for simulations of flow in complex geometries , 2001 .

[44]  V. C. Patel,et al.  Flow past a sphere up to a Reynolds number of 300 , 1999, Journal of Fluid Mechanics.

[45]  Feng Xiao,et al.  Unified formulation for compressible and incompressible flows by using multi-integrated moments I: one-dimensional inviscid compressible flow , 2004 .

[46]  K. Squires,et al.  LES and DES investigations of turbulent flow over a sphere , 2000 .

[47]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: one-dimensional systems , 1989 .

[48]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .