[논문 리뷰] A Chebyshev-based rectangular-polar integral solver for scattering by general geometries described by non-overlapping patches

Oscar P. Bruno, Emmanuel Garza|arXiv (Cornell University)|2018. 07. 04.

Electromagnetic Scattering and Analysis참고 문헌 19인용 수 9

한 줄 요약

이 논문은 일반적인 기하구조를 비중첩적인 논리적 사각형 패치로 기술하는 음향 산란 문제를 위한 고차정확도의 직각좌표-극좌표 적분 해법을 제시한다. 체비셰프 콜로케이션과 페지어 적분을 기반으로 하며, 표면에 따라 변하는 직각좌표-극좌표 변환, 사전 계산된 특이 적분, 체비셰프 전개를 결합함으로써, 파장 당 몇 점만으로도 스펙트럴 정확도를 달성한다. 이는 개방형 및 폐쇄형 표면, 복잡한 CAD 생성 기하구조를 포함한 강건하고 빠르며 고정밀도 해법을 가능하게 한다.

ABSTRACT

This paper introduces a high-order-accurate strategy for integration of singular kernels and edge-singular integral densities that appear in the context of boundary integral equation formulations of the problem of acoustic scattering. In particular, the proposed method is designed for use in conjunction with geometry descriptions given by a set of arbitrary non-overlapping logically-quadrilateral patches---which makes the algorithm particularly well suited for treatment of CAD-generated geometries. Fejér's first quadrature rule is incorporated in the algorithm, to provide a spectrally accurate method for evaluation of contributions from far integration regions, while highly-accurate precomputations of singular and near-singular integrals over certain "surface patches" together with two-dimensional Chebyshev transforms and suitable surface-varying "rectangular-polar" changes of variables, are used to obtain the contributions for singular and near-singular interactions. The overall integration method is then used in conjunction with the linear-algebra solver GMRES to produce solutions for sound-soft open- and closed-surface scattering obstacles, including an application to an aircraft described by means of a CAD representation. The approach is robust, fast, and highly accurate: use of a few points per wavelength suffices for the algorithm to produce far-field accuracies of a fraction of a percent, and slight increases in the discretization densities give rise to significant accuracy improvements.

연구 동기 및 목표

• 경계 적분방정식에 나타나는 특이 및 가까운 특이 핵을 다루는 고차정확도 적분 방법을 개발하기 위해.
• 비중첩적인 논리적 사각형 패치로 기술된 복잡한 기하구조에서 효율적이고 정확한 산란 문제 해법을 가능하게 하기 위해.
• 최소한의 메esh 점 수(파장 당)로 스펙트럴 수렴과 고정밀도를 달성하여 계산 비용을 감소시키기 위해.
• GMRES 기반 해법에 통합하여 개방형 및 폐쇄형 표면에서의 산란 문제를 강건하고 확장 가능한 방식으로 해결하기 위해.
• 실제 공학 기하구조, 예를 들어 항공기 및 글라이더 CAD 모델에 대한 적용 가능성을 입증하기 위해.

제안 방법

• 이 방법은 하이브리드 니스트롬-콜로케이션 전략을 사용한다: 원거리 영역의 매끄러운 적분자는 페지어의 첫 번째 적분 규칙을 사용하여 스펙트럴 정확도를 달성한다.
• 특이 및 가까운 특이 상호작용의 경우, 이중 체비셰프 변환을 사용하여 핵-통합 밀도의 고정밀 사전 계산을 수행한다.
• 특이점과 기하 모서리 근처에 적분점의 집중을 가능하게 하기 위해 표면에 따라 변하는 직각좌표-극좌표 좌표변환을 적용한다. 이는 해상도와 정확도를 향상시킨다.
• 좌표변환은 특이점에서 도함수가 0이 되게 하여, 특수 적분 규칙이 필요 없이도 고차수 수렴을 가능하게 한다.
• 밀도의 체비셰프 다항식 전개를 통해 효율적인 미분과 기하 양의 계산이 가능해지며, 이는 전자기 BIE에서 유용하다.
• 페지어 적분과 체비셰프 변환 간의 공유 노드를 활용하여 이산 변환 계산을 단순화하고 오버헤드를 감소시킨다.

실험 결과

연구 질문

• RQ1 일반적인 비중첩 패치 기반 기하구조에서 경계 적분방정식에 나타나는 특이 핵을 다루는 고차정확도 적분 방법을 개발할 수 있는가?
• RQ2 모서리 및 기하 특이점을 다룰 수 있으며, 파장 당 최소한의 점 수로도 스펙트럴 정확도를 달성할 수 있는가?
• RQ3 기하 재매핑이나 메esh 생성 없이도 복잡한 CAD 생성 기하구조에 얼마나 널리 적용될 수 있는가?
• RQ4 원거리 영역 정확도 및 해의 오차 측면에서 전진 연산자 및 전체 산란 해법의 수렴 행동은 어떠한가?
• RQ5 이 방법을 전자기 산란 및 음향 산란 이외의 적분방정식 설정으로 확장할 수 있는가?

주요 결과

• 파장 당 몇 점만으로도 원거리 영역 정확도가 백분율의 일부 수준으로 달성되어 높은 효율성과 강건성을 입증한다.
• 크기 2×2×2인 큐브에서 k=1일 때, 모서리 특이점에 맞는 좌표변환을 적용하면 원거리 오차가 수개월 수준으로 감소하며, p=2일 때 유의미한 향상이 이루어진다.
• 반지름 1인 원판에서 k=10π(λ=0.2)일 때, 이 방법은 포아송 반점 현상을 정확히 포착하여 물리적 신뢰성과 미세 구조의 해상도를 확인한다.
• 개방형 및 폐쇄형 표면 모두에서 고차수 수렴을 달성하며, 매끄럽고 특이한 적분자에 대해 이론적 기대와 일치하는 수렴 속도를 보인다.
• 87개의 패치로 구성된 CAD 기반 글라이더에 대해, 이 방법은 총장 및 강도 패턴을 성공적으로 계산하여 실제 공학 기하구조 적용 가능성을 입증한다.
• 체비셰프 기반 사전 계산과 직각좌표-극좌표 변환의 활용으로 성능이 향상되어, 과도한 메쉬 세분화 없이도 고정밀도를 달성할 수 있다.