BETI for acoustic and electromagnetic scattering
|
|
- Kristin Benson
- 5 years ago
- Views:
Transcription
1 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 and Finite Element Domain Decomposition Methods in Electromagnetics 1 / 38
2 Outline Model problem The Tearing and Interconnecting approach Preconditioning Outlook 2 / 38
3 Outline Model problem The Tearing and Interconnecting approach Preconditioning Outlook 3 / 38
4 Problem definition Find u : Ω R so that holds. Remark: u(x) k 2 u(x) = 0 for x Ω, γ D u(x) = g(x) for x Γ D, γ N u(x) = u n (x) = h(x) for x Γ N, ) (ikeγ D u(x) + γ N u(x) = k(x) for x Γ R. k(x) piecewise constant Extension to exterior boundary value problems 4 / 38
5 Domain: Geometric domain decomposition Ω = i I Ω i with Ω i Ω j = if i j. Ω i is bounded for i > 0, Ω 0 may be unbounded Boundary: Ω := Γ D Γ N, Γ I := Γ ij, i,j I Γ S := Γ i = Γ I Ω. i I Γ 23 Γ D Ω 2 Γ 12 Ω 1 Γ 34 Γ 14 Ω 3 Ω 4 Γ N 5 / 38
6 Dirichlet to Neumann map Helmholtz u k 2 u = 0 in Ω, u = g on Γ, t = Su BEM: (symmetric analytic version) S k := D k + ( 1 2 I + K k 1 )Vk ( 1 2 I + K k) FEM: (discrete) S k := A ΓΓ A ΓI A 1 II A IΓ Remark: These definitions do not allow an exact discretization of the Steklov-Poincare operator. If k 2 is not a Dirichlet eigenfrequency then the Steklov-Poincare operator is well defined, bounded and fulfills a Garding inequality. If k 2 is not a Neumann eigenfrequency then the Steklov-Poincare operator is invertible. 6 / 38
7 Formulation of the domain decomposition approach Find u H 1 (Ω) such that u i k 2 u i = 0 in Ω i, γ N u i = h on Γ N, γ N u i = γ N u j on Γ ij where u i = u Ωi. By using the Steklov-Poincare operator, it is possible to reduce the problem to functions defined on the boundary: Find γ D u H 1/2 (Γ S ) such that S i γ D u Γi + S j γ D u Γj =0 on Γ ij S i γ D u = h on Γ N Γ I. 7 / 38
8 Formulation of the domain decomposition approach Find γ D u i H 1/2 (Γ i ) such that S i γ D u i + S j γ D u j = 0 on Γ ij, γ D u i =γ D u j on Γ ij, If we rewrite the interface conditions we get: γ N u = h on Γ N. γ N u i + γ N u j =0 on Γ ij, γ D u i γ D u j =0 on Γ ij This can be replaced by a Dirichlet-Robin type coupling condition [Després, 1990, Gander, 2002]: (γ N u i + ikpγ D u i ) + (γ N u j ikpγ D u j ) =0 on Γ ij, γ D u i γ D u j =0 with an elliptic operator P : H 1/2 (Γ ij ) H 1/2 (Γ ij ). on Γ ij 8 / 38
9 Outline Model problem The Tearing and Interconnecting approach Preconditioning Outlook 9 / 38
10 Basic idea of Tearing and Interconnecting Ω 2 Ω 1 Ω 5 Ω 3 Ω 4 All degrees of freedom are global degrees of freedom (red) [Farhat and Roux, 1991] 10 / 38
11 Basic idea of Tearing and Interconnecting Ω 2 Ω 1 Ω 5 Ω 3 Ω 4 All degrees of freedom are global degrees of freedom (red) Now we tear the domains more degrees of freedom (but local ones) [Farhat and Roux, 1991] 10 / 38
12 Basic idea of Tearing and Interconnecting Ω 2 Ω 1 Ω 5 Ω 3 Ω 4 All degrees of freedom are global degrees of freedom (red) Now we tear the domains more degrees of freedom (but local ones) Now we connect corresponding d.o.f. again by Lagrangian multipliers [Farhat and Roux, 1991] 10 / 38
13 Tearing and Interconnecting Find γ D u H 1/2 (Γ S ) so that S i γ D u Γi + S j γ D u Γj = 0 on Γ ij, S i γ D u =h(x) on Γ N Γ ij. After a summation and discretization we end up with the linear system S h u = p A T i S h,i A i u = i=1 p A T i f i. i=1 Now we can apply the tearing and interconnecting idea and get: S h,1 B1 T u 1 f S h,p Bp T u p =. f p. B 1... B p 0 λ 0 11 / 38
14 Dirichlet-to-Robin map After eliminating the primal degrees of freedom we get Remarks: p i=1 B 1 S 1 h,i BT i λ = rhs. S i may be not well defined S h,i may be not invertible Idea: Using Robin boundary conditions instead of Neumann boundary conditions for coupling [Benamou and Desprès, 1997]: γ N u i + ikpγ D u i = g on Γ R This leads to the formulation (S i + ikp)γ D u Γi + (S j ikp)γ D u Γj =0 on Γ ij S i γ D u = h on Γ N Γ i. 12 / 38
15 Dirichlet-to-Robin map As in the Dirichlet to Neumann case, this leads to p A T i S h,i A i u = i=1 p A T i f i. But the discrete splitting leads to S h,1 ± ikp h,1 B1 T u 1 f S h,p ± ikp h,p Bp T u p =. f p. B 1... B p 0 λ 0 The local mapping operators S j ± ikp j are injective and fulfill a Garding inequality, so the local matrices S h,j ± ikp h,j are invertible if h is small enough. i=1 13 / 38
16 Robin-System In general we have to solve local systems like (S + ip)u = (D + ip)u + ( 1 2 I + K )V 1 ( 1 I + K)u = rhs (1) 2 but this operator is not well defined if k 2 is a local Dirichlet eigenvalue. Instead we can solve the system ( D + ip 1 2 I + K ) ( ( ) u rhs v) 1 2 I K V =. 0 Properties: Equivalent to system (1) if k 2 is not a Dirichlet eigenvalue Garding inequality + surjective uniquely solvable Covers natural properties of the boundary value problem 14 / 38
17 Robin-System Global problem: Find γ D u H 1/2 (Γ S ) and γ N u j H 1/2 (Γ j ) such that p ( Dj + ip j ( 1 2 I + K ) ( ) ( ) ) j γd u Γj γd v ( 1 2 I + K), Γj j V j γ N u j γ N v j j=1 = p j=1 ( ) ( ) rhsj γd v, Γj 0 γ N v j for all γ D v H 1/2 (Γ S ) and γ N v j H 1/2 (Γ j ) Properties: Fulfills a Garding inequality Injective 15 / 38
18 Robin-System After Tearing and Interconnecting (with 1 2 I +K = K and D +ip = D) we get ( ) ( ) ( ) D h K h B T 10 γd u ( ) 1 f K h V h 1 γ N u ( ) ( ). D h K h B T ( ) p γd u =. ( ). p f p K h V h 0 ( p B1 0 ) (... Bp 0 ) γ N u p 0 0 λ 0 After eliminating the primal degrees of freedom we get Fλ = p i=1 ( Bi 0 ) ( D h K h K h V h ) 1 i ( B T i 0 ) λ = rhs. 16 / 38
19 Example Neumann boundary value problem: Domain: u(x) k 2 u(x) = 0 for x Ω, γ N u(x) = h(x) for x Γ. Ω 0 Ω 1 Frequencies: k = 2.0 and k = (first eigenfrequency of the unit cube). 17 / 38
20 Robin-System Frequency: 2.0 Frequency: N i M i GMRES rel. error u 1 u 1,h L2(Γ 1) N i M i GMRES rel. error u 1 u 1,h L2(Γ 1) / 38
21 Example Neumann boundary value problem: Domain: Unit cube u(x) k 2 u(x) = 0 for x Ω, γ N u(x) = h(x) for x Γ. Frequencies: k = / 38
22 Robin-System p N i M i GMRES u 1 u 1,h L2 (Γ 1 ) / 38
23 Outline Model problem The Tearing and Interconnecting approach Preconditioning Outlook 21 / 38
24 Local preconditioning To precondition the system ( Dh + ip A = h ( 1 2 M ) h + K h ) ( 1 2 M h + K h ) V h we use a block diagonal preconditioner based on on the idea of operators of opposite order [Steinbach and Wendland, 1998]: with C 1 A ( M 1 = 1,h V h M 1 1,h M 1 0,h D h M 1 0,h ) V h [i, j] = Vψ i, ψ j D h[i, j] = Vφ i, φ j ψ i, ψ j S 1 h (Γ), φ i, φ i S 0 h(γ). 22 / 38
25 Local preconditioning We solve now C 1 A A = rhs on the unit cube. Iteration numbers in dependency of k and h. N i M i Remark: Iteration numbers depend on the domain Ω 23 / 38
26 Global Preconditioning Idea from FETI-H method [Farhat et al., 2000] The residual r = d Fλ should be orthogonal to a given vector space in each step: Q t r = 0. (2) This is also a solution constraint: r = d Fλ = n B i u i. i=1 So r is the jump of the solution on the interface. 24 / 38
27 Global Preconditioning Introduce new iterate λ = λ + µ = λ + Qγ (3) where γ is chosen such that Q t r = 0. Inserting (3) in (2) leads to Q t FQγ = Q t (d Fλ). Solving this equation and inserting it in (3) leads to λ = Pλ + λ with P = I Q(Q t FQ) 1 Q t F, λ = Q(Q t FQ) 1 Q t d. 25 / 38
28 Global Preconditioning Replacing λ by λ in Fλ = d leads to FPλ + Fλ = d By multiplying this equation with P t we get P t FPλ = P t d since P t Fλ = 0. Finally we have to solve this and the coarse system Q t FQγ = Q t d Remark: P t FP = FP 26 / 38
29 Local planar waves Planar waves: Column-wise construction of Q: Q j ( ) = e ikθt i x( ) with Q j j-th column of Q θ j j-th unitary direction vector x( ) vector of nodal point coordinates evaluated at ( ) Local waves: The same approach on subdomain level Q = [Q 1 Q s Q Ns ] 27 / 38
30 Domain: Unit cube (0, 1) 3 Frequency: k = 2.0 Numerical example Iteration number Number of subdomains Number of directions per subdomain Refinement level 28 / 38
31 Domain: Unit cube (0, 1) 3 Frequency: k = 4.0 Numerical example Iteration number Number of subdomains Number of directions per subdomain Refinement level 29 / 38
32 Domain: Unit cube (0, 1) 3 Frequency: k = 8.0 Numerical example Iteration number Number of subdomains Number of directions per subdomain Refinement level 30 / 38
33 Outline Model problem The Tearing and Interconnecting approach Preconditioning Outlook 31 / 38
34 Helmholtz-Maxwell Helmholtz Maxwell Equation scalar vectorial Operator k 2 I curlcurl k 2 I S k,d k Garding inequality generalized Garding inequality Energy space H 1 (Ω) H(curl,Ω) Dirichlet trace space H 1/2 (Γ) H 1/2 (curl Γ, Γ) Neumann trace space H 1/2 (Γ) H 1/2 (div Γ, Γ) Discretization Ph 0 and P1 h Raviart-Thomas elements 32 / 38
35 Very first Maxwell example Local problem Number of iterations to solve ( 1 D h 2 M ) ( ) h + K h x ( 1 2 M = rhs h + K h ) V h y with given right hand side rhs on the unit cube (0, 1) 3 for various frequencies / 38
36 Very first Maxwell example Neumann boundary value problem: curlcurlu(x) k 2 u(x) = 0 for x Ω, γ N u(x) = h(x) for x Γ. Domain: Unit cube (0, 1) 3, divided in p 3 subcubes it Table: 8 domains it Table: 27 domains No preconditioning No Robin interfaces 34 / 38
37 Very first Maxwell example Neumann boundary value problem: curlcurlu(x) k 2 u(x) = 0 for x Ω, γ N u(x) = h(x) for x Γ. Domain: Unit cube (0, 1) 3, divided in p 3 subcubes er Table: 8 domains er Table: 27 domains No preconditioning No Robin interfaces 35 / 38
38 Outlook Analysis for the Maxwell case Numerical Analysis for the Maxwell case Preconditioning for Maxwell Implementing remaining parts for Maxwell Preprint 2009/05. O. Steinbach, M. Windisch Stable Boundary element domain decomposition methods for the Helmholtz problem Proceedings DD19. O. Steinbach, M. Windisch Robust BE Domain Decomposition Methods in Acoustics 36 / 38
39 Benamou, J.-D. and Desprès, B. (1997). A domain decomposition method for the helmholtz equation and related optimal control problems. J. Comput. Phys., 136(1): Després, B. (1990). Décomposition de domaine et problème de Helmholtz. C. R. Acad. Sci. Paris Sér. I Math., 311(6): Farhat, C., Macedo, A., and Lesoinne, M. (2000). A two-level domain decomposition method for the iterative solution of high frequency exterior Helmholtz problems. Numer. Math., 85(2): Farhat, C. and Roux, F.-X. (1991). A method of finite element tearing and interconnecting and its parallel solution algorithm. 32(6): / 38
40 Gander, M. J. (2002). Optimized Schwarz methods for Helmholtz problems. In Domain decomposition methods in science and engineering (Lyon, 2000), Theory Eng. Appl. Comput. Methods, pages Internat. Center Numer. Methods Eng. (CIMNE), Barcelona. Steinbach, O. and Wendland, W. L. (1998). The construction of some efficient preconditioners in the boundary element method. Adv. Comput. Math., 9(1-2): Numerical treatment of boundary integral equations. 38 / 38
Technische Universität Graz
Technische Universität Graz Robust boundary element domain decomposition solvers in acoustics O. Steinbach, M. Windisch Berichte aus dem Institut für Numerische Mathematik Bericht 2009/9 Technische Universität
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 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 informationTechnische Universität Graz
Technische Universität Graz A note on the stable coupling of finite and boundary elements O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2009/4 Technische Universität Graz A
More informationInexact Data-Sparse BETI Methods by Ulrich Langer. (joint talk with G. Of, O. Steinbach and W. Zulehner)
Inexact Data-Sparse BETI Methods by Ulrich Langer (joint talk with G. Of, O. Steinbach and W. Zulehner) Radon Institute for Computational and Applied Mathematics Austrian Academy of Sciences http://www.ricam.oeaw.ac.at
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 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 informationDomain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions
Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions Ernst P. Stephan 1, Matthias Maischak 2, and Thanh Tran 3 1 Institut für Angewandte Mathematik, Leibniz
More informationA 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 informationCoupled FETI/BETI for Nonlinear Potential Problems
Coupled FETI/BETI for Nonlinear Potential Problems U. Langer 1 C. Pechstein 1 A. Pohoaţǎ 1 1 Institute of Computational Mathematics Johannes Kepler University Linz {ulanger,pechstein,pohoata}@numa.uni-linz.ac.at
More informationTechnische Universität Graz
Technische Universität Graz Modified combined field integral equations for electromagnetic scattering O. Steinbach, M. Windisch Berichte aus dem Institut für Numerische Mathematik Bericht 2007/6 Technische
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 informationTechnische Universität Graz
Technische Universität Graz The All-Floating Boundary Element Tearing and Interconnecting Method G. Of, O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2009/3 Technische Universität
More informationFast Multipole BEM for Structural Acoustics Simulation
Fast Boundary Element Methods in Industrial Applications Fast Multipole BEM for Structural Acoustics Simulation Matthias Fischer and Lothar Gaul Institut A für Mechanik, Universität Stuttgart, Germany
More informationScalable BETI for Variational Inequalities
Scalable BETI for Variational Inequalities Jiří Bouchala, Zdeněk Dostál and Marie Sadowská Department of Applied Mathematics, Faculty of Electrical Engineering and Computer Science, VŠB-Technical University
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 informationSharp condition number estimates for the symmetric 2-Lagrange multiplier method
Sharp condition number estimates for the symmetric -Lagrange multiplier method Stephen W. Drury and Sébastien Loisel Abstract Domain decomposition methods are used to find the numerical solution of large
More 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 informationFETI Methods for the Simulation of Biological Tissues
SpezialForschungsBereich F 32 Karl Franzens Universita t Graz Technische Universita t Graz Medizinische Universita t Graz FETI Methods for the Simulation of Biological Tissues Ch. Augustin O. Steinbach
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 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 informationThe Helmholtz Equation
The Helmholtz Equation Seminar BEM on Wave Scattering Rene Rühr ETH Zürich October 28, 2010 Outline Steklov-Poincare Operator Helmholtz Equation: From the Wave equation to Radiation condition Uniqueness
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 informationTechnische Universität Graz
Technische Universität Graz Inexact Data Sparse Boundary Element Tearing and Interconnecting Methods U. Langer, G. Of, O. Steinbach, W. Zulehner Berichte aus dem Institut für Mathematik D Numerik und Partielle
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 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 informationInexact Data Sparse Boundary Element Tearing and Interconnecting Methods
Inexact Data Sparse Boundary Element Tearing and Interconnecting Methods U. Langer 1,2, G. Of 3, O. Steinbach 4, W. Zulehner 2 1 Radon Institute for Computational and Applied Mathematics, Austrian Academy
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 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 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 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 informationAll-Floating Coupled Data-Sparse Boundary and Interface-Concentrated Finite Element Tearing and Interconnecting Methods
All-Floating Coupled Data-Sparse Boundary and Interface-Concentrated Finite Element Tearing and Interconnecting Methods Ulrich Langer 1,2 and Clemens Pechstein 2 1 Institute of Computational Mathematics,
More informationAn Iterative Domain Decomposition Method for the Solution of a Class of Indefinite Problems in Computational Structural Dynamics
An Iterative Domain Decomposition Method for the Solution of a Class of Indefinite Problems in Computational Structural Dynamics Charbel Farhat a, and Jing Li a a Department of Aerospace Engineering Sciences
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 informationNumerical Approximation Methods for Elliptic Boundary Value Problems
Numerical Approximation Methods for Elliptic Boundary Value Problems Olaf Steinbach Numerical Approximation Methods for Elliptic Boundary Value Problems Finite and Boundary Elements Olaf Steinbach Institute
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 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 informationPreconditioned space-time boundary element methods for the heat equation
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 Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More informationFETI-DPH: A DUAL-PRIMAL DOMAIN DECOMPOSITION METHOD FOR ACOUSTIC SCATTERING
Journal of Computational Acoustics, c IMACS FETI-DPH: A DUAL-PRIMAL DOMAIN DECOMPOSITION METHOD FOR ACOUSTIC SCATTERING Charbel Farhat, Philip Avery and Radek Tezaur Department of Mechanical Engineering
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 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 informationFEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment
FEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment Alastair J. Radcliffe Andreas Dedner Timo Betcke Warwick University, Coventry University College of London (UCL) U.K. Radcliffe
More informationNumerical approximation of output functionals for Maxwell equations
Numerical approximation of output functionals for Maxwell equations Ferenc Izsák ELTE, Budapest University of Twente, Enschede 11 September 2004 MAXWELL EQUATIONS Assumption: electric field ( electromagnetic
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 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 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 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 informationAn Efficient FETI Implementation on Distributed Shared Memory Machines with Independent Numbers of Subdomains and Processors
Contemporary Mathematics Volume 218, 1998 B 0-8218-0988-1-03024-7 An Efficient FETI Implementation on Distributed Shared Memory Machines with Independent Numbers of Subdomains and Processors Michel Lesoinne
More informationParallel Scalability of a FETI DP Mortar Method for Problems with Discontinuous Coefficients
Parallel Scalability of a FETI DP Mortar Method for Problems with Discontinuous Coefficients Nina Dokeva and Wlodek Proskurowski Department of Mathematics, University of Southern California, Los Angeles,
More 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 informationChapter 1 Mathematical Foundations
Computational Electromagnetics; Chapter 1 1 Chapter 1 Mathematical Foundations 1.1 Maxwell s Equations Electromagnetic phenomena can be described by the electric field E, the electric induction D, the
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 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 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 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 informationFrom Direct to Iterative Substructuring: some Parallel Experiences in 2 and 3D
From Direct to Iterative Substructuring: some Parallel Experiences in 2 and 3D Luc Giraud N7-IRIT, Toulouse MUMPS Day October 24, 2006, ENS-INRIA, Lyon, France Outline 1 General Framework 2 The direct
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 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 informationAn Adaptive Space-Time Boundary Element Method for the Wave Equation
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 An Adaptive Space-Time Boundary Element Method for the Wave Equation Marco Zank and Olaf Steinbach Institut für Numerische Mathematik AANMPDE(JS)-9-16,
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 informationarxiv: v2 [cs.ce] 22 Oct 2016
SIMULATION OF ELECTRICAL MACHINES A FEM-BEM COUPLING SCHEME LARS KIELHORN, THOMAS RÜBERG, JÜRGEN ZECHNER arxiv:1610.05472v2 [cs.ce] 22 Oct 2016 Abstract. Electrical machines commonly consist of moving
More informationTechnische Universität Graz
Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2 Technische
More informationADI iterations for. general elliptic problems. John Strain Mathematics Department UC Berkeley July 2013
ADI iterations for general elliptic problems John Strain Mathematics Department UC Berkeley July 2013 1 OVERVIEW Classical alternating direction implicit (ADI) iteration Essentially optimal in simple domains
More informationOn the choice of abstract projection vectors for second level preconditioners
On the choice of abstract projection vectors for second level preconditioners C. Vuik 1, J.M. Tang 1, and R. Nabben 2 1 Delft University of Technology 2 Technische Universität Berlin Institut für Mathematik
More informationIsogEometric Tearing and Interconnecting
IsogEometric Tearing and Interconnecting Christoph Hofer and Ludwig Mitter Johannes Kepler University, Linz 26.01.2017 Doctoral Program Computational Mathematics Numerical Analysis and Symbolic Computation
More informationNumerical Analysis of Electromagnetic Fields
Pei-bai Zhou Numerical Analysis of Electromagnetic Fields With 157 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Hong Kong Barcelona Budapest Contents Part 1 Universal Concepts
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 informationMultigrid and Iterative Strategies for Optimal Control Problems
Multigrid and Iterative Strategies for Optimal Control Problems John Pearson 1, Stefan Takacs 1 1 Mathematical Institute, 24 29 St. Giles, Oxford, OX1 3LB e-mail: john.pearson@worc.ox.ac.uk, takacs@maths.ox.ac.uk
More informationOverlapping Schwarz Preconditioners for Spectral. Problem in H(curl)
Overlapping Schwarz Preconditioners for Spectral Nédélec Elements for a Model Problem in H(curl) Technical Report TR2002-83 November 22, 2002 Department of Computer Science Courant Institute of Mathematical
More informationScientific Computing
Lecture on Scientific Computing Dr. Kersten Schmidt Lecture 4 Technische Universität Berlin Institut für Mathematik Wintersemester 2014/2015 Syllabus Linear Regression Fast Fourier transform Modelling
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
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 informationA Space-Time Boundary Element Method for the Wave Equation
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 Space-Time Boundary Element Method for the Wave Equation Marco Zank and Olaf Steinbach Institut für Numerische Mathematik Space-Time Methods for PDEs,
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 informationSOLVING ELLIPTIC PDES
university-logo SOLVING ELLIPTIC PDES School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 POISSON S EQUATION Equation and Boundary Conditions Solving the Model Problem 3 THE LINEAR ALGEBRA PROBLEM
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 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 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 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 informationA Non-Overlapping Domain Decomposition Method for the Exterior Helmholtz Problem
Contemporary Mathematics Volume 28, 998 B 0-828-0988--0300-6 A Non-Overlapping Domain Decomposition Method for the Exterior Helmholtz Problem Armel de La Bourdonnaye, Charbel Farhat, Antonini Macedo, Frédéric
More informationA coupled BEM and FEM for the interior transmission problem
A coupled BEM and FEM for the interior transmission problem George C. Hsiao, Liwei Xu, Fengshan Liu, Jiguang Sun Abstract The interior transmission problem (ITP) is a boundary value problem arising in
More informationA FETI-DP method for the parallel iterative solution of indefinite and complex-valued solid and shell vibration problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Published online in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/nme.1282 A FETI-DP method for the parallel iterative solution
More informationLecture on: Numerical sparse linear algebra and interpolation spaces. June 3, 2014
Lecture on: Numerical sparse linear algebra and interpolation spaces June 3, 2014 Finite dimensional Hilbert spaces and IR N 2 / 38 (, ) : H H IR scalar product and u H = (u, u) u H norm. Finite dimensional
More informationTechnische Universität Graz
Technische Universität Graz Boundary element methods for magnetostatic field problems: A critical view Z. Andjelic, G. Of, O. Steinbach, P. Urthaler Berichte aus dem Institut für Numerische Mathematik
More informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
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 informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 31, pp. 384-402, 2008. Copyright 2008,. ISSN 1068-9613. ETNA ON THE EQUIVALENCE OF PRIMAL AND DUAL SUBSTRUCTURING PRECONDITIONERS BEDŘICH SOUSEDÍK
More 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 informationLecture 18 Classical Iterative Methods
Lecture 18 Classical Iterative Methods MIT 18.335J / 6.337J Introduction to Numerical Methods Per-Olof Persson November 14, 2006 1 Iterative Methods for Linear Systems Direct methods for solving Ax = b,
More informationUncertainty analysis of large-scale systems using domain decomposition
Center for Turbulence Research Annual Research Briefs 2007 143 Uncertainty analysis of large-scale systems using domain decomposition By D. Ghosh, C. Farhat AND P. Avery 1. Motivation and objectives A
More informationWeighted Regularization of Maxwell Equations Computations in Curvilinear Polygons
Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial IRMAR, Université de Rennes, Campus de Beaulieu, Rennes,
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 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 informationPreprint Alexander Heinlein, Axel Klawonn, and Oliver Rheinbach Parallel Two-Level Overlapping Schwarz Methods in Fluid-Structure Interaction
Fakultät für Mathematik und Informatik Preprint 2015-15 Alexander Heinlein, Axel Klawonn, and Oliver Rheinbach Parallel Two-Level Overlapping Schwarz Methods in Fluid-Structure Interaction ISSN 1433-9307
More informationConstraint interface preconditioning for the incompressible Stokes equations Loghin, Daniel
Constraint interface preconditioning for the incompressible Stokes equations Loghin, Daniel DOI: 10.1137/16M1085437 License: Other (please specify with Rights Statement) Document Version Publisher's PDF,
More informationSpace-Time Domain Decomposition Methods for Transport Problems in Porous Media
1 / 49 Space-Time Domain Decomposition Methods for Transport Problems in Porous Media Thi-Thao-Phuong Hoang, Elyes Ahmed, Jérôme Jaffré, Caroline Japhet, Michel Kern, Jean Roberts INRIA Paris-Rocquencourt
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 informationIntroduction to Domain Decomposition Methods. Alfio Quarteroni
. Aachen, EU Regional School Series February 18, 2013 Introduction to Domain Decomposition Methods Alfio Quarteroni Chair of Modeling and Scientific Computing Mathematics Institute of Computational Science
More information