$hp$-a posteriori error estimates for hybrid high-order methods applied to biharmonic problems

문제 정의
본 논문은 클램프 경계 조건($u = \partial_n u = 0$ on $\partial\Omega$)을 갖는 비하모닉 문제($\Delta^2 u = f$)에 적용된 하이브리드 고차(Hybrid High-Order, HHO) 방법에 대한 사후 오차 추정(a posteriori error estimation)을 다룬다. 이 문제는 일반적인 위상(구멍을 포함할 수 있음)을 가진 2차원 및 3차원 다면체 리프시츠(Lipschitz) 영역에 대해 설정된다. 목표는 메쉬 크기($h$), 다항식 차수($p$ 또는 $k$), 그리고 도메인의 위상에 대해 강건한(robust) 잔차 기반 $hp$-오차 추정기를 유도하는 것이다.

방법론
저자들은 전체 오차를 일치(conforming) 성분과 비일치(nonconforming) 성분으로 분해함으로써 오차를 바운드하는 프레임워크를 개발한다.

1. 오차 분해: 오차는 일치 부분(재구성된 $H^2_0$ 함수의 에너지 노름으로 측정됨)과 비일치 부분(이산 해와 일치 공간 사이의 거리)으로 나뉜다.
2. 비일치 오차 제어:
   • 비일치 오차를 바운드하기 위해, 저자들은 $C^1$-분할 구(partition of unity)를 사용하여 이산 해의 일치하는 $H^2_0$ 재구성을 구성한다.
   • 이 분할은 $C^1$-일치 유한 요소 공간(구체적으로 [40]의 Alfeld 분할에 기반한 구성)을 사용하여 구축된다.
   • 결정적으로, 저자들은 정점 패치(vertex patches, stars) 상에서 **국소 헬름홀츠 분해(local Helmholtz decompositions)**를 채택한다. 기존의 접근 방식들이 도메인의 구멍 개수에 의존하는 상수를 도입하는 전역 헬름홀츠 분해에 의존했던 것과 달리, 이 국소적 접근 방식은 안정성 상수가 메쉬의 형상 정규성(shape-regularity)에만 의존하도록 보장한다.
3. 일치 오차 제어:
   • 두 가지 서로 다른 상한(upper bounds)이 각각 특정 보간 연산자를 사용하여 유도된다:
     • 설정 1: 정형 하이브리드 유한 요소 보간 연산자($C^{k+2}_T$)를 사용한다. 이 접근 방식은 이 연산자와 HHO 재구성의 합성이 $H^2$-타원 투영(elliptic projection)을 생성한다는 성질에 의존한다. 결과적인 바운드는 안정화 항과 데이터 진동(data oscillation)만을 포함하지만, 이는 보간 오차에 관한 수치적 가정(가정 3.4)에 달려 있다.
     • 설정 2: 수정된 Babuška–Suri $hp$-보간 연산자($I^{k+2}_{mBS}$)를 사용한다. 이 접근 방식은 추가적인 가정을 필요로 하지 않지만, 벌크 잔차(bulk residuals) 및 법선/접선 플럭스 점프(normal/tangential flux jumps)와 같은 추가적인 오차 지표 항들을 도입한다.
4. 잔차 추정기: 저자들은 적절한 $hp$-인자(factor)로 스케일링된 안정화 항, 벌크 잔차, 법선 플럭스 점프, 그리고 접선 점프를 포함하는 국소 오차 지표를 정의한다.

주요 기여

