Spectral (Finite) Volume Method for Conservation Laws on Unstructured Grids. Basic Formulation

A high-order, conservative, yet efficient method named the spectral volume (SV) method is developed for conservation laws on unstructured grids. The concept of a “spectral volume” is introduced to achieve high-order accuracy in an efficient manner similar to spectral element and multidomain spectral methods. Each spectral volume is further subdivided into control volumes, and cell-averaged data from these control volumes are used to reconstruct a high-order approximation in the spectral volume. Then Riemann solvers are used to compute the fluxes at spectral volume boundaries. Cell-averaged state variables in the control volumes are updated independently. Furthermore, total variation diminishing and total variation bounded limiters are introduced in the SV method to remove/reduce spurious oscillations near discontinuities. Unlike spectral element and multidomain spectral methods, the SV method can be applied to fully unstructured grids. A very desirable feature of the SV method is that the reconstruction is carried out analytically, and the reconstruction stencil is always nonsingular, in contrast to the memory and CPU-intensive reconstruction in a high-order k-exact finite volume method. Fundamental properties of the SV method are studied and high-order accuracy is demonstrated for several model problems with and without discontinuities.

[1]  H. Keller,et al.  Analysis of Numerical Methods , 1969 .

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

[3]  Jaime Peraire,et al.  An implicit finite element method for high speed flows , 1991 .

[4]  B. Leer,et al.  Flux-vector splitting for the Euler equations , 1997 .

[5]  Chi-Wang Shu TVB uniformly high-order schemes for conservation laws , 1987 .

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

[7]  Chi-Wang Shu,et al.  Monotonicity Preserving Weighted Essentially Non-oscillatory Schemes with Increasingly High Order of Accuracy , 2000 .

[8]  Dimitri J. Mavriplis,et al.  Multigrid solution of the Navier-Stokes equations on triangular meshes , 1989 .

[9]  Zhi J. Wang,et al.  Anisotropic Cartesian grid method for viscous turbulent flow , 2000 .

[10]  Chi-Wang Shu Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[11]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[12]  S. Osher Riemann Solvers, the Entropy Condition, and Difference , 1984 .

[13]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. II. Monotonicity and conservation combined in a second-order scheme , 1974 .

[14]  Yannis Kallinderis,et al.  Hybrid Prismatic/Tetrahedral Grid Generation for Viscous Flows Around Complex Geometries , 1996 .

[15]  P. Roe Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[16]  Thomas Sonar,et al.  On Families of Pointwise Optimal Finite Volume ENO Approximations , 1998 .

[17]  S. Rebay,et al.  High-Order Accurate Discontinuous Finite Element Solution of the 2D Euler Equations , 1997 .

[18]  Bryan E. Richards,et al.  High resolution schemes for steady flow computation , 1991 .

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

[20]  Jay P. Boris,et al.  Flux-corrected transport. I. SHASTA, a fluid transport algorithm that works , 1973 .

[21]  Harold L. Atkins,et al.  A Finite-Volume High-Order ENO Scheme for Two-Dimensional Hyperbolic Systems , 1993 .

[22]  Yen Liu,et al.  Quadratic reconstruction finite volume schemes on 3D arbitrary unstructured polyhedral grids , 1999 .

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

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

[25]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[26]  A. Harten ENO schemes with subcell resolution , 1989 .

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

[28]  J. Steger,et al.  Flux vector splitting of the inviscid gasdynamic equations with application to finite-difference methods , 1981 .

[29]  Harold L. Atkins,et al.  QUADRATURE-FREE IMPLEMENTATION OF DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC EQUATIONS , 1996 .

[30]  I. Bohachevsky,et al.  Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics , 1959 .

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

[32]  C. Angelopoulos High resolution schemes for hyperbolic conservation laws , 1992 .

[33]  V. N. Venkatakrishnan,et al.  Implicit Solvers for Unstructured Meshes , 1993 .

[34]  P. Roe The 'optimum' upwind advection on a triangular mesh , 1990 .

[35]  Rémi Abgrall,et al.  On essentially non-oscillatory schemes on unstructured meshes: analysis and implementation , 1994 .

[36]  A. Khawajay,et al.  Hybrid Prismatic/tetrahedral Grid Generation for Complex Geometries , 1996 .

[37]  David L. Book,et al.  Flux-corrected transport II: Generalizations of the method , 1975 .

[38]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[39]  Z. Wang A fast flux-splitting for all speed flow , 1997 .

[40]  Stanley Osher,et al.  Upwind schemes and boundary conditions with applications to Euler equations in general geometries , 1983 .

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

[42]  David A. Kopriva,et al.  Multidomain spectral solution of the Euler Gas-dynamics equations , 1991 .

[43]  Rainald Loehner,et al.  On the computation of compressible turbulent flows on unstructured grids , 2000 .

[44]  Pierre Sagaut,et al.  A problem-independent limiter for high-order Runge—Kutta discontinuous Galerkin methods , 2001 .

[45]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[46]  S. Osher,et al.  A new class of high accuracy TVD schemes for hyperbolic conservation laws. [Total Variation Diminishing] , 1985 .

[47]  A. Patera A spectral element method for fluid dynamics: Laminar flow in a channel expansion , 1984 .

[48]  Meng-Sing Liou,et al.  Mass Flux Schemes and Connection to Shock Instability , 2000 .

[49]  John H. Kolias,et al.  A CONSERVATIVE STAGGERED-GRID CHEBYSHEV MULTIDOMAIN METHOD FOR COMPRESSIBLE FLOWS , 1995 .