A Least Squares Finite Element Method (LSFEM) formulation for the detection of unknown boundary conditions in problems governed by the steady convection-diffusion equation will be presented. The method is capable of determining temperatures, and heat fluxes in location where such quantities are unknown provided such quantities are sufficiently over-specified in other locations. For the current formulation it is assumed the velocity field is known. The current formulation is unique in that it results in a sparse square system of equations even for partial differential equations that are not self-adjoint. Since this formulation always results in a symmetric positive-definite matrix, the solution can be found with standard sparse matrix solvers such as preconditioned conjugate gradient method. In addition, the formulation allows for equal order approximation of temperature and heat fluxes as it is not subject to the inf-sup condition. The formulation allow for a treatment of over-specified boundary conditions. Also, various forms of regularization can be naturally introduced within the formulation. Details of the discretization and sample results will be presented.