• 비하모닉 문제에 대한 HHO를 위한 첫 번째 $hp$-추정기: 본 연구는 일반적인 2D 및 3D 도메인에서의 비하모닉 문제를 위한 HHO 방법의 첫 번째 잔차 기반 $hp$-사후 오차 추정기를 제공한다.
• 위상적 강건성을 갖춘 국소 헬름홀츠 분해: 본 논문은 정점 스타(vertex stars) 상에서의 국소 헬름홀츠 분해에 대한 새로운 안정성 추정치를 확립한다. 저자들은 안정성 상수가 오직 메쉬의 형상 정규성에만 의존함을 증명함으로써, 전역 분해 접근 방식에서 나타나는 도메인의 구멍 개수에 따른 상수 저하 문제를 방지하였다.
• 이중 상한(Dual Upper Bounds): 일치 오차에 대한 두 가지 별개의 상한을 유도함으로써 유연성을 제공한다. 한 바운드는 "깔끔하며"(안정화 및 진동에만 의령함) 수치적으로 검증된 가정에 조건부인 반면, 다른 하나는 무조건적이지만 더 복잡한 잔차 항들을 포함한다.
• 일반 위상으로의 확장: 방법론은 임의의 위상(여러 개의 구멍)을 가진 도메인을 처리하도록 명시적으로 설계되었으며, 이는 4차 문제의 $hp$-추정에서 다루기 어려운 특징이다.

결과

• 이론적 바운드: 주요 정리(정리 5.6)는 전체 오차(broken Hessian norm 및 안정화 항)가 제안된 오차 지표들의 합에 의해 바운드됨을 확립한다. 이 바운드들은 $h$-최적(optimal)이며 $p$-서브옵티멀(suboptimal)하다(검증된 가정 하에서는 최대 1차 또는 잠재적으로 $1/2$ 차이만큼).
• 수치적 검증:
  • 가정 검증: 2D 및 3D에서의 수치 실험은 정형 하이브리드 연산자의 보간 오차에 관한 가정 3.4의 타당성을 확인한다.
  • 수렴 속도: 매끄러운 해에 대한 3D 테스트는 에너지 오차와 추정기가 최적의 속도 $O(\text{DoFs}^{-(k+1)/d})$로 수렴함을 보여준다.
  • 효과 지수(Effectivity Indices): 추정기는 다양한 다항식 차수와 메쉬 세분화에 걸쳐 완만한 효과 지수(대략 1.2에서 2.0 사이)를 나타낸다.
  • 위상적 강건성: 구멍의 개수가 증가하는(0, 1, 4, 8개) 2D 도메인에 대한 실험은 효과 지수가 안정적으로 유지됨을 보여주며, 이는 추정기가 도메인 위상에 대해 강건하다는 이론적 주장을 확인시켜 준다.
  • 적응적 성능: 추정기에 의해 구동되는 $h$-적응 알고리즘은 2D(L-자형 도메인) 및 3D(점 특이점)의 특이 해에 대해 균등 세분화보다 우수한 성능을 보이며 최적의 수렴 속도를 성공적으로 회복한다.

의의
저자들은 비하모닉 문제를 위한 $hp$-설정에서 비일치 오차를 바운드하기 위해$C^1$-분할 구와 국소 헬름홀츠 분해를 결합하여 사용하는 첫 번째 연구로 본 작업을 위치시킨다. 주요 의의는 도메인의 위상적 복잡성(구멍의 개수)에 따라 저하되는 안정성 상수 없이 비일치 오차를 제어할 수 있는 능력에 있다. 또한, 저자들은 이 접근 방식이 HHO 방법을 위해 개발되었으나, 국소 헬름홀츠 분해를 통해 비일치 오차를 바운드하는 기술은 Interior Penalty Discontinuous Galerkin (IPDG) 방법과 같은 다른 비일치 방법에도 잠재적으로 적용 가능하다고 언급한다. 본 논문은 하한(lower error bounds)을 제공한다고 주장하지 않으며, 그러한 하한은 기존 논거의 변형을 통해 도출될 수 있으나 특정 버블 함수(bubble functions) 및 확장 연산자(extension operators)를 구성하기 위한 현재의 $hp$-기술의 한계로 인해 비관적인$p$-서브옵티멀리티를 보일 가능성이 높다고 명시하였다.