We present an overview of our work on Hybridized Discontinuous Galerkin (HDG) method for a variety of steady and time dependent conservation laws. The essential ingredients are a local Galerkin projection of the underlying PDEs at the element level onto spaces of polynomials of degree k to parametrize the numerical solution in terms of the numerical trace; a judicious choice of the numerical flux to provide stability and consistency; and a global jump condition that enforces the continuity of the numerical flux to arrive at a global weak formulation in terms of the numerical trace. We will also present the applications of the HDG method to Plasmonics.