[논문 리뷰] Inverse estimates for elliptic boundary integral operators and their application to the adaptive coupling of FEM and BEM

Markus Aurada, Michael Feischl외 4명|arXiv (Cornell University)|2012. 11. 19.

Advanced Numerical Methods in Computational Mathematics참고 문헌 34인용 수 5

한 줄 요약

이 논문은 타원 경계요소법(BEM)에서 네 가지 고전적 경계적분연산자에 대한 역추정(inverse estimates)을 수립하여, h-적응형 유한요소법-BEM 결합 알고리즘의 수렴 분석을 가능하게 한다. 주요 기여는 가중 잔차 오차 추정자에 기반한 적응형 FEM-BEM 방법의 수렴성을 증명한 것으로, 특별히 특이적으로 경사진 문제에서 최적의 수렴 속도로 수치적으로 검증되었다.

ABSTRACT

We prove inverse-type estimates for the four classical boundary integral operators associated with the Laplace operator. These estimates are used to show convergence of an h-adaptive algorithm for the coupling of a finite element method with a boundary element method which is driven by a weighted residual error estimator.

연구 동기 및 목표

• 라플라스 연산자와 관련된 네 가지 고전적 경계적분연산자에 대한 역유형 추정을 수립하기 위해.
• 특이 해를 갖는 선형 전송 문제에 대한 h-적응형 FEM-BEM 결합의 수렴 분석을 가능하게 하기 위해.
• FEM-BEM 결합 프레임워크에 데이터 근사 항을 포함한 신뢰할 수 있는 사후 오차 추정기를 확장하기 위해.
• 가중 잔차 오차 추정자에 기반한 적응형 알고리즘의 수렴성을 추정자 감소 기법을 사용하여 증명하기 위해.
• 이론적 결과를 수치 실험을 통해 검증하여, 적응형 메esh 세분화 하에서 최적의 수렴 속도를 보여주기 위해.

제안 방법

• 국소 타원 정규성 이론을 사용하여 단일층, 이중층, 초특이, 병합 잠재력 연산자에 대한 새로운 역추정을 유도한다.
• 이 추정을 적용하여 FEM-BEM 결합 맥락에서 가중 잔차 오차 추정자의 효율성과 신뢰성을 증명한다.
• AFEM에서의 추정자 감소 프레임워크를 FEM-BEM 설정으로 확장하여 오차 추정자의 수축성을 확립한다.
• FEM-BEM 결합을 위해 Costabel의 대칭 결합 공식을 사용하여 대칭성과 안정성을 보장한다.
• 마킹 파라미터 θ = 0.25를 사용하여 적응형 알고리즘을 구현하고, 균일 세분화와 결과를 비교한다.
• 이산 경계적분연산자 조립과 선형 시스템의 직접 해법을 위해 MATLAB의 HILBERT BEM 라이브러리를 사용한다.

실험 결과

연구 질문

• RQ1 일반적인 다면체 표면에서 네 가지 고전적 경계적분연산자에 대한 역추정을 유도할 수 있는가?
• RQ2 모서리 특이성이 있는 문제에 대해 적응형 FEM-BEM 결합 알고리즘이 수렴하는가?
• RQ3 FEM-BEM 결합에서 가중 잔차 오차 추정기를 데이터 근사 항까지 확장할 수 있는가?
• RQ4 추정자 감소 기반의 적응형 알고리즘이 최적의 수렴 속도를 달성할 수 있는가?
• RQ5 특이 해를 갖는 문제에서 적응형 세분화와 균일 세분화 전략의 수렴 거동은 어떻게 비교되는가?

주요 결과

• 논문은 국소 타원 정규성 이론을 사용하여 다면체 표면에서 네 가지 고전적 경계적분연산자에 대한 역추정을 수립하였다.
• θ = 0.25를 사용한 적응형 FEM-BEM 결합 알고리즘은 ℓ → ∞ 일 때 오차 추정자 ϱℓ → 0 으로 수렴한다.
• 수치 실험 결과 최적의 수렴 속도를 확인하였다: FEM의 경우 α = p/2, BEM의 경우 β = 3/2 + q이며, 여기서 p와 q는 다항식 차수이다.
• 균일 세분화는 재entrant 모서리에서의 특이성으로 인해 다항식 차수와 관계없이 항상 하위 최적 수렴 속도 α = 2/7 을 보였다.
• 적응형 세분화는 최적의 수렴 속도를 회복하였으며, p = 1, q = 1일 경우 전체 오차는 FEM 부분이 지배적으로 작용하여 α = 1/2 를 얻었다.
• 오차 추정자 ϱℓ 는 에너지 오차에 대한 신뢰할 수 있는 상한을 제공하며, 체적 기여와 경계 기여 모두 최적의 속도로 수렴한다.