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

Adjoint-Based Error Estimation and Mesh Adaptation for Hybridized Discontinuous Galerkin Methods

We present a robust and efficient target-based mesh adaptation methodology, building on hybridized discontinuous Galerkin schemes for (nonlinear) convection-diffusion problems, including the compressible Euler and Navier-Stokes equations. Hybridization of finite element discretizations has the main advantage, that the resulting set of algebraic equations has globally coupled degrees of freedom only on the skeleton of the computational mesh. Consequently, solving for these degrees of freedom involves the solution of a potentially much smaller system. This not only reduces storage requirements, but also allows for a faster solution with iterative solvers. The mesh adaptation is driven by an error estimate obtained via a discrete adjoint approach. Furthermore, the computed target functional can be corrected with this error estimate to obtain an even more accurate value. The aim of this paper is twofold: Firstly, to show the superiority of adjoint-based mesh adaptation over uniform and residual-based mesh refinement, and secondly to investigate the efficiency of the global error estimate

A temporally adaptive hybridized discontinuous Galerkin method for time-dependent compressible flows

The potential of the hybridized discontinuous Galerkin (HDG) method has been recognized for the computation of stationary flows. Extending the method to time-dependent problems can, e.g., be done by backward difference formulae (BDF) or diagonally implicit Runge-Kutta (DIRK) methods. In this work, we investigate the use of embedded DIRK methods in an HDG solver, including the use of adaptive time-step control. Numerical results demonstrate the performance of the method for both linear and nonlinear (systems of) time-dependent convection-diffusion equations

Monolithic multiphysics simulation of hypersonic aerothermoelasticity using a hybridized discontinuous Galerkin method

This work presents implementation of a hybridized discontinuous Galerkin (DG) method for robust simulation of the hypersonic aerothermoelastic multiphysics system. Simulation of hypersonic vehicles requires accurate resolution of complex multiphysics interactions including the effects of high-speed turbulent flow, extreme heating, and vehicle deformation due to considerable pressure loads and thermal stresses. However, the state-of-the-art procedures for hypersonic aerothermoelasticity are comprised of low-fidelity approaches and partitioned coupling schemes. These approaches preclude robust design and analysis of hypersonic vehicles for a number of reasons. First, low-fidelity approaches limit their application to simple geometries and lack the ability to capture small scale flow features (e.g. turbulence, shocks, and boundary layers) which greatly degrades modeling robustness and solution accuracy. Second, partitioned coupling approaches can introduce considerable temporal and spatial inaccuracies which are not trivially remedied. In light of these barriers, we propose development of a monolithically-coupled hybridized DG approach to enable robust design and analysis of hypersonic vehicles with arbitrary geometries. Monolithic coupling methods implement a coupled multiphysics system as a single, or monolithic, equation system to be resolved by a single simulation approach. Further, monolithic approaches are free from the physical inaccuracies and instabilities imposed by partitioned approaches and enable time-accurate evolution of the coupled physics system. In this work, a DG method is considered due to its ability to accurately resolve second-order partial differential equations (PDEs) of all classes. We note that the hypersonic aerothermoelastic system is composed of PDEs of all three classes. Hybridized DG methods are specifically considered due to their exceptional computational efficiency compared to traditional DG methods. It is expected that our monolithic hybridized DG implementation of the hypersonic aerothermoelastic system will 1) provide the physical accuracy necessary to capture complex physical features, 2) be free from any spatial and temporal inaccuracies or instabilities inherent to partitioned coupling procedures, 3) represent a transition to high-fidelity simulation methods for hypersonic aerothermoelasticity, and 4) enable efficient analysis of hypersonic aerothermoelastic effects on arbitrary geometries