[논문 리뷰] Domain Decomposition Methods based on quasi-optimal transmission operators for the solution of Helmholtz transmission problems

Yassine Boubendir, Carlos Jerez-Hanckes외 2명|arXiv (Cornell University)|2017. 10. 07.

Electromagnetic Scattering and Analysis참고 문헌 22인용 수 3

한 줄 요약

이 논문은 딜리클레-투-노이만(DtN) 연산자의 근사치를 이용해 유도된 준최적 전달 연산자를 사용하는 헬름홀츠 전송 문제에 대한 겹치지 않는 도메인 분할 방법(DDM)을 제시한다. 복소수로 확장된 초특이 경계 적분 연산자 또는 제곱근 푸리에 다중연산자를 활용함으로써, 높은 대비, 고주파 영역에서도 거의 주파수, 대비 또는 하위 도메인 수에 민감하지 않게 몇 번의 반복만으로도 강건한 수렴성을 확보한다.

ABSTRACT

We present non-overlapping Domain Decomposition Methods (DDM) based on quasi-optimal transmission operators for the solution of Helmholtz transmission problems with piece-wise constant material properties. The quasi-optimal transmission boundary conditions incorporate readily available approximations of Dirichlet-to-Neumann operators. These approximations consist of either complexified hypersingular boundary integral operators for the Helmholtz equation or square root Fourier multipliers with complex wavenumbers. We show that under certain regularity assumptions on the closed interface of material discontinuity, the DDM with quasi-optimal transmission conditions are well-posed. We present a DDM framework based on Robin-to-Robin (RtR) operators that can be computed robustly via boundary integral formulations. More importantly, the use of quasi-optimal transmission operators results in DDM that converge in small numbers of iterations even in the challenging high-contrast, high-frequency regime of Helmholtz transmission problems. Furthermore, the DDM presented in this text require only minor modifications to handle the case of transmission problems in partially coated domains, while still maintaining excellent convergence properties. We also investigate the dependence of the DDM iterative performance on the number of subdomains.

연구 동기 및 목표

• 높은 대비, 고주파 영역에서 전통적인 도메인 분할 방법(DDM)의 느린 수렴 문제를 해결한다.
• 어려운 매개변수 영역에서 많은 케일리-크리스털 반복이 필요한 고전적 로빈 경계 조건과 표준 경계 적분 방정식(BIE) 설정의 한계를 극복한다.
• 부분적으로 코팅된 도메인을 포함한 복잡한 물질 구성에서도 수렴성이 강건하게 유지되는 DDM 프레임워크를 개발한다.
• 전달 연산자에 파동수를 복소수로 복소화하여 하위 도메인 문제의 잘 정의됨을 보장한다.
• 다양한 수의 하위 도메인과 고주파 매개변수에 걸쳐 반복 횟수를 최소화하여 확장 가능하고 효율적인 해법 전략을 확보한다.

제안 방법

• 디릴리클레-투-노이만(DtN) 연산자의 근사치를 기반으로 한 준최적 전달 연산자를 사용하는 겹치지 않는 DDM을 수립한다.
• 복소수 파동수를 갖는 복소수로 확장된 초특이 경계 적분 연산자(BIOs) 또는 제곱근 푸리에 다중연산자를 전달 연산자로 사용한다.
• 전달 연산자에 파동수를 복소수로 복소화하여 코ercivity 조건을 만족시킴으로써 하위 도메인 헬름홀츠 문제의 잘 정의됨을 보장한다.
• 경계 적분 수식을 통해 계산된 로빈-투-로빈(RtR) 연산자 프레임워크를 구현하여 강건하고 효율적인 반복를 확보한다.
• 전달 연산자를 수정함으로써 부분적으로 코팅된 도메인에 이 방법을 적응시키되, 수렴 성능을 저하시키지 않는다.
• 인터페이스 문제를 해결하기 위해 케일리-크리스털 부분공간 해법기를 사용하며, 반복 횟수와 원거리장 오차를 통해 수렴을 모니터링한다.

실험 결과

연구 질문

• RQ1 DtN 근사치를 기반으로 한 준최적 전달 연산자가 헬름홀츠 전송 문제의 DDM에서 반복 횟수를 크게 줄일 수 있는가?
• RQ2 준최적 전달 연산자를 갖는 DDM의 수렴 특성은 주파수와 물질 대비가 증가함에 따라 어떻게 변화하는가?
• RQ3 완전 전도체(PEC) 경계를 갖는 부분적으로 코팅된 도메인에 적용했을 때, DDM 프레임워크는 여전히 강건하고 효율적인가?
• RQ4 제안된 DDM에서 하위 도메인 수가 반복 횟수에 미치는 영향은 어떠한가?
• RQ5 높은 대비 헬름홀츠 문제에서 표준 CFIESK 수식보다 반복 횟수와 정확도 측면에서 더 우수한 성능을 보일 수 있는가?

주요 결과

• ε₁ = 16이고 ω = 64인 원형 도메인에서, 준최적 전달 연산자 $Z_j$를 사용한 DDM은 23회의 반복만으로도 충분했고, CFIESK 수식은 2271회의 반복을 필요로 하였다.
• 부분적으로 코팅된 원형 도메인의 경우, $Z_j$를 사용한 DDM은 ω = 32일 때 반복 횟수를 CFIESK의 2271회에서 23회로 줄였으며, 원거리장 오차는 8.5×10⁻³ 이하였다.
• 도메인이 네 개의 1/4 원형 하위 도메인으로 분할된 경우, ω = 32일 때 $Z_j$를 사용한 DDM은 173회의 반복을 필요로 하였고, 주파수 증가에 따라 반복 횟수의 증가가 거의 없었다.
• 두 개의 하위 도메인에서 $Z_j^{PS}$를 사용한 DDM은 ω = 32일 때 256회의 반복을 필요로 하였고, $Z_j^a$는 529회의 반복을 요구하여, 전달 연산자의 선택이 수렴 속도에 상당한 영향을 미친다.
• 세 개의 하위 도메인으로 나누어진 L자형 도메인의 경우, ω = 32일 때 $Z_j$를 사용한 DDM은 447회의 반복을 필요로 하였으며, 하위 도메인 수 증가에 따라도 수렴성이 유지됨을 보여주었다.
• 모든 수식에서 원거리장 오차는 8.5×10⁻³ 이하로 유지되었으며, 이는 적은 반복 횟수에도 불구하고 DDM의 정확도가 CFIESK 방법과 유사하거나 더 뛰어나다는 것을 시사한다.