Spectral difference method for unstructured grids I: Basic formulation

A new, high-order, conservative, and efficient method for conservation laws on unstructured grids is developed. It combines the best features of structured and unstructured grid methods to attain computational efficiency and geometric flexibility; it utilizes the concept of discontinuous and high-order local representations to achieve conservation and high accuracy; and it is based on the finite-difference formulation for simplicity. Universal reconstructions are obtained by distributing unknown and flux points in a geometrically similar manner for all unstructured cells. Placements of these points with various orders of accuracy are given for the triangular elements. Accuracy studies of the method are carried out with the two-dimensional linear wave equation and Burgers' equation, and each order of accuracy is verified numerically. Numerical solutions of plane electromagnetic waves incident on perfectly conducting circular cylinders are presented and compared with the exact solutions to demonstrate the capability of the method. Excellent agreement has been found. The method is much simpler than the discontinuous Galerkin and spectral volume methods for unstructured grids.

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

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

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

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

[5]  Jan S. Hesthaven,et al.  Regular Article: Spectral Collocation Time-Domain Modeling of Diffractive Optical Elements , 1999 .

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

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

[8]  R. F. Warming,et al.  An implicit finite-difference algorithm for hyperbolic systems in conservation-law form. [application to Eulerian gasdynamic equations , 1976 .

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

[10]  A. Jameson,et al.  Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes , 1981 .

[11]  David A. Kopriva,et al.  A Staggered-Grid Multidomain Spectral Method for the Compressible Navier-Stokes Equations , 1998 .

[12]  J. Hesthaven,et al.  Nodal high-order methods on unstructured grids , 2002 .

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

[14]  D. Gottlieb,et al.  Numerical analysis of spectral methods : theory and applications , 1977 .

[15]  D. Kopriva.  A spectral multidomain method for the solution of hyperbolic systems , 1986 .

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

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

[18]  Jan S. Hesthaven,et al.  From Electrostatics to Almost Optimal Nodal Sets for Polynomial Interpolation in a Simplex , 1998 .

[19]  R. Maccormack,et al.  Computational efficiency achieved by time splitting of finite difference operators. , 1972 .

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

[21]  Eric F Darve,et al.  Author ' s personal copy A hybrid method for the parallel computation of Green ’ s functions , 2009 .

[22]  Mark A. Taylor,et al.  An Algorithm for Computing Fekete Points in the Triangle , 2000, SIAM J. Numer. Anal..

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

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

[25]  T. A. Zang,et al.  Spectral methods for fluid dynamics , 1987 .

[26]  Ivo Babuška,et al.  The optimal symmetrical points for polynomial interpolation of real functions in the tetrahedron , 1995 .

[27]  G. Karniadakis,et al.  Spectral/hp Element Methods for CFD , 1999 .

[28]  D. Kopriva.  A Conservative Staggered-Grid Chebyshev Multidomain Method for Compressible Flows. II. A Semi-Structured Method , 1996 .

[29]  Stephen Wolfram,et al.  The Mathematica book (4th edition) , 1999 .

[30]  Ivo Babuška,et al.  Approximate optimal points for polynomial interpolation of real functions in an interval and in a triangle , 1995 .

[31]  Marcel Vinokur,et al.  Conservation equations of gasdynamics in curvilinear coordinate systems , 1974 .

[32]  S. Orszag.  Numerical Simulation of Incompressible Flows Within Simple Boundaries. I. Galerkin (Spectral) Representations , 1971 .

[33]  Marcel Vinokur,et al.  Exact Integrations of Polynomials and Symmetric Quadrature Formulas over Arbitrary Polyhedral Grids , 1998 .

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

[35]  Stephen Wolfram,et al.  The Mathematica Book , 1996 .

[36]  J. Benek,et al.  A flexible grid embedding technique with application to the Euler equations , 1983 .

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

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

[39]  K. C. Chung.  A GENERALIZED FINITE-DIFFERENCE METHOD FOR HEAT TRANSFER PROBLEMS OF IRREGULAR GEOMETRIES , 1981 .

[40]  Narayanaswamy Balakrishnan,et al.  An upwind finite difference scheme for meshless solvers , 2003 .

[41]  W. H. Reed,et al.  Triangular mesh methods for the neutron transport equation , 1973 .

[42]  John T. Batina,et al.  A gridless Euler/Navier-Stokes solution algorithm for complex-aircraft applications , 1993 .

[43]  Anthony T. Patera,et al.  An isoparametric spectral element method for solution of the Navier-Stokes equations in complex geometry , 1986 .