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.
Another class of high-order FEM is the discontinuous Galerkin (DG) methods. Unlike in isoparametric or spectral element methods, the basis functions in DG are not continuous between elements. Instead, a numerical flux is used at the element boundaries. DG is most conveniently applied to hyperbolic PDEs, such as the Navier-Stokes and compressible Euler equations @cite_6 @cite_31 , while extensions to parabolic and elliptic problems have also been developed @cite_29 . DG can achieve high-order accuracy and can be used on complex geometries @cite_28 , but high-quality curved elements are required for curved boundaries. Other methods related to DG include the interior penalty methods for elliptical problems, which can be unified with DG under one framework @cite_25 , the compact DG @cite_35 , which has a more compact stencil than local DG methods, and hybridized DG @cite_9 , which has a lower number of degrees of freedom.
@cite_35 @cite_28 @cite_29 @cite_9 @cite_6 @cite_31 @cite_25
140909251 791179977 2100903519 1971445355
Discontinuous Galerkin (DG) methods for solving partial differential equations, developed in the late 1990s, have become popular among computational scientists. This book covers both theory and computation as it focuses on three primal DG methods--the symmetric interior penalty Galerkin, incomplete interior penalty Galerkin, and nonsymmetric interior penalty Galerkin which are variations of interior penalty methods. The author provides the basic tools for analysis and discusses coding issues, including data structure, construction of local matrices, and assembling of the global matrix. Computational examples and applications to important engineering problems are also included. Discontinuous Galerkin Methods for Solving Elliptic and Parabolic Equations: Theory and Implementation is divided into three parts: Part I focuses on the application of DG methods to second order elliptic problems in one dimension and in higher dimensions. Part II presents the time-dependent parabolic problems without and with convection. Part III contains applications of DG methods to solid mechanics (linear elasticity), fluid dynamics (Stokes and Navier Stokes), and porous media flow (two-phase and miscible displacement). Appendices contain proofs and MATLAB code for one-dimensional problems for elliptic equations and routines written in C that correspond to algorithms for the implementation of DG methods in two or three dimensions. Audience: This book is intended for numerical analysts, computational and applied mathematicians interested in numerical methods for partial differential equations or who study the applications discussed in the book, and engineers who work in fluid dynamics and solid mechanics and want to use DG methods for their numerical results. The book is appropriate for graduate courses in finite element methods, numerical methods for partial differential equations, numerical analysis, and scientific computing. Chapter 1 is suitable for a senior undergraduate class in scientific computing. Contents: List of Figures; List of Tables; List of Algorithms; Preface; Part I: Elliptic Problems; Chapter 1: One-dimensional problem; Chapter 2: Higher dimensional problem; Part II: Parabolic Problems; Chaper 3: Purely parabolic problems; Chapter 4: Parabolic problems with convection; Part III: Applications; Chapter 5: Linear elasticity; Chapter 6: Stokes flow; Chapter 7: Navier-Stokes flow; Chapter 8: Flow in porous media; Appendix A: Quadrature rules; Appendix B: DG codes; Appendix C: An approximation result; Bibliography; Index. We present a time-explicit discontinuous Galerkin (DG) solver for the time-domain acoustic wave equation on hybrid meshes containing vertex-mapped hexahedral, wedge, pyramidal and tetrahedral elements. Discretely energy-stable formulations are presented for both Gauss-Legendre and Gauss-Legendre-Lobatto (Spectral Element) nodal bases for the hexahedron. Stable timestep restrictions for hybrid meshes are derived by bounding the spectral radius of the DG operator using order-dependent constants in trace and Markov inequalities. Computational efficiency is achieved under a combination of element-specific kernels (including new quadrature-free operators for the pyramid), multi-rate timestepping, and acceleration using Graphics Processing Units. We present an energy-stable low-storage high order discontinuous Galerkin method for the wave equation.We describe algorithms for efficient implementation on Graphics Processing Units (GPUs).We derive order-dependent global and local timestep restrictions for DG on hybrid meshes.We verify high order convergence and computational efficiency with numerical experiments. We present a compact discontinuous Galerkin (CDG) method for an elliptic model problem. The problem is first cast as a system of first order equations by introducing the gradient of the primal unknown, or flux, as an additional variable. A standard discontinuous Galerkin (DG) method is then applied to the resulting system of equations. The numerical interelement fluxes are such that the equations for the additional variable can be eliminated at the element level, thus resulting in a global system that involves only the original unknown variable. The proposed method is closely related to the local discontinuous Galerkin (LDG) method [B. Cockburn and C.-W. Shu, SIAM J. Numer. Anal., 35 (1998), pp. 2440-2463], but, unlike the LDG method, the sparsity pattern of the CDG method involves only nearest neighbors. Also, unlike the LDG method, the CDG method works without stabilization for an arbitrary orientation of the element interfaces. The computation of the numerical interface fluxes for the CDG method is slightly more involved than for the LDG method, but this additional complication is clearly offset by increased compactness and flexibility. Compared to the BR2 [F. Bassi and S. Rebay, J. Comput. Phys., 131 (1997), pp. 267-279] and IP [J. Douglas, Jr., and T. Dupont, in Computing Methods in Applied Sciences (Second Internat. Sympos., Versailles, 1975), Lecture Notes in Phys. 58, Springer, Berlin, 1976, pp. 207-216] methods, which are known to be compact, the present method produces fewer nonzero elements in the matrix and is computationally more efficient. We present a high-order discontinuous Galerkin (DG) method for solving the time-dependent Maxwell equations on non-conforming hybrid meshes. The hybrid mesh combines unstructured tetrahedra for the discretization of irregularly shaped objects with a hexahedral mesh for the rest of the computational domain. The transition between tetrahedra and hexahedra is completely non-conform, that is, no pyramidal or prismatic elements are introduced to link these elements. Within each mesh element, the electromagnetic field components are approximated by a arbitrary order nodal polynomial and a centered approximation is used for the evaluation of numerical fluxes at inter-element boundaries. The time integration of the associated semi-discrete equations is achieved by a fourth-order leap-frog scheme. The method is described and discussed, including algorithm formulation, stability, and practical implementation issues such as the hybrid mesh generation and the computation of flux matrices with cubature rules. We illustrate the performance of the proposed method on several two- and three-dimensional examples involving comparisons with DG methods on single element-type meshes. The results show that the use of non-conforming hybrid meshes in DG methods allows for a notable reduction in computing time without sacrificing accuracy.