A uniform additive Schwarz preconditioner for high order Discontinuous Galerkin approximations of elliptic problems
|
|
- Tiffany Norman
- 5 years ago
- Views:
Transcription
1 Noname manuscript No. (will be inserted by the editor) A uniform additive Schwarz preconditioner for high order Discontinuous Galerkin approximations of elliptic problems Paola F. Antonietti Marco Sarti Marco Verani Ludmil T. Zikatanov Received: date / Accepted: date Abstract In this paper we design and analyze a uniform preconditioner for a class of high order Discontinuous Galerkin schemes. The preconditioner is based on a space splitting involving the high order conforming subspace and results from the interpretation of the problem as a nearly-singular problem. We show that the proposed preconditioner exhibits spectral bounds that are uniform with respect to the discretization parameters, i.e., the mesh size, the polynomial degree and the penalization coefficient. The theoretical estimates obtained are supported by numerical tests. Keywords discontinuous Galerkin method high order discretizations uniform preconditioning. P. F. Antonietti MOX-Laboratory for Modeling and Scientific Computing, Dipartimento di Matematica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milano, Italy paola.antonietti@polimi.it M. Sarti MOX-Laboratory for Modeling and Scientific Computing, Dipartimento di Matematica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milano, Italy marco.sarti@polimi.it M. Verani MOX-Laboratory for Modeling and Scientific Computing, Dipartimento di Matematica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milano, Italy marco.verani@polimi.it L. T. Zikatanov Department of Mathematics, The Pennsylvania State University, University Park, PA 16802, U.S.A. and Institute of Mathematics and Informatics, Bulgarian Academy of Sciences, Acad. G. Bonchev str, bl. 8, 1113 Sofia, Bulgaria ludmil@psu.edu
2 2 Paola F. Antonietti et al. 1 Introduction In the last years, the design of efficient solution techniques for the system of equations arising from Discontinuous Galerkin (DG) discretizations of elliptic partial differential equations has become an increasingly active field of research. On the one hand, DG methods are characterized by a great versatility in treating a variety of problems and handling, for instance, non-conforming grids and hp-adaptive strategies. On the other hand, the main drawback of DG methods is the larger number of degrees of freedom compared to (standard) conforming discretizations. In this respect, the case of high order DG schemes is particularly representative, since the corresponding linear system of equations is very ill-conditioned: it can be proved that, for elliptic problems, the spectral condition number of the resulting stiffness matrix grows like h 2 p 4, h and p being the granularity of the underlying mesh and the polynomial approximation degree, respectively, cf. [6]. As a consequence, the design of effective tools for the solution of the linear system of equations arising from high order DG discretizations becomes particularly challenging. In the context of elliptic problems, Schwarz methods for low order DG schemes have been studied in [29], where overlapping and non-overlapping domain decomposition preconditioners are considered, and bounds of O(H/δ) and O(H/h), respectively, are obtained for the condition number of the preconditioned operator. Here H, h and δ stand for the granularity of the coarse and fine grids and the size of the overlap, respectively. Further extensions including inexact local solvers, and the extension of two-level Schwarz methods to advection-diffusion and fourth-order problems can be found in [38, 30, 1 3, 26, 12, 5]. In the field of Balancing Domain Decomposition (BDD) methods, a number of results exist in literature: exploiting a Neumann-Neumann type method, in [23, 24] a conforming discretization is used on each subdomain combined with interior penalty method on non-conforming boundaries, thus obtaining a bound for the condition number of the resulting preconditioner of O((1 log(h/h)) 2 ). In [21], using the unified framework of [11] a BDDC method is designed and analyzed for a wide range of DG methods. The auxiliary space method (ASM) (see e.g., [41,32,49,33]) is employed in the context of h-version DG methods to develop, for instance, the two-level preconditioners of [22] and the multilevel method of [15]. In both cases a stable splitting for the linear DG space is provided by a decomposition consisting of a conforming subspace and a correction, thus obtaining uniformly bounded preconditioners with respect to the mesh size. All the previous results focus on low order (i.e., linear) DG methods. In the context of preconditioning high order DG methods we mention [6, 8], where a class of non-overlapping Schwarz preconditioners is introduced, and [4], where a quasi-optimal (with respect to h and p) preconditioner is designed in the framework of substructuring methods for hp-nitsche-type discretizations. Multigrid methods for high order DG discretizations have been analyzed in [9], cf. also [7] for the extension to genereal polygonal/polyhedral meshes. A study of a BDDC scheme in the case of hp-spectral DG methods is addressed
3 A uniform additive Schwarz preconditioner for high order DG methods 3 in [19], where the DG framework is reduced to the conforming one via the ASM. The ASM framework is employed also in [14], where the high order conforming space is employed as auxiliary subspace, and a uniform multilevel preconditioner is designed for hp-dg spectral element methods in the case of locally varying polynomial degree. To the best of our knowledge, this preconditioner is the only uniform preconditioner designed for high order DG discretizations. We note that, in the framework of high order methods, the decomposition involving a conforming subspace was already employed in the case of a-posteriori error analysis, see for example [35, 17, 51]. In this paper, we address the issue of preconditioning high order DG methods by exploiting this kind of space splitting based on a high order conforming space and a correction. However, in our case the space decomposition is suggested by the interpretation of the high order DG scheme in terms of a nearly-singular problem, cf. [39]. Even though the space decomposition is similar to that of [14], the preconditioner and the analysis we present differs considerably since here we employ the abstract framework of subspace correction methods provided by [50]. More precisely, we are able to show that a simple pointwise Jacobi method paired with an overlapping additive Schwarz method for the conforming subspace, gives uniform convergence with respect to all the discretization parameters, i.e., the mesh size, the polynomial order and the penalization coefficient appearing in the DG bilinear form. The rest of the paper is organized as follows. In Section 2, we introduce the model problem and the corresponding discretization through a class of symmetric DG schemes. Section 3 is devoted to few auxiliary results regarding the Gauss-Legendre-Lobatto nodes, whose properties are fundamental to prove the stability of the space decomposition proposed in Section 4. The analysis of the preconditioner is presented in Section 5 and the theoretical results are supported by the numerical simulations of Section 6. We also test the performance of the proposed method in the case of more general diffusion equation, with isotropic discontinuous as well as anisotropic diffusion tensors. 2 Model problem and high order DG discretization In this section we introduce the model problem and its discretization through several Discontinuous Galerkin schemes, see also [11]. Throughout the paper, we use the notation x y and x y to denote the inequalities x Cy and x Cy, respectively, C being a positive constant independent of the discretization parameters. Moreover, x y means that there exist constants C 1, C 2 > 0 such that C 1 y x C 2 y. When needed, the constants are written explicitly. Given a convex polygonal/polyhedral domain Ω R d, d = 2, 3, and f L 2 (Ω), we consider the following weak formulation of the Poisson problem with ho-
4 4 Paola F. Antonietti et al. mogeneous Dirichlet boundary conditions: find u V := H0 1 (Ω), such that u v dx = fv dx v V. (1) Ω Ω Let T h denote a shape-regular, conforming, locally quasi-uniform partition of Ω into shape-regular elements κ of diameter h κ, and set h := max κ Th h κ. We also assume that each element κ T h results from the mapping, through an affine operator F κ, of a reference element ˆκ, which is the open, unit d-hypercube in R d, d = 2, 3. We denote by Fh I and F h B the set of internal and boundary faces (for d = 2 face means edge ) of T h, respectively, and define F h := Fh I F h B. We associate to any F F h a unit vector n F orthogonal to the face itself and also denote by n F,κ the outward normal vector to F κ with respect to κ. We observe that for any F Fh B, n F,κ = n F, since F belongs to a unique element. For any F Fh I, we assume F = κ+ κ, where κ + := {κ T h : F κ, n F n F,κ > 0}, κ := {κ T h : F κ, n F n F,κ < 0}. For regular enough vector-valued and scalar functions τ and v, we denote by τ ± and v ± the corresponding traces taken from the interior of κ ±, respectively, and define the jumps and averages across the face F Fh I as follows τ := τ + n F,κ + + τ n F,κ, v := v + n F,κ + + v n F,κ, {τ } := τ + + τ, 2 {v } := v+ + v, 2 For F F B h, the previous definitions reduce to v := vn F and {τ } := τ. We now associate to the partition T h, the high order Discontinuous Galerkin finite element space V hp defined as V hp := {v L 2 (Ω) : v F κ Q p (ˆκ) κ T h }, with Q p denoting the space of all tensor-product polynomials on ˆκ of degree p > 1 in each coordinate direction. We define the lifting operators R(τ ) := F F h r F (τ ) and L(v) := F F l h I F (v), where r F : [L 2 (F )] d [V hp ] d, r F (τ ) η dx := τ {η } ds F F h. Ω F l F : L 2 (F ) [V hp ] d, l F (v) η dx := v η ds F Fh, I Ω for any η [V hp ] d. We then introduce the DG finite element formulation: find u V hp such that A(u, v) = fv dx v V hp, (2) Ω F
5 A uniform additive Schwarz preconditioner for high order DG methods 5 with A(, ) : V hp V hp R defined as A(u, v) := u v dx + κ T h κ κ T h κ κ T h κ u (R( v ) + L(β v )) dx + (R( u ) + L(β u )) v dx + + θ Ω F F h (R( u ) + L(β u )) (R( v ) + L(β v )) dx, F σ u v ds(3) where θ = 0 for the SIPG method of [10] and θ = 1 for the LDG method of [20]. With regard to the vector function β, we have β = 0 for the SIPG method, while β R d is a uniformly bounded (and possibly null) vector for the LDG method. The penalization function σ L (F h ) is defined as p 2 σ F := α min(h κ +, h κ ), F F h, I σ F := α p2 F Fh B, h κ being α 1 and h κ ± the diameters of the neighboring elements κ ± T h sharing the face F Fh I. We endow the DG space V hp with the following norm v 2 DG := v 2 L 2 (κ) + σ 1/2 v 2 L 2 (F ), κ T h F F h and state the following result, cf. [36,44,6,45]. Lemma 1 The following results hold A(u, v) u DG v DG u, v V hp, A(u, u) u 2 DG u V hp. (4) For the SIPG formulation coercivity holds provided the penalization coefficient α is chosen large enough. From Lemma 1 and using the Poincaré inequality for piecewise H 1 functions of [13], the following spectral bounds hold, cf. [6]. Lemma 2 For any u V hp it holds that κ T h u 2 L 2 (κ) A(u, u) α p4 h 2 u 2 L 2 (κ). (5) κ T κ h 3 Gauss-Legendre-Lobatto nodes and quadrature rule In this section we provide some details regarding the choice of the basis functions spanning the space V hp and the corresponding degrees of freedom. On the reference d-hypercube [ 1, 1] d, we choose the basis obtained by the
6 6 Paola F. Antonietti et al. tensor product of the one-dimensional Lagrange polynomials on the reference interval [ 1, 1], based on Gauss-Legendre-Lobatto (GLL) nodes. We denote by N I (ˆκ) (N B (ˆκ)) the set of interior (boundary) nodes of ˆκ, and define N (ˆκ) := N I (ˆκ) N B (ˆκ). The analogous sets in the physical frame are denoted by N I (κ), N B (κ) and N (κ), where any ξ p N (κ) is obtained by applying the linear mapping F κ : ˆκ κ to the corresponding ˆξ p N (ˆκ). The choice of GLL points as degrees of freedom allow us to exploit the properties of the associated quadrature rule. We recall that, given (p + 1) d GLL quadrature nodes {ˆξ p } and weights {ŵ ξp }, we have ˆξ p N (ˆκ) v(ˆξ p )ŵ ξp = ˆκ v dx v Q 2p 1 (ˆκ), which implies that v(ˆξ p ) 2 ŵ ξp ˆξ p N (ˆκ) ˆκ v 2 dx v Q p (ˆκ). However, by defining, for v Q p (ˆκ), the following norm it can be proved that v 2 0,p,ˆκ := ξ p N (ˆκ) v(ˆξ p ) 2 ŵ ξp, v 2 0,p,ˆκ v 2 L 2 (ˆκ), (6) cf. [18, Section 5.3]. The same result holds for the physical frame κ, i.e., v 2 0,p,κ v 2 L 2 (κ). Considering the Lagrange basis {φ ξp }, ξ p κ T h N (κ), we can write any v V hp as v = v(ξ p )φ ξp = v ξp, (7) κ T h ξ p N (κ) where we note that v ξp = v(ξ p )φ ξp. κ T h ξ p N (κ) Lemma 3 For any v V hp, given the decomposition (7), the following equivalence holds v 2 L 2 (Ω) v ξp 2 L 2 (κ). κ T h ξ p N (κ) Proof The proof can be restricted to the case of a single element κ T h. We write v V hp as in (7), and observe that v ξp 2 0,p,κ = ξ p N (κ) v ξp (ξ p) 2 w ξp = v ξp (ξ p ) 2 w ξp,
7 A uniform additive Schwarz preconditioner for high order DG methods 7 hence, by (6), v 2 L 2 (κ) v(ξ p ) 2 w ξp = v ξp (ξ p ) 2 w ξp ξ p N (κ) = v ξp 2 0,p,κ ξ p N (κ) ξ p N (κ) ξ p N (κ) and the thesis follows summing over all κ T h. v ξp 2 L 2 (κ), 4 Space decomposition for high order DG methods The design of our preconditioner is based on a two-stage space decomposition: we first split the high order DG space as V hp = Vhp B + V hp C, with V hp B denoting a proper subspace of V hp, to be defined later, and Vhp C denoting the high order conforming subspace. As a second step, both spaces are further decomposed to build two corresponding additive Schwarz methods in each of the subspaces. The final preconditioner on V hp is then obtained by combining the two subspace preconditioners. The first space splitting is suggested by the interpretation of the high order DG formulation (2) as a nearly-singular problem. To present the motivation behind this choice, we briefly introduce the theoretical framework of [39] regarding space decomposition methods for this class of equations. Given a finite dimensional Hilbert space V, we consider the following problem: find U V such that AU = (A 0 + ɛa 1 )U = F, (8) where F is a given vector and where A 0 is symmetric and positive semi-definite and A 1 is symmetric and positive definite. As a consequence, if ɛ = 0, the problem is singular, but here we are interested in the case ɛ > 0 (with ɛ small, ɛ 1), i.e., (8) is nearly-singular. In general, the conditioning of problem (8) degenerates for decreasing ɛ, and this affects the performance of standard preconditioned iterative methods, unless proper initial guess are chosen. In the framework of space decomposition methods, in order to obtain a ɛ-uniform preconditioner, a key assumption on the space splitting V hp = N i=1 V i is needed. Assumption 1 ([39]) The decomposition V hp = N i=1 V i satisfies ker(a 0 ) = where ker(a 0 ) is the kernel of A 0. N (V i ker(a 0 )), i=1 We now turn to our DG framework, and show that a high order DG formulation can be indeed read as a nearly-singular problem with a suitable choice of ɛ. For the sake of simplicity, and without any loss of generality, we retrieve
8 8 Paola F. Antonietti et al. equation (8) working directly on a bilinear form that is spectrally equivalent to A(, ). To this aim, let the bilinear forms A (, ), A J (, ) and Ã(, ) be defined as A (u, v) := u v dx, κ T h A J (u, v) := F F h κ F u v ds, Ã(u, v) := A (u, v) + α p2 h A J(u, v), and let A, A J, and à be their corresponding operators. Clearly, A and A J are both symmetric and positive semi-definite, and à is symmetric and positive definite. Moreover, thanks to Lemma 1, and the local quasi-uniformity of the partition, the following spectral equivalence result holds A(u, u) u 2 DG Ã(u, u). From the spectral equivalence shown above, it is immediate to conclude that any optimal preconditioner for à is automatically an optimal preconditioner for A. Therefore, in the following we focus our attention on the following problem Ãũ = (A + 1 ɛ A J)ũ = f, (9) with ɛ := h/(αp 2 ) < 1. After some simple calculations, we can write (9) as [ɛ(a + A J ) + (1 ɛ)a J ] ũ = ɛf, which corresponds to (8) with A 1 = A + A J and A 0 = (1 ɛ)a J. In order to obtain a suitable space splitting satisfying Assumption 1, we observe that, according to the definition above, the kernel of A 0 is given by the space of continuous polynomial functions of degree p vanishing on the boundary Ω. We then derive the first space decomposition with V B hp := {v V hp : v(ξ p ) = 0 V hp = V B hp + V C hp, (10) ξ p V C hp := {v C 0 (Ω) : v F κ Q p (ˆκ) κ T h N I (κ)}, κ T h, v Ω = 0} H 1 0 (Ω), i.e., Vhp B consists of the functions in V hp that are null in any degree of freedom in the interior of any κ T h. Moreover, we observe that Vhp B V hp, and Vhp B V hp C V hp C, hence Assumption 1 is satisfied by decomposition (10), which is the basis to develop the analysis of our preconditioner for problem (2).
9 A uniform additive Schwarz preconditioner for high order DG methods Technical results In this subsection we present several results, which are fundamental for the forthcoming analysis. We introduce a suitable interpolation operator Q h : V hp V C hp, consisting of the Oswald operator, cf. [34,37,28,16,17]. For any v V hp, we can define on each κ T h the action of the operator Q h, by prescribing the value of Q h v in any ξ p N (κ): 0 if ξ p Ω, Q h v(ξ p ) := 1 v κ (ξ p ) otherwise, card(t ξp ) κ T ξp (11) with T ξp := {κ T h : ξ p κ }. Note that from the above definition it follows that v Q h v Vhp B, for any v V hp. In addition to the space of polynomials Q p (κ), we define Q p 0 (κ) as Q p 0 (κ) := {v Qp (κ) : v(ξ p ) = 0 ξ p N I (κ)}, and state the following trace and inverse trace inequalities. Lemma 4 ([17, Lemma 3.1]) The following trace and inverse trace inequalities hold v 2 L 2 ( κ) p2 h κ v 2 L 2 (κ) v Qp (κ), v 2 L 2 (κ) h κ p 2 v 2 L 2 ( κ) v Qp 0 (κ). (12) The next result is a keypoint for the forthcoming analysis, and can be found in [34] and [17, Lemma 3.2]. Lemma 5 ([34] and [17, Lemma 3.2]) For any v V hp, the following estimate holds v Q h v 2 L 2 (κ) h κ p 2 v 2 L 2 (F ), with F h (κ) := {F F h : F κ }. F F h (κ) Thanks to Lemma 5 we can prove the following theorem. Theorem 1 For any v V hp, it holds that A(v Q h v, v Q h v) + A(Q h v, Q h v) A(v, v), (13) where Q h v Vhp C in (10) is stable. is defined as in (11). Then the space decomposition defined
10 10 Paola F. Antonietti et al. Proof We observe that, from (5), the quasi-uniformity of the mesh and Lemma 5, we obtain A(v Q h v, v Q h v) α p4 h 2 κ T κ h v Q h v 2 L 2 (κ) α σ 1/2 v 2 L 2 (F ) A(v, v). F F h p4 h 2 κ T κ h h κ p 2 F F h (κ) v 2 L 2 (F ) The upper bound (13) follows from the triangle inequality and the above estimate A(Q h v, Q h v) A(v Q h v, v Q h v) + A(v, v) A(v, v). For any v V hp, we recall that v Q h v Vhp B, which implies inf v B V B hp,vc V C hp v B +v C =v A(v B, v B ) + A(v C, v C ) A(v Q h v, v Q h v) + A(Q h v, Q h v) A(v, v). 5 Construction and analysis of the preconditioner In this section we introduce our preconditioner and analyze the condition number of the preconditioned system. Employing the nomencalture of [48], the preconditioner is a parallel subspace correction method (also known as additive Schwarz preconditioner, see. e.g., [40, 47, 27]). Our construction uses a decomposition in two subspaces, cf. (10), and inexact subspace solvers. Each of the subspace solvers is a parallel subspace correction method itself. 5.1 Canonical representation of a parallel subspace correction method The main ingredients needed for the analysis of the parallel subspace correction (PSC) preconditioners are suitable space splittings and the corresponding subspace solvers (see [40,47,27,48,31,50,46]). In our analysis we use the notation and the general setting from [50]. We have the following abstract result. Lemma 6 ([50, Lemma 2.4]) Let V be a Hilbert space which is decomposed as V = N i=1 V i, V i V, i = 1,..., N, and T i : V V i, i = 1,..., N be operators whose restrictions on V i are symmetric and positive definite. For T := N i=1 T i the following identity holds A(T 1 v, v) = inf v i V i vi=v N i=1 A(T 1 i v i, v i ). (14)
11 A uniform additive Schwarz preconditioner for high order DG methods 11 According to the above lemma, to show a bound on the condition number of the preconditioned system we need to show that there exist positive constants c and C such that ca(v, v) A(T 1 v, v) CA(v, v). Remark 1 In many cases we have T i = P i, i = 1,..., N, where P i : V V i are the elliptic projections defined as follows: for v V, its projection P i v is the unique element of V i satisfying A(P i v, v i ) := A(v, v i ), for all v i V i. Note that by definition, P i is the identity on V i, namely, P i v i = v i = P 1 i v i V i. Hence, for T = N i=1 P i, the relation (14) gives A(T 1 v, v) = inf v i V i vi=v v i, for all N A(v i, v i ). (15) i=1 5.2 Space splitting and subspace solvers To fix the notation, let us point out that in what follows we use T (with subscript when necessary) to denote (sub)space solvers and preconditioners. Accordingly, P with subscript or superscript denotes elliptic projection on the corresponding subspace. We now define the space splitting and the corresponding subspace solvers. We recall the space decomposition from Section 4, V hp = Vhp B + V hp C, where V hp B are all functions in V hp for which the degrees of freedom in the interior of any κ T h vanish, and Vhp C is the space of high order continuous polynomials vanishing on Ω. Note that Vhp B V hp C {0}, and that V hp B contains non-smooth and oscillatory functions, while Vhp C contains the smooth part of the space V hp. Next, on each of these subspaces we define approximate solvers T B : V hp Vhp B and T C : V hp Vhp C. First, we decompose Vhp B as follows where V ξp := { V B hp = κ T h ξ p N B(κ) v V B hp : v(ξ p ) = 0 for any ξ p V ξp, (16) ( κ T h N B (κ) ) \ {ξ p } The approximate solver on V B then is a simple Jacobi method, defined as T B : V hp Vhp, B T B := P ξp P B. κ T h ξ p N B(κ) }.
12 12 Paola F. Antonietti et al. where P B and P ξp are the elliptic projections on Vhp B and V ξp, respectively. Note that T B is defined on all of V hp and is also an isomorphism when restricted to Vhp B, because the elliptic projection P B and P ξp are the identity on Vhp B and V ξp, respectively. In addition, the splitting is a direct sum, and, hence, any v Vhp B is uniquely represented as v = κ T h ξ p N B(κ) vξp, v ξp V ξp. Then, taking P i = P ξp P B : V hp V ξp, from (15), we have A(T 1 B vb, v B ) = κ T h ξ p N B(κ) A(v ξp, v ξp ), v B V B hp. (17) Next, we introduce the preconditioner T C on Vhp C. This is the two-level overlapping additive Schwarz method introduced in [42, 43] for high order conforming discretizations. If we denote by N V the number of interior vertices of T h, then this preconditioner corresponds to the following decomposition of Vhp C : N V Vhp C = Vi C. (18) i=0 Here V0 C is the (coarse) space of continuous piecewise linear functions on T h, and for i = 1,..., N V, Vi C := Vhp C H1 0 (Ω i ), where Ω i is the union of the elements sharing the i-th vertex (see Fig. 1 for a two-dimensional example). We recall that, in the case of Neumann and mixed boundary conditions, in order to obtain a uniform preconditioner, the decomposition (18) should be enriched with the subdomains associated to those vertices not lying on a Dirichlet boundary, see [42, 43] for details. Ω i Ω i Fig. 1: Examples of subdomains in a two-dimensional setting. Then, for any Vi C, i = 0,..., N V, we denote by Pi C : Vhp C V i C the elliptic projections on Vi C and define the two-level overlapping additive Schwarz operator as [ N V T C : V hp Vhp, C T C := P0 C + i=1 P C i ] P C = (P 0 + P V )P C, (19)
13 A uniform additive Schwarz preconditioner for high order DG methods 13 where P C is the elliptic projection on Vhp C. As in the case of T B, we have that the restriction of T C on Vhp C is an isomorphism. In addition, from (15) with P i = Pi CP C : V hp Vi C, we have A(T 1 C v, v) = inf v i V C i vi=v N V A(v i, v i ). (20) i=0 Observe that, especially for three dimensional problems, the dimension of the coarse space can be quite large, for example when the underlying mesh is very fine. In these cases, the solution of the corresponding coarse subproblem through a direct method can be computationally unfeasible. However, since the coarse level is the piece-wise linear conforming subspace, the associated linear system of equations can be efficiently solved by further preconditioning it with one of the efficient techniques already available for this class of problems, for example (Algebraic) Multigrid methods. 5.3 Definition of the global preconditioner Finally, we define the global preconditioner on V hp by setting T DG : V hp V hp, T DG := T B + T C, (21) We remark that from Lemma 6, with N = 2, T 1 = T B, V 1 = Vhp B, T 2 = T C, V 2 = Vhp C, we have A(T 1 DGv, v) = inf v B V B hp,vc V C hp v B +v C =v [ A(T 1 B vb, v B ) + A(T 1 C vc, v C ) ]. (22) 5.4 Condition number estimates: subspace solvers We now show the estimates on the conditioning of the subspace solvers needed to bound the condition number of T DG. The first result that we prove is on the conditioning of T B. Lemma 7 Let T B denote the Jacobi preconditioner defined in (5.2). Then there exist two positive constants C J 1 and C J 2, independent of the granularity of the mesh h, the polynomial approximation degree p and the penalization coefficient α, such that A(T 1 B vb, v B ) C J 1 A(v B, v B ) v B Vhp B (23) A(T 1 B (v Q hv), v Q h v) C J 2 A(v Q h v, v Q h v) v V hp, (24) with Q h v defined in (11).
14 14 Paola F. Antonietti et al. Proof We refer to the space decomposition (16) and write v B = v ξp. κ T h ξ p N B(κ) For the lower bound (23), we employ the eigenvalue estimate (5) and Lemma 3, thus obtaining A(v B, v B ) α p4 h 2 κ T κ h v B 2 L 2 (κ) We now observe that for any ξ p N B (κ), v ξp the inverse trace inequality (12) to obtain A(v B, v B ) α p4 h 2 κ T κ h ξ p N B(κ) v ξp 2 L 2 (κ) α p4 h 2 κ T κ h ξ p N B(κ) v ξp 2 L 2 (κ). Q p 0 (κ), and we can thus apply α p2 h κ κ T h ξ p N B(κ) Noting that v ξp 2 L 2 ( κ) = vξp 2 L 2 ( κ), it follows that A(v B, v B ) α p2 h κ κ T h ξ p N B(κ) κ T h ξ p N B(κ) v ξp 2 DG, v ξp 2 L 2 ( κ) κ T h ξ p N B(κ) v ξp 2 L 2 ( κ). σ 1/2 v ξp 2 L 2 ( κ) and the thesis follows from the coercivity bound (4) and (17). With regard to the upper bound (24), for the sake of simplicity we denote w = (I Q h )v, and observe that w = (I Q h )w. Since w Vhp B, we write w = w ξp, κ T h ξ p N B(κ) and, from (17), A(T 1 w, w) = B κ T h ξ p N B(κ) A(w ξp, w ξp ). Applying again the estimate (5) and Lemma 3, we obtain κ T h ξ p N B(κ) A(w ξp, w ξp ) κ T h ξ p N B(κ) α p4 h 2 κ T κ h α p4 h 2 w ξp 2 L 2 (κ) α p4 κ h 2 w 2 L 2 (κ) κ T κ h (I Q h )w 2 L 2 (κ) A(w, w), where the last steps follows from Lemma 5 and the quasi-uniformity of the mesh.
15 A uniform additive Schwarz preconditioner for high order DG methods 15 For the analysis of the additive preconditioner T C given in (19), we need several preliminary results (see [42, 43] for additional details). First of all, given the decomposition N V v = v 0 + i=1 v i v V C hp, v 0 V C 0, v i V C i, (25) we define the coarse function v 0 as the L 2 -projection on the space V C 0, i.e., v 0 := I 0 v with I 0 v satisfying v I 0 v 2 L 2 (Ω) h2 v 2 H 1 (Ω), (26) I 0 v 2 H 1 (Ω) v 2 H 1 (Ω), (27) for any v H 1 0 (Ω). For any i = 1,..., N V, the functions v i appearing in (25) are defined as v i := I p (θ i (v v 0 )), where θ i is a proper partition of unity and I p is an interpolation operator, described in the following. For any Ω i, i = 1,..., N V, the partition of unity θ i is such that θ i V C h1 and it can be defined by prescribing its values at the vertices {v} belonging to Ω i, and imposing it to be zero on Ω \ Ω i, see Fig. 2 for d = 2. More precisely, θ i (v) = { 1 if v is the internal vertex or Fv F B h, 0 otherwise, with F v := {F F h, F Ω i : v F } Ω Ω Fig. 2: Values of the partition of unity θ i for d = 2. It follows that: supp(θ i ) = Ω i, 0 θ i 1, N V θ i = 1, θ i 1 h. (28) i=1
16 16 Paola F. Antonietti et al. As interpolation operator I p, we make use of the operator defined in [42,43]: setting z := v v 0, we define I p (θ i z)(ξ p ) = (θ i z)(ξ p ) ξ p N (κ), κ Ω i. (29) Notice that, despite defined locally, I p (θ i z) belongs to Vi C since the interelement continuity is guaranteed by the fact that the (p + 1) d 1 GLL points on a face uniquely determine a tensor product polynomial of degree p defined on that face. The following result holds. Lemma 8 ([42, Lemma 3.1, Lemma 3.3]) The interpolation operator I p : Q p+1 (ˆκ) Q p (ˆκ), defined in (29), is bounded uniformly in the H 1 seminorm, i.e., I p (u) H 1 (ˆκ) u H 1 (ˆκ) u Q p+1 (ˆκ). (30) Once the partition of unity and the interpolation operator are defined, we are able to complete the analysis of T C. In analogy to Lemma 7, which is based on (17), we now use (20) and the above auxiliary results to show the following lemma. Lemma 9 Let T C denote the two-level overlapping additive Schwarz preconditioner defined in (19). Then there exist two positive constants C C 1 and C C 2, independent of the discretization parameters, i.e., the granularity of the mesh h and the polynomial approximation degree p, such that for any v V C hp. A(T 1 C v, v) CC 1 A(v, v) (31) A(T 1 C v, v) CC 2 A(v, v), (32) Proof We first prove the lower bound (31), and given the decomposition (25), we can write A(v, v) = N V i,j=0 A(v i, v j ) A(v 0, v 0 ) + N V i,j=1 A(v i, v j ). We now note that A(v i, v j ) 0 only if i = j and Ω i Ω j, and since each Ω i is overlapped by a limited number of neighboring subdomains, we conclude that A(v, v) A(v 0, v 0 ) + N V i,j=1 N V A(v i, v j ) A(v 0, v 0 ) + A(v i, v i ). Inequality (31) follows from the bound above and (20), denoting with C C 1 the hidden constant. Note that, from (20), the upper bound (32) is proved provided the following inequality holds N V A(v i, v i ) C C 2 A(v, v) v Vhp. C (33) i=0 i=1
17 A uniform additive Schwarz preconditioner for high order DG methods 17 We recall that v 0 = I 0 v, and from (27) it follows that A(v 0, v 0 ) = A(I 0 v, I 0 v) A(v, v). (34) For i = 1,..., N V, we have v i = I p (θ i z), with z = v v 0, and by (30), we obtain d 2 v i 2 H 1 (κ ) θ iz 2 H 1 (κ ) θ i z z + θ i x j x j, for any κ Ω i. By (28) it holds that j=1 L 2 (κ ) θ i 1 h, θ i L 1, hence, v i 2 H 1 (κ ) 1 h 2 z 2 L 2 (κ ) + d j=1 z x j 2 L 2 (κ ) 1 h 2 v v 0 2 L 2 (κ ) + v v 0 2 H 1 (κ ). On any element κ, a limited number of components v i are different from zero (at most four for d = 2, and eight for d = 3), which implies that we can sum over all the components v i, i = 1,..., N V, and then over all the elements, thus obtaining N V v i 2 H 1 (Ω) 1 h 2 v v 0 2 L 2 (Ω) + v v 0 2 H 1 (Ω) v 2 H 1 (Ω), i=1 where the last step follows from (26) and (27). The addition of the above result and (34), gives (33), denoting with C C 2 the resulting hidden constant. 5.5 Condition number estimates: global preconditioner We are now ready to prove the main result of the paper regarding the condition number of the preconditioned problem. Theorem 2 Let T DG be defined as in (21). Then, for any v V hp, it holds that A(v, v) A(T 1 DGv, v) A(v, v), (35) where the hidden constants are independent of the discretization parameters, i.e., the mesh size h, the polynomial approximation degree p, and the penalization coefficient α.
18 18 Paola F. Antonietti et al. Before proving Theorem 2, we note the following: As the operator T DG is symmetric and positive definite, from Theorem 2 it immediately follows that the spectral condition number of the preconditioned operator T DG κ(t DG ) := λ max(t DG ) λ min (T DG ) (36) is uniformly bounded, i.e. κ(t DG ) 1, where the hidden constant is independent of the the mesh size h, the polynomial approximation degree p, and the penalization coefficient α. In (36) λ max (T DG ) and λ min (T DG ) denotes the extremal eigenvalues of T DG defined as A(T DG v, v) A(T DG v, v) λ max (T DG ) := sup, λ min (T DG ) := inf. 0 v V hp A(v, v) 0 v V hp A(v, v) As a consequence, the conjugate gradient algorithm can be employed to solve the preconditioned linear system of equations and the number of iterations needed to reduce the norm of the (relative) residual below a given threshold is expected to be uniformly constant independently of the size of the underlying matrix, see, e.g., [46] for more details. Proof (Proof of Theorem 2). To prove the upper bound, we first consider the identity (22). Recalling that, by definition (11), v Q h v V B hp, for any v V hp, we obtain A(T 1 DGv, v) = inf v B V B hp,vc V C hp v B +v C =v A(T 1 B [ A(T 1 B vb, v B ) + A(T 1 C vc, v C ) ] (v Q hv), v Q h v) + A(T 1 Q hv, Q h v). From the bounds (24) and (32) for Q h v, it follows that C A(T 1 DG 1 v, v) A(TB (v Q hv), v Q h v) + A(T 1 C Q hv, Q h v) A(v Q h v, v Q h v) + A(Q h v, Q h v) A(v, v), where the last step follows from (13). The lower bound follows from (22), the bounds (23) and (31), and a triangle inequality A(T 1 DGv, v) = inf v B V B hp,vc V C hp v B +v C =v inf v B V B hp,vc V C hp v B +v C =v [ A(T 1 B vb, v B ) + A(T 1 C vc, v C ) ] [ A(v B, v B ) + A(v C, v C ) ] A(v, v).
19 A uniform additive Schwarz preconditioner for high order DG methods 19 6 Numerical experiments In this section we present some numerical tests to verify the theoretical estimates provided in Lemma 7, Lemma 9 and Theorem 2, as well as to test the performance of our preconditioner in the case of a variable diffusion tensor. To this aim we present several tests with anisotropies and/or heterogeneities in the diffusion coefficients following the lines of those proposed in [25]. As these tests show, the preconditioner is robust for problems with tensor coefficients, even when strong anisotropies are presented, as long as the variation of these coefficients throughout the domain is bounded. As expected, a deterioration in the condition number of the preconditioned system is observed when the coefficients have large jumps. This behavior is due to the use of a standard conforming coarse space. More complex constructions of coarse spaces depending on the coefficient variation are possible, but they are beyond the scope of our considerations here. Before presenting the numerical results assessing the practical performance of our preconditioner, we briefly provide some implementation details. The preconditioner is composed by the sum of two operators: i) a (block) Jacobi preconditioner acting on the degrees of freedom (dofs) of the DG space with the exclusion of the dofs associated to the interior of the elements (i.e., the bubble modes) and ii) a classical two-level overlapping Schwarz preconditioner acting on the space of high-order continuous polynomials vanishing on Ω. Implementing the action of the first operator is straightforward as it requires only a multiplication by a (block) diagonal matrix associated to the vertex/edge/face dofs. Regarding the action of the second operator, associated to the smooth part of the space V hp, it is easily seen that the computationally demanding part is the construction of the projector from the discontinuous space V hp to the conforming space. For example, in two-dimensions, it consists in constructing a matrix that averages the values at the vertex and edge dofs. Clearly, the projection matrix has a highly sparse structure. Once the conforming stiffness matrix is available, it is enough to apply the classical overlapping Schwarz preconditioner of [42, 43], with coarse space from piecewise continuous linear polynomials, for which efficient implementations are available. Finally, let us briefly comment on the actual cost of the proposed preconditioner relative to other strategies. Since the proposed approach is based on an overlapping partition of the computational domain, each application of the preconditioner is more expensive than the Schwarz preconditioners proposed in [6,8,1,2], that feature a high-level of scalability and load balancing because they are built using a non-overlapping approach. On the other hand, our preconditioner is uniform with respect to the polynomial approximation degree p whereas the preconditioner of [6, 8] exhibit spectral bounds that depend on p.
20 20 Paola F. Antonietti et al. 6.1 Example 1 We consider problem (2) in the two dimensional case with Ω = ( 1, 1) 2 and SIPG and LDG discretizations. For the first experiment, we set h = , the penalization parameter α = 10 and β = 1 for the LDG method. In Table 1, we show the numerical evaluation of the constants C J 1 and C J 2 of Lemma 7 and C C 1 and C C 2 of Lemma 9, as a function of the polynomial order employed in the discretization: the constants are independent of p, as expected from theory. With regard to the constants C C 1 and C C 2, we observe that the values are the same for both the SIPG and LDG methods, since the preconditioner on the conforming subspace reduces to the same operator regardless of the DG scheme employed. Table 1: Left and middle: numerical evaluation of the constants C J 1 and C J 2 of Lemma 7 as a function of p for the SIPG and LDG methods; right: numerical evaluation of the constants C C 1 and C C 2 of Lemma 9 as a function of p SIPG (α = 10, β = 0) LDG (α = 10, β = 1) C J 1 C J 2 C J 1 C J 2 C C 1 C C 2 p = p = p = p = p = Table 2 shows a comparison of the spectral condition number of the original system (K(A)) and of the preconditioned one (K(T DG )). While the former grows as p 4, cf. [6], the latter is constant with p, as stated in (35). The theoretical results are further confirmed by the number of iterations Niter P CG and of the Preconditioned Conjugate Gradient (PCG) and the Conjugate N CG iter Gradient (CG), respectively, needed to reduce the initial relative residual of a factor of Table 2: Condition number of the unpreconditioned (K(A)) and preconditioned (K(T DG )) linear systems of equations and corresponding CG (Niter CG ) and PCG ) iteration counts as a function of p for the SIPG and LDG methods (N P CG iter SIPG (α = 10, β = 0) LDG (α = 10, β = 1) K(A) Niter CG K(T DG ) Niter P CG K(A) Niter CG K(T DG ) Niter P CG p = p = p = p = p =
21 A uniform additive Schwarz preconditioner for high order DG methods 21 In the second numerical experiment, we consider the same test case presented above, but we now fix the polynomial approximation degree p = 6 and decrease the mesh-size h. The computed results obtained with the SIPG (β = 0, α = 10) and LDG (β = 1, α = 10) methods are shown in Table 3. For the sake of comparison, Table 3 also shows the computed spectral condition numbers of the unpreconditioned system as well as the number of CG iterations needed to reduced the relative residual below the given tolerance As predicted from our theoretical results, the condition number of the preconditioned system is insensitive to the the mesh-size, whereas K(A) grows quadratically as h tends to zero. Table 3: Condition number of the unpreconditioned (K(A)) and preconditioned (K(T DG )) linear systems of equations and corresponding CG (Niter CG ) and PCG (Niter P CG ) iteration counts as a function of the mesh-size h for the SIPG and LDG methods for p = 6. K(A) SIPG (α = 10, β = 0) LDG (α = 10, β = 1) Niter CG K(T DG ) Niter P CG K(A) Niter CG K(T DG ) Niter P CG h = 1/ h = 1/ h = 1/ The last numerical experiment of this section aims at verifying the uniformity of the proposed preconditioner with respect to the penalization coefficient α. In this case, we consider the same test case presented above, but we now fix the polynomial approximation degree p = 2 and increase α. The numerical data obtained are reported in Table 4: as done before, we compare the spectral condition numbers of the unpreconditioned and preconditioned systems and the iteration counts of the CG and PCG methods. As predicted from theory, while K(A) grows like α, the values of K(T DG ) are constant. Table 4: Condition number of the unpreconditioned (K(A)) and preconditioned (K(T DG )) linear systems of equations and corresponding CG (Niter CG ) and PCG ) iteration counts as a function of α for the SIPG and LDG methods (N P CG iter SIPG (p = 2, β = 0) LDG (p = 2, β = 1) K(A) Niter CG K(T DG ) Niter P CG K(A) Niter CG K(T DG ) Niter P CG α = α = α = α = α = α =
22 22 Paola F. Antonietti et al. 6.2 Example 2 We now address a more general model problem by introducing a diffusion tensor ρ in (1). We then rewrite the weak formulation (1) as follows: find u V, such that ρ u v dx = fv dx v V, Ω where, for simplicity we assume that [ ] ρ11 (x) 0 ρ ρ(x) =, 0 ρ 22 (x) Ω with ρ 11 (x) and ρ 22 (x) (possibly discontinuous) positive constants a.e. in Ω. We also assume that there exists an initial partition T h0 such that the (possible) discontinuities of ρ 11 (x) and ρ 22 (x) are aligned with T h0. We then modify the bilinear form (3) as follows, restricting, for simplicity, ourselves to the SIPG bilinear form, i.e., θ = 0 and β = 0, A ρ (u, v) := ρ u v dx+ ( u R( {ρ } A v )+R( {ρ } A u ) v) dx κ T h κ κ T h κ + F F h F σ {ρ } A u v ds where, on each face F F h, {ρ } A is defined as 2 (nt F,κ ρ + n + F,κ +)(n T F,κ ρ n F,κ ) {ρ } A := n T F,κ ρ + n + F,κ + + n T if F F F,κ ρ n h, I F,κ n T F ρn F if F Fh B, see, e.g., [25]. We then consider the following DG formulation: A ρ (u, v) = fv dx v V hp. (37) Ω The numerical results reported in this section have been obtained choosing piecewise quadratic elements, i.e., p = 2 and the penalization parameter α = 10. For the first test we choose an elementwise discontinuous isotropic diffusion coefficient, i.e. [ ] ρ11 (x) 0 ρ(x) =, 0 ρ 11 (x) where ρ 11 (x) follows an elementwise red-black partitioning of the computational domain, which means ρ 11 = 1 in red regions and ρ 11 = 10 γ, γ = 0, 1,..., 7, in the black ones. In Figure 3, we numerically evaluate the constant appearing in the inequality A ρ (Q h v, Q h v) C Qh A ρ (v, v), (38)
23 A uniform additive Schwarz preconditioner for high order DG methods 23 which is the analogous of (13), and observe that the presence of a discontinuous diffusion coefficient introduces a linear dependence of C Qh on the ratio between the maximum and minimum value of ρ. In Table 5 we report the estimated values of the constants C J 1 and C J 2 appearing in Lemma 7, and C C 1 and C C 2 appearing in Lemma 9, respectively, as a function of ρ 11. Here, we note that the two preconditioners, built on the corresponding subspaces, are robust with respect to the contrast in the diffusion tensor. From the results obtained by these numerical tests we can conclude that the condition number of the preconditioned operator T DG depends linearly on the jumps of the diffusion coefficient, as further confirmed by the numerical tests reported in Table 6 where the spectral condition number of T DG as well as the PCG iteration counts as a function of ρ 11 are shown. For the sake of comparison, the same quantities are also shown for the unpreconditioned system, showing that, as expected also the unpreconditioned operator has a spectral condition number that depends linearly on the contrast of the diffusion coefficient C Qh ρ ρ 11 Fig. 3: Numerical evaluation of the constant C Qh in (38) as a function of ρ 11. SIPG method (37) (p = 2, α = 10). Table 5: Left: numerical evaluation of the constants C J 1 and C J 2 of Lemma 7 as a function of ρ 11 for the SIPG method (37) (p = 2, α = 10); right: numerical evaluation of the constants C C 1 and C C 2 of Lemma 9 as a function of ρ 11. C J 1 C J 2 C C 1 C C 2 ρ 11 = ρ 11 = ρ 11 = ρ 11 = ρ 11 = ρ 11 = ρ 11 = ρ 11 =
24 24 Paola F. Antonietti et al. Table 6: Condition number of the unpreconditioned (K(A)) and preconditioned (K(T DG )) linear systems of equations and corresponding CG (Niter CG ) and PCG (Niter P CG) iteration counts as a function of ρ 11 for the SIPG method (37) (p = 2, α = 10). K(A) Niter CG K(T DG ) Niter P CG ρ 11 = ρ 11 = ρ 11 = ρ 11 = ρ 11 = ρ 11 = ρ 11 = ρ 11 = Next, we investigate the behavior of the proposed technique whenever the diffusion tensor features a strong anisotropy. For this set of experiments we set [ ] 1 0 ρ(x) =, 0 ρ 22 where ρ 22 is a positive constant over the whole computational domain. We have repeated the same set of experiments as before, again with p = 2. In Figure 4 we report the estimated constant C Qh in (38) as a function of ρ 22 = 10 γ, γ = 0, 1,..., 7. We observe that C Qh seems to be independent of the anisotropy of ρ(x). The computed constants C J 1 and C J 2 of Lemma 7 and C C 1 and C C 2 of C Qh ρ 22 Fig. 4: Numerical evaluation of the constant C Qh in (38) as a function of ρ 22. SIPG method (37) (p = 2, α = 10). Lemma 9 as a function of ρ 22 are reported in Table 7. We note that C J 2 C C 2 seem to be asymptotically independent of ρ 22. and
25 A uniform additive Schwarz preconditioner for high order DG methods 25 Table 7: Left: numerical evaluation of the constants C J 1 and C J 2 of Lemma 7 as a function of ρ 22 for the SIPG method in (37) (p = 2); right: numerical evaluation of the constants C C 1 and C C 2 of Lemma 9 as a function of ρ 22. C J 1 C J 2 C C 1 C C 2 ρ 22 = ρ 22 = ρ 22 = ρ 22 = ρ 22 = ρ 22 = ρ 22 = ρ 22 = The numerical data shown in Figure 4 and Table 7 again suggest the expected behavior of the final preconditioner. Indeed the condition number K(T DG ) seems to be asymptotically constant as ρ 22 grows, as confirmed by the results reported in Table 8. Notice however that also the unpreconditioned DG discretization seems to be itself asymptotically robust with respect to the anisotropy. Table 8: Condition number of the unpreconditioned (K(A)) and preconditioned (K(T DG )) linear systems of equations and corresponding CG (Niter CG ) and PCG ) iteration counts as a function of ρ 22 for the SIPG method (N P CG iter K(A) Niter CG K(T DG ) Niter P CG ρ 22 = ρ 22 = ρ 22 = ρ 22 = ρ 22 = ρ 22 = ρ 22 = ρ 22 = Acknowledgements Part of this work was developed during the visit of the second author at the Pennsylvania State University. Special thanks go to the Center for Computational Mathematics and Applications (CCMA) at the Mathematics Department, Penn State for the hospitality and support. The work of the fourth author was supported in part by NSF DMS , and NSF DMS The work of the first author was partially supported by SIR Project n. RBSI14VT0S PolyPDEs: Non-conforming polyhedral finite element methods for the approximation of partial differential equations funded by MIUR.
26 26 Paola F. Antonietti et al. References 1. P. F. Antonietti and B. Ayuso. Schwarz domain decomposition preconditioners for discontinuous Galerkin approximations of elliptic problems: non-overlapping case. M2AN Math. Model. Numer. Anal., 41(1):21 54, P. F. Antonietti and B. Ayuso. Multiplicative Schwarz methods for discontinuous Galerkin approximations of elliptic problems. M2AN Math. Model. Numer. Anal., 42(3): , P. F. Antonietti and B. Ayuso. Two-level Schwarz preconditioners for super penalty discontinuous Galerkin methods. Commun. Comput. Phys., 5(2-4): , P. F. Antonietti, B. Ayuso, S. Bertoluzza, and M. Pennacchio. Substructuring preconditioners for an hp domain decomposition method with Interior Penalty mortaring. Calcolo, 52(3): , P. F. Antonietti, B. Ayuso, S. C. Brenner, and L.-Y. Sung. Schwarz methods for a preconditioned WOPSIP method for elliptic problems. Comput. Meth. in Appl. Math., 12(3): , P. F. Antonietti and P. Houston. A class of domain decomposition preconditioners for hp-discontinuous Galerkin finite element methods. J. Sci. Comput., 46(1): , P. F. Antonietti, P. Houston, M. Sarti, and M. Verani. Multigrid algorithms for hpversion Interior Penalty Discontinuous Galerkin methods on polygonal and polyhedral meshes. MOX Report 55/2014, Submitted, P. F. Antonietti, P. Houston, and I. Smears. A note on optimal spectral bounds for nonoverlapping domain decomposition preconditioners for hp version discontinuous Galerkin methods. Int. J. Numer. Anal. Model., 4(13): , P. F. Antonietti, M. Sarti, and M. Verani. Multigrid algorithms for hp-discontinuous Galerkin discretizations of elliptic problems. SIAM J. Numer. Anal., 53(1): , D. N. Arnold. An interior penalty finite element method with discontinuous elements. SIAM J. Numer. Anal., 19(4): , D. N. Arnold, F. Brezzi, B. Cockburn, and L. D. Marini. Unified analysis of discontinuous Galerkin methods for elliptic problems. SIAM J. Numer. Anal., 39(5): , 2001/ A. T. Barker, S. C. Brenner, and L.-Y. Sung. Overlapping Schwarz domain decomposition preconditioners for the local discontinuous Galerkin method for elliptic problems. J. Numer. Math., 19(3): , S. C. Brenner. Poincaré-Friedrichs inequalities for piecewise H 1 functions. SIAM J. Numer. Anal., 41(1): , K. Brix, M. Campos Pinto, C. Canuto, and W. Dahmen. Multilevel Preconditioning of Discontinuous Galerkin Spectral Element Methods Part I: Geometrically Conforming Meshes. IMA J. Numer. Anal., To appear. 15. K. Brix, M. Pinto, and W. Dahmen. A multilevel preconditioner for the interior penalty discontinuous Galerkin method. SIAM J. Numer. Anal., 46(5): , cited By (since 1996) E. Burman. A unified analysis for conforming and nonconforming stabilized finite element methods using interior penalty. SIAM J. Numer. Anal., 43(5): (electronic), E. Burman and A. Ern. Continuous interior penalty hp-finite element methods for advection and advection-diffusion equations. Math. Comp., 76(259): (electronic), C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang. Spectral methods. Scientific Computation. Springer-Verlag, Berlin, Fundamentals in single domains. 19. C. Canuto, L. F. Pavarino, and A. B. Pieri. BDDC preconditioners for continuous and discontinuous Galerkin methods using spectral/hp elements with variable local polynomial degree. IMA J. Numer. Anal., 34(3): , B. Cockburn and C.-W. Shu. The local discontinuous Galerkin method for timedependent convection-diffusion systems. SIAM J. Numer. Anal., 35(6): (electronic), 1998.
ASM-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 informationarxiv: v2 [math.na] 29 Nov 2016
Noname manuscript No. (will be inserted by the editor) Multigrid algorithms for hp-version Interior Penalty Discontinuous Galerkin methods on polygonal and polyhedral meshes P. F. Antonietti P. Houston
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 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 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 informationDiscontinuous Galerkin Methods: Theory, Computation and Applications
Discontinuous Galerkin Methods: Theory, Computation and Applications Paola. Antonietti MOX, Dipartimento di Matematica Politecnico di Milano.MO. X MODELLISTICA E CALCOLO SCIENTIICO. MODELING AND SCIENTIIC
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 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 informationOptimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36
Optimal multilevel preconditioning of strongly anisotropic problems. Part II: non-conforming FEM. Svetozar Margenov margenov@parallel.bas.bg Institute for Parallel Processing, Bulgarian Academy of Sciences,
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 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 informationGeometric Multigrid Methods
Geometric Multigrid Methods Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University IMA Tutorial: Fast Solution Techniques November 28, 2010 Ideas
More 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 informationA Subspace Correction Method for Discontinuous Galerkin Discretizations of Linear Elasticity Equations
www.oeaw.ac.at A Subspace Correction Method for Discontinuous Galerkin Discretizations of Linear lasticity quations B. Ayuso de Dios, I. Georgiev, J. Kraus, L. Zikatanov RICAM-Report 2011-29 www.ricam.oeaw.ac.at
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationGALERKIN DISCRETIZATIONS OF ELLIPTIC PROBLEMS
SIAM J. NUMER. ANAL. Vol. 53, No. 1, pp. 598 618 c 2015 Society for Industrial and Applied Mathematics MULTIGRID ALGORITHMS FOR hp-discontinuous GALERKIN DISCRETIZATIONS OF ELLIPTIC PROBLEMS PAOLA F. ANTONIETTI,
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 information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More informationMULTIGRID 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 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 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 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 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 informationThe Discontinuous Galerkin Finite Element Method
The Discontinuous Galerkin Finite Element Method Michael A. Saum msaum@math.utk.edu Department of Mathematics University of Tennessee, Knoxville The Discontinuous Galerkin Finite Element Method p.1/41
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 informationA NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION
A NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION JOHNNY GUZMÁN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS Abstract. The analysis of finite-element-like Galerkin discretization techniques for the
More informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More informationDiscontinuous Galerkin Method for interface problem of coupling different order elliptic equations
Discontinuous Galerkin Method for interface problem of coupling different order elliptic equations Igor Mozolevski, Endre Süli Federal University of Santa Catarina, Brazil Oxford University Computing Laboratory,
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 informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
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 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] 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 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 informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
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 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 informationA note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations
A note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations Bernardo Cockburn Guido anschat Dominik Schötzau June 1, 2007 Journal of Scientific Computing, Vol. 31, 2007, pp.
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 informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More 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 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 informationPRECONDITIONING OF DISCONTINUOUS GALERKIN METHODS FOR SECOND ORDER ELLIPTIC PROBLEMS. A Dissertation VESELIN ASENOV DOBREV
PRECONDITIONING OF DISCONTINUOUS GALERKIN METHODS FOR SECOND ORDER ELLIPTIC PROBLEMS A Dissertation by VESELIN ASENOV DOBREV Submitted to the Office of Graduate Studies of Texas A&M University in partial
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 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 informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More 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 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 informationRemarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable?
Remarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable? Thomas Apel, Hans-G. Roos 22.7.2008 Abstract In the first part of the paper we discuss minimal
More informationPREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO
PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS
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 informationThe Virtual Element Method: an introduction with focus on fluid flows
The Virtual Element Method: an introduction with focus on fluid flows L. Beirão da Veiga Dipartimento di Matematica e Applicazioni Università di Milano-Bicocca VEM people (or group representatives): P.
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 information31. Successive Subspace Correction method for Singular System of Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 2003 DDM.org 31. Successive Subspace Correction method for Singular
More informationHybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
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 informationHigh order, finite volume method, flux conservation, finite element method
FLUX-CONSERVING FINITE ELEMENT METHODS SHANGYOU ZHANG, ZHIMIN ZHANG, AND QINGSONG ZOU Abstract. We analyze the flux conservation property of the finite element method. It is shown that the finite element
More informationEnergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations of the Navier-Stokes equations
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 19007; 1:1 [Version: 00/09/18 v1.01] nergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationA Two-grid Method for Coupled Free Flow with Porous Media Flow
A Two-grid Method for Coupled Free Flow with Porous Media Flow Prince Chidyagwai a and Béatrice Rivière a, a Department of Computational and Applied Mathematics, Rice University, 600 Main Street, Houston,
More informationComparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems
Comparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems Nikola Kosturski and Svetozar Margenov Institute for Parallel Processing, Bulgarian Academy of Sciences Abstract.
More informationNew Multigrid Solver Advances in TOPS
New Multigrid Solver Advances in TOPS R D Falgout 1, J Brannick 2, M Brezina 2, T Manteuffel 2 and S McCormick 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, P.O.
More informationNumerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods
Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods Contents Ralf Hartmann Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Lilienthalplatz 7, 3808
More informationKey words. preconditioned conjugate gradient method, saddle point problems, optimal control of PDEs, control and state constraints, multigrid method
PRECONDITIONED CONJUGATE GRADIENT METHOD FOR OPTIMAL CONTROL PROBLEMS WITH CONTROL AND STATE CONSTRAINTS ROLAND HERZOG AND EKKEHARD SACHS Abstract. Optimality systems and their linearizations arising in
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
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 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 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 informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationA Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators
A Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators Jeff Ovall University of Kentucky Mathematics www.math.uky.edu/ jovall jovall@ms.uky.edu Kentucky Applied and
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 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 informationAn a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element
Calcolo manuscript No. (will be inserted by the editor) An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element Dietrich Braess Faculty of Mathematics, Ruhr-University
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 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 informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
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 informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
More informationMultigrid and Domain Decomposition Methods for Electrostatics Problems
Multigrid and Domain Decomposition Methods for Electrostatics Problems Michael Holst and Faisal Saied Abstract. We consider multigrid and domain decomposition methods for the numerical solution of electrostatics
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 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 informationA local-structure-preserving local discontinuous Galerkin method for the Laplace equation
A local-structure-preserving local discontinuous Galerkin method for the Laplace equation Fengyan Li 1 and Chi-Wang Shu 2 Abstract In this paper, we present a local-structure-preserving local discontinuous
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 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 informationPaola F. Antonietti 1, Luca Formaggia 2, Anna Scotti 3, Marco Verani 4 and Nicola Verzotti 5. Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher MIMETIC FINITE DIFFERENCE APPROXIMATION OF FLOWS IN FRACTURED POROUS MEDIA Paola
More informationAdaptive algebraic multigrid methods in lattice computations
Adaptive algebraic multigrid methods in lattice computations Karsten Kahl Bergische Universität Wuppertal January 8, 2009 Acknowledgements Matthias Bolten, University of Wuppertal Achi Brandt, Weizmann
More informationA DISCONTINUOUS GALERKIN METHOD FOR ELLIPTIC INTERFACE PROBLEMS WITH APPLICATION TO ELECTROPORATION
A DISCONTINUOUS GALERIN METHOD FOR ELLIPTIC INTERFACE PROBLEMS WITH APPLICATION TO ELECTROPORATION GRÉGORY GUYOMARC H AND CHANG-OC LEE Abstract. We present a discontinuous Galerkin (DG) method to solve
More informationGoal. Robust A Posteriori Error Estimates for Stabilized Finite Element Discretizations of Non-Stationary Convection-Diffusion Problems.
Robust A Posteriori Error Estimates for Stabilized Finite Element s of Non-Stationary Convection-Diffusion Problems L. Tobiska and R. Verfürth Universität Magdeburg Ruhr-Universität Bochum www.ruhr-uni-bochum.de/num
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 informationLocal flux mimetic finite difference methods
Local flux mimetic finite difference methods Konstantin Lipnikov Mikhail Shashkov Ivan Yotov November 4, 2005 Abstract We develop a local flux mimetic finite difference method for second order elliptic
More 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 informationDomain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions
Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions Ernst P. Stephan 1, Matthias Maischak 2, and Thanh Tran 3 1 Institut für Angewandte Mathematik, Leibniz
More informationA SHORT NOTE COMPARING MULTIGRID AND DOMAIN DECOMPOSITION FOR PROTEIN MODELING EQUATIONS
A SHORT NOTE COMPARING MULTIGRID AND DOMAIN DECOMPOSITION FOR PROTEIN MODELING EQUATIONS MICHAEL HOLST AND FAISAL SAIED Abstract. We consider multigrid and domain decomposition methods for the numerical
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 informationSchwarz Preconditioner for the Stochastic Finite Element Method
Schwarz Preconditioner for the Stochastic Finite Element Method Waad Subber 1 and Sébastien Loisel 2 Preprint submitted to DD22 conference 1 Introduction The intrusive polynomial chaos approach for uncertainty
More informationA Posteriori Error Estimates For Discontinuous Galerkin Methods Using Non-polynomial Basis Functions
Lin Lin A Posteriori DG using Non-Polynomial Basis 1 A Posteriori Error Estimates For Discontinuous Galerkin Methods Using Non-polynomial Basis Functions Lin Lin Department of Mathematics, UC Berkeley;
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 informationBUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM ERIK BURMAN AND BENJAMIN STAMM Abstract. We propose a low order discontinuous Galerkin method for incompressible flows. Stability of the
More information