Convergence of parallel overlapping domain decomposition methods for the Helmholtz equation

We analyse parallel overlapping Schwarz domain decomposition methods for the Helmholtz equation, where the exchange of information between subdomains is achieved using first-order absorbing (impedance) transmission conditions, together with a partition of unity. We provide a novel analysis of this m...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:Numerische Mathematik 2022-10, Vol.152 (2), p.259-306
Hauptverfasser: Gong, Shihua, Gander, Martin J., Graham, Ivan G., Lafontaine, David, Spence, Euan A.
Format: Artikel
Sprache:eng
Schlagworte:
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
container_end_page 306
container_issue 2
container_start_page 259
container_title Numerische Mathematik
container_volume 152
creator Gong, Shihua
Gander, Martin J.
Graham, Ivan G.
Lafontaine, David
Spence, Euan A.
description We analyse parallel overlapping Schwarz domain decomposition methods for the Helmholtz equation, where the exchange of information between subdomains is achieved using first-order absorbing (impedance) transmission conditions, together with a partition of unity. We provide a novel analysis of this method at the PDE level (without discretization). First, we formulate the method as a fixed point iteration, and show (in dimensions 1, 2, 3) that it is well-defined in a tensor product of appropriate local function spaces, each with L 2 impedance boundary data. We then obtain a bound on the norm of the fixed point operator in terms of the local norms of certain impedance-to-impedance maps arising from local interactions between subdomains. These bounds provide conditions under which (some power of) the fixed point operator is a contraction. In 2-d, for rectangular domains and strip-wise domain decompositions (with each subdomain only overlapping its immediate neighbours), we present two techniques for verifying the assumptions on the impedance-to-impedance maps that ensure power contractivity of the fixed point operator. The first is through semiclassical analysis, which gives rigorous estimates valid as the frequency tends to infinity. At least for a model case with two subdomains, these results verify the required assumptions for sufficiently large overlap. For more realistic domain decompositions, we directly compute the norms of the impedance-to-impedance maps by solving certain canonical (local) eigenvalue problems. We give numerical experiments that illustrate the theory. These also show that the iterative method remains convergent and/or provides a good preconditioner in cases not covered by the theory, including for general domain decompositions, such as those obtained via automatic graph-partitioning software.
doi_str_mv 10.1007/s00211-022-01318-8
format Article
fullrecord <record><control><sourceid>proquest_hal_p</sourceid><recordid>TN_cdi_hal_primary_oai_HAL_hal_04284755v1</recordid><sourceformat>XML</sourceformat><sourcesystem>PC</sourcesystem><sourcerecordid>2717709873</sourcerecordid><originalsourceid>FETCH-LOGICAL-c397t-91cfbf0fc03c63c4be7b6a32ac65230a45e248f82ecd170513b8b3a3065c83a33</originalsourceid><addsrcrecordid>eNp9kEFLxDAUhIsouK7-AU8BTx6qL0nTpMdlUVdY8KIgeAhpmux2SZtu0hX019ta0ZuneQzfDI9JkksMNxiA30YAgnEKhKSAKRapOEpmUGQspSRjx8MNpEhZUbyeJmcx7gAwzzM8S96Wvn03YWNabZC3qFNBOWcc8oPrVNfV7QZVvlF1iyqjfdP5WPe1b1Fj-q2vIrI-oH5r0Mq4Zutd_4nM_qBG5Dw5scpFc_Gj8-Tl_u55uUrXTw-Py8U61bTgfVpgbUsLVgPVOdVZaXiZK0qUzhmhoDJmSCasIEZXmAPDtBQlVRRypsWgdJ5cT71b5WQX6kaFD-lVLVeLtRw9yIjIOGPveGCvJrYLfn8wsZc7fwjt8J4kHHMOheBjI5koHXyMwdjfWgxyHFxOg8thcPk9uBRDiE6hOMDtxoS_6n9SXy-8g7g</addsrcrecordid><sourcetype>Open Access Repository</sourcetype><iscdi>true</iscdi><recordtype>article</recordtype><pqid>2717709873</pqid></control><display><type>article</type><title>Convergence of parallel overlapping domain decomposition methods for the Helmholtz equation</title><source>SpringerNature Journals</source><creator>Gong, Shihua ; Gander, Martin J. ; Graham, Ivan G. ; Lafontaine, David ; Spence, Euan A.</creator><creatorcontrib>Gong, Shihua ; Gander, Martin J. ; Graham, Ivan G. ; Lafontaine, David ; Spence, Euan A.</creatorcontrib><description>We analyse parallel overlapping Schwarz domain decomposition methods for the Helmholtz equation, where the exchange of information between subdomains is achieved using first-order absorbing (impedance) transmission conditions, together with a partition of unity. We provide a novel analysis of this method at the PDE level (without discretization). First, we formulate the method as a fixed point iteration, and show (in dimensions 1, 2, 3) that it is well-defined in a tensor product of appropriate local function spaces, each with L 2 impedance boundary data. We then obtain a bound on the norm of the fixed point operator in terms of the local norms of certain impedance-to-impedance maps arising from local interactions between subdomains. These bounds provide conditions under which (some power of) the fixed point operator is a contraction. In 2-d, for rectangular domains and strip-wise domain decompositions (with each subdomain only overlapping its immediate neighbours), we present two techniques for verifying the assumptions on the impedance-to-impedance maps that ensure power contractivity of the fixed point operator. The first is through semiclassical analysis, which gives rigorous estimates valid as the frequency tends to infinity. At least for a model case with two subdomains, these results verify the required assumptions for sufficiently large overlap. For more realistic domain decompositions, we directly compute the norms of the impedance-to-impedance maps by solving certain canonical (local) eigenvalue problems. We give numerical experiments that illustrate the theory. These also show that the iterative method remains convergent and/or provides a good preconditioner in cases not covered by the theory, including for general domain decompositions, such as those obtained via automatic graph-partitioning software.</description><identifier>ISSN: 0029-599X</identifier><identifier>EISSN: 0945-3245</identifier><identifier>DOI: 10.1007/s00211-022-01318-8</identifier><language>eng</language><publisher>Berlin/Heidelberg: Springer Berlin Heidelberg</publisher><subject>Convergence ; Domain decomposition methods ; Eigenvalues ; Fixed points (mathematics) ; Function space ; Helmholtz equations ; Impedance ; Iterative methods ; Mathematical analysis ; Mathematical and Computational Engineering ; Mathematical and Computational Physics ; Mathematical Methods in Physics ; Mathematics ; Mathematics and Statistics ; Norms ; Numerical Analysis ; Numerical and Computational Physics ; Operators (mathematics) ; Simulation ; Tensors ; Theoretical</subject><ispartof>Numerische Mathematik, 2022-10, Vol.152 (2), p.259-306</ispartof><rights>The Author(s) 2022</rights><rights>The Author(s) 2022. This work is published under http://creativecommons.org/licenses/by/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.</rights><rights>Distributed under a Creative Commons Attribution 4.0 International License</rights><lds50>peer_reviewed</lds50><oa>free_for_read</oa><woscitedreferencessubscribed>false</woscitedreferencessubscribed><citedby>FETCH-LOGICAL-c397t-91cfbf0fc03c63c4be7b6a32ac65230a45e248f82ecd170513b8b3a3065c83a33</citedby><cites>FETCH-LOGICAL-c397t-91cfbf0fc03c63c4be7b6a32ac65230a45e248f82ecd170513b8b3a3065c83a33</cites></display><links><openurl>$$Topenurl_article</openurl><openurlfulltext>$$Topenurlfull_article</openurlfulltext><thumbnail>$$Tsyndetics_thumb_exl</thumbnail><linktopdf>$$Uhttps://link.springer.com/content/pdf/10.1007/s00211-022-01318-8$$EPDF$$P50$$Gspringer$$Hfree_for_read</linktopdf><linktohtml>$$Uhttps://link.springer.com/10.1007/s00211-022-01318-8$$EHTML$$P50$$Gspringer$$Hfree_for_read</linktohtml><link.rule.ids>230,314,780,784,885,27924,27925,41488,42557,51319</link.rule.ids><backlink>$$Uhttps://hal.science/hal-04284755$$DView record in HAL$$Hfree_for_read</backlink></links><search><creatorcontrib>Gong, Shihua</creatorcontrib><creatorcontrib>Gander, Martin J.</creatorcontrib><creatorcontrib>Graham, Ivan G.</creatorcontrib><creatorcontrib>Lafontaine, David</creatorcontrib><creatorcontrib>Spence, Euan A.</creatorcontrib><title>Convergence of parallel overlapping domain decomposition methods for the Helmholtz equation</title><title>Numerische Mathematik</title><addtitle>Numer. Math</addtitle><description>We analyse parallel overlapping Schwarz domain decomposition methods for the Helmholtz equation, where the exchange of information between subdomains is achieved using first-order absorbing (impedance) transmission conditions, together with a partition of unity. We provide a novel analysis of this method at the PDE level (without discretization). First, we formulate the method as a fixed point iteration, and show (in dimensions 1, 2, 3) that it is well-defined in a tensor product of appropriate local function spaces, each with L 2 impedance boundary data. We then obtain a bound on the norm of the fixed point operator in terms of the local norms of certain impedance-to-impedance maps arising from local interactions between subdomains. These bounds provide conditions under which (some power of) the fixed point operator is a contraction. In 2-d, for rectangular domains and strip-wise domain decompositions (with each subdomain only overlapping its immediate neighbours), we present two techniques for verifying the assumptions on the impedance-to-impedance maps that ensure power contractivity of the fixed point operator. The first is through semiclassical analysis, which gives rigorous estimates valid as the frequency tends to infinity. At least for a model case with two subdomains, these results verify the required assumptions for sufficiently large overlap. For more realistic domain decompositions, we directly compute the norms of the impedance-to-impedance maps by solving certain canonical (local) eigenvalue problems. We give numerical experiments that illustrate the theory. These also show that the iterative method remains convergent and/or provides a good preconditioner in cases not covered by the theory, including for general domain decompositions, such as those obtained via automatic graph-partitioning software.</description><subject>Convergence</subject><subject>Domain decomposition methods</subject><subject>Eigenvalues</subject><subject>Fixed points (mathematics)</subject><subject>Function space</subject><subject>Helmholtz equations</subject><subject>Impedance</subject><subject>Iterative methods</subject><subject>Mathematical analysis</subject><subject>Mathematical and Computational Engineering</subject><subject>Mathematical and Computational Physics</subject><subject>Mathematical Methods in Physics</subject><subject>Mathematics</subject><subject>Mathematics and Statistics</subject><subject>Norms</subject><subject>Numerical Analysis</subject><subject>Numerical and Computational Physics</subject><subject>Operators (mathematics)</subject><subject>Simulation</subject><subject>Tensors</subject><subject>Theoretical</subject><issn>0029-599X</issn><issn>0945-3245</issn><fulltext>true</fulltext><rsrctype>article</rsrctype><creationdate>2022</creationdate><recordtype>article</recordtype><sourceid>C6C</sourceid><recordid>eNp9kEFLxDAUhIsouK7-AU8BTx6qL0nTpMdlUVdY8KIgeAhpmux2SZtu0hX019ta0ZuneQzfDI9JkksMNxiA30YAgnEKhKSAKRapOEpmUGQspSRjx8MNpEhZUbyeJmcx7gAwzzM8S96Wvn03YWNabZC3qFNBOWcc8oPrVNfV7QZVvlF1iyqjfdP5WPe1b1Fj-q2vIrI-oH5r0Mq4Zutd_4nM_qBG5Dw5scpFc_Gj8-Tl_u55uUrXTw-Py8U61bTgfVpgbUsLVgPVOdVZaXiZK0qUzhmhoDJmSCasIEZXmAPDtBQlVRRypsWgdJ5cT71b5WQX6kaFD-lVLVeLtRw9yIjIOGPveGCvJrYLfn8wsZc7fwjt8J4kHHMOheBjI5koHXyMwdjfWgxyHFxOg8thcPk9uBRDiE6hOMDtxoS_6n9SXy-8g7g</recordid><startdate>20221001</startdate><enddate>20221001</enddate><creator>Gong, Shihua</creator><creator>Gander, Martin J.</creator><creator>Graham, Ivan G.</creator><creator>Lafontaine, David</creator><creator>Spence, Euan A.</creator><general>Springer Berlin Heidelberg</general><general>Springer Nature B.V</general><general>Springer Verlag</general><scope>C6C</scope><scope>AAYXX</scope><scope>CITATION</scope><scope>1XC</scope></search><sort><creationdate>20221001</creationdate><title>Convergence of parallel overlapping domain decomposition methods for the Helmholtz equation</title><author>Gong, Shihua ; Gander, Martin J. ; Graham, Ivan G. ; Lafontaine, David ; Spence, Euan A.</author></sort><facets><frbrtype>5</frbrtype><frbrgroupid>cdi_FETCH-LOGICAL-c397t-91cfbf0fc03c63c4be7b6a32ac65230a45e248f82ecd170513b8b3a3065c83a33</frbrgroupid><rsrctype>articles</rsrctype><prefilter>articles</prefilter><language>eng</language><creationdate>2022</creationdate><topic>Convergence</topic><topic>Domain decomposition methods</topic><topic>Eigenvalues</topic><topic>Fixed points (mathematics)</topic><topic>Function space</topic><topic>Helmholtz equations</topic><topic>Impedance</topic><topic>Iterative methods</topic><topic>Mathematical analysis</topic><topic>Mathematical and Computational Engineering</topic><topic>Mathematical and Computational Physics</topic><topic>Mathematical Methods in Physics</topic><topic>Mathematics</topic><topic>Mathematics and Statistics</topic><topic>Norms</topic><topic>Numerical Analysis</topic><topic>Numerical and Computational Physics</topic><topic>Operators (mathematics)</topic><topic>Simulation</topic><topic>Tensors</topic><topic>Theoretical</topic><toplevel>peer_reviewed</toplevel><toplevel>online_resources</toplevel><creatorcontrib>Gong, Shihua</creatorcontrib><creatorcontrib>Gander, Martin J.</creatorcontrib><creatorcontrib>Graham, Ivan G.</creatorcontrib><creatorcontrib>Lafontaine, David</creatorcontrib><creatorcontrib>Spence, Euan A.</creatorcontrib><collection>Springer Nature OA Free Journals</collection><collection>CrossRef</collection><collection>Hyper Article en Ligne (HAL)</collection><jtitle>Numerische Mathematik</jtitle></facets><delivery><delcategory>Remote Search Resource</delcategory><fulltext>fulltext</fulltext></delivery><addata><au>Gong, Shihua</au><au>Gander, Martin J.</au><au>Graham, Ivan G.</au><au>Lafontaine, David</au><au>Spence, Euan A.</au><format>journal</format><genre>article</genre><ristype>JOUR</ristype><atitle>Convergence of parallel overlapping domain decomposition methods for the Helmholtz equation</atitle><jtitle>Numerische Mathematik</jtitle><stitle>Numer. Math</stitle><date>2022-10-01</date><risdate>2022</risdate><volume>152</volume><issue>2</issue><spage>259</spage><epage>306</epage><pages>259-306</pages><issn>0029-599X</issn><eissn>0945-3245</eissn><abstract>We analyse parallel overlapping Schwarz domain decomposition methods for the Helmholtz equation, where the exchange of information between subdomains is achieved using first-order absorbing (impedance) transmission conditions, together with a partition of unity. We provide a novel analysis of this method at the PDE level (without discretization). First, we formulate the method as a fixed point iteration, and show (in dimensions 1, 2, 3) that it is well-defined in a tensor product of appropriate local function spaces, each with L 2 impedance boundary data. We then obtain a bound on the norm of the fixed point operator in terms of the local norms of certain impedance-to-impedance maps arising from local interactions between subdomains. These bounds provide conditions under which (some power of) the fixed point operator is a contraction. In 2-d, for rectangular domains and strip-wise domain decompositions (with each subdomain only overlapping its immediate neighbours), we present two techniques for verifying the assumptions on the impedance-to-impedance maps that ensure power contractivity of the fixed point operator. The first is through semiclassical analysis, which gives rigorous estimates valid as the frequency tends to infinity. At least for a model case with two subdomains, these results verify the required assumptions for sufficiently large overlap. For more realistic domain decompositions, we directly compute the norms of the impedance-to-impedance maps by solving certain canonical (local) eigenvalue problems. We give numerical experiments that illustrate the theory. These also show that the iterative method remains convergent and/or provides a good preconditioner in cases not covered by the theory, including for general domain decompositions, such as those obtained via automatic graph-partitioning software.</abstract><cop>Berlin/Heidelberg</cop><pub>Springer Berlin Heidelberg</pub><doi>10.1007/s00211-022-01318-8</doi><tpages>48</tpages><oa>free_for_read</oa></addata></record>
fulltext fulltext
identifier ISSN: 0029-599X
ispartof Numerische Mathematik, 2022-10, Vol.152 (2), p.259-306
issn 0029-599X
0945-3245
language eng
recordid cdi_hal_primary_oai_HAL_hal_04284755v1
source SpringerNature Journals
subjects Convergence
Domain decomposition methods
Eigenvalues
Fixed points (mathematics)
Function space
Helmholtz equations
Impedance
Iterative methods
Mathematical analysis
Mathematical and Computational Engineering
Mathematical and Computational Physics
Mathematical Methods in Physics
Mathematics
Mathematics and Statistics
Norms
Numerical Analysis
Numerical and Computational Physics
Operators (mathematics)
Simulation
Tensors
Theoretical
title Convergence of parallel overlapping domain decomposition methods for the Helmholtz equation
url https://sfx.bib-bvb.de/sfx_tum?ctx_ver=Z39.88-2004&ctx_enc=info:ofi/enc:UTF-8&ctx_tim=2024-12-25T00%3A23%3A27IST&url_ver=Z39.88-2004&url_ctx_fmt=infofi/fmt:kev:mtx:ctx&rfr_id=info:sid/primo.exlibrisgroup.com:primo3-Article-proquest_hal_p&rft_val_fmt=info:ofi/fmt:kev:mtx:journal&rft.genre=article&rft.atitle=Convergence%20of%20parallel%20overlapping%20domain%20decomposition%20methods%20for%20the%20Helmholtz%20equation&rft.jtitle=Numerische%20Mathematik&rft.au=Gong,%20Shihua&rft.date=2022-10-01&rft.volume=152&rft.issue=2&rft.spage=259&rft.epage=306&rft.pages=259-306&rft.issn=0029-599X&rft.eissn=0945-3245&rft_id=info:doi/10.1007/s00211-022-01318-8&rft_dat=%3Cproquest_hal_p%3E2717709873%3C/proquest_hal_p%3E%3Curl%3E%3C/url%3E&disable_directlink=true&sfx.directlink=off&sfx.report_link=0&rft_id=info:oai/&rft_pqid=2717709873&rft_id=info:pmid/&rfr_iscdi=true