text: The finite element methods are important and powerful tools for solving partial differential equations on complex geometries, but their high-order generalizations pose significant challenges in terms of robustness and mesh generation. In this paper, we introduce a high-order finite element method, which is insensitive to element quality and hence is more robust, and which requires only linear elements even for curved geometries and hence significantly simplifies meshing. Based on our recent work on Adaptive Extended Stencil FEM, or AES-FEM (Int. J. Num. Meth. Engrg., 2016, DOI:10.1002 nme.5246), the proposed method replaces the traditional interpolatory FEM basis functions with generalized Lagrange polynomial basis functions, constructed via local weighted least-squares approximations. In this work, we show that AES-FEM can achieve high-order accuracy using linear elements to discretize the geometry and using the standard FEM hat functions as test functions even for curved geometries. This is in contrast to most high-order generalizations of FEM, which require high-order curved elements for curved geometries. We describe the selection of stencils and the treatment of boundary conditions for high-order AES-FEM, and show that even-degree polynomial basis functions are preferred for even-order PDEs. We present numerical results in 2D and 3D for the Poisson equation and the convection-diffusion equation and demonstrate up to sixth order accuracy with AES-FEM. Overall, AES-FEM is much easier to use than FEM for curved boundaries, while being more accurate, more stable, and more efficient in terms of runtime versus error over relatively fine meshes.
  @cite_26: The implementation of @math hp-adaptivity is challenging as hanging nodes, edges, and faces have to be constrained to ensure compatibility of the shape functions. For this reason, most @math hp-code frameworks restrict themselves to @math 1-irregular meshes to ease the implementational effort. This work alleviates these difficulties by introducing a new formulation for high-order mesh adaptivity that provides full local @math hp-refinement capabilities at a comparably small implementational effort. Its main idea is the extension of the @math hp- @math d-method such that it allows for high-order overlay meshes yielding a hierarchical, multi-level @math hp-formulation of the Finite Element Method. This concept enables intuitive refinement and coarsening procedures, while linear independence and compatibility of the shape functions are guaranteed by construction. The proposed method is demonstrated to achieve exponential rates of convergence--both in terms of degrees of freedom and in run-time--for problems with non-smooth solutions. Furthermore, the scheme is used alongside the Finite Cell Method to simulate the heat flow around moving objects on a non-conforming background mesh and is combined with an energy-based refinement indicator for automatic @math hp-adaptivity.
  @cite_33: 
  @cite_24: The Finite Element Method: Its Basis and Fundamentals offers a complete introduction to the basis of the finite element method, covering fundamental theory and worked examples in the detail required for readers to apply the knowledge to their own engineering problems and understand more advanced applications. This edition sees a significant rearrangement of the book's content to enable clearer development of the finite element method, with major new chapters and sections added to cover: * Weak forms * Variational forms * Multi-dimensional field problems * Automatic mesh generation * Plate bending and shells * Developments in meshless techniques Focusing on the core knowledge, mathematical and analytical tools needed for successful application, The Finite Element Method: Its Basis and Fundamentals is the authoritative resource of choice for graduate level students, researchers and professional engineers involved in finite element-based engineering analysis. * A proven keystone reference in the library of any engineer needing to understand and apply the finite element method in design and development. * Founded by an influential pioneer in the field and updated in this seventh edition by an author team incorporating academic authority and industrial simulation experience. * Features reworked and reordered contents for clearer development of the theory, plus new chapters and sections on mesh generation, plate bending, shells, weak forms and variational forms.
  @cite_40: In this paper we present a new automatic adaptivity algorithm for the hp-FEM which is based on arbitrary-level hanging nodes and local element projections. The method is very simple to implement compared to other existing hp-adaptive strategies, while its performance is comparable or superior. This is demonstrated on several numerical examples which include the L-shape domain problem, a problem with internal layer, and the Girkmann problem of linear elasticity. With appropriate simplifications, the proposed technique can be applied to standard lower-order and spectral finite element methods.
  @cite_2: We present an algorithm, and a 2D implementation for a fully automatic hp-adaptive strategy for elliptic problems. Given a mesh, the next, optimally refined mesh, is determined by maximizing the rate of decrease of the hp-interpolation error for a reference solution. Numerical results confirm optimal, exponential convergence rates predicted by the theory of hp methods.
  @cite_16: 
summary: The isoparametric FEM in general uses the same-degree polynomial basis uniformly. The @math -adaptive FEM generalizes it by allowing the degrees of the basis functions to vary from element to element @cite_33 @cite_24 . The @math -adaptivity may be further combined with @math -adaptivity, leading to the @math -adaptive FEM or @math -FEM @cite_16 @cite_26 . In these methods, inter-element continuity is typically enforced as constraints @cite_40 . The @math -FEM can achieve high-order convergence, with error decreasing exponentially in the degree of basis functions, provided that the mesh is nearly optimal @cite_2 . Similar to isoparametric FEM, the @math -FEM also requires curved elements for curved geometries, requires good element shapes for accuracy, and is limited to only moderately high orders due to the use of equally spaced nodes.