[논문 리뷰] Natural domain decomposition algorithms for the solution of time-harmonic elastic waves

Romain Brunet, Victorita Dolean외 1명|arXiv (Cornell University)|2019. 04. 27.

Advanced Numerical Methods in Computational Mathematics참고 문헌 27인용 수 7

한 줄 요약

이 논문은 나비에 방정식에 의해 지배되는 시간 조화 탄성파 문제를 해결하기 위한 최적화된 슈바르츠 도메인 분할 방법을 도입하고 분석한다. 고전적 슈바르츠 방법이 이러한 방정식에서는 발산함을 증명하지만, 물리적으로 유용한 전달 조건—특히 국소적 저주파수 근사—를 사용할 경우 수렴하고, 켈러리 해법기의 전처리로써 뛰어난 성능을 발휘함을 보여주며, 이는 이중 및 다중 하위영역 문제, 특히 비균질 전달 문제에서 고전적 방법보다 크게 뛰어난 성능을 발휘한다.

ABSTRACT

We study for the first time Schwarz domain decomposition methods for the solution of the Navier equations modeling the propagation of elastic waves. These equations in the time harmonic regime are difficult to solve by iterative methods, even more so than the Helmholtz equation. We first prove that the classical Schwarz method is not convergent when applied to the Navier equations, and can thus not be used as an iterative solver, only as a preconditioner for a Krylov method. We then introduce more natural transmission conditions between the subdomains, and show that if the overlap is not too small, this new Schwarz method is convergent. We illustrate our results with numerical experiments, both for situations covered by our technical two subdomain analysis, and situations that go far beyond, including many subdomains, cross points, heterogeneous materials in a transmission problem, and Krylov acceleration. Our numerical results show that the Schwarz method with adapted transmission conditions leads systematically to a better solver for the Navier equations than the classical Schwarz method.

연구 동기 및 목표

• 시간 조화 탄성파 문제에 대해 나비에 방정식으로 모델링된 고전적 슈바르츠 방법의 수렴 행동을 조사하는 것.
• 탄성파 문제에 대해 효과적인 도메인 분할 해법기가 부족한 문제를 해결하는 것—이는 벡터 성격과 다수의 파동 유형으로 인해 헬름홀츠 문제보다 더 어렵기 때문이다.
• 나비에 방정식에 대해 슈바르츠 방법의 수렴을 보장하는 새로운 전달 조건을 개발하고 분석하는 것.
• 실용적 구현이 가능하면서도 수렴성을 유지하는 최적 전달 조건의 국소적 저주파수 근사를 설계하는 것.
• 복잡한 설정, 즉 많은 하위영역과 비균질 물질을 포함한 경우에 새로운 방법이 해법기 및 전처리기로써 고전적 방법보다 뛰어난 성능을 보임을 수치적으로 검증하는 것.

제안 방법

• 시간 조화 영역에서 두 겹침이 있는 하위영역에 대해 연속 수준에서 고전적 슈바르츠 알고리즘을 분석하여, 나비에 방정식에 대해 수렴하지 않음을 증명한다.
• 더 자연스럽고 물리적으로 유용한 하위영역 간 전달 조건을 도입하여, 이 조건이 유한한 단계 내에서 수렴하도록 할 수 있음을 보여준다.
• 최적 전달 조건의 국소적 저주파수 근사를 제안하여 비국소 연산자를 피하고 실용적 구현을 가능하게 한다.
• 하위영역 간 겹침이 너무 작지 않다면, 새로운 최적화된 슈바르츠 방법의 수렴성을 증명한다.
• 병렬 계산이 가능한 이산 문제에 대해 제한된 덧셈 슈바르츠(RAS) 및 최적화된 RAS(ORAS)를 사용하여 방법을 구현한다.
• 성능 평가를 위해 고전적 및 최적화된 슈바르츠 방법을 전처리기로 사용하는 GMRES를 켈러리 해법기로 활용한다.

실험 결과

연구 질문

• RQ1 왜 고전적 슈바르츠 방법은 시간 조화 탄성파에 대해 나비에 방정식에 적용될 때 수렴하지 못하는가?
• RQ2 P-파와 S-파의 물리적 거동을 반영하는 전달 조건이 나비에 방정식에 대해 수렴하는 슈바르츠 방법을 가능하게 할 수 있는가?
• RQ3 최적 전달 조건의 국소적 저주파수 근사를 구성할 수 있으며, 이는 나비에 방정식에 대해 수렴성을 유지할 수 있는가?
• RQ4 복잡한 다중 하위영역 및 비균질 문제에서 최적화된 슈바르츠 방법이 해법기 및 전처리기로써 고전적 방법과 비교해 어떻게 성능이 뛰어나게 되는가?
• RQ5 하위영역 간 겹침을 증가시키면 최적화된 슈바르츠 방법의 수렴성이 향상되는가? 만약 그렇다면 어떤 조건에서인가?

주요 결과

• 고전적 슈바르츠 방법은 겹침이 있더라도 나비에 방정식에 적용될 경우 발산하며, 반복적 해법기로 사용될 수 없다.
• 비국소 연산자를 포함하는 최적 전달 조건은 두 번의 반복 내에서 슈바르츠 방법의 수렴을 가능하게 하지만, 실용적 구현에는 부적합하다.
• 최적 전달 조건의 국소적 저주파수 근사는 겹침이 충분히 크다면 최적화된 슈바르츠 방법의 수렴성을 보장한다.
• 수치 실험 결과, 최적화된 슈바르츠 방법은 해법기로써 수렴하는 반면, 고전적 방법은 겹침이 증가함에도 불구하고 발산함을 확인하였다.
• GMRES의 전처리기로 사용할 경우, 최적화된 슈바르츠 방법은 특히 많은 하위영역, 교차점, 물성 비균질성이 있는 문제에서 고전적 방법보다 훨씬 적은 반복 수를 요구한다.
• 최적화된 슈바르츠 방법의 수렴성은 겹침이 증가함에 따라 향상되며, 반면 고전적 방법의 성능은 동일 조건에서 악화된다.