KMS Of Academy of mathematics and systems sciences, CAS
Trace transfer-based diagonal sweeping domain decomposition method for the Helmholtz equation: Algorithms and convergence analysis | |
Leng, Wei1; Ju, Lili2 | |
2022-04-15 | |
Source Publication | JOURNAL OF COMPUTATIONAL PHYSICS
![]() |
ISSN | 0021-9991 |
Volume | 455Pages:29 |
Abstract | By utilizing the perfectly matched layer (PML) and source transfer techniques, the diagonal sweeping domain decomposition method (DDM) was recently developed for solving the high-frequency Helmholtz equation in R-n, which uses 2(n) sweeps along respective diagonal directions with checkerboard domain decomposition. Although this diagonal sweeping DDM is essentially multiplicative, it is highly suitable for parallel computing of the Helmholtz problem with multiple right-hand sides when combined with the pipeline processing since the number of sequential steps in each sweep is much smaller than the number of subdomains. In this paper, we propose and analyze a trace transfer-based diagonal sweeping DDM. A major advantage of changing from source transfer to trace transfer for information passing between neighbor subdomains is that the resulting diagonal sweeps become easier to analyze and implement and more efficient, since the transferred traces have only 2n cardinal directions between neighbor subdomains while the transferred sources come from a total of 3(n) - 1 cardinal and corner directions. We rigorously prove that the proposed diagonal sweeping DDM not only gives the exact solution of the global PML problem in the constant medium case but also does it with at most one extra round of diagonal sweeps in the two-layered media case, which lays down the theoretical foundation of the method. Performance and parallel scalability of the proposed DDM as direct solver or preconditioner are also numerically demonstrated through extensive experiments in two and three dimensions. (C) 2022 Elsevier Inc. All rights reserved. |
Keyword | Domain decomposition method Diagonal sweeping Helmholtz equation Perfectly matched layer Trace transfer Parallel computing |
DOI | 10.1016/j.jcp.2022.110980 |
Indexed By | SCI |
Language | 英语 |
Funding Project | National Natural Science Foundation of China[12171464] ; National Natural Science Foundation of China[11771440] ; National Key R&D Program of China[2020YFA0711904] ; Strategic Priority Research Program of Chinese Academy of Sciences[XDB 41000000] ; National Center for Mathematics and Interdisciplinary Sciences of Chinese Academy of Sciences (NCMIS) ; U.S. National Science Foundation[DMS-1818438] ; U.S. National Science Foundation[DMS-2109633] |
WOS Research Area | Computer Science ; Physics |
WOS Subject | Computer Science, Interdisciplinary Applications ; Physics, Mathematical |
WOS ID | WOS:000762463300014 |
Publisher | ACADEMIC PRESS INC ELSEVIER SCIENCE |
Citation statistics | |
Document Type | 期刊论文 |
Identifier | http://ir.amss.ac.cn/handle/2S8OKBNM/60101 |
Collection | 中国科学院数学与系统科学研究院 |
Corresponding Author | Ju, Lili |
Affiliation | 1.Chinese Acad Sci, State Key Lab Sci & Engn Comp, Beijing 100190, Peoples R China 2.Univ South Carolina, Dept Math, Columbia, SC 29208 USA |
Recommended Citation GB/T 7714 | Leng, Wei,Ju, Lili. Trace transfer-based diagonal sweeping domain decomposition method for the Helmholtz equation: Algorithms and convergence analysis[J]. JOURNAL OF COMPUTATIONAL PHYSICS,2022,455:29. |
APA | Leng, Wei,&Ju, Lili.(2022).Trace transfer-based diagonal sweeping domain decomposition method for the Helmholtz equation: Algorithms and convergence analysis.JOURNAL OF COMPUTATIONAL PHYSICS,455,29. |
MLA | Leng, Wei,et al."Trace transfer-based diagonal sweeping domain decomposition method for the Helmholtz equation: Algorithms and convergence analysis".JOURNAL OF COMPUTATIONAL PHYSICS 455(2022):29. |
Files in This Item: | There are no files associated with this item. |
Items in the repository are protected by copyright, with all rights reserved, unless otherwise indicated.
Edit Comment