Optimized Schwarz methods with overlap for the Helmholtz equation
|
|
- Bathsheba McCarthy
- 5 years ago
- Views:
Transcription
1 Optimized Schwarz methods with overlap for the Helmholtz equation Martin J. Gander and Hui Zhang Introduction For the Helmholtz equation, simple absorbing conditions of the form n iω were proposed as transmission condition (TC) in Schwarz methods first without overlap in [4], and later also with overlap, see [3, ]. More advanced TCs can also be used, see e.g. [, 4, ]. Furthermore, parameters can be introduced into TCs and then optimized for rapid convergence, which led to the so called optimized Schwarz methods, see e.g. [6, 3] for elliptic equations. Without overlap, the parameters involved in some zero- and second-order TCs for the Helmholtz equation have been optimized in [8, 7]. With overlap, preliminary numerical studies of the parameters have been presented in [5, 9]. In this paper, we present the asymptotic solutions of the corresponding optimization problems with small overlap. We also compare the optimized parameters with other choices based on convergence factors and actual iteration numbers. We finally test for the first time Taylor second-order absorbing conditions for domain decomposition with overlap in the Helmholtz case. Schwarz Methods with Overlap As a model problem, we consider the Helmholtz equation in free space, (ω + )u = f (x,y), (x,y) R R d, equipped with the Sommerfeld radiation condition lim r d ( u d iω) = 0, r = x r r + y j. j= We decompose the domain into two overlapping subdomains Ω = (,L) R d and Ω = (0, ) R d with the overlap size L > 0. The Schwarz iteration reads ω u n+ + u n+ = f (x,y), (x,y) Ω, ( x + S )(u n+ )(L,y) = ( x + S )(u n )(L,y), y Rd, Section of Mathematics, University of Geneva, -4 rue du Liévre, Case postale 64, {martin.gander}{hui.zhang}@unige.ch
2 Martin J. Gander and Hui Zhang and ω u n+ + u n+ = f (x,y), (x,y) Ω, ( x + S )(u n+ )(0,y) = ( x + S )(u n )(0,y), y Rd, where S j, j =, are two linear operators in some trace spaces along {L} R d and {0} R d, respectively. For the analysis it suffices to consider by linearity the case f (x,y) = 0 and to analyze convergence to the zero solution. We take a Fourier transform in the y direction to obtain and (ω k )û n+ + xxû n+ = 0, x (,L), ( x + s )(û n+ )(L,k) = ( x + s )(û n )(L,k), (ω k )û n+ + xxû n+ = 0, x (0, ), ( x + s )(û n+ )(0,k) = ( x + s )(û n )(0,k), where k is the Fourier variable of y and s j denotes the symbol of S j. Since the Sommerfeld radiation condition excludes growing solutions as well as incoming modes at infinity we obtain the solutions û n+ (x,k) = û n+ (L,k)e λ(k)(x L), û n+ (x,k) = û n+ (0,k)e λ(k)x, where λ(k) denotes the root of the characteristic equation λ +(ω k ) = 0 with positive real part or negative imaginary part, { k λ(k) := ω for k > ω, i ω k for k < ω. () Substitution of the solutions into the transmission conditions yields û n+ (L,k) = s (k) λ(k) s (k)+λ(k) e λ(k)l û n (0,k), û n+ (0,k) = s (k) λ(k) s (k)+λ(k) e λ(k)l û n (L,k). By recursion we have û n+ (L,k) = ρ(k)û n (L,k) and û n+ (0,k) = ρ(k)û n (0,k), where the convergence factor ρ for a double iteration is defined by ρ(k) = s (k) λ(k) s (k) + λ(k) s(k) λ(k) s (k) + λ(k) e λ(k)l. () Setting the two complex parameters s = p iq and s = p iq, with p j,q j R, and inserting s, s and () into the convergence factor (), we find after simplifying
3 Optimized Schwarz methods with overlap for the Helmholtz equation 3 p +(q ω k ) p +(q ω k ), ρ(p,q, p,q,k) p = +(q + ω k ) p +(q + ω k ) k < ω, q +(p k ω ) q +(p k ω ) q +(p + k ω ) q +(p + k ω ) e 4λ(k)L, k > ω. As long as k = ω and p j,q j > 0, we have ρ <. Remark. It was shown in [6] that the two-sided operators S j = s j C can be transformed into the second-order operators S = S = r r y, with r = ω + s s, r =, (4) s + s s + s and the associated convergence factor for a single iteration is then given by which is just () with L replaced by L/. ρ(k) = s (k) λ(k) s (k) + λ(k) s(k) λ(k) s (k) + λ(k) e λ(k)l, (5) (3) 3 Optimized transmission conditions For simplicity, we consider p = q, p = q. Our goal is to find good parameters p, p such that the modulus of the squared convergence factor (3) is as small as possible over a range of frequencies k [k min,k ] [k +,k max ], where k < ω < k +. We require k to be away from ω because ρ = when k = ω, independently of what one chooses for the parameters p j and q j. Since in general we do not know how the Fourier coefficients of the initial error are distributed over the frequencies, we minimize ρ for the worst case, that is, we solve the min max problem ( ) argmin (p,p ) P max ρ(p, p, p, p,k), (6) k [k min,k ] [k +,k max ] where P is a certain search domain of the parameters. For well-posedness of the subdomain problems, we should choose P [0, ). The best approximation problem (6) is difficult to solve, and we only give asymptotic formulas for the parameters such that the convergence factor is as small as possible in different limiting processes in the mesh size h and the wave number ω. The proofs of the following theorems are beyond the scope of this short paper and will appear in [0]. Theorem. Let L = C L h, k max [C/h, ], C L,C,k min,k,k + and ω be positive and independent of h, k min < k < ω, k max > k + > ω and P = {(p, p ) 0 p p < }. Suppose h is small and k [k min,k ] [k +,k max ]. If we set
4 4 Martin J. Gander and Hui Zhang p = p = C/5 ω (4L) /5 / + o(h /5 ), p = p = C/5 ω (4L) 3/5 + o(h 3/5 ), (7) where C ω = min { ω k,k + ω }, then (3) is bounded by 4(4L C ω ) /5 + o(h /5 ). Moreover, any solution of (6) must satisfy (7). Theorem. Let L = C L h, h = C h /ω γ, γ, k max [C/h, ], δ ω = min{ω k,k + ω} with C L,C h,c,k and k + positive constants independent of ω, k min < k < ω, k max > k + > ω and P = {(p, p ) 0 p p < }. Suppose ω is large and k [k min,k ] [k +,k max ]. Then, for γ < 9/8 any solution of (6) satisfies p = δ 3/8 ω (ω/) 5/8 + o(ω 5/8 ), p = (δ ω) /8 ω 7/8 + o(ω 7/8 ), for which (3) is bounded by 4 /8 δω /8 ω /8 + o(ω /8 ). For γ > 9/8, any solution of (6) satisfies p = (δ ω ω) /5 L /5 / + o(ω /5+γ/5 ), p = (δ ω ω) /5 L 3/5 / + o(ω /5+3γ/5 ), and (3) is bounded by 4 (C h C L ) /5 δω /0 ω /0 γ/5 + o(ω /0 γ/5 ). Finally, for γ = 9/8, any solution of (6) satisfies p = (C hc L δ ω ) /3 ω 5/8 + o(ω 5/8 ), p = C hc L ω 7/8 + o(ω 7/8 ), and (3) is bounded by { 6C h C L ω /8 + o(ω /8 ), if 5/8 δω /8 C h C L, δω /6 C /3 h C /3 L ω /8 + o(ω /8 ), if 5/8 δω /8 C h C L. Remark. In the particular case γ = 9/8, the constant in front of the leading term of p can be an arbitrary number in the interval [ δ ω /(8C h C L ),3Ch 3C3 L ] in order to solve (6). But the choice in the above theorem is the best in the sense that it simultaneously minimizes the maximum of the other local but not global maxima. Remark 3. In practice, we use only the leading order terms of the optimized parameters. But it is also possible to extract higher order terms. Fig. shows the convergence factors of different Schwarz methods, obtained for the model problem in R, with ω = 0π and h = /00. The maximum of the convergence factors for double iterations over the chosen interval k [π,ω π] [ω + π, π/h] are.0, for the classical Schwarz method and Després method without overlap [4], for Després method with overlap [3, ], for the optimized Schwarz methods without overlap [7], and and for the same method with
5 Optimized Schwarz methods with overlap for the Helmholtz equation Classical Schwarz Despres method no overlap Despres method with overlap Optimized method no overlap Optimized method with overlap ρ ω k: y frequency Fig. Convergence factors of different Schwarz methods as functions of the Fourier parameter k, for ω = 0π, h = /00. The vertical lines indicate the ω, ω π and ω + π which are used to exclude a short interval for the optimized methods. overlap. The overlap size we chose is h, and we used the second-order formulation (4), (5) for the optimized methods. 4 Numerical experiments We used the ORAS formulation described in [3] for our implementation. As an alternative, one could also use a substructured formulation, see e.g. [9]. We implemented the second-order transmission conditions as indicated in Remark. We always solve the homogeneous equation with the zero solution and use a random initial guess to stimulate all frequencies. We use the domain decomposition Ω = (0, +h) (0,), Ω = ( h, ) (0,), and iterate until the relative residual is less than 0 8. We compare the overlapping Schwarz methods with optimized second-order transmission condition denoted by OO to those with the classical Dirichlet condition denoted by Cl, simple absorbing conditions of the form n iω (i.e. Després method with overlap, c.f. [3, ]) denoted by TO0, because it corresponds to a Taylor expansion of zero order of the symbol of the DtN operator, and the second-order low frequency absorbing condition, which is denoted by TO. Since the Schwarz methods can be used as a stationary iterative solver, or as a preconditioner for GMRES, both cases are tested, except for the classical Schwarz stationary iteration, which can not converge. We consider the open cavity problem with homogeneous Dirichlet boundary conditions on the top and the bottom of the unit square and the TO second-order absorbing conditions [] on the left and the right sides, and also the free space problem
6 6 Martin J. Gander and Hui Zhang Table Iteration numbers for the open cavity problem on the left and for the free space problem on the right, top half for ω = 9.5π, and below for ω = 0π. /h Stationary GMRES Stationary GMRES TO0 TO OO Cl. TO0 TO OO TO0 TO OO Cl. TO0 TO OO truncated to the unit square with the TO second-order absorbing conditions at the boundary. First, we fix ω = 9.5π (or ω = 0π) which are away from (or on) the sine frequencies at the continuous level in the y-direction. The corresponding iteration numbers are listed in Table. We can see that the minimum distance from ω to the frequencies at the discrete level in the y-direction plays an important role in all the stationary iterations while in the GMRES iterations this effect is only moderate. Fig. shows the asymptotic behavior of the different Schwarz methods as h 0, for the open cavity problem, and confirms our Fourier analysis results in Theorem. Now we fix hω or hω 3/ constant to see how the Schwarz methods behave for higher and higher wave numbers, which corresponds to Theorem. The iteration numbers are listed in Table. We see that the optimized method still converges faster than the others when used as a preconditioner for GMRES, while the stationary iterations Despres classical Krylov Despres Krylov optimized optimized Krylov h 5h 0.5 5h h Fig. Asymptotic behavior of the Schwarz methods for the open cavity, ω = 9.5π.
7 Optimized Schwarz methods with overlap for the Helmholtz equation 7 Table Iteration numbers for the open cavity problem on the left and for the free space problem on the right, top half for hω = π/5, and below for hω 3/ 3.5. /h Stationary GMRES Stationary GMRES TO0 TO OO Cl. TO0 TO OO TO0 TO OO Cl. TO0 TO OO iterations are again greatly affected by the discrete frequencies close to the wave number. The bars in the tables represent divergence. Next, we test the various Schwarz methods for an increasing number of subdomains. Since in most cases the stationary iterations diverge, we only show the GMRES iteration numbers in Table 3, where we use a bar to represent iteration numbers larger than We can see, neglecting the numbers in the parentheses for the moment, that all the methods deteriorate rapidly and the overlapping TO method outperforms the others eventually. Clearly the optimization of the two subdomain convergence factor does not predict well the optimal choice in the case of many subdomains for the Helmholtz equation. To partially improve the OO method, we introduce now two heuristics. First, we take δ ω = Nπ/ instead of δ ω = π in the former experiments, where N denotes the number of subdomains in the x-direction. Second, since the real parts of the parameters slow down the convergence for propagating modes, which becomes worse when the number of subdomains increases, we use s j = (/N i) p j ( j =,) instead of s j = ( i) p j. The new results are shown in the parentheses of Table 3, where the first numbers are obtained by using the two heuristics and the second numbers are from numerically optimized parameters based on a new many-subdomain Fourier analysis. But still, the low frequency Taylor conditions perform best in these exper- Table 3 Iteration numbers of GMRES, h = /56, ω = 5.π, overlap h. Sub. open cavity free space Cl. TO0 TO OO Cl. TO0 TO OO (45 40) (30 ) (9 88) (50 36) (83 88) (8 67) (77 66) (4 35) (638 60) (8 70) (859 54) (68 36)
8 8 Martin J. Gander and Hui Zhang iments. Our on-going work is to take a closer look at the multi-domain case and to seek better choices of parameters if it is possible. Acknowledgements The work was supported by the University of Geneva. The second author was also partially supported by the International Science and Technology Cooperation Program of China (00DFA4700). References. Bamberger, A., Joly, P., Roberts, J.E.: Second-order absorbing boundary conditions for the wave equation: a solution for the corner problem. SIAM J. Numer. Anal. 7, (990). Boubendir, Y., Antoine, X., Geuzaine, C.: A quasi-optimal non-overlapping domain decomposition algorithm for the Helmholtz equation. J. Comput. Phys. 3, 6 80 (0) 3. Cai, X.C., Casarin, M.A., Elliott Jr., F.W., Widlund, O.B.: Overlapping Schwarz algorithms for solving Helmholtz s equation. In: J. Mandel, C. Farhat, X.C. Cai (eds.) Domain Decomposition Methods 0, Contemporary Mathematics 8, Boulder, pp AMS (998) 4. Després, B.: Domain decomposition method and the Helmholtz problem. In: G.C. Cohen, L. Halpern, P. Joly (eds.) Mathematical and numerical aspects of wave propagation phenomena, Strasbourg, pp SIAM (99) 5. Gander, M.J.: Optimized Schwarz methods for Helmholtz problems. In: N. Debit, M. Garbey, R.H.W. Hoppe, et. al. (eds.) Domain Decomposition Methods in Science and Engineering, 3th International Conference on Domain Decomposition Methods, Barcelona, pp CIMNE (00) 6. Gander, M.J.: Optimized Schwarz methods. SIAM J. Numer. Anal. 44, (006) 7. Gander, M.J., Halpern, L., Magoulès, F.: An optimized Schwarz method with two-sided Robin transmission conditions for the Helmholtz equation. Int. J. Numer. Meth. Fluids 55, (007) 8. Gander, M.J., Magoulès, F., Nataf, F.: Optimized Schwarz methods without overlap for the Helmholtz equation. SIAM J. Sci. Comput. 4, (00) 9. Gander, M.J., Zhang, H.: Domain decomposition methods for the Helmholtz equation: a numerical investigation. In: R. Bank, M. Holst, J. Xu (eds.) Domain Decomposition Methods in Science and Engineering XX, Lecture Notes in Computational Science and Engineering, San Diego, p. 5. Springer-Verlag (0) 0. Gander, M.J., Zhang, H.: Optimized Schwarz methods with overlap for the Helmholtz equation. in preparation (0). Ghanemi, S.: A domain decomposition method for Helmholtz scattering problems. In: P.E. Bjorstad, M.S. Espedal, D.E. Keyes (eds.) Proceedings of the 9th Intl. Conf. on Domain Decomposition Methods, Ullensvang, pp. 05. ddm.org (998). Kimn, J.H., Sarkis, M.: Restricted overlapping balancing domain decomposition methods and restricted coarse problems for the Helmholtz problem. Comput. Methods Appl. Mech. Engrg. 96, (007) 3. St-Cyr, A., Gander, M.J., Thomas, S.J.: Optimized multiplicative, additive, and restricted additive Schwarz preconditioning. SIAM J. Sci. Comput. 9, (007) 4. Stupfel, B.: Improved transmission conditions for a one-dimensional domain decomposition method applied to the solution of the Helmholtz equation. J. Comput. Phys. 9, (00)
Convergence 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 informationOPTIMIZED SCHWARZ METHODS
SIAM J. NUMER. ANAL. Vol. 44, No., pp. 699 31 c 006 Society for Industrial and Applied Mathematics OPTIMIZED SCHWARZ METHODS MARTIN J. GANDER Abstract. Optimized Schwarz methods are a new class of Schwarz
More informationAdditive Schwarz method for scattering problems using the PML method at interfaces
Additive Schwarz method for scattering problems using the PML method at interfaces Achim Schädle 1 and Lin Zschiedrich 2 1 Zuse Institute, Takustr. 7, 14195 Berlin, Germany schaedle@zib.de 2 Zuse Institute,
More informationComparison of a One and Two Parameter Family of Transmission Conditions for Maxwell s Equations with Damping
Comparison of a One and Two Parameter Family of Transmission Conditions for Maxwell s Equations with Damping M. El Bouajaji 1, V. Dolean 2, M. J. Gander 3 and S. Lanteri 1 1 Introduction Transmission conditions
More informationUNCORRECTED PROOF. Domain Decomposition Methods for the Helmholtz 2 Equation: A Numerical Investigation 3. 1 Introduction 7
Domain Decomposition Methods for the Helmholtz 2 Equation: A Numerical Investigation 3 Martin J. Gander and Hui Zhang 4 University of Geneva martin.gander@unige.ch 5 2 hui.zhang@unige.ch 6 Introduction
More informationOptimized Schwarz Methods for Maxwell Equations with Discontinuous Coefficients
Optimized Schwarz Methods for Maxwell Equations with Discontinuous Coefficients Victorita Dolean, Martin Gander, Erwin Veneros To cite this version: Victorita Dolean, Martin Gander, Erwin Veneros. Optimized
More informationDG discretization of optimized Schwarz methods for Maxwell s equations
DG discretization of optimized Schwarz methods for Maxwell s equations Mohamed El Bouajaji, Victorita Dolean,3, Martin J. Gander 3, Stéphane Lanteri, and Ronan Perrussel 4 Introduction In the last decades,
More informationA Parallel Schwarz Method for Multiple Scattering Problems
A Parallel Schwarz Method for Multiple Scattering Problems Daisuke Koyama The University of Electro-Communications, Chofu, Japan, koyama@imuecacjp 1 Introduction Multiple scattering of waves is one of
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 informationAcceleration of a Domain Decomposition Method for Advection-Diffusion Problems
Acceleration of a Domain Decomposition Method for Advection-Diffusion Problems Gert Lube 1, Tobias Knopp 2, and Gerd Rapin 2 1 University of Göttingen, Institute of Numerical and Applied Mathematics (http://www.num.math.uni-goettingen.de/lube/)
More informationOptimal Interface Conditions for an Arbitrary Decomposition into Subdomains
Optimal Interface Conditions for an Arbitrary Decomposition into Subdomains Martin J. Gander and Felix Kwok Section de mathématiques, Université de Genève, Geneva CH-1211, Switzerland, Martin.Gander@unige.ch;
More 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 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 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 informationThe Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator
The Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator Martin J. Gander and Achim Schädle Mathematics Section, University of Geneva, CH-, Geneva, Switzerland, Martin.gander@unige.ch
More informationAdvection Diffusion Problems with Pure Advection Approximation in Subregions
Advection Diffusion Problems with Pure Advection Approximation in Subregions M. J. Gander, L. Halpern 2, C. Japhet 2, and V. Martin 3 Université de Genève, 2-4 rue du Lièvre, CP 64, CH-2 Genève. Switzerland.
More informationDirichlet-Neumann and Neumann-Neumann Methods
Dirichlet-Neumann and Neumann-Neumann Methods Felix Kwok Hong Kong Baptist University Introductory Domain Decomposition Short Course DD25, Memorial University of Newfoundland July 22, 2018 Outline Methods
More informationNew constructions of domain decomposition methods for systems of PDEs
New constructions of domain decomposition methods for systems of PDEs Nouvelles constructions de méthodes de décomposition de domaine pour des systèmes d équations aux dérivées partielles V. Dolean?? F.
More information2 CAI, KEYES AND MARCINKOWSKI proportional to the relative nonlinearity of the function; i.e., as the relative nonlinearity increases the domain of co
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2002; 00:1 6 [Version: 2000/07/27 v1.0] Nonlinear Additive Schwarz Preconditioners and Application in Computational Fluid
More informationOn the performance of the Algebraic Optimized Schwarz Methods with applications
On the performance of the Algebraic Optimized Schwarz Methods with applications Report 12-1-18 Department of Mathematics Temple University October 212 This report is available at http://www.math.temple.edu/~szyld
More informationAdditive Schwarz Methods with Nonreflecting Boundary Conditions for the Parallel Computation of Helmholtz Problems
Contemporary Mathematics Volume 18, 1998 B 0-818-0988-1-0305-9 Additive Schwarz Methods with Nonreflecting Boundary Conditions for the Parallel Computation of Helmholtz Problems Lois C. McInnes, Romeo
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 informationOptimized Schwarz Methods for Maxwell equations
Optimized Schwarz Methods for Maxwell equations Victorita Dolean, Martin Gander, Luca Gerardo-Giorda To cite this version: Victorita Dolean, Martin Gander, Luca Gerardo-Giorda. Optimized Schwarz Methods
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 informationTwo-level multiplicative domain decomposition algorithm for recovering the Lamé coefficient in biological tissues
Two-level multiplicative domain decomposition algorithm for recovering the Lamé coefficient in biological tissues Si Liu and Xiao-Chuan Cai Department of Applied Mathematics, University of Colorado at
More informationOptimized Schwarz methods in spherical geometry with an overset grid system
Optimized Schwarz methods in spherical geometry with an overset grid system J. Côté 1, M. J. Gander 2, L. Laayouni 3, and A. Qaddouri 4 1 Recherche en prévision numérique, Meteorological Service of Canada,
More informationA Quasi-Optimal Non-Overlapping Domain Decomposition Algorithm for the Helmholtz Equation
A Quasi-Optimal Non-Overlapping Domain Decomposition Algorithm for the Helmholtz Equation Y. Boubendir, X. Antoine, C. Geuzaine June 28, 2011 Abstract This paper presents a new non-overlapping domain decomposition
More informationA Quasi-Optimal Non-Overlapping Domain Decomposition Algorithm for the Helmholtz Equation
A Quasi-Optimal Non-Overlapping Domain Decomposition Algorithm for the Helmholtz Equation Y. Boubendir, X. Antoine, C. Geuzaine May 8, 2012 Abstract This paper presents a new non-overlapping domain decomposition
More informationSpace-Time Nonconforming Optimized Schwarz Waveform Relaxation for Heterogeneous Problems and General Geometries
Space-Time Nonconforming Optimized Schwarz Waveform Relaxation for Heterogeneous Problems and General Geometries Laurence Halpern, Caroline Japhet, and Jérémie Szeftel 3 LAGA, Université Paris XIII, Villetaneuse
More information27. A Non-Overlapping Optimized Schwarz Method which Converges with Arbitrarily Weak Dependence on h
Fourteent International Conference on Domain Decomposition Metods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 003 DDM.org 7. A Non-Overlapping Optimized Scwarz Metod wic Converges
More informationShifted Laplace and related preconditioning for the Helmholtz equation
Shifted Laplace and related preconditioning for the Helmholtz equation Ivan Graham and Euan Spence (Bath, UK) Collaborations with: Paul Childs (Schlumberger Gould Research), Martin Gander (Geneva) Douglas
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 informationTHE OPTIMISED SCHWARZ METHOD AND THE TWO-LAGRANGE MULTIPLIER METHOD FOR HETEROGENEOUS PROBLEMS. by Neil Greer
THE OPTIMISED SCHWARZ METHOD AND THE TWO-LAGRANGE MULTIPLIER METHOD FOR HETEROGENEOUS PROBLEMS by Neil Greer Submitted for the degree of Doctor of Philosophy Department of Mathematics School of Mathematical
More informationBETI for acoustic and electromagnetic scattering
BETI for acoustic and electromagnetic scattering O. Steinbach, M. Windisch Institut für Numerische Mathematik Technische Universität Graz Oberwolfach 18. Februar 2010 FWF-Project: Data-sparse Boundary
More 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 informationA Non-overlapping Quasi-optimal Optimized Schwarz 2 Domain Decomposition Algorithm for the Helmholtz 3 Equation 4 UNCORRECTED PROOF
1 A Non-overlapping Quasi-optimal Optimized Schwarz 2 Domain Decomposition Algorithm for the Helmholtz 3 Equation 4 Y. Boubendir 1, X. Antoine 2, and C. Geuzaine 3 5 1 Department of Mathematical Sciences
More informationTwo-level Domain decomposition preconditioning for the high-frequency time-harmonic Maxwell equations
Two-level Domain decomposition preconditioning for the high-frequency time-harmonic Maxwell equations Marcella Bonazzoli 2, Victorita Dolean 1,4, Ivan G. Graham 3, Euan A. Spence 3, Pierre-Henri Tournier
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 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 informationFrom the invention of the Schwarz method to the Best Current Methods for Oscillatory Problems: Part 2
and From the invention of the Schwarz method to the Best Current Methods for Oscillatory Problems: Part 2 martin.gander@unige.ch University of Geneva Woudschoten, October 2014 Are There Two Equations?
More informationElectromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems
Electromagnetic wave propagation ELEC 041-Modeling and design of electromagnetic systems EM wave propagation In general, open problems with a computation domain extending (in theory) to infinity not bounded
More informationA Nonoverlapping Subdomain Algorithm with Lagrange Multipliers and its Object Oriented Implementation for Interface Problems
Contemporary Mathematics Volume 8, 998 B 0-88-0988--03030- A Nonoverlapping Subdomain Algorithm with Lagrange Multipliers and its Object Oriented Implementation for Interface Problems Daoqi Yang. Introduction
More informationOverlapping Schwarz Waveform Relaxation for Parabolic Problems
Contemporary Mathematics Volume 218, 1998 B 0-8218-0988-1-03038-7 Overlapping Schwarz Waveform Relaxation for Parabolic Problems Martin J. Gander 1. Introduction We analyze a new domain decomposition algorithm
More informationDeflated Krylov Iterations in Domain Decomposition Methods
305 Deflated Krylov Iterations in Domain Decomposition Methods Y.L.Gurieva 1, V.P.Ilin 1,2, and D.V.Perevozkin 1 1 Introduction The goal of this research is an investigation of some advanced versions of
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 informationGenerating Equidistributed Meshes in 2D via Domain Decomposition
Generating Equidistributed Meshes in 2D via Domain Decomposition Ronald D. Haynes and Alexander J. M. Howse 2 Introduction There are many occasions when the use of a uniform spatial grid would be prohibitively
More informationIntroduction to Domain Decomposition Methods in the numerical approximation of PDEs. Luca Gerardo-Giorda
Introduction to Domain Decomposition Methods in the numerical approximation of PDEs Luca Gerardo-Giorda Plan 0 Motivation Non-overlapping domain decomposition methods Multi-domain formulation Variational
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 informationMulti-Domain Approaches for the Solution of High-Frequency Time-Harmonic Propagation Problems
Académie universitaire Wallonie Europe Université de Liège Faculté des Sciences Appliquées Collège de doctorat en Électricité, électronique et informatique Multi-Domain Approaches for the Solution of High-Frequency
More informationarxiv: v1 [math.na] 11 Jul 2018
An Additive Overlapping Domain Decomposition Method for the Helmholtz Equation Wei Leng Lili Ju arxiv:1807.04180v1 [math.na] 11 Jul 2018 July 12, 2018 Abstract In this paper, we propose and analyze an
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 informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
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 informationNewton additive and multiplicative Schwarz iterative methods
IMA Journal of Numerical Analysis (2008) 28, 143 161 doi:10.1093/imanum/drm015 Advance Access publication on July 16, 2007 Newton additive and multiplicative Schwarz iterative methods JOSEP ARNAL, VIOLETA
More informationTIME SPACE DOMAIN DECOMPOSITION AND REACTIVE TRANSPORT IN POROUS MEDIA
XVIII International Conference on Water Resources CMWR 2 J. Carrera (Ed) c CIMNE, Barcelona, 2 TIME SPACE DOMAIN DECOMPOSITION AND REACTIVE TRANSPORT IN POROUS MEDIA Florian Haeberlein and Anthony Michel
More informationLinear and Non-Linear Preconditioning
and Non- martin.gander@unige.ch University of Geneva June 2015 Invents an Iterative Method in a Letter (1823), in a letter to Gerling: in order to compute a least squares solution based on angle measurements
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 informationDomain decomposition methods via boundary integral equations
Domain decomposition methods via boundary integral equations G. C. Hsiao a O. Steinbach b W. L. Wendland b a Department of Mathematical Sciences, University of Delaware, Newark, Delaware 19716, USA. E
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 informationOptimized Schwarz Waveform Relaxation: Roots, Blossoms and Fruits
Optimized Schwarz Waveform Relaxation: Roots, Blossoms and Fruits Laurence Halpern 1 LAGA, Université Paris XIII, 99 Avenue J-B Clément, 93430 Villetaneuse, France, halpern@math.univ-paris13.fr 1 Introduction:
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 FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS
Proceedings of ALGORITMY 2005 pp. 222 229 A FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS ELENA BRAVERMAN, MOSHE ISRAELI, AND ALEXANDER SHERMAN Abstract. Based on a fast subtractional
More informationNonlinear Schwarz iterations with Reduced Rank Extrapolation. Sébastien Duminil, Hassane Sadok, and Daniel B. Szyld
Nonlinear Schwarz iterations with Reduced Rank Extrapolation Sébastien Duminil, Hassane Sadok, and Daniel B. Szyld Report 14-09-01 September 2014. Revised March 2015 Department of Mathematics Temple University
More informationOversampling for the partition of unity parallel finite element algorithm
International conference of Computational Mathematics and its application Urumqi, China, July 25-27, 203 Oversampling for the partition of unity parallel finite element algorithm HAIBIAO ZHENG School of
More informationarxiv: v1 [math.na] 14 May 2016
NONLINEAR PRECONDITIONING: HOW TO USE A NONLINEAR SCHWARZ METHOD TO PRECONDITION NEWTON S METHOD V. DOLEAN, M.J. GANDER, W. KHERIJI, F. KWOK, AND R. MASSON arxiv:65.449v [math.na] 4 May 6 Abstract. For
More informationPostprint.
http://www.diva-portal.org Postprint This is the accepted version of a chapter published in Domain Decomposition Methods in Science and Engineering XXI. Citation for the original published chapter: Gander,
More informationOpen-source finite element solver for domain decomposition problems
1/29 Open-source finite element solver for domain decomposition problems C. Geuzaine 1, X. Antoine 2,3, D. Colignon 1, M. El Bouajaji 3,2 and B. Thierry 4 1 - University of Liège, Belgium 2 - University
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 informationFinite Element Analysis of Acoustic Scattering
Frank Ihlenburg Finite Element Analysis of Acoustic Scattering With 88 Illustrations Springer Contents Preface vii 1 The Governing Equations of Time-Harmonic Wave Propagation, 1 1.1 Acoustic Waves 1 1.1.1
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 informationAltered Jacobian Newton Iterative Method for Nonlinear Elliptic Problems
Altered Jacobian Newton Iterative Method for Nonlinear Elliptic Problems Sanjay K Khattri Abstract We present an Altered Jacobian Newton Iterative Method for solving nonlinear elliptic problems Effectiveness
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 informationIntegration of Sequential Quadratic Programming and Domain Decomposition Methods for Nonlinear Optimal Control Problems
Integration of Sequential Quadratic Programming and Domain Decomposition Methods for Nonlinear Optimal Control Problems Matthias Heinkenschloss 1 and Denis Ridzal 2 1 Department of Computational and Applied
More informationSome Geometric and Algebraic Aspects of Domain Decomposition Methods
Some Geometric and Algebraic Aspects of Domain Decomposition Methods D.S.Butyugin 1, Y.L.Gurieva 1, V.P.Ilin 1,2, and D.V.Perevozkin 1 Abstract Some geometric and algebraic aspects of various domain decomposition
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 informationModelling in photonic crystal structures
Modelling in photonic crystal structures Kersten Schmidt MATHEON Nachwuchsgruppe Multiscale Modelling and Scientific Computing with PDEs in collaboration with Dirk Klindworth (MATHEON, TU Berlin) Holger
More informationANALYSIS AND NUMERICAL METHODS FOR SOME CRACK PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING, SERIES B Volume 2, Number 2-3, Pages 155 166 c 2011 Institute for Scientific Computing and Information ANALYSIS AND NUMERICAL METHODS FOR SOME
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 informationSOLVING MESH EIGENPROBLEMS WITH MULTIGRID EFFICIENCY
SOLVING MESH EIGENPROBLEMS WITH MULTIGRID EFFICIENCY KLAUS NEYMEYR ABSTRACT. Multigrid techniques can successfully be applied to mesh eigenvalue problems for elliptic differential operators. They allow
More informationA SIMPLE PRECONDITIONED DOMAIN DECOMPOSITION METHOD FOR ELECTROMAGNETIC SCATTERING PROBLEMS *
Journal of Computational Mathematics Vol.3, No., 203, 2. http://www.global-sci.org/jcm doi:0.4208/jcm.20-m3965 A SIMPLE PRECONDITIONED DOMAIN DECOMPOSITION METHOD FOR ELECTROMAGNETIC SCATTERING PROBLEMS
More informationDomain decomposition schemes with high-order accuracy and unconditional stability
Domain decomposition schemes with high-order accuracy and unconditional stability Wenrui Hao Shaohong Zhu March 7, 0 Abstract Parallel finite difference schemes with high-order accuracy and unconditional
More informationA Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation
A Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation Tao Zhao 1, Feng-Nan Hwang 2 and Xiao-Chuan Cai 3 Abstract In this paper, we develop an overlapping domain decomposition
More informationNew Approach for Accelerating Nonlinear Schwarz Iterations
Bol Soc Paran Mat (3s) v 00 0 (0000):???? c SPM ISSN-2175-1188 on line ISSN-00378712 in press SPM: wwwspmuembr/bspm doi:105269/bspmv38i437018 New Approach for Accelerating Nonlinear Schwarz Iterations
More informationNEAR FIELD REPRESENTATIONS OF THE ACOUSTIC GREEN S FUNCTION IN A SHALLOW OCEAN WITH FLUID-LIKE SEABED
Georgian Mathematical Journal Volume 4 2007, Number, 09 22 NEAR FIELD REPRESENTATIONS OF THE ACOUSTIC GREEN S FUNCTION IN A SHALLOW OCEAN WITH FLUID-LIKE SEAED ROERT GILERT AND MIAO-JUNG OU Abstract. In
More informationA PRECONDITIONER FOR THE HELMHOLTZ EQUATION WITH PERFECTLY MATCHED LAYER
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 A PRECONDITIONER FOR THE HELMHOLTZ EQUATION WITH PERFECTLY
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 informationA new method for the solution of scattering problems
A new method for the solution of scattering problems Thorsten Hohage, Frank Schmidt and Lin Zschiedrich Konrad-Zuse-Zentrum Berlin, hohage@zibde * after February 22: University of Göttingen Abstract We
More informationAn Introduction to Domain Decomposition Methods: algorithms, theory and parallel implementation
An Introduction to Domain Decomposition Methods: algorithms, theory and parallel implementation Victorita Dolean, Pierre Jolivet, Frédéric ataf To cite this version: Victorita Dolean, Pierre Jolivet, Frédéric
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 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 informationLinear and Non-Linear Preconditioning
and Non- martin.gander@unige.ch University of Geneva July 2015 Joint work with Victorita Dolean, Walid Kheriji, Felix Kwok and Roland Masson (1845): First Idea of After preconditioning, it takes only three
More informationNONCONFORMING IN TIME DOMAIN DECOMPOSITION METHOD FOR POROUS MEDIA APPLICATIONS
V European Conference on Computational Fluid Dynamics ECCOMAS CFD 2010 J. C. F. Pereira and A. Sequeira Eds Lisbon, Portugal,14-17 June 2010 NONCONFORMING IN TIME DOMAIN DECOMPOSITION METHOD FOR POROUS
More informationEfficient Numerical Simulation for Long Range Wave Propagation 1
Efficient Numerical Simulation for Long Range Wave Propagation 1 Kai Huang 2 George Papanicolaou 3 Knut Solna 2 Chrysoula Tsogka 4 Hongkai Zhao 2 1 The research is partially supported by ONR Grant N00014-02-1-0090,
More informationLeast-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations
Least-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations by Wilhelm Heinrichs Universität Duisburg Essen, Ingenieurmathematik Universitätsstr.
More 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 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 informationNumerical Methods for Large-Scale Nonlinear Systems
Numerical Methods for Large-Scale Nonlinear Systems Handouts by Ronald H.W. Hoppe following the monograph P. Deuflhard Newton Methods for Nonlinear Problems Springer, Berlin-Heidelberg-New York, 2004 Num.
More informationWHEN studying distributed simulations of power systems,
1096 IEEE TRANSACTIONS ON POWER SYSTEMS, VOL 21, NO 3, AUGUST 2006 A Jacobian-Free Newton-GMRES(m) Method with Adaptive Preconditioner and Its Application for Power Flow Calculations Ying Chen and Chen
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 informationThe Deflation Accelerated Schwarz Method for CFD
The Deflation Accelerated Schwarz Method for CFD J. Verkaik 1, C. Vuik 2,, B.D. Paarhuis 1, and A. Twerda 1 1 TNO Science and Industry, Stieltjesweg 1, P.O. Box 155, 2600 AD Delft, The Netherlands 2 Delft
More information