HYBRID FIRST-ORDER SYSTEM LEAST SQUARES FINITE ELEMENT METHODS WITH APPLICATION TO STOKES EQUATIONS

作者：Liu K^*; Manteuffel T A; McCormick S F; Ruge J W; Tang L

来源：SIAM Journal on Numerical Analysis, 2013, 51(4): 2214-2237.

DOI：10.1137/120868906

• least squares finite element method
• error measure
• nested iteration
• adaptive local refinement
• stokes equations
• mass conservation

This paper combines first-order system least squares (FOSLS) with first-order system LL* (FOSLL*) to create a Hybrid method. The FOSLS approach minimizes the error, e(h) = u(h) - u, over a finite element subspace, V-h, in the operator norm: min(uh is an element of Vh)parallel to L(u(h) - u)parallel to. The FOSLL* method looks for an approximation in the range of L*, setting uh = L*w(h) and choosing w(h) is an element of W-h, a standard finite element space. FOSLL* minimizes the L-2 norm of the error over L*(W-h): min(wh is an element of Wh) parallel to L*w(h) - u parallel to. FOSLS enjoys a locally sharp, globally reliable, and easily computable a posteriori error estimate, while FOSLL* does not. However, FOSLL* has the major advantage that it applies to problems that do not exhibit enough smoothness to enable the full advantages that the FOSLS approach otherwise provides. The Hybrid method attempts to retain the best properties of both FOSLS and FOSLL*. This is accomplished by combining the FOSLS functional, the FOSLL* functional, and an intermediate term that draws them together. The Hybrid method produces an approximation, u(h), that is nearly the optimal over V-h in the graph norm, parallel to e(h)parallel to(2)(G) :=1/2 parallel to e(h)parallel to(2)+parallel to Le(h)parallel to(2). The FOSLS and intermediate terms in the Hybrid functional provide a very effective a posteriori error measure. This paper establishes that the Hybrid functional is coercive and continuous in a graph-like norm with modest coercivity and continuity constants, c(0) = 1/3 and c(1) = 3; that both parallel to e(h)parallel to and parallel to Le(h)parallel to converge with rates based on standard interpolation bounds; and that if LL* has full H-2 regularity, the L-2 error, parallel to e(h)parallel to, converges with a full power of the discretization parameter, h, faster than the functional norm. Letting (u) over tilde (h) denote the optimum over V-h in the graph norm, this paper also shows that if superposition is used with nested iteration, then parallel to u(h) - u(h)parallel to(G) converges two powers of h faster than the functional norm. Numerical tests are provided to confirm the efficiency of the Hybrid method and effectiveness of the a posteriori error measure.

• 出版日期2013