  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.
