An Analysis of Solution Point Coordinates for Flux Reconstruction Schemes on Tetrahedral Elements

The flux reconstruction (FR) approach offers an efficient route to high-order accuracy on unstructured grids. In this work we study the effect of solution point placement on the stability and accuracy of FR schemes on tetrahedral grids. To accomplish this we generate a large number of solution point candidates that satisfy various criteria at polynomial orders $$\wp = 3,4,5$$℘=3,4,5. We then proceed to assess their properties by using them to solve the non-linear Euler equations on both structured and unstructured meshes. The results demonstrate that the location of the solution points is important in terms of both the stability and accuracy. Across a range of cases it is possible to outperform the solution points of Shunn and Ham for specific problems. However, there appears to be a degree of problem-dependence with regards to the optimal point set, and hence overall it is concluded that the Shunn and Ham points offer a good compromise in terms of practical utility.

[1]  Antony Jameson,et al.  Application of High-Order Energy Stable Flux Reconstruction Schemes to the Euler Equations , 2011 .

[2]  Len Bos,et al.  On certain configurations of points in R n which are unisolvent for polynomial interpolation , 1991 .

[3]  Antony Jameson,et al.  On the Non-linear Stability of Flux Reconstruction Schemes , 2012, J. Sci. Comput..

[4]  Antony Jameson,et al.  Insights from von Neumann analysis of high-order flux reconstruction schemes , 2011, J. Comput. Phys..

[5]  Frank E. Ham,et al.  Symmetric quadrature rules for tetrahedra based on a cubic close-packed lattice arrangement , 2012, J. Comput. Appl. Math..

[6]  H. T. Huynh,et al.  A Flux Reconstruction Approach to High-Order Schemes Including Discontinuous Galerkin Methods , 2007 .

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

[8]  Antony Jameson,et al.  Nodal Points and the Nonlinear Stability of High-Order Methods for Unsteady Flow Problems on Tetrahedral Meshes , 2013 .

[9]  Zhi J. Wang,et al.  High-Order Multidomain Spectral Difference Method for the Navier-Stokes Equations , 2006 .

[10]  I. Doležel,et al.  Higher-Order Finite Element Methods , 2003 .

[11]  Freddie D. Witherden,et al.  PyFR: An open source framework for solving advection-diffusion type problems on streaming architectures using the flux reconstruction approach , 2013, Comput. Phys. Commun..

[12]  Freddie D. Witherden,et al.  On the identification of symmetric quadrature rules for finite element methods , 2014, Comput. Math. Appl..

[13]  Chi-Wang Shu,et al.  Total variation diminishing Runge-Kutta schemes , 1998, Math. Comput..

[14]  Beth A. Wingate,et al.  Several new quadrature formulas for polynomial integration in the triangle , 2005 .

[15]  Antony Jameson,et al.  Symmetric quadrature rules for simplexes based on sphere close packed lattice arrangements , 2014, J. Comput. Appl. Math..

[16]  Freddie D. Witherden,et al.  An Analysis of Solution Point Coordinates for Flux Reconstruction Schemes on Triangular Elements , 2014, J. Sci. Comput..