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

eXtended hybridizable discontinuous Galerkin for incompressible flow problems with unfitted meshes and interfaces

The eXtended hybridizable discontinuous Galerkin (X-HDG) method is developed for the solution of Stokes problems with void or material interfaces. X-HDG is a novel method that combines the hybridizable discontinuous Galerkin (HDG) method with an eXtended finite element strategy, resulting in a high-order, unfitted, superconvergent method, with an explicit definition of the interface geometry by means of a level-set function. For elements not cut by the interface, the standard HDG formulation is applied, whereas a modified weak form for the local problem is proposed for cut elements. Heaviside enrichment is considered on cut faces and in cut elements in the case of bimaterial problems. Two-dimensional numerical examples demonstrate that the applicability, accuracy, and superconvergence properties of HDG are inherited in X-HDG, with the freedom of computational meshes that do not fit the interfacesPeer 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

HDG-NEFEM for two dimensional linear elasticity

This paper proposes a new methodology for the solution of two dimensional linear elastic problems in domains with curved boundaries. The proposed method exploits the advantages of the hybridisable discontinuous Galerkin method to obtain an accurate approximation of both the displacement and the stress fields by solving a global problem that only involves the displacement field on the element edges as unknown. In addition, the methodology incorporates the exact boundary representation of the domain by means of the so-called NURBS-enhanced finite element method. Numerical examples are used to illustrate the three main advantages of the proposed method, namely the reproducibility of polynomials in domains with curved boundaries, the super-convergence of the solution even for linear approximation and the effectiveness and reliability of degree adaptive processes driven by displacement or stresses

Adaptive low and high-order hybridized methods for unsteady incompressible flow simulations

