Some Domain Decomposition Methods for Discontinuous Coefficients
|
|
- Katrina Briggs
- 5 years ago
- Views:
Transcription
1 Some Domain Decomposition Methods for Discontinuous Coefficients Marcus Sarkis WPI RICAM-Linz, October 31, 2011 Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
2 Outline Discretizations P1 conforming, RT0, P1 nonconforming, DG, Mortar Some important differences between them Overlapping Schwarz Methods and/with Substructuring Methods COARSE SPACES Overlapping AS (additive Schwarz), ASHO (AS with harmonic overlap), RASHO (Restricted ASHO), OBBD (Overlapping BDD) Iterative Substructuring, BDDC First we consider discontinuous coefficients however constant inside substructures. Then, we consider more general discontinuities A general reference for DD: Andrea Toselli and Olof Widlund, Domain Decomposition Methods - Algorithms and Theory. Springer Series in Computational Mathematics, Vol. 34, Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
3 Geometrically Nonconforming Subdomain Partition Ω = N i=1 Ω i Ω i disjoint shaped regular polygonal subdomains of diameter O(H i ) T hi (Ω i ) shape regular triangulations In this talk: geometrically conforming partitions and matching meshes The Ω i are simplices (to drop later) with constant coefficient ρ i Subdomains, faces, edges are open sets Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
4 Problem of Interest Continuous PDE. Find u H0 1 (Ω) such that a(u, v) := N (ρ i u, v)dx = Ω i i=1 Ω fvdx =: f (v) v H 1 0 (Ω) Conforming FEM. V h (Ω) H0 1 (Ω) : continuous piecewise linear functions in Ω. Find u V h (Ω) such that a(u, v) = f (v) v V h (Ω) Nonconforming FEM. ˆV h (Ω): piecewise linear functions and continuous at the middle points of the edges of T h (Ω) and zero at the middle points of the edges of T h ( Ω). Find u ˆV h (Ω) such that â h (u, v) := (ρ u τ T h (Ω) τ Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
5 Conforming FEM P 1 conforming: V H (Ω) H 1 0 (Ω) Weighted L 2 approximation and weighted H 1 stability u I H u 2 L 2 ρ (Ω) H2 u 2 H 1 ρ (Ω) I H u 2 H 1 ρ (Ω) u 2 H 1 ρ (Ω) First inequality not always possible with constant independently of ρi u 2 L 2 ρ (Ω) := (ρu, u)dx u 2 Hρ 1 (Ω) := (ρ u, u)dx Ω Ω Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
6 Nonconforming FEM (both inequalities hold) Step 1: take averages where the coefficient is constant: average on elements, boundary of elements, faces (edges) in 3D (2D) Step 2: obtain Poincaré-Friedrich type inequalities inside each Ω i Step 3: establish local H 1 stability and L 2 approximation in Ω i Step 4: establish global weighted H 1 stab. and L 2 appr. in Ω Crouzeix-Raviart P1 nonconforming: VH (Ω) H 1 0 (Ω) Broken weighted H1 stability and weighted L 2 approximation ÎHu Ĥ1 ρ,h (Ω) u H 1 ρ (Ω) u ÎHu L 2 ρ (Ω) H u H 1 ρ (Ω) ÎH u averages u on faces (edges) of Ω i u 2 := N Ĥρ,H 1 (Ω) i=1 Ω i (ρ i u, u)dx Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
7 Nonconforming FEM (2D Local Analysis) m 1, m 2, m 3 middle points of edges E 1, E 2, E 3 of Ω i ū m1, ū m2, ū m3 edge averages ÎHu 2 H 1 (Ω i ) ū m 2 ū m1 2 + ū m3 ū m1 2 Note that ū m2 ū m1 2 E = 2 (u ū m1 )ds E 2 1ds 2 H 1 u ū m1 2 L 2 (E 2 ) H 2 u ū m1 2 L 2 (Ω i ) + u ū m 1 2 H 1 (Ω i ) By Friedrich s inequality, the H 1 stability in Ω i holds u ÎHu has zero average on edges of Ω i so Friedrich holds on Ω i u ÎHu 2 L 2 (Ω i ) H2 u ÎHu 2 H 1 (Ω i ) 2H2 ( u 2 H 1 (Ω i ) + ÎHu 2 H 1 (Ω i ) ) Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
8 Remark on CR elements CR system is spectrally equivalent to the reduced hybridizable RT0 RT0 has mass conservation properties For discontinuous coefficients, RT0 gives better results for the velocity and pressure (after postprocessing) compared to conforming P 1 In certain cases, discrete Poincaré-Friedrich-type inequalities with constants independently of coefficients are more difficult to obtain than for conforming elements Let V be a coarse vertex of T H (Ω). Let Ω V be the union of all Ω j touching V (let us assume there are 6 Ω j ) Coefficients (M, 1, 1, M, 1, 1) (anti-clockwise orientation). M very large Conforming P1 : u u(v ) L 2 ρ (Ω V ) H u H 1 ρ (Ω V ) for u V H (Ω V ) CR: broken weighted Poincaré does not hold in ÎH(Ω V ) Quasi-monotone coefficients: holds for ÎH(Ω V ) and H 1 (Ω V ) The preconditioning analysis for CR more complicated than conforming Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
9 CR: Boundary Element Averages How to define the values at the subdomain edges? Let us take an edge E ij := Ω i Ω j Let ū i, ū j element averages on Ω i and Ω j Define the ĨHu by ū ij,β = ρ β i ρβ i + ρ β j ū i + ρ β i ρβ j + ρ β j ū j, β [ 1 2, ] Denote m ij, m ik, m il middle points of E ij, E ik, E il Let ū mij, ū mik, ū mil be edge averages of u on E ij, E ik, E il Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
10 CR: Element Averages (cont.) Broken weighted H 1 stability ρ i ū ij,β ū ik,β 2 3ρ i ( ūij,β ū mij 2 + ū mij ū mik 2 + ū mik ū ik,β 2) ρ i ū ij,β ū mij 2 = ρ i ( ρ i ( ρ β i Weighted L 2 stability ρ β i ρ β j ρ β i + ρ β j + ρ β (ū i ū mij ) + ρβ j j ρ β i + ρ β (ū j ū mij ) j ) 2 ρ j if β 1/2 ) 2 ρ i u ĨHu 2 L 2 (Ω i ) ρ ih 2 u ĨHu 2 H 1 (Ω i ) + ρ ih u ĨHu 2 L 2 (E ij ) u ĨHu = ρ β i ρβ i + ρ β j (u ū i ) + ρβ j ρ β i + ρ β j (u ū j ) Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
11 CR: Face Coarse Spaces V H and ṼH V h Continuity at the middle points of the coarse edges ( V H ) does not imply continuity at the middle points of the fine edges ( V h ) Edge coarse space (in 2D) and Face coarse space (in 3D) Step 1: On each subdomain edge E ij : make Î E H u equal to (ÎHu)(m ij ) at every middle point of the fine triangulation on E ij Step 2: Discrete harmonic extension to define Î E H u inside Ω i The resulting coarse space can be constructed easily for general Ω i Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
12 CR-Face Coarse Basis Functions E ij := Ω i Ω j Face coarse basis functions θ Eij. On fine nodes, let: θ Eij : 1 on E ij θeij := 0 on all subdomain edges but E ij θeij : discrete harmonic extension inside Ω i and Ω j θeij := 0 elsewhere Note that the support of θ Eij is Ω i Ω j What is Î E H? Î E H u = Edges E T H (Ω) ū E θ E ˆV E H is the range of Î E H Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
13 CR: NN Coarse Basis Functions For each subdomain Ω i define θ i = E ij Ω i ρ β i β ρi + ρ β j θ Eij Support of θ i is Ω ext i Î NN H given by : union of all Ω j sharing an edge with Ω i Î NN H N = ū i θ i i=1 Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
14 CR: Face Coarse Space (cont.) Note that ÎHu and Î E H u have the same average on each coarse edge E ij It is possible to show: Weighted L 2 approximation u Î E H u 2 L 2 ρ (Ω i ) H2 u 2 Ĥ 1 ρ,h (Ω i ) Broken weighted H1 stability ( Î H E u log H ) u 2 Ĥρ,h 1 (Ω i ) h Ĥρ,h 1 (Ω i ) Broken weighted H1 norm u 2 Ĥ 1 ρ,h (Ω) := N i=1 τ T h (Ω i ) τ (ρ i u, u)dx ÎH NN : the bounds are similar, however, u 2 Ĥ 1 ρ,h (Ω i ) replaced by u 2 Ĥ 1 ρ,h (Ωext i ) Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
15 References (Two-level) Two-level: ˆV h = ˆV E(NN) H + N i=1 ˆV δ h (Ωδ i ) M. Sarkis. Nonstandard coarse spaces and Schwarz methods for elliptic problems with discontinuous coefficients using non-conforming elements. Numer. Math., 77(3), 1997, pp An earlier version in: Two-level Schwarz methods for P 1 nonconforming finite elements and discontinuous coefficients. Proceedings of the Sixth Copper Mountain Conference on Multigrid Methods, Number 3224, Part 2, pages , Hampton VA, NASA Important results on this paper: The β-weighting The face coarse spaces Iterative substructuring coarse space and overlapping Schwarz Tools to analyze preconditioners for nonconforming FEM: stable maps between conforming and nonconforming spaces Order of (1 + H/δ) (1 + log H/h) bounds independent of coefficients Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
16 References (Multi-level) Decomposition: ˆV hl = ˆV E(NN) h L + ˆV hl + L l=0 V h l (conforming spaces) E(NN) Bilinear forms: â h on ˆV h L, diag. of â h on ˆV hl, a-bpx on L l=0 V h l Order of (1 + log 2 H/h) bounds independent of coefficients M. Sarkis, Multilevel methods for P 1 nonconforming finite elements and discontinuous coefficients in three dimensions. Decomposition Methods in Scientific and Engineering (Univesity Park, PA, 1993), Contemp. Math., Vol. 180, Amer. Math. Soc., Providence, RI, 1994, pp M. Sarkis, Schwarz preconditioners for elliptic problems with discontinuous coefficients using conforming and non-conforming elements. PhD Dissertation, Technical Report 671, Courant Institute, September Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
17 P 1 Conforming FEM and Discrete Sobolev Inequalities V H (Ω) H0 1(Ω): P 1 conforming FEM in the coarse triangulation V h (Ω) V H (Ω): P 1 conforming FEM in the fine triangulation Pointwise interpolation I H : V h (Ω) V H (Ω), defined by (I H u)(x) = u(x), x vertices of the Ω i Discrete Sobolev inequalities u 2 L (Ω i ) (1 + log H h ) u 2 H 1 (Ω i ) in 2D u 2 L (Ω i ) H h u 2 H 1 (Ω i ) in 3D Same estimates hold for u u(x) 2 L (Ω i )... u 2 H 1 (Ω i ) x Ω i Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
18 Pointwise Interpolation: L 2 and H 1 Poincaré Inequalities L 2 approximation for 2D. For u V h (Ω i ) u I H u 2 L 2 (Ω i ) (1 + log H h )H2 u 2 H 1 (Ω i ), u V h(ω i ) H 1 stability for 2D I H u 2 H 1 (Ω i ) (1 + log H h ) u 2 H 1 (Ω i ), u V h(ω i ) In 3D, replace (1 + log H h ) by H h Estimates done locally in Ω i, hence, holds for weighted norms in Ω H/h factor in 3D (not so good) Not trivial to extend to general Ω i and discontinuous coefficients Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
19 Discrete Sobolev Inequalities for Conforming FEM (3D) Let u V be the value of u V h (Ω i ) at a vertex V of Ω i u u V 2 L 2 (Ω i ) (H h )H2 u 2 H 1 (Ω i ) Let ū E be the edge average of u V h (Ω i ) on an edge E of Ω i u ū E 2 L 2 (Ω i ) (1 + log H h )H2 u 2 H 1 (Ω i ) Let ū F be the face average of u V h (Ω i ) on a face F of Ω i u ū F 2 L 2 (Ω i ) H2 u 2 H 1 (Ω i ) Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
20 Coarse Basis Functions (3D) Similar construction as we did for the nonconforming case I B H u = V u V θ V + E ū E θ E + F ū F θ F θv := 1 at vertex V and zero at remaining interface nodes θ E := 1 at nodes on edge E and zero at remaining interface nodes θ F := 1 at nodes on face F and zero at remaining interface nodes θ V, θ E and θ F : discrete harmonic inside the subdomains H 1 norm of the basis functions θ V 2 H 1 (Ω i ) h θ E 2 H 1 (Ω i ) H θ F 2 H 1 (Ω i ) (1 + log H h )H Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
21 Stability and Approximation The θ V, θ E and θ F form a partition of unity on Ω i I B H u = V (ū u V )θ V + E (ū ū E )θ E + F (ū ū F )θ F Let us bound for instance (ū ū F )θ F H 1 (Ω i ), the others are similar (ū ū F )θ F 2 H 1 (Ω i ) ū ū F 2 (1 + log H h )H ū ū F 2 u ū 2 + u ū F 2 1 ) ( u H 3 ū 2 L 2 (Ω i ) + u ū F 2 L 2 (Ω i ) H 1 stability I B H u 2 H 1 (Ω i ) (1 + log H h ) u 2 H 1 (Ω i ) L 2 approximation u I B H u 2 L 2 (Ω i ) H2 u I B H u 2 H 1 (Ω i ) Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
22 Subspace Decomposition Global space V. Example: V = V h (Ω) H 1 0 (Ω) Global bilinear form a(u, v) for u, v V. Example: a(u, v) = (ρ u, v)dx Local spaces V i, i = 1 : N. Example: V i = V h (Ω) H0 1 (Ω δ i ) Ω Local bilinear forms a i (u i, v i ) for u i, v i V i. Example: a i (u i, v i ) = (ρ u i, v i )dx Ω δ i Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
23 Subspace Decomposition (cont) Coarse space V 0. Example: V 0 = V B h or V H h or V NN h Coarse bilinear form a0 (u 0, v 0 ) for u 0, v 0 V 0. Example: a 0 (u 0, v 0 ) = a(u 0, v 0 ) Extension operators R T i : V i V. Examples: Local spaces: u = R T i u i is the zero extension H 1 0 (Ωδ i ) to H1 0 (Ω) Coarse space: u = R T 0 u 0 is the identity since V 0 V Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
24 Abstract Theory of Schwarz Methods Additive Schwarz T as = ( N i=0 RT i A 1 i R i )A = (B)A Lower bound: C0 2 a(u, u) a(t asu, u) Lions Lemma: Find a C 0 > 0 such that, for any u V, there exists a decomposition u = N i=0 RT i u i such that N a i (u i, u i ) C0 2 a(u, u) i=0 Upper bound: a(t as u, u) ωρ(e)a(u, u) Inexact solvers: Find an ω where for any i = 0 : N and ui V i a(r T i u i, R T i u i ) ωa i (u i, u i ) Strengthened Cauchy-Schwarz: Find a upper bound for the spectral radius of E = {ɛ ij } i,j=0:n a(r T i u i, R T j u j ) ɛ ij a(r T i u i, R T i u i ) 1/2 a(r T j u j, R T j u j ) 1/2 Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
25 Decomposition Partition of Unity φ i (x) V h (Ω) for the overlapping subdomains Ω δ i supp(φ i ) Ω δ i 0 φ i (x) 1, x Ω δ i N φ i (x) = 1, i=1 φ i C/δ, Decomposition of u = u 0 + N i=1 u i x Ω 1 i N u 0 = I B h u. Let w = u u 0 ui = I h (φ i w) Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
26 Lower bound estimation We need to estimate a(u 0, u 0 ) + It follows from Hρ(Ω) 1 stability N a i (u i, u i ) i=1 a(u 0, u 0 ) = Ih B u 2 Hρ 1(Ω) (1 + log H h ) u 2 Hρ 1(Ω) = (1 + log H )a(u, u) h Easy to see: I h (φ i w) 2 w φ i 2 + φ i w 2 (ρ I h (φ i w), I h (φ i w))dx 1 δ 2 ρw 2 + ρ w 2 dx Ω δ i Ω δ i Ω δ i (1 + log H h )(H δ )2 u 2 Hρ(Ω 1 ext i ) + (1 + log H h ) u 2 Hρ(Ω 1 ext i ) Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
27 Small Overlapping Technique To improve the factor (H/δ) 2 to H/δ Note that the support of φ i is a δ-layer near Ω i Divide the support of φ i in pieces: the layer inside Ω i, and the layers Ω δ i Ω j. Note that in each of these layers, ρ is constant We can estimate the L 2 norm in each layer by the L 2 and H 1 seminorms on the corresponding subdomain the layer belongs 1D: Estimate w 2 L 2 (0,δ) in terms of w 2 L 2 (0,H) and w 2 H 1 (0,H) Use Fundamental Theorem of Calculus and Cauchy-Schwarz: w 2 (x) w 2 (y) + H w 2 H 1 (0,H) x, y (0, H) Integrate x in (0, δ) and y in (0, H) H w 2 L 2 (0,δ) δ w 2 L 2 (0,H) + H2 δ w 2 H 1 (0,H) Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
28 References Two-Level: C. R. Dohrmann, A. Klawonn, and O. B. Widlund. Domain Decomposition for Less Regular Subdomains: Overlapping Schwarz in Two Dimensions, March 2007, SIAM J. Numer. Anal., Vol. 46(4), 2008, Multi-level: M. Dryja, M. Sarkis and O. Widlund. Multilevel Schwarz methods for elliptic problems with discontinuous coefficients in three dimensions. Numer. Math., 72(3), 1996, pp M. Dryja and O. B. Widlund. Domain Decomposition Algorithms with Small Overlap, SIAM J. Sci. Stat. Comput., Vol. 15, No. 3, May 1994, M. Dryja, B. F. Smith, and O. B. Widlund. Schwarz Analysis of Iterative Substructuring Algorithms for Elliptic Problems in Three Dimensions, SIAM J. Numer. Anal., Vol. 31, No. 6, December 1994, Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
29 Primal Iterative Substructuring Methods Eliminate all the interior variables inside of the Ω i The idea: replace the local overlapping Dirichlet solvers by nonoverlapping Neumann solvers Local spaces V i := V h (Ω i ) with vertex, edge, face, or/and the whole boundary average constraint Local bilinear forms a i (u i, v i ) = (ρ i u i, v i )dx Ω i Extension operator R T i : u i u Nodes x Ωi define u(x) = u i (x)(ρ β i )/( j ρβ j ) Nodes x Ωj \ Ω i define u(x) = 0 Extend discrete harmonically inside the subdomains Note that this algorithm requires Neumann matrices, so, not algebraic Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
30 References M. Dryja and O. B. Widlund. Schwarz Methods of Neumann-Neumann Type for Three- Dimensional Elliptic Finite Element Problems, Comm. Pure Appl. Math., Vol. 48, No. 2, February 1995, M. Dryja, B. F. Smith, and O. B. Widlund. Schwarz Analysis of Iterative Substructuring Algorithms for Elliptic Problems in Three Dimensions, SIAM J. Numer. Anal., Vol. 31, No. 6, December 1994, Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
31 BDDC The idea is to use the same local bilinear forms to define the coarse basis functions The coarse basis functions are discontinuous across interfaces The local and global solvers are solved in parallel and the weights are used in the sum to make the iterative solutions continuous Reference: J. Mandel and C. R. Dohrmann, Convergence of a Balancing Domain Decomposition by Constraints and Energy Minimization, Numer. Lin. Alg. Appl. 10(2003) See also [TW] book Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
32 Partition of Unity Coarse Spaces Choose φ i as coarse basis functions ū i : Average of u in Ω δ i and define I PU H N u = ū i φ i i=1 For constant coefficients: H/δ bound. [S3] Obs: For φ i based on smoothed agglomeration techniques: (H/δ) 2 bound. [TW] For elasticity: coarse basis functions I h (φ i RBMs). RBMs are the rigid body motions. H/δ bound. [S2] Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
33 PU Coarse Spaces with Harmonic Overlap PU-ASHO: Restrict φ i to the Γ δ i := N i=1 Ωδ i and extend discrete harmonically elsewhere. We get H/δ bounds. [S3] Enhanced PU-Coarse Space: The bounds can be improved by considering coarse spaces based on I h (V Enh φ i ). In [S1]: For each subdomain Ω δ i we let the space V Enh be defined as the vector space generated by few lowest finite element eigenmodes associated to operator a Ω δ without assuming Dirichlet boundary condition on Ω δ i i PU-ASHO for discontinuous coefficients: For nodes x on Ω i define φ i (x) = ρ β i /( j ρβ j ). For nodes x on Γδ i, if inside Ω i define φ i (x) = 1, if outside φ i (x) = 0. The remaining nodes use discrete harmonic extension. [S1]. We get a (1 + H/h + (log H/h)(log δ/h)) bounds, the same bound as in RASHO [S4] PU-OBDD where the local problems are Neumann solvers on overlapping subdomais [S5] Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
34 References S1 M. Sarkis. Partition of unity coarse spaces: enhanced versions, discontinuous coefficients and applications to elasticity. Domain decomposition methods in science and engineering, Natl. Auton. Univ. Mex., Mexico, 2003, pp S2 M. Sarkis. A coarse space for elasticity: partition of unity rigid body motions coarse space. Applied Mathematics and Scientific Computing (Dubrovnik, 2001), Kluwer/Plenum, New York, 2003, pp S3 M. Sarkis. Partition of unity coarse spaces and Schwarz methods with harmonic overlap. Lect. Notes Comput. Sci. Eng. (Zurich, 2001), Springer-Verlag, Vol. 23, 2002, pp S4 X-C. Cai, M. Dryja and M. Sarkis. Restricted additive Schwarz preconditioner with harmonic overlap for symmetric positive definite linear systems. SIAM J. Numer. Anal., 41(4), 2003, pp S5 OBDD: overlapping balancing domain decomposition methods and generalizations to Helmholtz equations. LNCSE, vol. 55, Springer, 2006, pp Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
35 Galvis-Efendiev PU-Generalized Eigenvectors Overlapping regions Ω i : union of all Ω k sharing a vertex V i For each Ω i :, the enhanced space VEnh i : is defined as the span of all the eigenvectors associated to eigenvalues smaller than O(H) A i v (i,k) = λ (i,k) M i v (i,k) Coarse Space := Span (i,k) {I h (φ i v (i,k) )} Let PEnh i be the a-projection from V h(ω i ) to V Enh i (Ω i ) Key idea: (forany heterogeneous positive coefficients ρ) u P i Enh u 2 L 2 ρ(ω i ) H2 u 2 H 1 ρ(ω i ) P i Enh u 2 H 1 ρ(ω i ) u 2 H 1 ρ(ω i ) Variants: Mandel et al., Galvis-Efendiev et al., Nataf-Scheichl et al. Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
36 References M. Brezina, C. Heberton, J. Mandel, and P. Vanek, An iterative method with convergence rate chosen a priori, UCD/CCM Report 140, April An earlier version presented at the 1998 Copper Mountain Conference on Iterative Methods, April J. Galvis and Y. Efendiev, Domain Decomposition Preconditioners for Multiscale Flows in High-Contrast Media, Multiscale Model. Simul. Volume 8, Issue 4, pp (2010). J. Galvis and Y. Efendiev, Domain decomposition preconditioners for multiscale flows in high contrast media. Reduced dimension coarse spaces. Multiscale Model. Simul. Volume 8, Issue 5, pp (2010). N. Spillane, F. Nataf, V. Dolean, P. Hauret, C. Pechstein, and R. Scheichl and N. Spillane, Abstract Robust Coarse Spaces for Systems of PDEs via Generalized Eigenproblems in the Overlaps. Preprint NuMa Report No (2011). Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
37 A Few References on Heterogeneous Coefficients I.G. Graham, P. Lechner and R. Scheichl, Domain Decomposition for Multiscale PDEs, Numerische Mathematik 106: , 2007 R. Scheichl and E. Vainikko, Additive Schwarz with Aggregation-Based Coarsening for Elliptic Problems with Highly Variable Coefficients, Computing 80(4): , 2007 J. Van lent, R. Scheichl and I.G. Graham, Energy Minimizing Coarse Spaces for Two-level Schwarz Methods for Multiscale PDEs, Numerical Linear Algebra with Applications 16(10): , 2009 C. Pechstein and R. Scheichl, Weighted Poincare inequalities, submitted 10th December 2010, NuMa-Report , Institute of Computational Mathematics, Johannes Kepler University Linz, Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
38 References on Heterogeneous Coefficients (cont) A. Klawonn and O. Rheinbach, Robust FETI-DP methods for heterogeneous three dimensional elasticity problems, Comput. Meth. Appl. Mech. Engrg., Vol. 196, pp , January 2007 C. Pechstein and R. Scheichl, Analysis of FETI Methods for Multiscale PDEs - Part II: Interface Variation, Numerische Mathematik 118 (3): , 2011 M. Dryja and M. Sarkis. Boundary Layer Technical Tools for FETI-DP Methods to Heterogeneous Coefficients. In Domain Decomposition Methods in Science and Engineering XIX, Huang, Y.; Kornhuber, R.; Widlund, O.; Xu, J. (Eds.), Volume 78 of Lecture Notes in Computational Science and Engineering, Springer-Verlag, 2011, pp M. Dryja and M. Sarkis. Additive average Schwarz methods for discretization of elliptic problems with highly discontinuous coefficients. Computational Methods in Applied Mathematics, Vol 10 (2), pp , Marcus Sarkis (WPI) Domain Decomposition Methods RICAM / 38
Two 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 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 informationSubstructuring for multiscale problems
Substructuring for multiscale problems Clemens Pechstein Johannes Kepler University Linz (A) jointly with Rob Scheichl Marcus Sarkis Clark Dohrmann DD 21, Rennes, June 2012 Outline 1 Introduction 2 Weighted
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 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 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 informationAlgebraic Coarse Spaces for Overlapping Schwarz Preconditioners
Algebraic Coarse Spaces for Overlapping Schwarz Preconditioners 17 th International Conference on Domain Decomposition Methods St. Wolfgang/Strobl, Austria July 3-7, 2006 Clark R. Dohrmann Sandia National
More informationOn Iterative Substructuring Methods for Multiscale Problems
On Iterative Substructuring Methods for Multiscale Problems Clemens Pechstein Introduction Model Problem Let Ω R 2 or R 3 be a Lipschitz polytope with boundary Ω = Γ D Γ N, where Γ D Γ N = /0. We are interested
More informationDOMAIN DECOMPOSITION FOR LESS REGULAR SUBDOMAINS: OVERLAPPING SCHWARZ IN TWO DIMENSIONS TR
DOMAIN DECOMPOSITION FOR LESS REGULAR SUBDOMAINS: OVERLAPPING SCHWARZ IN TWO DIMENSIONS TR2007-888 CLARK R. DOHRMANN, AXEL KLAWONN, AND OLOF B. WIDLUND Abstract. In the theory of domain decomposition methods,
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 informationJOHANNES KEPLER UNIVERSITY LINZ. Weighted Poincaré Inequalities and Applications in Domain Decomposition
JOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics Weighted Poincaré Inequalities and Applications in Domain Decomposition Clemens Pechstein Institute of Computational Mathematics,
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 informationarxiv: v1 [math.na] 16 Dec 2015
ANALYSIS OF A NEW HARMONICALLY ENRICHED MULTISCALE COARSE SPACE FOR DOMAIN DECOMPOSITION METHODS MARTIN J. GANDER, ATLE LONELAND, AND TALAL RAHMAN arxiv:52.05285v [math.na] 6 Dec 205 Abstract. We propose
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 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 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 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 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 informationToward black-box adaptive domain decomposition methods
Toward black-box adaptive domain decomposition methods Frédéric Nataf Laboratory J.L. Lions (LJLL), CNRS, Alpines Inria and Univ. Paris VI joint work with Victorita Dolean (Univ. Nice Sophia-Antipolis)
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 informationAsymptotic expansions for high-contrast elliptic equations
Asymptotic expansions for high-contrast elliptic equations Victor M. Calo Yalchin Efendiev Juan Galvis 1 Abstract In this paper, we present a high-order expansion for elliptic equations in highcontrast
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 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 informationAdaptive Coarse Spaces and Multiple Search Directions: Tools for Robust Domain Decomposition Algorithms
Adaptive Coarse Spaces and Multiple Search Directions: Tools for Robust Domain Decomposition Algorithms Nicole Spillane Center for Mathematical Modelling at the Universidad de Chile in Santiago. July 9th,
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 informationDual-Primal Isogeometric Tearing and Interconnecting Solvers for Continuous and Discontinuous Galerkin IgA Equations
Dual-Primal Isogeometric Tearing and Interconnecting Solvers for Continuous and Discontinuous Galerkin IgA Equations Christoph Hofer and Ulrich Langer Doctoral Program Computational Mathematics Numerical
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 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 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 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 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 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 informationAccomodating Irregular Subdomains in Domain Decomposition Theory
Accomodating Irregular Subdomains in Domain Decomposition Theory Olof B. Widlund 1 Courant Institute, 251 Mercer Street, New York, NY 10012, USA widlund@cims.nyu.edu Summary. In the theory for domain decomposition
More informationTR THREE-LEVEL BDDC IN THREE DIMENSIONS
R2005-862 REE-LEVEL BDDC IN REE DIMENSIONS XUEMIN U Abstract BDDC methods are nonoverlapping iterative substructuring domain decomposition methods for the solution of large sparse linear algebraic systems
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 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 informationUne méthode parallèle hybride à deux niveaux interfacée dans un logiciel d éléments finis
Une méthode parallèle hybride à deux niveaux interfacée dans un logiciel d éléments finis Frédéric Nataf Laboratory J.L. Lions (LJLL), CNRS, Alpines Inria and Univ. Paris VI Victorita Dolean (Univ. Nice
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 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 informationFakultät für Mathematik und Informatik
Fakultät für Mathematik und Informatik Preprint 2016-09 Alexander Heinlein, Axel Klawonn, Jascha Knepper, Oliver Rheinbach Multiscale Coarse Spaces for Overlapping Schwarz Methods Based on the ACMS Space
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 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 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 informationMultigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids
Multigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids Long Chen 1, Ricardo H. Nochetto 2, and Chen-Song Zhang 3 1 Department of Mathematics, University of California at Irvine. chenlong@math.uci.edu
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 informationAlgebraic Adaptive Multipreconditioning applied to Restricted Additive Schwarz
Algebraic Adaptive Multipreconditioning applied to Restricted Additive Schwarz Nicole Spillane To cite this version: Nicole Spillane. Algebraic Adaptive Multipreconditioning applied to Restricted Additive
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 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 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 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 informationSOME PRACTICAL ASPECTS OF PARALLEL ADAPTIVE BDDC METHOD
Conference Applications of Mathematics 2012 in honor of the 60th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 2012 SOME PRACTICAL ASPECTS OF PARALLEL ADAPTIVE BDDC METHOD Jakub Šístek1,2,
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 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 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 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 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 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 informationBDDC deluxe for Isogeometric Analysis
BDDC deluxe for sogeometric Analysis L. Beirão da Veiga 1, L.. Pavarino 1, S. Scacchi 1, O. B. Widlund 2, and S. Zampini 3 1 ntroduction The main goal of this paper is to design, analyze, and test a BDDC
More informationAN ANALYSIS OF A FETI DP ALGORITHM ON IRREGULAR SUBDOMAINS IN THE PLANE TR
AN ANALYSIS OF A FETI DP ALGORITHM ON IRREGULAR SUBDOMAINS IN THE PLANE TR27 889 AXEL KLAWONN, OLIVER RHEINBACH, AND OLOF B. WIDLUND Abstract. In the theory for domain decomposition algorithms of the iterative
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 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 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 informationA Neumann-Dirichlet Preconditioner for FETI-DP 2 Method for Mortar Discretization of a Fourth Order 3 Problems in 2D 4 UNCORRECTED PROOF
1 A Neumann-Dirichlet Preconditioner for FETI-DP 2 Method for Mortar Discretization of a Fourth Order 3 Problems in 2D 4 Leszek Marcinkowski * 5 Faculty of Mathematics, University of Warsaw, Banacha 2,
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 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 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 informationSOME NONOVERLAPPING DOMAIN DECOMPOSITION METHODS
SIAM REV. c 1998 Society for Industrial and Applied Mathematics Vol. 40, No. 4, pp. 857 914, December 1998 004 SOME NONOVERLAPPING DOMAIN DECOMPOSITION METHODS JINCHAO XU AND JUN ZOU Abstract. The purpose
More informationJOHANNES KEPLER UNIVERSITY LINZ. Abstract Robust Coarse Spaces for Systems of PDEs via Generalized Eigenproblems in the Overlaps
JOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics Abstract Robust Coarse Spaces for Systems of PDEs via Generalized Eigenproblems in the Overlaps Nicole Spillane Frédérik Nataf Laboratoire
More informationApplication of Preconditioned Coupled FETI/BETI Solvers to 2D Magnetic Field Problems
Application of Preconditioned Coupled FETI/BETI Solvers to 2D Magnetic Field Problems U. Langer A. Pohoaţǎ O. Steinbach 27th September 2004 Institute of Computational Mathematics Johannes Kepler University
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 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 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 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 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 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 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 informationFrom the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes
From the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes Dylan Copeland 1, Ulrich Langer 2, and David Pusch 3 1 Institute of Computational Mathematics,
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 informationGRUPO DE GEOFÍSICA MATEMÁTICA Y COMPUTACIONAL MEMORIA Nº 8
GRUPO DE GEOFÍSICA MATEMÁTICA Y COMPUTACIONAL MEMORIA Nº 8 MÉXICO 2014 MEMORIAS GGMC INSTITUTO DE GEOFÍSICA UNAM MIEMBROS DEL GGMC Dr. Ismael Herrera Revilla iherrerarevilla@gmail.com Dr. Luis Miguel de
More informationDomain Decomposition solvers (FETI)
Domain Decomposition solvers (FETI) a random walk in history and some current trends Daniel J. Rixen Technische Universität München Institute of Applied Mechanics www.amm.mw.tum.de rixen@tum.de 8-10 October
More informationScalable Domain Decomposition Preconditioners For Heterogeneous Elliptic Problems
Scalable Domain Decomposition Preconditioners For Heterogeneous Elliptic Problems Pierre Jolivet, F. Hecht, F. Nataf, C. Prud homme Laboratoire Jacques-Louis Lions Laboratoire Jean Kuntzmann INRIA Rocquencourt
More informationA Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements
W I S S E N T E C H N I K L E I D E N S C H A F T A Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements Matthias Gsell and Olaf Steinbach Institute of Computational Mathematics
More informationFrom the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes
www.oeaw.ac.at From the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes D. Copeland, U. Langer, D. Pusch RICAM-Report 2008-10 www.ricam.oeaw.ac.at From the Boundary Element
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 informationJOHANNES KEPLER UNIVERSITY LINZ. A Robust FEM-BEM Solver for Time-Harmonic Eddy Current Problems
JOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics A Robust FEM-BEM Solver for Time-Harmonic Eddy Current Problems Michael Kolmbauer Institute of Computational Mathematics, Johannes
More informationarxiv: v2 [math.na] 17 Jun 2010
Numerische Mathematik manuscript No. (will be inserted by the editor) Local Multilevel Preconditioners for Elliptic Equations with Jump Coefficients on Bisection Grids Long Chen 1, Michael Holst 2, Jinchao
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 informationHybrid (DG) Methods for the Helmholtz Equation
Hybrid (DG) Methods for the Helmholtz Equation Joachim Schöberl Computational Mathematics in Engineering Institute for Analysis and Scientific Computing Vienna University of Technology Contributions by
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 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 informationNumerische Mathematik
umer. Math. 73: 149 167 (1996) umerische Mathematik c Springer-Verlag 1996 Overlapping Schwarz methods on unstructured meshes using non-matching coarse grids Tony F. Chan 1, Barry F. Smith 2, Jun Zou 3
More informationINTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR GRIDS. 1. Introduction
Trends in Mathematics Information Center for Mathematical Sciences Volume 9 Number 2 December 2006 Pages 0 INTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR
More 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 informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
More informationDISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS
DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS Sergey Korotov BCAM Basque Center for Applied Mathematics http://www.bcamath.org 1 The presentation is based on my collaboration with several
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 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 informationOn domain decomposition preconditioners for finite element approximations of the Helmholtz equation using absorption
On domain decomposition preconditioners for finite element approximations of the Helmholtz equation using absorption Ivan Graham and Euan Spence (Bath, UK) Collaborations with: Paul Childs (Emerson Roxar,
More informationElectronic Transactions on Numerical Analysis Volume 49, 2018
Electronic Transactions on Numerical Analysis Volume 49, 2018 Contents 1 Adaptive FETI-DP and BDDC methods with a generalized transformation of basis for heterogeneous problems. Axel Klawonn, Martin Kühn,
More informationCapacitance Matrix Method
Capacitance Matrix Method Marcus Sarkis New England Numerical Analysis Day at WPI, 2019 Thanks to Maksymilian Dryja, University of Warsaw 1 Outline 1 Timeline 2 Capacitance Matrix Method-CMM 3 Applications
More information