We design and analyze the first hybridizable discontinuous Galerkin methods for stationary, third-order linear equations in one-space dimension. The methods are defined as discrete versions of characterizations of the exact solution in terms of local problems and transmission conditions. They provide approximations to the exact solution u and its derivatives q := u' and p := u'' which are piecewise polynomials of degree ku, kq and kp, respectively. We consider the methods for which the difference between these polynomial degrees is at most two. We prove that all these methods have superconvergence properties which allows us to prove that their numerical traces converge at the nodes of the partition with order at least 2 k+1, where k is the minimum of ku, kq, and kp. This allows us to use an element-by-element postprocessing to obtain new approximations for u, q and p converging with order at least 2k +1 uniformly. Numerical results validating our error estimates are displayed.
Bibliographical noteFunding Information:
The research of the second author was partially supported by the National Science Foundation (grant DMS-1115331). The research of the third author was partially supported by the National Science Foundation (grant DMS-1419029).