[논문 리뷰] Analysis of parallel Schwarz algorithms for time-harmonic problems using block Toeplitz matrices. ETNA - Electronic Transactions on Numerical Analysis

Alexander Kyriakis, Jennifer Pestana외 2명|arXiv (Cornell University)|2021. 11. 09.

Matrix Theory and Algorithms참고 문헌 34인용 수 5

한 줄 요약

이 논문은 시간 조화파 Helmholtz 및 Maxwell 방정식에 대해 일차원 및 이차원 영역 분할 구조에서 Robin 전달 조건을 갖는 일수준 병렬 Schwarz 방법의 수렴성과 확장성에 대해 블록 토플리츠 행렬을 사용하여 분석한다. 흡수 조건(σ > 0) 하에서, 연속 수준에서 스트립형 영역 분할에 대해 방법이 약한 확장성을 확보하며, 하위영역 수 N이 클 경우 수렴 인자가 N에 독립적이게 된다. 이는 1차원 및 2차원 문제 모두에 대해 성립한다.

ABSTRACT

In this work we study the convergence properties of the one-level parallel Schwarz method with Robin transmission conditions applied to the one-dimensional and two-dimensional Helmholtz and Maxwell's equations. One-level methods are not scalable in general. However, it has recently been proven that when impedance transmission conditions are used in the case of the algorithm being applied to the equations with absorption, then, under certain assumptions, scalability can be achieved and no coarse space is required. We show here that this result is also true for the iterative version of the method at the continuous level for strip-wise decompositions into subdomains that are typically encountered when solving wave-guide problems. The convergence proof relies on the particular block Toeplitz structure of the global iteration matrix. Although non-Hermitian, we prove that its limiting spectrum has a near identical form to that of a Hermitian matrix of the same structure. We illustrate our results with numerical experiments.

연구 동기 및 목표

• 시간 조화파 문제에 대해 Robin 전달 조건을 갖는 일수준 병렬 Schwarz 방법의 확장성 연구.
• 파동 전파 맥락에서 Schwarz 방법에 대한 이론적 분석 부족 문제 해결, 특히 연속 수준에서의 분석.
• 1차원 및 2차원 Helmholtz 및 Maxwell 방정식에서 스트립형 분할에 대해 약한 확장성(하위영역 수 N에 독립적인 수렴 속도) 확립.
• 흡수 조건이 존재할 경우, 코arse 스페이스 없이도 확장성이 유지되는 조건을 규명.
• 반복적 Schwarz 스킴에서 발생하는 비에르미트 블록 토플리츠 행렬에 대해 스펙트럼 분석 기법을 확장.

제안 방법

• 연속 수준의 수식을 활용하여 Schwarz 반복을 블록 토플리츠 행렬이 경계 데이터에 작용하는 것으로 공식화.
• N → ∞의 극한에서 전체 반복 행렬의 스펙트럼을 분석하여, 비에르미트성에도 불구하고 거의 에르미트 행렬과 유사한 구조를 갖는다는 것을 증명.
• 푸리에 분석 및 스펙트럼 이론을 활용하여 다수의 하위영역에 대한 극한에서의 수렴 행동을 특성화.
• 문헌 [6]의 반복 행렬 형식을 파동 방정식과 흡수 조건을 고려한 스트립형 분할에 적용.
• 모든 σ > 0 에 대해 1차원 문제에서 수렴을 증명하고, 2차원 문제에서는 σ, δ, L에 대한 조건 하에서 수렴을 보장.
• 이론적 결과를 실용적 조건화 방법과 연결하여, 일수준 최적화된 제한된 덧셈 Schwarz(ORAS) 방법과의 연관성 규명.

실험 결과

연구 질문

• RQ1 Robin 조건을 갖는 일수준 Schwarz 방법이 코arse 스페이스 없이 시간 조화파 문제에 대해 약한 확장성을 달성할 수 있는가?
• RQ2 Helmholtz 및 Maxwell 방정식에서 스트립형 분할에 대해 Schwarz 반복 행렬의 극한 스펙트럼 행동은 어떻게 되는가?
• RQ3 흡수 조건(σ > 0)은 연속 수준의 공식화에서 수렴성과 확장성에 어떤 영향을 미치는가?
• RQ4 어떤 조건에서 방법이 k-로버스트성이 되며, 이는 영역 분할 파rameter 선택과 어떻게 관련되는가?
• RQ5 비에르미트 블록 토플리츠 행렬의 스펙트럼 분석을 통해 수렴성 및 확장성 결과를 도출할 수 있는가?

주요 결과

• 비에르미트 행렬임에도 불구하고, Schwarz 반복 행렬의 극한 스펙트럼은 동일한 블록 토플리츠 형식의 에르미트 행렬과 거의 동일한 구조를 갖는다.
• 1차원에서는 모든 σ > 0 에 대해 수렴하며, 충분히 큰 N에 대해 수렴 인자가 N에 독립적이게 되어 약한 확장성이 입증된다.
• 2차원에서는 감쇠 모드(˜k > k) 또는 σ, δ, L가 충분히 클 경우 확장성이 유지되며, σ ≥ k 일 때 수렴이 보장된다.
• 수치 실험 결과, GMRES 반복 횟수가 특정 N 이상에서 안정화되며, N과 파동 수 k에 독립적인 것으로 확인된다. (σ = 1일 때)
• 흡수 파rameter σ가 k에 비례할 경우, 방법은 k-로버스트성을 확보하며, 이는 이전 문헌의 결과와 일치한다.
• 이론적 프레임워크는 ORAS 방법과 같은 강건한 조건화 방법 설계를 지원하며, 연속 및 이산 수준 모두에서 분석 가능하게 한다.