Pressure-robust $hp$-a posteriori error estimates of $\boldsymbol{H}(\mathrm{div})$-conforming discontinuous Galerkin methods for the Stokes equations

문제 정의
본 논문은 2차원 및 3차원 다면체 립시츠(Lipschitz) 영역에서의 스토크스(Stokes) 방정식에 대한 사후 오차(a posteriori error) 추정을 다룬다. 구체적으로는 $H(\text{div})$-적합 불연속 갤러킨(dG) 방법(methods)에 초점을 맞춘다. 이러한 방법들에 대한 압력 강건한(pressure-robust) 사전 오차(a priori) 분석은 잘 확립되어 있는 반면, 그에 상응하는 사후 오차 추정기(error estimator)의 개발은 여전히 미개척 분야로 남아 있다. 고전적인 추정기들은 종종 점성($\mu$)의 역수($\mu^{-1}$)에 의해 스케일링된 압력 의존 항을 포함하며, 이는 압력이 지배적인 영역이나 무회전(irrotational) 포싱 항이 존재하는 경우 속도 오차를 크게 과대평가할 수 있다. 본 연구의 목표는 점성에 대해 강건하고 압력으로부터 독립적이며, 에너지 노름(energy norm)에서 속도 오차에 대한 전역 상한(global upper bound)과 국소 하한(local lower bound)을 모두 제공하는 잔차 기반(residual-based) $hp$-오차 추정기를 고안하는 것이다.

방법론
분석은 전체 오차를 적합(conforming) 부분과 비적합(nonconforming) 부분으로 분해하는 것에 기초한다.

1. 비적합 오차 제어: 저자들은 단위 분할(partition-of-unity) 프레임워크와 정점 패치(vertex patches)에서의 국소 헬름홀츠 분해(local Helmholtz decompositions)를 결합하여 사용한다. 핵심적인 기술적 혁신은 이러한 국소 헬름홀츠 분해에 대한 새로운 안정성 추정치를 유도한 것인데, 여기서 안정성 상수는 오직 메쉬의 형상 정규성(shape-regularity)에만 의존한다. 이 접근법은 $H(\text{div})$-적합 조각별 다항식과 해트 함수(hat function)의 곱을 발산(divergence)을 정확하게 보존하면서 $H^1_0(\omega_a)$로 들어 올리는(lift) 정교하게 설계된 패치별 스토크스 재구성(patchwise Stokes reconstruction)을 활용한다.
2. 적합 오차 제어: 적합 오차는 발산-무결성(divergence-free) 조건을 보존하는 보간 연산자를 구축하기 위해 일반화된 보고프스키(Bogovskiĭ) 연산자를 사용하여 분석된다. 저자들은 두 가지 구별된 접근 방식을 제안한다.
   • 접근 방식 1: 고전적인 바부슈카-수리(Babuška–Suri) 보간 연산자를 활용한다. 이는 다섯 개의 오차 지표와 데이터 진동(data oscillation) 항을 포함하는 상한을 산출한다. 네 개의 지표는 $p$-최적 스케일링을 보이지만, 점프 페널티(jump penalty) 항은 $p^{1/2}$만큼 차수가 낮은(suboptimal) 특성을 보인다.
   • 접근 방식 2: 교환(commuting) 및 직교(orthogonality) 성질을 갖는 네데렉(Nédélec) 투영 기반 보간 연산자에 의존한다. 이 접근 방식은 두 개의 잔차 지표와 데이터 진동 항을 갖는 상한을 산출하지만, 한 단계의 $p$-최적 차수를 희생한다.
3. 국소 하한: 4차 편미분 방정식 개발 기법에서 영감을 얻은 맞춤형 $C^1$-버블 함수(bubble functions)를 사용하여 압력 강건한 국소 하한을 확립한다.

주요 기여

• **최초의 압력 강건한 $hp$-추정기:** 저자들이 아는 바로는, 본 연구는 증명 가능한 상한과 국소 하한을 모두 특징으로 하는 스토크스 문제에 대한 최초의 엄밀한 압력 강건한$hp$-사후 오차 추정기이다.
• 국소 헬름홀츠 분해 안정성: 도메인의 구멍(hole) 개수에 의존하지 않는(전역 분해와 달리) 정점 패치에서의 국소 헬름홀츠 분해에 대한 안정성 추정치를 유도하였다.
• 새로운 스토크스 재구성: $H(\text{div})$-적합 다항식과 해트 함수의 곱을 발산을 정확히 보존하면서 $H^1_0$로 들어 올리는 패치별 스토크스 재구성을 제공한다.
• 두 가지 상한 전략: 적합 오차에 대해 오차 지표의 개수와 $p$-최적성 사이의 절충안을 제공하는 두 가지 구별된 압력 강건 상한을 제시한다.
• 개선된 하한: $H^2$-버블 함수를 사용하여 기존 결과([24]를 구체적으로 지칭)보다 개선된 압력 강건 국소 하한을 확립하였다.

결과
이론적 분석은 제안된 추정기가 다음을 제공함을 확인한다:

• 점성에 대해 강건하고 압력으로부터 독립적인 에너지 노름에서의 속도 오차에 대한 전역 상한.
• 압력 강건한 국소 하한.
• 2차원 및 3차원(매끄러운 해, 특이성을 가진 L-자형 도메인, 3D 코너 특이성 포함)에서의 수치 실험은 추정기의 신뢰성, 효율성 및 압력 강건성을 검증한다. 이 추정기에 의해 구동되는 적응형 스킴(adaptive schemes)은 특이성을 성공적으로 포착하고 최적의 수렴율을 달 수 있었던 반면, 비압력 강건 추정기는 압력 특이점 근처에서 적절하게 세분화(refine)하는 데 실패하였다.

의의 및 주장
본 논문은 제안된 방법론이 $H(\text{div})$-적합 dG 방법들에 국한되지 않고, 스코트-보겔리우스(Scott–Vogelius) 요소, 유리 함수(rational functions)를 사용하는 발산-무결 요소, 질량 보존 혼합 응력 방법, 격상(staggered) dG 방법 등 광범위한 압력 강건 이산화 방법들에 폭넓은 시사점을 준다고 주장한다. 저자들은 근저에 깔린 분해 및 재구성 기술이 독립적인 관심사이며, 하이브리드 불연속 갤러킨(HDG) 방법을 포함한 광범위한 비적합 유한 요소 방법들에 적용 가능함을 강조한다. 본 연구는 압력 강건한 스토크스 시뮬레이션을 위한 적응형 메쉬 세분화 전략의 이론적 토대를 제공함으로써 문헌의 중요한 공백을 메우며, 이를 통해 거대한 압력이나 무회전 힘이 지배하는 영역에서도 정확한 오차 제어가 보장되도록 한다.