Tesi en modalitat de cotutela: Universitat Politècnica de Catalunya i Università degli Studi di PaviaSimulations of incompressible flows are performed on a daily basis to solve problems of practical and industrial interest in several fields of engineering, including automotive, aeronautical, mechanical and biomedical applications. Although finite volume (FV) methods are still the preferred choice by the industry due to their efficiency and robustness, sensitivity to mesh quality and limited accuracy represent two main bottlenecks of these approaches. This is especially critical in the context of transient phenomena, in which FV methods show excessive numerical diffusion. In this context, there has been a growing interest towards high-order discretisation strategies in last decades. In this PhD thesis, a high-order adaptive hybidisable discontinuous Galerkin (HDG) method is proposed for the approximation of steady and unsteady laminar incompressible Navier-Stokes equations. Voigt notation for symmetric second-order tensors is exploited to devise an HDG method for the Cauchy formulation of the momentum equation with optimal convergence properties, even when low-order polynomial degrees of approximation are considered. In addition, a postprocessing strategy accounting for rigid translational and rotational modes is proposed to construct an element-by-element superconvergent velocity field. The discrepancy between the computed and postprocessed velocities is utilised to define a local error indicator to drive degree adaptivity procedures and accurately capture localised features of the flow. The resulting HDG solver is thus extended to the case of transient problems via high-order time integration schemes, namely the explicit singly diagonal implicit Runge-Kutta (ESDIRK) schemes. In this context, the embedded explicit step is exploited to define an inexpensive estimate of the temporal error to devise an efficient timestep control strategy. Finally, in order to efficiently solve the global problem arising from the HDG discretisation, a preconditioned iterative solver is proposed. This is critical in the context of high-order approximations in three-dimensional domains leading to large-scale problems, especially in transient simulations. A block diagonal preconditioner coupled with an inexpensive approximation of the Schur complement of the matrix is proposed to reduce the computational cost of the overall HDG solver. Extensive numerical validation of two and three-dimensional steady and unsteady benchmark tests of viscous laminar incompressible flows is performed to validate the proposed methodology.Simulaciones de flujo incompresible se emplean a diario para resolver problemas de interés práctico e industrial en varios campos de la ingeniería, p.ej. en aplicaciones automovilísticas, aeronáuticas, mecánicas y biomédicas. Aunque los métodos de volúmenes finitos (FV) siguen siendo la opción preferida por la industria debido a su eficiencia y robustez, la sensibilidad a la calidad de la malla y la baja precisión representan dos limitaciones importantes para estas técnicas. Estas limitaciones son todavía más críticas en el contexto de simulaciones de fenómenos transitorios, donde los FV están penalizados por su excesiva difusión numérica. En este contexto, las estrategias de discretización de alto orden han ganado una popularidad creciente en las últimas décadas para problemas transitorios dónde se necesitan soluciones precisas. Esta tesis propone un método de Galerkin discontinuo híbrido (HDG), de alto orden y adaptativo para la aproximación de las ecuaciones de Navier-Stokes incomprensible laminar, en el caso estacionario y transitorio en el entorno de aplicaciones ingenieriles. Para ello, la notación de Voigt para tensores simétricos de segundo orden (habituales en mecánica de los medios continuos) permite introducir un método HDG para la formulación de Cauchy de la ecuación de momento. La novedad de este resultado reside en la convergencia óptima alcanzada por el método, incluso para aproximaciones de orden polinómico bajo. Además, se desarrolla una estrategia de post-proceso local para construir elemento a elemento un campo de velocidad súper-convergente, tomando en cuenta los modos rígidos de traslación y rotación. La discrepancia entre el campo de velocidad calculado y el súper-convergente, obtenido a través del post-proceso, permite definir un indicador del error local. De esta forma, se desarrolla una estrategia para realizar adecuar elemento a elemento el grado de la aproximación polinómica y así mejorar la precisión adaptándose a las características localizadas del flujo. Seguidamente, se extiende el método HDG propuesto al tratamiento de problemas dependientes del tiempo. Más concretamente, se consideran los esquemas de integración temporal de alto orden explicit singly diagonal implicit Runge-Kutta (ESDIRK). En este contexto, se utiliza el paso explícito embedded para calcular una estimación computacionalmente eficiente del error temporal y definir una estrategia de adaptividad del paso de tiempo. Finalmente, se desarrolla un precondicionador adaptado a la estrategia HDG que acelera la convergencia del método iterativo empleado y, de esta forma, obtener resoluciones eficaces del problema global surgido de la discretización HDG. Es importante resaltar la importancia de una herramienta de resolución eficiente para problemas de gran escala en el contexto de aproximaciones de alto orden y en dominios tridimensionales. Estas herramientas se hacen aún más criticas en simulaciones transitorias. Más concretamente, se proponen un precondicionador diagonal por bloques y una aproximación eficiente del complemento Schur de la matriz para reducir el coste computacional del método HDG. Para validar la metodología propuesta, se realizan varias simulaciones numéricas de flujo incompresible laminar viscoso, para problemas estacionarios y transitorios, en dos y tres dimensiones.Postprint (published version

A superconvergent HDG method for stokes flow with strongly enforced symmetry of the stress tensor

This work proposes a superconvergent hybridizable discontinuous Galerkin (HDG) method for the approximation of the Cauchy formulation of the Stokes equation using same degree of polynomials for the primal and mixed variables. The novel formulation relies on the well-known Voigt notation to strongly enforce the symmetry of the stress tensor. The proposed strategy introduces several advantages with respect to the existing HDG formulations. First, it remedies the suboptimal behavior experienced by the classical HDG method for formulations involving the symmetric part of the gradient of the primal variable. The optimal convergence of the mixed variable is retrieved and an element-by-element postprocess procedure leads to a superconvergent velocity field, even for low-order approximations. Second, no additional enrichment of the discrete spaces is required and a gain in computational efficiency follows from reducing the quantity of stored information and the size of the local problems. Eventually, the novel formulation naturally imposes physical tractions on the Neumann boundary. Numerical validation of the optimality of the method and its superconvergent properties is performed in 2D and 3D using meshes of different element types

A superconvergent HDG method for stokes flow with strongly enforced symmetry of the stress tensor

The final publication is available at Springer via http://dx.doi.org/10.1007/s10915-018-0855-yThis work proposes a superconvergent hybridizable discontinuous Galerkin (HDG) method for the approximation of the Cauchy formulation of the Stokes equation using same degree of polynomials for the primal and mixed variables. The novel formulation relies on the well-known Voigt notation to strongly enforce the symmetry of the stress tensor. The proposed strategy introduces several advantages with respect to the existing HDG formulations. First, it remedies the suboptimal behavior experienced by the classical HDG method for formulations involving the symmetric part of the gradient of the primal variable. The optimal convergence of the mixed variable is retrieved and an element-by-element postprocess procedure leads to a superconvergent velocity field, even for low-order approximations. Second, no additional enrichment of the discrete spaces is required and a gain in computational efficiency follows from reducing the quantity of stored information and the size of the local problems. Eventually, the novel formulation naturally imposes physical tractions on the Neumann boundary. Numerical validation of the optimality of the method and its superconvergent properties is performed in 2D and 3D using meshes of different element types.Peer ReviewedPostprint (author's final draft

A superconvergent hybridisable discontinuous Galerkin method for linear elasticity

The first superconvergent hybridisable discontinuous Galerkin method for linear elastic problems capable of using the same degree of approximation for both the primal and mixed variables is presented. The key feature of the method is the strong imposition of the symmetry of the stress tensor by means of the well known and extensively used Voigt notation, circumventing the use of complex mathematical concepts to enforce the symmetry of the stress tensor either weakly or strongly. A novel procedure to construct element by element a superconvergent postprocessed displacement is proposed. Contrary to other hybridisable discontinuous Galerkin formulations, the methodology proposed here is able to produce a superconvergent displacement field for low-order approximations. The resulting method is robust and locking-free in the nearly incompressible limit. An extensive set of numerical examples is utilised to provide evidence of the optimality of the method and its superconvergent properties in two and three dimensions and for different element type

Tutorial on Hybridizable Discontinous Galerkin (HDG) for second-order elliptic problems

The HDG is a new class of discontinuous Galerkin (DG) methods that shares favorable properties with classical mixed methods such as the well known Raviart-Thomas methods. In particular, HDG provides optimal convergence of both the primal and the dual variables of the mixed formulation. This property enables the construction of superconvergent solutions, contrary to other popular DG methods. In addition, its reduced computational cost, compared to other DG methods, has made HDG an attractive alternative for solving problems governed by partial differential equations. A tutorial on HDG for the numerical solution of second-order elliptic problems is presented. Particular emphasis is placed on providing all the necessary details for the implementation of HDG methods.Peer ReviewedPreprin