[논문 리뷰] Corrected Trapezoidal Rules for Boundary Integral Equations in Three Dimensions

Bowei Wu, Per‐Gunnar Martinsson|arXiv (Cornell University)|2020. 07. 06.

Electromagnetic Scattering and Analysis참고 문헌 30인용 수 6

한 줄 요약

이 논문은 3D 표면을 균일한 격자로 매개변수화한 경우 경계 적분 방정식을 해결하기 위해 수정된 사다리꼴 적분 규칙을 소개한다. 이는 계수 행렬의 대각성분만 수정하여 $O(h^3)$ 수렴을 달성하고, 아울러 아홉 점 스텐실을 사용하면 $O(h^5)$ 수렴을 달성한다. 표준 사다리꼴 규칙의 단순성은 유지하면서도 특이핵을 보정함으로써 라플라스 및 헬름홀츠 핵에 대해 고차 정확도를 달성한다.

ABSTRACT

The manuscript describes a quadrature rule that is designed for the high order discretization of boundary integral equations (BIEs) using the Nyström method. The technique is designed for surfaces that can naturally be parameterized using a uniform grid on a rectangle, such as deformed tori, or channels with periodic boundary conditions. When a BIE on such a geometry is discretized using the Nyström method based on the Trapezoidal quadrature rule, the resulting scheme tends to converge only slowly, due to the singularity in the kernel function. The key finding of the manuscript is that the convergence order can be greatly improved by modifying only a very small number of elements in the coefficient matrix. Specifically, it is demonstrated that by correcting only the diagonal entries in the coefficient matrix, $O(h^{3})$ convergence can be attained for the single and double layer potentials associated with both the Laplace and the Helmholtz kernels. A nine-point correction stencil leads to an $O(h^5)$ scheme. The method proposed can be viewed as a generalization of the quadrature rule of Duan and Rokhlin, which was designed for the 2D Lippmann-Schwinger equation in the plane. The techniques proposed are supported by a rigorous error analysis that relies on Wigner-type limits involving the Epstein zeta function and its parametric derivatives.

연구 동기 및 목표

• 약한 특이핵으로 인해 표준 사다리꼴 적분의 수렴 속도가 느려지는 3D 표면에서 경계 적분 방정식의 뉴스트롬 이산화 문제를 해결한다.
• 핵의 특이성을 보정하면서도 사다리꼴 규칙의 단순성을 유지하는 고차 정확도의 적분 체계를 개발한다.
• 특히 대각성분을 수정함으로써 소수의 행렬 원소만을 수정하여 $O(h^3)$ 및 $O(h^5)$ 수렴을 달성한다.
• 2D의 Duan-Rokhlin 적분 규칙을 주기성 또는 토러스 기하학을 가진 3D 표면으로 일반화한다.
• 특이 적분의 점근적 행동 분석을 위해 Wigner 유형의 극한과 에프스타인 제타 함수 이론에 기반한 엄밀한 오차 분석을 제공한다.

제안 방법

• 3D 표면(예: 변형된 토러스 또는 주기적 채널)을 균일한 격자로 매개변수화하고, 표준 사다리꼴 적분을 사용한 뉴스트롬 방법을 적용한다.
• 비대각성분은 그대로 유지하면서 계수 행렬의 대각성분만 수정하여 핵의 특이성을 반영한다. 이로 인해 비대각성분은 $\mathbf{K}(i,j) = \mathcal{K}(\mathbf{x}_i,\mathbf{x}_j)w_j$ 형태로 유지된다.
• 아홉 점 보정 스텐실을 사용하여 $O(h^5)$ 수렴을 달성함으로써 최소한의 계산 부담으로 정확도를 향상시킨다.
• 특이 적분의 점근적 행동 분석을 위해 에프스타인 제타 함수와 그 매개변수 도함수를 포함하는 Wigner 유형의 극한을 사용하여 오차 한계를 유도한다.
• 부드러운 컷오프 함수 $\phi$ 및 $\eta$를 사용한 정규화 기법을 적용하여 핵을 정규성분과 특이성분으로 분해함으로써 정밀한 오차 추정이 가능해진다.
• 푸리에 도메인에서 오차를 추정하기 위해 $S^{(\beta-j)}(\mathbf{u},s)$를 포함하는 적분을 유계로 제한하고, $\eta^{(\ell)}(\mathbf{u}h)$의 테일러 전개를 사용하여 나머지 항을 제어한다.

실험 결과

연구 질문

• RQ1 표준 사다리꼴 적분을 사용한 뉴스트롬 방법의 수렴 속도를 약한 특이핵을 가진 3D 경계 적분 방정식에서 크게 향상시킬 수 있는가?
• RQ2 특히 대각성분을 수정함으로써 소수의 행렬 원소만으로도 고차 정확도를 얼마나 달성할 수 있는가?
• RQ3 제안된 보정 체계는 기존의 2D 적분 규칙(예: Duan-Rokhlin 방법)을 3차원 표면으로 어떻게 일반화하는가?
• RQ4 보정된 적분 규칙의 이론적 오차 한계는 무엇이며, 메쉬 간격 $h$와 핵의 매끄러움 정도에 따라 어떻게 달라지는가?
• RQ5 특이 적분의 맥락에서 에프스타인 제타 함수와 그 도함수를 사용하여 오차 분석을 엄밀하게 정당화할 수 있는가?

주요 결과

• 뉴스트롬 행렬의 대각성분만 보정하면, 라플라스 및 헬름홀츠 핵을 가진 단일 및 이중층 잠재함수에 대해 $O(h^3)$ 수렴이 달성된다.
• 아홉 점 보정 스텐실을 사용하면 수렴 속도가 $O(h^5)$로 향상되어 최소한의 시스템 행렬 수정으로 정확도가 크게 향상된다.
• 비대각성분에 대해 표준 사다리꼴 적분의 효율성은 유지되어 FMM 및 빠른 직접 해법기와의 호환성이 보장된다.
• 오차 분석은 Wigner 유형의 극한과 에프스타인 제타 함수에 기반하여 특이 적분의 수렴 행동에 대한 엄밀한 기초를 제공한다.
• 보정의 차수 $K$에 대해 이론적 오차 한계는 $O(h^{2K+1})$이며, 컷오프 함수 도함수의 테일러 전개를 통해 나머지 항을 명시적으로 제어할 수 있다.
• 강한 특이핵을 가진 핵에 대해서도 고차 정확도 수렴이 달성되며, 이는 $\left|\sum_{|\mathbf{k}|\neq 0}\hat{f}_{\phi}(\mathbf{k},h) - Z_1(s)\right| \leq C h^{2K+1}$로 표현되며, 여기서 $Z_1(s)$는 에프스타인 제타 함수 합이다.