Abstract:We present an overarching framework for stable spectral methods on a triangle, defined by a multivariate W-system and based on orthogonal polynomials on the triangle. Motivated by the Koornwinder orthogonal polynomials on the triangle, we introduce a Koornwinder W-system. Once discretised by this W-system, the resulting spatial differentiation matrix is skew symmetric, affording important advantages insofar as stability and conservation of structure are concerned. We analyse the construction of the differentiation matrix and matrix vector multiplication, demonstrating optimal computational cost. Numerical convergence is illustrated through experiments with different parameter choices. As a result, our method exhibits key characteristics of a practical spectral method, inclusive of rapid convergence, fast computation and the preservation of structure of the underlying partial differential equation.
What problem does this paper attempt to address?
This paper attempts to solve the problem of stably applying spectral methods to solve time - dependent partial differential equations (PDEs) on triangular regions. Specifically, the article proposes a framework based on multivariate W - systems, using orthogonal polynomials on triangles to construct spectral methods. By introducing the Koornwinder W - system, the discretized spatial differential matrix exhibits skew - symmetric characteristics, which provides important advantages for stability and structure preservation.
The following is a detailed interpretation of the main problems the paper attempts to solve:
1. **Constructing a stable spectral method**:
- The paper aims to develop a stable spectral method for time - dependent partial differential equations applicable to triangular regions. The application of traditional spectral methods on complex geometries (such as triangles) faces challenges, especially in ensuring numerical stability and structure preservation.
- To this end, the author introduces the Koornwinder W - system and proves that its corresponding differential matrix is skew - symmetric, which is crucial for stability.
2. **Improving computational efficiency**:
- The article explores the construction of differential matrices and matrix - vector multiplications, demonstrating the optimal computational cost. By introducing skew - symmetric differential matrices, the computational complexity can be significantly reduced and the algorithm efficiency can be improved.
- Numerical experiments show that this method has the key characteristics of fast convergence, fast computation, and preservation of the underlying partial differential equation structure.
3. **Handling different boundary conditions**:
- The paper not only considers zero Dirichlet boundary conditions but also discusses the treatment methods for general Dirichlet boundary conditions. By introducing ideas from computer - aided geometric design, the author proposes a method for handling arbitrary Dirichlet boundary conditions.
- This improvement enables the W - system to be applied to a wider range of boundary conditions, thus expanding its application scope in practical problems.
4. **Theoretical analysis and numerical verification**:
- The author conducts a detailed theoretical analysis, derives the specific expressions of differential matrices, and verifies the effectiveness of the method through numerical experiments.
- The experimental results show that the method exhibits good convergence and computational efficiency under different parameter selections.
In summary, this paper is committed to constructing an efficient and stable spectral method on triangular regions, solving the problems of stability, computational efficiency, and boundary condition handling faced by traditional spectral methods when applied to complex geometries. By introducing the Koornwinder W - system, the author successfully achieves these goals and provides new ideas and tools for future research.