  We propose a hybridizable discontinuous Galerkin (HDG) method to approximate
the solution of a distributed optimal control problem governed by an elliptic
convection diffusion PDE. We derive optimal a priori error estimates for the
state, adjoint state, their fluxes, and the optimal control. We present 2D and
3D numerical experiments to illustrate our theoretical results.
