Least-squares methods for the velocity-pressure-stress formulation of the Stokes equations

• Published in: Computer methods in applied mechanics and engineering Vol. 126; no. 3; pp. 267 - 287
• Main Authors: Bochev, Pavel B., Gunzburger, Max D.
• Format: Journal Article
• Language: English
• Published: Amsterdam Elsevier B.V 01.10.1995; Elsevier
• Subjects: Computational methods in fluid dynamics; Exact sciences and technology; Fluid dynamics; Fundamental areas of phenomenology (including applications); Physics; Finite element method; Error estimation; Stokes problem; Stress tensor; Incompressible fluid; Least square fit; Pressure; Velocity; Numerical convergence
• ISSN: 0045-7825; 1879-2138
• DOI: 10.1016/0045-7825(95)00826-M

We formulate and study finite element methods for the solution of the incompressible Stokes equations based on the application of least-squares minimization principle to an equivalent first order velocity-pressure-stress system. Our least-squares functional involves the L 2-norms of the residuals of each equation multiplied by a mesh-dependent weight. Each weight is determined according to the Agmon-Douglis-Nirenberg index of the corresponding equation. As a result, the approximating spaces are not subject to the LBB condition and conforming discretizations are possible with merely continuous finite element spaces. Moreover, the resulting discrete problems involve only symmetric, positive definite systems of linear equations, i.e. assembly of the discretization matrix is not required even at the element level. We prove that the least-squares approximations converge to the solutions of the Stokes problem at the best possible rate and then present some numerical examples illustrating our theoretical results. Among other things, these numerical examples indicate that the method is not optimal without the weights in the least-squares functional.