The direct discontinuous Galerkin (DDG) viscous flux scheme for the high order spectral volume method

Abstract The direct discontinuous Galerkin (DDG) method was developed by Liu and Yan to discretize the diffusion flux. It was implemented for the discontinuous Galerkin (DG) formulation. In this paper, we perform four tasks: (i) implement the direct discontinuous Galerkin (DDG) scheme for the spectral volume method (SV) method, (ii) design and implement two variants of DDG (called DDG2 and DDG3) for the SV method, (iii) perform a Fourier type analysis on both methods when solving the 1D diffusion equation and combine the above with a non-linear global optimizer, to obtain modified constants that give significantly smaller errors (in 1D), (iv) use the above coefficients as starting points in 2D. The dissipation properties of the above schemes were then compared with existing flux formulations (local discontinuous Galerkin, Penalty and BR2). The DDG, DDG2 and DDG3 formulations were found to be much more accurate than the above three existing flux formulations. The accuracy of the DDG scheme is heavily dependent on the penalizing coefficient for the odd ordered schemes. Hence a loss of accuracy was observed even for mildly non-uniform grids for odd ordered schemes. On the other hand, the DDG2 and DDG3 schemes were mildly dependent on the penalizing coefficient for both odd and even orders and retain their accuracy even on highly irregular grids. Temporal analysis was also performed and this yielded some interesting results. The DDG and its variants were implemented in 2D (on triangular meshes) for Navier–Stokes equations. Even the non-optimized versions of the DDG displayed lower errors than the existing schemes (in 2D). In general, the DDG and its variants show promising properties and it indicates that these approaches have a great potential for higher dimension flow problems.

[1]  Ravishekar Kannan An implicit LU‐SGS spectral volume method for moment models in device simulations II: Accuracy studies and performance enhancements using the penalty and BR2 formulations , 2011 .

[2]  Hailiang Liu,et al.  The Direct Discontinuous Galerkin (DDG) Methods for Diffusion Problems , 2008, SIAM J. Numer. Anal..

[3]  Zhi J. Wang,et al.  Extension of the spectral volume method to high-order boundary representation , 2006 .

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

[5]  Chi-Wang Shu,et al.  Runge–Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems , 2001, J. Sci. Comput..

[6]  F. Brezzi,et al.  Discontinuous Galerkin approximations for elliptic problems , 2000 .

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

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

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

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

[11]  Ravishekar Kannan,et al.  An implicit LU‐SGS spectral volume method for the moment models in device simulations: Formulation in 1D and application to a p‐multigrid algorithm , 2010 .

[12]  S. Rebay,et al.  A High-Order Accurate Discontinuous Finite Element Method for the Numerical Solution of the Compressible Navier-Stokes Equations , 1997 .

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

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

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

[16]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[17]  D. Arnold An Interior Penalty Finite Element Method with Discontinuous Elements , 1982 .

[18]  M. Liou,et al.  A New Flux Splitting Scheme , 1993 .

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

[20]  S. Dreyfus,et al.  Thermodynamical Approach to the Traveling Salesman Problem : An Efficient Simulation Algorithm , 2004 .

[21]  Yuzhi Sun,et al.  Spectral (finite) volume method for conservation laws on unstructured grids VI: Extension to viscous flow , 2006, J. Comput. Phys..

[22]  Marcel Vinokur,et al.  Spectral (finite) volume method for conservation laws on unstructured grids V: Extension to three-dimensional systems , 2006, J. Comput. Phys..

[23]  Chi-Wang Shu,et al.  The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems , 1998 .

[24]  R. Kannan,et al.  A Study of Viscous Flux Formulations for a p-Multigrid Spectral Volume Navier Stokes Solver , 2009, J. Sci. Comput..

[25]  M. Wheeler An Elliptic Collocation-Finite Element Method with Interior Penalties , 1978 .

[26]  Mengping Zhang,et al.  AN ANALYSIS OF THREE DIFFERENT FORMULATIONS OF THE DISCONTINUOUS GALERKIN METHOD FOR DIFFUSION EQUATIONS , 2003 .

[27]  Z. J. Wang,et al.  Efficient Implicit Non-linear LU-SGS Approach for Compressible Flow Computation Using High-Order Spectral Difference Method , 2008 .

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

[29]  Zhi J. Wang,et al.  Spectral (finite) volume method for conservation laws on unstructured grids IV: extension to two-dimensional systems , 2004 .

[30]  Chunlei Liang,et al.  A p-multigrid spectral difference method with explicit and implicit smoothers on unstructured triangular grids , 2009 .

[31]  Zhi Jian Wang,et al.  Spectral (Finite) Volume Method for Conservation Laws on Unstructured Grids III: One Dimensional Systems and Partition Optimization , 2004, J. Sci. Comput..

[32]  V. Cerný Thermodynamical approach to the traveling salesman problem: An efficient simulation algorithm , 1985 .