Search

We propose a new hybridizable discontinuous Galerkin (HDG) model order reduction technique based on proper orthogonal decomposition (POD). We consider the heat equation as a test problem and prove error bounds that converge to zero as the number of POD modes increases. We present 2D and 3D numerical results to illustrate the convergence analysis.

We begin an investigation of hybridizable discontinuous Galerkin (HDG) methods for approximating the solution of Dirichlet boundary control problems governed by elliptic PDEs. These problems can involve atypical variational formulations, and often have solutions with low regularity on polyhedral domains. These issues can provide challenges for numerical methods and the associated numerical analysis. We propose an HDG method for a Dirichlet boundary control problem for the Poisson equation, and obtain optimal a priori error estimates for the control. Specifically, under certain assumptions, for a 2D convex polygonal domain we show the control converges at a superlinear rate. We present 2D and 3D numerical experiments to demonstrate our theoretical results.

This paper presents a new hybridizable discontinuous Galerkin (HDG) method for linear elasticity on general polyhedral meshes, based on a strong symmetric stress formulation. The key feature of this new HDG method is the use of a special form of the numerical trace of the stresses, which makes the error analysis different from the projection-based error analyzes used for most other HDG methods. For arbitrary polyhedral elements, we approximate the stress by using polynomials of degree k>=1 and the displacement by using polynomials of degree k+1. In contrast, to approximate the numerical trace of the displacement on the faces, we use polynomials of degree k only. This allows for a very efficient implementation of the method, since the numerical trace of the displacement is the only globally-coupled unknown, but does not degrade the convergence properties of the method. Indeed, we prove optimal orders of convergence for both the stresses and displacements on the elements. In the almost incompressible case, we show the error of the stress is also optimal in the standard L2-norm. These optimal results are possible thanks to a special superconvergence property of the numerical traces of the displacement, and thanks to the use of a crucial elementwise Korn's inequality. Several numerical results are presented to support our theoretical findings in the end.