Discontinuous Galerkin approximations in computational mechanics: hybridization, exact geometry and degree adaptivity

Discontinuous Galerkin (DG) discretizations with exact representation of the geometry and local polynomial degree adaptivity are revisited. Hybridization techniques are employed to reduce the computational cost of DG approximations and devise the hybridizable discontinuous Galerkin (HDG) method. Exact geometry described by non-uniform rational B-splines (NURBS) is integrated into HDG using the framework of the NURBS-enhanced finite element method (NEFEM). Moreover, optimal convergence and superconvergence properties of HDG-Voigt formulation in presence of symmetric second-order tensors are exploited to construct inexpensive error indicators and drive degree adaptive procedures. Applications involving the numerical simulation of problems in electrostatics, linear elasticity and incompressible viscous flows are presented. Moreover, this is done for both high-order HDG approximations and the lowest-order framework of face-centered finite volumes (FCFV).Peer ReviewedPostprint (author's final draft

HDGlab: An Open-Source Implementation of the Hybridisable Discontinuous Galerkin Method in MATLAB

This paper presents HDGlab, an open source MATLAB implementation of the hybridisable discontinuous Galerkin (HDG) method. The main goal is to provide a detailed description of both the HDG method for elliptic problems and its implementation available in HDGlab. Ultimately, this is expected to make this relatively new advanced discretisation method more accessible to the computational engineering community. HDGlab presents some features not available in other implementations of the HDG method that can be found in the free domain. First, it implements high-order polynomial shape functions up to degree nine, with both equally-spaced and Fekete nodal distributions. Second, it supports curved isoparametric simplicial elements in two and three dimensions. Third, it supports non-uniform degree polynomial approximations and it provides a flexible structure to devise degree adaptivity strategies. Finally, an interface with the open-source high-order mesh generator Gmsh is provided to facilitate its application to practical engineering problems

An explicit hybridizable discontinuous Galerkin method for the 3D time-domain Maxwell equations

International audienceWe present an explicit hybridizable discontinuous Galerkin (HDG) method for numerically solving the system of three-dimensional (3D) time-domain Maxwell equations. The method is fully explicit similarly to classical so-called DGTD (Dis-continuous Galerkin Time-Domain) methods, is also high-order accurate in both space and time and can be seen as a generalization of the classical DGTD scheme based on upwind fluxes. We provide numerical results aiming at assessing its numerical convergence properties by considering a model problem and we present preliminary results of the superconvergence property on the H curl norm

A Class of Embedded DG Methods for Dirichlet Boundary Control of Convection Diffusion PDEs

We investigated an hybridizable discontinuous Galerkin (HDG) method for a convection diffusion Dirichlet boundary control problem in our earlier work [SIAM J. Numer. Anal. 56 (2018) 2262-2287] and obtained an optimal convergence rate for the control under some assumptions on the desired state and the domain. In this work, we obtain the same convergence rate for the control using a class of embedded DG methods proposed by Nguyen, Peraire and Cockburn [J. Comput. Phys. vol. 302 (2015), pp. 674-692] for simulating fluid flows. Since the global system for embedded DG methods uses continuous elements, the number of degrees of freedom for the embedded DG methods are smaller than the HDG method, which uses discontinuous elements for the global system. Moreover, we introduce a new simpler numerical analysis technique to handle low regularity solutions of the boundary control problem. We present some numerical experiments to confirm our theoretical results

Adaptive hybrid discontinuous methods for fluid and wave problems

This PhD thesis proposes a p-adaptive technique for the Hybridizable Discontinuous Galerkin method (HDG). The HDG method is a novel discontinuous Galerkin method (DG) with interesting characteristics. While retaining all the advantages of the common DG methods, such as the inherent stabilization and the local conservation properties, HDG allows to reduce the coupled degrees of freedom of the problem to those of an approximation of the solution de¿ned only on the faces of the mesh. Moreover, the convergence properties of the HDG solution allow to perform an element-by-element postprocess resulting in a superconvergent solution. Due to the discontinuous character of the approximation in HDG, p-variable computations are easily implemented. In this work the superconvergent postprocess is used to de¿ne a reliable and computationally cheap error estimator, that is used to drive an automatic adaptive process. The polynomial degree in each element is automatically adjusted aiming at obtaining a uniform error distribution below a user de¿ned tolerance. Since no topological modi¿cation of the discretization is involved, fast adaptations of the mesh are obtained. First, the p-adaptive HDG is applied to the solution of wave problems. In particular, the Mild Slope equation is used to model the problem of sea wave propagation is coastal areas and harbors. The HDG method is compared with the continuous Galerkin (CG) ¿nite element method, which is nowadays the common method used in the engineering practice for this kind of applications. Numerical experiments reveal that the e¿ciency of HDG is close to CG for uniform degree computations, clearly outperforming other DG methods such as the Compact Discontinuous Galerkin method. When p-adaptivity is considered, an important saving in computational cost is shown. Then, the methodology is applied to the solution of the incompressible Navier-Stokes equations for the simulation of laminar ¿ows. Both steady state and transient applications are considered. Various numerical experiments are presented, in 2D and 3D, including academic examples and more challenging applications of engineering interest. Despite the simplicity and low cost of the error estimator, high e¿ciency is exhibited for analytical examples. Moreover, even though the adaptive technique is based on an error estimate for just the velocity ¿eld, high accuracy is attained for all variables, with sharp resolution of the key features of the ¿ow and accurate evaluation of the ¿uid-dynamic forces. In particular, high degrees are automatically located along boundary layers, reducing the need for highly distorted elements in the computational mesh. Numerical tests show an important reduction in computational cost, compared to uniform degree computations, for both steady and unsteady computations

An explicit hybridizable discontinuous Galerkin method for the 3D time-domain Maxwell equations

International audienceWe present an explicit hybridizable discontinuous Galerkin (HDG) method for numerically solving the system of three-dimensional (3D) time-domain Maxwell equations. The method is fully explicit similarly to classical so-called DGTD (Dis-continuous Galerkin Time-Domain) methods, is also high-order accurate in both space and time and can be seen as a generalization of the classical DGTD scheme based on upwind fluxes. We provide numerical results aiming at assessing its numerical convergence properties by considering a model problem and we present preliminary results of the superconvergence property on the H curl norm