UZAWA TYPE ALGORITHMS FOR NONSYMMETRIC SADDLE POINT PROBLEMS
|
|
- Gavin Rose
- 6 years ago
- Views:
Transcription
1 MATHEMATICS OF COMPUTATION Volume 69, Number 230, Pages S (99) Article electronically published on May 17, 1999 UZAWA TYPE ALGORITHMS FOR NONSYMMETRIC SADDLE POINT PROBLEMS JAMES H BRAMBLE, JOSEPH E PASCIAK, AND APOSTOL T VASSILEV Abstract In this paper, we consider iterative algorithms of Uzawa type for solving linear nonsymmetric saddle point problems Specifically, we consider systems, written as usual in block form, where the upper left block is an invertible linear operator with positive definite symmetric part Such saddle point problems arise, for example, in certain finite element and finite difference discretizations of Navier Stokes equations, Oseen equations, and mixed finite element discretization of second order convection-diffusion problems We consider two algorithms, each of which utilizes a preconditioner for the operator in the upper left block Convergence results for the algorithms are established in appropriate norms The convergence of one of the algorithms is shown assuming only that the preconditioner is spectrally equivalent to the inverse of the symmetric part of the operator The other algorithm is shown to converge provided that the preconditioner is a sufficiently accurate approximation of the inverse of the upper left block Applications to the solution of steady-state Navier Stokes equations are discussed, and, finally, the results of numerical experiments involving the algorithms are presented 1 Introduction This paper provides an analysis for Uzawa type methods applied to the solution of linear nonsymmetric saddle point systems Such systems arise in certain discretizations of Navier Stokes equations and mixed discretizations of second order elliptic problems with convective terms (cf [9], [11], [14], [17]) The theory in this paper is an extension of that for symmetric saddle point problems developed in [4] Let H 1 and H 2 be finite dimensional Hilbert spaces with inner products which we shall denote by (, ) There is no ambiguity even though we use the same notation for the inner products on both of these spaces, since the particular inner product will be identified by the type of functions appearing We consider the system ( )( ) ( A B T X F (11) =, B 0 Y G) Received by the editor July 8, 1997 and, in revised form, June 23, Mathematics Subject Classification Primary 65N22, 65N30, 65F10 Key words and phrases Indefinite systems, iterative methods, preconditioners, saddle point problems, nonsymmetric saddle point systems, Navier-Stokes equations, Oseen equations, Uzawa algorithm This manuscript was written under contract number DE-AC02-76CH00016 with the US Department of Energy Accordingly, the US Government retains a non-exclusive, royalty-free license to publish or reproduce the published form of this contribution, or allow others to do so, for US Government purposes This work was also supported in part under the National Science Foundation Grant No DMS and by the Schlumberger Technology Corporation 667 c 2000 American Mathematical Society
2 668 J H BRAMBLE, J E PASCIAK, AND A T VASSILEV where F H 1 and G H 2 are given and X H 1 and Y H 2 are the unknowns Here A : H 1 H 1 is assumed to be a linear, nonsymmetric operator A T : H 1 H 1 is the adjoint of A with respect to the (, ) inner product In addition, the linear map B takes H 1 into H 2 and its adjoint, B T,takesH 2 into H 1 In general, (11) may not be solvable unless additional conditions on the operators A and B and the spaces H 1 and H 2 are imposed Throughout this paper we assume that A has a positive definite symmetric part Under this assumption, (11) is solvable if and only if the reduced problem (12) BA 1 B T Y = BA 1 F G is solvable In the case of a symmetric and positive definite operator A, the Ladyzhenskaya Babuška Brezzi (LBB) condition (cf [5]) is a necessary and sufficient condition for solvability of this problem As we shall see, the solvability of (11) in the nonsymmetric case is guaranteed provided that the LBB condition holds for the symmetric part of A The papers [7], [15] propose solving BA 1 B T by a preconditioned iteration One common problem with these approaches is that the evaluation of the action of the operator A 1 is required in each step of the iteration For many applications, this operation is expensive and is also implemented as an iteration The Uzawa method [1] is a particular implementation of a linear iterative method for solving (12) It is an exact algorithm in the sense that the action of A 1 is required for the implementation An alternative method which solves (11) by preconditioned iteration was proposed in [8] Their preconditioner also requires the evaluation of A 1 at each step of the iteration The Uzawa type methods studied here replace the exact inverse of A by an approximate evaluation of A 1 or a preconditioner for its symmetric part Such algorithms are defined in Sections 3 and 4 In this paper we distinguish two types of algorithms: (i) a linear one-step method, where the action of the inverse is replaced by a linear preconditioner such as one sweep of a multigrid procedure; (ii) a multistep method, where a sufficiently accurate approximation to A 1 is provided by some iterative method, eg, preconditioned GMRES [16] or preconditioned Lanczos [12] The Uzawa type algorithms applied to nonsymmetric problems are of interest because they are simple, efficient, and have minimal computer memory requirements They can be applied to the solution of difficult practical problems such as the Navier Stokes equation In addition, an exact Uzawa algorithm implemented as a double iteration can be easily modified to be an algorithm of the type studied here The paper is organized as follows In Section 2 we establish sufficient conditions for solvability of the abstract saddle point problem and analyze an exact Uzawa algorithm for solving it In Section 3 we define and analyze a linear one-step Uzawa type algorithm Next, a multistep inexact method is defined and analyzed in Section 4 Section 5 provides applications of the algorithms from Section 3 and Section 4 to the solution of indefinite systems of linear equations arising in finite element approximations of the steady-state Navier-Stokes equations Finally, the results of numerical computations involving the algorithms are given in Section 6
3 UZAWA TYPE ALGORITHMS II Analysis of the exact method In this section we establish sufficient conditions for solvability of (12) and analyze the exact Uzawa algorithm for the computation of its solution The analysis of this method and, in particular, the result of Theorem 22 below is important for the analysis of the algorithms defined in the subsequent sections The symmetric part A s of the operator A is defined by (21) A s = 1 2 (A + AT ) In the remainder of this paper a subscript s will be used to denote the symmetric part of various operators, defined as in (21) We assume that A s is positive definite and satisfies (22) (AX, Y ) α(a s X, X) 1/2 (A s Y,Y ) 1/2, for all X, Y H 1, for some number α Clearly, α 1 Moreover, since A s is positive definite, such an α always exists In many applications involving the numerical solution of partial differential equations, the constant α can be chosen independently of the mesh parameter In addition, the Ladyzhenskaya Babuška Brezzi condition is assumed to hold for the pair of spaces H 1 and H 2, ie (V,BU) 2 (23) sup U H 1 (A s U, U) c 0 V 2, for all V H 2, for some positive number c 0 Here denotes the norm in the space H 2 (or H 1 ) corresponding to the inner product (, ) As is well known, the condition (23) is sufficient to guarantee solvability of (11) when A is symmetric We will see that it also suffices in the case of nonsymmetric A To this end, we prove the following lemma Lemma 21 Suppose that A is an invertible linear operator with positive definite symmetric part A s that satisfies (22) Then(A 1 ) s is positive definite and satisfies (24) ((A 1 ) s W, W ) ((A s ) 1 W, W ) α 2 ((A 1 ) s W, W ), for all W H 1 Proof Clearly, ((A s ) 1 (W, U ) 2 W, W )= sup U H 1 (A s U, U) = sup ((A 1 ) T W, AU) 2 U H 1 (A s U, U) (25) α 2 (A 1 ) T W 2 A sup s U 2 A s U H 1 U 2 A s = α 2 ((A 1 ) s W, W ) = α 2 (A 1 ) T W 2 A s Here 2 A s =(A s, ) In the above inequalities we have used the Schwarz inequality, (22), and the fact that (26) (A s U, U) =(AU, U), for all U H 1 On the other hand, ((A 1 ) s U, U) =(A 1 U, U) =(A 1/2 s A 1 U, (A s ) 1/2 U) A 1 U As U (As) 1 =(A 1 U, U) 1/2 U (As) 1
4 670 J H BRAMBLE, J E PASCIAK, AND A T VASSILEV Therefore, (27) ((A 1 ) s U, U) ((A s ) 1 U, U) This completes the proof of the lemma It is now clear that Lemma 21 and (23) guarantee solvability of (12) Indeed, (BA 1 B T V,V )=((A 1 ) s B T V,B T V ) α 2 ((A s ) 1 B T V,B T V ) α 2 c 0 V 2 Thus, we have proved the following theorem Theorem 21 Suppose that the linear operator A is invertible and that (23) holds Then the reduced problem (12), or equivalently (11), is solvable Next, we turn to the analysis of the exact Uzawa algorithm applied to the solution of (12) The preconditioned variant of the exact Uzawa algorithm (cf [1, 4]) is defined as follows Algorithm 21 (Preconditioned exact Uzawa) For X 0 H 1 and Y 0 H 2 given, the sequence {(X i,y i )} is defined, for i =0, 1, 2,,by X i+1 = X i + A 1 ( F (AX i + B T Y i ) ), Y i+1 = Y i + τq 1 B (BX i+1 G) Here the preconditioner Q B : H 2 H 2 is a symmetric positive definite linear operator satisfying (28) γ(q B W, W ) (B(A s ) 1 B T W, W ) (Q B W, W ), for all W H 2, for some γ in the interval (0, 1], and τ is a positive parameter Notice that this condition implies appropriate scaling of Q B In many applications effective preconditioners that satisfy (28) with γ bounded away from zero are known Let (29a) and E X i = X X i (29b) E Y i = Y Y i be the iteration errors generated by the above method Note that Ei+1 Y =(I τq 1 B BA 1 B T )Ei Y Therefore, the convergence of Algorithm 21 is governed by the properties of the operator I τq 1 B BA 1 B T summarized in the following theorem Theorem 22 Suppose that A is invertible with positive definite symmetric part A s which satisfies (22) Suppose also that (23) holds In addition, let Q B be a symmetric and positive definite operator satisfying (28) Ifτ is a positive parameter with τ γ α 2,then (I τq 1 B BA 1 B T )U 2 Q B (1 γ ) (210) α 2 τ U 2 Q B, for all U H 2
5 UZAWA TYPE ALGORITHMS II 671 Remark 21 If A = A T,thenτ may be taken equal to one and (28) implies (cf [4]) that (I Q 1 B BA 1 B T )U 2 Q B (1 γ) 2 U 2 Q B Hence, Theorem 22 is not optimal in the limit when α 1 Proof of Theorem 22 The proof is based on Lemma 21 Let L = BA 1 B T Then, by (28) and Lemma 21, γ V 2 Q B ((A s ) 1 B T V,B T V ) (211) α 2 (LV,V ) In addition,using (24), (A 1 v, w) =((A s ) 1/2 A 1 v, (A s ) 1/2 w) (212) (A 1 v, v) 1/2 ((A s ) 1 w, w) 1/2 ((A s ) 1 v, v) 1/2 ((A s ) 1 w, w) 1/2 Taking v = B T V and w = B T W above gives (213) (LV,W) V QB W QB Next, (214) (I τq 1 B L)V 2 Q B = V 2 Q B 2τ(LV,V )+τ 2 (LV,Q 1 B LV ) By (212) and (211), the last term in the right hand side of (214) is estimated by (LV,Q 1 B LV ) V Q B Q 1 B LV Q B (215) = V QB (LV,Q 1 B LV )1/2 Using (211) and (215) in (214) yields ( (I τq 1 B L)V 2 Q B 1 2τγ ) α 2 + τ 2 V 2 Q B This concludes the proof of the theorem 3 Analysis of the linear one-step method In this section we define and analyze a linear one-step Uzawa type algorithm applied to (11) This section contains the main result of the paper We show that, under the minimal assumptions needed to guarantee solvability (cf Section 2), appropriately scaled linear preconditioners (cf (28) and (31) below) lead to an efficient and simple method for solving (11) The exact inverse of A is replaced by a preconditioner for the symmetric part of A Let A 0 : H 1 H 1 be a linear, symmetric, positive definite operator that satisfies (31) (A 0 V,V ) (A s V,V ) β(a 0 V,V ), for all V H 1, for some β 1 Remark 31 The inequalities (28) and (31) respectively imply scaling of Q B and A 0 In practice, the proper scaling these operators can be achieved using even crude estimates for the largest eigenvalues of à 1 0 A 1 s and Q B B(A s) 1 B T,where à 0 and Q B are unscaled preconditioners Usually, a few iterations of the power method are enough for obtaining such estimates
6 672 J H BRAMBLE, J E PASCIAK, AND A T VASSILEV The linear Uzawa type algorithm is defined as follows Algorithm 31 (Linear one-step method) For X 0 H 1 and Y 0 H 2 given, the sequence {(X i,y i )} is defined, for i =0, 1, 2,,by X i+1 = X i + δa 1 ( 0 F (AXi + B T Y i ) ), Y i+1 = Y i + τq 1 B (BX i+1 G) Here δ and τ are positive parameters We will assume that δ<1/β It then follows from (31) that A 0 δa s is positive definite The following theorem is the main result of this paper Theorem 31 Suppose that A has a positive definite symmetric part, A s, satisfying (22) Suppose also that Q B and A 0 are symmetric positive definite operators satisfying (28) and (31) Then Algorithm 31 converges if 0 <δ (3α 2 β 2 ) 1 and 0 <τ (4β) 1 Moreover, if (X, Y ) is the solution of (11) and (X i,y i ) is the approximation defined by Algorithm 31, then the iteration errors Ei X and Ei Y defined in (29) satisfy (32) {δ 1 E X i 2 A 0 δa s + τ 1 E Y i 2 Q B } 1/2 ρ i { δ 1 E X 0 2 A 0 δa s + τ 1 E Y 0 2 Q B } 1/2 for any i 1 Here ρ = δ/2 δτγ + (δ/2 δτγ) 2 +4(1 δ/2) 2 Remark 32 Convergence of Algorithm 31 follows from (32) Indeed, a simple algebraic manipulation using the fact that τγ 1/4 gives ρ δ/2 δτγ + (δ/2 δτγ) 2 +4(1 δ/2) < 1 δτγ 2 2 The quantity on the right hand side is clearly less than one Remark 33 The use of a symmetric preconditioner in Algorithm 31 results in a fundamental change in the convergence properties For example, if we take τ = γ α 2 in Algorithm 21, then Theorem 22 gives the reduction of 1 γ 2 /α 4 In contrast, if we set δ =(3α 2 β 2 ) 1 and τ =(4β) 1, then Theorem 31 gives a convergence rate which is bounded by γ 1 24α 2 β 3 This behaves significantly better in applications involving large α In order to analyze Algorithm 31 we reformulate it in terms of the iteration errors defined in (29) It is easy to see that Ei X and Ei Y satisfy the following equations: Ei+1 X = EX i δa 1 ( 0 AE X i + B T Ei Y )), Ei+1 Y = ( I δτq 1 ) B BA 1 0 BT Ei Y + τq 1 B B(I δa 1 0 A)EX i For convenience, these equations can be written in matrix form as ( E X ) ( i+1 (I δa 1A) δa 1 E X) i (33) E Y i+1 = τq 1 B 0 0 BT B(I δa 1 0 A) (I τδq 1 B BA 1 0 BT ) E Y i
7 UZAWA TYPE ALGORITHMS II 673 Straightforward manipulations of (33) give (34) N E i+1 = ME i, where ( ) E X E i = i Ei Y, ( ) δ N = 1 (A 0 δa T ) 0 0 τ 1, Q B and M = δ 1 (A 0 δa T )A 1 0 (A 0 δa) (A 0 δa T )A 1 0 BT BA 1 0 (A 0 δa) (τ 1 Q B δba 1 0 BT ) It is clear that we can study the convergence of Algorithm 31 by investigating the properties of the linear operators M and N We shall reduce this problem to estimation of the spectral radius of related symmetric operators Let M 1 be the symmetric matrix defined by M 1 = JM, where ( ) I 0 J = 0 I Our next lemma reduces the proof of the theorem to the estimation of the eigenvalues of the generalized eigenvalue problem (35) λn s ψ = M 1 ψ Since δ is less than 1/β, N s is positive definite and the above problem is well defined Because N s and M 1 are symmetric, the eigenvalues λ are real Lemma 31 The iteration error E i satisfies (N s E i+1,e i+1 ) 1/2 ρ(n s E i,e i ) 1/2, where ρ =max i λ i,with{λ i } the eigenvalues of (35) Proof Let {(λ i,ψ i )} be the eigenpairs for (35) Since N s is positive definite, {ψ i } spans the space H 1 H 2 We may choose the eigenvectors so that (N s ψ i,ψ j )=δ ij, where δ ij denotes the Kronecker delta Now any vectors v and w in H 1 H 2 can be represented as v = v i ψ i and w = w i ψ i Thus, i i (M 1 v, w) = ij v i w j (M 1 ψ i,ψ j )= i v i w i λ i (36) ρ ( ) 1/2 ( vi 2 i i = ρ v Ns w Ns w 2 i ) 1/2
8 674 J H BRAMBLE, J E PASCIAK, AND A T VASSILEV Since J 2 is the identity operator we have that M = JM 1 Therefore, using (34) we see that (N s E i+1,e i+1 )=(ME i,e i+1 )=(M 1 E i, J E i+1 ) ρ E i Ns J E i+1 Ns = ρ E i Ns E i+1 Ns The lemma immediately follows Our proof of Theorem 31 will require another lemma We need to provide some control on the convergence of the related linear iteration (37) to the solution U of U i+1 = U i + δa 1 0 (W AU i) AU = W Lemma 32 Let A 0 satisfy (31) and δ be a positive number with δ<1/β Then (38) where (I δa 1 0 A)V 2 A 0 δ ((A 0 δa s )V,V ), for all V H 1, δ =1 δ + α2 β 2 δ 2 1 δβ Remark 34 If, in addition (39) so that δ< 1 α 2 β 2 + β, α 2 β 2 δ 1 δβ < 1, then δ is less than one Proof of Lemma 32 By (31), (1 δβ)(a 0 V,V ) ((A 0 δa s )V,V ), for all V H 1 Hence, by (22) and (31), (310) (AV,W) α(a s V,V ) 1/2 (A s W, W ) 1/2 αβ (1 δβ) 1/2 (A 0V,V ) 1/2 ((A 0 δa s )W, W ) 1/2 On the other hand, (311) (I δa 1 0 A)V 2 A 0 = V 2 A 0 2δ(AV,V )+δ 2 (A 1 0 AV,AV ) =((A 0 δa s )V,V ) δ(av,v )+δ 2 (A 1 0 AV,AV ) In view of (31), we have (312) (AV,V ) (A 0 V,V ) ((A 0 δa s )V,V ) Also, (310) implies that (A 1 0 AV,AV ) αβ (A 1 (1 δβ) 1/2 0 AV,AV )1/2 ((A 0 δa s )V,V ) 1/2
9 UZAWA TYPE ALGORITHMS II 675 Thus, (A 1 0 AV,AV ) α2 β 2 (313) 1 δβ ((A 0 δa s )V,V ) Using (312) and (313) in (311) yields (38) Proof of Theorem 31 Let δ and τ satisfy the hypotheses of the theorem Because of Lemma 31 it suffices to bound the eigenvalues of (35) We begin with the negative eigenvalues Let (χ, ξ) be an eigenvector in H 1 H 2 with eigenvalue λ<0 Then multiplying the first equation of (35) by A 0 (A 0 δa T ) 1 gives (314) λδ 1 A 0 (A 0 δa T ) 1 (A 0 δa s )χ = δ 1 (A 0 δa)χ + B T ξ The second equation of (35) is (315) λτ 1 Q B ξ = BA 1 0 (A 0 δa)χ +(τ 1 Q B δba 1 0 BT )ξ Applying δba 1 0 to (314) and adding it to (315), we obtain (316) (1 λ)τ 1 Q B ξ = λb(a 0 δa T ) 1 (A 0 δa s )χ Eliminating ξ between (314) and (316), we find that 1 λ (A 0 δa)χ + δτ 1 λ BT Q 1 B B(A 0 δa T ) 1 (A 0 δa s )χ = A 0 (A 0 δa T ) 1 (A 0 δa s )χ Taking the inner product with (A 0 δa T ) 1 (A 0 δa s )χ yields (317) 1 λ ((A 0 δa s )χ, χ)+ δτ 1 λ B(A 0 δa T ) 1 (A 0 δa s )χ 2 Q 1 B = (A 0 δa T ) 1 (A 0 δa s )χ 2 A 0 For convenience, the last equation can be abbreviated as T 1 + T 2 = T 3 In order to bound T 2 we note that for any φ H 1, (318) (Q 1 (φ, B T ζ) 2 B Bφ, Bφ) = sup ζ H 2 (Q B ζ,ζ) = sup ((A s ) 1/2 φ, (A s ) 1/2 B T ζ) 2 ζ H 2 (Q B ζ,ζ) (A s φ, φ)(b(a s ) 1 B T ζ,ζ) sup (A s φ, φ) ζ H 2 (Q B ζ,ζ) For the last inequality above we used (28) Therefore, T 2 Using this in (317) gives ( (319) δτ 1 λ (A 0 δa T ) 1 (A 0 δa s )χ 2 A s βδτ 1 λ T 3 1 δτβ 1 λ ) T 3 1 λ ((A 0 δa s )χ, χ) By Lemma 32, for any φ H 1,wehave ( (320) (A0 δa T )A 1 0 (A 0 δa)φ, φ ) δ((a 0 δa s )φ, φ),
10 676 J H BRAMBLE, J E PASCIAK, AND A T VASSILEV which is equivalent to (321) ((A 0 δa s ) 1 φ, φ) δ((a 0 δa) 1 A 0 (A 0 δa T ) 1 φ, φ) for any φ H 1 Takingφ =(A 0 δa s )χ in (321), we obtain T 3 1 δ ((A 0 δa s )χ, χ) Combining this with (319) and using the fact that λ<0gives 1 λ ((A 0 δa s )χ, χ) 1 δ ( 1 δτβ ) ((A 0 δa s )χ, χ) 1 λ (322) 1 δτβ ((A 0 δa s )χ, χ) δ Now if χ were equal to zero then (314) would imply that B T ξ = 0, but (315), in turn, would imply that ξ = 0 Hence, since (χ, ξ) is an eigenvector, χ 0and therefore ((A 0 δa s )χ, χ) 0 Thus,from(322), (323) δ λ 1 δτβ Applying straightforward manipulations, we get δ =1 δ + α2 β 2 δ 2 ( 1 δβ 1 δ 1 1/3 ) (324) =1 δ 1 1/3 2 and (325) 1 1 δτβ 1 1 δ/4 Using (324) and (325) in (323) gives (326) λ 1 δ/2 1 δ/4 1 δ 4, which provides a bound for the negative part of the spectrum Next we obtain a bound for the positive eigenvalues of (35) To this end we factor M 1 as M 1 = D T M 2 D, where ( θ 1/2 (A 0 ) 1/2 ) (A 0 δa) 0 D =, 0 I M 2 = δ 1 θi θ 1/2 (A 0 ) 1/2 B T, θ 1/2 B(A 0 ) 1/2 τ 1 Q B δba 1 0 BT and θ =1 δ/2 The largest eigenvalue of (35) is given by (M 1 w, w) Λ= sup w H 1 H 2 (N s w, w) = sup (M 2 Dw, Dw) w H 1 H 2 (N s w, w) In order to obtain an upper bound for Λ, we first note that by (320) and (324) it follows that (327) θ 1 A 1/2 0 (A 0 δa)χ 2 ((A 0 δa s )χ, χ),
11 UZAWA TYPE ALGORITHMS II 677 and thus (328) ( ( ) ( θ 1 δ 1 A 1/2 0 (A 0 δa)χ 2 + τ 1 ξ 2 χ χ Q B N s, ξ ξ)) Thus, it suffices to show that for any vector (φ, ζ) H 1 H 2, ( ( ( ] φ φ M 2 D, D ρ [(δθ) ζ) ζ)) 1 A 1/2 0 (A 0 δa)φ 2 + τ 1 ζ 2 QB (329) = ρ (( δ 1 I 0 0 τ 1 Q B ) ( ( ) φ φ D, D ζ) ) ζ Then ρ will be an upper bound for Λ To this end let L = B(A 0 ) 1/2 Now M 2 may be written as ( ) δ M 2 = 1 θi θ 1/2 L T θ 1/2 L τ 1 Q B δll T The proof of (329) is now reduced to estimating the largest eigenvalue, λ, with eigenvector (χ, ξ) satisfying (330) θδ 1 χ + θ 1/2 L T ξ = λδ 1 χ and (331) θ 1/2 Lχ +(τ 1 Q B δll T )ξ = λτ 1 Q B ξ Solving for χ in (330), we get χ = δ(λ + θ) 1 θ 1/2 L T ξ Substituting this in (331) yields (1 λ)(λ + θ)q B ξ = δτλll T ξ Hence (332) (1 λ)(λ + θ)(q B ξ,ξ) =δτλ(l T ξ,l T ξ) Now ξ cannot be zero, since if it were, then (330) would imply that either χ =0 or λ = θ 0 Hence, since λ is a positive eigenvalue, it follows that ξ 0 In addition, by (31) and (28), (L T ξ,l T ξ)=(a 1 0 BT ξ,b T ξ) ((A s ) 1 B T ξ,b T ξ) γ(q B ξ,ξ) Using this in (332) gives (1 λ)(λ + θ) δτλγ, or equivalently λ 2 λ(1 θ δτγ) θ 0 From here we obtain that λ 1 θ δτγ + ((1 θ) δτγ) 2 +4θ 2 (333) = δ/2 δτγ + (δ/2 δτγ) 2 +4(1 δ/2) 2
12 678 J H BRAMBLE, J E PASCIAK, AND A T VASSILEV Finally, elementary inequalities imply that 1 δ 4 δ/2 δτγ + (δ/2 δτγ) 2 +4(1 δ/2), 2 which concludes the proof of the theorem 4 Analysis of the multistep inexact algorithm In this section we define and analyze an inexact Uzawa algorithm obtained by replacing A 1 with a sufficiently accurate approximation Such an algorithm is essentially different from the linear one-step method developed in the previous section for two main reasons First, achieving a certain accuracy of the approximation to A 1 typically requires more computational work than one evaluation of the action of a preconditioner Second, depending on the manner in which the accurate approximate inverse is computed, the resulting inexact Uzawa algorithm may not be linear In view of this, we shall approach the analysis of this method differently The approximate inverse may be described as a map Ψ : H 1 H 1, not necessarily linear In this section we shall assume that for any φ H 1,Ψ(φ) is close to the solution ξ of (41) Aξ = φ More precisely, we assume that (42) Ψ(φ) A 1 φ As ɛ A 1 φ As, for all φ H 1, for some positive ɛ with ɛ<1 Notice that for any ɛ (0, 1), (42) can be satisfied by taking sufficiently many steps of some iterative method for solving (41) which reduces the error in a norm equivalent to As For example, for an appropriate choice of the iteration parameter, the linear iteration (37) converges to the solution of the linear system (41) (cf Remark 34) Hence, an estimate of the type of (42) can be established for any ɛ<1, provided that sufficiently many iterations with (37) are performed Another example of Ψ results from a preconditioned generalized Lanczos procedure [12] In this case the resulting Uzawa algorithm will be nonlinear As an example, we consider the generalized minimal residual algorithm (GMRES) Specifically, let A 0 satisfy (31) and consider the GMRES algorithm applied in the inner product (, ) A0 (A 0, ) to the preconditioned equation A 1 0 Aξ = A 1 0 φ Using the initial iterate ξ 0 =0,wesetΨ(φ) =ξ n,whereξ n is the approximation obtained after n steps The GMRES method computes the best approximation to ξ (in the norm A 1 0 A A 0 )inthekrylovspacek n =span{(a 1 0 A)i ξ, i =1,,n} Thus, it follows from Lemma 32 that A 1 0 A(ξ n ξ) A0 (I δa 1 0 A)n A 1 0 Aξ A 0 δ n/2 A 1 0 Aξ A 0 Applying Lemma 21 and (31) gives A 1 0 A(ξ n ξ) A0 ξ n A 1 φ As Similarly, A 1 0 Aξ A 0 αβ 1/2 A 1 φ As
13 UZAWA TYPE ALGORITHMS II 679 Combining the above inequalities, we obtain ξ n A 1 φ As αβ 1/2 δn/2 A 1 φ As This shows that, given ɛ, thereexistsn such that Ψ(φ) =ξ n satisfies (42) provided that (39) holds The variant of the inexact Uzawa algorithm which we investigate in this section is defined as follows Algorithm 41 (Multistep inexact Uzawa) For X 0 H 1 and Y 0 H 2 given, the sequence {(X i,y i )} is defined, for i =0, 1, 2,,by X i+1 = X i +Ψ ( F ( AX i + B T )) Y i, Y i+1 = Y i + τq 1 B (BX i+1 G) Algorithm 41 reduces to Algorithm 21 if Ψ(φ) =A 1 φ for all φ H 1 The main result of this section is a bound for the rate of convergence of the multistep algorithm in terms of the factors α, γ, andɛ introduced in (22), (28), and (42) respectively The theorem below gives a sufficient condition on ɛ for convergence of the algorithm Theorem 41 Suppose that A has a positive definite symmetric part and satisfies (22), andq B is a symmetric positive definite operator satisfying (28) Assume that (42) holds and that the parameter τ is chosen so that 0 <τ γ α 2 Set θ = ( 1 τ γ α 2 ) 1/2 Then the multistep inexact Uzawa algorithm converges if (43) ɛ< 1 θ 1+2τ θ Moreover, if (X, Y ) is the solution of (11) and (X i,y i ) is the approximation defined by Algorithm 41, then the iteration errors Ei X and Ei Y defined in (29) satisfy ( ) ɛτ ɛτ (44) 1+ɛ EX i+1 2 A s + Ei+1 Y 2 Q B ρ 2(i+1) 1+ɛ EX 0 2 A s + E0 Y 2 Q B and (45) where E X i+1 2 A s ( ) ɛτ τ 1 (1 + ɛ)(1 + 2ɛ)ρ 2i 1+ɛ EX 0 2 A s + E0 Y 2 Q B, ρ = (1 + τ)ɛ + θ + ((1 + τ)ɛ + θ) 2 +4ɛ(τ θ) (46) < 1 2 Proof We start by deriving norm inequalities involving the errors Ei X and Ei Y Similarly to the approach in the previous section, we can write Ei+1 X = EX i Ψ ( AEi X + B T Ei Y ), (47) E Y i+1 = EY i + τq 1 B BEX i+1
14 680 J H BRAMBLE, J E PASCIAK, AND A T VASSILEV Thefirstequationabovecanberewritten (48) Ei+1 X =(A 1 Ψ) ( AEi X + B T Ei Y ) A 1 B T Ei Y It follows from the triangle inequality, (42), (24), and (28) that Ei+1 X As ɛ( Ei X As + A 1 B T Ei Y As )+ A 1 B T Ei Y As (49) = ɛ E X i A s +(1+ɛ) B T E Y i (A 1 ) s ɛ E X i A s +(1+ɛ) E Y i Q B Using (48) in the second equation of (47) gives Ei+1 Y =(I τq 1 B BA 1 B T )Ei Y Applying the QB inequality yields (410) E Y i+1 QB + τq 1 B B(A 1 Ψ)(AE X i + B T E Y i ) norm to both sides of the above equation and using the triangle (I τq 1 B BA 1 B T )E Y i QB + τ Q 1 B B(A 1 Ψ)(AE X i + B T E Y i ) QB Since τ γ, by (210) we have α2 ( (I τq 1 B BA 1 B T )Ei Y QB 1 τ γ ) 1/2 E Y (411) α 2 i QB = θ Ei Y QB Because of (318), (42), the triangle inequality, and (28), the second term in the right-hand side of (410) is bounded as follows: (412) Q 1 B B(A 1 Ψ)(AEi X + B T Ei Y ) QB ɛ( Ei X As + Ei Y QB ) Using (411) and (412) in (410) yields (413) Ei+1 Y Q B θ Ei Y Q B + ɛτ( Ei X A s + Ei Y Q B ) Combining (49) and (413) gives Ei+1 X A s ɛ Ei X A s +(1+ɛ) Ei Y Q B, (414) Ei+1 Y A s ɛτ Ei X A s +(θ + ɛτ) Ei Y Q B Let us adopt the notation ( ) ( ) x1 x2 y 1 y 2 for vectors of nonnegative numbers x 1,x 2,y 1,y 2 if x 1 x 2 and y 1 y 2 Hence, from (414) we obtain ( ) ( )( ) E X i+1 As ɛ 1+ɛ E X i As (415) Ei+1 Y Q B ɛτ θ + ɛτ Ei Y Q B Repeated application of (415) gives ( ) ( ) E X i+1 As E X Ei+1 Y M i+1 0 As (416) Q B E0 Y, Q B where M is given by ( ) ɛ 1+ɛ R = ɛτ θ + ɛτ
15 UZAWA TYPE ALGORITHMS II 681 We consider two dimensional Euclidean space with the inner product ( ) ( ) x1 x2, = ɛτ y 1 y 2 1+ɛ x 1x 2 + y 1 y 2 A trivial computation shows that R is symmetric with respect to the inner product It follows from (416) that ( ) ( ) ɛτ E X 1+ɛ EX i+1 2 A s + Ei+1 Y i+1 As E X i+1 As 2 Q B = Ei+1 Y, Q B Ei+1 Y Q B ( ) ( ) E X R i+1 0 As E X, R i+1 0 As E0 Y QB E0 Y QB ( ) ɛτ ρ 2(i+1) 1+ɛ EX 0 2 A s + E0 Y 2 Q B, where ρ is the norm of the matrix R with respect to the, -inner product Since R is symmetric in this inner product, its norm is equal to its spectral radius The eigenvalues of R are the roots of λ 2 ((1 + τ)ɛ + θ) λ ɛ(τ θ) =0 It is elementary to see that the root with largest absolute value is that given by (46) For any fixed positive τ and θ in the interval [0, 1], ρ is a function of ɛ only It is straightforward to check that ρ = 1 only if ɛ = 1 θ 1+2τ θ 1 θ 1+2τ θ ) Moreover, ρ = θ if ɛ =0 Thus,ρ<1forɛ [0, Finally, we prove (45) Multiplying both sides of the first inequality in (414) by τ 1/2 and using the fact that 0 <τ<1, we obtain τ 1/2 Ei+1 X A s τ 1/2 ɛ Ei X A s + τ 1/2 (1 + ɛ) Ei Y Q B τ 1/2 ɛ Ei X A s +(1+ɛ) Ei Y Q B We now apply the arithmetic-geometric mean inequality to the last inequality and get that, for any positive η, τ Ei+1 X 2 A s (1 + η)τɛ 2 Ei X 2 A s +(1+η 1 )(1 + ɛ) 2 Ei Y 2 Q B Inequality (45) follows by taking η = 1 + 1/ɛ and applying (44) This completes the proof of the theorem We conclude this section with the following remarks Remark 41 If we fix all parameters except ɛ, then provided that the assumptions of the theorem are satisfied, ρ [0, 1) Moreover, ρ is a continuous function of ɛ which is equal to θ at ɛ = 0, ie, the theorem reproduces the result of Theorem 22 when ɛ = 0 We clearly can achieve any convergence rate between one and θ with an appropriate choice of ɛ Remark 42 Theorem 41 is somewhat weaker than the result obtained in Section 3 for the linear case due to the threshold condition (43) on ɛ As discussed above, it is possible to take sufficiently many iterations n, so that (43) holds for any fixed γ, α, and τ In applications involving discretizations of partial differential equations,
16 682 J H BRAMBLE, J E PASCIAK, AND A T VASSILEV β or α may depend on the discretization parameter h If, however, these parameters can be bounded independently of h, then a fixed number of preconditioned GMRES iterations (independent of h) are sufficient to guarantee convergence of Algorithm 41 5 Application to Navier-Stokes problems In this section we consider an application of the algorithms developed in the previous sections to the problem of solving indefinite systems of linear equations arising from finite element approximations of the steady-state Navier-Stokes equations We consider the following model problem for the steady-state Navier-Stokes equations: (51a) ν u +(u )u p = f in Ω, (51b) u =0 inω, (51c) u =0 on Ω, (51d) pdx=0 Ω Here Ω is a the unit square in R 2, u is a vector valued function representing the fluid velocity, and ν is the kinematic viscosity of the flow The fluid pressure p is a scalar function The pressure is determined only up to an additive constant, so for uniqueness, we require (51d) Generalizations to more complex domains and nonhomogenious boundary conditions are possible For example, we shall consider a problem with nonzero Dirichet boundary conditions in the next section Let Π be the set of functions in L 2 (Ω) with zero mean value on Ω and let H 1 (Ω) denote the Sobolev space of order one on Ω ([6, 13]) The space H 1 0 (Ω) consists of those functions in H 1 (Ω) whose traces vanish on Ω Also, V =(H 1 0 (Ω))2 will denote the product space consisting of vector valued functions with each component in H 1 0 (Ω) In order to derive the weak formulation of (51) we multiply the first two equations of (51) by functions in V and Π respectively and integrate over Ω to get (52a) (52b) νd(u, v)+b(u, u, v)+(p, v) =(f, v), for all v V, ( u,q)=0, for all q Π Here (, ) isthel 2 (Ω) inner product and D(, ) denotes the Dirichlet form for vector functions on Ω defined by 2 D(v, w) = v i w i dx i=1 The trilinear form b(,, ) for vector functions on Ω is given by b(u, v, w) =((u )v, w) The existence of a solution to (52) has been shown (cf [17], [9]) It is well known that the Navier-Stokes equations may have more that one solution unless the data (the kinematic viscosity and the external forces) satisfy very stringent requirements (cf [9], [17]) On the other hand, it has been shown that in many practical cases these solutions are mostly isolated, ie there exists a neighborhood of ν and f in Ω
17 UZAWA TYPE ALGORITHMS II 683 which each solution is unique We refer the reader to [9] and [17] for additional discussion We next define our finite element approximation subspaces The discussion here is very closely related to the examples given in [3] and [2], where additional comments and other applications can be found We partition Ω into 2n 2n squares, where n is a positive integer, and we define h =1/2n Letx i = ih and y j = jh for i, j = 1,,2n Each of the squares is further partitioned into two triangles by its diagonal with positive slope Let S h be the space of functions that are continuous and piecewise linear with respect to the triangulation just defined and vanish on Ω We set V h S h S h V The space Π h of functions that are piecewise constant with respect to the square elements and have zero mean value on Ω together with V h as defined above form an unstable pair of approximation spaces This means that the inequality (53) p c 0 sup V V h ( V,p) D(V,V ) 1/2, for all p Π h, fails to hold In fact, there is a function p Π h such that ( V,p) = 0 for all V V h Here(, ) denotes the inner product in L 2 (Ω) and is the corresponding norm To get a divergence stable pair, we consider a smaller space defined as follows Let η kl for k, l = 1,,2n be the function that is 1 on the square element [x k 1,x k ] [y l 1,y l ] and vanishes elsewhere Define φ ij Π h for i, j = 1,,n by φ ij = η 2i 1,2j 1 η 2i,2j 1 η 2i 1,2j + η 2i,2j (see Figure 1) The space Π h is then defined by { Π h W Π } h : (W, φ ij )=0fori, j =1,,n Now (53), with Π h replaced by Π h, is satisfied with a constant c 0 independent of h [10] Moreover, the exclusion of the functions φ i,j does not change the order of approximation, since Π h still contains the piecewise constant functions on squares of size 2h Following Temam [17], we introduce a modification b(,, ) of the trilinear form b(,, ), given by b(u, v, w) =((u )v, w) ((u )w, v) The approximation to the solution of (52) is defined by the pair (X, Y ) V h Π h satisfying (54a) (54b) νd(x, V )+ b(x, X, V )+(Y, V)=(f,V), for all V V h, ( X, W) =0, for all W Π h Note that the use of b(,, ) above is justified by the observation that b(u,, ) = b(u,, ) for functions u which are divergence free The form of b(,, ) guarantees the existence of a solution to (54) (cf [17]) The uniqueness requires conditions on the data ν and f To solve (54) we apply a Picard iteration of the following type (cf [11]) Given an initial approximation X 0,wecompute(X i,y i ), for i =1, 2,, asthesolution
18 684 J H BRAMBLE, J E PASCIAK, AND A T VASSILEV of the linear system (55a) νd(x i,v)+ b(x i 1,X i,v)+(y i, V )=(f,v), for all V V h, (55b) ( X i,w)=0, for all W Π h It is shown in [11] that the algorithm converges under the assumption that ν 2 c 2 a >c b f 1, where c a and c b are the coercivity and boundedness constants of the trilinear form b(,, ) Such an assumption is enough to guarantee a unique solution of (52) The system (55) can be reformulated in the notation of the earlier sections Set H 1 = V h and H 2 =Π h Let B : H 1 H 2, (BU, W )=( U, W ), for all U H 1,W H 2, and B T : H 2 H 1, (B T W, V )=(W, V ), for all V H 1,W H 2 At each iterative step, X i 1 is fixed so that we can define A : H 1 H 1, (AU, V )=νd(u, V )+ b(x i 1,U,V), for all U, V H 1 It follows that the solution (X i,y i ) of (55) satisfies (11) with F equal to the L 2 (Ω) projection of f onto H 1 and G = 0 Notice also that b(u, v, w) = b(u, w, v) Therefore, (56) A s : H 1 H 1, (A s U, V )=νd(u, V ), for all U, V H 1 It is possible to show that (22) holds for A and A s with a constant α proportional to ν 1 (cf [17] and [9]) Moreover, it follows from (56) that (23) holds for A s, B, and B T as above with constant c 0 independent of the mesh size h Thisimpliesthat (28) is satisfied with Q B = ν 1 I and γ bounded away from zero independently of h Figure 1 Thesquaremeshusedfor H 2 ; the support (shaded) and values for a typical φ ij
19 UZAWA TYPE ALGORITHMS II 685 We still need to provide a preconditioner for A s Clearly, A s consists of two copies of the operator which results from a standard finite element discretization of Dirichlet s problem There has been an intensive effort focused on the development and analysis of preconditioners for such problems For the examples in Section 6, we will use a preconditioning operator which results from a V-cycle variational multigrid algorithm Such a preconditioner can be scaled so that (31) holds with β independent of the mesh parameter h Remark 51 By rescaling p, one can rewrite (51a) in the form u + Re(u )u p = Re f, where Re = ν 1 is the Reynolds number of the flow This results in a different scaling of the discrete problem (54) which is better suited for implementation on finite precision machines We use this scaling in our examples in the next section Remark 52 An alternative linearization of (54) can be defined by replacing b(x, X, V )with b(x i 1,X i 1,V) which provides a different Picard iteration We will call this an explicit Picard iteration, because the nonlinear term is handled in an explicit fashion This leads to a symmetric saddle point problem at each iteration The inexact Uzawa methods analyzed in [4] can be used here Even though the symmetric linear systems are easier to solve, this linearization is a less robust method for computing solutions to (54) than the implicit linearization defined above, because the explicit Picard iteration breaks down for values of ν where the implicit method converges We shall provide a comparison of these two methods in the next section 6 Numerical examples In this section we present the results from numerical computations that illustrate the theory developed in the earlier sections Our goals here are first to demonstrate the efficiency and the robustness of the algorithm, and also to provide a comparison between the implicit and the explicit Picard iteration applied to a Navier-Stokes problem with known analytic solution In addition, we show results from computations of a classical flow problem The finite element discretization defined in the previous section as well as the pressure rescaling according to Remark 51 are used in both cases The first computations are for the solution of (54) when the velocity X is given by ( ) x(1 x)y(1 y) (61) X =, x(1 x)y(1 y) and the pressure Y is given by (62) Y = x 1 2 Obviously, X 0, so that the right-hand side of (54b) has to be adjusted appropriately We first give some results which indicate the behavior of the linear iteration method Starting the implicit nonlinear iteration with zero initial guess, we report the number of iterations needed to solve the nonsymmetric linear system on the second nonlinear Picard iteration
20 686 J H BRAMBLE, J E PASCIAK, AND A T VASSILEV Table 1 Linear iteration numbers on the second nonlinear iteration for the implicit method h Re=1 Re=10 Re=100 Re=1000 1/ / / / In all of our examples, we use one V cycle sweep of variational multigrid with point Gauss-Seidel smoothing to define A 1 0 This results in (31) being satisfied with β independent of h Thus, wecanchooseτ independent of h The preconditioner Q B was taken to be the identity Table 1 gives the number of iterations of the linear one-step method for a 10 8 normalized reduction of the residual The nonsymmetric term was relatively weak for the case of Re = 1 and 10, so the discrete problem was dominated by its symmetric part Thus, we set δ =1/β and τ = 1 consistent with the theory for the symmetric problem (cf [4]) Note that in Algorithm 31, A 1 0 always appears multiplied by δ, so that using δ =1/β with a preconditioner A 0 satisfying (31) is equivalent to using δ = 1 with the preconditioner Ã0 = βa 0 satisfying 1 (63) β (Ã0V,V ) (A s V,V ) (Ã0V,V ), for all V H 1 The unscaled multigrid preconditioner Ã0 automatically satisfies (63) As the value of Re increased further, the values of τ and δ had to be adjusted to obtain stability of the iteration The value τ = 01 was chosen To obtain stability in the case of Re = 100, we also needed to reduce δ to 01/β Finally, for stability in thecaseofre = 1000, we needed to reduce δ to 0001/β This is in agreement with the α 2 behavior required by the theory The number of iterations was bounded independently of the mesh parameter h, as suggested by the theory The value of δ clearly depends on the strength of the nonsymmetric term This, in turn, depends on the solution in nonlinear applications For example, for the driven cavity results given below, the case of h =1/128 and Re = 1000 required δ =00001/β to remain stable, and diverged for δ =0001/β The next set of experiments illustrates the differences between the implicit and explicit methods described in the previous section Three conditions were common in all experiments First, at each Picard iteration, the corresponding linear problem was solved exactly (ie the L 2 norm of the normalized residual was reduced until less than 10 8 ) Second, the nonlinear iteration was considered to have converged when the L 2 norm of the difference U i U i 1 was less than 10 6 HereU consists of both velocity and pressure components Finally, the Picard iteration was started with zero initial iterate The numerical results from these experiments are shown in Tables 2 4 In the case of Re = 1 and 10, the nonlinearity is small, and both the explicit and implicit methods work well In contrast, the explicit method fails to converge for Re = 100 while the implicit method behaves quite well Our second numerical experiment is the calculation of the flow in a cavity The cavity domain Ω is the unit square and the flow is caused by a tangential velocity field applied to one of the square sides in the absence of other body forces Since
21 UZAWA TYPE ALGORITHMS II 687 Table 2 Errors and nonlinear iteration numbers for Re =1 for the implicit and explicit methods h Error (p) Error (u 1 ) Error (u 2 ) Implicit Explicit 1/8 102e-2 787e-4 886e /16 250e-3 193e-4 541e /32 618e-4 481e-5 299e /64 193e-4 126e-5 157e Table 3 Errors and nonlinear iteration numbers for Re =10for the implicit and explicit methods h Error (p) Error (u 1 ) Error (u 2 ) Implicit Explicit 1/8 106e-2 787e-4 886e /16 260e-3 193e-4 541e /32 643e-4 481e-5 299e /64 165e-4 125e-5 157e Table 4 Errors and nonlinear iteration numbers for Re = 100 for the implicit and explicit methods h Error (p) Error (u 1 ) Error (u 2 ) Implicit Explicit 1/8 315e-2 787e-4 885e /16 783e-3 193e-4 542e * 1/32 195e-3 480e-5 299e-3 8 ** 1/64 489e-4 120e-5 157e-3 8 ** * the algorithm converged to a different solution with corresponding errors (p, u 1, u 2 ) 781e-3, 205e-4, 537e-3 ** the algorithm could not converge to the solution all forces are independent of time, the flow in this case limits to a steady-state which is modeled by (51) with corresponding changes in the boundary conditions (51c) In particular, the solution u on the boundary is zero everywhere except on the boundary segment y =1,whereu =(1, 0) The corresponding discrete problem in the spaces H 1 and H 2 as defined in the previous section is similar to (54) and is given by D(X 0,V)+Re b(x, X 0,V)+(Y, V)=Re (f,v) D( ˆX,V ) Re b(x, ˆX,V ), ( X 0,W)=0, for all V H 1 and W H 2 Here X = X 0 + ˆX, withx 0 H 1 and ˆX satisfying the Dirichlet boundary conditions of u and vanishing at all interior vertex points from the triangulation of Ω Note that ˆX =0 The implicit Picard iteration for this nonlinear problem is given as follows Let ˆX be as defined above Then, given an initial iterate X 0 0, we compute (Xi 0,Yi ), for
22 688 J H BRAMBLE, J E PASCIAK, AND A T VASSILEV Figure 2 Streamlines for h =1/64, and Re = 1 (left); Re = 10 (right) Figure 3 Streamlines for h =1/64, Re = 100 (left); Re = 1000 (right) i =1, 2,, by D(X0 i,v)+re b(x i 1,X0 i,v)+(y i, V ) = Re (f,v) D( ˆX,V ) Re b(x i 1, ˆX,V ), ( X0,W)=0, i and set X i = X0 i + ˆX The streamlines of the velocity field X computed using this algorithm for a range of Reynolds numbers are shown in Figures 2 and 3 The effect of the Reynolds number on the flow pattern is clearly seen The flow for low Reynolds numbers (see Figure 2) has only one vortex center, located above the center of the domain (its location moves to the right as Re increases) As Re increases further, a second vortex center appears near the lower right corner (see Figure 3, the case of Re = 100), and, for even larger Reynolds numbers, a third vortex center develops near the lower left corner of the domain (see Figure 3, the case of Re = 1000)
23 UZAWA TYPE ALGORITHMS II 689 Again, the case of Re = 1000 was the most difficult problem, requiring a large amount of work in the linear solver for each Picard iteration The discretization with h =1/64 was sufficiently fine for resolving the essential flow behavior for all Reynolds numbers tested In contrast, the experimental results with h =1/16 and h =1/32 for Re = 100 did not show the vortex center near the lower right corner of the domain The experiment with h =1/128 and Re = 1000 resulted in a flow field whose streamlines were very similar to the ones from h =1/64 In conclusion, the implicit algorithm is a simple, robust and efficient method for solving Navier-Stokes equations for a wide range of Reynolds numbers For each nonlinear iteration it requires the solution of a nonsymmetric saddle point problem which can be solved effectively with the inexact Uzawa algorithm 31 References [1] K Arrow, L Hurwicz, and H Uzawa Studies in Linear and Nonlinear Programming Stanford University Press, Stanford, CA, 1958 MR 21:7115 [2] JH Bramble and JE Pasciak A preconditioning technique for indefinite systems resulting from mixed approximations of elliptic problems Math Comp, 50:1 18, 1988; 51: , 1988 MR 89m:65097 [3] JH Bramble and JE Pasciak Iterative techniques for time dependent Stokes problems Computers Math Applic, 33:13 30, 1997 MR 98e:65091 [4] JH Bramble, JE Pasciak, and AT Vassilev Analysis of the inexact Uzawa algorithm for saddle point problems SIAM J Numer Anal, 34: , 1997 MR 98c:65182 [5] FBrezziandMFortinMixed and Hybrid Finite Element Methods Springer-Verlag, New York, 1991 MR 92d:65187 [6] PG Ciarlet The Finite Element Method for Elliptic Problems North-Holland, New York, 1978 MR 58:25001 [7] H Elman Preconditioning for the steady-state Navier-Stokes equations with low viscosity Technical Report CS-TR-3712, Department of Computer Science, University of Maryland, College Park, MD 20742, 1996 [8] H Elman and D Silvester Fast nonsymmetric iterations and preconditioning for Navier- Stokes equations SIAM J Sci Comput, 17:33 46, 1996 MR 97e:65119 [9] V Girault and PA Raviart Finite Element Approximation of the Navier-Stokes Equations Lecture Notes in Math # 749, Springer-Verlag, New York, 1981 MR 83b:65122 [10] C Johnson and J Pitkäranta Analysis of some mixed finite element methods related to reduced integration Math Comp, 38: , 1982 MR 83d:65287 [11] OA Karakashian On a Galerkin-Lagrange multiplier method for the stationary Navier- Stokes equations SIAM J Numer Anal, 19: , 1982 MR 83j:65107 [12] C Lanczos An iteration method for the solution of the eigenvalue problem of linear differential and integral operators J Res National Bureau of Standards, 45: , 1950 MR 13:163d [13] JL Lions and E Magenes Problèmes aux Limites non Homogènes et Applications, volume1 Dunod, Paris, 1968 MR 40:512 [14] MM Liu, J Wang, and N-N Yan New error estimates for approximate solutions of convection-diffusion problems by mixed and discontinuous Galerkin methods SIAM J Numer Anal Submitted [15] MF Murphy and AJ Wathen On preconditioning for the Oseen equations Technical Report AM 95-07, Department of Mathematics, University of Bristol, 1995 [16] Y Saad and MH Schultz GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems SIAM J Sci Stat Comput, 7: , 1986 MR 87g:65064 [17] R Temam Navier-Stokes Equations North-Holland Publishing Co, New York, 1977 MR 58:29439 Department of Mathematics, Texas A&M University, College Station, Texas Schlumberger, 8311 N FM 620, Austin, Texas 78726
AMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationSchur Complements on Hilbert Spaces and Saddle Point Systems
Schur Complements on Hilbert Spaces and Saddle Point Systems Constantin Bacuta Mathematical Sciences, University of Delaware, 5 Ewing Hall 976 Abstract For any continuous bilinear form defined on a pair
More informationA Review of Preconditioning Techniques for Steady Incompressible Flow
Zeist 2009 p. 1/43 A Review of Preconditioning Techniques for Steady Incompressible Flow David Silvester School of Mathematics University of Manchester Zeist 2009 p. 2/43 PDEs Review : 1984 2005 Update
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationA Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations
A Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations Songul Kaya and Béatrice Rivière Abstract We formulate a subgrid eddy viscosity method for solving the steady-state
More informationA NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION
A NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION JOHNNY GUZMÁN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS Abstract. The analysis of finite-element-like Galerkin discretization techniques for the
More informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationA Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers
Applied and Computational Mathematics 2017; 6(4): 202-207 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20170604.18 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Robust Preconditioned
More informationFast Iterative Solution of Saddle Point Problems
Michele Benzi Department of Mathematics and Computer Science Emory University Atlanta, GA Acknowledgments NSF (Computational Mathematics) Maxim Olshanskii (Mech-Math, Moscow State U.) Zhen Wang (PhD student,
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationJournal of Computational and Applied Mathematics. Multigrid method for solving convection-diffusion problems with dominant convection
Journal of Computational and Applied Mathematics 226 (2009) 77 83 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
More informationMULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
More information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More informationMultigrid Methods for Saddle Point Problems
Multigrid Methods for Saddle Point Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University Advances in Mathematics of Finite Elements (In
More informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
More informationChapter 2 Finite Element Spaces for Linear Saddle Point Problems
Chapter 2 Finite Element Spaces for Linear Saddle Point Problems Remark 2.1. Motivation. This chapter deals with the first difficulty inherent to the incompressible Navier Stokes equations, see Remark
More informationFast solvers for steady incompressible flow
ICFD 25 p.1/21 Fast solvers for steady incompressible flow Andy Wathen Oxford University wathen@comlab.ox.ac.uk http://web.comlab.ox.ac.uk/~wathen/ Joint work with: Howard Elman (University of Maryland,
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationGeometric Multigrid Methods
Geometric Multigrid Methods Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University IMA Tutorial: Fast Solution Techniques November 28, 2010 Ideas
More informationSEMI-CONVERGENCE ANALYSIS OF THE INEXACT UZAWA METHOD FOR SINGULAR SADDLE POINT PROBLEMS
REVISTA DE LA UNIÓN MATEMÁTICA ARGENTINA Vol. 53, No. 1, 2012, 61 70 SEMI-CONVERGENCE ANALYSIS OF THE INEXACT UZAWA METHOD FOR SINGULAR SADDLE POINT PROBLEMS JIAN-LEI LI AND TING-ZHU HUANG Abstract. Recently,
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationPreconditioners for the incompressible Navier Stokes equations
Preconditioners for the incompressible Navier Stokes equations C. Vuik M. ur Rehman A. Segal Delft Institute of Applied Mathematics, TU Delft, The Netherlands SIAM Conference on Computational Science and
More information7.4 The Saddle Point Stokes Problem
346 CHAPTER 7. APPLIED FOURIER ANALYSIS 7.4 The Saddle Point Stokes Problem So far the matrix C has been diagonal no trouble to invert. This section jumps to a fluid flow problem that is still linear (simpler
More informationEfficient Solvers for the Navier Stokes Equations in Rotation Form
Efficient Solvers for the Navier Stokes Equations in Rotation Form Computer Research Institute Seminar Purdue University March 4, 2005 Michele Benzi Emory University Atlanta, GA Thanks to: NSF (MPS/Computational
More informationNumerische Mathematik
Numer. Math. (00) 93: 333 359 Digital Object Identifier (DOI) 10.1007/s00110100386 Numerische Mathematik Two new variants of nonlinear inexact Uzawa algorithms for saddle-point problems Qiya Hu 1,, Jun
More informationDomain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions
Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions Ernst P. Stephan 1, Matthias Maischak 2, and Thanh Tran 3 1 Institut für Angewandte Mathematik, Leibniz
More informationA STOKES INTERFACE PROBLEM: STABILITY, FINITE ELEMENT ANALYSIS AND A ROBUST SOLVER
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, K. Majava, and O. Pironneau (eds.) O. Nevanlinna and R. Rannacher (assoc. eds.) Jyväskylä,
More information20. A Dual-Primal FETI Method for solving Stokes/Navier-Stokes Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 23 DDM.org 2. A Dual-Primal FEI Method for solving Stokes/Navier-Stokes
More informationA mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N)
wwwijmercom Vol2, Issue1, Jan-Feb 2012 pp-464-472 ISSN: 2249-6645 A mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N) Jaouad El-Mekkaoui 1, Abdeslam Elakkad
More informationOn convergence rate of the augmented Lagrangian algorithm for nonsymmetric saddle point problems
Applied Numerical athematics 54 2005) 122 134 wwwelseviercom/locate/apnum On convergence rate of the augmented Lagrangian algorithm for nonsymmetric saddle point problems G Awanou a,,jlai b a Institute
More informationITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS
ITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS LONG CHEN In this chapter we discuss iterative methods for solving the finite element discretization of semi-linear elliptic equations of the form: find
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More informationKey words. preconditioned conjugate gradient method, saddle point problems, optimal control of PDEs, control and state constraints, multigrid method
PRECONDITIONED CONJUGATE GRADIENT METHOD FOR OPTIMAL CONTROL PROBLEMS WITH CONTROL AND STATE CONSTRAINTS ROLAND HERZOG AND EKKEHARD SACHS Abstract. Optimality systems and their linearizations arising in
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
More informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationA LEAST-SQUARES APPROACH BASED ON A DISCRETE MINUS ONE INNER PRODUCT FOR FIRST ORDER SYSTEMS
MATHEMATICS OF COMPUTATION Volume 66, Number 219, July 1997, Pages 935 955 S 0025-5718(97)00848-X A LEAST-SQUARES APPROACH BASED ON A DISCRETE MINUS ONE INNER PRODUCT FOR FIRST ORDER SYSTEMS JAMES H BRAMBLE,
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationMathematics and Computer Science
Technical Report TR-2007-002 Block preconditioning for saddle point systems with indefinite (1,1) block by Michele Benzi, Jia Liu Mathematics and Computer Science EMORY UNIVERSITY International Journal
More informationOn preconditioned Uzawa-type iterations for a saddle point problem with inequality constraints
On preconditioned Uzawa-type iterations for a saddle point problem with inequality constraints Carsten Gräser and Ralf Kornhuber FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationINTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR GRIDS. 1. Introduction
Trends in Mathematics Information Center for Mathematical Sciences Volume 9 Number 2 December 2006 Pages 0 INTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR
More informationIterative Methods for Solving A x = b
Iterative Methods for Solving A x = b A good (free) online source for iterative methods for solving A x = b is given in the description of a set of iterative solvers called templates found at netlib: http
More informationTermination criteria for inexact fixed point methods
Termination criteria for inexact fixed point methods Philipp Birken 1 October 1, 2013 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany Department of Mathematics/Computer
More informationPRECONDITIONING TECHNIQUES FOR NEWTON S METHOD FOR THE INCOMPRESSIBLE NAVIER STOKES EQUATIONS
BIT Numerical Mathematics 43: 961 974, 2003. 2003 Kluwer Academic Publishers. Printed in the Netherlands. 961 PRECONDITIONING TECHNIQUES FOR NEWTON S METHOD FOR THE INCOMPRESSIBLE NAVIER STOKES EQUATIONS
More informationXIAO-CHUAN CAI AND MAKSYMILIAN DRYJA. strongly elliptic equations discretized by the nite element methods.
Contemporary Mathematics Volume 00, 0000 Domain Decomposition Methods for Monotone Nonlinear Elliptic Problems XIAO-CHUAN CAI AND MAKSYMILIAN DRYJA Abstract. In this paper, we study several overlapping
More informationA NEW FAMILY OF STABLE MIXED FINITE ELEMENTS FOR THE 3D STOKES EQUATIONS
MATHEMATICS OF COMPUTATION Volume 74, Number 250, Pages 543 554 S 0025-5718(04)01711-9 Article electronically published on August 31, 2004 A NEW FAMILY OF STABLE MIXED FINITE ELEMENTS FOR THE 3D STOKES
More informationANALYSIS OF A FINITE ELEMENT PML APPROXIMATION FOR THE THREE DIMENSIONAL TIME-HARMONIC MAXWELL PROBLEM
MATHEMATICS OF COMPUTATION Volume 77, Number 261, January 2008, Pages 1 10 S 0025-5718(07)02037-6 Article electronically published on September 18, 2007 ANALYSIS OF A FINITE ELEMENT PML APPROXIMATION FOR
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
More informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationLocal flux mimetic finite difference methods
Local flux mimetic finite difference methods Konstantin Lipnikov Mikhail Shashkov Ivan Yotov November 4, 2005 Abstract We develop a local flux mimetic finite difference method for second order elliptic
More informationON A GENERAL CLASS OF PRECONDITIONERS FOR NONSYMMETRIC GENERALIZED SADDLE POINT PROBLEMS
U..B. Sci. Bull., Series A, Vol. 78, Iss. 4, 06 ISSN 3-707 ON A GENERAL CLASS OF RECONDIIONERS FOR NONSYMMERIC GENERALIZED SADDLE OIN ROBLE Fatemeh anjeh Ali BEIK his paper deals with applying a class
More informationUniversity of Maryland Department of Computer Science TR-5009 University of Maryland Institute for Advanced Computer Studies TR April 2012
University of Maryland Department of Computer Science TR-5009 University of Maryland Institute for Advanced Computer Studies TR-202-07 April 202 LYAPUNOV INVERSE ITERATION FOR COMPUTING A FEW RIGHTMOST
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationKey words. inf-sup constant, iterative solvers, preconditioning, saddle point problems
NATURAL PRECONDITIONING AND ITERATIVE METHODS FOR SADDLE POINT SYSTEMS JENNIFER PESTANA AND ANDREW J. WATHEN Abstract. The solution of quadratic or locally quadratic extremum problems subject to linear(ized)
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationGlobal Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations
Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations C. David Levermore Department of Mathematics and Institute for Physical Science and Technology
More informationA Finite Element Method for the Surface Stokes Problem
J A N U A R Y 2 0 1 8 P R E P R I N T 4 7 5 A Finite Element Method for the Surface Stokes Problem Maxim A. Olshanskii *, Annalisa Quaini, Arnold Reusken and Vladimir Yushutin Institut für Geometrie und
More informationFrom Completing the Squares and Orthogonal Projection to Finite Element Methods
From Completing the Squares and Orthogonal Projection to Finite Element Methods Mo MU Background In scientific computing, it is important to start with an appropriate model in order to design effective
More informationMODEL OF GROUNDWATER FLOW IN FRACTURED ENVIRONMENT
Proceedings of ALGORITMY 2002 Conference on Scientific Computing, pp. 138 145 MODEL OF GROUNDWATER FLOW IN FRACTURED ENVIRONMENT JIŘÍ MARYŠKA, OTTO SEVERÝN AND MARTIN VOHRALÍK Abstract. A stochastic discrete
More informationFIRST-ORDER SYSTEM LEAST SQUARES (FOSLS) FOR PLANAR LINEAR ELASTICITY: PURE TRACTION PROBLEM
SIAM J. NUMER. ANAL. c 1998 Society for Industrial Applied Mathematics Vol. 35, No. 1, pp. 320 335, February 1998 016 FIRST-ORDER SYSTEM LEAST SQUARES (FOSLS) FOR PLANAR LINEAR ELASTICITY: PURE TRACTION
More informationEfficient smoothers for all-at-once multigrid methods for Poisson and Stokes control problems
Efficient smoothers for all-at-once multigrid methods for Poisson and Stoes control problems Stefan Taacs stefan.taacs@numa.uni-linz.ac.at, WWW home page: http://www.numa.uni-linz.ac.at/~stefant/j3362/
More informationError estimates for the Raviart-Thomas interpolation under the maximum angle condition
Error estimates for the Raviart-Thomas interpolation under the maximum angle condition Ricardo G. Durán and Ariel L. Lombardi Abstract. The classical error analysis for the Raviart-Thomas interpolation
More informationSolving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners
Solving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University of Minnesota 2 Department
More informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 16-02 The Induced Dimension Reduction method applied to convection-diffusion-reaction problems R. Astudillo and M. B. van Gijzen ISSN 1389-6520 Reports of the Delft
More informationNew Discretizations of Turbulent Flow Problems
New Discretizations of Turbulent Flow Problems Carolina Cardoso Manica and Songul Kaya Merdan Abstract A suitable discretization for the Zeroth Order Model in Large Eddy Simulation of turbulent flows is
More informationNumerische Mathematik
Numer. Math. (2003) 94: 195 202 Digital Object Identifier (DOI) 10.1007/s002110100308 Numerische Mathematik Some observations on Babuška and Brezzi theories Jinchao Xu, Ludmil Zikatanov Department of Mathematics,
More informationFINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
More informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationA Two-grid Method for Coupled Free Flow with Porous Media Flow
A Two-grid Method for Coupled Free Flow with Porous Media Flow Prince Chidyagwai a and Béatrice Rivière a, a Department of Computational and Applied Mathematics, Rice University, 600 Main Street, Houston,
More informationResearch Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations
Abstract and Applied Analysis Volume 212, Article ID 391918, 11 pages doi:1.1155/212/391918 Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Chuanjun Chen
More informationA Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract
A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin
More informationMEYERS TYPE ESTIMATES FOR APPROXIMATE SOLUTIONS OF NONLINEAR ELLIPTIC EQUATIONS AND THEIR APPLICATIONS. Yalchin Efendiev.
Manuscript submitted to AIMS Journals Volume X, Number 0X, XX 200X Website: http://aimsciences.org pp. X XX MEYERS TYPE ESTIMATES FOR APPROXIMATE SOLUTIONS OF NONLINEAR ELLIPTIC EUATIONS AND THEIR APPLICATIONS
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationNEW ESTIMATES FOR RITZ VECTORS
MATHEMATICS OF COMPUTATION Volume 66, Number 219, July 1997, Pages 985 995 S 0025-5718(97)00855-7 NEW ESTIMATES FOR RITZ VECTORS ANDREW V. KNYAZEV Abstract. The following estimate for the Rayleigh Ritz
More information2 A Model, Harmonic Map, Problem
ELLIPTIC SYSTEMS JOHN E. HUTCHINSON Department of Mathematics School of Mathematical Sciences, A.N.U. 1 Introduction Elliptic equations model the behaviour of scalar quantities u, such as temperature or
More informationConvergence of a Class of Stationary Iterative Methods for Saddle Point Problems
Convergence of a Class of Stationary Iterative Methods for Saddle Point Problems Yin Zhang 张寅 August, 2010 Abstract A unified convergence result is derived for an entire class of stationary iterative methods
More informationAn Equal-order DG Method for the Incompressible Navier-Stokes Equations
An Equal-order DG Method for the Incompressible Navier-Stokes Equations Bernardo Cockburn Guido anschat Dominik Schötzau Journal of Scientific Computing, vol. 40, pp. 188 10, 009 Abstract We introduce
More informationNumerische Mathematik
Numer. Math. 000) 86: 443 470 Digital Object Identifier DOI) 0.007/s00000060 Numerische Mathematik Effective preconditioning of Uzawa type schemes for a generalized Stokes problem Georgij M. Kobelkov,
More informationNumerische Mathematik
Numer. Math. (2002) 90: 665 688 Digital Object Identifier (DOI) 10.1007/s002110100300 Numerische Mathematik Performance and analysis of saddle point preconditioners for the discrete steady-state Navier-Stokes
More informationSTOPPING CRITERIA FOR MIXED FINITE ELEMENT PROBLEMS
STOPPING CRITERIA FOR MIXED FINITE ELEMENT PROBLEMS M. ARIOLI 1,3 AND D. LOGHIN 2 Abstract. We study stopping criteria that are suitable in the solution by Krylov space based methods of linear and non
More informationVARIATIONAL AND NON-VARIATIONAL MULTIGRID ALGORITHMS FOR THE LAPLACE-BELTRAMI OPERATOR.
VARIATIONAL AND NON-VARIATIONAL MULTIGRID ALGORITHMS FOR THE LAPLACE-BELTRAMI OPERATOR. ANDREA BONITO AND JOSEPH E. PASCIAK Abstract. We design and analyze variational and non-variational multigrid algorithms
More informationOn the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations
On the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations Lutz Tobiska Institut für Analysis und Numerik Otto-von-Guericke-Universität
More informationA Note on Inverse Iteration
A Note on Inverse Iteration Klaus Neymeyr Universität Rostock, Fachbereich Mathematik, Universitätsplatz 1, 18051 Rostock, Germany; SUMMARY Inverse iteration, if applied to a symmetric positive definite
More informationAdaptive algebraic multigrid methods in lattice computations
Adaptive algebraic multigrid methods in lattice computations Karsten Kahl Bergische Universität Wuppertal January 8, 2009 Acknowledgements Matthias Bolten, University of Wuppertal Achi Brandt, Weizmann
More informationFast Solvers for Unsteady Incompressible Flow
ICFD25 Workshop p. 1/39 Fast Solvers for Unsteady Incompressible Flow David Silvester University of Manchester http://www.maths.manchester.ac.uk/~djs/ ICFD25 Workshop p. 2/39 Outline Semi-Discrete Navier-Stokes
More informationMultigrid absolute value preconditioning
Multigrid absolute value preconditioning Eugene Vecharynski 1 Andrew Knyazev 2 (speaker) 1 Department of Computer Science and Engineering University of Minnesota 2 Department of Mathematical and Statistical
More informationCourant Institute of Mathematical Sciences, New York University, New York, NY 10012,
EFFICIENT VARIABLE-COEFFICIENT FINITE-VOLUME STOKES SOLVERS MINGCHAO CAI, ANDY NONAKA, JOHN B. BELL, BOYCE E. GRIFFITH, AND ALEKSANDAR DONEV Abstract. We investigate several robust preconditioners for
More informationM.A. Botchev. September 5, 2014
Rome-Moscow school of Matrix Methods and Applied Linear Algebra 2014 A short introduction to Krylov subspaces for linear systems, matrix functions and inexact Newton methods. Plan and exercises. M.A. Botchev
More informationPreconditioners for reduced saddle point systems arising in elliptic PDE-constrained optimization problems
Zeng et al. Journal of Inequalities and Applications 205 205:355 DOI 0.86/s3660-05-0879-x RESEARCH Open Access Preconditioners for reduced saddle point systems arising in elliptic PDE-constrained optimization
More informationASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR THE POINTWISE GRADIENT ERROR ON EACH ELEMENT IN IRREGULAR MESHES. PART II: THE PIECEWISE LINEAR CASE
MATEMATICS OF COMPUTATION Volume 73, Number 246, Pages 517 523 S 0025-5718(0301570-9 Article electronically published on June 17, 2003 ASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR TE POINTWISE GRADIENT
More information