[논문 리뷰] Domain Decomposition with local impedance conditions for the Helmholtz equation with absorption

Ivan G. Graham, Euan A. Spence외 1명|arXiv (Cornell University)|2018. 06. 10.

Electromagnetic Scattering and Analysis참고 문헌 52인용 수 5

한 줄 요약

이 논문은 흡수성 헬름홀츠 방정식에 대해 겹치는 부분영역과 국소적 임피던스 조건을 사용하고, 유일성 기반의 확장 연산을 적용한 일차 수준의 애드티브 샤워츠 조절자(Preconditioner)를 제안한다. 수치 실험에서는 k-독립적인 GMRES 수렴을 달성하였고, 이는 ε/k 가 일정하거나 천천히 증가하는 흡수 문제에 대해 안정성을 보여주는 이론적 근거를 제공한다. 이는 순수 헬름홀츠 문제에 대한 조절자로서의 효과성을 뒷받침한다.

ABSTRACT

We consider one-level additive Schwarz preconditioners for a family of Helmholtz problems with absorption and increasing wavenumber $k$. These problems are discretized using the Galerkin method with nodal conforming finite elements of any (fixed) order on meshes with diameter $h = h(k)$, chosen to maintain accuracy as $k$ increases. The action of the preconditioner requires solution of independent (parallel) subproblems (with impedance boundary conditions) on overlapping subdomains of diameter $H$ and overlap $δ\leq H$. The solutions of these subproblems are linked together using prolongation/restriction operators defined using a partition of unity. In numerical experiments (with $δ\sim H$) for a model interior impedance problem, we observe robust (i.e. $k-$independent) GMRES convergence as $k$ increases. This provides a highly-parallel, $k-$robust one-level domain decomposition method. We provide supporting theory by studying the preconditioner applied to a range of absorptive problems, $k^2\mapsto k^2+ \mathrm{i} \varepsilon$, with absorption parameter $\varepsilon$. Working in the Helmholtz ``energy'' inner product, and using the underlying theory of Helmholtz boundary-value problems, we prove a $k-$independent upper bound on the norm of the preconditioned matrix, valid for all $\vert \varepsilon\vert \lesssim k^2$. We also prove a strictly-positive lower bound on the distance of the field of values of the preconditioned matrix from the origin which holds when $\varepsilon/k$ is constant or growing arbitrarily slowly with $k$. These results imply robustness of the preconditioner for the corresponding absorptive problem as k increases and give theoretical support for the observed robustness of the preconditioner for the pure Helmholtz problem.

연구 동기 및 목표

• 고주파 헬름홀츠 문제에 대해 병렬적이고 k-강건한 조절자를 개발한다.
• 겹치는 부분영역에서 국소적 임피던스 조건을 사용한 일차 수준의 애드티브 샤워츠 방법의 수렴 성질을 분석한다.
• 수치 실험에서 관측된 k-강건성에 대한 이론적 근거를 제공하기 위해, 조절된 행렬의 노름과 필드 값에 대한 경계를 도출한다.
• 흡수성 헬름홀츠 문제를 순수 헬름홀츠 문제에 대한 효과적인 조절자로 사용할 수 있음을 뒷받침한다.

제안 방법

• 방법은 지름 H와 겹침 δ ≤ H를 갖는 겹치는 부분영역을 사용하며, 부분영역 경계에서 임피던스 경계 조건을 적용한다.
• 유한요소 해석을 사용하여 부분문제를 병렬로 풀며, 노드 기반의 일관된 요소와 메쉬 크기 h = h(k)를 사용하여 k 증가에 따라 정확도를 유지한다.
• 확장 및 제약 연산자는 유일성 기반의 분할을 사용하여, 애드티브 샤워츠 프레임워크 내에서 일관성과 안정성을 보장한다.
• 조절자 분석은 헬름홀츠 에너지 내적을 기반으로 하며, 흡수성 헬름홀츠 경계값 문제 이론을 활용한다.
• |ε| ≤ Ck² 이고 ε/k 가 유계이거나 천천히 증가하는 조건 하에서, 조절된 행렬의 노름과 필드 값에 대한 이론적 경계를 유도한다.
• 수치 실험은 다양한 k, H, ε 및 부분영역 크기에서 수렴성을 시험하며, 고정된 H와 딜리클레 부분영역 조건을 포함한 경우도 포함한다.

실험 결과

연구 질문

• RQ1 국소적 임피던스 조건을 사용한 도메인 분할 조절자가 고파장수 헬름홀츠 문제에서 k-독립적인 GMRES 수렴을 달성할 수 있는가?
• RQ2 흡수성 헬름홀츠 설정에서 조절된 행렬의 노름과 필드 값에 대해 어떤 이론적 경계를 설정할 수 있는가?
• RQ3 조절자에서 흡수 매개변수 ε의 선택이 k 증가에 따라 수렴 강건성에 어떻게 영향을 미치는가?
• RQ4 부분영역에 임피던스 조건을 사용하는 것이 딜리클레 조건보다 조절자 효율성 측면에서 뚜렷하게 우수한가?
• RQ5 ε ~ k 인 흡수 문제를 순수 헬름홀츠 문제에 대한 강건한 조절자로 사용할 수 있으며, 이는 이론적으로도 뒷받침되는가?

주요 결과

• 수치 실험 결과, H ~ k^(-α) (α ∈ [0, 0.4]) 일 때 k-독립적인 GMRES 수렴이 관찰되어 흡수 문제에 대한 강건성을 확인한다.
• |ε| ≤ Ck² 일 때, 조절된 행렬의 노름에 대한 k-독립적인 상한 경계가 증명되었으며, 연구된 흡수 문제 전반에 걸쳐 유효하다.
• ε/k 가 일정하거나 임의로 천천히 증가할 경우, 필드 값의 원점으로부터의 거리에 대한 엄밀히 양수인 하한 경계가 확립되어 강건한 수렴을 보장한다.
• 조절자에서 흡수를 제거하더라도(ε_prec = 0) 조절자가 여전히 효과적이지만, ε_prec = k 인 경우보다 수렴이 약간 덜 강건하다.
• 부분영역에 딜리클레 조건을 사용하는 것보다 임피던스 조건을 사용할 경우, 특히 큰 k에서 수렴성이 뚜렷이 열등하여, 임피던스 조건의 핵심적 역할을 입증한다.
• 놀랍게도, GMRES에서 무작위 초깃값이 영 초깃값보다 더 적은 반복 수를 요구하는 것으로 나타났지만, 그 이유는 제공되지 않았다.