No abstract

We demonstrate the potential of collocation methods for efficient higher-order analysis on standard nodal finite element meshes. We focus on a collocation method that is variationally consistent and geometrically flexible, converges optimally, embraces concepts of reduced quadrature, and leads to symmetric stiffness and diagonal consistent mass matrices. At the same time, it minimizes the evaluation cost per quadrature point, thus reducing formation and assembly effort significantly with respect to standard Galerkin finite element methods. We provide a detailed review of all components of the technology in the context of elastodynamics, that is, weighted residual formulation, nodal basis functions on Gauss-Lobatto quadrature points, and symmetrization by averaging with the ultra-weak formulation. We quantify potential gains by comparing the computational efficiency of collocated and standard finite elements in terms of basic operation counts and timings. Our results show that collocation is significantly less expensive for problems dominated by the formation and assembly effort, such as higher-order elastostatic analysis. Furthermore, we illustrate the potential of collocation for efficient higher-order explicit dynamics. Throughout this work, we advocate a straightforward implementation based on simple modifications of standard finite element codes. We also point out the close connection to spectral element methods, where many of the key ideas are already established. AbstractWe demonstrate the potential of collocation methods for efficient higher-order analysis on standard nodal finite element meshes. We focus on a collocation method that is variationally consistent and geometrically flexible, converges optimally, embraces concepts of reduced quadrature, and leads to symmetric stiffness and diagonal consistent mass matrices. At the same time, it minimizes the evaluation cost per quadrature point, thus reducing formation and assembly effort significantly with respect to standard Galerkin finite element methods. We provide a detailed review of all components of the technology in the context of elastodynamics, i.e. weighted residual formulation, nodal basis functions on Gauss-Lobatto quadrature points, and symmetrization by averaging with the ultra-weak formulation. We quantify potential gains by comparing the computational efficiency of collocation and standard finite elements in terms of basic operation counts and timings. Our results show that collocation is significantly less expensive for problems dominated by the formation and assembly effort, such as higher-order elastostatic analysis. Furthermore we illustrate the potential of collocation for efficient higher-order explicit dynamics. Throughout this work, we advocate a straightforward implementation based on simple modifications of standard finite element codes. We also point out the close connection to spectral element methods, where many of the key ideas are already established.

Recently, a multi-level hp-version of the finite element method (FEM) was proposed to ease the difficulties of treating hanging nodes, while providing full hp-approximation capabilities. In the original paper, the refinement procedure made use of a-priori knowledge of the solution. However, adaptive procedures can produce discretizations which are more effective than an intuitive choice of element sizes h and polynomial degree distributions p. This is particularly prominent when a-priori knowledge of the solution is only vague or unavailable. The present contribution demonstrates that multi-level hp-adaptive schemes can be efficiently driven by an explicit a-posteriori error estimator. To this end, we adopt the classical residual-based error estimator. The main insight here is that its extension to multi-level hp-FEM is possible by considering the refined-most overlay elements as integration domains. We demonstrate on several two-and three-dimensional examples that exponential convergence rates can be obtained.

No abstract

No abstract