In this paper, a few dual least‐squares finite element methods and their application to scalar linear hyperbolic problems are studied. The purpose is to obtain L2‐norm approximations on finite element spaces of the exact solutions to hyperbolic partial differential equations of interest. This is approached by approximating the generally infeasible quadratic minimization that defines the L2‐orthogonal projection of the exact solution, by feasible least‐squares principles using the ideas of the original
method proposed in the context of elliptic equations. All methods in this paper are founded upon and extend the
approach that is rather general and applicable beyond the setting of elliptic problems. Error bounds are shown that point to the factors affecting the convergence and provide conditions that guarantee optimal rates. Furthermore, the preconditioning of the resulting linear systems is discussed. Numerical results are provided to illustrate the behavior of the methods on common finite element spaces.