1 Introduction In this paper, we propose and study a balancing Neumann-Neumann domain decomposition method for incompressible Stokes equations. Neuman
|
|
- Victor Dickerson
- 6 years ago
- Views:
Transcription
1 Balancing Neumann-Neumann methods for incompressible Stokes equations Luca F. Pavarino Λ Olof B. Widlund y February 7, 1 Abstract Balancing Neumann-Neumann methods are introduced and studied for incompressible Stokes equations discretized with mixed finite or spectral elements with discontinuous pressures. After decomposing the original domain of the problem into nonoverlapping subdomains, the interior unknowns, which are the interior velocity component and all except the constant pressure component, of each subdomain problem are implicitly eliminated. The resulting saddle point Schur complement is solved with a Krylov space method with a balancing Neumann-Neumann preconditioner based on the solution of a coarse Stokes problem with a few degrees of freedom per subdomain and on the solution of local Stokes problems with natural and essential boundary conditions on the subdomains. This preconditioner is of hybrid form in which the coarse problem is treated multiplicatively while the local problems are treated additively. The condition number of the preconditioned operator is independent of the number of subdomains and is bounded from above by the product of the square of the logarithm of the local number of unknowns in each subdomain and the inverse of the inf-sup constants of the discrete problem and of the coarse subproblem. Numerical results show that the method is quite fast; they are also fully consistent with the theory. Keywords. Stokes equations, mixed finite elements, spectral elements, preconditioned iterations, domain decomposition, balancing Neumann-Neumann methods AMS(MOS) subject classifications. 65N3, 65N35, 65N55 Λ Department of Mathematics, Universit a di Milano, Via Saldini 5, 133 Milano, ITALY. address: pavarino@mat.unimi.it, URL: This work was supported in part by the National Science Foundation under Grant NSF-CCR y Courant Institute of Mathematical Sciences, 51 Mercer Street, New York, NY 11. E- mail address: widlund@cs.nyu.edu. URL: This work was supported in part by the National Science Foundation under Grants NSF-CCR- 9738, and in part by the U.S. Department of Energy under contract DE-FG-9ER517. 1
2 1 Introduction In this paper, we propose and study a balancing Neumann-Neumann domain decomposition method for incompressible Stokes equations. Neumann-Neumann type algorithms, together with FETI type algorithms and overlapping Schwarz methods, have emerged in the last decade among the most efficient and thoroughly tested preconditioned iterative solvers for the parallel solution of very large linear systems arising from finite or spectral element discretizations of elliptic partial differential equations; see, e.g., Gropp, Kaushik, Keyes, and Smith [16] and Bhardwaj, Day, Farhat, Lesoinne, Pierson, and Rixen [5]. We refer to Smith, Bjfirstad, and Gropp [35] and Quarteroni and Valli [3] for a general introduction to domain decomposition methods. Most of the theoretical and numerical work for Neumann-Neumann methods has been carried out for second order elliptic problems; see Mandel [5], Mandel and Brezina [6], Cowsar, Mandel and Wheeler [9], Dryja and Widlund [1], and Pavarino [9]. More recently, this family of methods has been extended to plate and shell problems, see Le Tallec, Mandel, and Vidrascu [], to convection-diffusion problems, see Alart, Barboteu, Le Tallec, and Vidrascu [3] and Achdou, Le Tallec, Nataf, and Vidrascu [1], and to vector field problems, see Toselli [37]. We also note that the connection between Neumann-Neumann and FETI methods has been considered recently by Klawonn and Widlund [1]. In this paper, we extend the balancing Neumann-Neumann method to symmetric saddle point problems arising from incompressible Stokes equations. Previous work on domain decomposition methods for incompressible Stokes equations has been based on other iterative substructuring methods, see Bramble and Pasciak [6], Pasciak [8], Quarteroni [31], Marini and Quarteroni [7], Fischer and Rfinquist [1], Casarin [8], Rfinquist [33], Le Tallec and Patra [3], Pavarino and Widlund [3], and Ainsworth and Sherwin [], on overlapping Schwarz methods, see Gervasio [14], Fischer [11], Fischer, Miller, and Tufo [13], Klawonn and Pavarino [19], and Rfinquist [34], and on block preconditioners, see Klawonn [18, 17], Klawonn and Pavarino [] and the references therein. Our balancing Neumann-Neumann preconditioner will be built from a number of components. They include solvers of Dirichlet problems on the subdomains Ω i ; into which the given region Ω has been subdivided. Here we use the fact that these are well-posed problems which are solvable for any boundary values provided that the constraints corresponding to the pressure components which are constant in each subdomain are dropped. We also use Neumann solvers for the subregions to construct local components of the preconditioner. Finally, we use a coarse space correction which maps any velocity field given on the interface ; i.e., the union of the interior part of the subdomain boundaries, onto a balanced velocity field, i.e., a vector field with values on which allow a divergence free extension of the given boundary values into each subdomain. We must be careful and choose the coarse space so that it is stable, i.e., satisfies the appropriate infsup condition. It is well known that a coarse space correction is necessary, in addition to many local corrections, in the design of successful, scalable domain decompo-
3 sition preconditioners; see, e.g., Widlund [38]. The design and study of such coarse spaces are in fact at the core of many contributions to domain decomposition theory. For balancing NeumannNeumann methods, these issues have again come to the forefront in recent work by Achdou, Le Tallec, Nataf, and Vidrascu [1] and Alart, Barboteu, Le Tallec, and Vidrascu [3]. In the present study, a coarse space will also play a special and important role. Our coarse space correction will return the error vector of the iteration to the benign subspace of balanced velocities; the Stokes operator is symmetric, positive definite on this subspace and this allows us to use a preconditioned conjugate gradient method to achieve avery fast convergence rate. We also note that we are quite fortunate in that we can use a discontinuous pressure approximation in our finite element approximations of Stokes equations. This makes it possible to regard all except the piece-wise constant pressure components as local variables, together with the interior velocity modes, resulting in a Schur complement system which involves only the interface velocities and one pressure parameter per subdomain. This sets the stage for the construction of a coarse space which includes the piecewise constant pressure components, one for each subdomain, and just a few velocity degrees of freedom for each of them. In our bounds of the condition number of the preconditioned system, the infsup constant of this coarse problem, as well as the infsup constant of the entire finite element model, will enter; naturally, we will build our preconditioner from stable components representing solvers on subdomains and a coarse, global problem. The rest of this paper is organized as follows. In section, we briefly describe the Stokes system, certain finite and spectral element discretizations, and the resulting discrete system. The substructuring process, also known as static condensation, is described in section 3, both in matrix and variational form. The balancing Neumann-Neumann preconditioner is introduced in section 4, first in matrix form and then in the variational framework needed for the theoretical analysis. The main result of the paper and some auxiliary results are proven in section 5. Section 6 concludes the paper with some numerical results on a model Stokes problem in the plane. The Stokes system: the continuous and discrete problems Let Ω ρ R d ;d ; 3; be a polyhedral domain and L (Ω) be the subspace of L (Ω) of functions with zero mean value. The Stokes problem consists of finding the velocity u (H 1 (Ω)) d and the pressure p L (Ω) of an incompressible fluid 3
4 with viscosity ν by solving: 8 >< >: ν Z Z Ω Ω u j@ω g; ru : rvdx Z divv pdx Z Ω Ω f vdx 8v (H 1 (Ω)) d ; divu qdx 8q L (Ω); (1) R where g (H 1 (@Ω)) d ; f (H 1 (Ω)) d. We note that if g nds ; g is balanced, then the solution u of this Dirichlet problem will be divergence free. For simplicity, we will set g from now on. We also note that the techniques and results of this paper can be applied to the mixed formulation of the linear elasticity equations as well. We will also need to consider Stokes problems with natural boundary conditions: νru i n pn i r i i 1; ;d; () derived by using Green's formula. In contrast with the case of Dirichlet boundary conditions, the pressure is now uniquely determined and the pressure space should now be taken to be L (Ω). In this case, as for the Laplace operator, each velocity component is determined only up to a constant and there is a compatibility condition between r i and f i ; namely, Z dx Ω i1 f i dx + Z i1 r i ds : We note that if the boundary conditions are mixed, part essential, part natural, then there is a unique solution without any compatibility conditions. We will consider conforming discretizations of Stokes equations using finite as well as spectral finite elements, both with discontinuous pressures..1 Finite element methods with discontinuous pressures We assume that the domain Ω can be decomposed into N nonoverlapping subdomains Ω i of characteristic size H forming a quadrilateral (hexahedral) finite element mesh fi H. This coarse triangulation is further refined into a fine quadrilateral finite element triangulation fi h of characteristic size h. Among the many choices of mixed finite elements available for Stokes equations, we consider the following: a) Q (h) Q (h) mixed finite elements: the velocity space V is discretized by continuous, piecewise quadratic velocities: V h fv V : v kjt Q (T ); 8T fi h ; k 1; ;dg; while the pressure space is discretized by discontinuous piecewise constant functions on fi h U h fq U : q jt Q (T ); 8T fi h g: 4
5 These elements satisfy the uniform inf-sup condition (divv;q) sup fi h kqk vv kvk L 8q U h ; (3) h H 1 with fi h c> independent ofh, but they lead to nonoptimal error estimates; see Brezzi and Fortin [7, Ch. VI.4, p. 1]. b) Q (h)p 1 (h) mixed finite elements: the velocity space is as before, while the pressure space consists of piecewise linear discontinuous pressures: U h fq U : q jt P 1 (T ); 8T fi h g: These elements satisfy a uniform inf-sup condition (3) as well; there is also an optimal O(h ) error estimates for both velocities and pressures; see Brezzi and Fortin [7, Ch. VI, p. 16].. Spectral element methods: Q n Q n Let Ω ref be the reference square or cube (1; 1) d ;d ; 3, and let Q n (Ω ref )be the set of polynomials on Ω ref of degree n in each variable. We assume that the domain Ω can be decomposed into N nonoverlapping finite elements Ω i, each of which is an image Ω i ffi i (Ω ref ), with ffi i an affine mapping. V is discretized, component by component, by continuous, piecewise polynomials of degree n: V n fv V : v kjωi ffi ffi i Q n (Ω ref ); i 1; ;N; k 1; ;dg: The pressure space is discretized by piecewise polynomials of degree n, discontinuous across the boundaries of the elements Ω i : U n fq U : q jωi ffi ffi i Q n (Ω ref );i1; ;Ng: These mixed spectral elements are implemented using Gauss-Lobatto-Legendre (GLL(n)) quadrature, which also allows the construction of a very convenient nodal tensor-product basis for V n. Denote by fο i g n i the set of GLL(n) points of [1; 1], and by ff i the quadrature weight associated with ο i : Let l i (x) bethe Lagrange interpolating polynomial of degree n which vanishes at all the GLL(n) nodes except at ο i, where it equals one. Each element ofq n (Ω ref ) is expanded in the GLL(n) basis and each L inner product of two scalar components u and v is replaced, in the three dimensional case, by (u; v) n;ω nx s1 i;j;k (u ffi ffi s )(ο i ;ο j ;ο k )(v ffi ffi s )(ο i ;ο j ;ο k )jj s jff i ff j ff k ; where jj s j is the determinant of the Jacobian of ffi s. The mass matrix based on these basis elements and GLL(n) quadrature are diagonal. Similarly, a very convenient basis for U n consists of the tensor-product Lagrangian nodal basis functions associated with the internal GLL(n) nodes, i.e., the endpoints -1 and +1 are excluded. We will call these the pressure GLL(n) nodes. 5
6 The Q n Q n method satisfies a nonuniform inf-sup condition (divv;q) sup fi n kqk vv kvk L 8q U n ; (4) n H 1 d1 ( where fi n Cn ), d ; 3; and the constant C is independent ofn and q; see Maday, Meiron, Patera, and Rfinquist [4] and Stenberg and Suri [36]. However, numerical experiments, reported in [4], have also shown that for practical values of n, e.g., n 16, the inf-sup constant fi n of the Q n Q n method decays much slower than what would be first expected from the theoretical bound. An alternative, with a uniform bound on the inf-sup constant, is provided by the Q n P n1 method; see Bernardi and Maday [4]. However, this pressure space is less convenient than Q n as far as implementation is concerned..3 The discrete system Let e V and e U be the discrete velocity and pressure spaces. In the finite element case, we write e V e U V h U h, while in the spectral element case we write ev e U V n U n. Define the standard bilinear forms a(u; v) ν Z and the linear functional Ω ru : rvdx; b(v;p) F(v) Z Ω f vdx: Z Ω divv pdx The discrete system obtained from (1) using finite or spectral elements is: Find u e V and p e U such that ( a(u; v) + b(v;p) F(v) 8v e V b(u;q) 8q e U; (5) where we denote with the same letters the bilinear forms obtained using the appropriate quadrature rule described above. In matrix form, we have u A B T u b K : p B p On the subspace V B KerB fv e V : b(v;q) 8q e Ug; the Stokes problem (5) is equivalent to the positive definite problem: Find u V B such that a(u; v) F(v) 8v V B : 6
7 3 Substructuring for saddle point problems The domain Ω is decomposed into open, nonoverlapping quadrilateral (hexahedral) subdomains Ω i and the interface ; i.e., Ω[ N i1ω i [ : SN i1 i Each Ω i typically consists of one, or a few, spectral elements of degree n or of many finite elements. We denote by h h the set of nodes belonging to the interface respectively. The starting point of our algorithm is the implicit elimination of the interior degrees of freedom, i.e., the velocity component that is supported in the open subdomains and what we will call the interior pressure component which has zero average over the individual subdomains. This process, also known as static condensation, is carried out by solving decoupled local Stokes problems on each subdomain Ω i with Dirichlet boundary conditions for the velocities given i. We then obtain a saddle point Schur complement problem for the interface velocities and a constant pressure in each subdomain. This reduced problem will be solved by a preconditioned Krylov space iteration normally the preconditioned conjugate gradient method. For simplicity, we will use the same letters to denote both functions and their associated vector representations; the same convention will also be used for linear operators and their associated matrix forms. 3.1 Substructuring in matrix form In order to eliminate the interior degrees of freedom, we reorder the vector of unknowns as 6 4 u I p I u p interior velocities interior pressures with zero average interface velocities constant pressures in each Ω i : Then, after using the same permutation, the discrete Stokes system matrix can be written as 4 K II K T I K I K A II B T II A T I B II B I A I B T I A B T B : Eliminating the interior unknowns u I and p I by static condensation, we obtain the saddle point Schur complement system S u p ~ b ; (6) 7
8 where and A B T B ~ b b S K K I K 1 AI B T I S B T B AI B T I II KT I AII B T II B II ; AII B T II B II 1 A T I B I 1 bi : By using a second permutation that reorders the interior velocities and pressures subdomain by subdomain, we note that K 1 II represents the solution of N decoupled Stokes problems, one for each subdomain and all uniquely solvable, with Dirichlet data given i K 1 II 6 4 K (1) 1 II... K (N) 1 II : This is the matrix associated with the discrete Stokes extension operator SH described in the next section. The Schur complement S does not need to be explicitly assembled since only its action Sv on a vector v is needed in a Krylov iteration. This operation essentially only requires the action of K 1 II onavector, i.e., the solution of N decoupled Stokes problems. In other words, Sv is computed by subassembling the actions of the subdomain Schur complements S (i) defined for Ω i ; by " Once A (i) B (i)t B (i) u # S (i) K (i) K(i) I (K(i) II )1 K (i)t " is known, p p I ui We note that A (i) I B (i)t I " #" S (i) B (i)t B (i) p I ui AII B II T B II p I ui A (i) II B (i) # I B (i)t II II : # 1 " A (i)t I B (i) I can be found by back-substitution, 1 bi ( is independent ofp : 8 A T I B I u ): p #
9 3. Iterative substructuring in variational form The substructuring procedure described in the previous section is associated with the space decomposition where the interior spaces are defined as ev e U Φ N i1 V i U i Φ V U ; V i e V H 1 (Ω i ) U i e U L (Ω i ); and the spaces of interface velocities and coarse pressures, constant in each subdomain, are defined as V SH(e V)fv e V : vjωi SH(v j@ωi ); i 1; ;Ng; U fq e U : qjωi constant; i 1; ;Ng: Here SH : Vj e! V; e is the velocity component of the discrete Stokes harmonic extension operator that maps an interface velocity u Vj e onto the solution (~u; ~p) of the following homogeneous Stokes problem, defined on each subdomain separately: Find ~u e V and ~p P N i1 U i such that on each Ω i 8 < : a(~u; v) + b(v; ~p) 8v V i b(~u;q) 8q U i ~u u i : The following comparison of the energy of the discrete Stokes extensions SH and the discrete harmonic extensions H of each velocity component can be found in [15] and [6] for finite element discretizations and in [3] and [8] for spectral element discretizations. We note that the corresponding local bounds, for individual subdomains are equally valid and that the upper bound has an elementary proof. Lemma 1 cfia(shu ; SHu ) a(hu ; Hu ) a(shu ; SHu ) 8u V ; where fi is the inf-sup constant of the chosen mixed finite element space ~ V ~ U. If we define the interface inner product by s(u ; v )a(sh(u ); SH(v )) u T S v ; and by b (u ;p ) the restriction of the other bilinear form to the Stokes harmonic extensions, the variational formulation of the saddle point Schur complement problem (6) can be given by (see Pavarino and Widlund [3, Lemma 7.1]): 9
10 Find u V and p U such that, ρ s(u ; v ) + b (v ;p ) e F(v ) 8v V b (u ;q ) 8q U : (7) Problem (7) is equivalent to the positive definite problem: Find u V ;B such that s(u ; v )e F(v ) 8v V ;B ; (8) where V ;B is the subspace of balanced velocities defined by V ;B KerB fv V : B v g fv e V : vjωi SH(v j@ωi ) and i v n ; i 1; ;Ng: The equivalence of problems (7) and (8) follows from the fact that V ;B consists of the functions of SH(e V) that satisfy the zero flux constraints f Z@Ω i v n ; i 1; ;Ng. We obtain (7) by using Lagrange multipliers p (p ;i ) N i1 to enforce the constraints given in the definition of the balanced subspace. 4 A Neumann-Neumann preconditioner We will solve the saddle point Schur complement problem u S B S T u b ~ p B p (9) by a preconditioned Krylov space method such as GMRES or PCG. The latter can be applied to this indefinite problem because we will start and keep the iterates in the positive definite subspace V ;B. 4.1 Matrix form of the preconditioner We propose a balancing Neumann-Neumann preconditioner based on the solution of a coarse Stokes problem with a few degrees of freedom per subdomain and of local Stokes problems with natural and essential boundary conditions on each subdomain. This preconditioner is of hybrid form in which the coarse problem is treated multiplicatively while the local problems are treated additively; cf. [35, p. 15]. It is closely analogous to the balancing Neumann-Neumann preconditioner for the positive definite case, except that the coarse and local problems are saddle point problems. The matrix form of the preconditioner is Q Q H +(I Q H S) 1 i1 Q i (I SQ H );
11 where the coarse operator Q H and local operators Q i are defined below. The preconditioned operator is then T QS T +(I T ) i1 T i (I T ); where T Q H S and T i Q i S. Q can be written as a three-step preconditioner similar to Smith, Bjfirstad, and Gropp [35, p. 16]: Given a residual vector r the preconditioned vector ^r Qr is given by t ψ Q H r; t ψ i1 Q i (r St); ^r ψ t + Q H (r St): We note that the first step can be left out if r is a residual in the range of S(I T ); we will see that T is a projection and that therefore T (I T ): Coarse solver: Given a residual vector r, the coarse term Q H r is the solution of a coarse, global Stokes problem with a few velocity degrees of freedom and one constant pressure per subdomain Ω i : where and Q H RHS T 1 R H; L T R H ; I S R H SRH T L T S L L T B T B L : (1) We will consider four choices for the matrix L. Some of the columns of L are defined in terms of the Neumann-Neumann counting functions μ i associated with each subdomain Ω i : μ i is zero at the interface nodes i while its value at anynodeon@ω i equals the number of subdomains shared by that node. Its pseudo inverse μ y i is the function 1μ i(x) for all nodes where μ i (x) 6 ; and it vanishes at all other points of h h : We note that we use the function μ y i in all or almost all of the subdomains and for each velocity component. Then the columns of L are defined by one of the following four choices: ) the inverse counting functions μ y i ; 1) the μ y i and the continuous coarse piecewise bi- or tri-linear functions; ) the μ y i and the continuous coarse piecewise bi- or tri-quadratic functions; 3) the μ y i and the quadratic coarse edge/face bubble functions for the normal direction: 11
12 Choice ) provides a quite minimal coarse velocity space and it turns out not to be inf-sup stable; it is included since it is the standard choice for the laplacian. Choices 1), ), and 3) are enrichments of ) that all turn out to be inf-sup stable uniformly in N. In particular, ) is a very natural choice because it is based directly on a well known discretization of the Stokes problem on the coarse finite or spectral elements. We note that 1) and ) might not be easy to adapt to unstructured subdomain meshes produced by automatic mesh partitioners. This is one of the reasons why we also have considered choice 3), which should be relatively easy to adapt for general unstructured subdomains. In order to avoid linearly dependent μ y i functions, and hence a singular coarse space problem, we might have to drop all of the components of these functions for one subdomain. This depends on the coarse triangulation: on a regular hexahedral mesh in three dimensions, we should, e.g., include at most 3(N 1) such functions. We will use the same symbols μ i and μ y i for both the interface vectors and the functions obtained by extending them inside each subdomain by using the discrete Stokes harmonic extension operator SH. We can explicitly represent the velocity and pressure components of the coarse operator Q H S in matrix form. Lemma Let e S L T S L, e S B L e S 1 LT BT and consider v (I Q H S) u : Then, v (I P e P)u ; Gu ; where P and e P are projections defined by P L e S 1 LT S L e S 1 LT B T e S 1 B L e S 1 LT S ; ep L e S 1 LT B T e S 1 B ; and G e S 1 B L e S 1 LT S + e S 1 B : We note that e Pu ifu is balanced and that (I P e P)u is balanced because B (I P e P). Moreover, we note that v and are independent of, and that P P ; e P e P; and e PP P e P : Proof. These formulas follow by explicitly computing the inverse of the coarse matrix (1) by using the following explicit formula, proven by a direct computation. Let A be an invertible matrix and let B have full row rank. Then, A B T B 1 A 1 A 1 B T S 1 BA 1 A 1 B T S 1 S 1 BA 1 S 1 ; 1
13 where S BA 1 B T : Local solver: The local operators Q i will only be applied to residuals of balanced velocity fields and thus the second residual component will vanish. We have also shown in Lemma that the pressure components obtained in this step of the preconditioner plays no further role when we next apply the operator (I T ): Each local operator Q i is based on the solution of a local Stokes problem on Ω i with natural boundary condition. This local problem is nonsingular for any subdomain Ω i the boundary of which but it is singular otherwise, i.e., for the floating subdomains. In the latter case any constant velocity belongs to the null space, while the pressure is now uniquely determined because a pressure term is present in the Neumann boundary condition; see (). Toavoid possible complications with singular problems, we modify the local Stokes problems on the floating subdomains, by adding ffl times the velocity mass matrix to the local stiffness matrix K (i). We could also make these solutions unique by requiring that each velocity component has a zero average over Ω i ; the right hand sides will always be compatible. Given a residual vector with a first component r and a zero second component, Q i r is the weighted solutions of a local Stokes problem on subdomain Ω i with a natural boundary condition i Q i r R T i D 1 i " S (i) ;ffl B (i) B (i)t # 1 D 1 i R i r : (11) Here R i are ; 1 restriction matrices mapping r into r i and D i are diagonal matrices representing multiplication by the counting functions μ i. Moreover, S (i) ffl " S (i) ;ffl B (i) B (i)t is the local saddle point Schur complement, associated with subdomain Ω i,of the regularized local stiffness matrix where K (i) ffl 6 4 A (i) II;ffl B (i)t II B (i) II B (i) A (i) I;ffl # A (i)t I;ffl B (i)t I A (i) ;ffl B (i) A (i) ffl A (i) + fflm (i) : I B (i)t ffl is a positive parameter, and M (i) is the local velocity mass matrix ; 13
14 4. Variational formulation of the preconditioner The balancing Neumann-Neumann preconditioner introduced in the previous section is associated with further decomposing the interface space V U as V U V U + i1 V i U ;i ; where the coarse and local spaces are defined below. Coarse space: We will consider the following four choices for the coarse velocity space V V fv V : v spanfμ y i gg; V 1 fv V : v spanfμ y i g[q 1()g; V fv V : v spanfμ y i g[q ()g; V 3 fv V : v spanfμ y ig[spanfnormal direction edgeface bubble functionsgg: Here Q 1 () and Q () are the subspaces of V of piecewise linear and quadratic polynomials on, respectively. All these coarse spaces contain the kernel of the diagonally scaled local Stokes problems on the floating subdomains Ω i given in formula (13) that are used in the definition of the local components of the preconditioner. Coarse problem: Given u V U, define u V U as the solution of the coarse Stokes problem ρ s(u u ; v ) + b (v ;) 8v V b (u u ;q) 8q U : T u (1) We recall that the matrix form of this operator is T Q H S R T H S1 R HS where S is given by (1). Local spaces: The local spaces are defined by V ;i fv V : v(x) 8x h i;h g; U ;i spanfq ;i g: Local problems: We need to apply local Stokes solvers only to elements u w of Range(I T ). Since (I T ) is a projection on its range, we u have w (I T ) u (I T ) u ; according to Lemma, u (I P )u is balanced and Gu. T u We define T i w T i u i u V ;i U ;i as the solution of a local T p i u Stokes problem with natural boundary conditions: 8v i V ;i ; 8q U ;i, ρ si;ffl (μ i Ti u u ;μ i v i )+b ;i (μ i v i ;T p i u ) s(u ; v i )+b (v i ;Gu ) b ;i (μ i Ti u u ;q) b (u ;q) (13) 14
15 and Here ffl is a positive parameter, s i;ffl (u ; v )a i;ffl (SH i;ffl (u ); SH i;ffl (v )); a i;ffl (u; v) Z Ω i ru : rvdx + ffl Z Ω i u vdx; (14) and SH i;ffl is the velocity component of the discrete Stokes extension operator defined in terms of the regularized a i;ffl ( ; ) velocity bilinear form instead of the standard a( ; ) form. By μ i Ti u u and μ i v i we mean the finite or spectral element function with nodal values equal to the product of those of μ i and the other function. We recall that the matrix form of the operators T i is given by T i Q i S, with Q i given in (11). From the second equation in (13), we find that μ i T u i u is balanced i since u is balanced. If we also consider test functions in the same subspace, we see that the velocity component T u i u satisfies: 8v i V ;i such that μ i v i is balanced i : s i;ffl (μ i T u i u ;μ i v i )s(u ; v i )+b (v i ;Gu ): (15) Given that the right-hand side of (13) is compatible, even for ffl,we could alternatively choose to solve this singular system and enforce the constraints for each floating subdomain. Z Ω i SH(μ i T u i u )dx 5 Auxiliary and main results We need a few more auxiliary results before we can formulate and prove our main result. We will work with the S -inner product Du ; v E S s(u ; v )u T S v : On the balanced subspace V ;B U this expression coincides with the bilinear form defined by S fi u since B u B v. ; v fl S fi S u ; v fl ; 15
16 u Lemma 3 Let a) b) c) fi S u fi u S fi S u Range(I T ) and let ; v ; v ; v fl fl fl where (I P e P)v is the velocity component of u (I T ) v v be arbitrary. Then, fi fl S(I T ) u ; v fi fl u S ; (I T ) v ; fi fl u S ; v ; E DS u ; (I P P)v e DS (I P e P)u ; v E DS (I P )u ; v E ; (I T ) v Proof. a) Since Range(I T ) and I T is a projection onto that u subspace, we have (I T ) u : Since S and Q H are symmetric and T Q H S; we have fi S(I Q H S) u b) Since u ; v fi S u fl : fi (I Q H S) u ; (I Q H S) v fl : ;S v Range(I T ), u is balanced, and we find that fi u S ; v fi S u + B T fl ; v fi S u + B T fl fi S u ; v ; v fl fl : c) The first formula is obtained by applying a) and b) and Lemma : fi S u fi S u + B T ; v fl ; (IT ) v fi u S fl 16 ; (I T ) v fl fl u T S (IP e P)v + T B (IP e P)v
17 E u T S (I P P)v e DS u ; (I P P)v e : The second formula follows similarly and the third follows from the second since B u and therefore e Pu : Since V essentially is an enrichment ofq (); the following lemma follows from the stability of the well-known Q Q finite elements and the stability of the discrete Stokes extension; see Pavarino and Widlund [3, Lemma 7.] for a proof for the spectral element case and Casarin [8, Ch. 5] for related results. Lemma 4 The coarse space V U satisfies the inf-sup condition sup v V where fi is independent of h; n; N. (divsh(v );q ) a(sh(v ); SH(v )) fi kq k L 8q U ; Our numerical results, reported in section 6, indicate that a uniform inf-sup condition does not hold for the first coarse space V U. The results for all the others are quite satisfactory although we do not have a full theory. We note that the Q 1 Q elements by themselves are not inf-sup table but that we are using a richer velocity space which includes the μ y i functions. We also work in the somewhat different context of Stokes extensions of traces on : Lemma 5 ki P e PkS (1 + 1 fi ) where fi is the inf-sup constant of the coarse space. Proof. Any inf-sup stable problem satisfies a stability estimate given in Brezzi and Fortin [7, Th. 1.1, p. 4, eq. 1.9]): s(u ku ; v ) k S sup + b (u ;q) sup ; v V kv k S fi qu kqk L where s(u ; v ) and b (u ;q) are the two components of the right-hand side of the coarse problem (1). By the continuity of these two forms, it then follows By Lemma, u (P + e P)u ; hence and therefore ku k S (1 + fi )ku k S : ku k S k(p + e P)u k S (1 + fi )ku k S ; ki P e PkS kik S + kp + e PkS (1 + 1 fi ): 17
18 Theorem 1 On the balanced subspace V ;B U the balancing Neumann-Neumann operator T is symmetric positive definite with respect to the S bilinear form and where ff 8 < : cond(t ) C(1 + 1 fi ) 1 fi ff; (1 + log(hh)) for finite elements (1 + log n) for spectral elements; fi and fi are the inf-sup constants of the coarse problem and the original discrete Stokes problem, respectively. u Proof. Let u be balanced and let w. Then, D D D T w; w ES T w; w ES X + (I T ) i DX ES + D T w;t w D E T w; (I T )w i E T i (I T )w; w S E T i (I T )w; (I T )w ; (16) S because. Since T is an orthogonal projection on the S P balanced subspace, we are left with providing a lower and upper bound for i T i on D Range(I ET ); once such an inequality has been obtained we will simply add T w;t w to each side of the inequalities. We can limit ourselves to proving S lower and upper bounds for only P the velocity component i T i u: Since u is (I P )u balanced, we have (I T )w and the velocity component of Gu is X i T i (I T )w X Therefore, by Lemma 3c, DX i T i (IT )w; (IT )w DX i i P i Q i Q i [S (I P )u + B T Gu ] X i ES DX T u i u ; (I P )u ES i (I P )u S Gu T u i u : Q i [S (IP )u +B T Gu ];S (IP )u E D(I P e P) X i T u i u ; u ES : (17) 18
19 u Lower bound: Let be an element in Range(I T ); then u V ;B. Since the pseudoinverses P μ y i of the counting functions define a partition N of unity, wehave u u i1 i with u i μ y i u V ;i. Since u is balanced, we can use formula (15) and obtain because, i1 i1 Du ; u E s i;ffl (μ i T u i u ;μ i u i ) i1 S s(u ; u ) i1 b (u i ;)b ( b (u i ;) i1 i1 i1 s(u ; u i ) u i ;)b (u ;): By the Cauchy-Schwarz inequality and formula (15), s i;ffl (μ i T u i u ;μ i u i ); (18) X N 1 X N 1 s i;ffl (μ i Ti u u ;μ i u i ) s i;ffl (μ i Ti u u ;μ i Ti u u ) s i;ffl (μ i u i ;μ i u i ) i1 i1 N X i1 s(u ;T u i u )+ s(u ; i1 i1 T u i u )+b ( 1 X N 1 b (Ti u u ;Gu ) s i;ffl (u ; u ) i1 i1 The first factor in (19) is, by Lemma 3c, s(u ; i1 T u i u )+b ( i1 D u P N ; T i1 i u u Gu? T u i u ;Gu ) 1sffl (u ; u ) 1 : (19) D Ti u S u u ;Gu ) + B T Gu Here, the value of? makes no difference. estimated as follows: ES D u ; (I P P) ; i1 P N i1 T u i u? T u i u ES : () The second factor in (19) can be s ffl (u ; u )a ffl (SH ffl (u ); SH ffl (u )) a ffl (SH(u ); SH(u )) a(sh(u ); SH(u )) + fflksh(u )k L (Ω) (1 + ffl ffi )a(sh(u ); SH(u )) (1+ ffl ffi )s(u ; u ); (1) E 19
20 where a(sh(v ); SH(v )) ffi inf > : v ksh(v )k L (Ω) From (19), (), and (1), we then have i1 D s i;ffl (μ i Ti u u ;μ i u ) u ; (I P P) i1 Therefore, from this last estimate and (18), we find Du ; u ES (1 + ffl ffi ) D u ; (I P P) E 1 Ti u u (1 + ffl 1: S ffi )s(u ; u ) i1 T u i u ES ; and by (17) 1 D (1 + ffl ffi ) w; w ES DX i E T i (I T )w; (I T )w : () S Upper bound: We start by estimating the S -norm of an individual T u i u. We note that D T u i u ;T u i u E S s(t u i u ;T u i u )a(sh(t u i u ); SH(T u i u )): By comparing the energy of the discrete Stokes and harmonic extensions as in Lemma 1, we have a(sh(t u i u ); SH(T u i u )) C fi a(h(t u i u ); H(T u i u )): Since T u i u V ;i,we can now apply to each scalar component ofh(t u i u )a result proven in the scalar case in Dryja and Widlund [1, Lemma 4] for finite elements and in Pavarino [9, Lemma 6.] for spectral elements, and obtain a(h(t u i u ); H(T u i u )) C ffa i;ffl (H i;ffl (μ i T u i u ); H i;ffl (μ i T u i u )); where a i;ffl ( ; ) is given in formula (14) and ff 8 < : (1 + log(hh)) for finite elements (1 + log n) for spectral elements: By using a local variant of Lemma 1 for each subdomain Ω i,we can return to the discrete Stokes extensions: a i;ffl (H i;ffl (μ i T u i u ); H i;ffl (μ i T u i u )) a i;ffl (SH i;ffl (μ i T u i u ); SH i;ffl (μ i T u i u )) s i;ffl (μ i T u i u ;μ i T u i u ):
21 Therefore, D T u i u ;T u i u E S C ff fi s i;ffl(μ i T u i u ;μ i T u i u ): By formula (15), Lemma 3c, and the Cauchy-Schwarz inequality, we have Therefore s i;ffl (μ i T u i u ;μ i T u i u )s(u ;T u i u )+b(t u D u Gu ; T u i u? E 1 Du ; u D T u i u ;T u i u E i u ;Gu ) ES D u ; (I P e P)T u i u ES S ki P e PkS D S T u E 1 i u ;Ti u u S : C ff fi ki P e Pk S D u ; u ES : By a standard coloring argument (see, e.g., Dryja and Widlund [1, Th. 1, p. 18]), we obtain D N X i1 T u i u ; i1 Hence by the Cauchy-Schwarz inequality, D N X i1 T u i u ; u E E D Ti u u C ff S fi ki P Pk e S u ; u : ES S C ff fi k(i P e PkS D u ; u ES : Since u (I P )u,we find by using Lemma 3c, that and by (17) D N (I P P) e X E Ti u u ; u C ff D S fi k(i P PkS e u ; u ; ES i1 DX i T i (I T )w; (I T )w ES C ff D E fi k(i P PkS e w; w S : (3) In conclusion, we find by putting together the lower bound () and the upper bound (3) that 1 D (1 + ffl ffi ) w; w ES DX i E T i (IT )w; (IT )w C ki P PkS e D E ff w; w S fi ; S and we can conclude our proof by using Lemma 5 and equation (16). 1
22 6 Numerical results We report, in this last section, results of numerical experiments, carried out in Matlab 5.3 on Unix workstations, for a model Stokes problem on the unit square and with homogeneous Dirichlet boundary conditions. The problem was discretized p p with Q n Q n spectral elements and the domain Ω divided into N N square subdomains. After the implicit elimination of the interior unknowns, the saddle point Schur complement (9) is solved iteratively by either GMRES with restart every 5 iterations or PCG. The initial guess is always zero, the right hand side is random and uniformly distributed, and the stopping criterion is kr k k kr k 1 6, where r k is the residual at the kth iterate. As explained in section 4, the singularity of the local Neumann solves for the floating subdomains is avoided by shifting the diagonal of the local velocity stiffness matrices by ffl 1 5 : The Laplace operator. In order to have a reference to compare with, we first considered the Laplace instead of the Stokes operator. The results are reported in Table 1. In the upper half of the table the number of subdomains, N 3 3; is fixed, while the spectral degree n is increased from to 1; in the lower half the spectral degree n 4 is fixed and the number of subdomains N is increased from to1 1. The first three columns report the iteration count, and the maximum and minimum eigenvalues when CG is applied without any preconditioner. The next three columns report the same data when PCG is applied with a balancing Neumann-Neumann preconditioner with a coarse space spanned by only the μ y i functions associated with floating subdomains. The final three columns report the same data when PCG is applied with a balancing Neumann-Neumann preconditioner with a coarse space spanned by all the μ y i functions, except one in order to avoid a linearly dependent μ y i functions in this particular geometry. These results clearly show that CG without preconditioning is neither optimal nor scalable, since the iteration counts grow with n and N, while PCG with a balancing Neumann-Neumann preconditioner is quasi-optimal and scalable; the iteration counts appear to grow at most logarithmically with n and are independent ofn. The richer coarse space based on all μ y i functions but one yields a better preconditioner (the condition number is almost half of that of the other balancing preconditioner), but the iteration counts are essentially the same. The growth of the maximum eigenvalue of the preconditioned operator with n is also plotted in the left panel of Figure 3. The Stokes operator. In the next four tables, we report analogous sets of results for the Stokes operator. The saddle point Schur complement (9) for Stokes equations is solved by GMRES or PCG with and without our balancing Neumann-Neumann preconditioner. Each table corresponds to one of the four choices of coarse spaces described in section 4, V in Table, V 1 in Table 3, V in Table 4, and V 3 in Table 5. In the upper half of each table the number of subdomains, N 3 3; is fixed, while the spectral degree n is increased from to 1; in the lower half the spectral degree n 4 is fixed and the number of subdomains N is increased from to 1 1. In each table, the first column reports the iteration counts of GMRES(5) without precondition-
23 ing, the second column the iteration counts of GMRES(5) with our balancing Neumann-Neumann preconditioner. The next three columns report the iteration counts and extreme eigenvalues of PCG with our balancing Neumann-Neumann preconditioner. The last column reports the square of the inf-sup constant fi of the coarse Stokes problem, computed as the minimum nonzero eigenvalue of M 1 (B p L )(L T S L ) 1 (B L ) T, where M p is the coarse pressure mass matrix and (B L )(L T S L ) 1 (B L ) T is the pressure Schur complement for the coarse problem (1). The results of these tables are also plotted in Figure 1 (GM- RES iteration counts), Figure (PCG iteration counts), Figure 3 (maximum eigenvalue of the Laplacian and Stokes operators), Figure 4 (inverse coarse infsup constant1fi ). These results show that GMRES without preconditioning is neither optimal nor scalable, while both GMRES and PCG with our balancing Neumann-Neumann preconditioner are quasi-optimal and scalable, except with the first choice of coarse space V. In fact, the V coarse space does not seem to be inf-sup stable (last column of Table and right panel of Figure 4) and the iteration counts of both GMRES and PCG seem to grow linearly with N in that case (right panels of Figures 1 and ). The other three coarse spaces appear to be inf-sup stable uniformly in N and in fact the iteration counts of both GMRES and PCG are bounded from above independently of N. It is more difficult to detect numerically the dependence of the coarse inf-sup constant with respect to n; see right panel of Figure 4. Overall, these numerical results seem to be in good agreement with the estimates of our main result: since we are considering here the two dimensional case and spectral elements, Theorem 1 states that the maximum eigenvalue of the preconditioner operator T grows as O( p n(1 + log (n))), because fi n Cn 1 (see eq. (4)), while the minimum eigenvalue of T is always bounded from below by1(1 + ffl ffi ). We note that the condition numbers of T for the Stokes operator in Tables 3 and 4 are remarkably close to (almost always less than twice) the corresponding condition numbers of T for the Laplacian operator in Table 1. The situation would be even better if we had considered stable finite elements, since both the global and coarse inf-sup constants would be independent ofh and h: Therefore the maximum eigenvalue of T would grow only as O(1 + log (Hh)). 3
24 Table 1: Laplace solver with balancing Neumann-Neumann preconditioner. Fixed number of subdomains N 3 3 NO PREC. NN PREC. FLOAT ONLY NN PREC. ALL μ y i 1 n it. max min it. max min it. max min Fixed spectral degree: n 4 NO PREC. NN PREC. FLOAT ONLY NN PREC. ALL μ y i 1 N it. max min it. max min it. max min
25 Table : Stokes solver with balancing Neumann-Neumann preconditioner and first choice of coarse space V. Fixed number of subdomains N 3 3 GMRES GMRES NN PCG NN n it. it. it. max min fi e e e e e e e e e e e- Fixed spectral degree: n 4 GMRES GMRES NN PCG NN N it. it. it. max min fi e e e * e * e * e * e * e * e * e * e-3 5
26 Table 3: Stokes solver with balancing Neumann-Neumann preconditioner and second choice of coarse space V. 1 Fixed number of subdomains N 3 3 GMRES GMRES NN PCG NN n it. it. it. max min fi e e e e e e e e e e e-1 Fixed spectral degree: n 4 GMRES GMRES NN PCG NN N it. it. it. max min fi e e e * e * e * e * e * e * e * e * e-1 6
27 Table 4: Stokes solver with balancing Neumann-Neumann preconditioner and third choice of coarse space V. Fixed number of subdomains N 3 3 GMRES GMRES NN PCG NN n it. it. it. max min fi () () (1.3) (1.) e e e e e e e e e e-1 Fixed spectral degree: n 4 GMRES GMRES NN PCG NN N it. it. it. max min fi e e e * e * e * e * e * e * e * e * e-1 7
28 Table 5: Stokes solver with balancing Neumann-Neumann preconditioner and fourth choice of coarse space V. 3 Fixed number of subdomains N 3 3 GMRES GMRES NN PCG NN n it. it. it. max min fi e e e e e e e e e e e-1 Fixed spectral degree: n 4 GMRES GMRES NN PCG NN N it. it. it. max min fi e e * e * e * e * e * e * e * e * e-1 8
29 References [1] Y. Achdou, P. L. Tallec, F. Nataf, and M. Vidrascu, A domain decomposition preconditioner for an advection-diffusion problem, Comput. Methods Appl. Mech. Engrg, 184 (), pp [] M. Ainsworth and S. Sherwin, Domain decomposition preconditioners for p and hp finite element approximation of Stokes equations, Comput. Methods Appl. Mech. Engrg., 175 (1999), pp [3] P. Alart, M. Barboteu, P. L. Tallec, and M. Vidrascu, Méthode de Schwarz additive avec solveur grossier pour probl emes non symétriques, Comp. Rend. Acad. Sci. Paris, serie I, (), pp [4] C. Bernardi and Y. Maday, Uniform inf-sup conditions for the spectral discretization of the Stokes problem, Math. Models Methods Appl. Sci, 9 (1999), pp [5] M. Bhardwaj, D. Day, C. Farhat, M. Lesoinne, K. Pierson, and D. Rixen, Application of the FETI method to ASCI problems: Scalability results on one thousand processors and discussion of highly heterogeneous problems, Int. J. Numer. Meth. Engng., 47 (), pp [6] J. Bramble and J. Pasciak, A domain decomposition technique for Stokes problems, Appl. Numer. Math., 6 (1989/9), pp [7] F. Brezzi and M. Fortin, Mixed and Hybrid Finite Element Methods, Springer-Verlag, Berlin, [8] M. A. Casarin, Schwarz Preconditioners for Spectral and Mortar Finite Element Methods with Applications to Incompressible Fluids, PhD thesis, Dept. of Mathematics, Courant Institute of Mathematical Sciences, New York University, March [9] L. C. Cowsar, J. Mandel, and M. F. Wheeler, Balancing domain decomposition for mixed finite elements, Math. Comp., 64 (1995), pp [1] M. Dryja and O. B. Widlund, Schwarz methods of Neumann-Neumann type for three-dimensional elliptic finite element problems, Comm. Pure Appl. Math., 48 (1995), pp [11] P. Fischer, An overlapping Schwarz method for spectral element solution of the incompressible Navier-Stokes equations, J. Comp. Phys., 133 (1) (1997), pp [1] P. Fischer and E. Rfinquist, Spectral element methods for large scale parallel Navier-Stokes calculations, Comp. Meths. Appl. Mech. Eng., 116 (1994), pp
30 [13] P. F. Fischer, N. I. Miller, and H. M. Tufo, An overlapping Schwarz method for spectral element simulation of three-dimensional incompressible flows, in Parallel Solution of PDE, P. Bjfirstad and M. Luskin, eds., no. 1 in IMA Volumes in Mathematics and Its Applications, Springer-Verlag,, pp. 13. [14] P. Gervasio, Risoluzione delle equazioni alle derivate parziali con metodi spettrali in regioni partizionate in sottodomini, PhD thesis, Universit a di Milano, [15] V. Girault and P.-A. Raviart, Finite Element Methods for Navier- Stokes equations, Springer-Verlag, Berlin, [16] W. Gropp, D. Kaushik, D. E. Keyes, and B. F. Smith, High performance parallel implicit CFD, Tech. Rep. ANL/MCS-P863-1, Argonne National Laboratory,. [17] A. Klawonn, Block-triangular preconditioners for saddle point problems with a penalty term, SIAM J. Sci. Comp., 19 (1998), pp [18], An optimal preconditioner for a class of saddle point problems with a penalty term, SIAM J. Sci. Comp., 19 (1998), pp [19] A. Klawonn and L. F. Pavarino, Overlapping Schwarz methods for mixed linear elasticity and Stokes problems, Comp. Meths. Appl. Mech. Eng., 165 (1998), pp [], A comparison of overlapping Schwarz methods and block preconditioners for saddle point problems, Numer. Linear Alg. Appl., 7 (), pp. 15. [1] A. Klawonn and O. B. Widlund, FETI and Neumann-Neumann iterative substructuring methods: connections and new results, Comm. Pure Appl. Math., 54 (1), pp [] P. Le Tallec, J. Mandel, and M. Vidrascu, A Neumann-Neumann domain decomposition algorithm for solving plate and shell problems, SIAM J. Numer. Anal., 35 (1998), pp [3] P. Le Tallec and A. Patra, Non-overlapping domain decomposition methods for adaptive hp approximations of the Stokes problem with discontinuous pressure fields, Comp. Meths. Appl. Mech. Eng., 145 (1997), pp [4] Y. Maday, D. Meiron, A. Patera, and E. Rfinquist, Analysis of iterative methods for the steady and unsteady Stokes problem: application to spectral element discretizations, SIAM J. Sci. Comp., 14 (1993), pp
Overlapping 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 informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
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 information18. Balancing Neumann-Neumann for (In)Compressible Linear Elasticity and (Generalized) Stokes Parallel Implementation
Fourteenth nternational Conference on Domain Decomposition Methods Editors: smael Herrera, David E Keyes, Olof B Widlund, Robert Yates c 23 DDMorg 18 Balancing Neumann-Neumann for (n)compressible Linear
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY, USA. Dan_Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More informationDomain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions
Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Bernhard Hientzsch Courant Institute of Mathematical Sciences, New York University, 51 Mercer Street, New
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 informationOVERLAPPING SCHWARZ ALGORITHMS FOR ALMOST INCOMPRESSIBLE LINEAR ELASTICITY TR
OVERLAPPING SCHWARZ ALGORITHMS FOR ALMOST INCOMPRESSIBLE LINEAR ELASTICITY MINGCHAO CAI, LUCA F. PAVARINO, AND OLOF B. WIDLUND TR2014-969 Abstract. Low order finite element discretizations of the linear
More informationConstruction of a New Domain Decomposition Method for the Stokes Equations
Construction of a New Domain Decomposition Method for the Stokes Equations Frédéric Nataf 1 and Gerd Rapin 2 1 CMAP, CNRS; UMR7641, Ecole Polytechnique, 91128 Palaiseau Cedex, France 2 Math. Dep., NAM,
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 11, pp. 1-24, 2000. Copyright 2000,. ISSN 1068-9613. ETNA NEUMANN NEUMANN METHODS FOR VECTOR FIELD PROBLEMS ANDREA TOSELLI Abstract. In this paper,
More informationSelecting Constraints in Dual-Primal FETI Methods for Elasticity in Three Dimensions
Selecting Constraints in Dual-Primal FETI Methods for Elasticity in Three Dimensions Axel Klawonn 1 and Olof B. Widlund 2 1 Universität Duisburg-Essen, Campus Essen, Fachbereich Mathematik, (http://www.uni-essen.de/ingmath/axel.klawonn/)
More informationOverlapping Schwarz Preconditioners for Spectral. Problem in H(curl)
Overlapping Schwarz Preconditioners for Spectral Nédélec Elements for a Model Problem in H(curl) Technical Report TR2002-83 November 22, 2002 Department of Computer Science Courant Institute of Mathematical
More informationShort title: Total FETI. Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ Ostrava, Czech Republic
Short title: Total FETI Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ-70833 Ostrava, Czech Republic mail: zdenek.dostal@vsb.cz fax +420 596 919 597 phone
More informationCONVERGENCE ANALYSIS OF A BALANCING DOMAIN DECOMPOSITION METHOD FOR SOLVING A CLASS OF INDEFINITE LINEAR SYSTEMS
CONVERGENCE ANALYSIS OF A BALANCING DOMAIN DECOMPOSITION METHOD FOR SOLVING A CLASS OF INDEFINITE LINEAR SYSTEMS JING LI AND XUEMIN TU Abstract A variant of balancing domain decomposition method by constraints
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 information1. Fast Solvers and Schwarz Preconditioners for Spectral Nédélec Elements for a Model Problem in H(curl)
DDM Preprint Editors: editor1, editor2, editor3, editor4 c DDM.org 1. Fast Solvers and Schwarz Preconditioners for Spectral Nédélec Elements for a Model Problem in H(curl) Bernhard Hientzsch 1 1. Introduction.
More informationMultispace and Multilevel BDDC. Jan Mandel University of Colorado at Denver and Health Sciences Center
Multispace and Multilevel BDDC Jan Mandel University of Colorado at Denver and Health Sciences Center Based on joint work with Bedřich Sousedík, UCDHSC and Czech Technical University, and Clark R. Dohrmann,
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 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 informationMultilevel and Adaptive Iterative Substructuring Methods. Jan Mandel University of Colorado Denver
Multilevel and Adaptive Iterative Substructuring Methods Jan Mandel University of Colorado Denver The multilevel BDDC method is joint work with Bedřich Sousedík, Czech Technical University, and Clark Dohrmann,
More information33 RASHO: A Restricted Additive Schwarz Preconditioner with Harmonic Overlap
Thirteenth International Conference on Domain Decomposition ethods Editors: N. Debit,.Garbey, R. Hoppe, J. Périaux, D. Keyes, Y. Kuznetsov c 001 DD.org 33 RASHO: A Restricted Additive Schwarz Preconditioner
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 informationExtending the theory for domain decomposition algorithms to less regular subdomains
Extending the theory for domain decomposition algorithms to less regular subdomains Olof Widlund Courant Institute of Mathematical Sciences New York University http://www.cs.nyu.edu/cs/faculty/widlund/
More informationNew constructions of domain decomposition methods for systems of PDEs
New constructions of domain decomposition methods for systems of PDEs Nouvelles constructions de méthodes de décomposition de domaine pour des systèmes d équations aux dérivées partielles V. Dolean?? F.
More informationcfl Jing Li All rights reserved, 22
NYU-CS-TR83 Dual-Primal FETI Methods for Stationary Stokes and Navier-Stokes Equations by Jing Li A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy
More informationConvergence analysis of a balancing domain decomposition method for solving a class of indefinite linear systems
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2000; 00:1 6 [Version: 2002/09/18 v1.02] Convergence analysis of a balancing domain decomposition method for solving a class of indefinite
More informationASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM
ASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM C. Canuto 1, L. F. Pavarino 2, and A. B. Pieri 3 1 Introduction Discontinuous Galerkin (DG) methods for partial differential
More informationIN p-version AND SPECTRAL ELEMENT METHODS MARIO A. CASARIN
DIAGONAL EDGE PRECONDITIONERS IN p-version AND SPECTRAL ELEMENT METHODS MARIO A. CASARIN Abstract. Domain decomposition preconditioners for high-order Galerkin methods in two dimensions are often built
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 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 informationParallel scalability of a FETI DP mortar method for problems with discontinuous coefficients
Parallel scalability of a FETI DP mortar method for problems with discontinuous coefficients Nina Dokeva and Wlodek Proskurowski University of Southern California, Department of Mathematics Los Angeles,
More informationParallel Scalability of a FETI DP Mortar Method for Problems with Discontinuous Coefficients
Parallel Scalability of a FETI DP Mortar Method for Problems with Discontinuous Coefficients Nina Dokeva and Wlodek Proskurowski Department of Mathematics, University of Southern California, Los Angeles,
More informationAdaptive Coarse Space Selection in BDDC and FETI-DP Iterative Substructuring Methods: Towards Fast and Robust Solvers
Adaptive Coarse Space Selection in BDDC and FETI-DP Iterative Substructuring Methods: Towards Fast and Robust Solvers Jan Mandel University of Colorado at Denver Bedřich Sousedík Czech Technical University
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationNonoverlapping Domain Decomposition Methods with Simplified Coarse Spaces for Solving Three-dimensional Elliptic Problems
Nonoverlapping Domain Decomposition Methods with Simplified Coarse Spaces for Solving Three-dimensional Elliptic Problems Qiya Hu 1, Shi Shu 2 and Junxian Wang 3 Abstract In this paper we propose a substructuring
More informationDomain Decomposition Methods for Mortar Finite Elements
Domain Decomposition Methods for Mortar Finite Elements Dan Stefanica Courant Institute of Mathematical Sciences New York University September 1999 A dissertation in the Department of Mathematics Submitted
More informationConvergence Behavior of a Two-Level Optimized Schwarz Preconditioner
Convergence Behavior of a Two-Level Optimized Schwarz Preconditioner Olivier Dubois 1 and Martin J. Gander 2 1 IMA, University of Minnesota, 207 Church St. SE, Minneapolis, MN 55455 dubois@ima.umn.edu
More informationTwo new enriched multiscale coarse spaces for the Additive Average Schwarz method
346 Two new enriched multiscale coarse spaces for the Additive Average Schwarz method Leszek Marcinkowski 1 and Talal Rahman 2 1 Introduction We propose additive Schwarz methods with spectrally enriched
More informationA FETI-DP Method for Mortar Finite Element Discretization of a Fourth Order Problem
A FETI-DP Method for Mortar Finite Element Discretization of a Fourth Order Problem Leszek Marcinkowski 1 and Nina Dokeva 2 1 Department of Mathematics, Warsaw University, Banacha 2, 02 097 Warszawa, Poland,
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. A virtual overlapping Schwarz method for scalar elliptic problems in two dimensions
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES A virtual overlapping Schwarz method for scalar elliptic problems in two dimensions Juan Gabriel Calvo Preprint No. 25-2017 PRAHA 2017 A VIRTUAL
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 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 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 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 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 informationAn additive average Schwarz method for the plate bending problem
J. Numer. Math., Vol. 10, No. 2, pp. 109 125 (2002) c VSP 2002 Prepared using jnm.sty [Version: 02.02.2002 v1.2] An additive average Schwarz method for the plate bending problem X. Feng and T. Rahman Abstract
More informationThe All-floating BETI Method: Numerical Results
The All-floating BETI Method: Numerical Results Günther Of Institute of Computational Mathematics, Graz University of Technology, Steyrergasse 30, A-8010 Graz, Austria, of@tugraz.at Summary. The all-floating
More informationParallel Sums and Adaptive BDDC Deluxe
249 Parallel Sums and Adaptive BDDC Deluxe Olof B. Widlund 1 and Juan G. Calvo 2 1 Introduction There has recently been a considerable activity in developing adaptive methods for the selection of primal
More informationFETI-DP for Elasticity with Almost Incompressible 2 Material Components 3 UNCORRECTED PROOF. Sabrina Gippert, Axel Klawonn, and Oliver Rheinbach 4
1 FETI-DP for Elasticity with Almost Incompressible 2 Material Components 3 Sabrina Gippert, Axel Klawonn, and Oliver Rheinbach 4 Lehrstuhl für Numerische Mathematik, Fakultät für Mathematik, Universität
More informationfor three dimensional problems are often more complicated than the quite simple constructions that work well for problems in the plane; see [23] for a
DUAL-PRIMAL FETI METHODS FOR THREE-DIMENSIONAL ELLIPTIC PROBLEMS WITH HETEROGENEOUS COEFFICIENTS AEL KLAWONN Λ, OLOF B. WIDLUND y, AND MAKSYMILIAN DRYJA z Abstract. In this paper, certain iterative substructuring
More informationON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD. January 2000
ON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD JAN MANDEL AND RADEK TEZAUR January 2000 Abstract In the Dual-Primal FETI method, introduced by Farhat et al [5], the domain is decomposed into
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 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 informationLeast-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations
Least-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations by Wilhelm Heinrichs Universität Duisburg Essen, Ingenieurmathematik Universitätsstr.
More information2 CAI, KEYES AND MARCINKOWSKI proportional to the relative nonlinearity of the function; i.e., as the relative nonlinearity increases the domain of co
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2002; 00:1 6 [Version: 2000/07/27 v1.0] Nonlinear Additive Schwarz Preconditioners and Application in Computational Fluid
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 informationThe mortar element method for quasilinear elliptic boundary value problems
The mortar element method for quasilinear elliptic boundary value problems Leszek Marcinkowski 1 Abstract We consider a discretization of quasilinear elliptic boundary value problems by the mortar version
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 informationAn Efficient FETI Implementation on Distributed Shared Memory Machines with Independent Numbers of Subdomains and Processors
Contemporary Mathematics Volume 218, 1998 B 0-8218-0988-1-03024-7 An Efficient FETI Implementation on Distributed Shared Memory Machines with Independent Numbers of Subdomains and Processors Michel Lesoinne
More informationOn the Use of Inexact Subdomain Solvers for BDDC Algorithms
On the Use of Inexact Subdomain Solvers for BDDC Algorithms Jing Li a, and Olof B. Widlund b,1 a Department of Mathematical Sciences, Kent State University, Kent, OH, 44242-0001 b Courant Institute of
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
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 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 informationPreconditioned Mixed Spectral Element Methods for. Elasticity and Stokes Problems. Luca F. Pavarino. Department of Mathematics Pavia, ITALY
Preconditioned Mixed Spectral Element Methods for Elasticity and Stokes Problems Luca F. Pavarino Department of Mathematics Universita di Pavia, Via Abbiategrasso 209 27100 Pavia, ITALY pavarino@dragon.ian.pv.cnr.it
More informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
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 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 informationRobust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms
www.oeaw.ac.at Robust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms Y. Efendiev, J. Galvis, R. Lazarov, J. Willems RICAM-Report 2011-05 www.ricam.oeaw.ac.at
More informationNumerische Mathematik
Numer. Math. (1996) 75: 59 77 Numerische Mathematik c Springer-Verlag 1996 Electronic Edition A preconditioner for the h-p version of the finite element method in two dimensions Benqi Guo 1, and Weiming
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 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 informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 31, pp. 384-402, 2008. Copyright 2008,. ISSN 1068-9613. ETNA ON THE EQUIVALENCE OF PRIMAL AND DUAL SUBSTRUCTURING PRECONDITIONERS BEDŘICH SOUSEDÍK
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES On the approximation of a virtual coarse space for domain decomposition methods in two dimensions Juan Gabriel Calvo Preprint No. 31-2017 PRAHA 2017
More informationMultispace and Multilevel BDDC
Multispace and Multilevel BDDC Jan Mandel Bedřich Sousedík Clark R. Dohrmann February 11, 2018 arxiv:0712.3977v2 [math.na] 21 Jan 2008 Abstract BDDC method is the most advanced method from the Balancing
More informationANDREA TOSELLI. Abstract. Two-level overlapping Schwarz methods are considered for nite element problems
OVERLAPPING SCHWARZ METHODS FOR MAXWELL'S EQUATIONS IN THREE DIMENSIONS ANDREA TOSELLI Abstract. Two-level overlapping Schwarz methods are considered for nite element problems of 3D Maxwell's equations.
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 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 informationSharp condition number estimates for the symmetric 2-Lagrange multiplier method
Sharp condition number estimates for the symmetric -Lagrange multiplier method Stephen W. Drury and Sébastien Loisel Abstract Domain decomposition methods are used to find the numerical solution of large
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 informationMultipréconditionnement adaptatif pour les méthodes de décomposition de domaine. Nicole Spillane (CNRS, CMAP, École Polytechnique)
Multipréconditionnement adaptatif pour les méthodes de décomposition de domaine Nicole Spillane (CNRS, CMAP, École Polytechnique) C. Bovet (ONERA), P. Gosselet (ENS Cachan), A. Parret Fréaud (SafranTech),
More informationOptimal Interface Conditions for an Arbitrary Decomposition into Subdomains
Optimal Interface Conditions for an Arbitrary Decomposition into Subdomains Martin J. Gander and Felix Kwok Section de mathématiques, Université de Genève, Geneva CH-1211, Switzerland, Martin.Gander@unige.ch;
More informationAN ADAPTIVE CHOICE OF PRIMAL CONSTRAINTS FOR BDDC DOMAIN DECOMPOSITION ALGORITHMS
lectronic Transactions on Numerical Analysis. Volume 45, pp. 524 544, 2016. Copyright c 2016,. ISSN 1068 9613. TNA AN ADAPTIV CHOIC O PRIMAL CONSTRAINTS OR BDDC DOMAIN DCOMPOSITION ALGORITHMS JUAN G. CALVO
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 informationADDITIVE SCHWARZ FOR SCHUR COMPLEMENT 305 the parallel implementation of both preconditioners on distributed memory platforms, and compare their perfo
35 Additive Schwarz for the Schur Complement Method Luiz M. Carvalho and Luc Giraud 1 Introduction Domain decomposition methods for solving elliptic boundary problems have been receiving increasing attention
More informationA Nonoverlapping Subdomain Algorithm with Lagrange Multipliers and its Object Oriented Implementation for Interface Problems
Contemporary Mathematics Volume 8, 998 B 0-88-0988--03030- A Nonoverlapping Subdomain Algorithm with Lagrange Multipliers and its Object Oriented Implementation for Interface Problems Daoqi Yang. Introduction
More informationBETI for acoustic and electromagnetic scattering
BETI for acoustic and electromagnetic scattering O. Steinbach, M. Windisch Institut für Numerische Mathematik Technische Universität Graz Oberwolfach 18. Februar 2010 FWF-Project: Data-sparse Boundary
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 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 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 informationIntegration of Sequential Quadratic Programming and Domain Decomposition Methods for Nonlinear Optimal Control Problems
Integration of Sequential Quadratic Programming and Domain Decomposition Methods for Nonlinear Optimal Control Problems Matthias Heinkenschloss 1 and Denis Ridzal 2 1 Department of Computational and Applied
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 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 informationOptimal Left and Right Additive Schwarz Preconditioning for Minimal Residual Methods with Euclidean and Energy Norms
Optimal Left and Right Additive Schwarz Preconditioning for Minimal Residual Methods with Euclidean and Energy Norms Marcus Sarkis Worcester Polytechnic Inst., Mass. and IMPA, Rio de Janeiro and Daniel
More informationA LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES
A LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES P. HANSBO Department of Applied Mechanics, Chalmers University of Technology, S-4 96 Göteborg, Sweden E-mail: hansbo@solid.chalmers.se
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 posteriori error estimates for non conforming approximation of eigenvalue problems
A posteriori error estimates for non conforming approximation of eigenvalue problems E. Dari a, R. G. Durán b and C. Padra c, a Centro Atómico Bariloche, Comisión Nacional de Energía Atómica and CONICE,
More informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More informationAuxiliary space multigrid method for elliptic problems with highly varying coefficients
Auxiliary space multigrid method for elliptic problems with highly varying coefficients Johannes Kraus 1 and Maria Lymbery 2 1 Introduction The robust preconditioning of linear systems of algebraic equations
More informationSpectral element agglomerate AMGe
Spectral element agglomerate AMGe T. Chartier 1, R. Falgout 2, V. E. Henson 2, J. E. Jones 4, T. A. Manteuffel 3, S. F. McCormick 3, J. W. Ruge 3, and P. S. Vassilevski 2 1 Department of Mathematics, Davidson
More informationLuca F. Pavarino. Dipartimento di Matematica Pavia, Italy. Abstract
Domain Decomposition Algorithms for First-Order System Least Squares Methods Luca F. Pavarino Dipartimento di Matematica Universita dipavia Via Abbiategrasso 209 27100 Pavia, Italy pavarino@dragon.ian.pv.cnr.it.
More information