STOPPING CRITERIA FOR MIXED FINITE ELEMENT PROBLEMS
|
|
- Cameron Booker
- 5 years ago
- Views:
Transcription
1 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 linear systems of equations arising from the mixed and the mixed-hybrid finiteelement approximation of saddle point problems. Our approach is based on the equivalence between the Babuška and Brezzi conditions of stability which allows us to apply some of the results obtained in Arioli, Loghin and Wathen [1]. Our proposed criterion involves evaluating the residual in a norm defined on the discrete dual of the space where we seek a solution. We illustrate our approach using standard iterative methods such as MINRES and GMRES. We test our criteria on Stokes and Navier-Stokes problems both in a linear and nonlinear context. Key words. augmented systems, mixed and mixed-hybrid finite-element, stopping criteria, Krylov subspaces method AMS subject classifications. 65F10, 64F35, 65F50, 65N30 1. Introduction. Mixed and mixed-hybrid finite-element methods form a class of popular discretization methods designed to approximate systems of partial differential equations of saddle-point type arising in the modeling of a variety of physical phenomena in areas such as fluid-dynamics or linear elasticity. They generally give rise to large, nonsymmetric, indefinite linear and nonlinear systems for which the solution is typically sought via iterative approaches. An essential feature of such methods is the stopping criteria employed. This work aims to describe how to devise a suitable stopping procedure, given the well-defined theoretical context of variational formulation of partial differential equations and in particular the mixed finite element theory. The outline of the paper is as follows. In Section 2, we introduce the abstract formulation of a generic saddle-point problem as a system based on bilinear forms. Then, we describe a general framework in which we can formulate a stopping criterion based on the energy norm of the error between the exact solution of the continuous problem and the solution computed by an iterative method. Section 3 generalizes the stopping criterion derived in [1] to the case of mixed finite element formulations, discussing both the linear symmetric and nonsymmetric cases. We also propose a strategy for the extension of the stopping criteria to the nonlinear case. Finally, in Section 4, we present our class of test problems together with the convergence behaviour of some iterative algorithms showing the beneficial effect of our stopping criteria. 2. Mixed variational formulation. We start by summarizing the theoretical setting necessary to describe our problem. A comprehensive and exhaustive introduction can be found in the book of Brezzi and Fortin [6]. Let V, Q be Hilbert spaces with norms V, Q and duals V, Q, respectively. Consider the two real-valued bilinear forms a(, ) : V V, b(, ) : V Q and the two linear functionals f( ) V, g( ) Q. We are interested in the following abstract 1 Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire, OX11 0QX, UK. m.arioli@rl.ac.uk,phone: University of Birmingham, Watson Building, Birmingham B15 2TT, UK. loghind@for.mat.bham.ac.uk. Phone: The work of the first author was supported in part by EPSRC grants GR/R46641/01 and GR/S
2 2 M. ARIOLI and D. LOGHIN variational formulation Find (u, p) V Q such that for all (v, q) V Q (SP) a(u, v) + b(v, p) = f(v), b(u, q) = g(q). In the nonlinear case the bilinear form a(, ) is replaced by the nonlinear operator F : V V, as, for example, in the Navier-Stokes case. The variational formulation in this case reads Find (u, p) V Q such that for all (v, q) V Q (N SP) F (u), v (V,V) + b(v, p) = f(v) b(u, q) = g(q). Following [16], [12], and [26], we introduce the Hilbert space H = V Q with the norm graph: { H w = { u, q } w 2 H = v 2 V + q 2 Q, the bilinear form K : H H IR and the linear functional f : H IR, f H : K(u, p; v, q) = a(u, v) + b(v, p) + b(u, q), f(u, q) = f(v) + g(q), (2.1) where we equip H with the norm H given by f 2 H = f 2 V + g 2 Q. Problem SP can be reformulated as { Find u H such that for all v H K(u, v) = f(v). (2.2) Existence and uniqueness of solutions to problems of type (2.2) is guaranteed provided the following conditions hold for all u, v H sup v H\{0} K(w, v) C 1 w H v H (2.3a) K(w, v) C 2 w H, v H (2.3b) K(w, v) sup C 2 v H ; (2.3c) w H\{0} w H for some positive constants C 1, C 2. Remark 2.1. Requirements (2.3) are known as the Babuška conditions and can be shown to be equivalent to the Brezzi conditions which essentially are (i) continuity conditions (of type (2.3a)) on a(, ), b(, ), (ii) a condition of type (2.3b) for b(, ) and (iii) a coercivity condition on a(, ) [26, 12]. In the following, we find it convenient to work with the Babuška conditions. Consider now the finite dimensional spaces V h V and Q h Q with bases { ψ i }1 i n and { φ j }1 j m, respectively. Moreover, we denote by H h and its dual Hh the spaces H h = V h Q h, H h = V h Q h.
3 STOPPING CRITERIA FOR MIXED FINITE-ELEMENT PROBLEMS 3 Variational formulation (2.2) restricted to the finite dimensional space H h reads { Find uh H h such that for all v h H h (2.4) K h (u h, v h ) = f h (v h ). where K h (, ) is a bilinear form on H h H h and f h ( ) is a continuous linear form on H h. In the following, we assume that the Babuška conditions (2.3) hold for the bilinear form K h (, ). This allows us to derive the a priori error estimate u u h H ( 1 + C 1 C 2 ) min u v h H. (2.5) v h V h Remark 2.2. We shall be assuming that the variational formulations introduced above are weak formulations of a system of partial differential equations defined on some open subset Ω of IR d. Then the Hilbert spaces are spaces of real-valued functions defined on Ω, while V h, Q h are finite element spaces, spanned by basis functions defined on a subdivision Ω h of Ω. Replacing v h by the interpolant of u on Ω h and using standard interpolation error estimates, we can derive a priori bounds of the form u u h H C(u)C(h), which are very useful in informing our approach to designing stopping criteria. For the choice (2.1), the weak formulation (2.4) gives rise to a linear system of equations Ku = f, where the matrix K has the 2-by-2 block structure ( ) A B T K = B 0 with A ij = a(ψ j, ψ i ), B kj = b(ψ j, φ k ), i, j = 1 n, k = 1 m. Let us examine the discrete setting further. Note first that there is an isomorphism Π h between IR n+m and H h defined via ( ) ( n v Π h w = Π h = i=1 v ) ( ) m iψ i vh q j=1 q = = w jφ j q h. h In particular, since v h 2 V h = v T V v = v 2 V, q h 2 Q h = q t Qq = q 2 Q, where V IR n n and Q IR m m, the finite dimensional Hilbert spaces (V h, Vh ), (Q h, Qh ) are represented, respectively, by (IR n, V ), (IR m, Q ). Therefore, the space H h can be represented by IR n+m with norm H where H IR (n+m) (n+m) is given by H = [ V 0 0 Q ].
4 4 M. ARIOLI and D. LOGHIN The dual space H h can be shown to be represented by IRn+m with norm H 1. Finally, we have the following discrete representation K h (u h, v h ) = v T Ku u h, v h H h, which allows us to write the continuous stability conditions (2.3) as max max w R n \{0} v R n \{0} min max w R n \{0} v R n \{0} w T Kv w H v H C 1 w T Kv w H v H C 2 (2.6a) (2.6b) which is equivalent to uniform conditioning of K with respect to the norm induced by H: K H,H 1 C 1, K 1 H 1,H C 1 2, or, κ H (K) C 1 /C 2. We point out that both C 1 and C 2 are constants independent of h and, thus, independent of n and m. 3. Stopping criteria. Conditions (2.6) are sufficient for the main theorem in [1] to apply: Theorem 3.1. Let u be the solution of the weak formulation (2.2) and let u, u h = Π h u satisfy Ku = f; u u h H u h H C(h). Then ũ h = Π h ũ satisfies u ũ h H ũ h H C(h) = O(C(h)) if f Kũ H 1 ũ H ηc(h)c 2, (3.1) for some η (0, 1). Remark 3.1. This result means that one may replace the finite element solution u h by an approximation ũ h constructed by an iterative method provided the H 1 -norm of the residual r = f Kũ is of the same order as the finite element error. In the following, we consider in greater detail the application of the above criterion to saddle-point systems, both in a linear and nonlinear setting The linear case. Unlike the positive-definite case considered in [1], there is no obvious solution, or iterative method, that would allow for the approximation of r H 1 in an indefinite context. In fact, it appears that this may have to be computed by solving a linear system with coefficient matrix H. Fortunately, this is a procedure that is included already in some preconditioned iterative methods.
5 STOPPING CRITERIA FOR MIXED FINITE-ELEMENT PROBLEMS Symmetric indefinite problems. It is an established fact that symmetric saddle-point problems arising from the stable finite element discretization of a system of partial differential equations are rather amenable to iterative treatment in the sense that they come equipped with optimal preconditioners. We quote here a general result from [18] that expresses this fact. Theorem 3.2. Let (2.6) hold. Then H 1 K H = KH 1 H 1 C 1, (3.2a) K 1 H H = HK 1 H 1 C2 1 (3.2b) While the form of (3.2) is useful when we consider the nonsymmetric case, we note here that one can write the above bounds as a bound on the 2-norm condition number of K preconditioned centrally by the norm κ 2 (H 1/2 KH 1/2 ) C 1 C 2. This suggests that an iterative method such as the Minimum Residual method (MIN- RES) will converge in a number of steps independent of the size of the problem. Furthermore, the residual computed by this method is in fact measured in the right norm: H 1. Hence, one can easily incorporate in this approach bound (3.1). This we carry out in our numerics section. We note here that there is a significant amount of research devoted to the analysis of norm-based preconditioners for symmetric saddle-point problems and derivation of bounds of type (3.2). Some of the problems considered come from ground-water flow applications [4, 5, 10, 15, 22, 24], Stokes flow [8, 25, 23, 13, 9], elasticity [14, 2, 17, 7, 19], magnetostatics [20, 21] etc Nonsymmetric indefinite problems. While convergence of iterative methods for nonsymmetric problems is not fully understood, bounds such as (3.2) are clearly attractive in a preconditioning context. They guarantee that for H 1/2 KH 1/2 both the singular values and the absolute values of the eigenvalues are bounded from below and above. This means that the use of the norm as a preconditioner can be recommended also in the nonsymmetric case. The general approach, as suggested by the form of (3.2) is to employ an iterative solver in the H-inner product with left preconditioner H. The resulting algorithm is equivalent to employing an Euclidean inner-product and system matrix H 1/2 KH 1/2 and output a residual measured in the norm H 1 which is what we want to monitor. We carry out this kind of procedure in the case of the Generalized Minimum Residual method (GMRES) The nonlinear case. In the nonlinear case the approximation of the solutions by mixed or mixed-hybrid methods in combination with the linearization of the operator by a Newton method or a Picard approach, yields a sequence of finite dimensional problems of type (2.4), generally nonsymmetric, each of which satisfy the stability conditions (2.6). Writing the approximation of problem N SP as after linearization, we want to solve F(w) = 0 K k w k+1 = g k, (3.3)
6 6 M. ARIOLI and D. LOGHIN where w k+1 = ( ) uk+1 p k+1 ( fk and g k = 0 ). In practice, an inner-outer iteration process is the method of choice for large problems, with an inner linear solve, typically an iterative process, and an outer nonlinear update. A popular approach is the Newton-Krylov procedure, where the outer Newton iteration uses an inner iterative procedure of Krylov type to solve the linear system (3.3). The efficiency of such methods relies on the choice of inner stopping criteria. This was recognized by Dembo et al. [11]. We review briefly their result and adapt it to the finite element context Nonlinear stopping criteria. Let us assume that we seek the solution of problem (3.3) via an iterative routine in which at every step j we compute or estimate a residual r j = g k K k w k+1. Since initially system (3.3) approximates poorly the equation F(w) = 0 one can compute only a coarse approximation to the solution w k+1. As convergence of the outer iteration improves one would want to improve the quality of w k+1. One way of achieving this is through the use of of the following inner stopping criterion [11] r j F(w k ) c F(w k) q. (3.4) The norm employed in (3.4) is general, though in practice the standard Euclidean norm is employed. This is wasteful in a finite element context. In particular, given the result of Thm 3.1, we propose to evaluate the above criterion in the relevant norm H 1 r j H 1 F(w k ) H 1 c F(w k ) q H 1, (3.5) Criterion (3.5) is to be combined with criterion (3.1). Thus, while (3.1) is not satisfied, one employs at each nonlinear step an iterative method with criterion (3.5). Moreover, the choice of c, q in (3.5) needs to be related to C(h) in (3.1). Thus, if the problem is large, one can compute a satisfactory solution in just a few iterations, without the need to attain a very small order for the nonlinear residual, which in a finite element context has no relevant meaning. A typical algorithm for solving (3.3) is outlined below: k = 0, choose w k, r k = F k (w k ) tol out := ηc(h)c 2 while r k H 1/ w k H tol out end while w 0 = w n, r 0 = r k, tol in := c r 0 q H 1 w j = GMRES(K k, g k, w 0, tol in, H) k = k + 1, w k = w j, r k = F(w k )
7 STOPPING CRITERIA FOR MIXED FINITE-ELEMENT PROBLEMS 7 where the iterative routine GMRES(E, b, x 0, tol, H) applied to a matrix E computes an approximate solution x k such that b Ex k H 1 b Ex 0 H 1 tol. A GMRES routine which uses the H 1 -norm in its stopping criterion is a GMRES iteration in the H-inner product preconditioned from the left by H (cf. section 3.1.2) term GMRES. It was shown in the positive-definite case in [1] that the GMRES method in the H-inner product with left preconditioner H is a three-term recurrence provided H is the symmetric part of K. In this case, the preconditioned system matrix is a normal matrix H 1/2 KH 1/2 = I + S where S is a skew-symmetric matrix. Such an implementation of GMRES is storagefree, a desirable feature in an iterative solver. One would naturally want to extend this to the indefinite case, particularly for the case where we have to solve a long sequence of problems of type (3.3). We show how this can be achieved for a class of nonlinear saddle-point problems. Let K k in (3.3) have the form K k = ( Ak B T k B k 0 where A k are nonsymmetric positive-definite for all k, a standard assumption for a great variety of problems. Let us replace the sequence of problems (3.3) with the following sequence ( ) ( ) ( ) Ak Bk T uk+1 f = k. (3.6) τm p k+1 τmp k B k It is easy to see that this sequence converges to the same solution provided τ is sufficiently small (in fact, we require τ ρ(m 1 BA 1 k BT )). Multiplying the second set of equations by minus one, equation (3.6) becomes ( ) ( ) Ak Bk T uk+1 = B k τm p k+1 ) ( fk τmp k ) (3.7) and thus one can split the system matrix into a symmetric (positive-definite) and anti-symmetric part ( ) ( ) ( ) Ak Bk T Lk 0 Sk B = + k T. B k τm 0 τm B k 0 It is clear now that the three-term GMRES method devised for the scalar case [1] will work also in this case, provided we precondition the above matrix centrally with the hermitian part H k of the modified matrix K k H k = ( Lk 0 0 τm ).
8 8 M. ARIOLI and D. LOGHIN Residual convergence will then be automatically measured in the norm v 2 H 1 k = v 1 2 L 1 k + τ 1 v 2 2 M 1. It is clear that M will be chosen to be Q, while L k will be in general equivalent to V, when not identically equal to it. Thus, during the Arnoldi process, one can monitor strictly the H 1 -norm of the residual. Numerical experiments indicate that the method does not affect the convergence rate of the outer iteration for τ sufficiently small, the advantage being that one can employ a short-term recurrence for the inner iteration. 4. Experiments Test problems. We used two test problems suggested in [3]. The first is Stokes flow in the unit square while the second is 2D Navier-Stokes flow in a cavity, which tries to mimic the behaviour of the driven-cavity flow. The problems we solved are and u + p = f in Ω (4.1a) div u = 0 in Ω (4.1b) u(x) = u (x) on Γ, (4.1c) ε u + ( u ) u + p = f in Ω (4.2a) div u = 0 in Ω (4.2b) u(x) = u (x) on Γ, (4.2c) both of which have the exact solution ( u, p) = (u, v, p ) given by where u (x, y) = R 2 2π q 0(R 2, y) (1 cos 2πq(R 1, x)) sin 2πq(R 2, y) v (x, y) = R 1 2π q 0(R 1, x) (1 cos 2πq(R 2, y)) sin 2πq(R 1, x) p (x, y) = R 1 R 2 q 0 (R 1, x)q 0 (R 2, y) sin 2πq(R 1, x) sin 2πq(R 2, y). q(r, t) = ert 1 e R 1, q 0(R, t) = ert e R 1 and R 1, R 2 are two real constants that can be used to modify the flow behaviour. The pressure satisfies p dx = 0, (4.3) Ω and this is the type of condition one can use to ensure that equations (4.2) have a unique solution. Streamlines and pressure plots for various values of R 1, R 2 are given in Figs 4.1, 4.2. The problem used in our tests corresponds to the choice R 1 = 0.1, R 2 = 4.2. We solved (4.1), (4.2) using the discrete mixed formulation (2.4) where the space H h H with H = H 1 0 (Ω) L 2 (Ω). In particular, we chose to work with the norm v 2 H = ε v q 2 L 2 (Ω)
9 STOPPING CRITERIA FOR MIXED FINITE-ELEMENT PROBLEMS (a) Streamlines (b) Pressure Fig Exact solution for R 1 = 0.1, R 2 = (a) Streamlines (b) Pressure Fig Exact solution for R 1 = 1.2, R 2 = 0.1. where v = ( v, q) and v(x) H 1 0 (Ω) = v(x) 1 = α =1 Ω D α v(x) 2 dx Our discrete spaces V h, Q h were finite element spaces spanned by quadratic basis functions in the case of the velocity space and linear basis functions in the case of the pressure space The linear case. We chose to compare the stopping criterion (3.1) with both the exact finite element error and interpolation error measured in the norms inherited by the problem. We computed (i) : the exact relative errors between the solution at step k and the exact 1/2.
10 10 M. ARIOLI and D. LOGHIN continuous solution of either (4.1) or (4.2) F E := u uk h H u k h ; H (ii) FIE: the exact relative interpolation errors F IE := ui u k h H u k h ; H (iii) : the exact H 1 -norm criterion (3.1) with C 2 estimated on a coarse mesh; (iv) the standard 2-norm stopping criterion r k / r 0. We first display in Fig. 4.3 the results for MINRES preconditioned with the norm in the case of the Stokes problem. In Fig. 4.3 (a) we plot the value of the global error while in Fig. 4.3 (b), (c), and (d), we plot the values of, FIE,, and 2-norm of ε=1; N=2113; solver=minres ε=1; N=2113; solver=minres uvp u (a) Global (uvp) convergence (b) Component u convergence ε=1; N=2113; solver=minres ε=1; N=2113; solver=minres v p (c) Component v (d) Component p Fig Convergence criteria for preconditioned MINRES.
11 STOPPING CRITERIA FOR MIXED FINITE-ELEMENT PROBLEMS 11 the residual for each one of the components of the velocity and the pressure, which in this case appear distributed unevenly. The pressure component provides the largest error, while the velocity components appear to converge faster. This will not be the case for the nonsymmetric example. Moreover, the interpolation error seems to be higher than the energy in the case of the pressure and this is also reflected globally. Our guess is that this is to do with imposing condition (4.3) numerically. We remark here that Theorem 3.1 is only applicable to the global solution w = (u, v, p) and there is no reason to expect that the criterion should work componentwise. We also examined the convergence of the symmetrically norm-preconditioned GMRES applied to the Navier-Stokes problem. More precisely, we computed the solution to the nonlinear problem using a Picard iteration, but displayed only the results corresponding to the last linear system solve. As before, the convergence can be examined globally or separately for the different components of the solution. These curves are shown in Figs 4.4, 4.5 for two values of the diffusion parameter: ε = 0.1, ε = Again, the criterion works fine; moreover, it seems that in this case, the ε=0.1; N=2113; solver=sym PGMRES ε=0.1; N=2113; solver=sym PGMRES uvp u (a) Global (uvp) convergence (b) Component u convergence ε=0.1; N=2113; solver=sym PGMRES ε=0.1; N=2113; solver=sym PGMRES v p (c) Component v (d) Component p Fig Convergence criteria for symmetrically preconditioned GMRES for ε = 0.1.
12 12 M. ARIOLI and D. LOGHIN ε=0.01; N=2113; solver=sym PGMRES ε=0.01; N=2113; solver=sym PGMRES uvp u (a) Global (uvp) convergence (b) Component u convergence ε=0.01; N=2113; solver=sym PGMRES ε=0.01; N=2113; solver=sym PGMRES v p (c) Component v (d) Component p Fig Convergence criteria for symmetrically preconditioned GMRES for ε = component convergence can be described by components of our criterion, a feature which did not work for MINRES. More precisely, plots (b), (c), (d) seem to indicate that velocity and pressure residuals can provide respective bounds for the velocity and pressure forward errors The nonlinear case. In this section, we present the results for the fully nonlinear Navier-Stokes problem (4.2). First, we show the simple run of unpreconditioned GMRES as a solver for the Picard iteration applied to the nonlinear problem (4.2). The Picard iteration takes the form (3.3), but we used the modified iteration (3.7), in order to both exhibit its convergence properties and highlight the relevance of our stopping criteria. The choice c = 10 1, q = 1/4 in (3.4) does not affect the number of nonlinear iterations in our tests. Note that the norm in (3.4) is the Euclidean one. The purpose of this criterion, as described in [11] is to make GMRES work hard only when it matters (i.e., when the residual is sufficiently small). The nonlinear convergence history is displayed in Fig More precisely, the
13 STOPPING CRITERIA FOR MIXED FINITE-ELEMENT PROBLEMS 13 ε=0.1; N=2113; solver=gmres uvp (a) Nonlinear convergence and our stopping criterion (b) u difference (c) v difference (d) p difference Fig Convergence criteria for unpreconditioned GMRES for ε = 0.1 and the difference between the solution converged using (3.4) and using our proposed stopping criterion (3.5). 2-norm of the residuals r k, concatenated from all nonlinear iterations (7 in this case) is plotted together with the energy-error and the H 1 -norm of the residual, which is computed exactly for illustration purposes. Of course, in the case of unpreconditioned GMRES it is not clear how one could derive an approximation for r H 1, except through direct computation. We plotted also in Fig. 4.6 the difference between the final solution obtained using criterion (3.4) and that obtained by employing the H 1 -norm of the residual. The errors are of the order 10 3 for the velocities and for the pressure (given solutions of order one), which indeed are of the order of the M error. In Fig 4.7 we display the convergence of the 3-term GMRES method in the last Picard step of the nonlinear iteration of type (3.7) using the hermitian part to symmetrically precondition the system. We point out that the choice τ = in (3.6) did not change the number of nonlinear iterations. Finally, we present the results obtained using the 3-term GMRES algorithm suggested in Section for solving the full nonlinear problem. First, we note that different choices of c, q will lead to different nonlinear convergence curves. We present a typical example in Fig 4.8 (with c = 1, q = 0.5) and
14 14 M. ARIOLI and D. LOGHIN ε=0.1; N=2113; solver=three TERM PGMRES ε=0.01; N=2113; solver=three TERM PGMRES uvp uvp (a) Global convergence for ε = 0.1 (b) Global convergence for ε = 0.01 Fig Convergence criteria for symmetrically preconditioned 3-term GMRES. ε=0.1; N=2113; solver=sym PGMRES ε=0.01; N=2113; solver=sym PGMRES uvp uvp (a) ε = 0.1 (b) ε = 0.01 Fig Convergence criteria for the full nonlinear problem using symmetricallypreconditioned GMRES-Picard for ε = 0.1, 0.01 highlight convergence properties for several choices of parameters in Tables 4.1, 4.2. In particular, we chose to work with c = c(h), for three values of q. We set tol out= and highlighted the number of iterations needed for this classic criterion compared with that suggested in the algorithm above where tol out= ηc(h)c 2. We worked with η = C 2 = 1 and C(h) = h 2 for ε = 0.1 and C(h) = h 3/2 for ε = 0.01; we note that this leads to a robust stopping criterion which we highlight in the vertical lines across the convergence curves in Fig 4.8. As expected, a small value of c = c(h) leads to best performance in the H 1 -norm. In particular, the GMRES-Picard algorithm is most wasteful when we attempt to solve each iteration to the full M error level, i.e., when c(h) = h 2. On the other hand, when c(h) = h 1/2, the convergence is relatively robust with respect to the q parameter. 5. Conclusion. We showed how the results described in [1] can be used and extended in the framework of mixed and mixed-hybrid finite-element approximation of partial differential equation systems in saddle-point form.
15 STOPPING CRITERIA FOR MIXED FINITE-ELEMENT PROBLEMS 15 q = 0.25 q = 0.5 q = 0.75 c(h) = dual classic dual classic dual classic h 1/ h h Table 4.1 Total number of preconditioned GMRES iterations for the full nonlinear solution of our test problem with ε = 0.1 using both the classic (l2) and dual stopping criteria q = 0.25 q = 0.5 q = 0.75 c(h) = dual classic dual classic dual classic h 1/ h h Table 4.2 Total number of preconditioned GMRES iterations for the full nonlinear solution of our test problem with ε = 0.01 using both the classic (l2) and dual stopping criteria Moreover, we described how the dual norm of the residual can be easily used within classical Krylov methods to obtain reliable and efficient stopping criteria. Finally, we described how to generalize these techniques to the nonlinear case, thus obtaining a considerable gain in efficiency. In particular, we showed how the use of the dual norm of the residual (essentially, the energy norm of the error) can be successfully combined with a short term recurrence GMRES in order to solve nonlinear saddle-point problems such as Navier-Stokes equations. RERENCES [1] M. Arioli, D. Loghin, and A. Wathen, Stopping criteria for iterations in finite-element methods, Numer. Math., 99 (2005), pp [2] D. N. Arnold, R. S. Falk, and R. Winther, Preconditioning discrete approximations of the Reissner-Mindlin plate model, RAIRO Model. Math. Anal. Numer., 31 (1997), pp [3] S. Berrone, Adaptive discretization of stationary and incompressible navier stokes equations by stabilized finite element methods, Comput. Meth. Appl. Engrg., 190 (2001), pp [4] J. H. Bramble and J. E. Pasciak, A preconditioning technique for indefinite systems resulting from mixed approximations of elliptic problems, Math. Comp., 50 (1988), pp [5], Iterative techniques for time dependent Stokes problems, Computers Math. Applic., 33 (1997), pp [6] F. Brezzi and M. Fortin, Mixed and Hybrid Finite Element Methods, vol. 15, Springer-Verlag, Berlin, [7] B. M. Brown, P. K. Jimack, and M. D. Mihajlovic, An efficient direct solver for a class of mixed finite element problems, Applied Numerical Mathematics, 38 (2000), pp [8] J. Cahouet and J. P. Chabard, Some fast 3D finite element solvers for the generalized Stokes problem, Internat. J. Numer. Methods Fluids, 8 (1988), pp [9] H. H. Chen and J. C. Strikwerda, Preconditioning for regular elliptic systems, SIAM J. Num. Anal., 37 (1999), pp [10] Z. Chen, R. E. Ewing, and R. Lazarov, Domain decomposition algorithms for mixed methods for second order elliptic problems, Math. Comp., 65 (1996), pp [11] R. S. Dembo, S. C. Eisenstat, and T. Steihaug, Inexact Newton methods, SIAM J. Numer.
16 16 M. ARIOLI and D. LOGHIN Anal., 19 (1982), pp [12] L. Demkowicz, Babuška Brezzi??, Tech. Rep , ICES University of Texas at Austin, [13] B. Fischer, A. J. Wathen, and D. J. Silvester, The convergence of iterative solution methods for symmetric and indefinite linear systems, in Numerical analysis 1997 (Dundee), Longman, Harlow, 1998, pp [14] R. Glowinski and O. Pironneau, Numerical methods for the first biharmonic equation and for the two-dimensional Stokes problem, SIAM Rev., 21 (1979), pp [15] R. Glowinski and M. Wheeler, Domain decompostion and mixed finite element methods for elliptic problems, in First International Symposium on Domain Decomposition Methods for Partial Differential Equations, R. G. et. al., ed., SIAM, Philadelphia, 1988, pp [16] T. J. R. Hughes, L. Franca, and M. Balestra, A new finite element formulation for computational fluid dynamics: V. circumventing the Babuška-Brezzi condition: a stable Petrov-Galerkin formulation of the Stokes problem accommodating equal-order interpolations, Comput. methods Appl. Mech. Engrg., 59 (1986), pp [17] A. Klawonn, Block-triangular preconditioners for saddle-point problems with a penalty term, SIAM J. Sci. Comput., 19 (1998), pp [18] D. Loghin and A. J. Wathen, Analysis of preconditioners for saddle-point problems, SIAM J. Sci. Comput., 25 (2004), pp (electronic). [19] M. D. Mihajlovic and D. J. Silvester, A black-box multigrid preconditioner for the biharmonic equation, Tech. Rep. 362, UMIST, [20] I. Perugia and V. Simoncini, Optimal and quasi-optimal preconditioners for certain mixed finite element approximations, Tech. Rep. 1098, Istituto di Analisi Numerica del C.N.R., Pavia, [21] I. Perugia, V. Simoncini, and M. Arioli, Linear algebra methods in a mixed approximation of magnetostatic problems, SIAM J. Sci. Comput., 21 (1999), pp [22] T. Rusten and R. Winther, Substructure preconditioners for elliptic saddle point problems, Math. Comp., 60 (1993), pp [23] D. J. Silvester and A. J. Wathen, Fast iterative solution of stabilised Stokes systems. II. Using general block preconditioners., SIAM J. Num. Anal., 31 (1994), pp [24] P. S. Vassilevski and J. Wang, Multilevel iterative methods for mixed finite element discretizations of elliptic problems, Numer. Math., 63 (1992), pp [25] A. J. Wathen and D. J. Silvester, Fast iterative solution of stabilised Stokes systems. I. Using simple diagonal preconditioners., SIAM J. Num. Anal., 30 (1993), pp [26] J. Xu and L. Zikatanov, Some observations on Babuška and Brezzi theories, Numerische Mathematik, 94 (2003), pp
A 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 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 informationStopping criteria for iterations in finite element methods
Stopping criteria for iterations in finite element methods M. Arioli, D. Loghin and A. J. Wathen CERFACS Technical Report TR/PA/03/21 Abstract This work extends the results of Arioli [1], [2] on stopping
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 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 informationStructured Preconditioners for Saddle Point Problems
Structured Preconditioners for Saddle Point Problems V. Simoncini Dipartimento di Matematica Università di Bologna valeria@dm.unibo.it p. 1 Collaborators on this project Mario Arioli, RAL, UK Michele Benzi,
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 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 informationSome New Elements for the Reissner Mindlin Plate Model
Boundary Value Problems for Partial Differential Equations and Applications, J.-L. Lions and C. Baiocchi, eds., Masson, 1993, pp. 287 292. Some New Elements for the Reissner Mindlin Plate Model Douglas
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 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 informationNumerical behavior of inexact linear solvers
Numerical behavior of inexact linear solvers Miro Rozložník joint results with Zhong-zhi Bai and Pavel Jiránek Institute of Computer Science, Czech Academy of Sciences, Prague, Czech Republic The fourth
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 informationOn the accuracy of saddle point solvers
On the accuracy of saddle point solvers Miro Rozložník joint results with Valeria Simoncini and Pavel Jiránek Institute of Computer Science, Czech Academy of Sciences, Prague, Czech Republic Seminar at
More informationPreconditioned GMRES Revisited
Preconditioned GMRES Revisited Roland Herzog Kirk Soodhalter UBC (visiting) RICAM Linz Preconditioning Conference 2017 Vancouver August 01, 2017 Preconditioned GMRES Revisited Vancouver 1 / 32 Table of
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 informationUniform inf-sup condition for the Brinkman problem in highly heterogeneous media
Uniform inf-sup condition for the Brinkman problem in highly heterogeneous media Raytcho Lazarov & Aziz Takhirov Texas A&M May 3-4, 2016 R. Lazarov & A.T. (Texas A&M) Brinkman May 3-4, 2016 1 / 30 Outline
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 informationStructured Preconditioners for Saddle Point Problems
Structured Preconditioners for Saddle Point Problems V. Simoncini Dipartimento di Matematica Università di Bologna valeria@dm.unibo.it. p.1/16 Collaborators on this project Mario Arioli, RAL, UK Michele
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 informationAdaptive preconditioners for nonlinear systems of equations
Adaptive preconditioners for nonlinear systems of equations L. Loghin D. Ruiz A. Touhami CERFACS Technical Report TR/PA/04/42 Also ENSEEIHT-IRIT Technical Report RT/TLSE/04/02 Abstract The use of preconditioned
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 informationChebyshev semi-iteration in Preconditioning
Report no. 08/14 Chebyshev semi-iteration in Preconditioning Andrew J. Wathen Oxford University Computing Laboratory Tyrone Rees Oxford University Computing Laboratory Dedicated to Victor Pereyra on his
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 informationarxiv: v1 [math.na] 31 Oct 2014
arxiv:1411.0634v1 [math.na] 31 Oct 014 A NOTE ON STABILITY AND OPTIMAL APPROXIMATION ESTIMATES FOR SYMMETRIC SADDLE POINT SYSTEMS CONSTANTIN BACUTA Abstract. We establish sharp well-posedness and approximation
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 informationNumerical Methods for Large-Scale Nonlinear Equations
Slide 1 Numerical Methods for Large-Scale Nonlinear Equations Homer Walker MA 512 April 28, 2005 Inexact Newton and Newton Krylov Methods a. Newton-iterative and inexact Newton methods. Slide 2 i. Formulation
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 informationNatural preconditioners for saddle point systems
Natural preconditioners for saddle point systems Jennifer Pestana and Andrew J. Wathen August 29, 2013 Abstract The solution of quadratic or locally quadratic extremum problems subject to linear(ized)
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 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 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 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 informationINF-SUP CONDITION FOR OPERATOR EQUATIONS
INF-SUP CONDITION FOR OPERATOR EQUATIONS LONG CHEN We study the well-posedness of the operator equation (1) T u = f. where T is a linear and bounded operator between two linear vector spaces. We give equivalent
More informationAcceleration of a Domain Decomposition Method for Advection-Diffusion Problems
Acceleration of a Domain Decomposition Method for Advection-Diffusion Problems Gert Lube 1, Tobias Knopp 2, and Gerd Rapin 2 1 University of Göttingen, Institute of Numerical and Applied Mathematics (http://www.num.math.uni-goettingen.de/lube/)
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 informationNONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS
NONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS J.-L. GUERMOND 1, Abstract. This paper analyzes a nonstandard form of the Stokes problem where the mass conservation
More informationSpectral Properties of Saddle Point Linear Systems and Relations to Iterative Solvers Part I: Spectral Properties. V. Simoncini
Spectral Properties of Saddle Point Linear Systems and Relations to Iterative Solvers Part I: Spectral Properties V. Simoncini Dipartimento di Matematica, Università di ologna valeria@dm.unibo.it 1 Outline
More informationSTABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS
STABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS RAYTCHO LAZAROV AND XIU YE Abstract. In this paper, we derive two stabilized discontinuous finite element formulations, symmetric
More informationOn the interplay between discretization and preconditioning of Krylov subspace methods
On the interplay between discretization and preconditioning of Krylov subspace methods Josef Málek and Zdeněk Strakoš Nečas Center for Mathematical Modeling Charles University in Prague and Czech Academy
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 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 informationDomain Decomposition Algorithms for Saddle Point Problems
Contemporary Mathematics Volume 218, 1998 B 0-8218-0988-1-03007-7 Domain Decomposition Algorithms for Saddle Point Problems Luca F. Pavarino 1. Introduction In this paper, we introduce some domain decomposition
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationIterative solvers for saddle point algebraic linear systems: tools of the trade. V. Simoncini
Iterative solvers for saddle point algebraic linear systems: tools of the trade V. Simoncini Dipartimento di Matematica, Università di Bologna and CIRSA, Ravenna, Italy valeria@dm.unibo.it 1 The problem
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 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 informationIncomplete LU Preconditioning and Error Compensation Strategies for Sparse Matrices
Incomplete LU Preconditioning and Error Compensation Strategies for Sparse Matrices Eun-Joo Lee Department of Computer Science, East Stroudsburg University of Pennsylvania, 327 Science and Technology Center,
More informationJournal of Computational and Applied Mathematics. Optimization of the parameterized Uzawa preconditioners for saddle point matrices
Journal of Computational Applied Mathematics 6 (009) 136 154 Contents lists available at ScienceDirect Journal of Computational Applied Mathematics journal homepage: wwwelseviercom/locate/cam Optimization
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 informationThe antitriangular factorisation of saddle point matrices
The antitriangular factorisation of saddle point matrices J. Pestana and A. J. Wathen August 29, 2013 Abstract Mastronardi and Van Dooren [this journal, 34 (2013) pp. 173 196] recently introduced the block
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 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 informationINNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD
INNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD Abstract. Finite element methods for the Reissner Mindlin plate theory are discussed. Methods in which both the tranverse displacement and
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 informationAMS 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 informationTwo-scale Dirichlet-Neumann preconditioners for boundary refinements
Two-scale Dirichlet-Neumann preconditioners for boundary refinements Patrice Hauret 1 and Patrick Le Tallec 2 1 Graduate Aeronautical Laboratories, MS 25-45, California Institute of Technology Pasadena,
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 18 Outline
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 informationAn Accelerated Block-Parallel Newton Method via Overlapped Partitioning
An Accelerated Block-Parallel Newton Method via Overlapped Partitioning Yurong Chen Lab. of Parallel Computing, Institute of Software, CAS (http://www.rdcps.ac.cn/~ychen/english.htm) Summary. This paper
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 informationDepartment of Computer Science, University of Illinois at Urbana-Champaign
Department of Computer Science, University of Illinois at Urbana-Champaign Probing for Schur Complements and Preconditioning Generalized Saddle-Point Problems Eric de Sturler, sturler@cs.uiuc.edu, http://www-faculty.cs.uiuc.edu/~sturler
More informationFEM for Stokes Equations
FEM for Stokes Equations Kanglin Chen 15. Dezember 2009 Outline 1 Saddle point problem 2 Mixed FEM for Stokes equations 3 Numerical Results 2 / 21 Stokes equations Given (f, g). Find (u, p) s.t. u + p
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 12-10 Fast linear solver for pressure computation in layered domains P. van Slingerland and C. Vuik ISSN 1389-6520 Reports of the Department of Applied Mathematical
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 informationAN AUGMENTED LAGRANGIAN APPROACH TO LINEARIZED PROBLEMS IN HYDRODYNAMIC STABILITY
AN AUGMENTED LAGRANGIAN APPROACH TO LINEARIZED PROBLEMS IN HYDRODYNAMIC STABILITY MAXIM A. OLSHANSKII AND MICHELE BENZI Abstract. The solution of linear systems arising from the linear stability analysis
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 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 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 informationON THE ROLE OF COMMUTATOR ARGUMENTS IN THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS FOR STOKES CONTROL PROBLEMS
ON THE ROLE OF COUTATOR ARGUENTS IN THE DEVELOPENT OF PARAETER-ROBUST PRECONDITIONERS FOR STOKES CONTROL PROBLES JOHN W. PEARSON Abstract. The development of preconditioners for PDE-constrained optimization
More informationOUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative methods ffl Krylov subspace methods ffl Preconditioning techniques: Iterative methods ILU
Preconditioning Techniques for Solving Large Sparse Linear Systems Arnold Reusken Institut für Geometrie und Praktische Mathematik RWTH-Aachen OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative
More informationPreconditioning of Saddle Point Systems by Substructuring and a Penalty Approach
Preconditioning of Saddle Point Systems by Substructuring and a Penalty Approach Clark R. Dohrmann 1 Sandia National Laboratories, crdohrm@sandia.gov. Sandia is a multiprogram laboratory operated by Sandia
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 informationA MULTIGRID METHOD FOR THE PSEUDOSTRESS FORMULATION OF STOKES PROBLEMS
SIAM J. SCI. COMPUT. Vol. 29, No. 5, pp. 2078 2095 c 2007 Society for Industrial and Applied Mathematics A MULTIGRID METHOD FOR THE PSEUDOSTRESS FORMULATION OF STOKES PROBLEMS ZHIQIANG CAI AND YANQIU WANG
More informationTrust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization
Trust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization Denis Ridzal Department of Computational and Applied Mathematics Rice University, Houston, Texas dridzal@caam.rice.edu
More informationLeast-squares Finite Element Approximations for the Reissner Mindlin Plate
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl., 6, 479 496 999 Least-squares Finite Element Approximations for the Reissner Mindlin Plate Zhiqiang Cai, Xiu Ye 2 and Huilong Zhang
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 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 informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationPreconditioning for Nonsymmetry and Time-dependence
Preconditioning for Nonsymmetry and Time-dependence Andy Wathen Oxford University, UK joint work with Jen Pestana and Elle McDonald Jeju, Korea, 2015 p.1/24 Iterative methods For self-adjoint problems/symmetric
More informationA Review of Solution Techniques for Unsteady Incompressible Flow
Zeist 2009 p. 1/57 A Review of Solution Techniques for Unsteady Incompressible Flow David Silvester School of Mathematics University of Manchester Zeist 2009 p. 2/57 PDEs Review : 1966 1999 Update : 2000
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 37, pp. 166-172, 2010. Copyright 2010,. ISSN 1068-9613. ETNA A GRADIENT RECOVERY OPERATOR BASED ON AN OBLIQUE PROJECTION BISHNU P. LAMICHHANE Abstract.
More informationAdditive Average Schwarz Method for a Crouzeix-Raviart Finite Volume Element Discretization of Elliptic Problems
Additive Average Schwarz Method for a Crouzeix-Raviart Finite Volume Element Discretization of Elliptic Problems Atle Loneland 1, Leszek Marcinkowski 2, and Talal Rahman 3 1 Introduction In this paper
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 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 informationUnsteady Incompressible Flow Simulation Using Galerkin Finite Elements with Spatial/Temporal Adaptation
Unsteady Incompressible Flow Simulation Using Galerkin Finite Elements with Spatial/Temporal Adaptation Mohamed S. Ebeida Carnegie Mellon University, Pittsburgh, PA 15213 Roger L. Davis and Roland W. Freund
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 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 informationThe Nullspace free eigenvalue problem and the inexact Shift and invert Lanczos method. V. Simoncini. Dipartimento di Matematica, Università di Bologna
The Nullspace free eigenvalue problem and the inexact Shift and invert Lanczos method V. Simoncini Dipartimento di Matematica, Università di Bologna and CIRSA, Ravenna, Italy valeria@dm.unibo.it 1 The
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 informationNUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL
Proceedings of ALGORITMY 212 pp. 29 218 NUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL MARTIN HADRAVA, MILOSLAV FEISTAUER, AND PETR SVÁČEK Abstract. The present paper
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 informationTHE solution of the absolute value equation (AVE) of
The nonlinear HSS-like iterative method for absolute value equations Mu-Zheng Zhu Member, IAENG, and Ya-E Qi arxiv:1403.7013v4 [math.na] 2 Jan 2018 Abstract Salkuyeh proposed the Picard-HSS iteration method
More informationCONVERGENCE BOUNDS FOR PRECONDITIONED GMRES USING ELEMENT-BY-ELEMENT ESTIMATES OF THE FIELD OF VALUES
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 CONVERGENCE BOUNDS FOR PRECONDITIONED GMRES USING ELEMENT-BY-ELEMENT
More informationON AUGMENTED LAGRANGIAN METHODS FOR SADDLE-POINT LINEAR SYSTEMS WITH SINGULAR OR SEMIDEFINITE (1,1) BLOCKS * 1. Introduction
Journal of Computational Mathematics Vol.xx, No.x, 200x, 1 9. http://www.global-sci.org/jcm doi:10.4208/jcm.1401-cr7 ON AUGMENED LAGRANGIAN MEHODS FOR SADDLE-POIN LINEAR SYSEMS WIH SINGULAR OR SEMIDEFINIE
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 informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 13-10 Comparison of some preconditioners for the incompressible Navier-Stokes equations X. He and C. Vuik ISSN 1389-6520 Reports of the Delft Institute of Applied
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationLocal Mesh Refinement with the PCD Method
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 125 136 (2013) http://campus.mst.edu/adsa Local Mesh Refinement with the PCD Method Ahmed Tahiri Université Med Premier
More information