[논문 리뷰] Stable and efficient FEM-BEM coupling with OSRC regularisation for acoustic wave transmission

Elwin van ’t Wout|arXiv (Cornell University)|2021. 07. 20.

Numerical methods in engineering참고 문헌 75인용 수 10

한 줄 요약

이 논문은 표면 방사 조건(On-Surface Radiation Condition, OSRC) 정규화와 OSRC-ILU 조절자(preconditioning)를 사용하여 음향파 전파 문제에 대해 안정적이고 효율적인 FEM-BEM 결합 방법을 제안한다. 표준 정규화 방법 대신 노이만-디리클레(Neumann-to-Dirichlet, NtD) 연산자의 OSRC 근사치를 사용하고, FEM 블록에 대해 비완전한 LUP 분해(ILU)를 결합함으로써, 일반적인 문제와 대규모 문제에서 GMRES 반복 횟수를 크게 감소시켰으며, 공명 주파수를 포함한 모든 주파수에서 안정성을 유지한다.

ABSTRACT

The finite element method (FEM) and the boundary element method (BEM) can numerically solve the Helmholtz system for acoustic wave propagation. When an object with heterogeneous wave speed or density is embedded in an unbounded exterior medium, the coupled FEM-BEM algorithm promises to combine the strengths of each technique. The FEM handles the heterogeneous regions while the BEM models the homogeneous exterior. Even though standard FEM-BEM algorithms are effective, they do require stabilisation at resonance frequencies. One such approach is to add a regularisation term to the system of equations. This algorithm is stable at all frequencies but also brings higher computational costs. This study proposes a regulariser based on the on-surface radiation conditions (OSRC). The OSRC operators are also used to precondition the boundary integral operators and combined with incomplete LU factorisations for the volumetric weak formulation. The proposed preconditioning strategy improves the convergence of iterative linear solvers significantly, especially at higher frequencies.

연구 동기 및 목표

• 음향 전파 문제에서 공명 주파수에서 표준 FEM-BEM의 불안정성을 해결한다.
• 기존 정규화 방법을 사용하는 안정화된 FEM-BEM 형식의 높은 계산 비용을 감소시킨다.
• 특히 고주파수에서의 반복적 해법(GMRES)의 수렴성을 향상시킨다.
• BEM에 대한 OSRC 연산자와 FEM에 대한 ILU를 조합한 주파수에 강인하고 효율적인 조절 전략을 개발한다.
• 밀도와 파동 속도가 변하는 비균질 물질을 통과하는 파동 전파의 정확한 대규모 시뮬레이션을 가능하게 한다.

제안 방법

• 모든 주파수에서 안정성을 확보하기 위해 안정화된 FEM-BEM 시스템에 노이만-디리클레(NtD) 연산자의 OSRC 근사치를 정규화 요소로 도입한다.
• BEM 블록의 경계 적분 연산자에 OSRC 연산자를 사용한 연산자 조절 전략을 적용하여 조건수를 개선한다.
• 체적 FEM 부분에 대해 OSRC 조절 전략과 비완전한 LUP 분해(ILU)를 조합하여 반복 선형 해법의 수렴 속도를 가속화한다.
• FEM과 BEM이 단일 시스템 행렬 내에서 동시에 해석되는 강한 결합 형식을 사용하며, 상호 연결 조건 항을 포함한다.
• 공간 이산화에 대해 P1 유한요소를 사용하고, 선형 시스템의 반복적 해법에 대해 수렴 기준 10⁻⁵를 갖는 GMRES를 적용한다.
• 구형체, 고리형 기하구조와 같은 정준 기하구조 및 비균질 물성 상수를 가진 대규모 상호접촉 고리 형상에 대해 방법을 검증한다.

실험 결과

연구 질문

• RQ1 OSRC 기반 정규화는 공명 주파수를 포함한 모든 주파수에서 표준 정규화 방법을 대체하여 FEM-BEM 결합의 안정성을 확보할 수 있는가? 이 과정에서 임의의 비물리적 해(solution)가 발생하지는 않는가?
• RQ2 특히 고주파수 영역에서 OSRC 조절 전략이 결합된 FEM-BEM 시스템의 GMRES 수렴 속도를 현저히 향상시키는가?
• RQ3 기존 조절 전략 대비 OSRC-ILU 조합 전략은 반복 횟수와 계산 비용 측면에서 어떻게 비교되는가?
• RQ4 제안된 방법은 밀도와 속도가 매끄럽지만 불연속적인 비균질 물질을 통과하는 고주파수 파동 전파를 정확하게 시뮬레이션할 수 있는가?
• RQ5 대규모 3D 음향 전파 문제에서 OSRC 안정화 및 조절된 FEM-BEM 방법의 확장성과 강인성은 어떠한가?

주요 결과

• 표준 정규화 방법 대신 NtD 연산자의 OSRC 근사치를 사용함으로써, OSRC-정규화 FEM-BEM 형식은 공명 주파수를 포함한 모든 주파수에서 안정성을 확보한다.
• 비균질 밀도와 파동 속도를 가진 여섯 개의 상호접촉 고리에 대한 대규모 시뮬레이션에서, OSRC-ILU 조절 전략은 GMRES 반복 횟수를 1781회로 감소시켰다.
• 해당 방법은 밀도와 속도가 매끄럽지만 불연속적인 물성 조건을 가진 3D 고리 형상과 같은 다양한 시험 케이스에서 강인한 성능을 보였다.
• 계산 벤치마크 결과, OSRC 조절 전략은 기존 형식 대비 고주파수 영역에서 수렴 성능을 현저히 향상시킴을 확인하였다.
• 정규화로 인해 발생하는 약한 조건수 문제를 효과적으로 다루며, 정규화로 인해 추가로 증가한 자유도를 포함한 시스템 행렬에서도 효율성을 유지하였다.
• 대규모 시뮬레이션 결과, 비균질 내부 영역과 무한 외부 영역을 포함한 전체 영역에서 정확한 파동장 재구성이 이루어졌으며, 이는 복잡한 음향 전파 문제 해결 능력을 검증한다.