[논문 리뷰] Coupling of HDG with a double-layer potential BEM

Zhixing Fu, Norbert Heuer외 1명|arXiv (Cornell University)|2012. 12. 10.

Advanced Numerical Methods in Computational Mathematics참고 문헌 23인용 수 4

한 줄 요약

이 논문은 유한하지 않은 영역에서의 전달 문제를 해결하기 위해 이중층 잠재력 표현에 기반한 비대칭 결합 방법을 제안한다. 이는 하이브리드 불연속 갈레르킨(HDG) 방법과 경계 요소법(BEM) 간의 새로운 비대칭 결합 방법이다. 이 방법은 이산적으로 균일한 강제성(coercivity)을 보장하며, 모든 변수에 대해 최적 수렴률 $\mathcal{O}(h^{k+1})$을 달성하고, 특정 필드에서는 초수렴(superconvergence)을 $\mathcal{O}(h^{k+2})$로 달성한다. 이는 수치 실험을 통해 검증되었다.

ABSTRACT

In this paper we propose and analyze a new coupling procedure for the Hybridizable Discontinuous Galerkin Method with Galerkin Boundary Element Methods based on a double layer potential representation of the exterior component of the solution of a transmission problem. We show a discrete uniform coercivity estimate for the non-symmetric bilinear form and prove optimal convergence estimates for all the variables, as well as superconvergence for some of the discrete fields. Some numerical experiments support the theoretical findings.

연구 동기 및 목표

• 유계 영역 내에서 계수 변화가 있는 확산 문제와 외부 영역에서의 라플라스 방정식을 포함하는 전달 문제에 대해 안정적이고 비대칭적인 HDG와 BEM 간의 결합을 개발한다.
• 특히 계수 변화가 있는 확산계수에 대해 비대칭 HDG-BEM 형식에서 강제성을 증명하는 분석적 과제를 해결한다.
• HDG 및 BEM 변수 전역에 대해 최적 수렴률을 확립하며, 스칼라 및 유량 근사에서의 초수렴을 포함한다.
• HDG의 장점을 유지하는 프레임워크를 제공한다. 예를 들어, 자유도 수를 줄이고 후처리 가능성을 제공하며, BEM을 통해 비유한 영역로의 확장을 가능하게 한다.

제안 방법

• 외부 영역의 해를 경계면 상의 단일층 잠재력으로 표현하기 위해 이중층 잠재력 표현을 사용하여, 외부 문제를 경계 적분 방정식으로 축소한다.
• 스켈레톤 상에서 스칼라 필드, 유량, 트레이스 미지수에 대해 차수 $k \geq 0$의 다항식을 사용하여 내부 영역에 대해 HDG 방법을 적용한다.
• BEM 구성 요소는 이중층 잠재력 연산자 $\mathcal{W}$를 사용하며, 이는 $H^{1/2}(\Gamma)$에서 유계이자 강제성 있는 것으로 밝혀져 안정성 분석이 가능하다.
• HDG 약한 형태와 BEM 적분 방정식을 결합하여 비대칭 변분 형식을 구성함으로써 단일 비대칭 이차형식을 도출한다.
• 이산 강제성은 메쉬 크기 $h$와 다항식 차수 $k$에 대해 균일하게 증명되며, 이는 이중층 연산자의 성질과 트레이스 추정치에 기반한다.
• HDG 해에 대한 후처리를 통해 차수 $\mathcal{O}(h^{k+2})$의 초수렴 근사를 회복할 수 있다.

실험 결과

연구 질문

• RQ1 계수 변화가 있는 확산 문제에 대해 이산 강제성과 최적 수렴률을 보장하는 비대칭 HDG-BEM 간의 결합을 구성할 수 있는가?
• RQ2 내부 영역의 확산 계수 $\kappa$의 크기가 비대칭 HDG-BEM 형식의 안정성과 조건 수에 미치는 영향은 어떠한가?
• RQ3 비대칭 결합에서 다른 BEM 형식에 비해 이중층 잠재력 표현을 사용할 경우 안정성과 수렴성 향상이 이루어지는가?
• RQ4 이중층 잠재력에 의한 BEM과의 결합 시 HDG 변수에서 초수렴이 유지되는가?
• RQ5 이러한 결합 형식에서 최적 수렴을 보장하기 위해 필요한 도메인과 해에 대한 최소 정규성 가정은 무엇인가?

주요 결과

• 이중층 잠재력에 기반한 제안된 HDG-BEM 결합은 비대칭 이차형식에 대해 이산적으로 균일한 강제성을 확보하여 이산 시스템의 잘 정의됨을 보장한다.
• 모든 HDG 변수에 대해 최적 수렴률 $\mathcal{O}(h^{k+1})$이 증명되었으며, 이는 스칼라 필드 $u_h$, 유량 $\mathbf{q}_h$, 트레이스 미지수 $\widehat{u}_h$에 해당한다.
• HDG 스칼라 및 유량 근사에서 차수 $\mathcal{O}(h^{k+2})$의 초수렴이 달성되어 고정밀도 후처리가 가능하다.
• 확산 계수 $\kappa$의 크기와 관계없이 $\kappa^{-1} \in L^\infty(\Omega)$를 만족하는 한 방법은 안정적이고 수렴성이 보장된다. 이는 일부 대칭 결합에 비해 뚜렷한 이점이다.
• 수치 실험을 통해 이론적 수렴률이 확인되었으며, 다각형 및 다면체 영역에서도 방법의 효과성이 입증되었다.
• 이중층 잠재력의 사용은 외부 해가 자연스럽게 복사 조건을 만족시키며, 비凸성 및 비연결 영역을 쉽게 처리할 수 있다.