Abstract:To solve hyperbolic conservation laws, a new method is developed based on the spectral difference (SD) algorithm. The new scheme adopts hierarchical polynomials to represent the solution in each cell instead of Lagrange interpolation polynomials used by the original one. The degrees of freedom (DOFs) of the present scheme are the coefficients of these polynomials, which do not represent the states at the solution points like the original method. Therefore, the solution points defined in the original SD scheme are discarded, while the flux points are preserved to construct a Lagrange interpolation polynomial to approximate flux function in each cell. To update the DOFs, differential operators are applied to the governing equation as well as the Lagrange interpolation polynomial of flux function to evaluate first and higher order derivatives of both solution and flux at the centroid of the cell. The stability property of the current scheme is proved to be the same as the original SD method when the same solution space is adopted. One dimensional methods are always stable by the use of zeros of Legendre polynomials as inner flux points. For two dimensional problems, the introduction of Raviart–Thomas spaces for the interpolation of flux function proves stable schemes for triangles. Accuracy studies are performed with one- and two-dimensional problems. p-Multigrid algorithm is implemented with orthogonal hierarchical bases. The results verify the high efficiency and low memory requirements of implementation of p-multigrid algorithm with the proposed scheme.

The Development of Numerical Computing by Spectral Method

A New Spectral Difference Method Using Hierarchical Polynomial Bases for Hyperbolic Conservation Laws

Numerical evolutions of fields on the 2-sphere using a spectral method based on spin-weighted spherical harmonics

The ultraspherical spectral element method

A review of the application of spectral methods in computational ocean acoustics

On spectral numerical method for variable-order partial differential equations

A numerical method for solving elliptic equations on real closed algebraic curves and surfaces

Understanding the ultraspherical spectral method

Solving time-dependent PDEs with the ultraspherical spectral method

Spectral/hp element methods: recent developments, applications, and perspectives

A new spectral element method for numerical solution of partial differential equations on annular-type regions

A Review: Applications of the Spectral Finite Element Method

The spectral element method for the solution of Maxwell’s equations

Numerical Study of an Adaptive Domain Decomposition Algorithm Based on Chebyshev Tau Method for Solving Singular Perturbed Problems

Spectral solver for Cauchy problems in polar coordinates using discrete Hankel transforms

Efficient Spectral Methods for Quasi-Equilibrium Closure Approximations of Symmetric Problems on Unit Circle and Sphere

Application of spectral element method for solving Sobolev equations with error estimation

Benchmarking Numerical Algorithms for Harmonic Maps into the Sphere

A Spectral Method for Parabolic Differential Equations

A Domain Decomposition Based Spectral Collocation Method for Lane-Emden Equations