A diagonal-mass-matrix triangular-spectral-element method based on cubature points
Taylor, M. A.
MetadataShow full item record
The cornerstone of nodal spectral-element methods is the co-location of the interpolation and integration points, yielding a diagonal mass matrix that is efficient for time-integration. On quadrilateral elements, Legendre-Gauss-Lobatto points are both good interpolation and integration points; on triangles, analogous points have not yet been found. In this paper, a promising set of points for the triangle that were only available for polynomial degree N <= 5 are used. However, the derivation of these points is generalized to obtain degree N <= 7 points, which are referred to as cubature points because their selection is based on their integration accuracy. The diagonal-mass-matrix (DMM) triangular-spectral-element (TSE) method based on these points can be used for any set of equations and any type of domain. Because these cubature points integrate up to order 2N along the element boundaries and yield a diagonal mass matrix, they allow the triangular spectral elements to compete with quadrilateral spectral elements in terms of both accuracy and efficiency, while offering more geometric flexibility in the choice of grids. It is shown how to implement this DMM TSE for a variety of applications involving elliptic and hyperbolic equations on different domains. The DMM TSE method yields comparable accuracy to the exact integration (non-DMM) TSE method, while being far more efficient for time-dependent problems.
The article of record as published may be located at http://dx.doi.org/10.1007/s10665-006-9085-7
RightsThis publication is a work of the U.S. Government as defined in Title 17, United States Code, Section 101. Copyright protection is not available for this work in the United States.
Showing items related by title, author, creator and subject.
Giraldo, F.X.; Warburton, T. (2005);A nodal triangle-based spectral element (SE) method for the shallow water equations on the sphere is presented. The original SE method uses quadrilateral elements and high-order nodal Lagrange polynomials, constructed from ...
Efficient Construction of Unified Continuous and Discontinuous Galerkin Formulations for the 3D Euler Equations Abdi, Daniel S.; Giraldo, Francis X. (2016-05-28);A unified approach for the numerical solution of the 3D hyperbolic Euler equations using high order methods, namely continuous Galerkin (CG) and discontinuous Galerkin (DG) methods, is presented. First, we examine how ...
Yan, Hui; Ross, I. Michael; Fahroo, Fariba (The American Institute of Aeronautics and Astronautics (AIAA), 2001-06-25);We develop state feedback control laws for linear time-varying systems with quadratic cost criteria by an indirect Legendre pseudospectral method. This method approximates the linear two-point boundary value problem to a ...