[논문 리뷰] Estimation of quadrature errors in layer potential evaluation using quadrature by expansion

Ludvig af Klinteberg, Anna‐Karin Tornberg|arXiv (Cornell University)|2016. 03. 28.

Electromagnetic Scattering and Analysis참고 문헌 12인용 수 13

한 줄 요약

이 논문은 경계 적분 방법에서 거의 특이적일 때의 층 잠재함수를 평가하기 위한 적분 확장(QBX) 기법에서의 적분 오차를 추정하기 위해 잔여물 미적분 기반의 방법을 제시한다. QBX 계수 계산에서 발생하는 거의 특이적 적분을 분석함으로써, 2차원 및 3차원에서 라플라스 및 헬름홀츠 단일층 잠재함수에 대해 점점 더 정확한 오차 추정치를 유도한다. 이는 실무에서 최적의 매개변수 선택을 가능하게 한다.

ABSTRACT

In boundary integral methods it is often necessary to evaluate layer potentials on or close to the boundary, where the underlying integral is difficult to evaluate numerically. Quadrature by expansion (QBX) is a new method for dealing with such integrals, and it is based on forming a local expansion of the layer potential close to the boundary. In doing so, one introduces a new quadrature error due to nearly singular integration in the evaluation of expansion coefficients. Using a method based on contour integration and calculus of residues, the quadrature error of nearly singular integrals can be accurately estimated. This makes it possible to derive accurate estimates for the quadrature errors related to QBX, when applied to layer potentials in two and three dimensions. As examples we derive estimates for the Laplace and Helmholtz single layer potentials. These results can be used for parameter selection in practical applications.

연구 동기 및 목표

• 경계 적분 방법에서 발생하는 거의 특이적 적분에서의 적분 오차를 정확하게 추정하는 데 도전한다.
• QBX 전개 계수 계산 중 발생하는 오차를 추정하기 위한 체계적인 방법을 개발한다.
• 2차원 및 3차원에서 라플라스 및 헬름홀츠 단일층 잠재함수에 대해 실용적이고 점근적으로 정확한 오차 추정치를 제공한다.
• 최적의 정확도를 얻기 위해 QBX 매개변수(예: 전개 차수 및 적분 점 수)의 정보 기반 선택을 가능하게 한다.
• 해석적 근사와 수치적 검증을 통해 복잡한 기하 구조(예: 타원체)로 오차 추정을 확장한다.

제안 방법

• 근사적 특이적 적분에 대한 점근적 오차 추정치를 도출하기 위해 윤곽 적분과 잔여물 미적분을 사용한다.
• 경계 적분 방법에서 일반적으로 사용되는 가우스-레지앙드 및 트라프레지oidal 적분 규칙에 이 방법을 적용한다.
• 2차원(로그함수) 및 3차원(역거리)의 단일층 잠재함수 커널에 대해 급수 전개와 구면 조화함수를 사용하여 오차 추정치를 도출한다.
• 2차원의 경우, 복소 로그 전개를 통한 계수 계산 오차를 분석한다. 3차원의 경우, 구면 조화함수 전개를 사용한다.
• 잔여물 기반 오차 공식의 지수 수렴 항의 기저를 최대화함으로써 가장 나쁜 경우의 오차 추정치를 도입한다.
• 원형, 타원형, 타원체 등의 모델 기하 구조에서 추정치를 검증하고, 표면 전체에 걸쳐 국소 오차 추정치를 통합함으로써 복잡한 형상으로 확장한다.

실험 결과

연구 질문

• RQ1 QBX에서 거의 특이적 층 잠재함수를 평가할 때 적분 오차를 어떻게 정확하게 추정할 수 있는가?
• RQ2 가우스-레지앙드 및 트라프레지oidal 규칙를 사용하여 QBX 전개 계수를 계산할 때 오차의 점근적 행동은 어떻게 되는가?
• RQ3 오차 추정치는 경계에서의 거리, 전개 반경, 축률 등의 기하 매개변수에 어떻게 의존하는가?
• RQ4 잔여물 미적분이 QBX에서 발생하는 거의 특이적 적분의 오차 추정에 신뢰할 수 있는 프레임워크를 제공할 수 있는가?
• RQ5 유도된 오차 추정치는 QBX 응용에서 실용적인 매개변수 선택을 어떻게 안내할 수 있는가?

주요 결과

• 잔여물 미적분 방법은 QBX에서 거의 특이적 적분에 대해 점근적으로 정확한 오차 추정치를 제공하며, 중간 크기의 적분 노드 수에서도 유효하다.
• 2차원 라플라스 단일층 잠재함수의 경우, 가우스-레지앙드 규칙의 오차는 $ | ext{R}_{n,t}^{ ext{G}}[ ilde{ ho}_l]| \lesssim \frac{B_l}{\Gamma(l+1/2)} \frac{\pi^2 c b}{r^{3/2}} \left(\frac{n b}{r}\right)^{l-1/2} e^{-2b n} $로 추정되며, 여기서 $ b \approx 2r/\pi c $이다.
• 3차원 타원체에서의 총 QBX 적분 오차는 $ |e_Q(\bm{x}_t)| \lesssim |\sigma(\bm{x}_t)| \sum_{l=0}^p r^l \left( |\text{I}_t \text{R}_{n,s}^\text{T}[\psi_l]| + |\text{I}_s \text{R}_{n,t}^\text{G}[\psi_l]| \right) $로 추정되며, 각 성분에 대해 명시적인 추정치가 있다.
• 조건 $ a \leq c $ 하에서 3차원의 최종 오차 추정치는 $ |\text{I}_s \text{R}_{n,t}^\text{G}[\psi_l]| \approx \frac{B_l}{\Gamma(l+1/2)} \frac{\pi^{5/2} c}{n} \left( \frac{2n}{c\pi} \right)^l e^{-4nr/c\pi} $로, $ n $에 따라 지수적으로 감소함을 보여준다.
• 편평한 타원체와 긴 타원체에서의 수치적 검증을 통해 유도된 추정치의 정확성이 확인되었으며, $ l \leq 10 $ 범위에서 양호한 일치를 보였다.
• 이 방법은 QBX에서 실용적인 오차 예측과 매개변수 조정을 가능하게 하며, 기하 구조와 커널에 따라 오차가 약 $ \mathcal{O}(n^{-1}) $ 비례하고 $ n $에 따라 지수적으로 감소함을 보여준다.