Two new convection-dominated are derived under the approximate solutions least-squares mixed finite element procedures are formulated for solving Sobolev equations. Optimal H（div;Ω）×H1（Ω） norms error estimates standard mixed finite spaces. Moreover, these two schemes provide the with first-order and second-order accuracy in time increment, respectively.
