[논문 리뷰] Integral equations requiring small numbers of Krylov-subspace iterations for two-dimensional penetrable scattering problems

Yassine Boubendir, Oscar P. Bruno외 2명|arXiv (Cornell University)|2013. 10. 04.

Electromagnetic Scattering and Analysis참고 문헌 21인용 수 4

한 줄 요약

이 논문은 부드러운 경계를 가진 2차원 전자기 및 음향 산란 문제에 대해 새로운 종류의 정규화된 병합 경계 적분 방정식을 제안한다. 복소 파수수를 가진 층 전위 또는 푸리에 기반 연산자를 적용함으로써, 이 방법은 전이 문제를 고유하게 해를 가진 프레드홀름 방정식 제2종으로 변환하여, 주어진 정밀도에서 쿠리로-서브스페이스 반복 횟수와 전체 계산 비용을 크게 감소시킨다.

ABSTRACT

This paper presents a class of boundary integral equations for the solution of problems of electromagnetic and acoustic scattering by two dimensional homogeneous penetrable scatterers with smooth boundaries. The new integral equations, which, as is established in this paper, are uniquely solvable Fredholm equations of the second kind, result from representations of fields as combinations of single and double layer potentials acting on appropriately chosen regularizing operators. As demonstrated in this text by means of a variety of numerical examples (that resulted from a high-order Nystrom computational implementation of the new equations), these "regularized combined equations" can give rise to important reductions in computational costs, for a given accuracy, over those resulting from previous boundary integral formulations for transmission problems.

연구 동기 및 목표

• 투과성 산란 문제의 경계 적분 방정식에서 일반적으로 나타나는 높은 반복 횟수 문제를 해결한다.
• 유일하게 해를 가지는 프레드홀름 방정식 제2종이 되는 경계 적분 형식을 개발한다. 이는 스펙트럼 조건을 개선한다.
• 고정밀 시뮬레이션에서 수렴하기 위해 필요한 쿠리로 반복 횟수를 최소화하여 계산 비용을 감소시킨다.
• 복소 파수수를 가진 층 전위 및 푸리에 기반 기호를 기반으로 한 정규화 연산자를 도입하여 특이한 연산자들을 안정화하고 정규화한다.
• 고차수 뉴스트롬 이산화 및 부드러운 2차원 투과성 산란체에 대한 수치 실험을 통해 새로운 형식의 효과를 입증한다.

제안 방법

• 정규화 연산자를 적용한 단일 및 이중 층 전위의 조합을 사용하여 필드 표현을 구성함으로써 적분 방정식을 안정화시킨다.
• 복소 파수수를 가진 정규화 연산자(예: $ \tilde{S}_{\tilde{\nu}} $)를 적용하여 특이한 핵을 컴act하고 잘 조건화된 연산자로 변환한다.
• 복소 파수수 층 전위의 고주파 점근적 행동과 일치하는 기호를 가진 푸리에 기반 정규화 연산자를 사용한다.
• 적절한 소볼레프 공간에서 제2종 프레드홀름 방정식의 체계로 문제를 제시하여 고유한 해를 보장한다.
• 전역 삼각함수 다항식 근사, 핵을 특이 및 스무스 부분으로 분할, 로그 특이성에 대한 명시적 적분 공식을 사용한 고차수 뉴스트롬 방법을 구현한다.
• 베셀 함수 및 하이퍼헬름홀츠 함수의 알려진 점근 전개를 활용하여 정규화된 연산자의 부드러움 성질을 분석한다.

실험 결과

연구 질문

• RQ1 정규화 연산자를 사용하여 전이 문제의 특이한 적분 연산자를 컴팩트하고 잘 조건화된 제2종 방정식으로 변환할 수 있는가?
• RQ2 기존의 형식에 비해 새로운 정규화된 병합 적분 방정식이 쿠리로-서브스페이스 반복 횟수를 크게 줄이는가?
• RQ3 새로운 적분 연산자의 스펙트럼 행동은 어떻게 되며, 이는 반복 해법의 수렴 속도에 어떤 영향을 미치는가?
• RQ4 복소 파수수 층 전위 또는 푸리에 기반 기호를 기반으로 한 정규화 연산자가 고유하게 해를 가지는 제2종 방정식을 유도하고 조건을 개선하는가를 엄밀히 증명할 수 있는가?
• RQ5 고차수 뉴스트롬 이산화가 부드러운 산란체에 대해 최적의 수렴 속도를 달성하고 계산 비용을 줄이는 데 얼마나 기여하는가?

주요 결과

• 제안된 정규화된 병합 적분 방정식은 적절한 소볼레프 공간에서 고유하게 해를 가지는 프레드홀름 방정식 제2종이며, 이는 안정적이고 수렴 가능한 반복 해법의 적용을 보장한다.
• 복소 파수수 층 전위를 정규화 연산자로 사용함으로써 시스템 행렬의 스펙트럼 성질이 향상되어 쿠리로-서브스페이스 반복 횟수가 극적으로 감소한다.
• 수치 실험 결과, 동일한 정밀도에서 기존 형식에 비해 새로운 형식이 훨씬 적은 반복 횟수를 요구함을 입증하였으며, 일반적으로 2배에서 5배 정도 감소한다. 이는 전체 계산 비용을 감소시킨다.
• 파arametric 단일 층 전위 연산자와 그 푸리에 기반 정규화된 대응 연산자 간의 차이는 제3차 수준의 부드러움 연산자이며, 이는 정규화의 효과가 스펙트럼 조건을 향상시키는 데 기여함을 정당화한다.
• 고차수 뉴스트롬 구현은 부드러운 경계에서 스펙트럼 정밀도를 달성하며, 로그 특이성에 대한 명시적 적분 공식이 고정밀도를 가능하게 한다.
• 정규화 연산자가 유도하는 유리한 스펙트럼 성질 덕분에, 이전의 형식에 비해 특히 고주파 산란 문제에서 뛰어난 효율성을 달성한다.