“…Consequently, most finite element models of the Navier-Stokes equations based on the weak-form Galerkin procedure do not guarantee the minimization of the error in the solution or in the differential equation. Least-squares finite element models offers an appealing alternative to the commonly used weak-form Galerkin procedure for fluids and have received substantial attention in the academic literature in recent years (see, for example, [21,22,24,31,33,28,27,36,37,35,30]). The least-squares formulation allows for the construction of finite element models for fluids that, when combined with high-order finite element technology [22,4,5,38,17,29,31,31,49] possess many of the attractive qualities associated with the well-known Ritz method [43] such as global minimization, best approximation with respect to a well-defined norm, and symmetric positive-definiteness of the resulting finite element coefficient matrix [9].…”